Unveiling the Impact of Drug-Sensitive Mutations on HIV-1 Protease Dynamics: A Molecular Dynamics Simulation Study of the T12A, L63Q, and H69N Variants

HIV-1 protease (PR) is an essential enzyme in the viral life cycle and a primary target of antiretroviral therapies, particularly protease inhibitors (PIs). Understanding the dynamics of viral evolution and the factors governing the emergence or loss of resistance-associated mutations is critical for improving PI efficacy and managing drug resistance in HIV/AIDS treatment. In this study, we investigated the impact of three natural HIV-1 polymorphisms (T12A, L63Q, and H69N), whose prevalence varies depending on treatment status and viral subtype, on the structural stability and conformational dynamics of PR using molecular dynamics (MD) simulations. Three independent 500 ns MD simulations were performed for the native protease and each mutant system. Although none of the mutations disrupts the overall structural integrity of HIV-1 PR, they induce mutation-specific alterations in flexibility and residue interactions. In particular, T12A and H69N exhibit increased structural deviations, especially in the flap regions, along with enhanced conformational fluctuations. In contrast, the L63Q mutation shows a slight reduction in flap flexibility compared to both the native protease and the other mutants. Consistently, the fraction of time spent in open-flap conformations is higher for T12A and H69N and lower for L63Q relative to the native system. Moreover, mutations in the Fulcrum (T12A) and Cantilever (L63Q and H69N) regions do not disrupt the long-range network of correlated motions observed in the native protease, both inter- and intra-monomer, but instead increase the extent of correlated and anti-correlated motions in other regions of PR.