Accelerating Convergence of Free Energy Computations with Hamiltonian Simulated Annealing of Solvent (HSAS)

Coupling between binding of a ligand to a receptor and the displacement of a number of bound water molecules is a common event in molecular recognition processes. When the binding site is deeply buried and the exchange of water molecules with the bulk region is difficult to sample, the convergence and accuracy in free energy calculations can be severely compromised. Traditionally, Grand Canonical Monte Carlo (GCMC) based methods have been used to accelerate equilibration of water-at the expense, however, of lengthy trials before a molecular dynamics (MD) simulation. In this paper, a user-friendly and cost-efficient method, Hamiltonian simulated annealing of solvent in combination with lambda-exchange of free energy perturbation (FEP) is proposed to accelerate the sampling of water molecules in free energy calculations. As an illustrative example with reliable data from previous GCMC simulations, absolute binding affinity of camphor to cytochrome P450 was calculated. The simulated hydration state change in the buried binding pocket quantitatively agrees with GCMC simulations. It is shown that the new protocol significantly accelerates sampling of water in a buried binding pocket and the convergence of free energy, with negligible setup and computing costs compared to GCMC methods.