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.     

Crystal structures of a multidrug-resistant human immunodeficiency virus type 1 protease reveal an expanded active-site cavity

The goal of this study was to use X-ray crystallography to investigate the structural basis of resistance to human immunodeficiency virus type 1 (HIV-1) protease inhibitors. We overexpressed, purified, and crystallized a multidrug-resistant (MDR) HIV-1 protease enzyme derived from a patient failing on several protease inhibitor-containing regimens. This HIV-1 variant contained codon mutations at positions 10, 36, 46, 54, 63, 71, 82, 84, and 90 that confer drug resistance to protease inhibitors. The 1.8-angstrom (A) crystal structure of this MDR patient isolate reveals an expanded active-site cavity. The active-site expansion includes position 82 and 84 mutations due to the alterations in the amino acid side chains from longer to shorter (e.g., V82A and I84V). The MDR isolate 769 protease "flaps" stay open wider, and the difference in the flap tip distances in the MDR 769 variant is 12 A. The MDR 769 protease crystal complexes with lopinavir and DMP450 reveal completely different binding modes. The network of interactions between the ligands and the MDR 769 protease is completely different from that seen with the wild-type protease-ligand complexes. The water molecule-forming hydrogen bonds bridging between the two flaps and either the substrate or the peptide-based inhibitor are lacking in the MDR 769 clinical isolate. The S1, S1', S3, and S3' pockets show expansion and conformational change. Surface plasmon resonance measurements with the MDR 769 protease indicate higher k(off) rates, resulting in a change of binding affinity. Surface plasmon resonance measurements provide k(on) and k(off) data (K(d) = k(off)/k(on)) to measure binding of the multidrug-resistant protease to various ligands. This MDR 769 protease represents a new antiviral target, presenting the possibility of designing novel inhibitors with activity against the open and expanded protease forms.     

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.     

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.     

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.     

An in silico pharmacological approach toward the discovery of potent inhibitors to combat drug resistance HIV-1 protease variants

Protease inhibitors (PIs) are crucial drugs in highly active antiretroviral therapy for human immunodeficiency virus-1 (HIV-1) infections. However, resistance owing to mutations challenge the long-term efficacy in the medication of HIV-1-infected individuals. Lopinavir (LPV) and darunavir (DRV), two second-generation drugs are the most potent among PIs, hustling the drug resistance when mutations occur in the active and nonactive site of the protease (PR). Herein, we strive for compounds that can stifle the function of wild-type (WT) HIV-1 PR along with four major single mutants (I54M, V82T, I84V, and L90M) instigating resistance to the PIs using in silico approach. Six common compounds are retrieved from six databases using combined pharmacophore-based and structure-based virtual screening methodology. LPV and DRV are docked and the binding free energy is calculated to set the cut-off value for selecting compounds. Further, to gain insight into the stability of the complexes the molecular dynamics simulation (MDS) is carried out, which uncovers two lead molecules namely NCI-524545 and ZINC12866729. Both the lead molecules connect with WT and mutant HIV-1 PRs through strong and stable hydrogen bond interactions when compared with LPV and DRV throughout the trajectory analysis. Interestingly, NCI-524545 and ZINC12866729 exhibit direct interactions with I50/50′ by replacing the conserved water molecule as evidenced by MDS, which indicates the credible potency of these compounds. Hence, we concluded that NCI-524545 and ZINC12866729 have great puissant to restrain the role of drug resistance HIV-1 PR variants, which can also show better activity through in vivo and in vitro conditions.     

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 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.     

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.     

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.     

Genotypic changes in human immunodeficiency virus type 1 protease associated with reduced susceptibility and virologic response to the protease inhibitor tipranavir

Tipranavir is a novel, nonpeptidic protease inhibitor of human immunodeficiency virus type 1 (HIV-1) with activity against clinical HIV-1 isolates from treatment-experienced patients. HIV-1 genotypic and phenotypic data from phase II and III clinical trials of tipranavir with protease inhibitor-experienced patients were analyzed to determine the association of protease mutations with reduced susceptibility and virologic response to tipranavir. Specific protease mutations were identified based on stepwise multiple-regression analyses of phase II study data sets. Validation included analyses of phase III study data sets to determine if the same mutations would be selected and to assess how these mutations contribute to multiple-regression models of tipranavir-related phenotype and of virologic response. A tipranavir mutation score was developed from these analyses, which consisted of a unique string of 16 protease positions and 21 mutations (10V, 13V, 20M/R/V, 33F, 35G, 36I, 43T, 46L, 47V, 54A/M/V, 58E, 69K, 74P, 82L/T, 83D, and 84V). HIV-1 isolates displaying an increasing number of these tipranavir resistance-associated mutations had a reduced phenotypic susceptibility and virologic response to tipranavir. Regression models for predicting virologic response in phase III trials revealed that each point in the tipranavir score was associated with a 0.16-log10 copies/ml-lower virologic response to tipranavir at week 24 of treatment. A lower number of points in the tipranavir score and a greater number of active drugs in the background regimen were predictive of virologic success. These analyses demonstrate that the tipranavir mutation score is a potentially valuable tool for predicting the virologic response to tipranavir in protease inhibitor-experienced patients.     

X-ray crystal structures of human immunodeficiency virus type 1 protease mutants complexed with atazanavir

Atazanavir, which is marketed as REYATAZ, is the first human immunodeficiency virus type 1 (HIV-1) protease inhibitor approved for once-daily administration. As previously reported, atazanavir offers improved inhibitory profiles against several common variants of HIV-1 protease over those of the other peptidomimetic inhibitors currently on the market. This work describes the X-ray crystal structures of complexes of atazanavir with two HIV-1 protease variants, namely, (i) an enzyme optimized for resistance to autolysis and oxidation, referred to as the cleavage-resistant mutant (CRM); and (ii) the M46I/V82F/I84V/L90M mutant of the CRM enzyme, which is resistant to all approved HIV-1 protease inhibitors, referred to as the inhibitor-resistant mutant. In these two complexes, atazanavir adopts distinct bound conformations in response to the V82F substitution, which may explain why this substitution, at least in isolation, has yet to be selected in vitro or in the clinic. Because of its nearly symmetrical chemical structure, atazanavir is able to make several analogous contacts with each monomer of the biological dimer.     

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.     

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.     

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.     

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.     

Ultra-high resolution crystal structure of HIV-1 protease mutant reveals two binding sites for clinical inhibitor TMC114

TMC114 (darunavir) is a promising clinical inhibitor of HIV-1 protease (PR) for treatment of drug resistant HIV/AIDS. We report the ultra-high 0.84 A resolution crystal structure of the TMC114 complex with PR containing the drug-resistant mutation V32I (PR(V32I)), and the 1.22 A resolution structure of a complex with PR(M46L). These structures show TMC114 bound at two distinct sites, one in the active-site cavity and the second on the surface of one of the flexible flaps in the PR dimer. Remarkably, TMC114 binds at these two sites simultaneously in two diastereomers related by inversion of the sulfonamide nitrogen. Moreover, the flap site is shaped to accommodate the diastereomer with the S-enantiomeric nitrogen rather than the one with the R-enantiomeric nitrogen. The existence of the second binding site and two diastereomers suggest a mechanism for the high effectiveness of TMC114 on drug-resistant HIV and the potential design of new inhibitors.     

Multi-drug resistance profile of PR20 HIV-1 protease is attributed to distorted conformational and drug binding landscape: molecular dynamics insights

The PR20 HIV-1 protease, a variant with 20 mutations, exhibits high levels of multi-drug resistance; however, to date, there has been no report detailing the impact of these 20 mutations on the conformational and drug binding landscape at a molecular level. In this report, we demonstrate the first account of a comprehensive study designed to elaborate on the impact of these mutations on the dynamic features as well as drug binding and resistance profile, using extensive molecular dynamics analyses. Comparative MD simulations for the wild-type and PR20 HIV proteases, starting from bound and unbound conformations in each case, were performed. Results showed that the apo conformation of the PR20 variant of the HIV protease displayed a tendency to remain in the open conformation for a longer period of time when compared to the wild type. This led to a phenomena in which the inhibitor seated at the active site of PR20 tends to diffuse away from the binding site leading to a significant change in inhibitor–protein association. Calculating the per-residue fluctuation (RMSF) and radius of gyration, further validated these findings. MM/GBSA showed that the occurrence of 20 mutations led to a drop in the calculated binding free energies (ΔG_bind) by ~25.17 kcal/mol and ~5 kcal/mol for p2-NC, a natural peptide substrate, and darunavir, respectively, when compared to wild type. Furthermore, the residue interaction network showed a diminished inter-residue hydrogen bond network and changes in inter-residue connections as a result of these mutations. The increased conformational flexibility in PR20 as a result of loss of intra- and inter-molecular hydrogen bond interactions and other prominent binding forces led to a loss of protease grip on ligand. It is interesting to note that the difference in conformational flexibility between PR20 and WT conformations was much higher in the case of substrate-bound conformation as compared to DRV. Thus, developing analogues of DRV by retaining its key pharmacophore features will be the way forward in the search for novel protease inhibitors against multi-drug resistant strains.     

Structural evidence for effectiveness of darunavir and two related antiviral inhibitors against HIV-2 protease

No drug has been targeted specifically for HIV-2 (human immunodeficiency virus type 2) infection despite its increasing prevalence worldwide. The antiviral HIV-1 (human immunodeficiency virus type 1) protease (PR) inhibitor darunavir and the chemically related GRL98065 and GRL06579A were designed with the same chemical scaffold and different substituents at P2 and P2′ to optimize polar interactions for HIV-1 PR (PR1). These inhibitors are also effective antiviral agents for HIV-2-infected cells. Therefore, crystal structures of HIV-2 PR (PR2) complexes with the three inhibitors have been solved at 1.2-Å resolution to analyze the molecular basis for their antiviral potency. Unusually, the crystals were grown in imidazole and zinc acetate buffer, which formed interactions with the PR2 and the inhibitors. Overall, the structures were very similar to the corresponding inhibitor complexes of PR1 with an RMSD of 1.1 Å on main-chain atoms. Most hydrogen-bond and weaker C-H…O interactions with inhibitors were conserved in the PR2 and PR1 complexes, except for small changes in interactions with water or disordered side chains. Small differences were observed in the hydrophobic contacts for the darunavir complexes, in agreement with relative inhibition of the two PRs. These near-atomic-resolution crystal structures verify the inhibitor potency for PR1 and PR2 and will provide the basis for the development of antiviral inhibitors targeting PR2.