Effectiveness of commercial inhibitors against subtype F HIV-1 protease

Subtype F wild type HIV protease has been kinetically characterized using six commercial inhibitors (amprenavir, indinavir, lopinavir, nelfinavir, ritonavir and saquinavir) commonly used for HIV/AIDS treatment, as well as inhibitor TL-3 and acetyl-pepstatin. We also obtained kinetic parameters for two multi-resistant proteases (one of subtype B and one of subtype F) harboring primary and secondary mutations selected by intensive treatment with ritonavir/nelfinavir. This newly obtained biochemical data shows that all six studied commercially available protease inhibitors are significantly less effective against subtype F HIV proteases than against HIV proteases of subtype B, as judged by increased K(i) and biochemical fitness (vitality) values. Comparison with previously reported kinetic values for subtype A and C HIV proteases show that subtype F wild type proteases are significantly less susceptible to inhibition. These results demonstrate that the accumulation of natural polymorphisms in subtype F proteases yields catalytically more active enzymes with a large degree of cross-resistance, which thus results in strong virus viability.     

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 crystal structure of protease Sapp1p from Candida parapsilosis in complex with the HIV protease inhibitor ritonavir

Secreted aspartic proteases (Saps) are extracellular proteolytic enzymes that enhance the virulence of Candida pathogens. These enzymes therefore represent possible targets for therapeutic drug design. Saps are inhibited by nanomolar concentrations of the classical inhibitor of aspartic proteases pepstatin A and also by the inhibitors of the HIV protease, but with the K_i of micromolar values or higher. To contribute to the discussion regarding whether HIV protease inhibitors can act against opportunistic mycoses by the inhibition of Saps, we determined the structure of Sapp1p from Candida parapsilosis in complex with ritonavir (RTV), a clinically used inhibitor of the HIV protease. The crystal structure refined at resolution 2.4 Å proved binding of RTV into the active site of Sapp1p and provided the structural information necessary to evaluate the stability and specificity of the protein-inhibitor interaction.     

The contribution of naturally occurring polymorphisms in altering the biochemical and structural characteristics of HIV-1 subtype C protease

Fourteen subtype B and C protease variants have been engineered in an effort to study whether the preexistent baseline polymorphisms, by themselves or in combination with drug resistance mutations, differentially alter the biochemical and structural features of the subtype C protease when compared with those of subtype B protease. The kinetic studies performed in this work showed that the preexistent polymorphisms in subtype C protease, by themselves, do not provide for a greater level of resistance. Inhibition analysis with eight clinically used protease inhibitors revealed that the natural polymorphisms found in subtype C protease, in combination with drug resistance mutations, can influence enzymatic catalytic efficiency and inhibitor resistance. Structural analyses of the subtype C protease bound to nelfinavir and indinavir showed that these inhibitors form similar interactions with the residues in the active site of subtype B and C proteases. It also revealed that the naturally occurring polymorphisms could alter the position of the outer loops of the subtype C protease, especially the 60's loop.     

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.     

Structural characterization of B and non-B subtypes of HIV-protease: insights into the natural susceptibility to drug resistance development

Although a majority of HIV-1 infections in Brazil are caused by the subtype B virus (also prevalent in the United States and Western Europe), viral subtypes F and C are also found very frequently. Genomic differences between the subtypes give rise to sequence variations in the encoded proteins, including the HIV-1 protease. The current anti-HIV drugs have been developed primarily against subtype B and the effects arising from the combination of drug-resistance mutations with the naturally existing polymorphisms in non-B HIV-1 subtypes are only beginning to be elucidated. To gain more insights into the structure and function of different variants of HIV proteases, we have determined a 2.1 A structure of the native subtype F HIV-1 protease (PR) in complex with the protease inhibitor TL-3. We have also solved crystal structures of two multi-drug resistant mutant HIV PRs in complex with TL-3, from subtype B (Bmut) carrying the primary mutations V82A and L90M, and from subtype F (Fmut) carrying the primary mutation V82A plus the secondary mutation M36I, at 1.75 A and 2.8 A resolution, respectively. The proteases Bmut, Fwt and Fmut exhibit sevenfold, threefold, and 54-fold resistance to TL-3, respectively. In addition, the structure of subtype B wild type HIV-PR in complex with TL-3 has been redetermined in space group P6(1), consistent with the other three structures. Our results show that the primary mutation V82A causes the known effect of collapsing the S1/S1' pockets that ultimately lead to the reduced inhibitory effect of TL-3. Our results further indicate that two naturally occurring polymorphic substitutions in subtype F and other non-B HIV proteases, M36I and L89M, may lead to early development of drug resistance in patients infected with non-B HIV subtypes.     

Structural and kinetic analysis of pyrrolidine-based inhibitors of the drug-resistant Ile84Val mutant of HIV-1 protease

Human immunodeficiency virus (HIV) protease is a well-established drug target in HIV chemotherapy. However, continuously increasing resistance towards approved drugs inevitably requires the development of new inhibitors preferably showing no susceptibility against resistant HIV protease strains. Recently, symmetric pyrrolidine-3,4-bis-N-benzyl-sulfonamides have been developed as a new class of HIV-1 protease inhibitors. The most promising candidate exhibited a K_i of 74 nM towards a wild-type protease. Herein, we report the influence of the active-site mutations Ile50Val and Ile84Val on these inhibitors by structural and kinetic analysis. Although the Ile50Val mutation leads to a significant decrease in affinity for all compounds in this series, they retain or even show increased affinity towards the important Ile84Val mutation. By detailed analysis of the crystal structures of two representatives in complex with wild-type and mutant proteases, we were able to elucidate the structural basis of this phenomenon.     

Novel arylsulfonamides possessing sub-picomolar HIV protease activities and potent anti-HIV activity against wild-type and drug-resistant viral strains

A novel series of P1' chain-extended arylsufonamides was synthesized and evaluated for wild-type HIV protease inhibitory activity and in vitro antiviral activity against wild type virus and two protease inhibitor-resistant mutant viruses. All of the compounds showed dramatic increases in enzyme activity as compared to the currently marketed HIV protease inhibitors amprenavir, indinavir, and nelfinavir. In addition, significant improvements in antiviral potencies against wild type and the two mutant viruses were also realized.     

Relation between flexibility and positively selected HIV‐1 protease mutants against inhibitors

The antiretroviral chemotherapy helps to reduce the mortality of HIVs infected patients. However, RNA dependant virus replication has a high mutation rate. Human immunodeficiency virus Type 1 protease plays an essential role in viral replication cycle. This protein is an important target for therapy with viral protein inhibitors. There are few works using normal mode analysis to investigate this problem from the structural changes viewpoint. The investigation of protein flexibility may be important for the study of processes associated with conformational changes and state transitions. The normal mode analysis allowed us to investigate structural changes in the protease (such as flexibility) in a straightforward way and try to associate these changes with the increase of fitness for each positively selected HIV‐1 mutant protease of patients treated with several protease inhibitors (saquinavir, indinavir, ritonavir, nelfinavir, lopinavir, fosamprenavir, atazanavir, darunavir, and tripanavir) in combination or separately. These positively selected mutations introduce significant flexibility in important regions such as the active site cavity and flaps. These mutations were also able to cause changes in accessible solvent area. This study showed that the majority of HIV‐1 protease mutants can be grouped into two main classes of protein flexibility behavior. We presented a new approach to study structural changes caused by positively selected mutations in a pathogen protein, for instance the HIV‐1 protease and their relationship with their resistance mechanism against known inhibitors. The method can be applied to any pharmaceutically relevant pathogen proteins and could be very useful to understand the effects of positively selected mutations in the context of structural changes. Proteins 2012; © 2012 Wiley Periodicals, Inc.     

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.     

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.     

The crystal structure of protease Sapp1p from Candida parapsilosis in complex with the HIV protease inhibitor ritonavir

Secreted aspartic proteases (Saps) are extracellular proteolytic enzymes that enhance the virulence of Candida pathogens. These enzymes therefore represent possible targets for therapeutic drug design. Saps are inhibited by nanomolar concentrations of the classical inhibitor of aspartic proteases pepstatin A and also by the inhibitors of the HIV protease, but with the K(i) of micromolar values or higher. To contribute to the discussion regarding whether HIV protease inhibitors can act against opportunistic mycoses by the inhibition of Saps, we determined the structure of Sapp1p from Candida parapsilosis in complex with ritonavir (RTV), a clinically used inhibitor of the HIV protease. The crystal structure refined at resolution 2.4 ? proved binding of RTV into the active site of Sapp1p and provided the structural information necessary to evaluate the stability and specificity of the protein-inhibitor interaction.     

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.     

Pharmacological eEF2K activation promotes cell death and inhibits cancer progression

Activation of the elongation factor 2 kinase (eEF2K) leads to the phosphorylation and inhibition of the elongation factor eEF2, reducing mRNA translation rates. Emerging evidence indicates that the regulation of factors involved in protein synthesis may be critical for controlling diverse biological processes including cancer progression. Here we show that inhibitors of the HIV aspartyl protease (HIV‐PIs), nelfinavir in particular, trigger a robust activation of eEF2K leading to the phosphorylation of eEF2. Beyond its anti‐viral effects, nelfinavir has antitumoral activity and promotes cell death. We show that nelfinavir‐resistant cells specifically evade eEF2 inhibition. Decreased cell viability induced by nelfinavir is impaired in cells lacking eEF2K. Moreover, nelfinavir‐mediated anti‐tumoral activity is severely compromised in eEF2K‐deficient engrafted tumors in vivo. Our findings imply that exacerbated activation of eEF2K is detrimental for tumor survival and describe a mechanism explaining the anti‐tumoral properties of HIV‐PIs.     

Kinetic and thermodynamic characterization of HIV-1 protease inhibitors

Interaction kinetic and thermodynamic analyses provide information beyond that obtained in general inhibition studies, and may contribute to the design of improved inhibitors and increased understanding of molecular interactions. Thus, a biosensor-based method was used to characterize the interactions between HIV-1 protease and seven inhibitors, revealing distinguishing kinetic and thermodynamic characteristics for the inhibitors. Lopinavir had fast association and the highest affinity of the tested compounds, and the interaction kinetics were less temperature-dependent as compared with the other inhibitors. Amprenavir, indinavir and ritonavir showed non-linear temperature dependencies of the kinetics. The free energy, enthalpy and entropy (DeltaG, DeltaH, DeltaS) were determined, and the energetics of complex association (DeltaG(on), DeltaH(on), DeltaS(on)) and dissociation (DeltaG(off), DeltaH(off), DeltaS(off)) were resolved. In general, the energetics for the studied inhibitors was in the same range, with the negative free energy change (DeltaG < 0) due primarily to increased entropy (DeltaS > 0). Thus, the driving force of the interaction was increased degrees of freedom in the system (entropy) rather than the formation of bonds between the enzyme and inhibitor (enthalpy). Although the DeltaG(on) and DeltaG(off) were in the same range for all inhibitors, the enthalpy and entropy terms contributed differently to association and dissociation, distinguishing these phases energetically. Dissociation was accompanied by positive enthalpy (DeltaH(off) > 0) and negative entropy (DeltaS(off) < 0) changes, whereas association for all inhibitors except lopinavir had positive entropy changes (DeltaS(on) > 0), demonstrating unique energetic characteristics for lopinavir. This study indicates that this type of data will be useful for the characterization of target-ligand interactions and the development of new inhibitors of HIV-1 protease.     

Artificial neural network method for predicting HIV protease cleavage sites in protein

Knowledge of the polyprotein cleavage sites by HIV protease will refine our understanding of its specificity, and the information thus acquired will be useful for designing specific and efficient HIV protease inhibitors. The search for inhibitors of HIV protease will be greatly expedited if one can find and accurate, robust, and rapid method for predicting the cleavage sites in proteins by HIV protease. In this paper, Kohonen’s self-organization model, which uses typical artificial neural networks, is applied to predict the cleavability of oligopeptides by proteases with multiple and extended specificity subsites. We selected HIV-1 protease as the subject of study. We chose 299 oligopeptides for the training set, and another 63 oligopeptides for the test set. Because of its high rate of correct prediction (58/63=92.06%) and stronger fault-tolerant ability, the neural network method should be a useful technique for finding effective inhibitors of HIV protease, which is one of the targets in designing potential drugs against AIDS. The principle of the artificial neural network method can also be applied to analyzing the specificity of any multisubsite enzyme.     

Artificial neural network method for predicting HIV protease cleavage sites in protein

Knowledge of the polyprotein cleavage sites by HIV protease will refine our understanding of its specificity, and the information thus acquired will be useful for designing specific and efficient HIV protease inhibitors. The search for inhibitors of HIV protease will be greatly expedited if one can find and accurate, robust, and rapid method for predicting the cleavage sites in proteins by HIV protease. In this paper, Kohonen’s self-organization model, which uses typical artificial neural networks, is applied to predict the cleavability of oligopeptides by proteases with multiple and extended specificity subsites. We selected HIV-1 protease as the subject of study. We chose 299 oligopeptides for the training set, and another 63 oligopeptides for the test set. Because of its high rate of correct prediction (58/63=92.06%) and stronger fault-tolerant ability, the neural network method should be a useful technique for finding effective inhibitors of HIV protease, which is one of the targets in designing potential drugs against AIDS. The principle of the artificial neural network method can also be applied to analyzing the specificity of any multisubsite enzyme.     

Structural insights into the South African HIV-1 subtype C protease: impact of hinge region dynamics and flap flexibility in drug resistance

The HIV protease plays a major role in the life cycle of the virus and has long been a target in antiviral therapy. Resistance of HIV protease to protease inhibitors (PIs) is problematic for the effective treatment of HIV infection. The South African HIV-1 subtype C protease (C-SA PR), which contains eight polymorphisms relative to the consensus HIV-1 subtype B protease, was expressed in Escherichia coli, purified, and crystallized. The crystal structure of the C-SA PR was resolved at 2.7?Å, which is the first crystal structure of a HIV-1 subtype C protease that predominates in Africa. Structural analyses of the C-SA PR in comparison to HIV-1 subtype B proteases indicated that polymorphisms at position 36 of the homodimeric HIV-1 protease may impact on the stability of the hinge region of the protease, and hence the dynamics of the flap region. Molecular dynamics simulations showed that the flap region of the C-SA PR displays a wider range of movements over time as compared to the subtype B proteases. Reduced stability in the hinge region resulting from the absent E35-R57 salt bridge in the C-SA PR, most likely contributes to the increased flexibility of the flaps which may be associated with reduced susceptibility to PIs.

An animated interactive 3D complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:36