Cantilever-centric mechanism of cooperative non-active site mutations in HIV protease: Implications for flap dynamics

Mechanism of drug resistance in HIV-1 protease subtype C in the presence of Atazanavir

AIDS is one of the deadliest diseases in the history of humankind caused by HIV. Despite the technological development, curtailing the viral infection inside human host still remains a challenge. Therapies such as HAART uses a combination of drugs to inhibit the viral activity. One of the important targets includes HIV protease and inhibiting its activity will minimize the production of mature structural proteins. However, the genetic diversity and the occurrence of drug resistant mutations adds complexity to effective drug design. In this study, we aimed at understanding the drug binding mechanism of one such subtype, namely subtype C and its insertion variant L38HL. We performed multiple molecular dynamics simulations along with binding free energy analysis of wild-type and L38HL bound to Atazanavir (ATV). From the analysis, we revealed that the insertion alters the hydrogen bond and hydrophobic interaction networks. The alterations in the interaction networks increase flexibility at the hinge-fulcrum interface. Further, the effects of these changes affect flap tip curling. Moreover, the changes in the hinge-fulcrum-cantilever interface alters the concerted motion of the functional regions leading to change in the direction of flap movement thus causing a subtle change in the active site volume. Additionally, formation of intramolecular hydrogen bonds in the ATV docked to L38HL restricted the movement of R1 and R2 groups thereby altering the interactions. Overall, the changes in the flexibility of flap together with the changes in the active site volume and compactness of the ligand provide insights for increased binding affinity of ATV with L38HL.

FMO-guided design of darunavir analogs as HIV-1 protease inhibitors

The prevalence of HIV-1 infection continues to pose a significant global public health issue, highlighting the need for antiretroviral drugs that target viral proteins to reduce viral replication. One such target is HIV-1 protease (PR), responsible for cleaving viral polyproteins, leading to the maturation of viral proteins. While darunavir (DRV) is a potent HIV-1 PR inhibitor, drug resistance can arise due to mutations in HIV-1 PR. To address this issue, we developed a novel approach using the fragment molecular orbital (FMO) method and structure-based drug design to create DRV analogs. Using combinatorial programming, we generated novel analogs freely accessible via an on-the-cloud mode implemented in Google Colab, Combined Analog generator Tool (CAT). The designed analogs underwent cascade screening through molecular docking with HIV-1 PR wild-type and major mutations at the active site. Molecular dynamics (MD) simulations confirmed the assess ligand binding and susceptibility of screened designed analogs. Our findings indicate that the three designed analogs guided by FMO, 19-0-14-3, 19-8-10-0, and 19-8-14-3, are superior to DRV and have the potential to serve as efficient PR inhibitors. These findings demonstrate the effectiveness of our approach and its potential to be used in further studies for developing new antiretroviral drugs.

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 Drug Resistance of Wild-Type and L38HL Insertion Mutant of HIV-1 C Protease to Saquinavir

Acquired immunodeficiency syndrome (AIDS) is one of the most challenging infectious diseases to treat on a global scale. Understanding the mechanisms underlying the development of drug resistance is necessary for novel therapeutics. HIV subtype C is known to harbor mutations at critical positions of HIV aspartic protease compared to HIV subtype B, which affects the binding affinity. Recently, a novel double-insertion mutation at codon 38 (L38HL) was characterized in HIV subtype C protease, whose effects on the interaction with protease inhibitors are hitherto unknown. In this study, the potential of L38HL double-insertion in HIV subtype C protease to induce a drug resistance phenotype towards the protease inhibitor, Saquinavir (SQV), was probed using various computational techniques, such as molecular dynamics simulations, binding free energy calculations, local conformational changes and principal component analysis. The results indicate that the L38HL mutation exhibits an increase in flexibility at the hinge and flap regions with a decrease in the binding affinity of SQV in comparison with wild-type HIV protease C. Further, we observed a wide opening at the binding site in the L38HL variant due to an alteration in flap dynamics, leading to a decrease in interactions with the binding site of the mutant protease. It is supported by an altered direction of motion of flap residues in the L38HL variant compared with the wild-type. These results provide deep insights into understanding the potential drug resistance phenotype in infected individuals.

Non-active site mutations in the HIV protease: Diminished drug binding affinity is achieved through modulating the hydrophobic sliding mechanism

The global HIV/AIDS epidemic still currently affects approximately 38 million individuals globally. The protease enzyme of the human immunodeficiency virus is a major drug target in antiviral therapy, however, under the influence of reverse transcriptase and in the context of drug pressure, the rapid PR mutation rate contributes significantly to clinical failure. The set of cooperative non-active site mutations, I13V/I62V/V77I, have been associated with reduced inhibitor susceptibility and are the focus of the current study. When compared to the wild-type protease the mutant protease exhibited decreased binding affinities towards ATV and DRV by 64- and 12-fold, respectively, and decreased the overall favourable Gibbs free energy for ATV, DRV, RTV and SQV. Moreover, these mutations decreased the thermal stability of the protease when in complex with ATV and DRV by approximately 6.4 and 4.2 °C, respectively. The crystal structure of the mutant protease revealed that the location of these mutations and their effect on the hydrophobic sliding mechanism may be crucial in their role in resistance.

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.