Simple, Intuitive Calculations of Free Energy of Binding for Protein-Ligand Complexes. 3. Including the Free Energy Contribution of Structural Water Molecules

Structural water molecules within protein active sites are relevant for ligand-protein recognition because they modify the active site geometry and contribute to binding affinity. In this work an analysis of the interactions between 23 ligands and dimeric HIV-1 protease is reported. The X-ray structures of these complexes show the presence of four types of structural water molecules: water 301 (on the symmetry axis), water 313, water 313bis, and peripheral waters. Except for water 301, these are generally complemented with a symmetry-related set. The GRID program was used both for checking water locations and for placing water molecules that appear to be missing from the complexes due to crystallographic uncertainty. Hydropathic analysis of the energetic contributions using HINT indicates a significant improvement of the correlation between HINT scores and the experimentally determined binding constants when the appropriate bridging water molecules are taken into account. In the absence of water r(2) = 0.30 with a standard error of +/- 1.30 kcal mol(-1) and when the energetic contributions of the constrained waters are included r(2) = 0.61 with a standard error of 0,98 kcal mol(-1). HINT was shown to be able to map quantitatively the contribution of individual structural waters to binding energy. The order of relevance for the various types of water is water 301 > water 313 > water 313bis > peripheral waters. Thus, to obtain the most reliable free energy predictions, the contributions of structural water molecules should be included. However, care must be taken to include the effects of water molecules that add information value and not just noise.