Molecular dynamics and ligand docking of a hinge region variant of South African HIV-1 subtype C protease

Abstract

HIV-1 protease is an important antiretroviral drug target due to its key role in viral maturation. Computational models have been successfully used in the past to understand the dynamics of HIV-1 protease variants. We performed molecular dynamics simulations and induced fit docking on a wild-type South African HIV-1 subtype C protease and an N37T↑V hinge region variant. The simulations were initiated in a cubic cell universe and run in explicit solvent, with the wild-type and variant proteases in the fully closed conformation and under periodic boundary conditions. The trajectory for each simulation totalled 20 ns. The results indicate that the N37T↑V hinge region mutation and insertion alter the molecular dynamics of the flap and hinge regions when compared to the wild-type protease. Specifically, the destabilisation of the hinge region allowed a larger and protracted opening of the flap region due to the formation of two key hinge/cantilever salt-bridges, which are absent in the wild-type protease. Domain-domain anti-correlation was observed between the flap and hinge region for both models. However, the N37T↑V variant protease displayed a lower degree of anti-correlation. The mutations affected the thermodynamic landscape of inhibitor binding as there were fewer observable chemical contacts between the N37T↑V variant protease and lopinavir, atazanavir and darunavir, respectively. These data elucidate the biophysical basis for the selection of hinge region insertion mutations by the HI virus.