The prediction of binding energies from the three-dimensional (3D) structure of a protein–ligand complex is an important goal of biophysics and structural biology. Here, we critically assess the use of empirical, solvent-accessible surface area-based calculations for the prediction of the binding of Src-SH2 domain with a series of tyrosyl phosphopeptides based on the high-affinity ligand from the hamster middle T antigen (hmT), where the residue in the pY+3 position has been changed. Two other peptides based on the C-terminal regulatory site of the Src protein and the platelet-derived growth factor receptor (PDGFR) are also investigated. Here, we take into account the effects of proton linkage on binding, and test five different surface area-based models that include different treatments for the contributions to conformational change and protein solvation. These differences relate to the treatment of conformational flexibility in the peptide ligand and the inclusion of proximal ordered solvent molecules in the surface area calculations. This allowed the calculation of a range of thermodynamic state functions (ΔCp, ΔS, ΔH, and ΔG) directly from structure. Comparison with the experimentally derived data shows little agreement for the interaction of SrcSH2 domain and the range of tyrosyl phosphopeptides. Furthermore, the adoption of the different models to treat conformational change and solvation has a dramatic effect on the calculated thermodynamic functions, making the predicted binding energies highly model dependent. While empirical, solvent-accessible surface area based calculations are becoming widely adopted to interpret thermodynamic data, this study highlights potential problems with application and interpretation of this type of approach. There is undoubtedly some agreement between predicted and experimentally determined thermodynamic parameters; however, the tolerance of this approach is not sufficient to make it ubiquitously applicable.

The crystal structures of four active site-directed thrombin inhibitors, 1–4, in a complex with human α-thrombin have been determined and refined at up to 2.0 Å resolution using X-ray crystallography. These compounds belong to a structurally novel family of inhibitors based on a 2,3-disubstituted benzo[b]thiophene structure. Compared to traditional active-site directed inhibitors, the X-ray crystal structures of these complexes reveal a novel binding mode. Unexpectedly, the lipophilic benzo[b]thiophene nucleus of the inhibitor appears to bind in the S1 specificity pocket. At the same time, the basic amine of the C-3 side chain of the inhibitor interacts with the mostly hydrophobic proximal, S2, and distal, S3, binding sites. The second, basic amine side chain at C-2 was found to point away from the active site, occupying a location between the S1 and S′1 sites. Together, the aromatic rings of the C-2 and C-3 side chains sandwich the indole ring of Trp60D contained in the thrombin S2 insertion loop defined by the sequence “Tyr-Pro-Pro-Trp.” [The thrombin residue numbering used in this study is equivalent to that reported for chymotrypsinogen (Hartley BS, Shotton DM, 1971, The enzymes, vol. 3. New York: Academic Press. pp 323–373).] In contrast to the binding mode of more classical thrombin inhibitors (d-Phe-Pro-Arg-H, NAPAP, Argatroban), this novel class of benzo[b]thiophene derivatives does not engage in hydrogen bond formation with Gly216 of the thrombin active site. A detailed analysis of the three-dimensional structures not only provides a clearer understanding of the interaction of these agents with thrombin, but forms a foundation for rational structure-based drug design. The use of the data from this study has led to the design of derivatives that are up to 2,900-fold more potent than the screening hit 1.