Binding cooperativity between a ligand carbonyl group and a hydrophobic side chain can be enhanced by additional H-bonds in a distance dependent manner: A case study with thrombin inhibitors

One of the underappreciated non-covalent binding factors, which can significantly affect ligand-protein binding affinity, is the cooperativity between ligand functional groups. Using four different series of thrombin inhibitors, we reveal a strong positive cooperativity between an H-bond accepting carbonyl functionality and the adjacent P3 hydrophobic side chain. Adding an H-bond donating amine adjacent to the P3 hydrophobic side chain further increases this positive cooperativity thereby improving the K_i by as much as 546-fold. In contrast, adding an amidine multiple H-bond/salt bridge group in the distal S1 pocket does not affect this cooperativity. An analysis of the crystallographic B-factors of the ligand groups inside the binding site indicates that the strong cooperativity is mainly due to a significant mutual reduction in the residual mobility of the hydrophobic side chain and the H-bonding functionalities that is absent when the separation distance is large. This type of cooperativity is important to encode in binding affinity prediction software, and to consider in SAR studies.