Analysis of HIV-1 CRF_01 A/E protease inhibitor resistance: structural determinants for maintaining sensitivity and developing resistance to atazanavir

A series of HIV-1 protease mutants has been designed in an effort to analyze the contribution to drug resistance provided by natural polymorphisms as well as therapy-selective (active and non-active site) mutations in the HIV-1 CRF_01 A/E (AE) protease when compared to that of the subtype B (B) protease. Kinetic analysis of these variants using chromogenic substrates showed differences in substrate specificity between pretherapy B and AE proteases. Inhibition analysis with ritonavir, indinavir, nelfinavir, amprenavir, saquinavir, lopinavir, and atazanavir revealed that the natural polymorphisms found in A/E can influence inhibitor resistance. It was also apparent that a high level of resistance in the A/E protease, as with B protease, is due to it aquiring a combination of active site and non-active site mutations. Structural analysis of atazanavir bound to a pretherapy B protease showed that the ability of atazanavir to maintain its binding affinity for variants containing some resistance mutations is due to its unique interactions with flap residues. This structure also explains why the I50L and I84V mutations are important in decreasing the binding affinity of atazanavir.     

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.     

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.     

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.     

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.     

Exploring the drug resistance mechanism of active site,non-active site mutations and their cooperative effects in CRF01_AE HIV-1 protease: molecular dynamics simulations and free energy calculations

Human immunodeficiency virus type 1 protease is essential for virus replication and maturation and has been considered as one of the important drug target for the antiretroviral treatment of HIV infection. The majority of HIV infections are caused due to non-B subtypes in developing countries. Subtype AE is spreading rapidly and infecting huge population worldwide. Understanding the interdependence of active and non-active site mutations in conferring drug resistance is crucial for the development effective inhibitors in subtype AE protease. In this work, we have investigated the mechanism of resistance against indinavir (IDV) due to therapy selected active site mutation V82F, non-active site mutations PF82V and their cooperative effects PV82F in subtype AE-protease using molecular dynamics simulations and binding free energy calculations. The simulations suggested all the three complexes lead to decrease in binding affinity of IDV, whereas the PF82V complex resulted in an enhanced binding affinity compared to V82F and PV82F complexes. Large positional deviation of IDV was observed in V82F complex. The preservation of hydrogen bonds of IDV with active site Asp25/Asp25′ and flap residue Ile50/50′ via a water molecule is crucial for effective binding. Owing to the close contact of 80s loop with Ile50′ and Asp25, the alteration between residues Thr80 and Val82, further induces conformational change thereby resulting in loss of interactions between IDV and the residues in the active site cavity, leading to drug resistance. Our present study shed light on the effect of active, non-active site mutations and their cooperative effects in AE protease.

Communicated by Ramaswamy H. Sarma     

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.     

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.     

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.     

Secondary mutations M36I and A71V in the human immunodeficiency virus type 1 protease can provide an advantage for the emergence of the primary mutation D30N

Development of resistance mutations in enzymatic targets of human immunodeficiency virus 1 (HIV-1) hampers the ability to provide adequate therapy. Of special interest is the effect mutations outside the active site of HIV-1 protease have on inhibitor binding and virus viability. We engineered protease mutants containing the active site mutation D30N alone and with the nonactive site polymorphisms M36I and/or A71V. We determined the K(i) values for the inhibitors nelfinavir, ritonavir, indinavir, KNI272, and AG1776 as well as the catalytic efficiency of the mutants. Single and double mutation combinations exhibited a decrease in catalytic efficiency, while the triple mutant displayed catalytic efficiency greater than that of the wild type. Variants containing M36I or A71V alone did not display a significant change in binding affinities to the inhibitors tested. The variant containing mutation D30N displayed a 2-6-fold increase in K(i) for all inhibitors tested, with nelfinavir showing the greatest increase. The double mutants containing a combination of mutations D30N, M36I, and A71V displayed -0.5-fold to +6-fold changes in the K(i) of all inhibitors tested, with ritonavir and nelfinavir most affected. Only the triple mutant showed a significant increase (>10-fold) in K(i) for inhibitor nelfinavir, ritonavir, or AG-1776 displaying 22-, 19-, or 15-fold increases, respectively. Our study shows that the M36I and A71V mutations provide a greater level of inhibitor cross-resistance combined with active site mutation D30N. M36I and A71V, when present as natural polymorphisms, could aid the virus in developing active site mutations to escape inhibitor binding while maintaining catalytic efficiency.     

Non-active site changes elicit broad-based cross-resistance of the HIV-1 protease to inhibitors.

Three high level, cross-resistant variants of the HIV-1 protease have been analyzed for their ability to bind four protease inhibitors approved by the Food and Drug Administration (saquinavir, ritonavir, indinavir, and nelfinavir) as AIDS therapeutics. The loss in binding energy (DeltaDeltaG(b)) going from the wild-type enzyme to mutant enzymes ranges from 2.5 to 4.4 kcal/mol, 40-65% of which is attributed to amino acid substitutions away from the active site of the protease and not in direct contact with the inhibitor. The data suggest that non-active site changes are collectively a major contributor toward engendering resistance against the protease inhibitor and cannot be ignored when considering cross-resistance issues of drugs against the HIV-1 protease.     

Resistance mechanism revealed by crystal structures of unliganded nelfinavir-resistant HIV-1 protease non-active site mutants N88D and N88S

Nelfinavir is an inhibitor of HIV-1 protease, and is used for treatment of patients suffering from HIV/AIDS. However, treatment results in drug resistant mutations in HIV-1 protease. N88D and N88S are two such mutations which occur in the non-active site region of the enzyme. We have determined crystal structures of unliganded N88D and N88S mutants of HIV-1 protease to resolution of 1.65 Å and 1.8 Å, respectively. These structures refined against synchrotron data lead to R-factors of 0.1859 and 0.1780, respectively. While structural effects of N88D are very subtle, the mutation N88S has caused a significant conformational change in D30, an active site residue crucial for substrate and inhibitor 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.     

Accessory mutations balance the marginal stability of the HIV-1 protease in drug resistance

The HIV-1 protease is a major target of inhibitor drugs in AIDS therapies. The therapies are impaired by mutations of the HIV-1 protease that can lead to resistance to protease inhibitors. These mutations are classified into major mutations, which usually occur first and clearly reduce the susceptibility to protease inhibitors, and minor, accessory mutations that occur later and individually do not substantially affect the susceptibility to inhibitors. Major mutations are predominantly located in the active site of the HIV-1 protease and can directly interfere with inhibitor binding. Minor mutations, in contrast, are typically located distal to the active site. A central question is how these distal mutations contribute to resistance development. In this article, we present a systematic computational investigation of stability changes caused by major and minor mutations of the HIV-1 protease. As most small single-domain proteins, the HIV-1 protease is only marginally stable. Mutations that destabilize the folded, active state of the protease therefore can shift the conformational equilibrium towards the unfolded, inactive state. We find that the most frequent major mutations destabilize the HIV-1 protease, whereas roughly half of the frequent minor mutations are stabilizing. An analysis of protease sequences from patients in treatment indicates that the stabilizing minor mutations are frequently correlated with destabilizing major mutations, and that highly resistant HIV-1 proteases exhibit significant fractions of stabilizing mutations. Our results thus indicate a central role of minor mutations in balancing the marginal stability of the protease against the destabilization induced by the most frequent major mutations.     

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.     

Structural implications of drug-resistant mutants of HIV-1 protease: high-resolution crystal structures of the mutant protease/substrate analogue complexes

Emergence of drug-resistant mutants of HIV-1 protease is an ongoing problem in the fight against AIDS. The mechanisms governing resistance are both complex and varied. We have determined crystal structures of HIV-1 protease mutants, D30N, K45I, N88D, and L90M complexed with peptide inhibitor analogues of CA-p2 and p2-NC cleavage sites in the Gag-pol precursor in order to study the structural mechanisms underlying resistance. The structures were determined at 1.55-1.9-A resolution and compared with the wild-type structure. The conformational disorder seen for most of the hydrophobic side-chains around the inhibitor binding site indicates flexibility of binding. Eight water molecules are conserved in all 9 structures; their location suggests that they are important for catalysis as well as structural stability. Structural differences among the mutants were analyzed in relation to the observed changes in protease activity and stability. Mutant L90M shows steric contacts with the catalytic Asp25 that could destabilize the catalytic loop at the dimer interface, leading to its observed decreased dimer stability and activity. Mutant K45I reduces the mobility of the flap and the inhibitor and contributes to an enhancement in structural stability and activity. The side-chain variations at residue 30 relative to wild-type are the largest in D30N and the changes are consistent with the altered activity observed with peptide substrates. Polar interactions in D30N are maintained, in agreement with the observed urea sensitivity. The side-chains of D30N and N88D are linked through a water molecule suggesting correlated changes at the two sites, as seen with clinical inhibitors. Structural changes seen in N88D are small; however, water molecules that mediate interactions between Asn88 and Thr74/Thr31/Asp30 in other complexes are missing in N88D.     

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.     

HIV-1 protease inhibitors: enthalpic versus entropic optimization of the binding affinity

Existing experimental as well as computational screening methods select potential ligands or drug candidates on the basis of binding affinity. Since the binding affinity is a function of the enthalpy (DeltaH) and entropy (DeltaS) changes, it is apparent that improved binding can be achieved in different ways: by optimizing DeltaH, DeltaS, or a combination of both. However, the behavior of enthalpically or entropically optimized inhibitors is fundamentally different, including their response to mutations that may elicit drug resistance. In the design of HIV-1 protease inhibitors, high binding affinity has usually been achieved by preshaping lead compounds to the geometry of the binding site and by incorporating a high degree of hydrophobicity. The thermodynamic consequence of that approach is that the binding affinity of the resulting inhibitors becomes entropically favorable but enthalpically unfavorable. Specifically, the resulting high binding affinity is due to an increased solvation entropy (hydrophobic effect) combined with a reduced loss of conformational entropy of the inhibitor upon binding (structural rigidity). Here we report that tripeptide inhibitors derived from the transframe region of Gag-Pol (Glu-Asp-Leu and Glu-Asp-Phe) bind to the HIV-1 protease with a favorable enthalpy change. This behavior is qualitatively different from that of known inhibitors and points to new strategies for inhibitor design. Since the binding affinities of enthalpically favorable and enthalpically unfavorable inhibitors have opposite temperature dependence, it is possible to design fast screening protocols that simultaneously select inhibitors on the basis of affinity and enthalpy.     

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

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