Theoretical Studies on the Selectivity Mechanisms of Glycogen Synthase Kinase 3β (GSK3β) with Pyrazine ATP-competitive Inhibitors by 3DQSAR, Molecular Docking, Molecular Dynamics Simulation and Free Energy Calculations

Background: Glycogen synthase kinase-3 (GSK3) is associated with various key biological processes and has been considered as an important therapeutic target for the treatment of many diseases. Great efforts have been made on the development of GSK3 inhibitors, especially ATP-competitive GSK3β inhibitor, but it is still a great challenge to develop selective GSK3β inhibitors because of the high sequence homology with other kinases. Objective: In order to reveal the selectivity mechanisms of GSK3β inhibition at the molecular level, a series of ATP-competitive GSK3β inhibitor was analyzed by a systematic computational method, combining 3DQSAR, molecular docking, molecular dynamic simulations and free energy calculations. Methods: Firstly, 3D-QSAR with CoMFA was built to explore the general structure activity relationships. Secondly, CDOCKER and Flexible docking were employed to predicted the reasonable docking poses of all studied inhibitors. And then, both GSK3β and CDK2 complexes were selected to conduct molecular dynamics simulations. Finally, the free energy calculations were employed to find the key selective-residues. Results: CoMFA model suggested the steric, hydrophobic fields play key roles in the bioactivities of inhibitors, and the binding mechanisms were well analyzed through molecular docking. The binding free energies predicted are in good agreement with the experimental bioactivities and the free energy calculations showed that the binding of GSK3β/inhibitors was mainly contributed from hydrogen bonding and hydrophobic interaction. Conclusion: Some key residues for selective binding were highlighted, which may afford important guidance for the rational design of novel ATP-competitive GSK3β inhibitors.