Application of free energy simulations to the binding of a transition-state-analogue inhibitor to HIV protease.

Free energy simulations (slow-change method) have been used to estimate quantitatively the ratio of the binding constants of (S) and (R) isomers of a novel HIV protease inhibitor, JG365. As a starting geometry, we used the X-ray crystallographic structure of a complex of HIV protease and JG365 provided by A. Wlodawer. According to our results the (S) configuration, i.e. the form previously identified experimentally, binds considerably more tightly to the protease (delta delta G degrees = 2.9 kcal/mol). When the (S) inhibitor is bound, there is a very strong preference for protonation of the Asp125 (rather than the Asp25) residue of the protease. This study is the first to apply a new method for quantitatively assessing the precision of free energies calculated by the slow-change method.