Investigating Drug-Target Residence Time in Kinases through Enhanced Sampling Simulations

It is widely accepted that drug-target association and dissociation rates directly affect drug efficacy and safety. To rationally optimize drug binding kinetics, one must know the atomic arrangement of the protein-ligand complex during the binding/unbinding process in order to detect stable and metastable states. Whereas experimental approaches can determine kinetic constants with fairly good accuracy, computational approaches based on molecular dynamics (MD) simulations can deliver the atomistic details of the unbinding process. Furthermore, they can also be utilized prospectively to predict residence time (i.e. the inverse of unbinding kinetics constant, koff) with acceptable level of accuracy. Here, we report on a novel method based on adiabatic bias MD with an electrostatics-like collective variable (dubbed elABMD) for sampling protein-ligand dissociation events in two kinases. elABMD correctly ranked a ligand series on glucokinase, in agreement with experimental data and previous calculations. Subseq...