Molecular dynamic simulations to investigate the structural impact of known drug resistance mutations on HIV-1C Integrase-Dolutegravir binding

Abstract

Resistance associated mutations (RAMs) threaten the long-term success of combination antiretroviral therapy (cART) outcomes for HIV-1 treatment. HIV-1 Integrase (IN) strand transfer inhibitors (INSTIs) have proven to be a viable option for highly specific HIV-1 therapy. The INSTI, Dolutegravir is recommended by the World Health Organization for use as first-line cART. This study aims to understand how RAMs affect the stability of IN, as well as the binding of the drug Dolutegravir to the catalytic pocket of the protein. Molecular modelling of HIV-1C IN was performed using the SWISS-MODEL webserver; with quality assessment performed using internal methods and external software tools. The site directed mutator webserver was used to predict destabilizing and/or stabilizing effects of known RAMs while FoldX confirmed any changes in protein energy upon introduction of mutation. Interaction analysis between neighbouring residues was done using PyMOL. Three randomly selected mutations were chosen for molecular dynamic simulation studies using Gromacs. Trajectory analysis included Root mean square deviation and fluctuation, Radius of gyration, Principal component analysis and Interaction analysis between Dolutegravir and protein residues. The structural quality assessment indicated high reliability of the HIV-1C IN tetrameric structure, with more than 90% confidence in modelled regions. Change in free energy for the G140S mutant indicated a stabilizing effect and simulation analysis showed it to affect structural stability and flexibility of the protein structure. This was further supported by the drug being expelled from the G140S mutant active site, as indicated by interaction analysis. Our findings suggest the G140S mutant has a strong effect on the HIV-1C IN protein structure and Dolutegravir binding and should be validated using laboratory-based experiments. This approach can be applied to determine the effect of other mutations/variants on HIV-1C integrase drug binding.