The 80's loop (residues 78 to 85) is important for the differential activity of retroviral proteases

Quantum mechanical analysis of newly synthesized HIV-1 protease inhibitors: evaluation of wild-type and resistant strain binding interactions

Inhibition of HIV-1 protease is a cornerstone of antiretroviral therapy. However, the notorious ability of HIV-1 to develop resistance to protease inhibitors (PIs), particularly darunavir (DRV), poses a major challenge. Using quantum chemistry and computer simulations, this study aims to investigate the interactions between two novel PIs, GRL-004 and GRL-063, as well as a wild-type (WT) HIV strain and a DRV-resistant mutant strain. To do this, we used molecular docking, molecular dynamics simulations, and quantum mechanical calculations to check how well GRL-004 and GRL-063 bound to both WT and DRV-resistant proteases. The results show that GRL-004 and GRL-063 bind very well to ASP29 in the WT strain. ASP29 is an important amino acid in the HIV protease dimer. Remarkably, amino acids such as ILE50 in the WT strains showed substantial binding energies to both drugs. Quantum energy calculations showed a slight reduction in the energy affinity of the interaction between the MUT strain and the ligand GRL-063, compared to the WT strain. GRL-004 showed similar interaction energy for both strains, suggesting that it has greater plasticity than GRL-063 despite its lower interaction affinity. Furthermore, GLY49B demonstrated strong binding energies regardless of mutations. Other relevant residues with strong binding energies include GLY49B, PHE82A, PRO81A, ASP29A, ASP25A and ALA28B. This study improves our understanding of receptor-ligand dynamics and the adaptability of new protease inhibitors (PIs), which has profound implications for the innovation of future antiretroviral drugs.

C-5-Modified Tetrahydropyrano-Tetrahydofuran-Derived Protease Inhibitors (PIs) Exert Potent Inhibition of the Replication of HIV-1 Variants Highly Resistant to Various PIs, including Darunavir

We identified three nonpeptidic HIV-1 protease inhibitors (PIs), GRL-015, -085, and -097, containing tetrahydropyrano-tetrahydrofuran (Tp-THF) with a C-5 hydroxyl. The three compounds were potent against a wild-type laboratory HIV-1 strain (HIV-1(WT)), with 50% effective concentrations (EC50s) of 3.0 to 49 nM, and exhibited minimal cytotoxicity, with 50% cytotoxic concentrations (CC50) for GRL-015, -085, and -097 of 80, >100, and >100 μM, respectively. All the three compounds potently inhibited the replication of highly PI-resistant HIV-1 variants selected with each of the currently available PIs and recombinant clinical HIV-1 isolates obtained from patients harboring multidrug-resistant HIV-1 variants (HIVMDR). Importantly, darunavir (DRV) was >1,000 times less active against a highly DRV-resistant HIV-1 variant (HIV-1DRV(R) P51); the three compounds remained active against HIV-1DRV(R) P51 with only a 6.8- to 68-fold reduction. Moreover, the emergence of HIV-1 variants resistant to the three compounds was considerably delayed compared to the case of DRV. In particular, HIV-1 variants resistant to GRL-085 and -097 did not emerge even when two different highly DRV-resistant HIV-1 variants were used as a starting population. In the structural analyses, Tp-THF of GRL-015, -085, and -097 showed strong hydrogen bond interactions with the backbone atoms of active-site amino acid residues (Asp29 and Asp30) of HIV-1 protease. A strong hydrogen bonding formation between the hydroxyl moiety of Tp-THF and a carbonyl oxygen atom of Gly48 was newly identified. The present findings indicate that the three compounds warrant further study as possible therapeutic agents for treating individuals harboring wild-type HIV and/or HIVMDR.

IMPORTANCE

Darunavir (DRV) inhibits the replication of most existing multidrug-resistant HIV-1 strains and has a high genetic barrier. However, the emergence of highly DRV-resistant HIV-1 strains (HIVDRV(R) ) has recently been observed in vivo and in vitro. Here, we identified three novel HIV-1 protease inhibitors (PIs) containing a tetrahydropyrano-tetrahydrofuran (Tp-THF) moiety with a C-5 hydroxyl (GRL-015, -085, and -097) which potently suppress the replication of HIVDRV(R) . Moreover, the emergence of HIV-1 strains resistant to the three compounds was considerably delayed compared to the case of DRV. The C-5 hydroxyl formed a strong hydrogen bonding interaction with the carbonyl oxygen atom of Gly48 of protease as examined in the structural analyses. Interestingly, a compound with Tp-THF lacking the hydroxyl moiety substantially decreased activity against HIVDRV(R) . The three novel compounds should be further developed as potential drugs for treating individuals harboring wild-type and multi-PI-resistant HIV variants as well as HIVDRV(R) .

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.

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.

Narrow substrate specificity and sensitivity toward ligand-binding site mutations of human T-cell Leukemia virus type 1 protease

Human T-cell leukemia virus type 1 (HTLV-1) is associated with a number of human diseases; therefore, its protease is a potential target for chemotherapy. To compare the specificity of HTLV-1 protease with that of human immunodeficiency virus type 1 (HIV-1) protease, oligopeptides representing naturally occurring cleavage sites in various retroviruses were tested. The number of hydrolyzed peptides as well as the specificity constants suggested a substantially broader specificity of the HIV protease. Amino acid residues of HTLV-1 protease substrate-binding sites were replaced by equivalent ones of HIV-1 protease. Most of the single and multiple mutants had altered specificity and a dramatically reduced folding and catalytic capability, suggesting that mutations are not well tolerated in HTLV-1 protease. The catalytically most efficient mutant was that with the flap residues of HIV-1 protease. The inhibition profile of the mutants was also determined for five inhibitors used in clinical practice and inhibitor analogs of HTLV-1 cleavage sites. Except for indinavir, the HIV-1 protease inhibitors did not inhibit wild type and most of the mutant HTLV-1 proteases. The wild type HTLV-1 protease was inhibited by the reduced peptide bond-containing substrate analogs, whereas the mutants showed various degrees of weakened binding capability. Most interesting, the enzyme with HIV-1-like residues in the flap region was the most sensitive to the HIV-1 protease inhibitors and least sensitive to the HTLV-1 protease inhibitors, indicating that the flap plays an important role in defining the specificity differences of retroviral proteases.

Drug resistance in HIV-1 protease: Flexibility-assisted mechanism of compensatory mutations

The emergence of drug-resistant variants is a serious side effect associated with acquired immune deficiency syndrome therapies based on inhibition of human immunodeficiency virus type 1 protease (HIV-1 PR). In these variants, compensatory mutations, usually located far from the active site, are able to affect the enzymatic activity via molecular mechanisms that have been related to differences in the conformational flexibility, although the detailed mechanistic aspects have not been clarified so far. Here, we perform multinanosecond molecular dynamics simulations on L63P HIV-1 PR, corresponding to the wild type, and one of its most frequently occurring compensatory mutations, M46I, complexed with the substrate and an enzymatic intermediate. The quality of the calculations is established by comparison with the available nuclear magnetic resonance data. Our calculations indicate that the dynamical fluctuations of the mutated enzyme differ from those in the wild type. These differences in the dynamic properties of the adducts with the substrate and with the gem-diol intermediate might be directly related to variations in the enzymatic activity and therefore offer an explanation of the observed changes in catalytic rate between wild type and mutated enzyme. We anticipate that this "flexibility-assisted" mechanism might be effective in the vast majority of compensatory mutations, which do not change the electrostatic properties of the enzyme.

Role of conformational fluctuations in the enzymatic reaction of HIV-1 protease

The emergence of compensatory drug-resistant mutations in HIV-1 protease challenges the common view of the reaction mechanism of this enzyme. Here, we address this issue by performing classical and ab initio molecular dynamics simulations (MD) on a complex between the enzyme and a peptide substrate. The classical MD calculation reveals large-scale protein motions involving the flaps and the cantilever. These motions modulate the conformational properties of the substrate at the cleavage site. The ab initio calculations show in turn that substrate motion modulates the activation free energy barrier of the enzymatic reaction dramatically. Thus, the catalytic power of the enzyme does not arise from the presence of a pre-organized active site but from the protein mechanical fluctuations. The implications of this finding for the emergence of drug-resistance are discussed.

Altered substrate specificity of drug-resistant human immunodeficiency virus type 1 protease

Resistance to human immunodeficiency virus type 1 protease (HIV PR) inhibitors results primarily from the selection of multiple mutations in the protease region. Because many of these mutations are selected for the ability to decrease inhibitor binding in the active site, they also affect substrate binding and potentially substrate specificity. This work investigates the substrate specificity of a panel of clinically derived protease inhibitor-resistant HIV PR variants. To compare protease specificity, we have used positional-scanning, synthetic combinatorial peptide libraries as well as a select number of individual substrates. The subsite preferences of wild-type HIV PR determined by using the substrate libraries are consistent with prior reports, validating the use of these libraries to compare specificity among a panel of HIV PR variants. Five out of seven protease variants demonstrated subtle differences in specificity that may have significant impacts on their abilities to function in viral maturation. Of these, four variants demonstrated up to fourfold changes in the preference for valine relative to alanine at position P2 when tested on individual peptide substrates. This change correlated with a common mutation in the viral NC/p1 cleavage site. These mutations may represent a mechanism by which severely compromised, drug-resistant viral strains can increase fitness levels. Understanding the altered substrate specificity of drug-resistant HIV PR should be valuable in the design of future generations of protease inhibitors as well as in elucidating the molecular basis of regulation of proteolysis in HIV.

1.9 A x-ray study shows closed flap conformation in crystals of tethered HIV-1 PR

Three-dimensional structure of an asymmetrically mutated (C95M) tethered human immunodeficiency virus type 1 protease enzyme (HIV-1 PR) has been determined in an unliganded form using X-ray diffraction data to 1.9 A resolution. The structure, refined using X-PLOR to an R factor of 19.5%, is unexpectedly similar to the ligand-bound native enzyme, rather than to the ligand-free native enzyme. In particular, the two flaps in the tethered dimer are in a closed configuration. The environments around M95 and C1095 are identical, showing no structural effect of this asymmetric mutation at position 95. Oxidation of Cys1095 has been observed for the first time. There is one well-defined water molecule that hydrogen bonds to both carboxyl groups of the essential aspartic acids in the active site. Proteins 2001;43:57-64.

Alteration of substrate and inhibitor specificity of feline immunodeficiency virus protease

Feline immunodeficiency virus (FIV) protease is structurally very similar to human immunodeficiency virus (HIV) protease but exhibits distinct substrate and inhibitor specificities. We performed mutagenesis of subsite residues of FIV protease in order to define interactions that dictate this specificity. The I37V, N55M, M56I, V59I, and Q99V mutants yielded full activity. The I37V, N55M, V59I, and Q99V mutants showed a significant increase in activity against the HIV-1 reverse transcriptase/integrase and P2/nucleocapsid junction peptides compared with wild-type (wt) FIV protease. The I37V, V59I, and Q99V mutants also showed an increase in activity against two rapidly cleaved peptides selected by cleavage of a phage display library with HIV-1 protease. Mutations at Q54K, I98P, and L101I dramatically reduced activity. Mutants containing a I35D or I57G substitution showed no activity against either FIV or HIV substrates. FIV proteases all failed to cut HIV-1 matrix/capsid, P1/P6, P6/protease, and protease/reverse transcriptase junctions, indicating that none of the substitutions were sufficient to change the specificity completely. The I37V, N55M, M56I, V59I, and Q99V mutants, compared with wt FIV protease, all showed inhibitor specificity more similar to that of HIV-1 protease. The data also suggest that FIV protease prefers a hydrophobic P2/P2' residue like Val over Asn or Glu, which are utilized by HIV-1 protease, and that S2/S2' might play a critical role in distinguishing FIV and HIV-1 protease by specificity. The findings extend our observations regarding the interactions involved in substrate binding and aid in the development of broad-based inhibitors.