Nine Crystal Structures Determine the Substrate Envelope of the MDR HIV-1 Protease

Under drug selection pressure, emerging mutations render HIV-1 protease drug resistant, leading to the therapy failure in anti-HIV treatment. It is known that nine substrate cleavage site peptides bind to wild type (WT) HIV-1 protease in a conserved pattern. However, how the multidrug-resistant (MDR) HIV-1 protease binds to the substrate cleavage site peptides is yet to be determined. MDR769 HIV-1 protease (resistant mutations at residues 10, 36, 46, 54, 62, 63, 71, 82, 84, and 90) was selected for present study to understand the binding to its natural substrates. MDR769 HIV-1 protease was co-crystallized with nine substrate cleavage site hepta-peptides. Crystallographic studies show that MDR769 HIV-1 protease has an expanded substrate envelope with wide open flaps. Furthermore, ligand binding energy calculations indicate weaker binding in MDR769 HIV-1 protease-substrate complexes. These results help in designing the next generation of HIV-1 protease inhibitors by targeting the MDR HIV-1 protease.     

The higher barrier of darunavir and tipranavir resistance for HIV-1 protease

Darunavir and tipranavir are two inhibitors that are active against multi-drug resistant (MDR) HIV-1 protease variants. In this study, the invitro inhibitory efficacy was tested against a MDR HIV-1 protease variant, MDR 769 82T, containing the drug resistance mutations of 46L/54V/82T/84V/90M. Crystallographic and enzymatic studies were performed to examine the mechanism of resistance and the relative maintenance of potency. The key findings are as follows: (i) The MDR protease exhibits decreased susceptibility to all nine HIV-1 protease inhibitors approved by the US Food and Drug Administration (FDA), among which darunavir and tipranavir are the most potent; (ii) the threonine 82 mutation on the protease greatly enhances drug resistance by altering the hydrophobicity of the binding pocket; (iii) darunavir or tipranavir binding facilitates closure of the wide-open flaps of the MDR protease; and (iv) the remaining potency of tipranavir may be preserved by stabilizing the flaps in the inhibitor-protease complex while darunavir maintains its potency by preserving protein main chain hydrogen bonds with the flexible P2 group. These results could provide new insights into drug design strategies to overcome multi-drug resistance of HIV-1 protease variants.     

Crystal structures of multidrug-resistant HIV-1 protease in complex with two potent anti-malarial compounds

Two potent inhibitors (compounds 1 and 2) of malarial aspartyl protease, plasmepsin-II, were evaluated against wild type (NL4-3) and multidrug-resistant clinical isolate 769 (MDR) variants of human immunodeficiency virus type-1 (HIV-1) aspartyl protease. Enzyme inhibition assays showed that both 1 and 2 have better potency against NL4-3 than against MDR protease. Crystal structures of MDR protease in complex with 1 and 2 were solved and analyzed. Crystallographic analysis revealed that the MDR protease exhibits a typical wide-open conformation of the flaps (Gly48 to Gly52) causing an overall expansion in the active site cavity, which, in turn caused unstable binding of the inhibitors. Due to the expansion of the active site cavity, both compounds showed loss of direct contacts with the MDR protease compared to the docking models of NL4-3. Multiple water molecules showed a rich network of hydrogen bonds contributing to the stability of the ligand binding in the distorted binding pockets of the MDR protease in both crystal structures. Docking analysis of 1 and 2 showed a decrease in the binding affinity for both compounds against MDR supporting our structure-function studies. Thus, compounds 1 and 2 show promising inhibitory activity against HIV-1 protease variants and hence are good candidates for further development to enhance their potency against NL4-3 as well as MDR HIV-1 protease variants.     

Conserved hydrogen bonds and water molecules in MDR HIV-1 protease substrate complexes

The success of highly active antiretroviral therapy (HAART) in anti-HIV therapy is severely compromised by the rapidly developing drug resistance. HIV-1 protease inhibitors, part of HAART, are losing their potency and efficacy in inhibiting the target. Multi-drug resistant (MDR) 769 HIV-1 protease (resistant mutations at residues 10, 36, 46, 54, 62, 63, 71, 82, 84, 90) was selected for the present study to understand the binding to its natural substrates. The nine crystal structures of MDR769 HIV-1 protease substrate hepta-peptide complexes were analyzed in order to reveal the conserved structural elements for the purpose of drug design against MDR HIV-1 protease. Our structural studies demonstrated that highly conserved hydrogen bonds between the protease and substrate peptides, together with the conserved crystallographic water molecules, played a crucial role in the substrate recognition, substrate stabilization and protease stabilization. In addition, the absence of the key flap-ligand bridging water molecule might imply a different catalytic mechanism of MDR769 HIV-1 protease compared to that of wild type (WT) HIV-1 protease.     

The L33F darunavir resistance mutation acts as a molecular anchor reducing the flexibility of the HIV-1 protease 30s and 80s loops

HIV-1 protease (PR) is a 99 amino acid protein responsible for proteolytic processing of the viral polyprotein – an essential step in the HIV-1 life cycle. Drug resistance mutations in PR that are selected during antiretroviral therapy lead to reduced efficacy of protease inhibitors (PI) including darunavir (DRV). To identify the structural mechanisms associated with the DRV resistance mutation L33F, we performed X-ray crystallographic studies with a multi-drug resistant HIV-1 protease isolate that contains the L33F mutation (MDR769 L33F). In contrast to other PR L33F DRV complexes, the structure of MDR769 L33F complexed with DRV reported here displays the protease flaps in an open conformation. The L33F mutation increases noncovalent interactions in the hydrophobic pocket of the PR compared to the wild-type (WT) structure. As a result, L33F appears to act as a molecular anchor, reducing the flexibility of the 30s loop (residues 29–35) and the 80s loop (residues 79–84). Molecular anchoring of the 30s and 80s loops leaves an open S1/S1′ subsite and distorts the conserved hydrogen-bonding network of DRV. These findings are consistent with previous reports despite structural differences with regards to flap conformation.     

Structural and thermodynamic basis for the binding of TMC114, a next-generation human immunodeficiency virus type 1 protease inhibitor

TMC114, a newly designed human immunodeficiency virus type 1 (HIV-1) protease inhibitor, is extremely potent against both wild-type (wt) and multidrug-resistant (MDR) viruses in vitro as well as in vivo. Although chemically similar to amprenavir (APV), the potency of TMC114 is substantially greater. To examine the basis for this potency, we solved crystal structures of TMC114 complexed with wt HIV-1 protease and TMC114 and APV complexed with an MDR (L63P, V82T, and I84V) protease variant. In addition, we determined the corresponding binding thermodynamics by isothermal titration calorimetry. TMC114 binds approximately 2 orders of magnitude more tightly to the wt enzyme (K(d) = 4.5 x 10(-12) M) than APV (K(d) = 3.9 x 10(-10) M). Our X-ray data (resolution ranging from 2.2 to 1.2 A) reveal strong interactions between the bis-tetrahydrofuranyl urethane moiety of TMC114 and main-chain atoms of D29 and D30. These interactions appear largely responsible for TMC114's very favorable binding enthalpy to the wt protease (-12.1 kcal/mol). However, TMC114 binding to the MDR HIV-1 protease is reduced by a factor of 13.3, whereas the APV binding constant is reduced only by a factor of 5.1. However, even with the reduction in binding affinity to the MDR HIV protease, TMC114 still binds with an affinity that is more than 1.5 orders of magnitude tighter than the first-generation inhibitors. Both APV and TMC114 fit predominantly within the substrate envelope, a property that may be associated with decreased susceptibility to drug-resistant mutations relative to that of first-generation inhibitors. Overall, TMC114's potency against MDR viruses is likely a combination of its extremely high affinity and close fit within the substrate envelope.     

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.     

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.     

Thermodynamic basis of resistance to HIV-1 protease inhibition: calorimetric analysis of the V82F/I84V active site resistant mutant

One of the most serious side effects associated with the therapy of HIV-1 infection is the appearance of viral strains that exhibit resistance to protease inhibitors. The active site mutant V82F/I84V has been shown to lower the binding affinity of protease inhibitors in clinical use. To identify the origin of this effect, we have investigated the binding thermodynamics of the protease inhibitors indinavir, ritonavir, saquinavir, and nelfinavir to the wild-type HIV-1 protease and to the V82F/I84V resistant mutant. The main driving force for the binding of all four inhibitors is a large positive entropy change originating from the burial of a significant hydrophobic surface upon binding. At 25 degrees C, the binding enthalpy is unfavorable for all inhibitors except ritonavir, for which it is slightly favorable (-2.3 kcal/mol). Since the inhibitors are preshaped to the geometry of the binding site, their conformational entropy loss upon binding is small, a property that contributes to their high binding affinity. The V82F/I84V active site mutation lowers the affinity of the inhibitors by making the binding enthalpy more positive and making the entropy change slightly less favorable. The effect on the enthalpy change is, however, the major one. The predominantly enthalpic effect of the V82F/I84V mutation is consistent with the idea that the introduction of the bulkier Phe side chain at position 82 and the Val side chain at position 84 distort the binding site and weaken van der Waals and other favorable interactions with inhibitors preshaped to the wild-type binding site. Another contribution of the V82F/I84V to binding affinity originates from an increase in the energy penalty associated with the conformational change of the protease upon binding. The V82F/I84V mutant is structurally more stable than the wild-type protease by about 1.4 kcal/mol. This effect, however, affects equally the binding affinity of substrate and inhibitors.     

Binding kinetics of darunavir to human immunodeficiency virus type 1 protease explain the potent antiviral activity and high genetic barrier

The high incidence of cross-resistance between human immunodeficiency virus type 1 (HIV-1) protease inhibitors (PIs) limits their sequential use. This necessitates the development of PIs with a high genetic barrier and a broad spectrum of activity against PI-resistant HIV, such as tipranavir and darunavir (TMC114). We performed a surface plasmon resonance-based kinetic study to investigate the impact of PI resistance-associated mutations on the protease binding of five PIs used clinically: amprenavir, atazanavir, darunavir, lopinavir, and tipranavir. With wild-type protease, the binding affinity of darunavir was more than 100-fold higher than with the other PIs, due to a very slow dissociation rate. Consequently, the dissociative half-life of darunavir was much higher (>240 h) than that of the other PIs, including darunavir's structural analogue amprenavir. The influence of protease mutations on the binding kinetics was tested with five multidrug-resistant (MDR) proteases derived from clinical isolates harboring 10 to 14 PI resistance-associated mutations with a decreased susceptibility to various PIs. In general, all PIs bound to the MDR proteases with lower binding affinities, caused mainly by a faster dissociation rate. For amprenavir, atazanavir, lopinavir, and tipranavir, the decrease in affinity with MDR proteases resulted in reduced antiviral activity. For darunavir, however, a nearly 1,000-fold decrease in binding affinity did not translate into a weaker antiviral activity; a further decrease in affinity was required for the reduced antiviral effect. These observations provide a mechanistic explanation for darunavir's potent antiviral activity and high genetic barrier to the development of resistance.     

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.     

Lack of synergy for inhibitors targeting a multi-drug-resistant HIV-1 protease

The three-dimensional structures of indinavir and three newly synthesized indinavir analogs in complex with a multi-drug-resistant variant (L63P, V82T, I84V) of HIV-1 protease were determined to approximately 2.2 A resolution. Two of the three analogs have only a single modification of indinavir, and their binding affinities to the variant HIV-1 protease are enhanced over that of indinavir. However, when both modifications were combined into a single compound, the binding affinity to the protease variant was reduced. On close examination, the structural rearrangements in the protease that occur in the tightest binding inhibitor complex are mutually exclusive with the structural rearrangements seen in the second tightest inhibitor complex. This occurs as adaptations in the S1 pocket of one monomer propagate through the dimer and affect the conformation of the S1 loop near P81 of the other monomer. Therefore, structural rearrangements that occur within the protease when it binds to an inhibitor with a single modification must be accounted for in the design of inhibitors with multiple modifications. This consideration is necessary to develop inhibitors that bind sufficiently tightly to drug-resistant variants of HIV-1 protease to potentially become the next generation of therapeutic agents.     

Ligand modifications to reduce the relative resistance of multi-drug resistant HIV-1 protease

Proper proteolytic processing of the HIV-1 Gag/Pol polyprotein is required for HIV infection and viral replication. This feature has made HIV-1 protease an attractive target for antiretroviral drug design for the treatment of HIV-1 infected patients. To examine the role of the P1 and P1′positions of the substrate in inhibitory efficacy of multi-drug resistant HIV-1 protease 769 (MDR 769), we performed a series of structure–function studies. Using the original CA/p2 cleavage site sequence, we generated heptapeptides containing one reduced peptide bond with an L to F and A to F double mutation at P1 and P1′ (F-r-F), and an A to F at P1′ (L-r-F) resulting in P1/P1′ modified ligands. Here, we present an analysis of co-crystal structures of CA/p2 F-r-F, and CA/p2 L-r-F in complex with MDR 769. To examine conformational changes in the complex structure, molecular dynamic (MD) simulations were performed with MDR769–ligand complexes. MD trajectories show the isobutyl group of both the lopinavir analog and the CA/p2 L-r-F substrate cause a conformational change of in the active site of MDR 769. IC₅₀ measurements suggest the non identical P1/P1′ ligands (CA/p2 L-r-F and lopinavir analog) are more effective against MDR proteases as opposed to identical P1/P1′ligands. Our results suggest that a non identical P1/P1′composition may be more favorable for the inhibition of MDR 769 as they induce conformational changes in the active site of the enzyme resulting in disruption of the two-fold symmetry of the protease, thus, stabilizing the inhibitor in the active site.     

Active-site mutations in the South african human immunodeficiency virus type 1 subtype C protease have a significant impact on clinical inhibitor binding: kinetic and thermodynamic study

Human immunodeficiency virus (HIV) infections in sub-Saharan Africa represent about 56% of global infections. Study of active-site mutations (the V82A single mutation and the V82F I84V double mutation) in the less-studied South African HIV type 1 subtype C (C-SA) protease indicated that neither mutation had a significant impact on the proteolytic functioning of the protease. However, the binding affinities of, and inhibition by, saquinavir, ritonavir, indinavir, and nelfinavir were weaker for each variant than for the wild-type protease, with the double mutant exhibiting the most dramatic change. Therefore, our results show that the C-SA V82F I84V double mutation decreased the binding affinities of protease inhibitors to levels significantly lower than that required for effective inhibition.     

Amplification of the effects of drug resistance mutations by background polymorphisms in HIV-1 protease from African subtypes

The vast majority of HIV-1 infections worldwide are caused by the C and A viral subtypes rather than the B subtype prevalent in the United States and Western Europe. Genomic differences between subtypes give rise to sequence variations in the encoded proteins, including those identified as targets for antiretroviral therapies. In the case of the HIV-1 protease, we reported earlier [Velazquez-Campoy et al. (2001) Proc. Natl. Acad. Sci. U.S.A. 98, 6062-6067] that proteases from the C and A subtypes exhibit a higher biochemical fitness in the presence of widely prescribed protease inhibitors. In this paper we present a complete thermodynamic dissection of the differences between proteases from different subtypes and the effects of the V82F/I84V drug-resistant mutation within the framework of the B, C, and A subtypes. These studies involved four inhibitors in clinical use (indinavir, saquinavir, ritonavir, and nelfinavir) and a second-generation protease inhibitor (KNI-764). Naturally occurring amino acid polymorphisms found in proteases from the C and A subtypes lower the binding affinities of existing clinical inhibitors by factors ranging between 2 and 7.5 which by themselves are not enough to cause drug resistance. The preexisting lower affinity in the C and A subtypes, however, significantly amplifies the effects of the drug-resistant mutation. Relative to the wild-type B subtype protease, the V82F/I84V drug-resistant mutation within the C and A subtypes lowers the binding affinity of inhibitors by factors ranging between 40 and 3000. When the enzyme kinetic properties (k(cat) and K(m)) are included in the analysis, the biochemical fitness of the C and A subtype drug-resistant mutants can be up to 1000-fold higher than that of the wild-type B subtype protease in the presence of the studied inhibitors. From a thermodynamic standpoint, the combined effects of the drug-resistant mutations and the natural amino acid polymorphisms on the Gibbs energy are additive and involve significant alterations in the enthalpy and entropy changes associated with inhibitor binding. At the biochemical level, the combined effects of naturally existing polymorphisms and drug-resistant mutations might have important consequences on the long-term viability of current HIV-1 protease inhibitors.     

The binding energetics of first- and second-generation HIV-1 protease inhibitors: implications for drug design

KNI-764 is a powerful HIV-1 protease inhibitor with a reported low susceptibility to the effects of protease mutations commonly associated with drug resistance. In this paper the binding thermodynamics of KNI-764 to the wild-type and drug-resistant mutant V82F/I84V are presented and the results compared to those obtained with existing clinical inhibitors. KNI-764 binds to the wild-type HIV-1 protease with very high affinity (3.1 x 10(10) M(-1) or 32 pM) in a process strongly favored by both enthalpic and entropic contributions to the Gibbs energy of binding (Delta G = -RTlnK(a)). When compared to existing clinical inhibitors, the binding affinity of KNI-764 is about 100 fold higher than that of indinavir, saquinavir, and nelfinavir, but comparable to that of ritonavir. Unlike the existing clinical inhibitors, which bind to the protease with unfavorable or only slightly favorable enthalpy changes, the binding of KNI-764 is strongly exothermic (-7.6 kcal/mol). The resistant mutation V82F/I84V lowers the binding affinity of KNI-764 26-fold, which can be accounted almost entirely by a less favorable binding enthalpy to the mutant. Since KNI-764 binds to the wild type with extremely high affinity, even after a 26-fold decrease, it still binds to the resistant mutant with an affinity comparable to that of other inhibitors against the wild type. These results indicate that the effectiveness of this inhibitor against the resistant mutant is related to two factors: extremely high affinity against the wild type achieved by combining favorable enthalpic and entropic interactions, and a mild effect of the protease mutation due to the presence of flexible structural elements at critical locations in the inhibitor molecule. The conclusions derived from the HIV-1 protease provide important thermodynamic guidelines that can be implemented in general drug design strategies.     

"Wide-open" 1.3 A structure of a multidrug-resistant HIV-1 protease as a drug target

This report examines structural changes in a highly mutated, clinical multidrug-resistant HIV-1 protease, and the crystal structure has been solved to 1.3 A resolution in the absence of any inhibitor. This protease variant contains codon mutations at positions 10, 36, 46, 54, 62, 63, 71, 82, 84, and 90 that confer resistance to protease inhibitors. Major differences between the wild-type and the variant include a structural change initiated by the M36V mutation and amplified by additional mutations in the flaps of the protease, resulting in a "wide-open" structure that represents an opening that is 8 A wider than the "open" structure of the wild-type protease. A second structural change is triggered by the L90M mutation that results in reshaping the 23-32 segment. A third key structural change of the protease is due to the mutations from longer to shorter amino acid side chains at positions 82 and 84.     

A structural and thermodynamic escape mechanism from a drug resistant mutation of the HIV-1 protease

The efficacy of HIV-1 protease inhibition therapies is often compromised by the appearance of mutations in the protease molecule that lower the binding affinity of inhibitors while maintaining viable catalytic activity and substrate affinity. The V82F/I84V double mutation is located within the binding site cavity and affects all protease inhibitors in clinical use. KNI-764, a second-generation inhibitor currently under development, maintains significant potency against this mutation by entropically compensating for enthalpic losses, thus minimizing the loss in binding affinity. KNI-577 differs from KNI-764 by a single functional group critical to the inhibitor response to the protease mutation. This single difference changes the response of the two inhibitors to the mutation by one order of magnitude. Accordingly, a structural understanding of the inhibitor response will provide important guidelines for the design of inhibitors that are less susceptible to mutations conveying drug resistance. The structures of the two compounds bound to the wild type and V82F/I84V HIV-1 protease have been determined by X-ray crystallography at 2.0 A resolution. The presence of two asymmetric functional groups, linked by rotatable bonds to the inhibitor scaffold, allows KNI-764 to adapt to the mutated binding site cavity more readily than KNI-577, with a single asymmetric group. Both inhibitors lose about 2.5 kcal/mol in binding enthalpy when facing the drug-resistant mutant protease; however KNI-764 gains binding entropy while KNI-577 loses binding entropy. The gain in binding entropy by KNI-764 accounts for its low susceptibility to the drug-resistant mutation. The heat capacity change associated with binding becomes more negative when KNI-764 binds to the mutant protease, consistent with increased desolvation. With KNI-577, the opposite effect is observed. Structurally, the crystallographic B factors increase for KNI-764 when it is bound to the drug-resistant mutant. The opposite is observed for KNI-577. Consistent with these observations, it appears that KNI-764 is able to gain binding entropy by a two-fold mechanism: it gains solvation entropy by burying itself deeper within the binding pocket and gains conformational entropy by losing interaction with the protease.     

Why is Substrate Peptide Binding Unsusceptible to Multidrug-Resistant Mutations in HIV-1 Protease? A Structural and Energetic Analysis

Understanding of the molecular mechanism and biological implication underlying the difference in binding of substrate peptides and small-molecule inhibitors to multidrug-resistant mutants of HIV-1 protease would help to develop new anti-HIV agents combating drug resistance. Here, an integration of rigorous quantum mechanics/molecular mechanics (QM/MM) analysis and empirical Poisson–Boltzmann/surface area (PB/SA) model is described to investigate the structural basis and energetic property of wild-type HIV-1 protease and its mutants in recognizing and binding with a wide variety of ligands, including the peptides derived from its cognate cleavage sites and the cleavage site variants as well as a number of FDA-approved protease inhibitors, attempting to explain why is substrate binding unsusceptible to most observed HIV-1 protease mutations. A preliminary test study demonstrates that the combined QM/MM–PB/SA scheme is able to effectively reproduce the relative ligand binding energy changes upon protease single- and double-mutations, albeit the absolute values appear to be different significantly between the calculated and experimental results. With the QM/MM–PB/SA calculations a complete mutation energy map of HIV-1 protease–ligand interactions is created, which unravels distinct affinity pictures of wild-type substrates, substrate variants and, particularly, the protease inhibitors bound to HIV-1 protease mutants, suggesting that, on the one hand, the evaluation pressure under anti-HIV chemotherapies addresses site-directed protease mutations that impair and undermine the intermolecular interactions specific to inhibitors but not substrates; on the other hand, co-evaluation of protease and its substrate peptides provides a more effective mechanism to avoid therapeutic surveillance. Further, nonbonded interaction analysis and computational alanine scanning reveal 12 key residues that is critical for substrate binding, from which the Asn25, Gly27, Ala28, Asp29 and Pro81 are identified that have not yet been found to cause drug resistance and hence would be the promising sites targeted by new protease inhibitors.     

Structure at 2.5-A resolution of chemically synthesized human immunodeficiency virus type 1 protease complexed with a hydroxyethylene-based inhibitor

The crystal structure of a complex between chemically synthesized human immunodeficiency virus type 1 (HIV-1) protease and an octapeptide inhibitor has been refined to an R factor of 0.138 at 2.5-A resolution. The substrate-based inhibitor, H-Val-Ser-Gln-Asn-Leu psi [CH(OH)CH2]Val-Ile-Val-OH (U-85548e) contains a hydroxyethylene isostere replacement at the scissile bond that is believed to mimic the tetrahedral transition state of the proteolytic reaction. This potent inhibitor has Ki less than 1 nM and was developed as an active-site titrant of the HIV-1 protease. The inhibitor binds in an extended conformation and is involved in beta-sheet interactions with the active-site floor and flaps of the enzyme, which form the substrate/inhibitor cavity. The inhibitor diastereomer has the S configuration at the chiral carbon atom of the hydroxyethylene insert, and the hydroxyl group is within H-bonding distance of the two active-site carboxyl groups in the enzyme dimer. The two subunits of the enzyme are related by a pseudodyad, which superposes them at a 178 degrees rotation. The main difference between the subunits is in the beta turns of the flaps, which have different conformations in the two monomers. The inhibitor has a clear preferred orientation in the active site and the alternative conformation, if any, is a minor one (occupancy of less than 30%). A new model of the enzymatic mechanism is proposed in which the proteolytic reaction is viewed as a one-step process during which the nucleophile (water molecule) and electrophile (an acidic proton) attack the scissile bond in a concerted manner.