Prediction of the Binding Mode Between GSK3β and a Peptide Derived from GSKIP Using Molecular Dynamics Simulation

GSK3 beta plays an important role in many physiological functions; dysregulated GSK3 beta is involved in human diseases such as diabetes, cancer, and Alzheimer's disease. This study uses MD simulations to determine the interaction between GSK3 beta and a peptide derived from GSKIP, a novel GSK3 beta interacting protein. Results show that GSKIPtide is inlaid in a binding pocket consisting of an alpha-helix and an extended loop near the carboxy-terminal end. This binding pocket is hydrophobic, and is responsible for the protein-protein interaction of two other GSK3 beta interacting proteins: FRAT and Axin. The GSKIPtide binding mode is closer to that of AxinGID (in the Axin-GSK3-interacting domain). The single-point mutations of V267G and Y288F in GSK3 beta differentiate the binding modes between GSK3 and GSKIPtide, AxinGID, and FRATide. The V2677G mutation of GSK3 beta reduces the GSKIPtide binding affinity by 70% and abolishes the binding affinity with AxinGID, but has no effect on FRATide. However, GSK3 beta Y288F completely abolishes the FRATide binding without affecting GSKIPtide or AxinGID binding. An analysis of the GSK3 beta-GSKIPtide complex structure and the X-ray crystal structures of GSK3 beta-FRATide and GSK3 beta-AxinGID complexes suggests that the hydroxyl group of Y288 is crucial to maintaining a hydrogen bond network in GSK3 beta-FRATide. The hydrophobic side chain of V267 maintains the integrity of helix-helix ridge-groove hydrophobic interaction for GSK3 beta-GSKIPtide and GSK3 beta-AxinGID. This study simulates these two mutant systems to provide atomic-level evidence of the aforementioned experimental results and validate the wild-type complex structure prediction. (C) 2011 Wiley Periodicals, Inc. Biopolymers 95: 461-471, 2011.