Home  |  About JAPTR |  Editorial board  |  Search |  Ahead of print  |  Current issue  |  Archives |  Submit article  |  Instructions  |  Subscribe  |  Advertise  |  Contacts  |Login 
Users Online: 1201   Home Print this page Email this page Small font sizeDefault font sizeIncrease font size
     

 Table of Contents  
ORIGINAL ARTICLE
Year : 2015  |  Volume : 6  |  Issue : 1  |  Page : 13-18  

Three-dimensional quantitative structure-activity relationships and docking studies of some structurally diverse flavonoids and design of new aldose reductase inhibitors


1 Department of Chemistry, Tripura University, Agartala, Tripura, India
2 Department of Chemistry, MBB College, Agartala, Tripura, India
3 Department of Pharmacy, Bengal Institute of Pharmaceutical Sciences, Kalyani, Nadia, West Bengal, India

Date of Web Publication30-Jan-2015

Correspondence Address:
Sudhan Debnath
Department of Chemistry, MBB College, West Tripura, Agartala, Tripura - 799 004
India
Login to access the Email id

Source of Support: Financial Support from DBT, New Delhi., Conflict of Interest: None


DOI: 10.4103/2231-4040.150366

Rights and Permissions
  Abstract 

Aldose reductase (AR) plays an important role in the development of several long-term diabetic complications. Inhibition of AR activities is a strategy for controlling complications arising from chronic diabetes. Several AR inhibitors have been reported in the literature. Flavonoid type compounds are shown to have significant AR inhibition. The objective of this study was to perform a computational work to get an idea about structural insight of flavonoid type compounds for developing as well as for searching new flavonoid based AR inhibitors. The data-set comprising 68 flavones along with their pIC 50 values ranging from 0.44 to 4.59 have been collected from literature. Structure of all the flavonoids were drawn in Chembiodraw Ultra 11.0, converted into corresponding three-dimensional structure, saved as mole file and then imported to maestro project table. Imported ligands were prepared using LigPrep option of maestro 9.6 version. Three-dimensional quantitative structure-activity relationships and docking studies were performed with appropriate options of maestro 9.6 version installed in HP Z820 workstation with CentOS 6.3 (Linux). A model with partial least squares factor 5, standard deviation 0.2482, R 2 = 0.9502 and variance ratio of regression 122 has been found as the best statistical model.

Keywords: Aldose reductase, flavonoids, three-dimensional quantitative structure-activity relationships, extra precision glide docking, designed inhibitor


How to cite this article:
De UC, Debnath T, Sen D, Debnath S. Three-dimensional quantitative structure-activity relationships and docking studies of some structurally diverse flavonoids and design of new aldose reductase inhibitors. J Adv Pharm Technol Res 2015;6:13-8

How to cite this URL:
De UC, Debnath T, Sen D, Debnath S. Three-dimensional quantitative structure-activity relationships and docking studies of some structurally diverse flavonoids and design of new aldose reductase inhibitors. J Adv Pharm Technol Res [serial online] 2015 [cited 2019 Nov 20];6:13-8. Available from: http://www.japtr.org/text.asp?2015/6/1/13/150366


  Introduction Top


Diabetes mellitus has become a common disease in many countries around the world, currently 387 million people worldwide suffering from this disease and expected to affect about 592 million peoples by 2035 [Figure 1]. [1],[2] Under normal glycemic conditions, the majority of glucose is metabolized through glycolysis and only a small fraction of glucose is metabolized through the polyol pathway [3],[4] . However, in diabetic condition glucose concentrations are elevated in tissues leading to increased flux of glucose through the polyol pathway where aldose reductase (AR) plays an important role as a key enzyme. AR being present in the lens, retina, Schwann cells of peripheral nerves may lead to blindness associated with increased risk for painful neuropathy, heart disease and kidney failure due to increased sorbitol flux and NADH/NAD + ratio associated with decreased NADPH/NADP + ratio [Figure 1], which plays an important role on other enzymes. Thus, AR is primarily important in the development of degenerative long-term diabetic complications such as cataract, neuropathy, retinopathy, and nephropathy. [5],[6],[7],[8],[9] Inhibition of AR activities is, therefore, a useful strategy for prevention and treatment of complications arising from chronic diabetes. Several classes of compounds have been reported in the literature [10],[11],[12],[13],[14],[15],[16] as AR inhibitors such as flavonoids, isoflavonoids, coumarins, stilbenes, rosmarinic acid derivatives, thiazolidinediones, [17] etc., Among these, flavonoids and 2, 4-thiazolidinedione derivatives have been found to be more potent. These facts prompted us to perform a computational work with a view to get an idea regarding the structural insight of flavonoids to develop and searching for new selective and effective flavonoid based AR inhibitors.
Figure 1: Glucose metabolism in glycolysis and polyol pathway

Click here to view



  Materials and methods Top


Data-set

The data-set is comprising 68 flavones along with their pIC 50 values ranging from 0.44 to 4.59 have been collected from literature. [18] It is hypothesized that active compounds share all or most of the required features for binding with the active site of the target molecule while the inactive compounds experience steric hindrance and other disfavored interactions. In view of these, the whole data-set of 68 molecules were clustered, and 54 flavonoid molecules were selected for the computational study. The selected subset randomly divided into a training set of 38 compounds (~70%) and a test set of 16 compounds [Table 1].
Table 1: List of the flavonoid derivatives with their observed and predicted biological activities

Click here to view


Ligand preparation

Structure of all the flavonoids were drawn in Chembiodraw Ultra 11.0 (Cambridge Soft), converted into corresponding three-dimensional structure, saved as mole file and then imported to maestro project table. Imported ligands were prepared using LigPrep option of Schrodinger maestro 9.6 version [19] where ligand minimization was done applying orthogonal partial least squares_2005 force field. Maximum of 100 conformers was generated by 100 steps preprocess and 50 steps postprocess minimization. Conformers were filtered for 11.4 kcal/mol energy and 2Ǻ atom deviation with chiralities form three-dimensional structure and retaining original states of ionization.

Building three-dimensional-quantitative structure- activity relationships model

Ligand data-set was aligned [Figure 2] by atom type macro model option under shape screen. In Build quantitative structure-activity relationships (QSAR) option randomly selected 70% molecules were kept as training set and model was generated with 1Ε grid spacing and partial least squares (PLS) factor 7. A model with PLS factor 5, 0.2482 as a standard deviation, 0.9502 as R 2 and 122 as F (variance ratio for the regression) has been found as the best statistical model.
Figure 2: Structural alignments of all molecules in the data-set

Click here to view


Docking study

The docking studies were carried out with human AR (PDB ID: 3M0I, resolution 1.07 Ε7 retrieved from RCSB protein data bank. The root mean square deviation (0.20 Ε0 of superimposition between docking and co-crystalline pose of co-ligand was calculated. This AR was mutated with Val-113, refined by protein preparation wizard of Schrodinger Maestro version 9.6. A grid box with maximum 12Ε edge lengths was generated centering the co-ligand and retaining other default settings. We have virtually replaced Val-113 (mutated residue) of 3M0I by Thr-113 (normal residue) and compared the ligand interactions with the receptor. Extra precision (XP) glide docking [19],[20],[21] was performed with both AR separately for the entire data-set [Table 1] enabling the XP descriptor option and other default values in glide module [22] of maestro 9.6 were retained. The co-crystal ligand structure of 3M0I was superimposed [Figure 3] on the docking poses of the same with mutated (Val-113) and nonmutated AR.
Figure 3: Superimposition of the coligand of 3M0I with it's docked poses computed by extra precision glide (a) without mutation (b) with mutation

Click here to view



  Results and Discussion Top


Atom-based three-dimensional QSAR models were generated for a series of 54 structurally diverse flavonoid derivatives with pIC 50 ranging from 3.96 to 6.59 [Table 1]. A statistically significant three-dimensional -QSAR model was generated, and over-fitting was avoided with the results of different PLS factor and relevant statistical parameters are summarized in [Table 2].
Table 2: Summary of statistical parameters of the more significant models with different PLS factors

Click here to view


Three-dimensional quantitative structure-activity relationships analysis

In this study, AR inhibitory activity of the selected ligands were tried to visualize by a developed QSAR model. Different statistical parameters of QSAR models so far generated are presented in [Table 2]. Results of different structural features viz. hydrogen bonding, hydrophobic interaction and electron withdrawing effect were applied to the most active compound 1 and the least active compound 54 and are shown [Figure 4],[Figure 5] and [Figure 6]. Inspection of three-dimensional QSAR visualization [Figure 4] generated by PHASE [23] revealed that 5, 7, 3, 4 positions as indicated by blue cubes are favorable for hydrogen bond donor groups in flavonoid skeleton and the orange cubes in the vicinity of 6, 3 and 2' positions are unfavorable. These observations are corroborated by the activity of compound 1 (highest active) in which the hydrogen bond donor groups are available at position 5, 7, 3', 4' while in least active compound 54, the 7-OH and 3'-OH groups are replaced by-OCH 3 .
Figure 4: Hydrogen bonding interactions: Blue implies positive and orange implies negative effect; 1: Highest active, 54: Lowest active

Click here to view
Figure 5: Hydrophobic interaction: Green implies positive and purple implies negative effect; 1: Highest active, 54: Lowest active

Click here to view
Figure 6: Electron withdrawing effect: Green implies positive and orange implies negative effect; 1: Highest active, 54: Lowest active

Click here to view


Similarly presence of hydrophobic groups at 3, 6 and 8 positions as indicated by green color [Figure 5] are favorable for better activity while presence of similar group at 7, 3' and 4' indicated by purple color are detrimental for activity. It is seen that in compound 1,-OCH 3 group occupied 3 and 6 positions and had favorable contribution to its enhanced activity. In compound 54 presence of - OH group at position 6,-OCH 3 group at 7 and absence of the hydrophobic group at 3 positions greatly reduced the activity.

Three-dimensional quantitative structure-activity relationships visualization [Figure 6] indicates that electron withdrawing group at 3, 3 and 4n positions that is, in the vicinity of green region are favorable but position 8 should be devoid of such groups for better activity. This is further corroborated by the activity of compound 1 in which electron withdrawing oxygen atom is present at position 3', 4' (as -OH), and 3 (as -OCH 3 ). On the other hand in compound 54 one of the favorable positions (3') is occupied by weak electron withdrawing -OCH 3 group (-OH>-OCH 3 due to + I effect of -CH 3 ) and least favorable '8' position is occupied by - OCH 3 group represented by orange color in [Figure 6].

The experimental and predicted pIC 50 values of the training and test sets are given in [Table 1] while regression graph of experimental versus predicted pIC 50 values of the training and test sets are illustrated [Figure 7].
Figure 7: Correlation diagram of observed and predicted activity (a: Training set b: Test set)

Click here to view


Docking analysis

The docking results (XP Glide docking parameters) of the selected flavonoids with mutated, as well as normal AR, are summarized in [Table 3].
Table 3: XP glide docking parameters of some selected (pIC50>7) flavonoids of the dataset

Click here to view


Ligand-receptor interactions in mutated as well as in normal AR are found similar. This may be due to because mutation of the AR (Thr-113 ↔X) is not in the active site. The highest active compound 1 and the co-ligand show similar interactions with residues viz. Leu-300 [Figure 8]a. This result also corroborated with the three-dimensional-QSAR model that developed in Phase module. The docking study also indicates that the 5-OH act as favorable hydrogen bond acceptor in flavones. This feature is supported by the activity (within parenthesis) of 1 (7.59), 7 (7.24), 17 (6.69), 26 (6.52), 32 (6.02) that interacts with Leu-300 of protein (AR) that is, interaction of 5-OH with Leu-300 is crucial for potential activity. Hydrogen bond donor interaction with Val-47 and 5-OH of flavones were found in 31 (6.09), 33 (6.07), 34 (5.92), 35 (5.92), 37 (6.54), 42 (5.09), 43 (5.08) and attributed the observed activity of these compounds. The ðhe staking interaction between Trp-111 and R1 of flavones in compound 1 (7.59), 7 (7.24), 17 (6.69), 32 (6.02) are also plays an important role to have good activity. The interaction of Trp-111, an amino acid residues in the active site of AR with Ri [Figure 9] of flavones are also found in maximum number (41) of compounds such as 2 (7.49), 3 (7.44), 4 (7.47), 5 (7.41), 6 (7.35), 8 (7.19), 9 (7.13), 10 (7.11), 11 (7.04), 12 (6.92), 13 (6.85), 14 (6.79), 16 (6.77), 18 (6.92), 19 (6.64), 20 (6.62), 22 (6.57), 23 (6.55), 24 (6.55), 26 (6.52), 27 (6.52), 28 (6.46), 29 (6.39), 30 (6.27), 31 (6.09), 33 (6.07), 34 (5.92), 35 (5.92), 36 (5.85), 37 (6.54), 38 (5.35), 39 (5.20), 40 (5.17), 41 (5.14), 42 (5.09), 43 (5.08), 44 (5.05), 49 (4.74), 51 (4.67), 52 (4.42), 54 (3.96). Out of these 41 interacting flavones twenty-seven have activity > 6, ten have activity in the range of 5-6 and only four have activity < 5. Therefore, interaction between Trp-111 and R2 of flavones play an important role in the activity against AR. Another ðno staking interaction with Trp-20 and Rn of flavones plays a significant role in ligand interaction and are also observed in 1 (7.59), 7 (7.24), 15 (6.79), 17 (6.69), 21 (6.60), 25 (6.52), 28 (6.46), 32 (6.02), 53 (4.44). Considering altogether the favorable and unfavorable features we have designed large number of AR inhibitors and performed similar docking studies. Two-dimensional views of ligand-receptor interaction diagram of co-ligand, highest active ligand and two best designed new molecules [Figure 8]a and b. Since we are practically not in a position to study the AR inhibition and hence we did not attempt to synthesize the highly active designed molecule. However, it can be tried in the future.
Figure 8: (a) Two-dimensional ligand interaction diagram of coligand and highest active ligand (b) Two-dimensional ligand interaction diagram of two best designed new molecules

Click here to view
Figure 9: Skeletal structure of flavone

Click here to view



  Conclusion Top


This atom based three-dimensional-QSAR model and docking studies could help in understanding the relationship between structure and activity of flavonoid type compounds as AR inhibitors and gives us various options to design a novel and potent inhibitors for AR. However, other conditions like transport properties and steric factors are also needed to be considered while further new molecules are designed.


  Acknowledgment Top


The authors are thankful to Department of Biotechnology, New Delhi, Government of India for providing financial support [F.No.BT/327/NE/TBP/2012] for this research work.

 
  References Top

1.
Zimmet P, Alberti KG, Shaw J. Global and societal implications of the diabetes epidemic. Nature 2001;414:782-7.  Back to cited text no. 1
    
2.
IDF Diabetes Atlas Sixth Edition Poster Update 2014. Available from: http://www.idf.org/diabetesatlas/update-2014.  Back to cited text no. 2
    
3.
Gabbay KH, Merola LO, Field RA. Sorbitol pathway: Presence in nerve and cord with substrate accumulation in diabetes. Science 1966;151:209-10.  Back to cited text no. 3
[PUBMED]    
4.
Lindstad RI, McKinley-McKee JS. Methylglyoxal and the polyol pathway. Three-carbon compounds are substrates for sheep liver sorbitol dehydrogenase. FEBS Lett 1993;330:31-5.  Back to cited text no. 4
    
5.
Gabbay KH, O′Sullivan JB. The sorbitol pathway. Enzyme localization and content in normal and diabetic nerve and cord. Diabetes 1968;17:239-43.  Back to cited text no. 5
[PUBMED]    
6.
Robison WG Jr, Kador PF, Kinoshita JH. Retinal capillaries: Basement membrane thickening by galactosemia prevented with aldose reductase inhibitor. Science 1983;221:1177-9.  Back to cited text no. 6
[PUBMED]    
7.
Lipinski CA, Hutson NJ. Aldose reductase inhibitors as a new approach to the treatment of diabetic complications. Annu Rep Med Chem 1984;19:169-77.  Back to cited text no. 7
    
8.
Kinoshita JH, Dvornik D, Kraml M, Gabbay KH. The effect of an aldose reductase inhibitor on the galactose-exposed rabbit lens. Biochim Biophys Acta 1968;158:472-5.  Back to cited text no. 8
[PUBMED]    
9.
Cogan DG, Kinoshita JH, Kador PF, Robison G, Datilis MB, Cobo LM, et al. NIH conference. Aldose reductase and complications of diabetes. Ann Intern Med 1984;101:82-91.  Back to cited text no. 9
[PUBMED]    
10.
Kador PF, Robison WG Jr, Kinoshita JH. The pharmacology of aldose reductase inhibitors. Annu Rev Pharmacol Toxicol 1985;25:691-714.  Back to cited text no. 10
[PUBMED]    
11.
Tomlinson DR, Willars GB, Carrington AL. Aldose reductase inhibitors and diabetic complications. Pharmacol Ther 1992;54:151-94.  Back to cited text no. 11
    
12.
Kador PF, Kinoshita JH, Sharpless NE. Aldose reductase inhibitors: A potential new class of agents for the pharmacological control of certain diabetic complications. J Med Chem 1985;28:841-9.  Back to cited text no. 12
[PUBMED]    
13.
Sarges R, Oates PJ. Aldose reductase inhibitors: Recent developments. Prog Drug Res 1993;40:99-161.  Back to cited text no. 13
    
14.
Humber LG. The medicinal chemistry of aldose reductase inhibitors. Prog Med Chem 1987;24:299-343.  Back to cited text no. 14
[PUBMED]    
15.
Liu HY, Liu SS, Qin LT, Mo LY. CoMFA and CoMSIA analysis of 2, 4-thiazolidinediones derivatives as aldose reductase inhibitors. J Mol Model 2009;15:837-45.  Back to cited text no. 15
    
16.
Miyamoto S. Recent advances in aldose reductase inhibitors: Potential agents for the treatment of diabetic complications. Expert Opin Ther Pat 2002;12:621-31.  Back to cited text no. 16
    
17.
Schemmel KE, Padiyara RS, D′Souza JJ. Aldose reductase inhibitors in the treatment of diabetic peripheral neuropathy: A review. J Diabetes Complications 2010;24:354-60.  Back to cited text no. 17
    
18.
Caballero J. 3D-QSAR (CoMFA and CoMSIA) and pharmacophore (GALAHAD) studies on the differential inhibition of aldose reductase by flavonoid compounds. J Mol Graph Model 2010;29:363-71.  Back to cited text no. 18
    
19.
Friesner RA, Banks JL, Murphy RB, Halgren TA, Klicic JJ, Mainz DT, et al. Glide: A new approach for rapid, accurate docking and scoring 1. Method and assessment of docking accuracy. J Med Chem 2004;47:1739-49.  Back to cited text no. 19
    
20.
Halgren TA, Murphy RB, Friesner RA, Beard HS, Frye LL, Pollard WT, et al. Glide: A new approach for rapid, accurate docking and scoring 2. Enrichment factors in database screening. J Med Chem 2004;47:1750-9.  Back to cited text no. 20
    
21.
Friesner RA, Murphy RB, Repasky MP, Frye LL, Greenwood JR, Halgren TA, et al. Extra precision glide: Docking and scoring incorporating a model of hydrophobic enclosure for protein-ligand complexes. J Med Chem 2006;49:6177-96.  Back to cited text no. 21
    
22.
Glide. Ver. 6.0. New York: Schrodinger, LLC; 2013.  Back to cited text no. 22
    
23.
Maestro. Ver. 9.6. New York: Schrodinger, LLC; 2013.  Back to cited text no. 23
    


    Figures

  [Figure 1], [Figure 2], [Figure 3], [Figure 4], [Figure 5], [Figure 6], [Figure 7], [Figure 8], [Figure 9]
 
 
    Tables

  [Table 1], [Table 2], [Table 3]


This article has been cited by
1 In vitro Anti-Cancer Effect of Wellness-Compound (Ochnaflavone)
Jae-Sook Lee,Hwa-Jung Choi,Myung-Ju Kim,Jang-Soon Park
Journal of Digital Convergence. 2015; 13(5): 337
[Pubmed] | [DOI]



 

Top
 
 
  Search
 
    Similar in PUBMED
   Search Pubmed for
   Search in Google Scholar for
 Related articles
    Access Statistics
    Email Alert *
    Add to My List *
* Registration required (free)  

 
  In this article
   Abstract
  Introduction
   Materials and me...
   Results and Disc...
  Conclusion
  Acknowledgment
   References
   Article Figures
   Article Tables

 Article Access Statistics
    Viewed2135    
    Printed114    
    Emailed0    
    PDF Downloaded347    
    Comments [Add]    
    Cited by others 1    

Recommend this journal