Revealing the drug resistance mechanism of saquinavir due to G48V and V82F mutations in subtype CRF01_AE HIV-1 protease: molecular dynamics simulation and binding free energy calculations

Abstract

Human immunodeficiency virus-1 (HIV-1) protease is one of the important targets in AIDS therapy. The majority of HIV infections are caused due to non-B subtypes in developing countries. The co-occurrence of mutations along with naturally occurring polymorphisms in HIV-1 protease cause resistance to the FDA approved drugs, thereby posing a major challenge in the treatment of antiretroviral therapy. In this work, the resistance mechanism against SQV due to active site mutations G48V and V82F in CRF01_AE (AE) protease was explored. The binding free energy calculations showed that the direct substitution of valine at position 48 introduces a bulkier side chain, directly impairing the interaction with SQV in the binding pocket. Also, the intramolecular hydrogen bonding network of the neighboring residues is altered, indirectly affecting the binding of SQV. Interestingly, the substitution of phenylalanine at position 82 induces conformational changes in the 80’s loop and the flap region, thereby favoring the binding of SQV. The V82F mutant structure also maintains similar intramolecular hydrogen bond interactions as observed in AE-WT. Communicated by Ramaswamy H. Sarma