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 hotspot residues in the two groups of systems are Trp50,