Crystal structures of HIV protease V82A and L90M mutants reveal changes in the indinavir-binding site

Hybrid QM/MM Free-Energy Evaluation of Drug-Resistant Mutational Effect on the Binding of an Inhibitor Indinavir to HIV-1 Protease

A human immunodeficiency virus-1 (HIV-1) protease is a homodimeric aspartic protease essential for the replication of HIV. The HIV-1 protease is a target protein in drug discovery for antiretroviral therapy, and various inhibitor molecules of transition state analogues have been developed. However, serious drug-resistant mutants have emerged. For understanding the molecular mechanism of the drug resistance, an accurate examination of the impacts of the mutations on ligand binding and enzymatic activity is necessary. Here, we present a molecular simulation study on the ligand binding of indinavir, a potent transition state analogue inhibitor, to the wild-type protein and a V82T/I84V drug-resistant mutant of the HIV-1 protease. We employed a hybrid ab initio quantum mechanical/molecular mechanical (QM/MM) free-energy optimization technique which combines a highly accurate QM description of the ligand molecule and its interaction with statistically ample conformational sampling of the MM protein environment by long-time molecular dynamics simulations. Through the free-energy calculations of protonation states of catalytic groups at the binding pocket and of the ligand-binding affinity changes upon the mutations, we successfully reproduced the experimentally observed significant reduction of the binding affinity upon the drug-resistant mutations and elucidated the underlying molecular mechanism. The present study opens the way for understanding the molecular mechanism of drug resistance through the direct quantitative comparison of ligand binding and enzymatic reaction with the same accuracy.

Comparative analysis of the unbinding pathways of antiviral drug Indinavir from HIV and HTLV1 proteases by supervised molecular dynamics simulation

Determining the unbinding pathways of potential small molecule compounds from their target proteins is of great significance for designing efficacious treatment solutions. One of these potential compounds is the approved HIV-1 protease inhibitor, Indinavir, which has a weak effect on the HTLV-1 protease. In this work, by employing the SuMD method, we reconstructed the unbinding pathways of Indinavir from HIV and HTLV-1 proteases to compare and understand the mechanism of the unbinding and to discover the reasons for the lack of inhibitory activity of Indinavir against the HTLV-1 protease. We achieved multiple unbinding events from both HIV and HTLV-1 proteases in which the RMSD values of Indinavir reached over 40 Å. Also, we found that the mobility and fluctuations of the flap region are higher in the HTLV-1 protease, making the drug less stable. We realized that critically positioned aromatic residues such as Trp98/Trp98' and Phe67/Phe67' in the HTLV-1 protease could make strong π-Stacking interactions with Indinavir in the unbinding pathway, which are unfavorable for the stability of Indinavir in the active site. The details found in this study can make a reasonable explanation for the lack of inhibitory activity of this drug against HTLV-1 protease. We believe the details discovered in this work can help design more effective and selective inhibitors for the HTLV-1 protease.

Prediction of HIV drug resistance based on the 3D protein structure: Proposal of molecular field mapping

A method for predicting HIV drug resistance by using genotypes would greatly assist in selecting appropriate combinations of antiviral drugs. Models reported previously have had two major problems: lack of information on the 3D protein structure and processing of incomplete sequencing data in the modeling procedure. We propose obtaining the 3D structural information of viral proteins by using homology modeling and molecular field mapping, instead of just their primary amino acid sequences. The molecular field potential parameters reflect the physicochemical characteristics associated with the 3D structure of the proteins. We also introduce the Bayesian conditional mutual information theory to estimate the probabilities of occurrence of all possible protein candidates from an incomplete sequencing sample. This approach allows for the effective use of uncertain information for the modeling process. We applied these data analysis techniques to the HIV-1 protease inhibitor dataset and developed drug resistance prediction models with reasonable performance.

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.

Structural Impacts of Drug-Resistance Mutations Appearing in HIV-2 Protease

The use of antiretroviral drugs is accompanied by the emergence of HIV-2 resistances. Thus, it is important to elucidate the mechanisms of resistance to antiretroviral drugs. Here, we propose a structural analysis of 31 drug-resistant mutants of HIV-2 protease (PR2) that is an important target against HIV-2 infection. First, we modeled the structures of each mutant. We then located structural shifts putatively induced by mutations. Finally, we compared wild-type and mutant inhibitor-binding pockets and interfaces to explore the impacts of these induced structural deformations on these two regions. Our results showed that one mutation could induce large structural rearrangements in side-chain and backbone atoms of mutated residue, in its vicinity or further. Structural deformations observed in side-chain atoms are frequent and of greater magnitude, that confirms that to fight drug resistance, interactions with backbone atoms should be favored. We showed that these observed structural deformations modify the conformation, volume, and hydrophobicity of the binding pocket and the composition and size of the PR2 interface. These results suggest that resistance mutations could alter ligand binding by modifying pocket properties and PR2 stability by impacting its interface. Our results reinforce the understanding of the effects of mutations that occurred in PR2 and the different mechanisms of PR2 resistance.

An Overview of Antiretroviral Agents for Treating HIV Infection in Paediatric Population

Paediatric Acquired ImmunoDeficiency Syndrome (AIDS) is a life-threatening and infectious disease in which the Human Immunodeficiency Virus (HIV) is mainly transmitted through Mother-To- Child Transmission (MTCT) during pregnancy, labour and delivery, or breastfeeding. This review provides an overview of the distinct therapeutic alternatives to abolish the systemic viral replication in paediatric HIV-1 infection. Numerous classes of antiretroviral agents have emerged as therapeutic tools for downregulation of different steps in the HIV replication process. These classes encompass Non- Nucleoside Analogue Reverse Transcriptase Inhibitors (NNRTIs), Nucleoside/Nucleotide Analogue Reverse Transcriptase Inhibitors (NRTIs/NtRTIs), INtegrase Inhibitors (INIs), Protease Inhibitors (PIs), and Entry Inhibitors (EIs). Co-administration of certain antiretroviral drugs with Pharmacokinetic Enhancers (PEs) may boost the effectiveness of the primary therapeutic agent. The combination of multiple antiretroviral drug regimens (Highly Active AntiRetroviral Therapy - HAART) is currently the standard therapeutic approach for HIV infection. So far, the use of HAART offers the best opportunity for prolonged and maximal viral suppression, and preservation of the immune system upon HIV infection. Still, the frequent administration of high doses of multiple drugs, their inefficient ability to reach the viral reservoirs in adequate doses, the development of drug resistance, and the lack of patient compliance compromise the complete HIV elimination. The development of nanotechnology-based drug delivery systems may enable targeted delivery of antiretroviral agents to inaccessible viral reservoir sites at therapeutic concentrations. In addition, the application of Computer-Aided Drug Design (CADD) approaches has provided valuable tools for the development of anti-HIV drug candidates with favourable pharmacodynamics and pharmacokinetic properties.

Novel radial distribution function approach in the study of point mutations: the HIV-1 protease case study

Background: Mutations are one of the engines of evolution. Under constant stress pressure, mutations can lead to the emergence of unwanted, drug-resistant entities. Methodology: The radial distribution function weighted by the number of valence shell electrons is used to design quantitative structure-activity relationship (QSAR) model relating descriptors with the inhibition constant for a series of wild-type HIV-1 protease inhibitor complexes. The residuals of complexes with mutant HIV-1 protease were correlated with the energy of the highest occupied molecular orbitals of the residues introduced to enzyme via point mutations. Conclusion: Successful identification of residues Ile3, Asp25, Val32 and Ile50 as the one whose substitution influences the inhibition constant the most, demonstrates the potential of the proposed methodology for the study of the effects of point mutations.

Halogen Bond Interactions of Novel HIV-1 Protease Inhibitors (PI) (GRL-001-15 and GRL-003-15) with the Flap of Protease Are Critical for Their Potent Activity against Wild-Type HIV-1 and Multi-PI-Resistant Variants

We generated two novel nonpeptidic HIV-1 protease inhibitors (PIs), GRL-001-15 and GRL-003-15, which contain unique crown-like tetrahydropyranofuran (Crn-THF) and P2'-cyclopropyl-aminobenzothiazole (Cp-Abt) moieties as P2 and P2' ligands, respectively. GRL-001-15 and GRL-003-15 have meta-monofluorophenyl and para-monofluorophenyl at the P1 site, respectively, exert highly potent activity against wild-type HIV-1 with 50% effective concentrations (EC50s) of 57 and 50 pM, respectively, and have favorable cytotoxicity profiles with 50% cytotoxic concentrations (CC50s) of 38 and 11 μM, respectively. The activity of GRL-001-15 against multi-PI-resistant HIV-1 variants was generally greater than that of GRL-003-15. The EC50 of GRL-001-15 against an HIV-1 variant that was highly resistant to multiple PIs, including darunavir (DRV) (HIV-1DRVRP30), was 0.17 nM, and that of GRL-003-15 was 3.3 nM, while DRV was much less active, with an EC50 of 216 nM. The emergence of HIV-1 variants resistant to GRL-001-15 and GRL-003-15 was significantly delayed compared to that of variants resistant to selected PIs, including DRV. Structural analyses of wild-type protease (PRWT) complexed with the novel PIs revealed that GRL-001-15's meta-fluorine atom forms halogen bond interactions (2.9 and 3.0 Å) with Gly49 and Ile50, respectively, of the protease flap region and with Pro81' (2.7 and 3.2 Å), which is located close to the protease active site, and that two fluorine atoms of GRL-142-13 form multiple halogen bond interactions with Gly49, Ile50, Pro81', Ile82', and Arg8'. In contrast, GRL-003-15 forms halogen bond interactions with Pro81' alone, suggesting that the reduced antiviral activity of GRL-003-15 is due to the loss of the interactions with the flap region.

Highly Drug-Resistant HIV-1 Protease Mutant PRS17 Shows Enhanced Binding to Substrate Analogues

We report the structural analysis of highly drug-resistant human immunodeficiency virus protease (PR) variant PR^S17, rationally selected by machine learning, in complex with substrate analogues. Crystal structures were solved of inhibitor-free inactive PR^S17-D25N, wild-type PR/CA-p2 complex, and PR^S17 in complex with substrate analogues, CA-p2 and p2-NC. Peptide analogues p2-NC and CA-p2 exhibit inhibition constants of 514 and 22 nM, respectively, for PR^S17 or approximately 3-fold better than for PR. CA-p2 is a better inhibitor of PR^S17 than are clinical inhibitors (K _i = 50-8390 nM) except for amprenavir (K _i = 11 nM). G48V resistance mutation induces curled flap tips in PR^S17-D25N structure. The inner P2-P2' residues of substrate analogues in PR^S17 complexes maintain similar conformations to those of wild-type complex, while significant conformational changes are observed in the peripheral residues P3, P4' of CA-p2 and P3, P4, and P3' of p2-NC. The loss of β-branched side chain by V82S mutation initiates a shift in 80's loop and reshapes the S3/S3' subsite, which enhances substrate binding with new hydrogen bonds and van der Waals interactions that are absent in the wild-type structures. The steric hindrance caused by G48V mutation in the flap of PR^S17 contributes to altered binding interactions of P3 Arg, P4' norleucine of CA-p2, and P4 and P3' of p2-NC with the addition of new hydrogen bonds and van der Waals contacts. The enhanced interaction of PR^S17 with substrate analogues agrees with their relative inhibition, suggesting that this mutant improves substrate binding while decreasing affinity for clinical inhibitors.

Inference of Epistatic Effects Leading to Entrenchment and Drug Resistance in HIV-1 Protease

Understanding the complex mutation patterns that give rise to drug resistant viral strains provides a foundation for developing more effective treatment strategies for HIV/AIDS. Multiple sequence alignments of drug-experienced HIV-1 protease sequences contain networks of many pair correlations which can be used to build a (Potts) Hamiltonian model of these mutation patterns. Using this Hamiltonian model, we translate HIV-1 protease sequence covariation data into quantitative predictions for the probability of observing specific mutation patterns which are in agreement with the observed sequence statistics. We find that the statistical energies of the Potts model are correlated with the fitness of individual proteins containing therapy-associated mutations as estimated by in vitro measurements of protein stability and viral infectivity. We show that the penalty for acquiring primary resistance mutations depends on the epistatic interactions with the sequence background. Primary mutations which lead to drug resistance can become highly advantageous (or entrenched) by the complex mutation patterns which arise in response to drug therapy despite being destabilizing in the wildtype background. Anticipating epistatic effects is important for the design of future protease inhibitor therapies.

Structural Studies of a Rationally Selected Multi-Drug Resistant HIV-1 Protease Reveal Synergistic Effect of Distal Mutations on Flap Dynamics

We report structural analysis of HIV protease variant PR^S17 which was rationally selected by machine learning to represent wide classes of highly drug-resistant variants. Crystal structures were solved of PR^S17 in the inhibitor-free form and in complex with antiviral inhibitor, darunavir. Despite its 17 mutations, PR^S17 has only one mutation (V82S) in the inhibitor/substrate binding cavity, yet exhibits high resistance to all clinical inhibitors. PR^S17 has none of the major mutations (I47V, I50V, I54ML, L76V and I84V) associated with darunavir resistance, but has 10,000-fold weaker binding affinity relative to the wild type PR. Comparable binding affinity of 8000-fold weaker than PR is seen for drug resistant mutant PR20, which bears 3 mutations associated with major resistance to darunavir (I47V, I54L and I84V). Inhibitor-free PR^S17 shows an open flap conformation with a curled tip correlating with G48V flap mutation. NMR studies on inactive PR^S17_D25N unambiguously confirm that the flaps adopt mainly an open conformation in solution very similar to that in the inhibitor-free crystal structure. In PR^S17, the hinge loop cluster of mutations, E35D, M36I and S37D, contributes to the altered flap dynamics by a mechanism similar to that of PR20. An additional K20R mutation anchors an altered conformation of the hinge loop. Flap mutations M46L and G48V in PR^S17/DRV complex alter the Phe53 conformation by steric hindrance between the side chains. Unlike the L10F mutation in PR20, L10I in PR^S17 does not break the inter-subunit ion pair or diminish the dimer stability, consistent with a very low dimer dissociation constant comparable to that of wild type PR. Distal mutations A71V, L90M and I93L propagate alterations to the catalytic site of PR^S17. PR^S17 exhibits a molecular mechanism whereby mutations act synergistically to alter the flap dynamics resulting in significantly weaker binding yet maintaining active site contacts with darunavir.

Recent Progress in the Development of HIV-1 Protease Inhibitors for the Treatment of HIV/AIDS

HIV-1 protease inhibitors continue to play an important role in the treatment of HIV/AIDS, transforming this deadly ailment into a more manageable chronic infection. Over the years, intensive research has led to a variety of approved protease inhibitors for the treatment of HIV/AIDS. In this review, we outline current drug design and medicinal chemistry efforts toward the development of next-generation protease inhibitors beyond the currently approved drugs.

Computational Studies of a Mechanism for Binding and Drug Resistance in the Wild Type and Four Mutations of HIV-1 Protease with a GRL-0519 Inhibitor

Drug resistance of mutations in HIV-1 protease (PR) is the most severe challenge to the long-term efficacy of HIV-1 PR inhibitor in highly active antiretroviral therapy. To elucidate the molecular mechanism of drug resistance associated with mutations (D30N, I50V, I54M, and V82A) and inhibitor (GRL-0519) complexes, we have performed five molecular dynamics (MD) simulations and calculated the binding free energies using the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method. The ranking of calculated binding free energies is in accordance with the experimental data. The free energy spectra of each residue and inhibitor interaction for all complexes show a similar binding model. Analysis based on the MD trajectories and contribution of each residues show that groups R2 and R3 mainly contribute van der Waals energies, while groups R1 and R4 contribute electrostatic interaction by hydrogen bonds. The drug resistance of D30N can be attributed to the decline in binding affinity of residues 28 and 29. The size of Val50 is smaller than Ile50 causes the residue to move, especially in chain A. The stable hydrophobic core, including the side chain of Ile54 in the wild type (WT) complex, became unstable in I54M because the side chain of Met54 is flexible with two alternative conformations. The binding affinity of Ala82 in V82A decreases relative to Val82 in WT. The present study could provide important guidance for the design of a potent new drug resisting the mutation inhibitors.

Molecular Dynamics Simulations of HIV-1 Protease Suggest Different Mechanisms Contributing to Drug Resistance

A major problem in the antiretroviral treatment of HIV-infections with protease-inhibitors is the emergence of resistance, resulting from the occurrence of distinct mutations within the protease molecule. In the present work molecular dynamics simulations of an active-site mutation (D30N) and a nonactive-site mutation (N88S) of HIV-1 protease that both directly confer resistance to the protease inhibitor Nelfinavir but not to Amprenavir were performed and compared to wild-type HIV-protease. A decreased interaction energy between protease and Nelfinavir was observed for the D30N mutant giving a plausible explanation for resistance, while the N88S mutation did not significantly affect the interaction energies in the bound form. Structural analysis including both ligand-bound and unliganded HIV-1 proteases revealed that the free N88S mutant protease shows significant differences in its hydrogen bonding pattern compared to free or Nelfinavir-bound wild-type protease. In particular, Asp30 forms more frequently a hydrogen bond with Ser88 in the unbound N88S mutant thus interfering with the Asp30-Nelfinavir interaction. These findings suggest that different molecular mechanisms contribute to resistance in active-site and nonactive-site mutants and propose a mechanism for the N88S mutant that is based on a shift of the conformational equilibrium of the unbound protease.

Structural basis for HTLV-1 protease inhibition by the HIV-1 protease inhibitor indinavir

HTLV-1 protease (HTLV-1 PR) is an aspartic protease which represents a promising drug target for the discovery of novel anti-HTLV-1 drugs. The X-ray structure of HTLV-1 PR in complex with the well-known and approved HIV-1 PR inhibitor Indinavir was determined at 2.40 Å resolution. In this contribution, we describe the first crystal structure in complex with a nonpeptidic inhibitor that accounts for rationalizing the rather moderate affinity of Indinavir against HTLV-1 PR and provides the basis for further structure-guided optimization strategies.

How conformational changes can affect catalysis, inhibition and drug resistance of enzymes with induced-fit binding mechanism such as the HIV-1 protease

A central question is how the conformational changes of proteins affect their function and the inhibition of this function by drug molecules. Many enzymes change from an open to a closed conformation upon binding of substrate or inhibitor molecules. These conformational changes have been suggested to follow an induced-fit mechanism in which the molecules first bind in the open conformation in those cases where binding in the closed conformation appears to be sterically obstructed such as for the HIV-1 protease. In this article, we present a general model for the catalysis and inhibition of enzymes with induced-fit binding mechanism. We derive general expressions that specify how the overall catalytic rate of the enzymes depends on the rates for binding, for the conformational changes, and for the chemical reaction. Based on these expressions, we analyze the effect of mutations that mainly shift the conformational equilibrium on catalysis and inhibition. If the overall catalytic rate is limited by product unbinding, we find that mutations that destabilize the closed conformation relative to the open conformation increase the catalytic rate in the presence of inhibitors by a factor exp(ΔΔGC/RT) where ΔΔGC is the mutation-induced shift of the free-energy difference between the conformations. This increase in the catalytic rate due to changes in the conformational equilibrium is independent of the inhibitor molecule and, thus, may help to understand how non-active-site mutations can contribute to the multi-drug-resistance that has been observed for the HIV-1 protease. A comparison to experimental data for the non-active-site mutation L90M of the HIV-1 protease indicates that the mutation slightly destabilizes the closed conformation of the enzyme. This article is part of a Special Issue entitled: The emerging dynamic view of proteins: Protein plasticity in allostery, evolution and self-assembly.

Novel P2 tris-tetrahydrofuran group in antiviral compound 1 (GRL-0519) fills the S2 binding pocket of selected mutants of HIV-1 protease

GRL-0519 (1) is a potent antiviral inhibitor of HIV-1 protease (PR) possessing tris-tetrahydrofuran (tris-THF) at P2. The high resolution X-ray crystal structures of inhibitor 1 in complexes with single substitution mutants PR(R8Q), PR(D30N), PR(I50V), PR(I54M), and PR(V82A) were analyzed in relation to kinetic data. The smaller valine side chain in PR(I50V) eliminated hydrophobic interactions with inhibitor and the other subunit consistent with 60-fold worse inhibition. Asn30 in PR(D30N) showed altered interactions with neighboring residues and 18-fold worse inhibition. Mutations V82A and I54M showed compensating structural changes consistent with 6-7-fold lower inhibition. Gln8 in PR(R8Q) replaced the ionic interactions of wild type Arg8 with hydrogen bond interactions without changing the inhibition significantly. The carbonyl oxygen of Gly48 showed two alternative conformations in all structures likely due to the snug fit of the large tris-THF group in the S2 subsite in agreement with high antiviral efficacy of 1 on resistant virus.

Binding of single walled carbon nanotube to WT and mutant HIV-1 proteases: analysis of flap dynamics and binding mechanism

Most of the currently treated HIV-1 protease (HIV-PR) inhibitors have been prone to suffer from the mutations associated drug resistance. Therefore, it is necessary to search for potent alternatives against the drug resistance. In the current study we have tested the single-walled carbon nanotube (SWCNT) as an inhibitor in wild type (WT) as well as in three primary mutants (I50V(PR), V82A(PR) and I84V(PR)) of the HIV-1-PR through docking the SWCNT in the active site region, and then performed all-atom MD simulations for the complexes. The conformational dynamics of HIV-PR with a 20 ns trajectory reveals that the SWCNT can effectively bind to the HIV-1-PR active site and regulate the flap dynamics such as maintaining the flap-flap closed. To gain an insight into the binding affinity, we also performed the MM-PBSA based binding free energy calculations for the four HIV-PR/SWCNT complexes. It was observed that, although the binding between the SWCNT and the HIV-PR decreases due to the mutations, the SWCNTs bind to the HIV-PRs 3-5 folds stronger than the most potent HIV-1-PR inhibitor, TMC114. Remarkably, the significant interactions with binding energy higher than 1kcal/mol focus on the flap and active regions, which favors closing flap-flap and deactivating the active residues of the HIV-PR. The flap dynamics and binding strength information for HIV-PR and SWCNTs can help design SWCNT-based HIV-1-PR inhibitors.

Potent antiviral HIV-1 protease inhibitor GRL-02031 adapts to the structures of drug resistant mutants with its P1'-pyrrolidinone ring

GRL-02031 (1) is an HIV-1 protease (PR) inhibitor containing a novel P1' (R)-aminomethyl-2-pyrrolidinone group. Crystal structures at resolutions of 1.25-1.55 Å were analyzed for complexes of 1 with the PR containing major drug resistant mutations, PR(I47V), PR(L76V), PR(V82A), and PR(N88D). Mutations of I47V and V82A alter residues in the inhibitor-binding site, while L76V and N88D are distal mutations having no direct contact with the inhibitor. Substitution of a smaller amino acid in PR(I47V) and PR(L76V) and the altered charge of PR(N88D) are associated with significant local structural changes compared to the wild-type PR(WT), while substitution of alanine in PR(V82A) increases the size of the S1' subsite. The P1' pyrrolidinone group of 1 accommodates to these local changes by assuming two different conformations. Overall, the conformation and interactions of 1 with PR mutants resemble those of PR(WT) with similar inhibition constants in good agreement with the antiviral potency on multidrug resistant HIV-1.

HIV-1 protease with 20 mutations exhibits extreme resistance to clinical inhibitors through coordinated structural rearrangements

The escape mutant of HIV-1 protease (PR) containing 20 mutations (PR20) undergoes efficient polyprotein processing even in the presence of clinical protease inhibitors (PIs). PR20 shows >3 orders of magnitude decreased affinity for PIs darunavir (DRV) and saquinavir (SQV) relative to PR. Crystal structures of PR20 crystallized with yttrium, substrate analogue p2-NC, DRV, and SQV reveal three distinct conformations of the flexible flaps and diminished interactions with inhibitors through the combination of multiple mutations. PR20 with yttrium at the active site exhibits widely separated flaps lacking the usual intersubunit contacts seen in other inhibitor-free dimers. Mutations of residues 35-37 in the hinge loop eliminate interactions and perturb the flap conformation. Crystals of PR20/p2-NC contain one uninhibited dimer with one very open flap and one closed flap and a second inhibitor-bound dimer in the closed form showing six fewer hydrogen bonds with the substrate analogue relative to wild-type PR. PR20 complexes with PIs exhibit expanded S2/S2' pockets and fewer PI interactions arising from coordinated effects of mutations throughout the structure, in agreement with the strikingly reduced affinity. In particular, insertion of the large aromatic side chains of L10F and L33F alters intersubunit interactions and widens the PI binding site through a network of hydrophobic contacts. The two very open conformations of PR20 as well as the expanded binding site of the inhibitor-bound closed form suggest possible approaches for modifying inhibitors to target extreme drug-resistant HIV.

Enhancing protein backbone binding--a fruitful concept for combating drug-resistant HIV

The evolution of drug resistance is one of the most fundamental problems in medicine. In HIV/AIDS, the rapid emergence of drug-resistant HIV-1 variants is a major obstacle to current treatments. HIV-1 protease inhibitors are essential components of present antiretroviral therapies. However, with these protease inhibitors, resistance occurs through viral mutations that alter inhibitor binding, resulting in a loss of efficacy. This loss of potency has raised serious questions with regard to effective long-term antiretroviral therapy for HIV/AIDS. In this context, our research has focused on designing inhibitors that form extensive hydrogen-bonding interactions with the enzyme's backbone in the active site. In doing so, we limit the protease's ability to acquire drug resistance as the geometry of the catalytic site must be conserved to maintain functionality. In this Review, we examine the underlying principles of enzyme structure that support our backbone-binding concept as an effective means to combat drug resistance and highlight their application in our recent work on antiviral HIV-1 protease inhibitors.

Flexible catalytic site conformations implicated in modulation of HIV-1 protease autoprocessing reactions

BACKGROUND

The HIV-1 protease is initially synthesized as part of the Gag-Pol polyprotein in the infected cell. Protease autoprocessing, by which the protease domain embedded in the precursor catalyzes essential cleavage reactions, leads to liberation of the free mature protease at the late stage of the replication cycle. To examine autoprocessing reactions in transfected mammalian cells, we previously described an assay using a fusion precursor consisting of the mature protease (PR) along with its upstream transframe region (p6*) sandwiched between GST and a small peptide epitope.

RESULTS

In this report, we studied two autoprocessing cleavage reactions, one between p6* and PR (the proximal site) and the other in the N-terminal region of p6* (the distal site) catalyzed by the embedded protease, using our cell-based assay. A fusion precursor carrying the NL4-3 derived protease cleaved both sites, whereas a precursor with a pseudo wild type protease preferentially autoprocessed the proximal site. Mutagenesis analysis demonstrated that several residues outside the active site (Q7, L33, N37, L63, C67 and H69) contributed to the differential substrate specificity. Furthermore, the cleavage reaction at the proximal site mediated by the embedded protease in precursors carrying different protease sequences or C-terminal fusion peptides displayed varied sensitivity to inhibition by darunavir, a catalytic site inhibitor. On the other hand, polypeptides such as a GCN4 motif, GFP, or hsp70 fused to the N-terminus of p6* had a minimal effect on darunavir inhibition of either cleavage reaction.

CONCLUSIONS

Taken together, our data suggest that several non-active site residues and the C-terminal flanking peptides regulate embedded protease activity through modulation of the catalytic site conformation. The cell-based assay provides a sensitive tool to study protease autoprocessing reactions in mammalian cells.

Accessory mutations maintain stability in drug-resistant HIV-1 protease

The underlying mechanisms driving the evolution of drug resistance in human immunodeficiency virus (HIV) are only partially understood. We investigated the evolutionary cost of the major resistance mutations in HIV-1 protease in terms of protein stability. The accumulation of resistance mutations destabilizes the protease, limiting further adaptation. From an analysis of clinical isolates, we identified specific accessory mutations that were able to restore the stability of the protease or even increase it beyond the wild-type baseline. Resistance mutations were also found to decrease the activity of HIV protease near neutral pH values, while incorporating stabilizing mutations improved the enzyme's pH tolerance. These findings help us to explain the prevalence of mutations located far from the active site and underscore the importance of protein stability during the evolution of drug resistance in HIV.

The L76V drug resistance mutation decreases the dimer stability and rate of autoprocessing of HIV-1 protease by reducing internal hydrophobic contacts

The mature HIV-1 protease (PR) bearing the L76V drug resistance mutation (PR(L76V)) is significantly less stable, with a >7-fold higher dimer dissociation constant (K(d)) of 71 ± 24 nM and twice the sensitivity to urea denaturation (UC(50) = 0.85 M) relative to those of PR. Differential scanning calorimetry showed decreases in T(m) of 12 °C for PR(L76V) in the absence of inhibitors and 5-7 °C in the presence of inhibitors darunavir (DRV), saquinavir (SQV), and lopinavir (LPV), relative to that of PR. Isothermal titration calorimetry gave a ligand dissociation constant of 0.8 nM for DRV, ∼160-fold higher than that of PR, consistent with DRV resistance. Crystal structures of PR(L76V) in complexes with DRV and SQV were determined at resolutions of 1.45-1.46 Å. Compared to the corresponding PR complexes, the mutated Val76 lacks hydrophobic interactions with Asp30, Lys45, Ile47, and Thr74 and exhibits closer interactions with Val32 and Val56. The bound DRV lacks one hydrogen bond with the main chain of Asp30 in PR(L76V) relative to PR, possibly accounting for the resistance to DRV. SQV shows slightly improved polar interactions with PR(L76V) compared to those with PR. Although the L76V mutation significantly slows the N-terminal autoprocessing of the precursor TFR-PR(L76V) to give rise to the mature PR(L76V), the coselected M46I mutation counteracts the effect by enhancing this rate but renders the TFR-PR(M46I/L76V) precursor less responsive to inhibition by 6 μM LPV while preserving inhibition by SQV and DRV. The correlation of lowered stability, higher K(d), and impaired autoprocessing with reduced internal hydrophobic contacts suggests a novel molecular mechanism for drug resistance.

Effect of structural stress on the flexibility and adaptability of HIV-1 protease

Resistance remains a major issue with regards to HIV-1 protease, despite the availability of numerous HIV-1 protease inhibitors and copious amounts of structural and binding data. In an effort to improve our understanding of how HIV-1 protease is able to "outsmart" new drugs, we have investigated the flexibility of HIV-1 protease and in particular how it adapts to different structural stresses. Our analysis has highlighted the effects of space group on the variability between structures of HIV-1 protease and suggests that consideration of multiple structures and appropriate consideration of different conformations of the Ile50 residue is necessary in any structural analysis. Calculation of the root-mean-square deviation on a per-residue basis has been used to identify 'natural variation', while mutational and ligand analyses have been carried out to identify the effect on structure as a result of specific stresses. It was observed that mutations readily cause changes to occur at sites both close to and distant from a mutation site, with changes more likely to occur at residues that are sites of other major mutations. It is also revealed that HIV-1 protease adaption is dependent on the type and the structure of any bound ligand. Identification of the specific changes that occur due to these stresses will aid in the understanding of resistance and also aid in the design of new drugs.

Dynamic flaps in HIV-1 protease adopt unique ordering at different stages in the catalytic cycle

The flexibility of HIV-1 protease flaps is known to be essential for the enzymatic activity. Here we attempt to capture a multitude of conformations of the free and substrate-bound HIV-1 protease that differ drastically in their flap arrangements. The substrate binding process suggests the opening of active site gate in conjunction with a reversal of flap tip ordering, from the native semiopen state. The reversed-flap, open-gated enzyme readily transforms to a closed conformation after proper placement of the substrate into the binding cleft. After substrate processing, the closed state protease which possessed opposite flap ordering relative to the semiopen state, encounters another flap reversal via a second open conformation that facilitates the evolution of native semiopen state of correct flap ordering. The complicated transitional pathway, comprising of many high and low energy states, is explored by combining standard and activated molecular dynamics (MD) simulation techniques. The study not only complements the existing findings from X-ray, NMR, EPR, and MD studies but also provides a wealth of detailed information that could help the structure-based drug design process.

Amprenavir complexes with HIV-1 protease and its drug-resistant mutants altering hydrophobic clusters

The structural and kinetic effects of amprenavir (APV), a clinical HIV protease (PR) inhibitor, were analyzed with wild-type enzyme and mutants with single substitutions of V32I, I50V, I54V, I54M, I84V and L90M that are common in drug resistance. Crystal structures of the APV complexes at resolutions of 1.02-1.85 Å reveal the structural changes due to the mutations. Substitution of the larger side chains in PR(V32I) , PR(I54M) and PR(L90M) resulted in the formation of new hydrophobic contacts with flap residues, residues 79 and 80, and Asp25, respectively. Mutation to smaller side chains eliminated hydrophobic interactions in the PR(I50V) and PR(I54V) structures. The PR(I84V)-APV complex had lost hydrophobic contacts with APV, the PR(V32I)-APV complex showed increased hydrophobic contacts within the hydrophobic cluster and the PR(I50V) complex had weaker polar and hydrophobic interactions with APV. The observed structural changes in PR(I84V)-APV, PR(V32I)-APV and PR(I50V)-APV were related to their reduced inhibition by APV of six-, 10- and 30-fold, respectively, relative to wild-type PR. The APV complexes were compared with the corresponding saquinavir complexes. The PR dimers had distinct rearrangements of the flaps and 80's loops that adapt to the different P1' groups of the inhibitors, while maintaining contacts within the hydrophobic cluster. These small changes in the loops and weak internal interactions produce the different patterns of resistant mutations for the two drugs.

Prediction of potency of protease inhibitors using free energy simulations with polarizable quantum mechanics-based ligand charges and a hybrid water model

Reliable and robust prediction of the binding affinity for drug molecules continues to be a daunting challenge. We simulated the binding interactions and free energy of binding of nine protease inhibitors (PIs) with wild-type and various mutant proteases by performing GBSA simulations in which each PI's partial charge was determined by quantum mechanics (QM) and the partial charge accounts for the polarization induced by the protease environment. We employed a hybrid solvation model that retains selected explicit water molecules in the protein with surface-generalized Born (SGB) implicit solvent. We examined the correlation of the free energy with the antiviral potency of PIs with regard to amino acid substitutions in protease. The GBSA free energy thus simulated showed strong correlations (r > 0.75) with antiviral IC(50) values of PIs when amino acid substitutions were present in the protease active site. We also simulated the binding free energy of PIs with P2-bis-tetrahydrofuranylurethane (bis-THF) or related cores, utilizing a bis-THF-containing protease crystal structure as a template. The free energy showed a strong correlation (r = 0.93) with experimentally determined anti-HIV-1 potency. The present data suggest that the presence of selected explicit water in protein and protein polarization-induced quantum charges for the inhibitor, compared to lack of explicit water and a static force-field-based charge model, can serve as an improved lead optimization tool and warrants further exploration.

Detecting and understanding combinatorial mutation patterns responsible for HIV drug resistance

We propose a systematic approach for a better understanding of how HIV viruses employ various combinations of mutations to resist drug treatments, which is critical to developing new drugs and optimizing the use of existing drugs. By probabilistically modeling mutations in the HIV-1 protease or reverse transcriptase (RT) isolated from drug-treated patients, we present a statistical procedure that first detects mutation combinations associated with drug resistance and then infers detailed interaction structures of these mutations. The molecular basis of our statistical predictions is further studied by using molecular dynamics simulations and free energy calculations. We have demonstrated the usefulness of this systematic procedure on three HIV drugs, (Indinavir, Zidovudine, and Nevirapine), discovered unique interaction features between viral mutations induced by these drugs, and revealed the structural basis of such interactions.

HIV-1 Protease: Structural Perspectives on Drug Resistance

Antiviral inhibitors of HIV-1 protease are a notable success of structure-based drug design and have dramatically improved AIDS therapy. Analysis of the structures and activities of drug resistant protease variants has revealed novel molecular mechanisms of drug resistance and guided the design of tight-binding inhibitors for resistant variants. The plethora of structures reveals distinct molecular mechanisms associated with resistance: mutations that alter the protease interactions with inhibitors or substrates; mutations that alter dimer stability; and distal mutations that transmit changes to the active site. These insights will inform the continuing design of novel antiviral inhibitors targeting resistant strains of HIV.

Visualizing transient events in amino-terminal autoprocessing of HIV-1 protease

HIV-1 protease processes the Gag and Gag-Pol polyproteins into mature structural and functional proteins, including itself, and is therefore indispensable for viral maturation. The mature protease is active only as a dimer with each subunit contributing catalytic residues. The full-length transframe region protease precursor appears to be monomeric yet undergoes maturation via intramolecular cleavage of a putative precursor dimer, concomitant with the appearance of mature-like catalytic activity. How such intramolecular cleavage can occur when the amino and carboxy termini of the mature protease are part of an intersubunit beta-sheet located distal from the active site is unclear. Here we visualize the early events in N-terminal autoprocessing using an inactive mini-precursor with a four-residue N-terminal extension that mimics the transframe region protease precursor. Using paramagnetic relaxation enhancement, a technique that is exquisitely sensitive to the presence of minor species, we show that the mini-precursor forms highly transient, lowly populated (3-5%) dimeric encounter complexes that involve the mature dimer interface but occupy a wide range of subunit orientations relative to the mature dimer. Furthermore, the occupancy of the mature dimer configuration constitutes a very small fraction of the self-associated species (accounting for the very low enzymatic activity of the protease precursor), and the N-terminal extension makes transient intra- and intersubunit contacts with the substrate binding site and is therefore available for autocleavage when the correct dimer orientation is sampled within the encounter complex ensemble.

Prediction of HIV-1 protease inhibitor resistance by Molecular Modeling Protocols (MMPs) using GenMol software

This paper investigates the contribution of Molecular Modeling to (i) predict and (ii) understand more fundamentally HIV drug resistance. Based on a new automated GenMol module, these goals are approached by Molecular Modeling Protocols (MMPs), respectively, (i) the Molecular Modeling Phenotype Protocol (MMPP) and (ii) the Molecular Modeling Phenotype-Genotype Protocol (MMGPP). Section 2 recalls clinical practice with a reference case study and Section 3 presents atomistic simulation tools. Section 4 is the heart of the paper. In Section 4.1, MMPP drug resistance prediction is based on correlations between fold resistances versus binding energies on 2959 HIV-1 complexes with 6 protease inhibitors. Based on a drug sensitivity twofold criterion, modeling prediction is able to replace long and costly phenotype tests. In Section 4.2, MMGPP enlightens drug resistance by investigating steric and energetic residues/inhibitor interaction. Section 5 gives a synthesis on modeling contribution to drug resistance prediction. In conclusion, the most promising trend consists of MMP automats that are able to suggest a real time diagnosis taking into account the history of each patient, to enrich databases and to develop therapy strategy and new drugs.

Effect of the active site D25N mutation on the structure, stability, and ligand binding of the mature HIV-1 protease

All aspartic proteases, including retroviral proteases, share the triplet DTG critical for the active site geometry and catalytic function. These residues interact closely in the active, dimeric structure of HIV-1 protease (PR). We have systematically assessed the effect of the D25N mutation on the structure and stability of the mature PR monomer and dimer. The D25N mutation (PR(D25N)) increases the equilibrium dimer dissociation constant by a factor >100-fold (1.3 +/- 0.09 microm) relative to PR. In the absence of inhibitor, NMR studies reveal clear structural differences between PR and PR(D25N) in the relatively mobile P1 loop (residues 79-83) and flap regions, and differential scanning calorimetric analyses show that the mutation lowers the stabilities of both the monomer and dimer folds by 5 and 7.3 degrees C, respectively. Only minimal differences are observed in high resolution crystal structures of PR(D25N) complexed to darunavir (DRV), a potent clinical inhibitor, or a non-hydrolyzable substrate analogue, Ac-Thr-Ile-Nle-r-Nle-Gln-Arg-NH(2) (RPB), as compared with PR.DRV and PR.RPB complexes. Although complexation with RPB stabilizes both dimers, the effect on their T(m) is smaller for PR(D25N) (6.2 degrees C) than for PR (8.7 degrees C). The T(m) of PR(D25N).DRV increases by only 3 degrees C relative to free PR(D25N), as compared with a 22 degrees C increase for PR.DRV, and the mutation increases the ligand dissociation constant of PR(D25N).DRV by a factor of approximately 10(6) relative to PR.DRV. These results suggest that interactions mediated by the catalytic Asp residues make a major contribution to the tight binding of DRV to PR.

Mutational patterns and correlated amino acid substitutions in the HIV-1 protease after virological failure to nelfinavir- and lopinavir/ritonavir-based treatments

Human immunodeficiency virus type 1 (HIV-1) antiviral drug resistance is a major consequence of therapy failure and compromises future therapeutic options. Nelfinavir and lopinavir/ritonavir-based therapies have been widely used in the treatment of HIV-infected patients, in combination with reverse transcriptase inhibitors. The aim of this observational study was the identification and characterization of mutations or combinations of mutations associated with resistance to nelfinavir and lopinavir/ritonavir in treated patients. Nucleotide sequences of 1,515 subtype B HIV-1 isolates from 1,313 persons with different treatment histories (including naïve and treated patients) were collected in 31 Spanish hospitals over the years 2002-2005. Chi-square contingency tests were performed to detect mutations associated with failure to protease inhibitor-based therapies, and correlated mutations were identified using statistical methods. Virological failure to nelfinavir was associated with two different mutational pathways. D30N and N88D appeared mostly in patients without previous exposure to protease inhibitors, while K20T was identified as a secondary resistance mutation in those patients. On the other hand, L90M together with L10I, I54V, A71V, G73S, and V82A were selected in protease inhibitor-experienced patients. A series of correlated mutations including L10I, M46I, I54V, A71V, G73S, and L90M appeared as a common cluster of amino acid substitutions, associated with failure to lopinavir/ritonavir-based treatments. Despite the relatively high genetic barrier of some protease inhibitors, a relatively small cluster of mutations, previously selected under drug pressure, can seriously compromise the efficiency of nelfinavir- and lopinavir/ritonavir-based therapies.

Molecular analysis of the HIV-1 resistance development: enzymatic activities, crystal structures, and thermodynamics of nelfinavir-resistant HIV protease mutants

Human immunodeficiency virus (HIV) encodes an aspartic protease (PR) that cleaves viral polyproteins into mature proteins, thus leading to the formation of infectious particles. Protease inhibitors (PIs) are successful virostatics. However, their efficiency is compromised by antiviral resistance. In the PR sequence of viral variants resistant to the PI nelfinavir, the mutations D30N and L90M appear frequently. However, these two mutations are seldom found together in vivo, suggesting that there are two alternative evolutionary pathways leading to nelfinavir resistance. Here we analyze the proteolytic activities, X-ray structures, and thermodynamics of inhibitor binding to HIV-1 PRs harboring the D30N and L90M mutations alone and in combination with other compensatory mutations. Vitality values obtained for recombinant mutant proteases and selected PR inhibitors confirm the crucial role of mutations in positions 30 and 90 for nelfinavir resistance. The combination of the D30N and L90M mutations significantly increases the enzyme vitality in the presence of nelfinavir, without a dramatic decrease in the catalytic efficiency of the recombinant enzyme. Crystal structures, molecular dynamics simulations, and calorimetric data for four mutants (D30N, D30N/A71V, D30N/N88D, and D30N/L90M) were used to augment our kinetic data. Calorimetric analysis revealed that the entropic contribution to the mutant PR/nelfinavir interaction is less favorable than the entropic contribution to the binding of nelfinavir by wild-type PR. This finding is supported by the structural data and simulations; nelfinavir binds most strongly to the wild-type protease, which has the lowest number of protein-ligand hydrogen bonds and whose structure exhibits the greatest degree of fluctuation upon inhibitor binding.

Computational simulations of HIV-1 proteases--multi-drug resistance due to nonactive site mutation L90M

Human immunodeficiency virus type 1 protease (HIV-1 PR) is one of the proteins that currently available anti-HIV-1 drugs target. Inhibitors of HIV-1 PR have become available, and they have lowered the rate of mortality from acquired immune deficiency syndrome (AIDS) in advanced countries. However, the rate of emergence of drug-resistant HIV-1 variants is quite high because of their short retroviral life cycle and their high mutation rate. Serious drug-resistant mutations against HIV-1 PR inhibitors (PIs) frequently appear at the active site of PR. Exceptionally, some other mutations such as L90M cause drug resistance, although these appear at nonactive sites. The mechanism of resistance due to nonactive site mutations is difficult to explain. In this study, we carried out computational simulations of L90M PR in complex with each of three kinds of inhibitors and one typical substrate, and we clarified the mechanism of resistance. The L90M mutation causes changes in interaction between the side chain atoms of the 90th residue and the main chain atoms of the 25th residue, and a slight dislocation of the 25th residue causes rotation of the side chain at the 84th residue. The rotation of the 84th residue leads to displacement of the inhibitor from the appropriate binding location, resulting in a collision with the flap or loop region. The difference in levels of resistance to the three inhibitors has been explained from energetic and structural viewpoints, which provides the suggestion for promising drugs keeping its efficacy even for the L90M mutant.

Potent new antiviral compound shows similar inhibition and structural interactions with drug resistant mutants and wild type HIV-1 protease

The potent new antiviral inhibitor GRL-98065 (1) of HIV-1 protease (PR) has been studied with PR variants containing the single mutations D30N, I50V, V82A, and I84V that provide resistance to the major clinical inhibitors. Compound 1 had inhibition constants of 17-fold, 8-fold, 3-fold, and 3-fold, respectively, for PR(D30N), PR(I50V), PR(V82A), and PR(I84V) relative to wild type PR. The chemically related darunavir had similar relative inhibition, except for PR(D30N), where inhibitor 1 was approximately 3-fold less potent. The high resolution (1.11-1.60 Angstrom) crystal structures of PR mutant complexes with inhibitor 1 showed small changes relative to the wild type enzyme. PR(D30N) and PR(V82A) showed compensating interactions with inhibitor 1 relative to those of PR, while reduced hydrophobic contacts were observed with PR(I50V) and PR(I84V). Importantly, inhibitor 1 complexes showed fewer changes relative to wild type enzyme than reported for darunavir complexes. Therefore, inhibitor 1 is a valuable addition to the antiviral inhibitors with high potency against resistant strains of HIV.