Glycogen synthase kinase-3 beta (GSK-3β) is an enzyme pertinently linked to neurodegenerative diseases since it is associated with the regulation of key neuropathological features in the central nervous system. Among the different kinds of inhibitors of this kinase, the allosteric ones stand out due to their selective and subtle modulation, lowering the chance of producing side effects. The mechanism of GSK-3β allosteric modulators may be considered still vague in terms of elucidating a well-defined binding pocket and a bioactive pose for them. In this context, we propose to reinvestigate and reinforce such knowledge by the application of an extensive set of in silico methodologies, such as cavity detection, ligand 3D shape analysis and docking (with robust validation of corresponding protocols), and molecular dynamics. The results here obtained were consensually consistent in furnishing new structural data, in particular by providing a solid bioactive pose of one of the most representative GSK-3β allosteric modulators. We further applied this to the prospect for new compounds by ligand-based virtual screening and analyzed the potential of the two obtained virtual hits by quantum chemical calculations. All potential hits achieved will be subsequently tested by in vitro assays in order to validate our approaches as well as to unveil novel chemical entities as GSK-3β allosteric modulators.
Keywords: GSK-3β; allosteric modulators; binding pose; cavity detection; computer-aided drug design; docking; neurodegenerative diseases; shape similarity.