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.

Hydrogen atoms play a key role in protein-ligand recognition. They determine the quality of established H-bonding networks and define the protonation of bound ligands. Structural visualization of H atoms by X-ray crystallography is rarely possible. We used neutron diffraction to determine the positions of the hydrogen atoms in the ligands aniline and 2-aminopyridine bound to the archetypical serine protease trypsin. The resulting structures show the best resolution so far achieved for proteins larger than 100 residues and allow an accurate description of the protonation states and interactions with nearby water molecules. Despite its low pK of 4.6 and a large distance of 3.6 Å to the charged Asp189 at the bottom of the S1 pocket, the amino group of aniline becomes protonated, whereas in 2-aminopyridine, the pyridine nitrogen picks up the proton although its amino group is 1.6 Å closer to Asp189. Therefore, apart from charge-charge distances, tautomer stability is decisive for the resulting binding poses, an aspect that is pivotal for predicting correct binding.

Thermodynamics and kinetics of protein−ligand binding are both important aspects for the design of novel drug molecules. Presently, thermodynamic data are collected with isothermal titration calorimetry, while kinetic data are mostly derived from surface plasmon resonance. The new method of kinITC provides both thermodynamic and kinetic data from calorimetric titration measurements. The present study demonstrates the convenient collection of calorimetric data suitable for both thermodynamic and kinetic analysis for two series of congeneric ligands of human carbonic anhydrase II and correlates these findings with structural data obtained by macromolecular crystallography to shed light on the importance of shape complementarity for thermodynamics and kinetics governing a protein−ligand binding event. The study shows how minute chemical alterations change preferred ligand conformation and can be used to manipulate thermodynamic and kinetic signatures of binding. They give rise to the observation that analogous n-alkyl and nalkyloxy derivatives of identical chain length swap their binding kinetic properties at unchanged binding affinity.

Trypsin and thrombin, structurally similar serine proteases, recognize different substrates; thrombin cleaves after Arg, whereas trypsin cleaves after Lys/Arg. Both recognize basic substrate headgroups via Asp189 at the bottom of the S1 pocket. By crystallography and isothermal titration calorimetry (ITC), we studied a series of D-Phe/D-DiPhe-Pro-(amino)pyridines. Identical ligand pairs show the same binding poses. Surprisingly, one ligand binds to trypsin in protonated state and to thrombin in unprotonated state at P1 along with differences in the residual solvation pattern. While trypsin binding is mediated by an ordered water molecule, in thrombin, water is scattered over three hydration sites. Although having highly similar S1 pockets, our results suggest different electrostatic properties of Asp189 possibly contributing to the selectivity determinant. Thrombin binds a specific Na + ion next to Asp189, which is absent in trypsin. The electrostatic properties across the S1 pocket are further attenuated by charged Glu192 at the rim of S1 in thrombin, which is replaced by uncharged Gln192 in trypsin.

The fluorination of lead-like compounds is a common tool in medicinal chemistry to alter molecular properties in various ways and with different goals. We herein present a detailed study of the binding of fluorinated benzenesulfonamides to human Carbonic Anhydrase II by complementing macromolecular X-ray crystallographic observations with thermodynamic and kinetic data collected with the novel method of kinITC. Our findings comprise so far unknown alternative binding modes in the crystalline state for some of the investigated compounds as well as complex thermodynamic and kinetic structure-activity relationships. They suggest that fluorination of the benzenesulfonamide core is especially advantageous in one position with respect to the kinetic signatures of binding and that a higher degree of fluorination does not necessarily provide for a higher affinity or more favorable kinetic binding profiles. Lastly, we propose a relationship between the kinetics of binding and ligand acidity based on a small set of compounds with similar substitution patterns.