PHARMACOPHORE MAPPING AND VIRTUAL SCREENING ON HUMAN ALPHA AMYLASE INHIBITORS
HTML Full TextPHARMACOPHORE MAPPING AND VIRTUAL SCREENING ON HUMAN ALPHA AMYLASE INHIBITORS
- Santhosh
Department of Bioinformatics, School of Chemical and Biotechnology, SASTRA University, Thirumalaisamudram, Thanjavur - 613401, India.
ABSTRACT: Diabetes is a non communicable chronic metabolic disease characterized by high blood glucose level in the body. Targeting alpha amylase to control diabetes has long been studied. This study focused on developing pharmacophore model for set of inhibitor that possess anti amylase activity based on previous studies. Developing a pharmacophore model for these compound yields a hypothesis of ADDRR. Furthermore 3D QSAR model is carried out with good correlation activity for training set as .98 and test set as .92. Database screening for the developed model list out top compounds. High throughput virtual screening is performed for this top compound set. Two ligands quercetin 3-glucosyl rhamnosyl galactoside and rhodionin is found to have good binding affinity over alpha amylase
Keywords: |
Diabetes, Alpha amylase, Pharmacophore, Virtual screening
INTRODUCTION: Non communicable chronic diseases have become the major causes of mortality and disease burden worldwide. According to the latest estimates by the International Diabetes Federation (IDF), South East-Asia (SEA) Region account for more than 72 million adults with diabetes in 2013 and is expected to exceed 123 million in 2035 out of which nearly 95% of people have type 2 diabetes. 1 Thus there is a need for the development of a utile approach to the diabetes prevention as well as understanding its underlying pathophysiology. 2 The improved life style involving a healthy diet, reduced weight and increased physical activity have been observed to have long lasting effects in the prevention of type 2diabetes. 3
The most practical approach involves targeting the carbohydrate metabolism by the inhibition of enzyme alpha amylase and alpha glucosidase. 4 Alpha amylases catalyze the hydrolysis of a-(1,4) glycosidic linkages in starch resulting in the formation di and tri saccharides which are further broken into monomeric glucose by α-glucosidase. 5 Currently several synthetic drugs are available such as acarbose, voglibiose and miglitol which act as inhibitors of digestive enzymes in type 2 diabetic patients. 6 However various side effects, such as flatulence, abdominal distension, diarrhea, and several reports of hepatitis have been observed with the acarbose use, which was not observed after the discontinuation of medication. 7, 8 Therefore attempts have been aimed towards finding natural and comparatively safer α-amylase and α-glucosidase inhibitors.
The present study was focused for the identification of novel inhibitors for the α-amylase from medicinal plants. Pharmacophore analysis was done to identify the active functional groups responsible for the binding of inhibitor in the active site of amylase. QSAR studies and regression coefficient analysis were performed to identify the reliability of model. Finally database screening was done against 1,60,000 natural compounds from database like pubchem, ZINC to identify best hits. Molecular docking was done for the top hits to calculate its binding affinity over α-amylase.
MATERIALS AND METHODS:
Data set:
A set of 30 alpha amylase inhibitors with known IC50 are collected based on the previous reports by 8, 9, 10 11, 12 13 14, 15 shown in Table 1. The 2D structure of 30 compounds are shown in Fig 1. IC50 values were converted into PIC50 which indicates the negative logarithm of IC50 values. Pharmacophore mapping were made using PHASE (2009) and compounds are partitioned into active and inactive set.
FIG.1: STRUCTURE OF COMPOUNDS USED FOR PHARMACOPHORE MODELLING.
TABLE 1: A SET OF 30 COMPOUNDS PARTITIONED INTO ACTIVE AND INACTIVE BASED ON PIC50 FOR PHARMACOPHORE MODELING
Compound No | IC50(μM) | PIC50 | Compound Name |
1 | 3.420 | 5.4659 | Myricetin |
2 | 3.301 | 5.4813 | Quercetin |
3 | 2.275 | 5.643 | Kaempferol |
4 | 3.221 | 5.492 | Fisetin |
5 | 3.443 | 5.463 | Luteolin |
6 | 3.107 | 5.571 | Cyanidin-3-glucoside |
7 | 3.154 | 5.501 | Apigenin |
8 | 3.187 | 5.496 | Baicalein |
9 | 3.096 | 5.509 | Naringenin |
10 | 3.517 | 5.454 | Catechin |
11 | 3.045 | 5.516 | Hesperetin |
12 |
3.240 |
5.489 |
Daidzein |
13 | 3.110 | 5.507 | Genistein |
14 | 3.517 | 5.454 | Epicatechin |
15 | 3.015 | 5.521 | Trans-chalcone |
16 | 3.612 | 5.442 | Cyanidin-3-rutinoside |
17 | 3.665 | 5.436 | Beta-sitosterol |
18 | 3.853 | 5.414 | Lupeol |
19 | 2.229 | 5.652 | Avicularin |
20 | 4.096 | 5.388 | Cyanidin |
21 | 2.251 | 5.648 | Guaijaverin |
22 | 2.214 | 5.659 | Hyperin |
23 | 3.190 | 5.496 | Cyanidin 3-galactoside |
24 | 2.752 | 5.560 | Cinnamyl cinnamate |
25 | 2.853 | 5.545 | Chlorogenic acid |
26 | 2.896 | 5.538 | Chrysoeriol |
27 | 3.420 | 5.466 | Caffeic acid |
28 | 4.602 | 5.337 | Bisdemethoxycurcumin |
29 | 1.576 | 5.8202 | Quinic acid |
30 | 3.275 | 5.845 | Diphlorethohydroxycarmalol |
Ligand preparation:
Ligand preparations were done using which lead to the development of energy minimized structures. It results in the addition of hydrogen and Neutralization of charged acidic or basic groups 16. For conformer generation, Confgen module was used 17. The output of the Confgen may have high energy compounds which are processed to remove those conformers before going into the next step. Then total number of conformation per rotatable bond was kept 100 using OPLS_2005 Force field. Distance dependent dielectrics were used for solvation treatment. The RMSD cutoff was set to 1Å to remove possible duplicate entries.
Pharmacophore Hypothesis Generation:
Site points are generated based on six pharmacophore features such as hydrogen bond donors, hydrogen bond acceptors, aromatic ring, hydrophobic groups, cations and anions. Site points are the active groups in the ligand that binds to receptor atoms. Common pharmacophore were identified using tree based partitioning technique. Common pharmacophore hypotheses (CPH) are calculated from highly active set of conformations which were examined using scoring function to get best alignment with RMSD box size of 1.2 Å.
Generation of QSAR models
Quantitative structure activity relationship (QSAR) was used to predict the biological activity of the molecule based on its structure. Models were developed by partitioning the compound into training set and tests set. All hypotheses were generated and scoring was done to build pharmacophore-based QSAR models while setting the grid spacing value at 1.0Å. QSAR model are represented as cubes. Cubes are colored according to their coefficient values. Positive coefficient represented in blue color indicates increase in the pharmacophore activity. Negative coefficients which indicate low activity were represented in red color.
Database screening:
Virtual screening of databases was used to identify novel compounds with highest activity. Screening were done 1,60,000 natural compounds from PubChem 18 and ZINC database 19 to obtain best ligand that matched the selected hypothesis.
Docking studies
Molecular docking studies were performed for ligands identified with highest activity to determine the binding pose and affinity in the active site of α-amylase. Docking was done using Glide 5.5. Crystal structure of α-amylase was obtained from PDB with PDB ID of 1HNY 20 which was used for protein preparation using Protein preparation wizard a module of Glide 5.5. Site map was used to determine the binding site in the protein and ligands were docked. Ligands were ranked based on their binding energy and interaction pattern was analyzed for docked structures.
RESULTS AND DISCUSSION:
Pharmacophore models were generated keeping a terminal box size of 1Å with different active molecules using a tree based partition algorithm. To find the common pharmacophore hypothesis, the dataset was divided into active and inactive set depending on the observed activity. Activities of the compounds are shown in Table 1. The top model was found to be associated with the five-point hypotheses ADDRR, which consists of one hydrogen bond acceptor (A), two hydrogen bond donors (D) and two aromatic rings (R). This is denoted as A7, D20, D17, R21 and R22. The pharmacophoric sites obtained are shown in Fig 2.
FIG. 1: COMMON PHARMACOPHORIC DESCRIPTORS OBTAINED FROM PHASE. PINK: HYDROGEN-BOND ACCEPTOR. BLUE: HYDROGEN-BOND DONOR. RED: AROMATIC RINGS.
These hypothesis was subjected to PLS factor of 3 for QSAR model development. The best hypotheses fitted with an active ligand are shown in Fig 3.
FIG. 2: PHARMACOPHORE MODEL FITTED WITH ACTIVE LIGAND HYPEROSIDE.
This active ligand was found to be hyperoside. Angle and distance geometry of pharmacophoric sites are shown in Fig 4.
FIG. 3: GEOMETRY OF PHARMACOPHORE HYPOTHESIS a) DISTANCE GEOMETRY b) ANGLE GEOMETRY.
For the QSAR models generation, inactive molecules in the dataset were then aligned, based on matching with at least three pharmacophore features. These dataset were randomly divided into training set of 19 compounds and the test set of 10 compounds for QSAR study. The best pharmacophore model resulted are ADDRR. Pictorial representation of contours generated from QSAR fitted with active and inactive ligand are shown in Fig 5.
FIG.4: PICTORIAL REPRESENTATION OF CONTOURS GENERATED USING QSAR MODEL WITH BLUE COLOR REPRESENTING POSITIVE CORRELATION AND RED COLOR REPRESENTING NEGATIVE CORRELATION. a) QSAR MODEL FITTED WITH ACTIVE LIGAND HYPEROSIDE b) QSAR MODEL FITTED WITH INACTIVE LIGAND BISDEMETHOXYCURCUMIN
Regression plots of the training and test set showed a good regression values of r2=.98 and r2=.92 respectively. The regression plot of the training set and test set are shown in Fig 6.
FIG. 5: SCATTER PLOT OBTAINED FOR QSAR MODEL THAT ARE GENERATED a) TRAINING DATA b) TEST DATA.
Virtual screening was done for 1, 60,000 compounds. Since the flexibility of ligand is limited in case of virtual screening docking was performed using Glide. Quercetin 3-glucosyl rhamnosyl galactoside (QGRG) and rhodionin shown in Fig 7 were identified as best ligand from the database.
FIG.6: TOP TWO LIGAND OBTAINED FROM VIRTUAL SCREENING FITTED WITH PHARMACOPHORE HYPOTHESIS. a) QUERCETIN 3-GLUCOSYL RHAMNOSYL GALACTOSIDE b) RHODIONIN.
Interaction of ligand with protein are shown in ligand interaction diagram Fig 8 represent many favourable interactions.
QGRG has been reported to be present in green tea by 21 whereas rhodionin is present in rhodiola species 22. Both green tea and rhodiola are reported to have anti-amylase activity 23, 24, 25, 26, 27, 28, 29 Inclusion of above two plants in normal culinary activity can help us to prevent diabetes.
FIG. 7: LIGAND INTERACTION DIAGRAM FOR TOP TWO LIGAND OBTAINED AFTER DOCKING WITH ALPHA-AMYLASE a) QUERCETIN 3-GLUCOSYL RHAMNOSYL GALACTOSIDE b) RHODIONIN.
CONCLUSION: The purpose of this study was to construct a pharmacophore model for alpha amylase inhibitor using PHASE module from schrodinger. A five point pharmacophore hypothesis was developed ADDRR which constitute one hydrogen bond acceptor, two hydrogen bond donor and two aromatic rings. 3D QSAR study was carried out on the series of phytochemicals as an alpha amylase inhibitor by generating volume occupied maps, which demonstrated that the activity may be increased by substituting the groups in the pharmacophore features.
Accuracy of QSAR model was evaluated based on regression plot obtained for training set as .98 and test set as .92. Database screening against natural compounds also provide top hits with good binding affinity against alpha amylase. These compounds were also found to be present in medicinal plants and hence inclusion of these plants in normal diet can help prevent diabetes.
ACKNOWLEDGEMENT: I would like to thank SASTRA University for providing the necessary facilities to carry out this research project.
CONFLICT OF INTEREST: There exists no conflict of interest.
REFERENCES:
- Xu Y, Wang L, He J, Bi Y, Li M, Wang T, Wang L, Jiang Y, Dai M, Lu J, Xu M, Li Y, Hu N, Li J, Mi S, Chen CS, Li G, Mu Y, Zhao J, Kong L, Chen J, Lai S, Wang W, Zhao W, Ning G: Prevalence and Control of Diabetes in Chinese Adults. JAMA 2013; 310(9):948-958.
- Alberti KG, Zimmet P; Epidemiology:Global burden of disease-where does diabetes mellitus fit in?. Nat Rev Endocrinol 2013; 9:258-60.
- Amarowicz R, Pegg RB, Rahimi-Moghaddam P, Barl B and Weil JA: Free-radical scavenging capacity and antioxidant activity of selected plant species from the Canadian prairies. Food Chem 2004; 84:551-562.
- Shim YJ, Doo HK, Ahn SY, Kim YS, Seong JK, Park IS and Min BH: Inhibitory effect of aqueous extract from the gal of Rhuz chinesis on alpha-glucosidase activities and postprandial blood J Ethnopharmacol 2003; 85: 283–287.
- Qin X, Ren L, Yang X, Bai F, Wang L, Geng P, Bai G and Shen Y: Structures of human pancreatic α-amylase in complex with acarviostatins: Implications for drug design against type II diabetes. J Struct Biol 2011; 174:196-202.
- Baker JR, Zyzak DV, Thorpe SR and Baynes JW: Chemistry of the fructosamine assay: D-glucosone is the product of oxidation of Amadori compounds during the fructosamine assay. Clin Chem 1994; 40:1950-1955.
- Baker JR, Zyzak DV, Thorpe SR and Baynes JW: Mechanism of Fructosamine Assay: Evidence against Role of Superoxide as Intermediate in Nitroblue Tetrazolium Reduction. Clin Chem 1993; 39:2460-2465.
- Bendich A, Machlin LJ, Scandurra O, Burton GW and Wayner DDM: The antioxidant role of vitamin C. Free Radical Bio Med 1986; 2:419-444.
- Tadera K, Minami Y, Takamatsu K and Matsuoka T: Inhibition of α-glucosidase and α-amylase by flavonoids. J Nutr Sci Vitaminol 2006; 52:149-153.
- Najafian M, Ebrahim-Habibi A, Hazareh N, Yaghmaei P, Parivar K and Larijani B: Trans-chalcone: a novel small molecule inhibitor of mammalian α–amylase. Mol Biol Rep 2011; 38:1617-1620.
- Yu L, Zhao M, Yang B, Zhao Q and Jiyang Y: α-glucosidase and α-amylase inhibitory activities and antioxidant activities of a compound isolated from Terminalia sericea. Food chem 2007; 104:176-181.
- Bhandari MR, Anurakkun NJ, Hong G and Kawabata J: α-Glucosidase and α-amylase inhibitory activities of Nepalese medicinal herb Pakhanbhed. Food chem 2008; 106:247-252.
- Wang H, Du YJ and Song HC: α-Glucosidase and α-amylase inhibitory activities of guava leaves, Food chem 2010; 123:6-13.
- Funke I and Melzig M: Traditionally used plants in diabetes therapy - phytotherapeutics as inhibitors of α-amylase activity. Braz J Pharmacog 2006; 16:1-6.
- Rohn S, Rawel HM and Kroll J: Inhibitory effects of plant phenols on the activity of selected enzymes. J Agric Food Chem 2002; 50:3566-3571.
- Akkarachiyasit S, Charoenlertkul P, Yibchok-anun S and Adisakwattana S: Inhibitory Activities of Cyanidin and Its Glycosides and Synergistic Effect with Acarbose against Intestinal α-Glucosidase and Pancreatic α-Amylase. Int J Mol Sci 2010; 11:3387-3396.
- Sales PMD, Souza PMD, Simeoni LA, Magalhães PDO and Silveira D: α-Amylase Inhibitors: A Review of Raw Material and Isolated Compounds from Plant Source. Journal of Pharmacology and Pharmaceutical Sciences 2012; 15:141-183.
- Phase, version 3.1, Schrödinger, LLC, New York, NY, 2009.
- LigPrep, version 2.3, Schrödinger, LLC, New York, NY, 2009.
- ConfGen, version 2.1, Schrödinger, LLC, New York, NY, 2009.
- Bolton E, Wang Y, Thiessen PA and Bryant SH: PubChem: Integrated Platform of Small Molecules and Biological Activities. Chapter 12 INAnnual Reports in Computational Chemistry 2008; Volume 4, American Chemical Society.
- Irwin JJ and Shoichet BK: ZINC - A Free Database of Commercially Available Compounds for Virtual Screening. J Chem Inf Model 2005; 45(1):177-182.
- Brayer GD, Luo Y and Withers SG: The structure of human pancreatic alpha-amylase at 1.8 A resolution and comparisons with related enzymes. Protein Sci 1995; 4:1730-1742.
- Protein Preparation Wizard, version 2.0, Schrödinger, LLC, New York, NY, 2009.
- Glide, version 5.5, Schrödinger, Inc., New York, NY, 2009.
- Kim Y, Goodner KL, Park JD, Choi J and Talcott ST: Changes in antioxidant phytochemicals and volatile composition of Camellia sinensis by oxidation during tea fermentation. Food Chem 2011; 129:1331-1342.
- Li T and Zhang H: Identification and Comparative Determination of Rhodionin in Traditional Tibetan Medicinal Plants of Fourteen Rhodiola Species by High-Performance Liquid Chromatography-Photodiode Array Detection and Electrospray Ionization-Mass Spectrometry. Chem Pharm Bull 2008; 56: 807-814.
- Gao J, Xu P, Wang Y, Wang Y and Hochstetter D: Combined Effects of Green Tea Extracts, Green Tea Polyphenols or Epigallocatechin Gallate with Acarbose on Inhibition against α-Amylase and α-Glucosidase in-vitro. Molecules 2013; 18:11614-1623.
- Kwon YI, Jang HD and Shetty K: Evaluation of Rhodiola crenulata and Rhodiola rosea for management of Type II diabetes and hypertension. Asia Pac J Clin Nutr 2006; 15: 425-432
How to cite this article:
Santhosh S: Pharmacophore Mapping and Virtual Screening on Human Alpha Amylase Inhibitors. Int J Pharm Sci Res 2015; 6(5): 2127-32.doi: 10.13040/IJPSR.0975-8232.6(5).2127-32.
All © 2013 are reserved by International Journal os Pharmaceutical Sciences and Research. This Journal licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 3.0 Unported License.
Article Information
38
2127-32
782
1564
English
Ijpsr
S. Santhosh
Department of Bioinformatics, School of Chemical and Biotechnology, SASTRA University, Thirumalaisamudram, Thanjavur - 613401, India.
santhosharni@gmail.com
26 September, 2014
22 November, 2014
20 January, 2015
10.13040/IJPSR.0975-8232.6(5).2127-32
01 May, 2015