Glycogen synthase kinase 3β (GSK-3β) is a potential therapeutic target for cancer, type-2 diabetes, and Alzheimer’s disease. This paper proposes a new lead identification protocol that predicts new GSK-3β ATP competitive inhibitors with topologically diverse scaffolds. First, three-dimensional quantitative structure–activity relationship (3D QSAR) models were built and validated. These models are based upon known GSK-3β inhibitors, benzofuran-3-yl-(indol-3-yl) maleimides, by means of comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). Second, 28 826 maleimide derivatives were selected from the PubChem database. After filtration via Lipinski’s rules, 10 429 maleimide derivatives were left. Third, the FlexX-dock program was employed to virtually screen the 10 429 compounds against GSK-3β. This resulted in 617 virtual hits. Fourth, the 3D QSAR models predicted that from the 617 virtual hits, 93 compounds would have GSK-3β inhibition values of less than 15 nM. Finally, from the 93 predicted active hits, 23 compounds were confirmed as GSK-3β inhibitors from literatures; their GSK-3β inhibition ranged from 1.3 to 480 nM. Therefore, the hits rate of our virtual screening protocol is greater than 25%. The protocol combines ligand- and structure-based approaches and therefore validates both approaches and is capable of identifying new hits with topologically diverse scaffolds.