Glycogen synthase kinase-3 (GSK-3) isoforms α and β have diverse roles within cell biology, and have been linked with multiple diseases that include prominent CNS conditions such as Alzheimer’s disease and several psychiatric disorders. In this study, motivated by computation, we aimed to identify novel ATP-binding site inhibitors of GSK-3 with CNS-active potential. A ligand screening (docking) protocol against GSK-3β was first optimized, employing an active/decoy benchmarking set, with the final protocol selected based on statistical performance analysis. The optimized protocol involved pre-filtering of ligands using a three-point 3D-pharmacophore, followed by Glide-SP docking applying hinge region hydrogen bonding constraints. Using this approach, the Biogenic subset of the ZINC15 compound database was screened, focused on compounds with potential for CNS-activity. Twelve compounds (generation I) were selected for experimental validation usingin vitroGSK-3β binding assays. Two hit compounds,1and2, with 6-amino-7H-benzo[e]perimidin-7-one and 1-(phenylamino)-3H-naphtho[1,2,3-de]quinoline-2,7-dione type scaffolds were identified with IC50values of 1.63 µM and 20.55 µM, respectively. Ten analogues of2(generation II) were selected for structure activity relationship (SAR) analysis and revealed four low micromolar inhibitors (<10 µM), with19(IC50= 4.1 µM)~five times more potent than initial hit compound2. Selectivity screening of low micromolar inhibitors14and19(comparing aryl- and alkyl-substituents) against 10 homologous kinases revealed unique selectivity profiles, with both compounds more potent against the GSK-3α isoform (IC50s~2 µM) and, additionally, inhibitors of PKBβ (IC50s < 25 µM). Compound14also inhibited ERK2 and19, PKCγ, but generally good selectivity for GSK-3 isoforms over the other kinases was observed. The compounds had excellent predicted oral bioavailability and CNS-activity profiles, presenting promising candidates for future testing in cellular models of disease.