Effect of set up protocols on the accuracy of alchemical free energy calculation over a set of ACK1 inhibitors

José M Granadino-Roldán, Antonia S J S Mey, Juan J Pérez González, Stefano Bosisio, Jaime Rubio-Martinez, Julien Michel

Research output: Contribution to journalArticlepeer-review

Abstract / Description of output

Hit-to-lead virtual screening frequently relies on a cascade of computational methods that starts with rapid calculations applied to a large number of compounds and ends with more expensive computations restricted to a subset of compounds that passed initial filters. This work focuses on set up protocols for alchemical free energy (AFE) scoring in the context of a Docking-MM/PBSA-AFE cascade. A dataset of 15 congeneric inhibitors of the ACK1 protein was used to evaluate the performance of AFE set up protocols that varied in the steps taken to prepare input files (using previously docked and best scored poses, manual selection of poses, manual placement of binding site water molecules). The main finding is that use of knowledge derived from X-ray structures to model binding modes, together with the manual placement of a bridging water molecule, improves the R2 from 0.45 ± 0.06 to 0.76 ± 0.02 and decreases the mean unsigned error from 2.11 ± 0.08 to 1.24 ± 0.04 kcal mol-1. By contrast a brute force automated protocol that increased the sampling time ten-fold lead to little improvements in accuracy. Besides, it is shown that for the present dataset hysteresis can be used to flag poses that need further attention even without prior knowledge of experimental binding affinities.

Original languageEnglish
Article numbere0213217
JournalPLoS ONE
Issue number3
Publication statusPublished - 12 Mar 2019


Dive into the research topics of 'Effect of set up protocols on the accuracy of alchemical free energy calculation over a set of ACK1 inhibitors'. Together they form a unique fingerprint.

Cite this