Pinpointing top inhibitors for GSK3β from pool of indirubin derivatives using rigorous computational workflow and their validation using molecular dynamics (MD) simulations

Abstract

Glycogen synthase kinase-3β (GSK3β) is a pivotal protein kinase implicated in a spectrum of debilitating diseases, encompassing cancer, diabetes, and neurodegenerative disorders. While the therapeutic potential of GSK3β inhibition is widely recognized, there remains an unmet need for a rigorous, systematic analysis probing the theoretical inhibition dynamics of a comprehensive library of indirubin derivatives against GSK3β using advanced computational methodologies. Addressing this gap, this study embarked on an ambitious endeavor, leveraging indirubin—a renowned scaffold—as a template to curate a vast library of 1000 indirubin derivatives from PubChem. These were enriched with varied substitutions and modifications, identified via a structure similarity search with a Tanimoto similarity threshold of 85%. Harnessing a robust virtual screening workflow, we meticulously identified the top 10 contenders based on XP docking scores. Delving deeper, we gauged the binding free energy differentials (ΔGBind) of these hits, spotlighting the top three compounds that showcased unparalleled binding prowess. A comparative pharmacophore feature mapping with the reference inhibitor OH8, co-crystallized with GSK3β (PDB ID: 6Y9R), was undertaken. The binding dynamics of these elite compounds were further corroborated with 100 ns molecular dynamics simulations, underlining their stable and potent interactions with GSK3β. Remarkably, our findings unveil that these indirubin derivatives not only match but, in certain scenarios, surpass the binding affinity and specificity of OH8. By bridging this research chasm, our study amplifies the therapeutic promise of indirubin derivatives, positioning them as frontrunners in the quest for groundbreaking GSK3β inhibitors, potentially revolutionizing treatments for a myriad of ailments.