Open AccessAccounting for the Central Role of Interfacial Water in Protein–Ligand Binding Free Energy Calculations
Author(s) -
Ido Ben-Shalom,
Zhihua Lin,
Brian K. Radak,
Charles Lin,
Woody Sherman,
Michael K. Gilson
Publication year - 2020
Publication title -
journal of chemical theory and computation
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 2.001
H-Index - 185
eISSN - 1549-9626
pISSN - 1549-9618
DOI - 10.1021/acs.jctc.0c00785
Subject(s) - free energy perturbation , molecular dynamics , context (archaeology) , binding energy , thermodynamic integration , energy landscape , chemistry , ligand (biochemistry) , chemical physics , molecule , umbrella sampling , solvent models , computational chemistry , metadynamics , water model , monte carlo method , statistical physics , solvation , physics , thermodynamics , atomic physics , paleontology , biochemistry , receptor , organic chemistry , statistics , mathematics , biology
Rigorous binding free energy methods in drug discovery are growing in popularity because of a combination of methodological advances, improvements in computer hardware, and workflow automation. These calculations typically use molecular dynamics (MD) to sample from the Boltzmann distribution of conformational states. However, when part or all of the binding sites is inaccessible to the bulk solvent, the time needed for water molecules to equilibrate between bulk solvent and the binding site can be well beyond what is practical with standard MD. This sampling limitation is problematic in relative binding free energy calculations, which compute the reversible work of converting ligand 1 to ligand 2 within the binding site. Thus, if ligand 1 is smaller and/or more polar than ligand 2, the perturbation may allow additional water molecules to occupy a region of the binding site. However, this change in hydration may not be captured by standard MD simulations and may therefore lead to errors in the computed free energy. We recently developed a hybrid Monte Carlo/MD (MC/MD) method, which speeds up the equilibration of water between bulk solvent and buried cavities, while sampling from the intended distribution of states. Here, we report on the use of this approach in the context of alchemical binding free energy calculations. We find that using MC/MD markedly improves the accuracy of the calculations and also reduces hysteresis between the forward and reverse perturbations, relative to matched calculations using only MD with or without the crystallographic water molecules. The present method is available for use in AMBER simulation software.