Calculation of protein-ligand binding free energy using smooth reaction path generation (SRPG) method: a comparison of the explicit water model, gb/sa model and docking score function.
We compared the protein-ligand binding free energies (G) obtained by the explicit water model, the MM-GB/SA (molecular-mechanics generalized Born surface area) model, and the docking scoring function. The free energies by the explicit water model and the MM-GB/SA model were calculated by the previously developed Smooth Reaction Path Generation (SRPG) method. In the SRPG method, a smooth reaction path was generated by linking two coordinates, one a bound state and the other an unbound state. The free energy surface along the path was calculated by a molecular dynamics (MD) simulation, and the binding free energy was estimated from the free energy surface. We applied these methods to the streptavidin-and-biotin system. The G value by the explicit water model was close to the experimental value. The G value by the MM-GB/SA model was overestimated and that by the scoring function was underestimated. The free energy surface by the explicit water model was close to that by the GB/SA model around the bound state (distances of < 6 A), but the discrepancy appears at distances of > 6 A. Thus, the difference in long-range Coulomb interaction should cause the error in G. The scoring function cannot take into account the entropy change of the protein. Thus, the error of G could depend on the target protein.
Mitomo
Daisuke
D
Japan Biological Informatics Consortium (JBIC), 2-41-6, Aomi, Koto-ku, Tokyo 135-0064, Japan. d-mitomo@aist.go.jp
Fukunishi
Yoshifumi
Y
Higo
Junichi
J
Nakamura
Haruki
H
D008962:Models, Theoretical; D011485:Protein Binding; D013816:Thermodynamics; D014867:Water
D014867:Water