Abstract
OBJECTIVES: Vascular endothelial growth factor receptor 2 (VEGFR2) is a critical therapeutic target in hepatocellular carcinoma (HCC) due to its role in angiogenesis and tumor progression. While several inhibitors are currently used, clinical utility is often limited by resistance and adverse effects, necessitating the discovery of novel therapeutic agents. The aim of this study was to identify and characterize novel, highly effective VEGFR2 inhibitors using an integrated computational pipeline to advance the development of new HCC treatments. METHODS: A comprehensive dataset from the ChEMBL database was curated and standardized for Quantitative Structure-Activity Relationship (QSAR) modeling. A binary classification framework was employed, where a Light Gradient Boosting Machine (LGBM) model demonstrated superior predictive performance. Two lead compounds and a reference were selected for in-depth molecular modeling. Their binding poses were predicted via molecular docking and subsequently subjected to 200 ns Molecular Dynamics (MD) simulations to assess stability and conformational dynamics. Thermodynamic binding affinities were calculated using the Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) method. RESULTS: The LGBM model achieved high accuracy and a robust Matthews Correlation Coefficient (MCC) on an independent test set. MD analysis, including Root Mean Square Deviation (RMSD) and Radius of Gyration (Rg), confirmed stable binding throughout the 200 ns trajectory. MMPBSA calculations validated the binding affinities, identifying van der Waals and electrostatic interactions as the primary driving forces for complex stability. CONCLUSION: This study successfully bridges machine learning with advanced molecular simulations, offering a validated workflow for the rational design and optimization of novel small-molecule VEGFR2 inhibitors.