Absolute FKBP binding affinities obtained via nonequilibrium unbinding simulations

We compute the absolute binding affinities for two ligands bound to the FKBP protein using nonequilibrium unbinding simulations. The methodology is straightforward requiring little or no modification to many modern molecular simulation packages. The approach makes use of a physical pathway, eliminating the need for complicated alchemical decoupling schemes. We compare our nonequilibrium results to those obtained via a fully equilibrium approach and to experiment. The results of this study suggest that to obtain accurate results using nonequilibrium approaches one should use the stiff-spring approximation with the second cumulant expansion. From this study we conclude that nonequilibrium simulation could provide a simple means to estimate protein-ligand binding affinities.