In this study, 2 groups of 10 modified ligand systems with modified P3 and P2 side chains are used to study the binding mechanism with thrombin. Experimental results show that the binding affinity is enhanced by complex ligand side chains. The binding free energy obtained from the polarized protein-specific charge (PPC) force field combined with the newly developed interaction entropy (IE) method is consistent with the experimental values with a high correlation coefficient. On the contrary, poor correlation is obtained using the traditional normal mode (Nmode) method for calculating the entropy change. Furthermore, the binding free energy and hot-spot residue energy are decomposed, and the common hot-spot residues in the two groups of systems are Trp50, Leu96, Ile179, Asp199, Cyx201, Ser226, Trp227, Gly228, and Gly230. The electrostatic and van der Waals interaction energies were found to be the main contributors in the binding energy difference. CH-π and CH-CH interactions of Leu96 ligands are significantly related to the energy change due to the modified side chain, and the hydrogen bond between Asp199 and the ligand provides a strong electrostatic interaction, contributing to the binding free energy. Investigating the B-factor, principal component, and binding pocket also explains the change in the binding affinity caused by the modified side chains in ligands from the viewpoint of conformational change. This study demonstrates that the new IE method is superior to the Nmode method in the predicting binding free energy and emphasizes the importance of electronic polarization in molecular dynamics simulation. Moreover, from the viewpoint of energy and structure analysis, this study reveals the origin of the change in binding free energy in modified ligands with different binding sites.
Read full abstract