Large scale relative protein ligand binding affinities using non-equilibrium alchemy

2020 | journal article

Jump to:Cite & Linked | Documents & Media | Details | Version history

Cite this publication

​Large scale relative protein ligand binding affinities using non-equilibrium alchemy​
Gapsys, V. ; Pérez-Benito, L.; Aldeghi, M.; Seeliger, D.; van Vlijmen, H.; Tresadern, G. & de Groot, B. L. ​ (2020) 
Chemical Science11(4) pp. 1140​-1152​.​ DOI: https://doi.org/10.1039/C9SC03754C 

Documents & Media

document.pdf1.42 MBAdobe PDF

License

GRO License GRO License

Details

Authors
Gapsys, Vytautas ; Pérez-Benito, Laura; Aldeghi, Matteo; Seeliger, Daniel; van Vlijmen, Herman; Tresadern, Gary; de Groot, Bert L. 
Abstract
Relative ligand binding affinity calculations based on molecular dynamics (MD) simulations and non-physical (alchemical) thermodynamic cycles have shown great promise for structure-based drug design.
Ligand binding affinity calculations based on molecular dynamics (MD) simulations and non-physical (alchemical) thermodynamic cycles have shown great promise for structure-based drug design. However, their broad uptake and impact is held back by the notoriously complex setup of the calculations. Only a few tools other than the free energy perturbation approach by Schrödinger Inc. (referred to as FEP+) currently enable end-to-end application. Here, we present for the first time an approach based on the open-source software pmx that allows to easily set up and run alchemical calculations for diverse sets of small molecules using the GROMACS MD engine. The method relies on theoretically rigorous non-equilibrium thermodynamic integration (TI) foundations, and its flexibility allows calculations with multiple force fields. In this study, results from the Amber and Charmm force fields were combined to yield a consensus outcome performing on par with the commercial FEP+ approach. A large dataset of 482 perturbations from 13 different protein–ligand datasets led to an average unsigned error (AUE) of 3.64 ± 0.14 kJ mol −1 , equivalent to Schrödinger's FEP+ AUE of 3.66 ± 0.14 kJ mol −1 . For the first time, a setup is presented for overall high precision and high accuracy relative protein–ligand alchemical free energy calculations based on open-source software.
Relative ligand binding affinity calculations based on molecular dynamics (MD) simulations and non-physical (alchemical) thermodynamic cycles have shown great promise for structure-based drug design.
Ligand binding affinity calculations based on molecular dynamics (MD) simulations and non-physical (alchemical) thermodynamic cycles have shown great promise for structure-based drug design. However, their broad uptake and impact is held back by the notoriously complex setup of the calculations. Only a few tools other than the free energy perturbation approach by Schrödinger Inc. (referred to as FEP+) currently enable end-to-end application. Here, we present for the first time an approach based on the open-source software pmx that allows to easily set up and run alchemical calculations for diverse sets of small molecules using the GROMACS MD engine. The method relies on theoretically rigorous non-equilibrium thermodynamic integration (TI) foundations, and its flexibility allows calculations with multiple force fields. In this study, results from the Amber and Charmm force fields were combined to yield a consensus outcome performing on par with the commercial FEP+ approach. A large dataset of 482 perturbations from 13 different protein–ligand datasets led to an average unsigned error (AUE) of 3.64 ± 0.14 kJ mol −1 , equivalent to Schrödinger's FEP+ AUE of 3.66 ± 0.14 kJ mol −1 . For the first time, a setup is presented for overall high precision and high accuracy relative protein–ligand alchemical free energy calculations based on open-source software.
Issue Date
2020
Journal
Chemical Science 
ISSN
2041-6520
eISSN
2041-6539
Language
English

Reference

Citations


Social Media