Molecular docking of citrus flavo- noids with some targets related to diabetes

Citrus flavonoids isolated from citrus peel (flavedo and albedo) were taken as ligands for molecular docking. The molecular targets, (i.e. glucokinase, glycogen synthase kinases 3 , peroxisome proliferator-activated receptor gamma, and dipeptidyl peptidase IV) whose crystallographic structures are available on the PDB database as 1V4S, 1Q4L, 2PRG, 2ONC respectively, were used for the docking analysis using the Autodock tool V 4.2 and ADT v1.5.4 programs. The docking studies of the ligands citrus flavonoids with four different target proteins showed that citrus flavonoids are good molecules which dock well with various targets related to diabetes mellitus. The above results demonstrate that citrus flavonoids might be potentially used for blood glucose regulation. Article Info Received: 16 March 2013 Accepted: 22 March 2013 Available Online: 19 April 2013 DOI: 10.3329/bjp.v8i2.14240 Cite this article: Shen W. Lu YH. Molecular docking of citrus flavonoids with some targets related to diabetes. Bangladesh J Pharmacol. 2013; 8: 156-70. This work is licensed under a Creative Commons Attribution 3.0 License. You are free to copy, distribute and perform the work. You must attribute the work in the manner specified by the author or licensor. step in glycolysis, the main producer of ATP in the body (Balamurugan et al., 2012; Kamata et al., 2004). Therefore, GK would be an ideal drug target for type 2 diabetes (T2D) diseases because of its high impact in glucose homeostasis, and its activation results in lower blood glucose level irrespective of the cause of hyperglycemia. Glycogen synthase kinases 3 (GSK 3 , PDB: 1Q4L) belongs to the super family of mitogen-activated protein kinase. GSK 3 has been implicated in the development of insulin resistance and regulation of glycogen synthesis (Osolodkin et al., 011). It is one of the important targets in the treatment of T2D. Inhibitors of GSK-3 have antidiabetic properties because they improve insulin sensitivity, glycogen synthesis, and glucose metabolism in skeletal muscles of diabetic patients (Akhtar and Bharatam, 2012; Johnson et al., 2011; Khanfar et al., 2010; Wauwe, 2003). Peroxisome proliferator-activated receptor gamma (PPAR , PDB: 2PRG) belongs to the nuclear receptor super family of transcription factors and is an important regulator of target genes involved in glucose and lipid homeostasis (Choi et al., 2011; Nolte et al., 1998). PPAR and their targets have been invested as attractive therapeutic targets for T2D (Bruning et al., 2007; Cho et al., 2011; Maltarollo and Honório, 2012). Dipeptidyl peptidase IV (DPP IV, PDB: 2ONC) has become an attractive target of drug discovery and diabetes treatment (Drucker and Nauck, 2006; Verspohl, 2009). DPP IV is a membrane-bound, serine protease ectoenzyme found in numerous sites, including the kidney, intestine, and capillary endothelium. DPP IV is responsible for the degradation of a number of biological peptides including GLP-1 and GIP, which are incretins released from the gut in response to food and play an essential role in maintaining glucose homeostasis (Abu-Hamdah et al., 2009; Baggio and Drucker, 2007; Deacon, 2004; Verspohl, 2009). The above mentioned targets were subjected to molecular docking with a view to identify how citrus flavonoids play an important role in the process of hyperglycemia. Methods and Materials Retrieval of the three-dimensional structure of target proteins The structures of the target receptor binding sites of human glucokinase (PDB:1V4S), glycogen synthase kinases 3 (PDB:1Q4L), peroxisome proliferator-activated receptor gamma (PDB:2PRG), and dipeptidyl peptidase IV (PDB:2ONC) were obtained from the RCSB protein Data Bank, http://www.pdb.org/pdb/ home/home.do.


Introduction
Diabetes mellitus has reached epidemic proportions and affects more than 170 million individuals worldwide (Stumvoll et al., 2005).Compared with synthetic compounds, natural small molecules with special bioactivity have become the major resource of bioactive agents and played a key role in diabetes therapy (Liu et al., 2010).Therefore, management of diabetes without any side effects is still a challenge to the medical system.This leads to increasing demand for natural products with antidiabetic activity with less side effects.
Citrus flavonoids have received much attention in recent years, for its potential therapeutic qualities and relatively low toxicity to animals (Benavente-Garcia and Castillo, 2008;Choi and Ahn, 2008;Manthey et al., 2001).In our previous study, hesperidin, naringin, neohesperidin and nobiletin exhibited antidiabetic activities partly by binding to starch, delaying the starch digestion (Shen et al., 2012).And these flavonoids showed week inhibitory activity against digestive enzymes (e.g., pancreatic -amylase and -glucosidase).
Likewise, hesperidin, naringin and nobiletin also showed hypoglycemic effects by improving insulin sensitivity in diabetic animals (Akiyama et al., 2009;Akiyama et al., 2010;Jung et al., 2006;Jung et al., 2004;Lee et al., 2010).However, whether these citrus flavonoids regulated blood glucose via other target proteins or genes related to diabetes were unclear.
Bioinformatics tools have become very important to pinpoint the targets for different ligands (Osguthorpe et al., 2012).Using bioinformatics tools we tried to evaluate whether citrus flavonoids are good ligands to some of the target proteins or gene related to diabetes such as glucokinase, glycogen synthase kinases 3 , peroxisome proliferator-activated receptor gamma, and dipeptidyl peptidase IV.
Glucokinase (GK, PDB ID: 1V4S) is a monomeric cytoplasmic enzyme found in the liver and pancreas.Its main function is regulation of glucose levels in these organs.Through phosphorylation glucokinase is able to increase the metabolism of glucose.In the liver it increases the synthesis of glycogen and it is the first step in glycolysis, the main producer of ATP in the body (Balamurugan et al., 2012;Kamata et al., 2004).Therefore, GK would be an ideal drug target for type 2 diabetes (T2D) diseases because of its high impact in glucose homeostasis, and its activation results in lower blood glucose level irrespective of the cause of hyperglycemia.
Glycogen synthase kinases 3 (GSK 3 , PDB: 1Q4L) belongs to the super family of mitogen-activated protein kinase.GSK 3 has been implicated in the development of insulin resistance and regulation of glycogen synthesis (Osolodkin et al.,011).It is one of the important targets in the treatment of T2D.Inhibitors of GSK-3 have antidiabetic properties because they improve insulin sensitivity, glycogen synthesis, and glucose metabolism in skeletal muscles of diabetic patients (Akhtar and Bharatam, 2012;Johnson et al., 2011;Khanfar et al., 2010;Wauwe, 2003).
Peroxisome proliferator-activated receptor gamma (PPAR , PDB: 2PRG) belongs to the nuclear receptor super family of transcription factors and is an important regulator of target genes involved in glucose and lipid homeostasis (Choi et al., 2011;Nolte et al., 1998).PPAR and their targets have been invested as attractive therapeutic targets for T2D (Bruning et al., 2007;Cho et al., 2011;Maltarollo and Honório, 2012).
Dipeptidyl peptidase IV (DPP IV, PDB: 2ONC) has become an attractive target of drug discovery and diabetes treatment (Drucker and Nauck, 2006;Verspohl, 2009).DPP IV is a membrane-bound, serine protease ectoenzyme found in numerous sites, including the kidney, intestine, and capillary endothelium.DPP IV is responsible for the degradation of a number of biological peptides including GLP-1 and GIP, which are incretins released from the gut in response to food and play an essential role in maintaining glucose homeostasis (Abu-Hamdah et al., 2009;Baggio and Drucker, 2007;Deacon, 2004;Verspohl, 2009).
The above mentioned targets were subjected to molecular docking with a view to identify how citrus flavonoids play an important role in the process of hyperglycemia.

Ligand selection
Citrus flavonoids (i.e.hesperidin, naringin, neohesperidin, nobiletin) and other positive drugs were chosen from the National Centre for Biotechnology Informaton (NCBI) PubChem compound database.These molecules were downloaded in Structure Date File (SDF) format and converted to Protein Data Bank (PDB) coordinates by using Open Babel (http:// openbabel.org)converter.All the chemical structures of ligand compounds used in the study were shown in Figure 1.

Receptor and ligand optimization
PDB coordinates of the target receptor proteins and ligands molecules were optimized using Gromacs 4.0 suite force field analysis and UCSF Chimera (http:// www.cgl.ucsf.edu/chimera)tools, respectively.The optimized structures had minimum energy confirmation, which provided stability to the structure.These optimized receptors and ligands molecules were used for the docking study.

Docking analysis
The docking analysis of citrus flavonoids were carried out by means of the Autodock tools (ADT) v1.5.4 and autodock v4.2 program; (Autodock, Autogrid, Autotors, Copyright-1991e2000) from the Scripps Research Institute, http://www.scripps.edu/mb/olson/doc/autodock.To run autodock, we used a searching grid extended over ligand moieties, Kollman charges were assigned and atomic solvation parameteres were added.Polar hydrogen charges of the Gasteiger-type were assigned and nonpolar hydrogens were merged with the carbons and the internal degrees of freedom and torsions were set.Citrus flavonoids were docked to all the target protein complexes with the molecule considered as a rigid body.The search was extended over the whole receptor protein used as blind docking.Affinity maps for all the atom types present, as well as an electrostatic map, were computed with a grid spacing of 0.375.The search was carried out with the Lamarckian Genetic Algorithm; populations of 100 individuals with a mutation rate of 0.02 have been evolved for 10 generations.The remaining parameters were set as default.A root mean square deviation (RMSD) tolerance for each docking was set at 2.0.Evaluation of the results was done by sorting the different complexes with respect to the predicted binding energy.A cluster analysis based on root mean square deviation values, with reference to the starting geometry, was subsequently performed and the lowest energy conformation of the more populated cluster was

Results and Discussion
The docking simulations in the active sites of 1V4S, 1Q4L, 2PRG, and 2ONC were performed by the Auto dock program, which has been shown to successfully reproduce experimentally observed binding modes in terms of lowest docking energy.The target proteins structures of 1V4S, 1Q4L, 2PRG, and 2ONC were docked with citrus flavonoids, which provided excellent results as were seen by the least values of the binding energy.
The best possible binding modes of all the ligands at target protein's active sites are displayed in Figure 2 Figure 2 shows the result of docking analysis of human glucokinase (1V4S) with citrus flavonoids.Figure 2A showed the binding site of protein and metformin.Glucokinase protein residues Glu 216, Tyr 215, Glu 96 were formed H-bond with metformin.Figure 2B showed that hesperidin exhibited strong interaction with glucokinase protein.2D showed that nobiletin formed H-bonds with Thr 65, Tyr 215, Leu 451 residues.Figure 2E depicted that neohesperidin exhibited weak interaction with glucokinase via H-bond interaction with Ser 64, Tyr 215 residues.The glucokinase resultant binding energy, Ki, H-bond and other interaction with the active site residues are given below in the Table I.
Glucokinase is expressed only in liver and pancreatic beta cells and plays a key role in the regulation of glucose homeostasis.In the hepatocyte, the phosphorylation of glucose by glucokinase facilitates the uptake and metabolism of glucose by maintaining a gradient for glucose transport into these cells thereby regulating hepatic glucose disposal.In the beta cells, GK is believed to be part of the glucose-sensing mechanism and to be involved in the regulation of insulin release (Stoffel et al., 1992).During diabetes condition total or partial deficiency of insulin causes derangements in carbohydrate metabolism that decreases activity of several key enzymes including glucokinase resulting in the impaired glucose utilization and augmented hepatic glucose production.Chandramohan et al. reported that diabetic rats treated with 3-HMX active principle from plant increased GK activity (Chandramohan et al., 2008).In the same way citrus flavonoids increases GK activity by docking into GK's active sites, thereby increasing the utilization of glucose leading to decreased blood glucose level.3E showed the interaction between GSK-3 and nobiletin.Two Hbonds were formed between GSK-3 and nobiletin.The GSK-3 resultant binding energy, Ki, H-bond and other interaction with the active site residues are given below in the Table II.
GSK-3 can be inhibited through three distinct mechanisms: i), ATP non-competitive (in substrate interaction domain), ii), ATP competitive (in ATP binding pocket), and iii) metal ion competitive (in Mg2+ binding site).GSK-3 has been identified as an important kinase in the intercellular signaling pathway downstream from the insulin receptor.GSK-3 inactivates GS by phosphorylation resulting in glycogen synthesis inhibition (Khanfar et al., 2010;Wauwe, 2003).It has been reported that some flavonoids, such as luteolin, rutin, narirutin, etc., could bind with B-ring hydroxyls stabilized by hydrogen bonding with Arg 141 and Tyr 134 in the hinge; A-ring hydroxyls stabilized by hydrogen bonding with Asn 64, Gly 65, Lys 85, and Asp 200 residues in the glycine-rich loop (Johnson et al., 2011).Our results indicated that four citrus flavonoids could form H-bonds or via other interactions with the above amino acid residues.Therefore, citrus flavonoids showed GSK-3 inhibitory activity.Moreover, among the four tested citrus flavonoids, hesperidin and naringin present the better inhibitory activity against GSK-3 .
Figure 4 depicts the docking analysis of PPAR .Rosiglitazone is full agonists of PPAR .Figure 4A showed the hydrogen-bonding network of the rosiglitazone head group to conserved PPAR residues.Figure 4B showed the hydrogen-bonding network of the hesperidin to protein residues Tyr 327, Tyr 473, Ser 342, Glu 291, His 449, Ile 326, Arg 288. Figure 4C depicted the orientation of naringin bound in the active site of the PPAR crystal structure.Figure 4D illustrated hydrogen-bonding interaction of nobiletin to PPAR protein Arg 288, Glu 343 residues.Figure 4E showed the neohesperidin formed H-bonds with protein residues Leu 340, Glu 291, Gly 284.The PPAR resultant binding energy, Ki, H-bond and other interaction with the active site   Hit list the interacting residues, binding energy and inhibit contant of docked ligands to GSK-3 receptor       III.

Hit list the interacting residues, binding energy and inhibit contant of docked ligands to GK receptor
It has been reported that rosiglitazone could form hydrogen bonds with residues His 323, Tyr 473, His 449, Ser 289, Gln286, Tyr 327 of PPAR (Bruning et al., 2007;Choi et al., 2011).Our result was accordance with the reported result.Compared with rosiglitazone, four citrus flavonoids exhibited PPAR agonists.Bruning and his coworkers reported that BVT.13 and nTZDpa were part agonists of PPAR .They could contact the surface area of H3, making different hydrogen bonds with H3, as well as making more hydrophobic contacts overall.Both structures made several contacts with Arg 288, Ser 342 found on H3 (Bruning et al., 2007).Four tested citrus flavonoids exhibited hydrogen-bond interaction with H3 residues.Therefore, citrus flavonoids belong to part agonists.
Ligands docking to DPP IV protein active sites were shown in Figure 5. Figure 5A (Parmar et al., 2012).Hesperidin, naringin, neohesperidin obtained higher docking scores than sitagliptin.Interestingly, these three citrus flavonoids formed more hydrogen-bonds with DPP IV protein than sitagliptin.However, nobiletin formed only three Hydrogen-bonds with DPP IV protein and obtained the highest docking scores.These results indicated that the chemistry structure of citrus flavonoids might be contribute to the inhibitory activity of DPP IV.Hesperidin, naringin, neohesperidin have plenty of OH functional groups in the structures, which might promote the formation of H-bonds between flavonoids and protein residues.

Conclusion
Docking studies of the ligand citrus flavonoids with four different target proteins showed that citrus flavonoids are good molecules which dock well with various targets related to diabetes mellitus.Therefore, citrus flavonoids play important roles in blood glucose regulation, might via activation of GK and PPAR , whereas inhibition of GSK-3 and DPP IV.

Figure 1 :Figure 1 :
Figure 1: Structures of complexes used in this study -5 by using PYMOL tool v1.1.Ligands hydrogen-bonding to four target proteins and their corresponding energy values are listed in TableI-IV.

Table IV Hit list the interacting residues, binding energy and inhibit contant of docked ligands to DPPⅣ receptor
showed sitagliptin formed H-bonds with residues Glu 205, Glu 206, Tyr 662, Ser 209. Figure 5B exhibited the strong interaction between hesperidin and DPP IV protein residues Glu 205, Ser 209, Cys 551, Lys 554, Tyr 585, Trp 629, Tyr 631, Tyr 662. Figure 5C illustrated the naringin formed H-bonds with protein residues Val 546, Ser 209, Tyr 662, Trp 629, Glu 206, Arg 125, Lys 554.Figure5Ddepicted the interaction between neohesperidin and protein active sites.Figure5Eillustrated that nobiletin exhibited weak interaction with DPP IV protein via residues Trp 629, Ser 630, Tyr 631.The DPP IV resultant binding energy, Ki, H-bond and other interaction with the active site residues are given below in the TableIV.Parmar et al.indicated that sitagliptin exhibited good DPP IV inhibitory activity might be due to the its chemistry structure contains more electronegative groups F and N