Hydrogen bonds are key interactions determining protein-ligand binding affinity and therefore fundamental to any biological process. Unfortunately, explicit structural information about hydrogen positions and thus H-bonds in protein-ligand complexes is extremely rare and similarly the important role of water during binding remains poorly understood. Here, we report on neutron structures of trypsin determined at very high resolutions ≤1.5 Å in uncomplexed and inhibited state complemented by X-ray and thermodynamic data and computer simulations. Our structures show the precise geometry of H-bonds between protein and the inhibitors N-amidinopiperidine and benzamidine along with the dynamics of the residual solvation pattern. Prior to binding, the ligand-free binding pocket is occupied by water molecules characterized by a paucity of H-bonds and high mobility resulting in an imperfect hydration of the critical residue Asp189. This phenomenon likely constitutes a key factor fueling ligand binding via water displacement and helps improving our current view on water influencing protein–ligand recognition.

Biophysical parameters can accelerate drug development; e.g., rigid ligands may reduce entropic penalty and improve binding affinity. We studied systematically the impact of ligand rigidification on thermodynamics using a series of fasudil derivatives inhibiting protein kinase A by crystallography, isothermal titration calorimetry, nuclear magnetic resonance, and molecular dynamics simulations. The ligands varied in their internal degrees of freedom but conserve the number of heteroatoms. Counterintuitively, the most flexible ligand displays the entropically most favored binding. As experiment shows, this cannot be explained by higher residual flexibility of ligand, protein, or formed complex nor by a deviating or increased release of water molecules upon complex formation. NMR and crystal structures show no differences in flexibility and water release, although strong ligand-induced adaptations are observed. Instead, the flexible ligand entraps more efficiently water molecules in solution prior to protein binding, and by release of these waters, the favored entropic binding is observed.

Insights into the thermodynamic and kinetic signature of the transient opening of a protein-binding pocket resulting from accommodation of suitable substituents attached to a given parent ligand scaffold are presented. As a target, we selected human aldose reductase, an enzyme involved in the development of late-stage diabetic complications. To recognize a large scope of substrate molecules, this reductase opens a transient specificity pocket. The pocket-opening step was studied by X-ray crystallography, microcalorimetry, and surface plasmon resonance using a narrow series of 2-carbamoyl-phenoxy-acetic acid derivatives. Molecular dynamics simulations suggest that pocket opening occurs only once an appropriate substituent is attached to the parent scaffold. Transient pocket opening of the uncomplexed protein is hardly recorded. Hydration-site analysis suggests that up to five water molecules entering the opened pocket cannot stabilize this state. Sole substitution with a benzyl group stabilizes the opened state, and the energetic barrier for opening is estimated to be ∼5 kJ/mol. Additional decoration of the pocket-opening benzyl substituent with a nitro group results in a huge enthalpy-driven potency increase; on the other hand, an isosteric carboxylic acid group reduces the potency 1000-fold, and binding occurs without pocket opening. We suggest a ligand induced-fit mechanism for the pocket-opening step, which, however, does not represent the rate-determining step in binding kinetics.

Drug binding involves changes of the local water structure around proteins including water rearrangements across surface-solvation layers around protein and ligand portions exposed to the newly formed complex surface. For a series of thermolysin-binding phosphonamidates, we discovered that variations of the partly exposed P2'-substituents modulate binding affinity up to 10 kJ mol(-1) with even larger enthalpy/entropy partitioning of the binding signature. The observed profiles cannot be completely explained by desolvation effects. Instead, the quality and completeness of the surface water network wrapping around the formed complexes provide an explanation for the observed structure-activity relationship. We used molecular dynamics to compute surface water networks and predict solvation sites around the complexes. A fairly good correspondence with experimental difference electron densities in high-resolution crystal structures is achieved; in detail some problems with the potentials were discovered. Charge-assisted contacts to waters appeared as exaggerated by AMBER, and stabilizing contributions of water-to-methyl contacts were underestimated.

A ligand-binding study is presented focusing on thermodynamics of fragment expansion. The binding of four compounds with increasing molecular weight to protein kinase A (PKA) was analyzed. The ligands display affinities between low-micromolar to nanomolar potency despite their low molecular weight. Binding free energies were measured by isothermal titration calorimetry, revealing a trend toward more entropic and less enthalpic binding with increase in molecular weight. All protein-ligand complexes were analyzed by crystallography and solution NMR spectroscopy. Crystal structures and solution NMR data are highly consistent, and no major differences in complex dynamics across the series are observed that would explain the differences in the thermodynamic profiles. Instead, the thermodynamic trends result either from differences in the solvation patterns of the conformationally more flexible ligand in aqueous solution prior to protein binding as molecular dynamics simulations suggest, or from local shifts of the water structure in the ligand-bound state. Our data thus provide evidence that changes in the solvation pattern constitute an important parameter for the understanding of thermodynamic data in protein-ligand complex formation.