Abstract
The characteristics of NH…O, OH…O, and CH…O hydrogen bonds are examined in a group of 28 high-resolution crystal structures of protein–ligand complexes from the Protein Data Bank and compared with interactions found in small-molecule crystal structures from the Cambridge Structural Database. It is found that both strong and weak hydrogen bonds are involved in ligand binding. Because of the prevalence of multifurcation, the restrictive geometrical criteria set up for hydrogen bonds in small-molecule crystal structures may need to be relaxed in macromolecular structures. For example, there are definite deviations from linearity for the hydrogen bonds in protein–ligand complexes. The formation of CH…O hydrogen bonds is influenced by the activation of the CαH atoms and by the flexibility of the side-chain atoms. In contrast to small-molecule structures, anticooperative geometries are common in the macromolecular structures studied here, and there is a gradual lengthening as the extent of furcation increases. CH…O bonds formed by Gly, Phe, and Tyr residues are noteworthy. The numbers of hydrogen bond donors and acceptors agree with Lipinski's “rule of five” that predicts drug-like properties. Hydrogen bonds formed by water are also seen to be relevant in ligand binding. Ligand CH…Ow interactions are abundant when compared to NH…Ow and OH…Ow. This suggests that ligands prefer to use their stronger hydrogen bond capabilities for use with the protein residues, leaving the weaker interactions to bind with water. In summary, the interplay between strong and weak interactions in ligand binding possibly leads to a satisfactory enthalpy–entropy balance. The implications of these results to crystallographic refinement and molecular dynamics software are discussed. Proteins 2004. ©2003 Wiley-Liss, Inc.