Abstract
The degree of sequence variation in HIV-1 integrase genes among infected patients and their impact on clinical response to Anti retroviral therapy (ART) is of interest. Therefore, we collected plasma samples from 161 HIV-1 infected individuals for subsequent integrase gene amplification (1087 bp). Thus, 102 complete integrase gene sequences identified as HIV-1 subtype-C was assembled. This sequence data was further used for sequence analysis and multiple sequence alignment (MSA) to assess position specific frequency of mutations within pol gene among infected individuals. We also used biophysical geometric optimization technique based molecular modeling and docking (Schrodinger suite) methods to infer differential function caused by position specific sequence mutations towards improved inhibitor selection. We thus identified accessory mutations (usually reduce susceptibility) leading to the resistance of some known integrase inhibitors in 14% of sequences in this data set. The Stanford HIV-1 drug resistance database provided complementary information on integrase resistance mutations to deduce molecular basis for such observation. Modeling and docking analysis show reduced binding by mutants for known compounds. The predicted binding values further reduced for models with combination of mutations among subtype C clinical strains. Thus, the molecular basis implied for the consequence of mutations in different variants of integrase genes of HIV-1 subtype C clinical strains from South India is reported. This data finds utility in the design, modification and development of a representative yet an improved inhibitor for HIV-1 integrase.