Thermodynamic contributions of the ordered water molecule in HIV-1 protease

Binding between biomolecules is usually accompanied by the formation of direct interactions with displacement of water from the binding sites. In some cases, however, the interactions are mediated by ordered water molecules, whose effect on binding affinity and the other thermodynamic functions is unclear. In this work, we compute the contribution of one such water molecule, the strongly bound water molecule at the binding site of HIV-1 protease, to the thermodynamic properties using statistical mechanical formulas for the energy and entropy. The requisite correlation functions are obtained by molecular dynamics simulations. We find that the entropic penalty of ordering is large but is outweighed by the favorable water-protein interactions. We also find a large negative contribution from this water molecule to the heat capacity. This approach could be useful in rational drug design by estimating which bound water molecules would be most favorable to displace.