Complexity in modeling and understanding protonation states: computational titration of HIV-1 protease inhibitor complexes

The computational-titration (CT) algorithm based on the 'natural' Hydropathic INTeractions (HINT) force field is described. The HINT software model is an empirical, non-Newtonian force field derived from experimentally measured partition coefficients for solvent transfer between octanol and H2O (log P-O/W). The CT algorithm allows the identification, modeling, and optimization of multiple protonation states of residues and ligand functional groups at the protein-ligand active site. The importance of taking into account pH and ionization states of residues, which strongly affect the process of ligand binding, for correctly predicting binding free energies is discussed. The application of the CT protocol to a set of six cyclic inhibitors in their complexes with HIV-1 protease is presented, and the advance of HINT as a virtual-screening tool is outlined.