Figures
Abstract
Background
Natural extracts have played an important role in the prevention and treatment of diseases and are important sources for drug discovery. However, to be effectively used in these processes, natural extracts must be characterized through the identification of their active compounds and their modes of action.
Methodology/Principal Findings
From an initial set of 29,779 natural products that are annotated with their natural source and using a previously developed virtual screening procedure (carefully validated experimentally), we have predicted as potential peroxisome proliferators-activated receptor gamma (PPARγ) partial agonists 12 molecules from 11 extracts known to have antidiabetic activity. Six of these molecules are similar to molecules with described antidiabetic activity but whose mechanism of action is unknown. Therefore, it is plausible that these 12 molecules could be the bioactive molecules responsible, at least in part, for the antidiabetic activity of the extracts containing them. In addition, we have also identified as potential PPARγ partial agonists 10 molecules from 16 plants with undescribed antidiabetic activity but that are related (i.e., they are from the same genus) to plants with known antidiabetic properties. None of the 22 molecules that we predict as PPARγ partial agonists show chemical similarity with a group of 211 known PPARγ partial agonists obtained from the literature.
Conclusions/Significance
Our results provide a new hypothesis about the active molecules of natural extracts with antidiabetic properties and their mode of action. We also suggest plants with undescribed antidiabetic activity that may contain PPARγ partial agonists. These plants represent a new source of potential antidiabetic extracts. Consequently, our work opens the door to the discovery of new antidiabetic extracts and molecules that can be of use, for instance, in the design of new antidiabetic drugs or functional foods focused towards the prevention/treatment of type 2 Diabetes Mellitus.
Citation: Guasch L, Sala E, Mulero M, Valls C, Salvadó MJ, Pujadas G, et al. (2013) Identification of PPARgamma Partial Agonists of Natural Origin (II): In Silico Prediction in Natural Extracts with Known Antidiabetic Activity. PLoS ONE 8(2): e55889. https://doi.org/10.1371/journal.pone.0055889
Editor: Vladimir N. Uversky, University of South Florida College of Medicine, United States of America
Received: September 16, 2011; Accepted: January 7, 2013; Published: February 6, 2013
Copyright: © 2013 Guasch et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This study was supported by grants AGL2008-00387/ALI and AGL2011-25831/ALI from the “Ministerio de Educación y Ciencia” of the Spanish Government and the ACC1Ó (TECCT11-1-0012) program and XRQTC grant from the “Generalitat de Catalunya”. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Since ancient times, natural products (NPs) have played an important role in the treatment of type 2 diabetes mellitus (T2DM) [1]. Plants are one of the most important sources of antidiabetic compounds. Thus, 656 species from 437 genera, representing 111 plant families, with antidiabetic properties have been identified [1]. The plant families most studied as a result of their confirmed antidiabetic effects include Leguminoseae, Lamiaceae, Liliaceae, Cucurbitaceae, Asteraceae, Moraceae, Rosaceae, Euphorbiaceae and Araliaceae [2].
Although plant extracts have been used for the treatment of T2DM for hundreds of years in India [3], [4], China and other parts of the world, more research is needed for the identification of their active compounds and their mode of action. Some of the active principles associated with the antidiabetic activity of plant extracts are alkaloids, saponins, xanthones, flavonoids and nonstarch polysaccharides [1]. Despite the wide array of these active principles with a demonstrated antidiabetic activity, to date, metformin is the only drug approved for treatment of T2DM derived from a medicinal plant [5]. Therefore, the identification of the active compounds and the modes of action from plants traditionally used in the treatment of T2DM is an important issue for the discovery of new antidiabetic drugs and for the validation, standardization and rational use of traditional herbal remedies [1].
Numerous mechanisms of antidiabetic action have been proposed for several plant extracts [1], [6] and some hypotheses relate their effects to the increase of the insulin-stimulated glucose uptake. One target of interest for antidiabetic drugs is peroxisome proliferators-activated receptor gamma (PPARγ). PPARγ is a member of the nuclear receptor superfamily that regulate the gene expression of proteins involved in the control of glucose and lipid metabolism [7]. Indeed, the importance of PPARγ in regulating the insulin sensitivity has motivated research groups in both academia and the pharmaceutical industry to devote increasing efforts toward developing synthetic PPARγ agonists, which could be of therapeutic use in patients affected by T2DM [8]. Thiazolidinediones (TZDs) are one important class of synthetic agonists of PPARγ. TZDs are antidiabetic agents currently used in the treatment of T2DM that target adipose tissue and improve insulin sensitivity. Despite the clinical benefit of these drugs, the use of TZDs has been associated with adverse effects including weight gain, increased adipogenesis, renal fluid retention and possible increased incidence of cardiovascular events [9], [10]. Therefore, new PPARγ ligands with enhanced therapeutic efficacy and reduced adverse effects are needed. A promising new class of such ligands is selective PPARγ modulators (i.e., SPPARγMs) [9], [10]. These compounds act as partial agonists of PPARγ and display different binding properties in comparison to full agonists [11], [12]. Several natural products or plant extracts have been found to increase insulin-stimulated glucose uptake through the action of PPARγ with no or little effect on adipocyte differentiation [6], [13]–[15]. Thus, PPARγ partial agonists from natural extracts are promising candidates for the treatment of T2DM. There are successful examples of the application of structure-based drug design methods to discover new PPARγ partial agonists from natural products [16], [17].
Based on the hypothesis that it would be possible to identify PPARγ partial agonists among medicinal extracts previously used as hypoglycemic agents, the goal of the present work was to find natural extracts with known antidiabetic activity that contain at least one molecule that we predict as a PPARγ partial agonist through a virtual screening (VS) workflow that has previously been carefully validated experimentally [18]. Our results provide new information about potential active molecules of natural extracts with antidiabetic properties and their mode of action, i.e., the increase of the insulin-stimulated glucose uptake through the action of PPARγ. We also suggest plants with undescribed antidiabetic activity that may contain PPARγ partial agonists and are related to plants with known antidiabetic activity. These plants represent a potential new source of antidiabetic extracts. In addition, the new PPARγ partial agonists that we have predicted are chemically different from known PPARγ partial agonists and could be used as lead-hopping candidates for the development of new antidiabetic drugs.
Results and Discussion
Virtual Screening Description, Validation and Application
We used a slightly modified version of a VS workflow that was previously developed and validated experimentally [18] to identify PPARγ partial agonists from a large in-house database of compounds. Briefly, the VS used consists of a combination of two pharmacophore modeling methods (i.e., one of them to discard potential PPARγ full agonists and the second one to identify PPARγ partial agonists), a protein-ligand docking and an electrostatic and shape similarity search. The discriminatory power of the VS workflow to identify PPARγ partial agonists was evaluated by applying it to a group of 211 known PPARγ partial agonists obtained from the literature and to 3,122 decoys obtained from the DUD database [19]. See Table 1 for data about how many of these molecules survived each VS step. Because we were interested in discovering novel PPARγ partial agonists but not full agonists, we developed an initial structure-based pharmacophore, called the antipharmacophore, to exclude possible full agonists. We used this strategy because full agonists present more clearly defined features than partial agonists. Although both types of agonists interact with the ligand-binding domain of PPARγ through several hydrophobic contacts, their mode of binding, and thus their effects, are different [11], [12]. Full agonists are characterized by making a hydrogen-bond network with Ser289, Tyr473, His323 and His449 PPARγ residues, but most partial agonists form a hydrogen bond with Ser342 [11]. In total, 135 known PPARγ partial agonists and 2,204 decoys survived the antipharmacophore step, i.e., they were not identified as potential PPARγ full agonists and served as the input molecules in the next step (Table 1). From the molecules that survived the antipharmacophore step, 111 known PPARγ partial agonists and 964 decoys were identified as PPARγ partial agonists by our partial agonist pharmacophore. This represents an enrichment factor (EF) of 1.79 (Table 1). To find docking poses that were compatible with the partial agonist pharmacophore, the compounds that had at least one conformer, generated in vacuo, that matched with the partial agonist pharmacophore were also docked to the PPARγ structure from 2Q5S. The best docking poses were then matched again to the partial agonist pharmacophore, identifying that 72 out of 111 partial agonists and 382 out 964 decoys that survived the previous step have at least one docked pose that simultaneously accomplished the following: (a) compatibility with the PPARγ ligand-binding site; and (b) possession of functional groups that match the 3D location of the sites of the partial agonists pharmacophore. Finally an electrostatic and shape similarity analysis was applied. Using the experimental poses of five known PPARγ partial agonists as queries, 65 out of 72 partial agonists and 102 out 382 decoys were identified as partial agonist candidates by this VS step (Table 1). In terms of sensitivity (Se), specificity (Sp) and EF, the electrostatic/shape similarity analysis was the best step of the VS (Table 1). Overall, our VS workflow identified as partial agonists 65 and 102 out of the initial 211 and 3,122 molecules labeled as partial agonists and decoys, respectively. Therefore, the EF of the process was 6.15 (a 38.92% of 15.80 that would correspond to the highest possible EF value) and the Se and the Sp were 30.81% and 96.73%, respectively. The high Sp and moderate Se of our procedure reflect the correct assignment of inactive compounds and the loss of potential partial agonists, respectively. However, because of the high number of initial compounds and the difficulties in differentiating partial from full agonists, we preferred a specific, but less sensible, VS workflow. This VS workflow therefore seems adequate to identify molecules with antidiabetic properties that could act as PPARγ partial agonists.
Once the VS workflow was validated, it was applied to an in-house database formed by 29,779 NPs that contained an annotation of their natural source. After applying the VS workflow described above, a group of 65 PPARγ partial agonist candidates were ultimately identified (see Table 1 for viewing the number of molecules that survived each step of the VS workflow).
Virtual Screening Hits in Natural Extracts with Known Antidiabetic Activity
According to the information available in our in-house NP database, the 65 molecules that were predicted by the VS workflow as potential PPARγ partial agonists have been isolated from 74 different natural sources. Interestingly, a systematic bibliographic search of PubMed (http://www.pubmed.org) revealed that 11 out of these 74 natural extracts were described previously as having antidiabetic activity (Table 2). These 11 extracts contained 12 molecules that we predict to be PPARγ partial agonists (see Table 2 and Figure S1), therefore, it is plausible that they could contribute to the observed antidiabetic activity of their corresponding extracts. In fact, a search with SciFinder (http://www.cas.org/products/sfacad) revealed that 6 out of these 12 natural compounds are extremely similar to molecules for which antidiabetic properties have already been described (Table 2 and Figure 1), although no mechanism of action has been suggested for them. This finding validates our methodology and suggests that the mode of action of these molecules could be through PPARγ. The remaining 6 natural compounds not identified previously as antidiabetic molecules represent new molecules with this activity. The most significant compounds found in these 11 antidiabetic extracts will be discussed below:
Each row represents the comparison of the 2D chemical structure between a molecule predicted as a PPARγ partial agonist through our VS workflow and a similar molecule that has been described to present antidiabetic activity.
- The genus name Salvia derives from the Latin salvere meaning “to save” perhaps referring to the healing properties of plants from this genus. Leafs, roots or flowers from species of Salvia, like Salvia officinalis [20], Salvia miltiorrhiza (a regional Chinese variety) [21], Salvia fruticosa [22] and Salvia lavandulifolia [23] have been used traditionally worldwide to treat diabetes [1]. The molecules deoxyneocryptotanshinone and miltionone I, which are found in S. miltiorrhiza extracts (Table 2), were predicted by our VS as PPARγ partial agonists, and they are extremely similar to the main lipophilic diterpene compounds from Danshen (i.e., the dried root of S. miltiorrhiza), and in particular to tanshinone IIA (see Figure 1A for a comparison of the three structures). Tanshinone IIA enhances low-dose insulin-mediated tyrosine auto-phosphorylation of the insulin receptor β-subunit [24]. Although S. miltiorrhiza extracts have been shown to have anti-atherosclerotic and antidiabetic properties [21], there is not any evidence that relates the antidiabetic action of the extracts from S. miltiorrhiza with PPARγ. However, it is known that extracts from the leaves of S. officinalis activate PPARγ [6]. Deoxyneocryptotanshinone and miltionone I molecules may be useful for the development of a new class of specific insulin receptor activators that combine this action with the action of PPARγ partial agonists. In addition, we have predicted as PPARγ partial agonists three extra molecules from other extracts of Salvia whose species have never been described as antidiabetic: (a) sanigerone from Salvia lanigera; (b) 12-hydroxysapriparaquinone from Salvia prionitis and Salvia eriophora and (c) prionitin from S. prionitis (Table 3). These molecules are new candidates of PPARγ partial agonists.
- The 2′,5,6′-trihydroxy-6,7,8-trimethoxyflavone that is isolated from the roots of four species of plants of the genus Scutellaria (Scutellaria baicalensis, Scutellaria adenostegia, Scutellaria alpina and Scutellaria ramosissima), was also identified as a PPARγ partial agonist in our VS procedure (Table 2). Extracts from S. baicalensis are prescribed in Kampo medicines in Japan [1], and they are reported to enhance the antidiabetic activity of metformin [25]. Baicalein (5,6,7-trihydroxyflavone), a related compound of the flavone hit, isolated from the roots of S. baicalensis (see Figure 1B for a comparison of both structures), is an α-glucosidase inhibitor [26]. S. baicalensis extracts may therefore contain more than one active component with different modes of antidiabetic action.
- Cryptolepis sanguinolenta, a shrub indigenous to West Africa, has been employed by traditional healers in the treatment of various fevers, including malaria [27]. Cryptolepine, an indoloquinolone alkaloid isolated from C. sanguinolenta, significantly lowers glucose when given orally in a mouse model of diabetes [5], and its antihyperglycemic activity has been demonstrated by several cryptolepine analogs [28]. Cryptolepicarboline is a cryptolepine analog isolated from C. sanguinolenta [27] (see Figure 1C for a comparison of both structures) that we predict as a PPARγ partial agonist (see Table 2). This result suggests that the increase of glucose uptake caused by cryptolepine and analogous compounds could therefore be mediated by the action of PPARγ.
- Dysidine is a sesquiterpene quinone from the marine sponge Dysidea villosa that greatly promotes glucose uptake in 3T3-L1 cells and shows strong insulin-sensitizing activity [29]. The results of our VS procedure suggest that an analog of dysidine isolated from D. villosa (see Figure 1D for a comparison of both structures) may be a PPARγ partial agonist (Table 2). Although it has been suggested that dysidine exhibits its cellular effects through the activation of the insulin pathway, possibly through the inhibition of protein tyrosine phosphatases [29], it is possible that the mode of action of dysidine and analogous molecules could also be through the action of PPARγ, or that different components of a D. villosa extract show antidiabetic activity through different mechanisms. Dysidine and analogous molecules are therefore potential lead compounds for the discovery of new antidiabetic compounds.
- Xanthoangelol F from the Japanese plant Angelica keiskei significantly enhances glucose uptake without activating the transactivation activity of PPARγ [30]. This agrees with the results of our VS workflow that suggest that xanthoangelol F acts as a PPARγ partial agonist (Table 2). This compound may therefore belong to the interesting group of PPARγ partial agonists that stimulate glucose uptake without promoting the transactivation activity of PPARγ and avoid some of the problematic side effects of PPARγ full agonists [31], [32].
- The remaining 6 molecules predicted as PPARγ partial agonists through our VS workflow that belong to extracts with described antidiabetic properties (see Table 2) are: (a) 7-hydroxy-3,5,8-trimethoxyflavone from Achyrocline satureoidesm, a species widely used as medicinal plant in South America; (b) 5-hydroxy-7,8,2′,3′-tetramethoxyflavone from Andrographis paniculata; (c) 2,4,6,2′,6′-pentamethoxybiphenyl isolated from Fucus vesiculosus; (d) hericerin from Hericium erinaceum; (e) the molecule with CAS number 78279-81-9 from Aspergillus terreus; and (f) bazouanthrone from Harungana madagascariensis. Our results suggest that these molecules could be PPARγ partial agonists and that extracts containing these molecules could stimulate glucose uptake through the action of PPARγ. This information is novel and relevant because it is the first time that antidiabetic properties for these molecules have been suggested.
Taking into account the fact that extracts from closer species of the same genus may share a high number of components, we also look for species that contain a molecule that we predict as a PPARγ partial agonist and, although they have been not described previously as antidiabetic, they are related (i.e., they belong to the same genus) to species with known antidiabetic properties. Thus, we identified 10 molecules isolated from 16 different plants, such as Acradenia franklinii, Annona purpurea, Artocarpus gomezianus, Euodia lunuankenda, Evodia elleryana, Helichrysum mixtum, Helichrysum odoratissimum, Helichrysum stenopterum, Melicope ptelefolia, Melicope simplex, Murraya paniculata, Salvia eriophora, Salvia lanigera, Salvia prionitis, Swertia hookeri and Tephrosia watsoniana (Table 3 and Figure S1), whose extracts could show antidiabetic properties mediated by the action of PPARγ. For example, the 1,3,5,8-tetramethoxy-9H-xanthen-9-one from Swertia hookeri was identified as a PPARγ partial agonist by our VS (see Table 3). However, neither the molecule nor an extract from this species has been identified previously as an antidiabetic agent. Nevertheless, the whole plants of Swertia japonica and Swertia chirayita have been reported to exhibit hypoglycemic effects by oral administration, and the xanthone constituents, bellidifolin and methylswertianin, have been isolated as active constituents [33]–[35]. Methylswertianin and bellidifolin are molecules highly similar to the tetramethoxyxanthon from Swertia hookeri that we identified as a PPARγ partial agonist (see Figure 1E for a comparison of their structures). Our results therefore suggest that the antidiabetic action of Swertia species could be mediated at least in part by PPARγ.
To compare the 22 molecules from Tables 2 and 3 that we predict to be PPARγ partial agonists with known PPARγ partial agonists, we merged their structures with 211 structures of known PPARγ partial agonists obtained from the literature. The resulting set was classified into 26 clusters according to structure similarity. The 22 NP hits of our VS were classified into 12 clusters. None of these clusters contained any of the 211 known PPARγ partial agonists. Thus, our 22 predicted PPARγ partial agonists represent 12 different chemical scaffolds that are different from the ones present in known synthetic PPARγ partial agonists. Therefore, these scaffolds are lead-hoping candidates for searching for new PPARγ partial agonists.
Conclusions
We have applied an experimentally validated VS workflow based on (a) two structure-based pharmacophores, (b) protein-ligand docking and (c) an electrostatic/shape similarity analysis to identify NPs that may be novel scaffolds for the discovery of new PPARγ partial agonists. Thus, from an initial set of 29,779 NPs that are annotated with their natural source, we predict 22 molecules to be potential PPARγ partial agonists. A subset of 12 of these molecules are present in 11 natural extracts with known antidiabetic activity and 10 of them are present in extracts related (i.e., they are from species of the same genus) to plants with known antidiabetic activity. None of the 22 hits show chemical similarity with 211 known PPARγ partial agonists obtained from the literature and, therefore, are new chemical scaffold candidates for the development of PPARγ partial agonists. Moreover, our results provide a new hypothesis about the active molecules of natural extracts with antidiabetic properties and their mode of action, i.e., the insulin-stimulated glucose uptake is increased through the action of PPARγ. We also suggest plants with undescribed antidiabetic activity that may contain PPARγ partial agonists and are related to plants with known antidiabetic activity. These plants represent a new source of potential antidiabetic extracts. Consequently, our work opens the door to the discovery of new antidiabetic extracts and molecules that can be of use, for instance, in the design of new antidiabetic drugs or functional foods focused towards the prevention/treatment of T2DM.
Materials and Methods
Initial Dataset of Natural Compounds Used
The initial in-house dataset of natural compounds that was filtered through the VS contained 29,779 compounds annotated with the natural sources from which they were obtained and the bibliographic references that describe how to extract them from each natural source. Moreover, according to the FAF-Drugs2 program [36], all of these molecules (a) show good ADME properties according to the Lipinski rule of five [37] (i.e., only one violation of this rule was allowed) and (b) are not potentially toxic (i.e., they lack “warhead” chelators, frequent hitters, promiscuous inhibitors and other undesirable functional groups). Conformations and sites for the 3D structures of these 29,779 compounds were determined during the generation of the corresponding Phase v3.1 (Schrödinger LLC., Portland, USA; http://www.schrodinger.com) [38] databases with the Generate Phase Database graphic front-end. Conformers are generated using the ConfGen facility. ConfGen carefully and systematically selects which conformations to produce, based upon an examination of the structure of the ligand being processed. During conformation generation, the ligand is first divided into a core region and a periphery. The conformational search generates all core configurations and then varies the peripheral configurations. The parameter values used during this conformer generation were the default values (i.e., Rapid conformational sampling method and energy threshold of 25 kcal/mol for discarding conformers), with the exception of the maximum number of conformers per structure, which was increased from 100 (the default value) to 200. The conformer sites were generated with definitions made by adding the ability to consider aromatic rings as hydrophobic groups to the default built-in Phase definitions.
Virtual Screening Workflow
The VS workflow used in this work is a slightly modified version of a VS workflow developed previously (that was also validated experimentally) to identify PPARγ partial agonists in chemical databases [18]. Briefly, the VS workflow consisted of several steps that must be applied sequentially (i.e., the output molecules of one step were the input molecules for the next step). Thus, the filters applied (and sorted according their usage) were the following: (1) a structure-based antipharmacophore screening; (2) a structure-based pharmacophore screening (called partial agonist pharmacophore); and (3) an electrostatic/shape similarity analysis (the previously developed VS workflow was altered for the current work with lower threshold values for the electrostatic and shape comparisons; see below for more details). We have previously used a similar VS workflow to identify novel IKK-2 inhibitors [39], [40] and DPP-IV inhibitors [41], [42]. All of the PDB files used in that work were superposed with the DeepView v3.7 program (http://spdbv.vital-it.ch/) [43] to ensure that all of them had the same relative orientation. From then on, only the resulting re-oriented coordinates for these PDB files were used during the subsequent structure-based pharmacophore generation and in the steps of the VS workflow where spatial orientation is crucial (i.e., pharmacophore-based searches, protein-ligand docking studies and shape and electrostatic-potential comparisons).
The initial set of compounds was filtered by a structure-based antipharmacophore with the aim of discarding potential PPARγ full agonists. This pharmacophore is formed by 5 sites (two hydrogen-bond acceptors and three hydrophobic sites) that are present in most of the validated 19 complexes of full agonists (where validated means that the coordinates for the ligand and the PPARγ active site are reliable according to their corresponding electron density map) and is completed with receptor-based excluded volumes obtained from the PDB file coded as 1FM9. Thus, this filter removed from the sample those molecules that had at least one in vacuo-generated conformer that matched at least 4 out of 5 sites of the antipharmacophore. The fitting between the molecules and the pharmacophore was analyzed with Phase v3.1 [38]. The subset of molecules that did not match the antipharmacophore was then used to identify possible partial agonists. To accomplish this task, a second pharmacophore obtained from the common sites of 12 validated complexes between PPARγ, and a partial agonist was used. It consists of one hydrogen-bond acceptor and three hydrophobic sites with receptor-based excluded volumes obtained from the PDB file coded as 2Q5S. Molecules that had at least one in vacuo-generated conformer and that matched with the 4 sites of the partial agonist pharmacophore were initially identified as putative PPARγ partial agonists. To find docking poses that were compatible with the partial agonist pharmacophore, those molecules identified as putative PPARγ partial agonists were then docked (using a rigid protein and flexible ligand) to the ligand-binding site of 2Q5S (ligands and water molecules from this structure were removed prior to docking). Thus, the best 32 docked poses predicted by the eHiTS v2009 program (SimBioSys Inc., Toronto, Canada; http://www.simbiosys.ca/ehits) [44] were filtered again with Phase through the partial agonist pharmacophore, using the same filtering options as the first pharmacophore matching, except that no re-orientation of the poses was allowed during the search. Default docking conditions were selected with the exception of the size of the sides of the cubic box encompassing the PPARγ binding site, which was increased from 10 to 15 Å. By default, eHiTS systematically evaluates all possible protonation states for the receptor and ligand, automatically for every ligand pose. Then during the docking algorithm each state is evaluated and scored and the most favorable state is selected [44].
The aligned poses that passed the pharmacophore and docking screenings were submitted to an electrostatic/shape similarity analysis, using the PPARγ partial agonists crystallized in the structures 2G0H, 4PRG, 2Q5S, 2FVJ and 2Q6S as a queries. Both electrostatic and shape similarity analyses was performed with EON v2.0.1 (OpenEye Scientific Software, Inc., Santa Fe, New Mexico, USA; http://www.eyesopen.com). EON is an electrostatics comparison program that compares electrostatic potential maps of pre-aligned molecules and determines the Electrostatic Tanimoto combo (ET_combo) score as the similarity criteria. In our case, the molecules were pre-aligned at the pharmacophore and docking steps. As the electrostatic potential of molecules is not strongly linked to the molecular graph, EON can screen molecular databases for electrostatic similarity to a query compound. The ET_combo score is the sum of two calculations: (a) the Shape Tanimoto (ST) score, which is a quantitative measure of three-dimensional overlap (where 1 corresponds to a perfect overlap, i.e., the same shape) and (b) the Possion-Boltzman Electrostatic Tanimoto (ET_pb) score that compares the electrostatic potential of two small molecules and ranges from 1 (identical potential) to negative values that results from the overlap of positive and negative charges. In this work, we selected the EON thresholds taking into account the results of the comparison between a group of experimental poses for PPARγ partial agonists in their complexes with PPARγ. Applying the five query poses against twelve other experimental poses of previously described PPARγ partial agonists, the lowest values for the ET_pb score and ST were 0.2 and 0.4, respectively. Therefore, these values were used as thresholds during the VS electrostatic/shape similarity analysis. The thresholds used in the original VS workflow [18] were more restrictive (0.3 and 0.5 for ET_pb and ST, respectively).
Virtual Screening Workflow Validation
The ability of the VS workflow to identify PPARγ partial agonists was tested by applying it to a group of 211 known PPARγ partial agonists obtained from the literature and 3,122 decoys obtained from the DUD database [19]. The structures of the 211 partial agonists were built with ChemDraw Ultra v11.0 (CambridgeSoft Corporation, Cambridge, MA, USA; http://www.cambridgesoft.com/) [45] and were cleaned using LigPrep v2.3 (Schrödinger LLC., Portland, USA; http://www.schrodinger.com). We calculated an EF and values for sensitivity (Se) and specificity (Sp) for the global VS workflow and each step [46]. The EF was obtained as the quotient between the fraction of actives in the sample that survived a particular VS step and the fraction of actives that were in the sample before applying this step. The EF therefore represents the ratio of the number of actives actually retrieved by a method compared to the number expected purely by chance. Se describes how well the model correctly identifies active compounds and it is calculated as the ratio between the number of active molecules that survived a particular VS step and the number of all active compounds that were in the sample before applying the VS step. Sp measures the correct assignment of inactive compounds. It is calculated as the ratio between the number of inactive compounds that were discarded at a particular VS step and the number of all inactive molecules that were in the sample before applying the VS step. Because the aim of the antipharmacophore step was to remove full agonists of PPARγ, for this step an EF and values of Se and Sp could not be estimated.
Structural Similarity Analysis
To obtain new scaffolds for PPARγ partial agonists, the VS hits were merged with the 211 PPARγ partial agonists previously used for validating the VS workflow and clustered with Canvas v1.2 (Schrödinger LLC., Portland, USA; http://www.schrodinger.com). Using a fingerprint precision of 32 bits, MOLPRINT2D fingerprints [47] were calculated for each molecule, and then a hierarchical clustering, based on Tanimoto similarities, was obtained. The number of clusters obtained was defined using the Kelley criterion [48].
Supporting Information
Figure S1.
2D-structures, product name and CAS number of the molecules predicted to be PPARγ partial agonists.
https://doi.org/10.1371/journal.pone.0055889.s001
(PDF)
Acknowledgments
This manuscript was edited for English language fluency by American Journal Experts. The authors wish to thank the “Servei de Disseny de Fàrmacs” (Pharmaceutical Design Service) of the Center for Scientific and Academic Services of Catalonia (CESCA) for providing access to the Schrödinger suite and OpenEye for providing access to its software.
Author Contributions
Conceived and designed the experiments: LG GP SG-V. Performed the experiments: LG ES. Analyzed the data: LG MM CV MJS GP SG-V. Wrote the paper: LG GP SG-V.
References
- 1.
Simmonds MSJ, Howes M-JR (2005) Plants Used in the Treatment of Diabetes. In: Soumyanath A, editor. Traditional Medicines for Modern Times. Antidiabetic Plants. Taylor & Francis. pp. 19–82.
- 2. Bnouham M, Ziyyat A, Mekhfi H, Tahri A, Legssyer A (2006) Medicinal plants with potential antidiabetic activity - A review of ten years of herbal medicine research (1990–2000). Int J Diabetes & Metab 14: 1–25.
- 3. Modak M, Dixit P, Londhe J, Ghaskadbi S, Devasagayam TP (2007) Indian herbs and herbal drugs used for the treatment of diabetes. J Clin Biochem Nutr 40: 163–173.
- 4. Grover JK, Yadav S, Vats V (2002) Medicinal plants of India with anti-diabetic potential. J Ethnopharmacol 81: 81–100.
- 5. Luo J, Fort DM, Carlson TJ, Noamesi BK, nii-Amon-Kotei D, et al. (1998) Cryptolepis sanguinolenta: an ethnobotanical approach to drug discovery and the isolation of a potentially useful new antihyperglycaemic agent. Diabet Med 15: 367–374 .
- 6. Christensen KB, Minet A, Svenstrup H, Grevsen K, Zhang H, et al. (2009) Identification of plant extracts with potential antidiabetic properties: effect on human peroxisome proliferator-activated receptor (PPAR), adipocyte differentiation and insulin-stimulated glucose uptake. Phytother Res 23: 1316–1325 .
- 7. Berger J, Moller DE (2002) The mechanisms of action of PPARs. Annu Rev Med 53: 409–435 .
- 8. Shearer BG, Billin AN (2007) The next generation of PPAR drugs: do we have the tools to find them? Biochim Biophys Acta 1771: 1082–1093 .
- 9. Feldman PL, Lambert MH, Henke BR (2008) PPAR modulators and PPAR pan agonists for metabolic diseases: the next generation of drugs targeting peroxisome proliferator-activated receptors? Curr Top Med Chem 8: 728–749.
- 10. Pourcet B, Fruchart J-C, Staels B, Glineur C (2006) Selective PPAR modulators, dual and pan PPAR agonists: multimodal drugs for the treatment of type 2 diabetes and atherosclerosis. Expert Opin Emerg Drugs 11: 379–401 .
- 11. Guasch L, Sala E, Valls C, Blay M, Mulero M, et al. (2011) Structural insights for the design of new PPARgamma partial agonists with high binding affinity and low transactivation activity. J Comput Aided Mol Des 25: 717–728 .
- 12. Guasch L, Sala E, Valls C, Mulero M, Pujadas G, et al. (2012) Development of docking-based 3D-QSAR models for PPARgamma full agonists. J Mol Graph Model 36: 1–9 .
- 13. Dey D, Medicherla S, Neogi P, Gowri M, Cheng J, et al. (2003) A novel peroxisome proliferator-activated gamma (PPAR gamma) agonist, CLX-0921, has potent antihyperglycemic activity with low adipogenic potential. Metabolism 52: 1012–1018.
- 14. Puhl AC, Bernardes A, Silveira RL, Yuan J, Campos JLO, et al. (2012) Mode of peroxisome proliferator-activated receptor γ activation by luteolin. Mol Pharmacol 81: 788–799 .
- 15. Weidner C, De Groot JC, Prasad A, Freiwald A, Quedenau C, et al. (2012) Amorfrutins are potent antidiabetic dietary natural products. Proc Natl Acad Sci U S A 109: 7257–7262 .
- 16. Fakhrudin N, Ladurner A, Atanasov AG, Heiss EH, Baumgartner L, et al. (2010) Computer-aided discovery, validation, and mechanistic characterization of novel neolignan activators of peroxisome proliferator-activated receptor gamma. Mol Pharmacol 77: 559–566 .
- 17. Petersen RK, Christensen KB, Assimopoulou AN, Fretté X, Papageorgiou VP, et al. (2011) Pharmacophore-driven identification of PPARγ agonists from natural sources. J Comput Aided Mol Des 25: 107–116 .
- 18. Guasch L, Sala E, Castell-Auví A, Cedó L, Liedl KR, et al. (2012) Identification of PPARgamma Partial Agonists of Natural Origin (I): Development of a Virtual Screening Procedure and In Vitro Validation. PLoS ONE 7: e50816 .
- 19. Huang N, Shoichet BK, Irwin JJ (2006) Benchmarking sets for molecular docking. J Med Chem 49: 6789–6801 .
- 20. Eidi M, Eidi A, Zamanizadeh H (2005) Effect of Salvia officinalis L. leaves on serum glucose and insulin in healthy and streptozotocin-induced diabetic rats. J Ethnopharmacol 100: 310–313 .
- 21. Tan Y, Kamal MA, Wang Z-Z, Xiao W, Seale JP, et al. (2011) Chinese herbal extracts (SK0506) as a potential candidate for the therapy of the metabolic syndrome. Clin Sci (Lond) 120: 297–305 .
- 22. Perfumi M, Arnold N, Tacconi R (1991) Hypoglycemic activity of Salvia fruticosa Mill. from Cyprus. J Ethnopharmacol 34: 135–140.
- 23. Jimenez J, Risco S, Ruiz T, Zarzuelo A (1986) Hypoglycemic activity of Salvia lavandulifolia. Planta Med 260–262.
- 24. Jung SH, Seol HJ, Jeon SJ, Son KH, Lee JR (2009) Insulin-sensitizing activities of tanshinones, diterpene compounds of the root of Salvia miltiorrhiza Bunge. Phytomedicine 16: 327–335 .
- 25. Waisundara VY, Hsu A, Huang D, Tan BK-H (2008) Scutellaria baicalensis enhances the anti-diabetic activity of metformin in streptozotocin-induced diabetic Wistar rats. Am J Chin Med 36: 517–540.
- 26. Nishioka T, Kawabata J, Aoyama Y (1998) Baicalein, an alpha-glucosidase inhibitor from Scutellaria baicalensis. J Nat Prod 61: 1413–1415 .
- 27. Sharaf MHM, Schiff PL, Tackie AN, Phoebe CH, Howard L, et al. (1995) Submicromole structure elucidation: Cryptolepicarboline—a novel dimeric alkaloid from Cryptolepis sanguinolenta. Magn Reson Chem 33: 767–778.
- 28. Bierer DE, Dubenko LG, Zhang P, Lu Q, Imbach PA, et al. (1998) Antihyperglycemic activities of cryptolepine analogues: an ethnobotanical lead structure isolated from Cryptolepis sanguinolenta. J Med Chem 41: 2754–2764 .
- 29. Zhang Y, Li Y, Guo Y, Jiang H, Shen X (2009) A sesquiterpene quinone, dysidine, from the sponge Dysidea villosa, activates the insulin pathway through inhibition of PTPases. Acta Pharmacol Sin 30: 333–345 .
- 30. Enoki T, Ohnogi H, Nagamine K, Kudo Y, Sugiyama K, et al. (2007) Antidiabetic activities of chalcones isolated from a Japanese Herb, Angelica keiskei. J Agric Food Chem 55: 6013–6017 .
- 31. Choi JH, Banks AS, Estall JL, Kajimura S, Boström P, et al. (2010) Anti-diabetic drugs inhibit obesity-linked phosphorylation of PPARgamma by Cdk5. Nature 466: 451–456 .
- 32. Choi JH, Banks AS, Kamenecka TM, Busby SA, Chalmers MJ, et al. (2011) Antidiabetic actions of a non-agonist PPARγ ligand blocking Cdk5-mediated phosphorylation. Nature 477: 477–481 .
- 33. Tian L-Y, Bai X, Chen X-H, Fang J-B, Liu S-H, et al. (2010) Anti-diabetic effect of methylswertianin and bellidifolin from Swertia punicea Hemsl. and its potential mechanism. Phytomedicine 17: 533–539 .
- 34. Basnet P, Kadota S, Shimizu M, Namba T (1994) Bellidifolin: a potent hypoglycemic agent in streptozotocin (STZ)-induced diabetic rats from Swertia japonica. Planta Med 60: 507–511 .
- 35. Chandrasekar B, Bajpai MB, Mukherjee SK (1990) Hypoglycemic activity of Swertia chirayita (Roxb ex Flem) Karst. Indian J Exp Biol 28: 616–618.
- 36. Lagorce D, Sperandio O, Galons H, Miteva MA, Villoutreix BO (2008) FAF-Drugs2: free ADME/tox filtering tool to assist drug discovery and chemical biology projects. BMC Bioinformatics 9: 396 .
- 37. Lipinski CA, Lombardo F, Dominy BW, Feeney PJ (2001) Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings. Adv Drug Deliv Rev 46: 3–26.
- 38. Dixon SL, Smondyrev AM, Knoll EH, Rao SN, Shaw DE, et al. (2006) PHASE: a new engine for pharmacophore perception, 3D QSAR model development, and 3D database screening: 1. Methodology and preliminary results. J Comput Aided Mol Des 20: 647–671 .
- 39. Sala E, Guasch L, Iwaszkiewicz J, Mulero M, Salvadó M-J, et al. (2011) Identification of human IKK-2 inhibitors of natural origin (part I): modeling of the IKK-2 kinase domain, virtual screening and activity assays. PLoS ONE 6: e16903 .
- 40. Sala E, Guasch L, Iwaszkiewicz J, Mulero M, Salvadó M-J, et al. (2011) Identification of human IKK-2 inhibitors of natural origin (Part II): in Silico prediction of IKK-2 inhibitors in natural extracts with known anti-inflammatory activity. Eur J Med Chem 46: 6098–6103 .
- 41. Guasch L, Ojeda MJ, González-Abuín N, Sala E, Cereto-Massagué A, et al. (2012) Identification of novel human dipeptidyl peptidase-IV inhibitors of natural origin (part I): virtual screening and activity assays. PLoS ONE 7: e44971 .
- 42. Guasch L, Sala E, Ojeda MJ, Valls C, Bladé C, et al. (2012) Identification of novel human dipeptidyl peptidase-IV inhibitors of natural origin (Part II): in silico prediction in antidiabetic extracts. PLoS ONE 7: e44972 .
- 43. Guex N, Peitsch MC (1997) SWISS-MODEL and the Swiss-PdbViewer: an environment for comparative protein modeling. Electrophoresis 18: 2714–2723 .
- 44. Zsoldos Z, Reid D, Simon A, Sadjad SB, Johnson AP (2007) eHiTS: a new fast, exhaustive flexible ligand docking system. J Mol Graph Model 26: 198–212 .
- 45. Cousins KR (2011) Computer review of ChemDraw Ultra 12.0. J Am Chem Soc 133: 8388 .
- 46. Schuster D, Wolber G (2010) Identification of bioactive natural products by pharmacophore-based virtual screening. Curr Pharm Des 16: 1666–1681.
- 47. Duan J, Dixon SL, Lowrie JF, Sherman W (2010) Analysis and comparison of 2D fingerprints: insights into database screening performance using eight fingerprint methods. J Mol Graph Model 29: 157–170 .
- 48. Kelley LA, Gardner SP, Sutcliffe MJ (1996) An automated approach for clustering an ensemble of NMR-derived protein structures into conformationally related subfamilies. Protein Eng 9: 1063–1065.
- 49. Mesquita AAL, Corrêa DDB, De Pádua AP, Guedes MLO, Gottlieb OR (1986) Flavonoids from four compositae species. Phytochemistry 25: 1255–1256 .
- 50. Arredondo MF, Blasina F, Echeverry C, Morquio A, Ferreira M, et al. (2004) Cytoprotection by Achyrocline satureioides (Lam) D.C. and some of its main flavonoids against oxidative stress. J Ethnopharmacol 91: 13–20 .
- 51. Koteswara Rao Y, Vimalamma G, Rao CV, Tzeng Y-M (2004) Flavonoids and andrographolides from Andrographis paniculata. Phytochemistry 65: 2317–2321 .
- 52. Jahan N, Malik A, Muhammad P (2001) New Flavonoid from Mentha longifolia. Heterocycles 55: 1951–1955.
- 53. Wibudi A, Kiranadi B, Manalu W, winarto A, Suyono S (2008) The traditional plant, Andrographis paniculata (Sambiloto), exhibits insulin-releasing actions in vitro. Acta Med Indones 40: 63–68.
- 54. Nishimura R, Tabata K, Arakawa M, Ito Y, Kimura Y, et al. (2007) Isobavachalcone, a chalcone constituent of Angelica keiskei, induces apoptosis in neuroblastoma. Biol Pharm Bull 30: 1878–1883.
- 55. Matsuura M, Kimura Y, Nakata K, Baba K, Okuda H (2001) Artery relaxation by chalcones isolated from the roots of Angelica keiskei. Planta Med 67: 230–235 .
- 56. Arai K, Shimizu S, Yamamoto Y (1981) Metabolic Products of Aspregillus terreus. VI. Metabolites of the Strain IFO 8835. (3). The Isolation and Chemical Structures of Colorless Metabolites. Chem Pharm Bull 29: 1005–1012.
- 57. Dewi RT, Iskandar YM, Hanafi M, Kardono LBS, Angelina M, et al. (2007) Inhibitory effect of koji Aspergillus terreus on alpha-glucosidase activity and postprandial hyperglycemia. Pak J Biol Sci 10: 3131–3135.
- 58. Sharaf MHM, Schiff PL, Tackie AN, Phoebe CH, Johnson RL, et al. (1996) The isolation and structure determination of cryptomisrine, a novel indolo[3,2- b ]quinoline dimeric alkaloid from cryptolepis sanguinolenta. J Heterocycl Chem 33: 789–797 .
- 59. Stewart M, Fell PM, Blunt JW, Munro MHG (1997) Avarol and Related Compounds from the New Zealand Marine Sponge Dysidea sp. Aust J Chem 50: 341–348.
- 60. Glombitza K-W, Lentz G (1981) Antibiotics from algae—XXVIII: Cleavage of high molecular phlorotannin derivatives from the brown alga fucus vesiculosus L. Tetrahedron 37: 3861–3866 .
- 61. Lamela M, Anca J, Villar R, Otero J, Calleja JM (1989) Hypoglycemic activity of several seaweed extracts. J Ethnopharmacol 27: 35–43.
- 62. Lenta B, Ngouela S, Boyom F, Tantangmo F, Tchouya G, et al. (2007) Anti-plasmodial activity of some constituents of the root bark of Harungana madagascariensis LAM. (Hypericaceae). Chem Pharm Bull 55: 464–467.
- 63. Agbor GA, Kuate D, Oben JE (2007) Medicinal plants can be good source of antioxidants: case study in Cameroon. Pak J Biol Sci 10: 537–544.
- 64. Kimura Y, Nishibe M, Nakajima H, Hamasaki T, Shimada A, et al. (1991) Hericerin, a New Pollen Growth Inhibitor from the Mushroom Hericium erinaceum. Agric Biol Chem 55: 2673–2674.
- 65. Michihiro F, Nobuaki E, Ken'ichi M (2006) Utility and Function of Mushrooms Neuroprotective Effect of Hericium erinaceum. Foods & Food Ingred J Jpn 211: 141–147.
- 66. Ikeshiro Y, Hashimoto I, Iwamoto Y, Mase I, Tomita Y (1991) Diterpenoids from Salvia miltiorrhiza. Phytochemistry 30: 2791–2792 .
- 67. Ikeshiro Y, Mase I, Tomita Y (1989) Abietane type diterpenoids from Salvia miltiorrhiza. Phytochemistry 28: 3139–3141 .
- 68. Chemesova II, Iinuma M, Budantsev AL (1993) Investigation of the flavonoid composition of Scutellaria adenostegia. Chem Nat Compd 29: 133–134 .
- 69. Prabhakar PK, Doble M (2011) Mechanism of action of natural products used in the treatment of diabetes mellitus. Chin J Integr Med 17: 563–574 .
- 70. Chang FR, Wei JL, Teng CM, Wu YC (1998) Two new 7-dehydroaporphine alkaloids and antiplatelet action aporphines from the leaves of Annona purpurea. Phytochemistry 49: 2015–2018.
- 71. Shirwaikar A, Rajendran K, Dinesh Kumar C, Bodla R (2004) Antidiabetic activity of aqueous leaf extract of Annona squamosa in streptozotocin-nicotinamide type 2 diabetic rats. J Ethnopharmacol 91: 171–175 .
- 72. Likhitwitayawuid K, Sritularak B, De-Eknamkul W (2000) Tyrosinase inhibitors from Artocarpus gomezianus. Planta Med 66: 275–277 .
- 73. Fernando MR, Thabrew MI, Karunanayake EH (1990) Hypoglycaemic activity of some medicinal plants in Sri-Lanka. Gen Pharmacol 21: 779–782.
- 74. Jakupovic J, Kuhnke J, Schuster A, Metwally MA, Bohlmann F (1986) Phloroglucinol derivatives and other constituents from South African Helichrysum species. Phytochemistry 25: 1133–1142 .
- 75. Aslan M, Deliorman Orhan D, Orhan N, Sezik E, Yesilada E (2007) In vivo antidiabetic and antioxidant potential of Helichrysum plicatum ssp. plicatum capitulums in streptozotocin-induced-diabetic rats. J Ethnopharmacol 109: 54–59 .
- 76. Hänsel R, Cybulski E-M, Çubukçu B, Meriçli AH, Bohlmann F, et al. (1980) Neue pyron-derivate aus Helichrysum-arten. Phytochemistry 19: 639–644 .
- 77. Kamperdick C, Van NH, Sung TV, Adam G (1997) Benzopyrans from Melicope ptelefolia leaves. Phytochemistry 45: 1049–1056 .
- 78. Yeo J, Kang Y-M, Cho S-I, Jung M-H (2011) Effects of a multi-herbal extract on type 2 diabetes. Chin Med 6: 10 .
- 79. Briggs LH, Locker RH (1950) 486. Chemistry of New Zealand Melicope species. Part IV. Constituents of the bark of Melicope simplex. J Chem Soc 2376–2379 .
- 80. Kinoshita T, Wu J-B, Ho F-C (1996) Prenylcoumarins from Murraya paniculata var. omphalocarpa (Rutaceae): The Absolute Configuration of Sibiricin, Mexoticin and Omphamurin. Chem Pharm Bull 44: 1208–1211.
- 81. El-Lakany AM (2003) Two new diterpene quinones from the roots of Salvia lanigera Poir. Pharmazie 58: 75–76.
- 82. Li M, Zhang JS, Chen MQ (2001) A novel dimeric diterpene from Salvia prionitis. J Nat Prod 64: 971–972.
- 83. Ulubelen A, Birman H, Oksüz S, Topçu G, Kolak U, et al. (2002) Cardioactive diterpenes from the roots of Salvia eriophora. Planta Med 68: 818–821 .
- 84. Blasko G, Lin L-Z, Cordell GA (n.d.) Determination of a new tetracyclic diterpene skeleton through selective INEPT spectroscopy. J Org Chem 53: 6113–6115.
- 85. Ghosal S, Biswas K, Jaiswal DK (1980) Xanthone and flavonol constituents of Swertia hookeri. Phytochemistry 19: 123–126 .
- 86. Basnet P, Kadota S, Shimizu M, Takata Y, Kobayashi M, et al. (1995) Bellidifolin stimulates glucose uptake in rat 1 fibroblasts and ameliorates hyperglycemia in streptozotocin (STZ)-induced diabetic rats. Planta Med 61: 402–405 .
- 87. Negi JS, Singh P, Pant GJ, Rawat MSM (2010) RP-HPLC analysis and antidiabetic activity of Swertia paniculata. Nat Prod Commun 5: 907–910.
- 88. Gómez F, Quijano L, Calderón JS, Rodríquez C, Ríos T (1985) Prenylflavans from Tephrosia watsoniana. Phytochemistry 24: 1057–1059 .
- 89. Pavana P, Manoharan S, Renju GL, Sethupathy S (2007) Antihyperglycemic and antihyperlipidemic effects of Tephrosia purpurea leaf extract in streptozotocin induced diabetic rats. J Environ Biol 28: 833–837.