Double trouble? Gag in conjunction with double insert in HIV protease contributes to reduced DRV susceptibility

HIV Protease Hinge Region Insertions at Codon 38 Affect Enzyme Kinetics, Conformational Stability and Dynamics

HIV-1 protease is essential for the production of mature, infectious virions and is a major target in antiretroviral therapy. We successfully purified a HIV-1 subtype C variant, L38↑N↑L- 4, containing an insertion of asparagine and leucine at position 38 without the four background mutations - K20R, E35D, R57K, V82I using a modified purification protocol. Isothermal titration calorimetry indicated that 50% of the variant protease sample was in the active conformation compared to 62% of the wild type protease. The secondary structure composition of the variant protease was unaffected by the double insertion. The specific activity and kcat values of the variant protease were approximately 50% lower than the wild type protease values. The variant protease also exhibited a 1.6-fold increase in kcat/KM when compared to the wild type protease. Differential scanning calorimetry showed a 5 °C increase in Tm of the variant protease, indicating the variant was more stable than the wild type. Molecular dynamics simulations indicated the variant was more stable and compact than the wild type protease. A 3-4% increase in the flexibility of the hinge regions of the variant protease was observed. In addition, increased flexibility of the flaps, cantilever and fulcrum regions of the variant protease B chain was observed. The variant protease sampled only the closed flap conformation indicating a potential mechanism for drug resistance. The present study highlights the direct impact of a double amino acid insertion in hinge region on enzyme kinetics, conformational stability and dynamics of an HIV-1 subtype C variant protease.

Understanding the co-evolutionary molecular mechanisms of resistance in the HIV-1 Gag and protease

Due to high human immunodeficiency virus type 1 (HIV-1) subtype C infections coupled with increasing antiretroviral treatment failure, the elucidation of complex drug resistance mutational patterns arising through protein co-evolution is required. Despite the inclusion of potent protease inhibitors Lopinavir (LPV) and Darunavir (DRV) in second- and third-line therapies, many patients still fail treatment due to the accumulation of mutations in protease (PR) and recently, Gag. To understand the co-evolutionary molecular mechanisms of resistance in the HIV-1 PR and Gag, we performed 100 ns molecular dynamic simulations on multidrug resistant PR's when bound to LPV, DRV or a mutated A431V NC|p1 Gag cleavage site (CS). Here we showed that distinct changes in PR's active site, flap and elbow regions due to several PR resistance mutations (L10F, M46I, I54V, L76V, V82A) were found to alter LPV and DRV drug binding. However, binding was significantly exacerbated when the mutant PRs were bound to the NC|p1 Gag CS. Although A431V was shown to coordinate several residues in PR, the L76V PR mutation was found to have a significant role in substrate recognition. Consequently, a greater binding affinity was observed when the mutated substrate was bound to an L76V-inclusive PR mutant (G_bind: -62.46 ± 5.75 kcal/mol) than without (Gbind: -50.34 ± 6.28 kcal/mol). These data showed that the co-selection of resistance mutations in the enzyme and substrate can simultaneously constrict regions in PR's active site whilst flexing the flaps to allow flexible movement of the substrate and multiple, complex mechanisms of resistance to occur. Communicated by Ramaswamy H. Sarma.

Elasticity-Associated Functionality and Inhibition of the HIV Protease

HIV protease plays a critical role in the life cycle of the virus through the generation of mature and infectious virions. Detailed knowledge of the structure of the enzyme and its substrate has led to the development of protease inhibitors. However, the development of resistance to all currently available protease inhibitors has contributed greatly to the decreased success of antiretroviral therapy. When therapy failure occurs, multiple mutations are found within the protease sequence starting with primary mutations, which directly impact inhibitor binding, which can also negatively impact viral fitness and replicative capacity by decreasing the binding affinity of the natural substrates to the protease. As such, secondary mutations which are located outside of the active site region accumulate to compensate for the recurrently deleterious effects of primary mutations. However, the resistance mechanism of these secondary mutations is not well understood, but what is known is that these secondary mutations contribute to resistance in one of two ways, either through increasing the energetic penalty associated with bringing the protease into the closed conformation, or, through decreasing the stability of the protein/drug complex in a manner that increases the dissociation rate of the drug, leading to diminished inhibition. As a result, the elasticity of the enzyme-substrate complex has been implicated in the successful recognition and catalysis of the substrates which may be inferred to suggest that the elasticity of the enzyme/drug complex plays a role in resistance. A realistic representation of the dynamic nature of the protease may provide a more powerful tool in structure-based drug design algorithms.

Acquired HIV-1 Protease Conformational Flexibility Associated with Lopinavir Failure May Shape the Outcome of Darunavir Therapy after Antiretroviral Therapy Switch

Understanding the underlying molecular interaction during a therapy switch from lopinavir (LPV) to darunavir (DRV) is essential to achieve long-term virological suppression. We investigated the kinetic and structural characteristics of multidrug-resistant South African HIV-1 subtype C protease (HIV-1 PR) during therapy switch from LPV to DRV using enzyme activity and inhibition assay, fluorescence spectroscopy, and molecular dynamic simulation. The HIV-1 protease variants were from clinical isolates with a combination of drug resistance mutations; MUT-1 (M46I, I54V, V82A, and L10F), MUT-2 (M46I, I54V, L76V, V82A, L10F, and L33F), and MUT-3 (M46I, I54V, L76V, V82A, L90M, and F53L). Enzyme kinetics analysis shows an association between increased relative resistance to LPV and DRV with the progressive decrease in the mutant HIV-1 PR variants' catalytic efficiency. A direct relationship between high-level resistance to LPV and intermediate resistance to DRV with intrinsic changes in the three-dimensional structure of the mutant HIV-1 PR as a function of the multidrug-resistance mutation was observed. In silico analysis attributed these structural adjustments to the multidrug-resistance mutations affecting the LPV and DRV binding landscape. Though DRV showed superiority to LPV, as a lower concentration was needed to inhibit the HIV-1 PR variants, the inherent structural changes resulting from mutations selected during LPV therapy may dynamically shape the DRV treatment outcome after the therapy switch.

Concerted hydrolysis mechanism of HIV-1 natural substrate against subtypes B and C-SA PR: insight through molecular dynamics and hybrid QM/MM studies

It is well known that understanding the catalytic mechanism of HIV-1 PR is the rationale on which its inhibitors were developed; therefore, a better understanding of the mechanism of natural substrate hydrolysis is important. Herein, the reaction mechanism of HIV-1 natural substrates with subtypes B and common mutant in South Africa (subtype C-SA) protease were studied through transition state modelling, using a general acid-general base (GA-GB) one-step concerted process. The activation free energies of enzyme-substrate complexes were compared based on their rate of hydrolysis using a two-layered ONIOM (B3LYP/6-31++G(d,p):AMBER) method. We expanded our computational model to obtain a better understanding of the mechanism of hydrolysis as well as how the enzyme recognises or chooses the cleavage site of the scissile bonds. Using this model, a potential substrate-based inhibitor could be developed with better potency. The calculated activation energies of natural substrates in our previous study correlated well with experimental data. A similar trend was observed for the Gag and Gag-Pol natural substrates in the present work for both enzyme complexes except for the PR-RT substrate. Natural bond orbital (NBO) analysis was also applied to determine the extent of charge transfer within the QM part of both enzymes considered and the PR-RT natural substrate. The result of this study shows that the method can be utilized as a dependable computational technique to rationalize lead compounds against specific targets.