Prediction and validation of HIV-1 gp41 ecto-transmembrane domain post-fusion trimeric structure using molecular modeling

The glycoproteins on the surface of human immunodeficiency virus (HIV) undergoes cascade of conformational transitions to evade the human immune system. The virus replicates inside the host and infects the T-cells instigating acquired immunodeficiency syndrome (AIDS). The glycoprotein 41 (gp41) of HIV helps to mediate the fusion of virus and host membranes. The detailed mechanism of host cell invasion by virus remains obscure due to the unavailability of experimental structure of complete gp41. In the current study, the post-fusion (PoF) trimeric structure of ecto-domain including transmembrane domain of gp41 was modeled using multiple homologous templates of Simian immunodeficiency virus (SIV) and HIV-1. In order to validate the gp41 model, interactions of three peptide inhibitors: T20, C37 and C34; were studied using all-atom molecular dynamics (MD) simulations, binding free-energy calculation and per-residue energy decomposition analysis. The binding free energy calculated using MM-PBSA (Molecular Mechanics Poisson-Boltzmann surface area) method predicts maximum affinity for C34 and minimum by T20 for gp41, which is in good agreement with the available computational and experimental studies. The van der Waals interaction is a dominant contributor for the peptide-gp41 complexes. The per-residue decomposition of energy confirmed the role of Trp117, Trp120 and Ile124, present in C34 and C37, for the strong hydrophobic interactions with the deep pocket localized around the N-terminal of gp41, which is lacking in T20. The HIV-1 gp41 structure developed in this work can be used in future study to gain insight into the mechanism of virus invasion and probing potent inhibitor to eliminate AIDS