Resistance mechanism of human immunodeficiency virus type-1 protease to inhibitors: A molecular dynamic approach

Human immunodeficiency virus type 1 (HIV-1) protease inhibitors comprise an important class of drugs used in HIV treatments. However, mutations of protease genes accelerated by low fidelity of reverse transcriptase yield drug resistant mutants of reduced affinities for the inhibitors. This problem is considered to be a serious barrier against HIV treatment for the foreseeable future. In this study, molecular dynamic simulation method was used to examine the combinational and additive effects of all known mutations involved in drug resistance against FDA approved inhibitors. Results showed that drug resistant mutations are not randomly distributed along the protease sequence; instead, they are localized on flexible or hot points of the protein chain.  Substitution of more hydrophobic residues in flexible points of protease chains tends to increase the folding, lower the flexibility and decrease the active site area of the protease. The reduced affinities of HIV-1 protease for inhibitors seemed to be due to substantial decrease in the size of the active site and flap mobility. A correlation was found between the binding energy of inhibitors and their affinities for each mutant suggesting the distortion of the active site geometry in drug resistance by preventing effective fitting of inhibitors into the enzymes' active site. To overcome the problem of drug resistance of HIV-1 protease, designing inhibitors of variable functional groups and configurations is proposed.