Neurodegenerative diseases (NDDs), including Alzheimer's disease (AD) and Parkinson's disease (PD), are a growing global health concern, especially among the elderly, posing significant challenges to well-being and survival. GSK3β, a serine/threonine kinase, is a key molecular player in the pathogenesis of NDDs. Dysregulated activity of GSK3β has been linked to neurodegenerative complications. Targeting GSK3β with active-site-specific inhibitors presents a promising therapeutic strategy for mitigating its pathological effects and potentially intercepting NDD progression. This study aimed to identify potential GSK3β inhibitors through an integrated in silico approach combining machine learning (ML)-based virtual screening, molecular docking, molecular dynamics (MD) simulations, and MM/GBSA binding free energy calculations. ML models were trained using known GSK3β inhibitors from BindingDB. Among all models, the Random Forest (RF) algorithm had the best prediction accuracy, with a value of 0.6832 on the test set and 0.7432 on the training set, and was employed to screen the target library of 11,032 phytochemicals. The ML-based screening identified 2,898 compounds with potential inhibitory action against GSK3β. Further screening based on Lipinski's Rule of Five gave 221 drug-like candidates. These compounds were further evaluated for GSK3β interaction via molecular docking. The analyses found ZINC136900288, ZINC7267, and ZINC519549 bind strongly and interact well with key residues in GSK3β active site with their binding scores being - 9.9, -8.8, and - 8.7 kcal/mol, respectively. MD simulations were conducted for both ligand-bound and apo GSK3β to assess structural stability. The simulation results showed that the ligand bound complexes were structurally stable with less fluctuations and higher conformational stability. In addition, (MM/GBSA) binding free energy calculations were carried out to quantify the affinity of the candidate compounds, and the candidate compound ZINC136900288 has the strongest binding affinity (-24.86 kcal/mol) of the three. Notably, these identified compounds feature novel chemical scaffolds that are structurally distinct from previously reported GSK3β inhibitors, emphasizing the originality and therapeutic potential of this study. These results show that ZINC136900288 may serve as suitable GSK3β inhibitors. Nevertheless, the efficacy and safety of these compounds need to be further validated experimentally and further studied in vivo for possible therapeutic application in NDDs.