Unbinding Pathways of VEGFR2 Inhibitors Revealed by Steered Molecular Dynamics

A detailed atomistic description of the unbinding process of sorafenib and sunitinib, two known VEGFR2 inhibitors clinically used to treat renal cell carcinoma, was unraveled by using steered molecular dynamics (SMD) simulations. While sunitinib is a fast-dissociating binder, sorafenib exhibits quite a long residence time at this enzyme, which might impact its duration of action in vivo. In order to gain insights into the kinetically different behaviors of the two inhibitors, an SMD study was carried out, which involved a careful optimization of the force and velocity parameters. We were able to identify two different binding pathways for the two inhibitors, as sunitinib exited the ATP binding site from the cavity entrance without a rupture point while sorafenib moved opposite to the ATP binding site entrance. Furthermore, the calculated Delta G(off) values clearly reflect on a qualitative level the distinct off-rates of the two inhibitors, thus suggesting that this protocol could be tried on other VEGFR2 ligands to assess its robustness and then used to rank structural analogues of these derivatives.