Computational exploration of FOXM1 inhibitors for glioblastoma: an integrated virtual screening and molecular dynamics simulation study

In this study, a comprehensive investigation of a set of phytochemicals to identify potential inhibitors for the Forkhead box protein M1 (FOXM1) was conducted. FOXM1 is overexpressed in glioblastoma (GBM) cells and plays a crucial role in cell cycle progression, proliferation, and invasion. FOXM1 inhibitors have shown promising results in preclinical studies, and ongoing clinical trials are assessing their efficacy in GBM patients. However, there are limited studies on the identification of novel compounds against this attractive therapeutic target. To address this, the NPACT database containing 1,574 phytochemicals was used, employing a hierarchical multistep docking approach, followed by an estimation of relative binding free energy. By fixing user-defined XP-dock and MM-GBSA cut-off scores of −6.096 and −37.881 kcal/mol, the chemical space was further narrowed. Through exhaustive analysis of molecular binding interactions and various pharmacokinetics profiles, we identified four compounds, namely NPACT00002, NPACT01454, NPACT00856, and NPACT01417, as potential FOXM1 inhibitors. To assess the stability of protein-ligand binding in dynamic conditions, 100 ns Molecular dynamics (MD) simulations studies were performed. Furthermore, Molecular mechanics with generalized Born and surface area solvation (MM-GBSA) based binding free energy estimations of the entire simulation trajectories revealed a strong binding affinity of all identified compounds towards FOXM1, surpassing that of the control drug Troglitazone. Based on extensively studied multistep docking approaches, we propose that these molecules hold promise as FOXM1 inhibitors for potential therapeutic applications in GBM. However, experimental validation will be necessary to confirm their efficacy as targeted therapies.