Abstract
Explicit water molecules in the binding site of proteins play a crucial role for protein-ligand association. Recent advances in computer-aided drug discovery methodology allow for an accurate prediction of the localized position and thermodynamic profile of water molecules (i.e., hydration sites) in the binding site. The underlying calculations are based on MD simulations of explicit water molecules in a restrained protein structure. However, the ligand-binding process is typically associated with protein conformational change that influences the position and thermodynamic properties of the hydration site. In this manuscript, we present the developments of two methods to incorporate the influence of protein conformational change on hydration sites either by following the conformational transition step-by-step (method I) or to match the hydration sites of the two transition end states using local coordinate systems (method II). Using these methods, we highlight the difference in the estimated protein desolvation free energy with and without inclusion of protein flexibility. To the best of our knowledge, this is the first study that explicitly studies the influence of protein conformational change on the position and thermodynamic profiles of water molecules and provides methodology to incorporate protein flexibility into the estimation of the desolvation free energy.