Prediction of Protein-Ligand Binding Pose and Affinity Using the gREST plus FEP Method

The accurate prediction of protein-ligand binding affinity is a central challenge in computational chemistry and in-silico drug discovery. The free energy perturbation (FEP) method based on molecular dynamics (MD) simulation provides reasonably accurate results only if a reliable structure is available via high-resolution X-ray crystallography. To overcome the limitation, we propose a sequential prediction protocol using generalized replica exchange with solute tempering (gREST) and FEP. At first, ligand binding poses are predicted using gREST, which weakens protein-ligand interactions at high temperatures to sample multiple binding poses. To avoid ligand dissociation at high temperatures, a flat-bottom restraint potential centered on the binding site is applied in the simulation. The binding affinity of the most reliable pose is then calculated using FEP. The protocol is applied to the bindings of ten ligands to FK506 binding proteins (FKBP), showing the excellent agreement between the calculated and experimental binding affinities. The present protocol, which is referred to as the gREST+FEP method, would help to predict the binding affinities without high-resolution structural information on the ligand-bound state.