Role of structural water molecule in HIV protease-inhibitor complexes: a QM/MM study

Structural water molecule 301 found at the interface of HIV protease-inhibitor complexes function as a hydrogen bond (H-bond) donor to carbonyl groups of the inhibitor as well as H-bond acceptor to amide/amine groups of the flap region of the protease. In this study, six systems of HIV protease-inhibitor complexes were analyzed, which have the presence of this "conserved" structural water molecule using a two-layer QM/MM ONIOM method. The combination of QM/MM and QM method enabled the calculation of strain energies of the bound ligands as well as the determination of their binding energies in the ligand-water and ligand-water-protease complexes. Although the ligand experiences considerable strain in the protein bound structure, the H-bond interactions through the structural water overcomes this strain effect to give a net stability in the range of 16-24 kcal/mol. For instance, in 1HIV system, the strain energy of the ligand was 12.2 kcal/mol, whereas the binding energy associated with the structural water molecule was 20.8 kcal/mol. In most of the cases, the calculated binding energy of structural water molecule showed the same trend as that of the experimental binding free energy values. Further, the classical MD simulations carried out on 1HVL system with and without structural water 301 showed that this conserved water molecule enhances the H-bond dynamics occurring at the Asp-bound active site region of the protease-inhibitor system, and therefore it will have a direct influence on the mechanism of drug action.