A folding inhibitor of the HIV-1 protease

Because the human immunodeficiency virus type 1 protease (HIV-1-PR) is an essential enzyme in the viral life cycle, its inhibition can control AIDS. The folding of single-domain proteins, like each of the monomers forming the HIV-1-PR homodimer, is controlled by local elementary structures (LES, folding units stabilized by strongly interacting, highly conserved, as a rule hydrophobic, amino acids). These LES have evolved over myriad generations to recognize and strongly attract each other, so as to make the protein fold fast and be stable in its native conformation. Consequently, peptides displaying a sequence identical to those segments of the monomers associated with LES are expected to act as competitive inhibitors and thus destabilize the native structure of the enzyme. These inhibitors are unlikely to lead to escape mutants as they bind to the protease monomers through highly conserved amino acids, which play an essential role in the folding process. The properties of one of the most promising inhibitors of the folding of the HIV-1-PR monomers found among these peptides are demonstrated with the help of spectrophotometric assays and circular dichroism spectroscopy.     

Suppression of HIV-1 protease inhibitor resistance by phosphonate-mediated solvent anchoring

The introduction of human immunodeficiency virus type 1 (HIV-1) protease inhibitors (PIs) markedly improved the clinical outcome and control of HIV-1 infection. However, cross-resistance among PIs due to a wide spectrum of mutations in viral protease is a major factor limiting their broader clinical use. Here we report on the suppression of PI resistance using a covalent attachment of a phosphonic acid motif to a peptidomimetic inhibitor scaffold. The resulting phosphonate analogs maintain high binding affinity to HIV-1 protease, potent antiretroviral activity, and unlike the parent molecules, display no loss of potency against a panel of clinically important PI-resistant HIV-1 strains. As shown by crystallographic analysis, the phosphonate moiety is highly exposed to solvent with no discernable interactions with any of the enzyme active site or surface residues. We term this effect "solvent anchoring" and demonstrate that it is driven by a favorable change in the inhibitor binding entropy upon the interaction with mutant enzymes. This type of thermodynamic behavior, which was not found with the parent scaffold fully buried in the enzyme active site, is a result of the increased degeneracy of inhibitor binding states, allowing effective molecular adaptation to the expanded cavity volume of mutant proteases. This strategy, which is applicable to various PI scaffolds, should facilitate the design of novel PIs and potentially other antiviral therapeutics.     

Molecular mechanisms of resistance: free energy calculations of mutation effects on inhibitor binding to HIV-1 protease

The changes in the inhibitor binding constants due to the mutation of isoleucine to valine at position 84 of HIV-1 protease are calculated using molecular dynamics simulations. The calculations are done for three potent inhibitors--KNI-272, L-735,524 (indinavir or MK-639), and Ro 31-8959 (saquinavir). The calculations agree with the experimental data both in terms of an overall trend and in the magnitude of the resulting free energy change. HIV-1 protease is a homodimer, so each mutation causes two changes in the enzyme. The decrease in the binding free energy from each mutated side chain differs among the three inhibitors and correlates well with the size of the cavities induced in the protein interior near the mutated residue. The cavities are created as a result of a mutation to a smaller side chain, but the cavities are less than would be predicted from the wild-type structures, indicating that there is significant relaxation to partially fill the cavities.     

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.     

Molecular mechanics analysis of inhibitor binding to HIV-1 protease.

Crystallographic structures of HIV protease with three different peptide-mimetic inhibitors were subjected to energy minimization using molecular mechanics, the minimized structures analyzed and the inhibitor binding energies calculated. Partial charge assignment for the hydrogen bonded catalytic aspartic acids, Asp25 and -25', was in good agreement with charge calculations using semi-empirical molecular orbital methods. Root mean square deviations on minimization were small and similar for both subunits in the protease dimer. The surface loops, which had the largest B factors, changed most on minimization; the hydrophobic core and the inhibitor binding site showed little change. The distance-dependent dielectric of D(r) = 4r was found to be preferable to D(r) = r. Distance restraints were applied for the intermolecular hydrogen bonds to maintain the conformation of the inhibitor binding site. Using the dielectric of D(r) = 4r, the calculated interaction energy of the three inhibitors with the protease ranged from -53 to -56 kcal/mol. The psi groups of the inhibitors were changed to add or remove a 'transition state analogue' hydroxyl group, and the loss in energy on the removal of this group was calculated to be 0.9-1.7 kcal/mol. This would represent 19-36% of the total measured difference in binding energy between the inhibitors JG365 and MVT-101.     

Evaluating the potency of HIV-1 protease drugs to combat resistance

HIV-1 protease has been an important drug target for the antiretroviral treatment of HIV infection. The efficacy of protease drugs is impaired by the rapid emergence of resistant virus strains. Understanding the molecular basis and evaluating the potency of an inhibitor to combat resistance are no doubt important in AIDS therapy. In this study, we first identified residues that have significant contributions to binding with six substrates using molecular dynamics simulations and Molecular Mechanics Generalized Born Surface Area calculations. Among the critical residues, Asp25, Gly27, Ala28, Asp29, and Gly49 are well conserved, with which the potent drugs should form strong interactions. We then calculated the contribution of each residue to binding with eight FDA approved drugs. We analyzed the conservation of each protease residue and also compared the interaction between the HIV protease and individual residues of the drugs and substrates. Our analyses showed that resistant mutations usually occur at less conserved residues forming more favorable interactions with drugs than with substrates. To quantitatively integrate the binding free energy and conservation information, we defined an empirical parameter called free energy/variability (FV) value, which is the product of the contribution of a single residue to the binding free energy and the sequence variability at that position. As a validation, the FV value was shown to identify single resistant mutations with an accuracy of 88%. Finally, we evaluated the potency of a newly approved drug, darunavir, to combat resistance and predicted that darunavir is more potent than amprenavir but may be susceptible to mutations on Val32 and Ile84.     

Combining mutations in HIV-1 protease to understand mechanisms of resistance

HIV-1 develops resistance to protease inhibitors predominantly by selecting mutations in the protease gene. Studies of resistant mutants of HIV-1 protease with single amino acid substitutions have shown a range of independent effects on specificity, inhibition, and stability. Four double mutants, K45I/L90M, K45I/V82S, D30N/V82S, and N88D/L90M were selected for analysis on the basis of observations of increased or decreased stability or enzymatic activity for the respective single mutants. The double mutants were assayed for catalysis, inhibition, and stability. Crystal structures were analyzed for the double mutants at resolutions of 2.2-1.2 A to determine the associated molecular changes. Sequence-dependent changes in protease-inhibitor interactions were observed in the crystal structures. Mutations D30N, K45I, and V82S showed altered interactions with inhibitor residues at P2/P2', P3/P3'/P4/P4', and P1/P1', respectively. One of the conformations of Met90 in K45I/L90M has an unfavorably close contact with the carbonyl oxygen of Asp25, as observed previously in the L90M single mutant. The observed catalytic efficiency and inhibition for the double mutants depended on the specific substrate or inhibitor. In particular, large variation in cleavage of p6(pol)-PR substrate was observed, which is likely to result in defects in the maturation of the protease from the Gag-Pol precursor and hence viral replication. Three of the double mutants showed values for stability that were intermediate between the values observed for the respective single mutants. D30N/V82S mutant showed lower stability than either of the two individual mutations, which is possibly due to concerted changes in the central P2-P2' and S2-S2' sites. The complex effects of combining mutations are discussed.     

The effect of inhibitor binding on the structural stability and cooperativity of the HIV-1 protease.

The effects of the peptide inhibitor acetyl pepstatin on the structural stability of the HIV-1 protease have been measured by high sensitivity calorimetric techniques. At 25 degrees C and pH 3.6, acetyl pepstatin binds to HIV-1 protease with an affinity of 1.6 x 10(7 )M-1 and an enthalpy of 7.3 +/- 0.5 kcal/mol, indicating that binding is not favored enthalpically and that the favorable Gibbs energy originates from a large positive entropy. Since the binding of acetyl pepstatin is associated with a negative change in heat capacity (-450 cal/K*mol) the association reaction becomes enthalpically favored at temperatures higher than 40 degrees C. The presence of the inhibitor stabilizes the dimeric structure of the protease in a fashion that can be quantitatively described by a set of thermodynamic linkage equations. The combination of titration and differential scanning calorimetry provides an accurate way of determining binding constants for high affinity inhibitors that cannot be determined by titration calorimetry alone. A structure-based thermodynamic analysis of the binding process indicates that the stabilization effect is not distributed uniformly throughout the protease molecule. The binding of the inhibitor selectively stabilizes those conformational states in which the binding site is formed, triggering a redistribution of the state probabilities in the ensemble of conformations populated under native conditions. As a result, the stability constants for individual residues do not exhibit the same change in magnitude upon inhibitor binding. Residues in certain areas of the protein are affected significantly whereas residues in other areas are not affected at all. In particular, inhibitor binding has a significant effect on those regions that define the binding site, especially the flap region which becomes structurally stable as a result of the additional binding free energy. The induced stabilization propagates to regions not in direct contact with the inhibitor, particularly to the strand between residues Pro9 and Ala22 and the helix between Arg87 and Gly94. On the other hand, the stability of the strand between Asp60 and Leu76 is not significantly affected by inhibitor binding. The structural distribution of binding effects define cooperative pathways within the protease molecule. Proteins 1999;36:147-156.     

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

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

Comparison of protease genotype and phenotype in HIV-1 infected patients exposed to more than one protease inhibitor

Many patterns of mutations selected by HIV-1 protease inhibitors have been described, but in most cases isolates with these patterns have been obtained from pre-clinical studies or after failures of monotherapies. We compared genotype and phenotype in HIV-1 infected patients who have failed more than one PI-including regimen. Phenotypic resistance could arise also in the absence of specific primary mutations and in the presence of different substitutions among those known to confer resistance to ritonavir, indinavir or nelfinavir. The number of secondary mutations was significantly associated with phenotypic resistance for each protease inhibitor. Thus, more study of mutational patterns in heavily pretreated patients is warranted; in the mean time treatment choices might be optimized if phenotyping could integrate genotyping within this setting.     

A new structural theme in C2-symmetric HIV-1 protease inhibitors: ortho-substituted P1/P1' side chains

In this report, the rapid syntheses of 24 novel C2-symmetric HIV-1 protease inhibitors are described. Two ortho-iodobenzyloxy containing C-terminal duplicated inhibitors served as starting materials for microwave-enhanced palladium(0)-catalyzed carbon-carbon bond forming reactions (Suzuki, Sonogashira, Heck, and Negishi). Highly potent inhibitors equipped with ortho-functionalized P1/P1' side chains as the structural theme were identified. Computational efforts were applied to study the binding mode of this class of inhibitors and to establish structure-activity relationships. The overall orientation of the inhibitors in the active site was reproduced by docking which suggested three possible conformations of the P1/P1' groups of which two seem more plausible.     

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

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

Comparative sensitivity of Vibrio cholerae 01 E1 Tor and Escherichia coli to disinfectants

The sensitivity of two strains of Vibrio cholerae to disinfectant compounds used in food processing, kitchen and personal hygiene has been compared with the sensitivity of a 'disinfectant-test' strain of Escherichia coli. In a suspension test, both strains of V. cholerae were slightly more sensitive than E. coli to all the compounds. When used in a hard-surface disinfection test, the vibrios died rapidly during the initial drying phase. Disinfectant products which are effective to control the risks from pathogenic enterobacteriaceae should also be appropriate for V. cholerae.     

A target within the target: probing cruzain's P1' site to define structural determinants for the Chagas' disease protease

BACKGROUND: Cysteine proteases of the papain superfamily are present in nearly all groups of eukaryotes and play vital roles in a wide range of biological processes and diseases, including antigen and hormone processing, bacterial infection, arthritis, osteoporosis, Alzheimer's disease and cancer-cell invasion. Because they are critical to the life-cycle progression of many pathogenic protozoa, they represent potential targets for selective inhibitors. Chagas' disease, the leading cause of death due to heart disease in Latin American countries, is transmitted by Trypanosoma cruzi. Cruzain is the major cysteine protease of T cruzi and has been the target of extensive structure-based drug design. RESULTS: High-resolution crystal structures of cruzain bound to a series of potent phenyl-containing vinyl-sulfone, sulfonate and sulfonamide inhibitors have been determined. The structures show a consistent mode of interaction for this family of inhibitors based on a covalent Michael addition formed at the enzyme's active-site cysteine, hydrophobic interactions in the S2 substrate-binding pocket and a strong constellation of hydrogen bonding in the S1' region. CONCLUSIONS: The series of vinyl-sulfone-based inhibitors examined in complex with cruzain was designed to probe recognition and binding potential of an aromatic-rich region of the enzyme. Analysis of the interactions formed shows that aromatic interactions play a less significant role, whereas the strength and importance of hydrogen bonding in the conformation adopted by the inhibitor upon binding to the enzyme was highlighted. A derivative of one inhibitor examined is currently under development as a therapeutic agent against Chagas' disease.     

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

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