A multiscale simulation approach to modelling drug-protein binding kinetics

Susanta Haldar, Federico Comitani, Giorgio Saladino, Christopher Woods, Marc Van der Kamp, Adrian Mulholland, Francesco Luigi Gervasio

Research output: Contribution to journalArticle (Academic Journal)peer-review

26 Citations (Scopus)
361 Downloads (Pure)


Drug–target binding kinetics has recently emerged as a sometimes critical determinant of in vivo efficacy and toxicity. Its rational optimization to improve potency or reduce side effects of drugs is, however, extremely difficult. Molecular simulations can play a crucial role in identifying features and properties of small ligands and their protein targets affecting the binding kinetics, but significant challenges include the long time scales involved in (un)binding events and the limited accuracy of empirical atomistic force fields (lacking, e.g., changes in electronic polarization). In an effort to overcome these hurdles, we propose a method that combines state-of-the-art enhanced sampling simulations and quantum mechanics/molecular mechanics (QM/MM) calculations at the BLYP/VDZ level to compute association free energy profiles and characterize the binding kinetics in terms of structure and dynamics of the transition state ensemble. We test our combined approach on the binding of the anticancer drug Imatinib to Src kinase, a well-characterized target for cancer therapy with a complex binding mechanism involving significant conformational changes. The results indicate significant changes in polarization along the binding pathways, which affect the predicted binding kinetics. This is likely to be of widespread importance in binding of ligands to protein targets.
Original languageEnglish
Number of pages9
JournalJournal of Chemical Theory and Computation
Early online date13 Sept 2018
Publication statusE-pub ahead of print - 13 Sept 2018

Cite this