Molecular modeling and identification of novel glucokinase activators through stepwise virtual screening

The glucose phosphorylating enzyme glucokinase (GK) is a 50kD monomeric protein having 465 amino acids. It maintains glucose homeostasis inside cells, acts as a glucose sensor in pancreatic β-cells and as a rate controlling enzyme for hepatic glucose clearance and glycogen synthesis. It has two binding sites, one for binding d-glucose and the other for a putative allosteric activator named glucokinase activator (GKA). The GKAs interact with the same region of the GK enzyme that is commonly affected by naturally occurring mutations in humans. However, many GKAs do not bind to GK in the absence of glucose. Recently, it has been reported that GKAs are highly effective in patients with type 2 diabetes mellitus. In this milieu a molecular modeling study has been carried out on three natural variants of GK that lie in the GKA binding site and are known to cause maturity onset diabetes of young (MODY). Additionally, a 10ns molecular dynamics simulation was done on each of the modeled variant in order to explore the flexibility of this site. Subsequently, a systematic virtual screening study was done to identify compounds which can bind with high affinity at GKA binding site of mutant GK.