Estimation of Protein-Ligand Unbinding Kinetics Using Non-Equilibrium Targeted Molecular Dynamics Simulations

We here report on non-equilibrium targeted Molecular Dynamics simulations as tool for the estimation of protein-ligand unbinding kinetics. With this method, we furthermore investigate the molecular basis determining unbinding rates, correlating simulations with experimental data from SPR kinetics measurements and X-ray crystallography on two small molecule compound libraries bound to the N-terminal domain of the chaperone Hsp90. Within the investigated libraries, we find ligand conformational changes and protein-ligand nonbonded interactions as discriminators for unbinding rates. Ligands with flexible chemical scaffold may remain longer at the protein target if they need to pass through extended conformations upon unbinding, or if they exhibit strong electrostatic and/or van der Waals interactions with the target. Ligands with rigid chemical scaffold can exhibit longer residence times if they need to perform any kind of conformational change for unbinding, while electrostatic interactions with the protein can facilitate unbinding. Our results show that understanding the unbinding pathway and the protein-ligand interactions along this path is crucial for the prediction of small molecule ligands with defined unbinding kinetics.