Transforming growth factor β type 1 receptor (TGFβR1), a crucial serine-threonine kinase, is central to the TGFβ/Smad signaling pathway, governing cellular processes like growth, differentiation, apoptosis, and immune response. This pathway is closely linked to the epithelial-mesenchymal transition (EMT) process, which plays an important role in the metastasis of hepatocellular carcinoma (HCC). To date, only limited inhibitors targeting TGFβR1 have entered clinical trials, yet they encounter challenges, notably high toxicity, in clinical applications. Herein, an efficient virtual screening pipeline was developed. Eighty compounds were screened from a pool of over 17 million molecules based on docking scores and binding free energy. Four compounds were manually selected with the assistance of enhanced sampling method BPMD (binding pose metadynamics). The binding stability of these four compounds complexed with TGFβR1 was subsequently studied through long-timescale conventional molecular dynamics simulations. The three most promising compounds were subjected to