|Year : 2015 | Volume
| 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
Utpal Chandra De1, Tanusree Debnath2, Debanjan Sen3, Sudhan Debnath2
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 Publication||30-Jan-2015|
Department of Chemistry, MBB College, West Tripura, Agartala, Tripura - 799 004
Source of Support: Financial Support from DBT, New Delhi., Conflict of Interest: None
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|| |
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]. , 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 , . 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. ,,,, 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 ,,,,,, as AR inhibitors such as flavonoids, isoflavonoids, coumarins, stilbenes, rosmarinic acid derivatives, thiazolidinediones,  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.
| Materials and methods|| |
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.  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
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  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.
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 ,, 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  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|| |
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  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
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
| Conclusion|| |
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|| |
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|| |
Zimmet P, Alberti KG, Shaw J. Global and societal implications of the diabetes epidemic. Nature 2001;414:782-7.
IDF Diabetes Atlas Sixth Edition Poster Update 2014. Available from: http://www.idf.org/diabetesatlas/update-2014.
Gabbay KH, Merola LO, Field RA. Sorbitol pathway: Presence in nerve and cord with substrate accumulation in diabetes. Science 1966;151:209-10.
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.
Gabbay KH, O′Sullivan JB. The sorbitol pathway. Enzyme localization and content in normal and diabetic nerve and cord. Diabetes 1968;17:239-43.
Robison WG Jr, Kador PF, Kinoshita JH. Retinal capillaries: Basement membrane thickening by galactosemia prevented with aldose reductase inhibitor. Science 1983;221:1177-9.
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.
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.
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.
Kador PF, Robison WG Jr, Kinoshita JH. The pharmacology of aldose reductase inhibitors. Annu Rev Pharmacol Toxicol 1985;25:691-714.
Tomlinson DR, Willars GB, Carrington AL. Aldose reductase inhibitors and diabetic complications. Pharmacol Ther 1992;54:151-94.
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.
Sarges R, Oates PJ. Aldose reductase inhibitors: Recent developments. Prog Drug Res 1993;40:99-161.
Humber LG. The medicinal chemistry of aldose reductase inhibitors. Prog Med Chem 1987;24:299-343.
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.
Miyamoto S. Recent advances in aldose reductase inhibitors: Potential agents for the treatment of diabetic complications. Expert Opin Ther Pat 2002;12:621-31.
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.
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.
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.
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.
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.
Glide. Ver. 6.0. New York: Schrodinger, LLC; 2013.
Maestro. Ver. 9.6. New York: Schrodinger, LLC; 2013.
[Figure 1], [Figure 2], [Figure 3], [Figure 4], [Figure 5], [Figure 6], [Figure 7], [Figure 8], [Figure 9]
[Table 1], [Table 2], [Table 3]