Insights into a mutation-assisted lateral drug escape mechanism from the HIV-1 protease active site

We provide insight into the first stages of a kinetic mechanism of lateral drug expulsion from the active site of HIV-1 protease, by conducting all atom molecular dynamics simulations with explicit solvent over a time scale of 24 ns for saquinavir bound to the wildtype, G48V, L90M and G48V/L90M mutant proteases. We find a consistent escape mechanism associated with the G48V mutation. First, increased hydrophilic and hydrophobic flap coupling and water mediated disruption of catalytic dyad hydrogen bonding induce drug motion away from the dyad and promote protease flap transition to the semi-open form. Conversely, flap-inhibitor motion is decoupled in the wildtype. Second, the decrease of total interactions causes unidirectional lateral inhibitor translation by up to 4 A toward the P3 subsite exit of the active site, increased P3 subsite exposure to solvent and a complete loss of hydrophobic interactions with the opposite end of the active site. The P1 subsite moves beyond the hydrophobic active site side pocket, the only remaining steric barrier to complete expulsion being the "breathable" residue, P81. Significant inhibitor deviation is reported over 24 ns, and subsequent complete expulsion, implemented using steered molecular dynamics simulations, is shown to occur most easily for the G48V-containing mutants. Our simulations thus provide compelling support for lateral drug escape from a protease in a semi-open flap conformation. It is likely that some mutations take advantage of this escape mechanism to increase the rate of inhibitor dissociation from the protease. Finally, unidirectional translation may be countered by designing inhibitors with terminal subsites that provide sufficient anchoring to the flaps, thus increasing the steric barrier for translation in either direction.     

Structural and kinetic analysis of drug resistant mutants of HIV-1 protease.

Mutants of HIV-1 protease that are commonly selected on exposure to different drugs, V82S, G48V, N88D and L90M, showed reduced catalytic activity compared to the wild-type protease on cleavage site peptides, CA-p2, p6pol-PR and PR-RT, critical for viral maturation. Mutant V82S is the least active (2-20% of wild-type protease), mutants N88D, R8Q, and L90M exhibit activities ranging from 20 to 40% and G48V from 50 to 80% of the wild-type activity. In contrast, D30N is variable in its activity on different substrates (10-110% of wild-type), with the PR-RT site being the most affected. Mutants K45I and M46L, usually selected in combination with other mutations, showed activities that are similar to (60-110%) or greater than (110-530%) wild-type, respectively. No direct relationship was observed between catalytic activity, inhibition, and structural stability. The mutants D30N and V82S were similar to wild-type protease in their stability toward urea denaturation, while R8Q, G48V, and L90M showed 1.5 to 2.7-fold decreased stability, and N88D and K45I showed 1.6 to 1.7-fold increased stability. The crystal structures of R8Q, K45I and L90M mutants complexed with a CA-p2 analog inhibitor were determined at 2.0, 1.55 and 1.88 A resolution, respectively, and compared to the wild-type structure. The intersubunit hydrophobic contacts observed in the crystal structures are in good agreement with the relative structural stability of the mutant proteases. All these results suggest that viral resistance does not arise by a single mechanism.     

Replicative fitness of protease inhibitor-resistant mutants of human immunodeficiency virus type 1

The relative replicative fitness of human immunodeficiency virus type 1 (HIV-1) mutants selected by different protease inhibitors (PIs) in vivo was determined. Each mutant was compared to wild type (WT), NL4-3, in the absence of drugs by several methods, including clonal genotyping of cultures infected with two competing viral variants, kinetics of viral antigen production, and viral infectivity/virion particle ratios. A nelfinavir-selected protease D30N substitution substantially decreased replicative capacity relative to WT, while a saquinavir-selected L90M substitution moderately decreased fitness. The D30N mutant virus was also outcompeted by the L90M mutant in the absence of drugs. A major natural polymorphism of the HIV-1 protease, L63P, compensated well for the impairment of fitness caused by L90M but only slightly improved the fitness of D30N. Multiply substituted indinavir-selected mutants M46I/L63P/V82T/I84V and L10R/M46I/L63P/V82T/I84V were just as fit as WT. These results indicate that the mutations which are usually initially selected by nelfinavir and saquinavir, D30N and L90M, respectively, impair fitness. However, additional mutations may improve the replicative capacity of these and other drug-resistant mutants. Hypotheses based on the greater fitness impairment of the nelfinavir-selected D30N mutant are suggested to explain observations that prolonged responses to delayed salvage regimens, including alternate PIs, may be relatively common after nelfinavir failure.     

Unique thermodynamic response of tipranavir to human immunodeficiency virus type 1 protease drug resistance mutations

Drug resistance is a major problem affecting the clinical efficacy of antiretroviral agents, including protease inhibitors, in the treatment of infection with human immunodeficiency virus type 1 (HIV-1)/AIDS. Consequently, the elucidation of the mechanisms by which HIV-1 protease inhibitors maintain antiviral activity in the presence of mutations is critical to the development of superior inhibitors. Tipranavir, a nonpeptidic HIV-1 protease inhibitor, has been recently approved for the treatment of HIV infection. Tipranavir inhibits wild-type protease with high potency (K(i) = 19 pM) and demonstrates durable efficacy in the treatment of patients infected with HIV-1 strains containing multiple common mutations associated with resistance. The high potency of tipranavir results from a very large favorable entropy change (-TDeltaS = -14.6 kcal/mol) combined with a favorable, albeit small, enthalpy change (DeltaH = -0.7 kcal/mol, 25 degrees C). Characterization of tipranavir binding to wild-type protease, active site mutants I50V and V82F/I84V, the multidrug-resistant mutant L10I/L33I/M46I/I54V/L63I/V82A/I84V/L90M, and the tipranavir in vitro-selected mutant I13V/V32L/L33F/K45I/V82L/I84V was performed by isothermal titration calorimetry and crystallography. Thermodynamically, the good response of tipranavir arises from a unique behavior: it compensates for entropic losses by actual enthalpic gains or by sustaining minimal enthalpic losses when facing the mutants. The net result is a small loss in binding affinity. Structurally, tipranavir establishes a very strong hydrogen bond network with invariant regions of the protease, which is maintained with the mutants, including catalytic Asp25 and the backbone of Asp29, Asp30, Gly48 and Ile50. Moreover, tipranavir forms hydrogen bonds directly to Ile50, while all other inhibitors do so by being mediated by a water molecule.     

Molecular mechanics analysis of drug-resistant mutants of HIV protease.

Drug-resistant mutants of HIV-1 protease limit the long-term effectiveness of current anti-viral therapy. In order to study drug resistance, the wild-type HIV-1 protease and the mutants R8Q, V32I, M46I, V82A, V82I, V82F, I84V, V32I/I84V and M46I/I84V were modeled with the inhibitors saquinavir and indinavir using the program AMMP. A new screen term was introduced to reproduce more correctly the electron distribution of atoms. The atomic partial charge was represented as a delocalized charge distribution instead of a point charge. The calculated protease-saquinavir interaction energies showed the highly significant correlation of 0.79 with free energy differences derived from the measured inhibition constants for all 10 models. Three different protonation states of indinavir were evaluated. The best indinavir model included a sulfate and gave a correlation coefficient of 0.68 between the calculated interaction energies and free energies from inhibition constants for nine models. The exception was R8Q with indinavir, probably due to differences in the solvation energy. No significant correlation was found using the standard molecular mechanics terms. The incorporation of the new screen correction resulted in better prediction of the effects of inhibitors on resistant protease variants and has potential for selecting more effective inhibitors for resistant virus.     

Multidrug resistance to HIV-1 protease inhibition requires cooperative coupling between distal mutations

The appearance of viral strains that are resistant to protease inhibitors is one of the most serious problems in the chemotherapy of HIV-1/AIDS. The most pervasive drug-resistant mutants are those that affect all inhibitors in clinical use. In this paper, we have characterized a multiple-drug-resistant mutant of the HIV-1 protease that affects indinavir, nelfinavir, saquinavir, ritonavir, amprenavir, and lopinavir. This mutant (MDR-HM) contains six amino acid mutations (L10I/M46I/I54V/V82A/I84V/L90M) located within and outside the active site of the enzyme. Microcalorimetric and enzyme kinetic measurements indicate that this mutant lowers the affinity of all inhibitors by 2-3 orders of magnitude. By comparison, the multiiple-drug-resistant mutant only increased the K(m) of the substrate by a factor of 2, indicating that the substrate is able to adapt to the changes caused by the mutations and maintain its binding affinity. To understand the origin of resistance, three submutants containing mutations in specific regions were also studied, i.e., the active site (V82A/I84V), flap region (M46I/I54V), and dimerization region (L10I/L90M). None of these sets of mutations by themselves lowered the affinity of inhibitors by more than 1 order of magnitude, and additionally, the sum of the effects of each set of mutations did not add up to the overall effect, indicating the presence of cooperative effects. A mutant containing only the four active site mutations (V82A/I84V/M46I/I54V) only showed a small cooperative effect, suggesting that the mutations at the dimer interface (L10I/L90M) play a major role in eliciting a cooperative response. These studies demonstrate that cooperative interactions contribute an average of 1.2 +/- 0.7 kcal/mol to the overall resistance, most of the cooperative effect (0.8 +/- 0.7 kcal/mol) being mediated by the mutations at the dimerization interface. Not all inhibitors in clinical use are affected the same by long-range cooperative interactions between mutations. These interactions can amplify the effects of individual mutations by factors ranging between 2 and 40 depending on the inhibitor. Dissection of the energetics of drug resistance into enthalpic and entropic components provides a quantitative account of the inhibitor response and a set of thermodynamic guidelines for the design of inhibitors with a lower susceptibility to this type of mutations.     

Identification of Genotypic Changes in Human Immunodeficiency Virus Protease That Correlate with Reduced Susceptibility to the Protease Inhibitor Lopinavir among Viral Isolates from Protease Inhibitor-Experienced Patients

The association of genotypic changes in human immunodeficiency virus (HIV) protease with reduced in vitro susceptibility to the new protease inhibitor lopinavir (previously ABT-378) was explored using a panel of viral isolates from subjects failing therapy with other protease inhibitors. Two statistical tests showed that specific mutations at 11 amino acid positions in protease (L10F/I/R/V, K20M/R, L24I, M46I/L, F53L, I54L/T/V, L63P, A71I/L/T/V, V82A/F/T, I84V, and L90M) were associated with reduced susceptibility. Mutations at positions 82, 54, 10, 63, 71, and 84 were most closely associated with relatively modest (4- and 10-fold) changes in phenotype, while the K20M/R and F53L mutations, in conjunction with multiple other mutations, were associated with >20- and >40-fold-reduced susceptibility, respectively. The median 50% inhibitory concentrations (IC(50)) of lopinavir against isolates with 0 to 3, 4 or 5, 6 or 7, and 8 to 10 of the above 11 mutations were 0.8-, 2.7-, 13.5-, and 44.0-fold higher, respectively, than the IC(50) against wild-type HIV. On average, the IC(50) of lopinavir increased by 1.74-fold per mutation in isolates containing three or more mutations. Each of the 16 viruses that displayed a >20-fold change in susceptibility contained mutations at residues 10, 54, 63, and 82 and/or 84, along with a median of three mutations at residues 20, 24, 46, 53, 71, and 90. The number of protease mutations from the 11 identified in these analyses (the lopinavir mutation score) may be useful for the interpretation of HIV genotypic resistance testing with respect to lopinavir-ritonavir (Kaletra) regimens and may provide insight into the genetic barrier to resistance to lopinavir-ritonavir in both antiretroviral therapy-naive and protease inhibitor-experienced patients.     

Combining mutations in HIV-1 protease to understand mechanisms of resistance

HIV-1 develops resistance to protease inhibitors predominantly by selecting mutations in the protease gene. Studies of resistant mutants of HIV-1 protease with single amino acid substitutions have shown a range of independent effects on specificity, inhibition, and stability. Four double mutants, K45I/L90M, K45I/V82S, D30N/V82S, and N88D/L90M were selected for analysis on the basis of observations of increased or decreased stability or enzymatic activity for the respective single mutants. The double mutants were assayed for catalysis, inhibition, and stability. Crystal structures were analyzed for the double mutants at resolutions of 2.2-1.2 A to determine the associated molecular changes. Sequence-dependent changes in protease-inhibitor interactions were observed in the crystal structures. Mutations D30N, K45I, and V82S showed altered interactions with inhibitor residues at P2/P2', P3/P3'/P4/P4', and P1/P1', respectively. One of the conformations of Met90 in K45I/L90M has an unfavorably close contact with the carbonyl oxygen of Asp25, as observed previously in the L90M single mutant. The observed catalytic efficiency and inhibition for the double mutants depended on the specific substrate or inhibitor. In particular, large variation in cleavage of p6(pol)-PR substrate was observed, which is likely to result in defects in the maturation of the protease from the Gag-Pol precursor and hence viral replication. Three of the double mutants showed values for stability that were intermediate between the values observed for the respective single mutants. D30N/V82S mutant showed lower stability than either of the two individual mutations, which is possibly due to concerted changes in the central P2-P2' and S2-S2' sites. The complex effects of combining mutations are discussed.     

Comparing the accumulation of active- and nonactive-site mutations in the HIV-1 protease

Protease inhibitor resistance still poses one of the greatest challenges in treating HIV. To better design inhibitors able to target resistant proteases, a deeper understanding is needed of the effects of accumulating mutations and the contributions of active- and nonactive-site mutations to the resistance. We have engineered a series of variants containing the nonactive-site mutations M46I and I54V and the active-site mutation I84V. These mutations were added to a protease clone (V6) isolated from a pediatric patient on ritonavir therapy. This variant possessed the ritonavir-resistance-associated mutations in the active-site (V32I and V82A) and nonactive-site mutations (K20R, L33F, M36I, L63P, A71V, and L90M). The I84V mutation had the greatest effect on decreasing catalytic efficiency, 10-fold when compared to the pretherapy clone LAI. The decrease in catalytic efficiency was partially recovered by the addition of mutations M46I and I54V. The M46I and I54V were just as effective at decreasing inhibitor binding as the I84V mutation when compared to V6 and LAI. The V6(54/84) variant showed over 1000-fold decrease in inhibitor-binding strength to ritonavir, indinavir, and nelfinavir when compared to LAI and V6. Crystal-structure analysis of the V6(54/84) variant bound to ritonavir and indinavir shows structural changes in the 80's loops and active site, which lead to an enlarged binding cavity when compared to pretherapy structures in the Protein Data Bank. Structural changes are also seen in the 10's and 30's loops, which suggest possible changes in the dynamics of flap opening and closing.     

Analysis of resistance to human immunodeficiency virus type 1 protease inhibitors by using matched bacterial expression and proviral infection vectors.

There are already reports, from clinical trials with human immunodeficiency virus type 1 protease inhibitors, of the emergence of drug-resistant mutants which have one or more point mutations in their protease genes. To examine roles of individual and multiple amino acid substitutions in terms of altered enzyme and virus drug sensitivities, we have produced matched vectors for bacterial expression and virus production. Both vectors accept the same restriction enzyme fragment, produced by PCR or PCR-mutagenesis of the protease gene, allowing parallel expression of mutant enzymes in Escherichia coli and in recombinant viruses. The utility of this vector system was demonstrated by using protease variants glycine to valine at amino acid 48 (G48V) and leucine to methionine at amino acid 90 (L90M) identified after passage of HIV-1 in the Roche phase II clinical trial protease inhibitor Ro 31-8959 (H. Jacobsen, K. Yasargil, D. L. Winslow, J. C. Craig, A. Krohn, I. B. Duncan, and J. Mous, Virology 206:527, 1995). G48V, L90M, and G48V/L90M exhibited successively less processing in vitro than the wild-type enzyme, and the purified enzymes were 220-, 20-, and 720-fold, respectively, less sensitive to Ro 31-8959. The reduced enzyme sensitivity correlated directly with the sensitivities of the matched recombinant viruses, in that individual mutations L90M and G48V conferred 2-fold and 4- to 6-fold increases in 50% inhibitory concentration, respectively, whereas G48V/L90M was 8 to 10 times less sensitive to Ro 31-8959. A proviral vector with the entire protease gene deleted was constructed for use as an in vivo recombination target for an overlapping protease PCR fragment, generating wild-type infectious virus. Finally, direct ligation of restriction fragments, generated from random PCR mutagenesis, into the proviral vector should provide a library of protease mutations that allow extremely rapid selection of highly resistant viral variants.     

Crystal structure of an in vivo HIV-1 protease mutant in complex with saquinavir: insights into the mechanisms of drug resistance

Saquinavir is a widely used HIV-1 protease inhibitor drug for AIDS therapy. Its effectiveness, however, has been hindered by the emergence of resistant mutations, a common problem for inhibitor drugs that target HIV-1 viral enzymes. Three HIV-1 protease mutant species, G48V, L90M, and G48V/L90M double mutant, are associated in vivo with saquinavir resistance by the enzyme (Jacobsen et al., 1996). Kinetic studies on these mutants demonstrate a 13.5-, 3-, and 419-fold increase in Ki values, respectively, compared to the wild-type enzyme (Ermolieff J, Lin X, Tang J, 1997, Biochemistry 36:12364-12370). To gain an understanding of how these mutations modulate inhibitor binding, we have solved the HIV-1 protease crystal structure of the G48V/L90M double mutant in complex with saquinavir at 2.6 A resolution. This mutant complex is compared with that of the wild-type enzyme bound to the same inhibitor (Krohn A, Redshaw S, Richie JC, Graves BJ, Hatada MH, 1991, J Med Chem 34:3340-3342). Our analysis shows that to accommodate a valine side chain at position 48, the inhibitor moves away from the protease, resulting in the formation of larger gaps between the inhibitor P3 subsite and the flap region of the enzyme. Other subsites also demonstrate reduced inhibitor interaction due to an overall change of inhibitor conformation. The new methionine side chain at position 90 has van der Waals interactions with main-chain atoms of the active site residues resulting in a decrease in the volume and the structural flexibility of S1/S1' substrate binding pockets. Indirect interactions between the mutant methionine side chain and the substrate scissile bond or the isostere part of the inhibitor may differ from those of the wild-type enzyme and therefore may facilitate catalysis by the resistant mutant.     

Drug resistance mutations can effect dimer stability of HIV-1 protease at neutral pH.

The monomer-dimer equilibrium for the human immunodeficiency virus type 1 (HIV-1) protease has been investigated under physiological conditions. Dimer dissociation at pH 7.0 was correlated with a loss in beta-sheet structure and a lower degree of ANS binding. An autolysis-resistant mutant, Q7K/L33I/L63I, was used to facilitate sedimentation equilibrium studies at neutral pH where the wild-type enzyme is typically unstable in the absence of bound inhibitor. The dimer dissociation constant (KD) of the triple mutant was 5.8 microM at pH 7.0 and was below the limit of measurement (approximately 100 nM) at pH 4.5. Similar studies using the catalytically inactive D25N mutant yielded a KD value of 1.0 microM at pH 7.0. These values differ significantly from a previously reported value of 23 nM obtained indirectly from inhibitor binding measurements (Darke et al., 1994). We show that the discrepancy may result from the thermodynamic linkage between the monomer-dimer and inhibitor binding equilibria. Under conditions where a significant degree of monomer is present, both substrates and competitive inhibitors will shift the equilibrium toward the dimer, resulting in apparent increases in dimer stability and decreases in ligand binding affinity. Sedimentation equilibrium studies were also carried out on several drug-resistant HIV-1 protease mutants: V82F, V82F/I84V, V82T/I84V, and L90M. All four mutants exhibited reduced dimer stability relative to the autolysis-resistant mutant at pH 7.0. Our results indicate that reductions in drug affinity may be due to the combined effects of mutations on both dimer stability and inhibitor binding.     

HIV-1 protease molecular dynamics of a wild-type and of the V82F/I84V mutant: possible contributions to drug resistance and a potential new target site for drugs

The protease from type 1 human immunodeficiency virus (HIV-1) is a critical drug target against which many therapeutically useful inhibitors have been developed; however, the set of viral strains in the population has been shifting to become more drug-resistant. Because indirect effects are contributing to drug resistance, an examination of the dynamic structures of a wild-type and a mutant could be insightful. Consequently, this study examined structural properties sampled during 22 nsec, all atom molecular dynamics (MD) simulations (in explicit water) of both a wild-type and the drug-resistant V82F/I84V mutant of HIV-1 protease. The V82F/I84V mutation significantly decreases the binding affinity of all HIV-1 protease inhibitors currently used clinically. Simulations have shown that the curling of the tips of the active site flaps immediately results in flap opening. In the 22-nsec MD simulations presented here, more frequent and more rapid curling of the mutant's active site flap tips was observed. The mutant protease's flaps also opened farther than the wild-type's flaps did and displayed more flexibility. This suggests that the effect of the mutations on the equilibrium between the semiopen and closed conformations could be one aspect of the mechanism of drug resistance for this mutant. In addition, correlated fluctuations in the active site and periphery were noted that point to a possible binding site for allosteric inhibitors.     

Molecular basis of telaprevir resistance due to V36 and T54 mutations in the NS3-4A protease of the hepatitis C virus

Background

The inhibitor telaprevir (VX-950) of the hepatitis C virus (HCV) protease NS3-4A has been tested in a recent phase 1b clinical trial in patients infected with HCV genotype 1. This trial revealed residue mutations that confer varying degrees of drug resistance. In particular, two protease positions with the mutations V36A/G/L/M and T54A/S were associated with low to medium levels of drug resistance during viral breakthrough, together with only an intermediate reduction of viral replication fitness. These mutations are located in the protein interior and far away from the ligand binding pocket.

Results

Based on the available experimental structures of NS3-4A, we analyze the binding mode of different ligands. We also investigate the binding mode of VX-950 by protein-ligand docking. A network of non-covalent interactions between amino acids of the protease structure and the interacting ligands is analyzed to discover possible mechanisms of drug resistance. We describe the potential impact of V36 and T54 mutants on the side chain and backbone conformations and on the non-covalent residue interactions. We propose possible explanations for their effects on the antiviral efficacy of drugs and viral fitness. Molecular dynamics simulations of T54A/S mutants and rotamer analysis of V36A/G/L/M side chains support our interpretations. Experimental data using an HCV V36G replicon assay corroborate our findings.

Conclusion

T54 mutants are expected to interfere with the catalytic triad and with the ligand binding site of the protease. Thus, the T54 mutants are assumed to affect the viral replication efficacy to a larger degree than V36 mutants. Mutations at V36 and/or T54 result in impaired interaction of the protease residues with the VX-950 cyclopropyl group, which explains the development of viral breakthrough variants.     

Effect of flap mutations on structure of HIV-1 protease and inhibition by saquinavir and darunavir

HIV-1 (human immunodeficiency virus type 1) protease (PR) and its mutants are important antiviral drug targets. The PR flap region is critical for binding substrates or inhibitors and catalytic activity. Hence, mutations of flap residues frequently contribute to reduced susceptibility to PR inhibitors in drug-resistant HIV. Structural and kinetic analyses were used to investigate the role of flap residues Gly48, Ile50, and Ile54 in the development of drug resistance. The crystal structures of flap mutants PR_I50V (PR with I50V mutation), PR_I54V (PR with I54V mutation), and PR_I54M (PR with I54M mutation) complexed with saquinavir (SQV) as well as PR_G48V (PR with G48V mutation), PR_I54V, and PR_I54M complexed with darunavir (DRV) were determined at resolutions of 1.05-1.40 Å. The PR mutants showed changes in flap conformation, interactions with adjacent residues, inhibitor binding, and the conformation of the 80s loop relative to the wild-type PR. The PR contacts with DRV were closer in PR_G48V-DRV than in the wild-type PR-DRV, whereas they were longer in PR_I54M-DRV. The relative inhibition of PR_I54V and that of PR_I54M were similar for SQV and DRV. PR_G48V was about twofold less susceptible to SQV than to DRV, whereas the opposite was observed for PR_I50V. The observed inhibition was in agreement with the association of G48V and I50V with clinical resistance to SQV and DRV, respectively. This analysis of structural and kinetic effects of the mutants will assist in the development of more effective inhibitors for drug-resistant HIV.     

Overcoming drug resistance in HIV-1 chemotherapy: the binding thermodynamics of Amprenavir and TMC-126 to wild-type and drug-resistant mutants of the HIV-1 protease

Amprenavir is one of six protease inhibitors presently approved for clinical use in the therapeutic treatment of AIDS. Biochemical and clinical studies have shown that, unlike other inhibitors, Amprenavir is severely affected by the protease mutation I50V, located in the flap region of the enzyme. TMC-126 is a second-generation inhibitor, chemically related to Amprenavir, with a reported extremely low susceptibility to existing resistant mutations including I50V. In this paper, we have studied the thermodynamic and molecular origin of the response of these two inhibitors to the I50V mutation and the double active-site mutation V82F/I84V that affects all existing clinical inhibitors. Amprenavir binds to the wild-type HIV-1 protease with high affinity (5.0 x 10(9) M(-1) or 200 pM) in a process equally favored by enthalpic and entropic contributions. The mutations I50V and V82F/I84V lower the binding affinity of Amprenavir by a factor of 147 and 104, respectively. TMC-126, on the other hand, binds to the wild-type protease with extremely high binding affinity (2.6 x 10(11) M(-1) or 3.9 pM) in a process in which enthalpic contributions overpower entropic contributions by almost a factor of 4. The mutations I50V and V82F/I84V lower the binding affinity of TMC-126 by only a factor of 16 and 11, respectively, indicating that the binding affinity of TMC-126 to the drug-resistant mutants is still higher than the affinity of Amprenavir to the wild-type protease. Analysis of the data for TMC-126 and KNI-764, another second-generation inhibitor, indicates that their low susceptibility to mutations is caused by their ability to compensate for the loss of interactions with the mutated target by a more favorable entropy of binding.     

Studies on adaptability of binding residues and flap region of TMC-114 resistance HIV-1 protease mutants

Drug resistant mutations have severely restricted the success of HIV therapy. These mutations frequently involve the aspartic protease encoded by the virus. Knowledge of the molecular mechanisms underlying the conformational changes of HIV-1 protease mutants may be useful in developing more effective and longer lasting treatment regimes. The flap regions of the protease are the target of a particular type of mutations occurring far from the active site, which are able to produce significant resistance against the anti-HIV drug TMC-114. We provide insight into the molecular basis of TMC-114 resistance major flap mutations (I50V and I54M) in HIV-1 protease. It reports the shape complementarity and receptor-ligand interaction analysis supported by unrestrained all-atom molecular dynamics simulations of wild and major flap mutants of HIV-1 protease that sample large conformational changes of the flaps and active site binding residues. Both resistant flap mutants showed less atomic interaction toward TMC-114 and more structural deviation compared to wild HIV-protease. It is due to increasing flexibility at TMC-114 binding cavity and deviation of binding residues in 3-D space. Distortion in binding cavity and deviation in binding residues are the result of alteration in hydrogen bonding. Flap region also exhibited similar behaviour due to changes in number of hydrogen bonds during simulations.     

Structural insights into the mechanisms of drug resistance in HIV-1 protease NL4-3

The development of resistance to anti-retroviral drugs targeted against HIV is an increasing clinical problem in the treatment of HIV-1-infected individuals. Many patients develop drug-resistant strains of the virus after treatment with inhibitor cocktails (HAART therapy), which include multiple protease inhibitors. Therefore, it is imperative that we understand the mechanisms by which the viral proteins, in particular HIV-1 protease, develop resistance. We have determined the three-dimensional structure of HIV-1 protease NL4-3 in complex with the potent protease inhibitor TL-3 at 2.0 A resolution. We have also obtained the crystal structures of three mutant forms of NL4-3 protease containing one (V82A), three (V82A, M46I, F53L) and six (V82A, M46I, F53L, V77I, L24I, L63P) point mutations in complex with TL-3. The three protease mutants arose sequentially under ex vivo selective pressure in the presence of TL-3, and exhibit fourfold, 11-fold, and 30-fold resistance to TL-3, respectively. This series of protease crystal structures offers insights into the biochemical and structural mechanisms by which the enzyme can overcome inhibition by TL-3 while recovering some of its native catalytic activity.