A computational exploration of global and temporal dynamics of selection pressure on HIV-1 Vif polymorphism

Virion infectivity factor (Vif), an accessory protein of HIV-1 (human immunodeficiency virus type 1), antagonizes host APOBEC3 protein (apolipoprotein B mRNA editing enzyme, catalytic polypeptide 3) or A3 via proteasomal degradation, facilitating viral replication. HLA (Human leukocyte antigens) alleles, host restriction factors, and error-prone reverse transcription contribute to the global polymorphic dynamics of HIV, impacting effective vaccine design. Our computational analysis of over 50,000 HIV-1 M vif sequences from the Los Alamos National Laboratory (LANL) database (1998-2021) revealed positive selection pressure on the vif gene (nonsynonymous to synonymous ratio, dn/ds=1.58) and an average entropy score of 0.372 in protein level. Interestingly, over the years (1998-2021), a decreasing trend of dn/ds (1.68 to 1.47) and an increasing trend of entropy (0.309 to 0.399) was observed. The predicted mutational frequency against Vif consensus sequence decreased over time (slope = -0.00024, p < 0.0001). Sequence conservation was observed in Vif functional motifs F1, F2, F3, G, BC box, and CBF β binding region, while variability was observed mainly in N- and C- terminal and Zinc finger region, which were dominantly under immune pressure by host HLA-I-restricted CD8+ T cell. Computational analysis of ∆∆G^stability through protein stability prediction tools suggested that missense mutation may affect Vif stability, especially in the Vif-A3 binding interface. Notably, mutations R17K and Y44F in F1 and G box were predicted to destabilize the Vif-A3 binding interface by altering bond formations with adjacent amino acids. Therefore, our analysis demonstrates Vif adaptation with host physiology by maintaining sequence conservation, especially in A3 interacting functional motifs, highlighting important therapeutic candidate regions of Vif against HIV-1 infections.