Resistant mechanism against nelfinavir of subtype C human immunodeficiency virus type 1 proteases

Human immunodeficiency virus type 1 protease (HIV-1 PR) is one of the major targets of anti-AIDS drug discovery, and the subtype C HIV-1 is infecting more and more humans. In this work, we executed computational simulations of subtype B (labeled B(WT)) HIV-1 PR, D30N mutant B (labeled B(D30N)) HIV-1 PR, subtype C (labeled C(Ref)) HIV-1 PR, D30N mutant C (labeled C(D30N)) HIV-1 PR, and D30N/N83T mutant C (labeled C(D30N/N83T)) HIV-1 PR with drug nelfinavir (NFV), aiming at clarifying (1) the resistant mechanism against NFV due to D30N mutation in subtype C HIV-1 PR; (2) the reason that the emergence rate of N83T mutation in C(D30N) HIV-1 PR is higher than that in B(D30N) HIV-1 PR; (3) the affinity of NFV with C(D30N/N83T) HIV-1 PR is higher than B(D30N) and C(D30N) HIV-1 PRs. The results indicate: (1) D30N mutation in subtype C HIV-1 PR reduces the hydrogen bond between the 30th residue and NFV, and the binding free energy contributions of some residues decrease; (2) the hydrogen bonds between the 83th/83′th residue and the 34th/34′th residue exist in both B(D30N) and C(D30N/N83T) HIV-1 PRs, while they disappear in C(D30N) HIV-1 PR. Meanwhile, the binding free energy contribution of N30 in C(D30N) is lower than that in B(D30N) and C(D30N/N83T); (3) N83T mutation makes some residues dislocate, and the contributions of these residues to binding free energy in C(D30N/N83T) increase comparing to those in B(D30N) and C(D30N). Our findings suggest that despite the nonactive site mutations, the polymorphisms regulate the emergence rates of these drug-resistant mutants.