Designing of sorafenib analogs to target c-Raf for the management of hepatocellular carcinoma: Molecular dynamics and mmPBSA analysis.

Sorafenib remains the only approved treatment for advanced hepatocellular carcinoma (HCC), yet its clinical use is hindered by toxicity and the emergence of drug resistance. Sorafenib's anticancer effects are largely attributed to its inhibition of multiple kinases, including c-Raf, a key player in the Ras-Raf-MEK-ERK signaling cascade that promotes cell growth and survival. Given the critical role of c-Raf in tumor progression, targeting this kinase offers a promising strategy for improving therapeutic outcomes. Developing new analogs with stronger c-Raf inhibition, better pharmacokinetics, and reduced side effects could help address the current limitations of sorafenib. This study aimed to design novel sorafenib analogs with enhanced binding affinity and favorable pharmacokinetic profiles, specifically targeting the c-Raf kinase to increase therapeutic efficacy against HCC. By using a fragment replacement approach combined with computational methods, the goal was to identify candidates capable of forming stronger, more stable interactions with c-Raf, potentially overcoming resistance linked to sorafenib treatment. A total of 84 sorafenib analogs (A1-A84) were generated by modifying key functional groups, including the 2-picolinamide and substituted phenyl moieties known to influence kinase binding and anticancer activity. These analogs were evaluated through chemoinformatics and pharmacokinetic screening to assess their drug-likeness and safety. Molecular docking was performed to estimate their binding affinity toward c-Raf. Six top-performing analogs (A2, A6, A9, A20, A22, A63) were selected for further analysis. To evaluate their dynamic behavior, 100[Formula: see text]ns all-atom molecular dynamics simulations were conducted, followed by Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) calculations to determine binding free energies. Principal component analysis (PCA) was carried out to explore key motion patterns within the protein-ligand complexes. Molecular docking showed that the selected analogs exhibited stronger binding affinities (-11.6 to -10.9[Formula: see text]kcal/mol) compared to sorafenib (-9.3[Formula: see text]kcal/mol) and regorafenib (-9.5[Formula: see text]kcal/mol). Molecular dynamics simulations substantiated the docking results. MM-PBSA results revealed that at 100[Formula: see text]ns, the binding free energy for the c-Raf-sorafenib complex was 86.751[Formula: see text]kJ/mol, while the c-Raf complexes with A2, A6, A9, A20, A22, and A63 demonstrated significantly lower free energies of -129.114, -135.637, -136.242, -127.178, -94.25, and -123.176[Formula: see text]kJ/mol, respectively, indicating stronger and more stable binding. PCA further confirmed the stability and favorable dynamic profiles of these analogs trajectory with c-Raf. The improved binding affinities and lower free energies of the top analogs indicate that specific structural changes to sorafenib can enhance its effectiveness against c-Raf. Molecular dynamics and MM-PBSA results suggest the stability and strength of these interactions, particularly for A2, A6, and A9. This study identified six promising sorafenib analogs with improved binding affinity, favorable pharmacokinetic characteristics, and stable interactions with c-Raf. By focusing on c-Raf inhibition, the combined use of computational modeling, molecular simulations and mmPBSA analysis provided valuable insights for drug design. Among the candidates, A2, A6, and A9 emerged as promising drug candidates for further development, supporting the potential of targeting c-Raf to enhance therapeutic strategies against HCC.