Machine Learning-Enhanced Calculation of Quantum-Classical Binding Free Energies.

Journal: Journal of chemical theory and computation
Published Date:

Abstract

Binding free energies are key elements in understanding and predicting the strength of protein-drug interactions. While classical free energy simulations yield good results for many purely organic ligands, drugs, including transition metal atoms, often require quantum chemical methods for an accurate description. We propose a general and automated workflow that samples the potential energy surface with hybrid quantum mechanics/molecular mechanics (QM/MM) calculations and trains a machine learning (ML) potential on the QM/MM energies and forces to enable efficient alchemical free energy simulations. To represent systems including many different chemical elements efficiently and to account for the different descriptions of QM and MM atoms, we propose an extension of element-embracing atom-centered symmetry functions for QM/MM data as an ML descriptor. The ML potential approach takes electrostatic embedding and long-range electrostatics into account. We demonstrate the applicability of the workflow on the well-studied protein-ligand complex of myeloid cell leukemia 1 and the inhibitor 19G and on the anticancer drug NKP1339 acting on the glucose-regulated protein 78.

Authors

  • Moritz Bensberg
    ETH Zurich Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland.
  • Marco Eckhoff
    ETH Zurich Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland.
  • F Emil Thomasen
    University of Copenhagen, Department of Biology, Linderstrøm-Lang Centre for Protein Science, Ole Maaløes Vej 5, DK-2200 Copenhagen N, Denmark.
  • William Bro-Jørgensen
    Department of Chemistry and Nano-Science Center, University of Copenhagen, Universitetsparken 5, DK-2100, Copenhagen Ø, Denmark. gsolomon@chem.ku.dk.
  • Matthew S Teynor
    University of Copenhagen Department of Chemistry and Nano-Science Center, Universitetsparken 5, DK-2100 Copenhagen Ø, Denmark.
  • Valentina Sora
    Department of Computer Science, University of Copenhagen, Copenhagen, Denmark.
  • Thomas Weymuth
    ETH Zurich Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland.
  • Raphael T Husistein
    ETH Zurich Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland.
  • Frederik E Knudsen
    Department of Biology, Linderstrøm-Lang Centre for Protein Science, University of Copenhagen, Ole Maaløes Vej 5, Copenhagen N DK-2200, Denmark.
  • Anders Krogh
    Department of Computer Science, University of Copenhagen, Copenhagen, Denmark.
  • Kresten Lindorff-Larsen
    Structural Biology and NMR Laboratory & Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen. Ole Maaløes Vej 5, DK-2200 Copenhagen N, Denmark. Electronic address: lindorff@bio.ku.dk.
  • Markus Reiher
    ETH Zurich Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland.
  • Gemma C Solomon
    Department of Chemistry and Nano-Science Center, University of Copenhagen, Universitetsparken 5, DK-2100, Copenhagen Ø, Denmark. gsolomon@chem.ku.dk.

Keywords

No keywords available for this article.