ENERGY CALCULATIONS AND ANALYSIS OF HIV-1 PROTEASE INHIBITOR CRYSTAL-STRUCTURES

The interactions between HIV-1 protease and its bound inhibitors have been investigated by molecular mechanics calculations and by analysis of crystal structures of the complexes in order to determine general rules for inhibitor and substrate binding to the protease. Fifteen crystal structures of HIV-1 protease with different peptidomimetic inhibitors showed conservation of hydrogen bond interactions between the main chain C=O and NH groups of the inhibitors and the C=O and NH groups of the protease extending from P3 C=O to P3' NH. The mean length of the hydrogen bonds between the inhibitor and the flexible flaps and the conserved water molecule (2.9 Angstrom) is slightly shorter than the mean length of hydrogen bonds between the inhibitor and the more rigid active site region (3.1 Angstrom) of the protease. The two hydrogen bonds between the conserved water and P2 and P1' carbonyl oxygen atoms of the inhibitor are the shortest and are predicted to be important for the tight binding of inhibitors. Molecular mechanics analysis of three crystal structures of HIV-1 protease with different inhibitors with independent calculations using the programs Discover and Brugel gave an estimate of 56-68% for the contribution of all the inhibitor main chain atoms to the total calculated protease-inhibitor interaction energy. The contribution of individual inhibitor residues to the interaction energy was calculated using Brugel. The main chain atoms of residue P2 had a consistently large favorable contribution to the total interaction energy, probably due to the presence of the two short hydrogen bonds to the flexible flap. The contribution of individual inhibitor side chains depended on the size of the side chain and the presence of specific hydrogen bond interactions with the protease.