Mutational and Structural Studies Aimed at Characterizing the Monomer of HIV-1 Protease and Its Precursor

Interplay between protease and reverse transcriptase dimerization in a model HIV-1 polyprotein

The Gag-Pol polyprotein in human immunodeficiency virus type I (HIV-1) encodes enzymes that are essential for virus replication: protease (PR), reverse transcriptase (RT), and integrase (IN). The mature forms of PR, RT and IN are homodimer, heterodimer and tetramer, respectively. The precise mechanism underlying the formation of dimer or tetramer is not yet understood. Here, to gain insight into the dimerization of PR and RT in the precursor, we prepared a model precursor, PR-RT, incorporating an inactivating mutation at the PR active site, D25A, and including two residues in the p6* region, fused to a SUMO-tag, at the N-terminus of the PR region. We also prepared two mutants of PR-RT containing a dimer dissociation mutation either in the PR region, PR(T26A)-RT, or in the RT region, PR-RT(W401A). Size exclusion chromatography showed both monomer and dimer fractions in PR-RT and PR(T26A)-RT, but only monomer in PR-RT(W401A). SEC experiments of PR-RT in the presence of protease inhibitor, darunavir, significantly enhanced the dimerization. Additionally, SEC results suggest an estimated PR-RT dimer dissociation constant that is higher than that of the mature RT heterodimer, p66/p51, but slightly lower than the premature RT homodimer, p66/p66. Reverse transcriptase assays and RT maturation assays were performed as tools to assess the effects of the PR dimer-interface on these functions. Our results consistently indicate that the RT dimer-interface plays a crucial role in the dimerization in PR-RT, whereas the PR dimer-interface has a lesser role.

Rational design, synthesis and biological evaluation of novel HIV-1 protease inhibitors containing 2-phenylacetamide derivatives as P2 ligands with potent activity against DRV-Resistant HIV-1 variants

A novel kind of potent HIV-1 protease inhibitors, containing diverse hydroxyphenylacetic acids as the P2-ligands and 4-substituted phenyl sulfonamides as the P2' ligands, were designed, synthesized and evaluated in this work. Majority of the target compounds exhibited good to excellent activity against HIV-1 protease with IC50 values below 200 nM. In particular, compound 18d with a 2-(3,4-dihydroxyphenyl) acetamide as the P2 ligand and a 4- methoxybenzene sulfonamide P2' ligand exhibited inhibitory activity IC50 value of 0.54 nM, which was better than that of the positive control darunavir (DRV). More importantly, no significant decline of the potency against HIV-1DRVRS (DRV-resistant mutation) and HIV-1NL4_3 variant (wild type) for 18d was detected. The molecular docking study of 18d with HIV-1 protease (PDB-ID: 1T3R, www.rcsb.org) revealed possible binding mode with the HIV-1 protease. These results suggested the validity of introducing phenol-derived moieties into the P2 ligand and deserve further optimization which was of great value for future discovery of novel HIV-1 protease.

Mechanism of darunavir binding to monomeric HIV-1 protease: a step forward in the rational design of dimerization inhibitors

HIV protease (HIVPR) is a key target in AIDS therapeutics. All ten FDA-approved drugs that compete with substrates in binding to this dimeric enzyme's active site have become ineffective due to the emergence of drug resistant mutants. Blocking the dimerization interface of HIVPR is thus being explored as an alternate strategy. The latest drug, darunavir (DRV), which exhibited a high genetic barrier to viral resistance, is said to have a dual mode of action - (i) binding to the dimeric active site, and (ii) preventing the dimerization by binding to the HIVPR monomer. Despite several reports on DRV complexation with dimeric HIVPR, the mode and mechanism of the binding of DRV to the HIVPR monomer are poorly understood. In this study, we utilized all-atomic MD simulations and umbrella sampling techniques to identify the best possible binding mode of DRV to the monomeric HIVPR and its mechanism of association. The results suggest that DRV binds between the active site and the flap of the monomer, and the flap plays a crucial role in directing the drug to bind and driving the other protein domains to undergo induced fit changes for stronger complexation. The obtained binding mode of DRV was validated by comparing with various mutational data from clinical isolates to reported in vitro mutations. The identified binding pose was also able to successfully reproduce the experimental K_i value in the picomolar range. The residue-level information extracted from this study could accelerate the structure-based drug designing approaches targeting HIVPR dimerization.

Amino acid substitutions at the HIV-1 transframe region significantly impair virus infectivity

A transframe region within HIV-1 Gag-Pol (referred to as p6* or p6pol), directly linked to the protease (PR) N-terminus, plays a pivotal role in modulating PR activation. To identify specific p6* residues involved in PR activation, we created a series of p6* mutants by making substitutions for conserved p6* residues. Our results indicate that some p6* mutants were defective in terms of virus infectivity, despite displaying a wild-type virus particle processing pattern. Mutations at p6* F8 reduced virus infectivity associated with insufficient virus processing, due in part to impaired PR maturation and RT packaging. Our data strongly suggest that conserved Phe (F) residues at position 8 of p6* are involved in the PR maturation process.

Conditional activation of an HIV-1 protease attenuated mutant by a leucine zipper dimerization motif

Mature HIV-1 protease (PR) functions as a dimer. Changes in HIV-1 PR activation can block virus assembly via premature or enhanced Gag cleavage. HIV-1 PR precursor contains N terminal-linked p6*, a possible modulating factor in PR activation. We found that p6* replacement with a leucine zipper (LZ) dimerization motif (creating a DWzPR construct) or an LZ insertion at the PR C-terminus significantly reduced virus yields due to enhanced Gag cleavage, suggesting that an LZ insertion promotes PR activation by facilitating PR dimer formation. However, introducing T26S (a PR activity-attenuated mutation) into DWzPR strongly impaired Gag cleavage, except when the native C-terminal p6* tetrapeptide remained at the LZ/PR junction. LZ insertion at the PR C-terminus still strongly enhanced PR T26S Gag cleavage. Our data suggest that in addition to p6* mutations, a single amino acid substitution within PR can impair PR activation, likely due to conformational changes triggered by the PR precursor.

Effect of Lysyl-tRNA Synthetase on the Maturation of HIV-1 Reverse Transcriptase

In human immunodeficiency virus-1 (HIV-1), reverse transcriptase (RT) is encoded as a 66 kDa protein, p66, in the Gag-Pol polyprotein. This protein is proteolytically cleaved by HIV-1 protease (PR) to finally generate a mature RT that is a heterodimer, composed of a p66 subunit and a p66-derived 51 kDa subunit, p51. In our prior work, we demonstrated that tRNA^Lys3 binding to p66/p66 facilitates efficient cleavage of p66 to p51 by PR. However, tRNA^Lys3 is known to be recruited to the virus by forming a complex with lysyl-tRNA synthetase (LysRS). Herein, we tested whether LysRS can have an effect on RT maturation in vitro. Importantly, our data show no significant differences in RT maturation in the presence of LysRS. Furthermore, no apparent p66/66 interaction with LysRS was observed. Although PR cleaved LysRS, it did not immediately release tRNA^Lys3 from LysRS. Thus, we conclude that a free fraction of tRNA^Lys3, which is in equilibrium with a LysRS-bound form, interacts with p66/p66 without any additional mechanism involving release of tRNA^Lys3 from LysRS. Given that only transient tRNA^Lys3-p66/p66 interaction is needed for efficient RT maturation, a small amount of free tRNA may be sufficient for this process. These studies reveal molecular level insights into RT maturation and will be useful for the design of cellular/viral experiments to better understand the role of tRNA in HIV-1 replication.

Diverse Folding Pathways of HIV-1 Protease Monomer on a Rugged Energy Landscape

The modern energy landscape theory of protein folding predicts multiple folding pathways connecting a myriad of unfolded conformations and a well-defined folded state. However, direct experimental observation of heterogeneous folding pathways is difficult. Naturally evolved proteins typically exhibit a smooth folding energy landscape for fast and efficient folding by avoiding unfavorable kinetic traps. In this case, rapid fluctuations between unfolded conformations result in apparent two-state behavior and make different pathways indistinguishable. However, the landscape roughness can be different, depending on the selection pressures during evolution. Here, we characterize the unusually rugged folding energy landscape of human immunodeficiency virus-1 protease monomer using single-molecule Förster resonance energy transfer spectroscopy. Our data show that fluctuations between unfolded conformations are slow, which enables the experimental observation of heterogeneous folding pathways as predicted by the landscape theory. Although the landscape ruggedness is sensitive to the mutations and fluorophore locations, the folding rate is similar for various protease constructs. The natural evolution of the protease to have a rugged energy landscape likely results from intrinsic pressures to maintain robust folding when human immunodeficiency virus-1 mutates frequently, which is essential for its survival.

Inhibition of the precursor and mature forms of HIV-1 protease as a tool for drug evaluation

HIV-1 protease (PR) is a homodimeric enzyme that is autocatalytically cleaved from the Gag-Pol precursor. Known PR inhibitors bind the mature enzyme several orders of magnitude more strongly than the PR precursor. Inhibition of PR at the precursor level, however, may stop the process at its rate-limiting step before the proteolytic cascade is initiated. Due to its structural heterogeneity, limited solubility and autoprocessing, the PR precursor is difficult to access by classical methods, and limited knowledge regarding precursor inhibition is available. Here, we describe a cell-based assay addressing precursor inhibition. We used a reporter molecule containing the transframe (TFP) and p6* peptides, PR, and N-terminal fragment of reverse transcriptase flanked by the fluorescent proteins mCherry and EGFP on its N- and C- termini, respectively. The level of FRET between EGFP and mCherry indicates the amount of unprocessed reporter, allowing specific monitoring of precursor inhibition. The inhibition can be quantified by flow cytometry. Additionally, two microscopy techniques confirmed that the reporter remains unprocessed within individual cells upon inhibition. We tested darunavir, atazanavir and nelfinavir and their combinations against wild-type PR. Shedding light on an inhibitor’s ability to act on non-mature forms of PR may aid novel strategies for next-generation drug design.

HIV-1 protease with leucine zipper fused at N-terminus exhibits enhanced linker amino acid-dependent activity

Background

HIV-1 protease (PR) activation is triggered by Gag-Pol dimerization. Premature PR activation results in reduced virion yields due to enhanced Gag cleavage. A p6* transframe peptide located directly upstream of protease is believed to play a modulating role in PR activation. Previous reports indicate that the C-terminal p6* tetra-peptide prevents premature PR activation triggered by a leucine zipper (LZ) dimerization motif inserted in the deleted p6* region. To clarify the involvement of C-terminal p6* residues in mitigating enhanced LZ-incurred Gag processing, we engineered constructs containing C-terminal p6* residue substitutions with and without a mutation blocking the p6*/PR cleavage site, and created other Gag or p6* domain-removing constructs. The capabilities of these constructs to mediate virus maturation were assessed by Western blotting and single-cycle infection assays.

Results

p6*-PR cleavage blocking did not significantly reduce the LZ enhancement effect on Gag cleavage when only four amino acid residues were present between the p6* and PR. This suggests that the potent LZ dimerization motif may enhance PR activation by facilitating PR dimer formation, and that PR precursors may trigger sufficient enzymatic activity without breaking off from the PR N-terminus. Enhanced LZ-induced activation of PR embedded in Gag-Pol was found to be independent of the Gag assembly domain. In contrast, the LZ enhancement effect was markedly reduced when six amino acids were present at the p6*-PR junction, in part due to impaired PR maturation by substitution mutations. We also observed that a proline substitution at the P3 position eliminated the ability of p6*-deleted Gag-Pol to mediate virus maturation, thus emphasizing the importance of C-terminal p6* residues to modulating PR activation.

Conclusions

The ability of HIV-1 C-terminal p6* amino acid residues to modulate PR activation contributes, at least in part, to their ability to counteract enhanced Gag cleavage induced by a leucine zipper substituted for a deleted p6*. Changes in C-terminal p6* residues between LZ and PR may affect PR-mediated virus maturation, thus providing a possible method for assessing HIV-1 protease precursor activation in the context of virus assembly.

Probing Structural Changes among Analogous Inhibitor-Bound Forms of HIV-1 Protease and a Drug-Resistant Mutant in Solution by Nuclear Magnetic Resonance

In the era of state-of-the-art inhibitor design and high-resolution structural studies, detection of significant but small protein structural differences in the inhibitor-bound forms is critical to further developing the inhibitor. Here, we probed differences in HIV-1 protease (PR) conformation among darunavir and four analogous inhibitor-bound forms and compared them with a drug-resistant mutant using nuclear magnetic resonance chemical shifts. Changes in amide chemical shifts of wild-type (WT) PR among these inhibitor-bound forms, ΔCSP, were subtle but detectable and extended >10 Å from the inhibitor-binding site, asymmetrically between the two subunits of PR. Molecular dynamics simulations revealed differential local hydrogen bonding as the molecular basis of this remote asymmetric change. Inhibitor-bound forms of the drug-resistant mutant also showed a similar long-range ΔCSP pattern. Differences in ΔCSP values of the WT and the mutant (ΔΔCSPs) were observed at the inhibitor-binding site and in the surrounding region. Comparing chemical shift changes among highly analogous inhibitors and ΔΔCSPs effectively eliminated local environmental effects stemming from different chemical groups and enabled exploitation of these sensitive parameters to detect subtle protein conformational changes and to elucidate asymmetric and remote conformational effects upon inhibitor interaction.

C-Terminal HIV-1 Transframe p6* Tetrapeptide Blocks Enhanced Gag Cleavage Incurred by Leucine Zipper Replacement of a Deleted p6* Domain

HIV-1 protease (PR) functions as a homodimer mediating virus maturation following virus budding. Gag-Pol dimerization is believed to trigger embedded PR activation by promoting PR dimer formation. Early PR activation can lead to markedly reduced virus yields due to premature Gag cleavage. The p6* peptide, located between Gag and PR, is believed to ensure virus production by preventing early PR maturation. Studies aimed at finding supporting evidence for this proposal are limited due to a reading frame overlap between p6* and the p6gag budding domain. To determine if p6* affects virus production via the modulation of PR activation, we engineered multiple constructs derived from Dp6*PR (an assembly- and processing-competent construct with Pol fused at the inactivated PR C terminus). The data indicated that a p6* deletion adjacent to active PR significantly impaired virus processing. We also observed that the insertion of a leucine zipper (LZ) dimerization motif in the deleted region eliminated virus production in a PR activity-dependent manner, suggesting that the LZ insertion triggered premature PR activation by facilitating PR dimer formation. As few as four C-terminal p6* residues remaining at the p6*/PR junction were sufficient to restore virus yields, with a Gag processing profile similar to that of the wild type. Our study provides supporting evidence in a virus assembly context that the C-terminal p6* tetrapeptide plays a role in preventing premature PR maturation.IMPORTANCE Supporting evidence for the assumption that p6* retards PR maturation in the context of virus assembly is lacking. We found that replacing p6* with a leucine zipper peptide abolished virus assembly due to the significant enhancement of Gag cleavage. However, as few as four C-terminal p6* residues remaining in the deleted region were sufficient for significant PR release, as well as for counteracting leucine zipper-incurred premature Gag cleavage. Our data provide evidence that (i) p6* ensures virus assembly by preventing early PR activation and (ii) four C-terminal p6* residues are critical for modulating PR activation. Current PR inhibitor development efforts are aimed largely at mature PR, but there is a tendency for HIV-1 variants that are resistant to multiple protease inhibitors to emerge. Our data support the idea of modulating PR activation by targeting PR precursors as an alternative approach to controlling HIV-1/AIDS.

Structural Studies of a Rationally Selected Multi-Drug Resistant HIV-1 Protease Reveal Synergistic Effect of Distal Mutations on Flap Dynamics

We report structural analysis of HIV protease variant PR^S17 which was rationally selected by machine learning to represent wide classes of highly drug-resistant variants. Crystal structures were solved of PR^S17 in the inhibitor-free form and in complex with antiviral inhibitor, darunavir. Despite its 17 mutations, PR^S17 has only one mutation (V82S) in the inhibitor/substrate binding cavity, yet exhibits high resistance to all clinical inhibitors. PR^S17 has none of the major mutations (I47V, I50V, I54ML, L76V and I84V) associated with darunavir resistance, but has 10,000-fold weaker binding affinity relative to the wild type PR. Comparable binding affinity of 8000-fold weaker than PR is seen for drug resistant mutant PR20, which bears 3 mutations associated with major resistance to darunavir (I47V, I54L and I84V). Inhibitor-free PR^S17 shows an open flap conformation with a curled tip correlating with G48V flap mutation. NMR studies on inactive PR^S17_D25N unambiguously confirm that the flaps adopt mainly an open conformation in solution very similar to that in the inhibitor-free crystal structure. In PR^S17, the hinge loop cluster of mutations, E35D, M36I and S37D, contributes to the altered flap dynamics by a mechanism similar to that of PR20. An additional K20R mutation anchors an altered conformation of the hinge loop. Flap mutations M46L and G48V in PR^S17/DRV complex alter the Phe53 conformation by steric hindrance between the side chains. Unlike the L10F mutation in PR20, L10I in PR^S17 does not break the inter-subunit ion pair or diminish the dimer stability, consistent with a very low dimer dissociation constant comparable to that of wild type PR. Distal mutations A71V, L90M and I93L propagate alterations to the catalytic site of PR^S17. PR^S17 exhibits a molecular mechanism whereby mutations act synergistically to alter the flap dynamics resulting in significantly weaker binding yet maintaining active site contacts with darunavir.

Evolution under Drug Pressure Remodels the Folding Free-Energy Landscape of Mature HIV-1 Protease

Using high-pressure NMR spectroscopy and differential scanning calorimetry, we investigate the folding landscape of the mature HIV-1 protease homodimer. The cooperativity of unfolding was measured in the absence or presence of a symmetric active site inhibitor for the optimized wild type protease (PR), its inactive variant PRD25N, and an extremely multidrug-resistant mutant, PR20. The individual fit of the pressure denaturation profiles gives rise to first order, ∆GNMR, and second order, ∆VNMR (the derivative of ∆GNMR with pressure); apparent thermodynamic parameters for each amide proton considered. Heterogeneity in the apparent ∆VNMR values reflects departure from an ideal cooperative unfolding transition. The narrow to broad distribution of ∆VNMR spanning the extremes from inhibitor-free PR20D25N to PR-DMP323 complex, and distinctively for PRD25N-DMP323 complex, indicated large variations in folding cooperativity. Consistent with this data, the shape of thermal unfolding transitions varies from asymmetric for PR to nearly symmetric for PR20, as dimer-inhibitor ternary complexes. Lack of structural cooperativity was observed between regions located close to the active site, including the hinge and tip of the glycine-rich flaps, and the rest of the protein. These results strongly suggest that inhibitor binding drastically decreases the cooperativity of unfolding by trapping the closed flap conformation in a deep energy minimum. To evade this conformational trap, PR20 evolves exhibiting a smoother folding landscape with nearly an ideal two-state (cooperative) unfolding transition. This study highlights the malleability of retroviral protease folding pathways by illustrating how the selection of mutations under drug pressure remodels the free-energy landscape as a primary mechanism.

Binding of Clinical Inhibitors to a Model Precursor of a Rationally Selected Multidrug Resistant HIV-1 Protease Is Significantly Weaker Than That to the Released Mature Enzyme

We have systematically validated the activity and inhibition of a HIV-1 protease (PR) variant bearing 17 mutations (PR(S17)), selected to represent high resistance by machine learning on genotype-phenotype data. Three of five mutations in PR(S17) correlating with major drug resistance, M46L, G48V, and V82S, and five of 11 natural variations differ from the mutations in two clinically derived extreme mutants, PR20 and PR22 bearing 19 and 22 mutations, respectively. PR(S17), which forms a stable dimer (<10 nM), is ∼10- and 2-fold less efficient in processing the Gag polyprotein than the wild type and PR20, respectively, but maintains the same cleavage order. Isolation of a model precursor of PR(S17) flanked by the 56-amino acid transframe region (TFP-p6pol) at its N-terminus, which is impossible upon expression of an analogous PR20 precursor, allowed systematic comparison of inhibition of TFP-p6pol-PR(S17) and mature PR(S17). Resistance of PR(S17) to eight protease inhibitors (PIs) relative to PR (Ki) increases by 1.5-5 orders of magnitude from 0.01 to 8.4 μM. Amprenavir, darunavir, atazanavir, and lopinavir, the most effective of the eight PIs, inhibit precursor autoprocessing at the p6pol/PR site with IC50 values ranging from ∼7.5 to 60 μM. Thus, this process, crucial for stable dimer formation, shows inhibition ∼200-800-fold weaker than that of the mature PR(S17). TFP/p6pol cleavage, which occurs faster, is inhibited even more weakly by all PIs except darunavir (IC50 = 15 μM); amprenavir shows a 2-fold increase in IC50 (∼15 μM), and atazanavir and lopinavir show increased IC50 values of >42 and >70 μM, respectively.

Mutations Proximal to Sites of Autoproteolysis and the α-Helix That Co-evolve under Drug Pressure Modulate the Autoprocessing and Vitality of HIV-1 Protease

N-Terminal self-cleavage (autoprocessing) of the HIV-1 protease precursor is crucial for liberating the active dimer. Under drug pressure, evolving mutations are predicted to modulate autoprocessing, and the reduced catalytic activity of the mature protease (PR) is likely compensated by enhanced conformational/dimer stability and reduced susceptibility to self-degradation (autoproteolysis). One such highly evolved, multidrug resistant protease, PR20, bears 19 mutations contiguous to sites of autoproteolysis in retroviral proteases, namely clusters 1-3 comprising residues 30-37, 60-67, and 88-95, respectively, accounting for 11 of the 19 mutations. By systematically replacing corresponding clusters in PR with those of PR20, and vice versa, we assess their influence on the properties mentioned above and observe no strict correlation. A 10-35-fold decrease in the cleavage efficiency of peptide substrates by PR20, relative to PR, is reflected by an only ∼4-fold decrease in the rate of Gag processing with no change in cleavage order. Importantly, optimal N-terminal autoprocessing requires all 19 PR20 mutations as evaluated in vitro using the model precursor TFR-PR20 in which PR is flanked by the transframe region. Substituting PR20 cluster 3 into TFR-PR (TFR-PR(PR20-3)) requires the presence of PR20 cluster 1 and/or 2 for autoprocessing. In accordance, substituting PR clusters 1 and 2 into TFR-PR20 affects the rate of autoprocessing more drastically (>300-fold) compared to that of TFR-PR(PR20-3) because of the cumulative effect of eight noncluster mutations present in TFR-PR20(PR-12). Overall, these studies imply that drug resistance involves a complex synchronized selection of mutations modulating all of the properties mentioned above governing PR regulation and function.

Gag-Pol Transframe Domain p6* Is Essential for HIV-1 Protease-Mediated Virus Maturation

HIV-1 protease (PR) is encoded by pol, which is initially translated as a Pr160^gag-pol polyprotein by a ribosomal frameshift event. Within Gag-Pol, truncated p6gag is replaced by a transframe domain (referred to as p6* or p6pol) located directly upstream of PR. p6* has been proposed as playing a role in modulating PR activation. Overlapping reading frames between p6* and p6gag present a challenge to researchers using genetic approaches to studying p6* biological functions. To determine the role of p6* in PR activation without affecting the gag reading frame, we constructed a series of Gag/Gag-Pol expression vectors by duplicating PR with or without p6* between PR pairs, and observed that PR duplication eliminated virus production due to significant Gag cleavage enhancement. This effect was mitigated when p6* was placed between the two PRs. Further, Gag cleavage enhancement was markedly reduced when either one of the two PRs was mutationally inactivated. Additional reduction in Gag cleavage efficiency was noted following the removal of p6* from between the two PRs. The insertion of a NC domain (wild-type or mutant) directly upstream of PR or p6*PR did not significantly improve Gag processing efficiency. With the exception of those containing p6* directly upstream of an active PR, all constructs were either noninfectious or weakly infectious. Our results suggest that (a) p6* is essential for triggering PR activation, (b) p6* has a role in preventing premature virus processing, and (c) the NC domain within Gag-Pol is not a major determinant of PR activation.

Inhibitor and substrate binding induced stability of HIV-1 protease against sequential dissociation and unfolding revealed by high pressure spectroscopy and kinetics

High-pressure methods have become an interesting tool of investigation of structural stability of proteins. They are used to study protein unfolding, but dissociation of oligomeric proteins can be addressed this way, too. HIV-1 protease, although an interesting object of biophysical experiments, has not been studied at high pressure yet. In this study HIV-1 protease is investigated by high pressure (up to 600 MPa) fluorescence spectroscopy of either the inherent tryptophan residues or external 8-anilino-1-naphtalenesulfonic acid at 25°C. A fast concentration-dependent structural transition is detected that corresponds to the dimer-monomer equilibrium. This transition is followed by a slow concentration independent transition that can be assigned to the monomer unfolding. In the presence of a tight-binding inhibitor none of these transitions are observed, which confirms the stabilizing effect of inhibitor. High-pressure enzyme kinetics (up to 350 MPa) also reveals the stabilizing effect of substrate. Unfolding of the protease can thus proceed only from the monomeric state after dimer dissociation and is unfavourable at atmospheric pressure. Dimer-destabilizing effect of high pressure is caused by negative volume change of dimer dissociation of -32.5 mL/mol. It helps us to determine the atmospheric pressure dimerization constant of 0.92 μM. High-pressure methods thus enable the investigation of structural phenomena that are difficult or impossible to measure at atmospheric pressure.

Conformation of inhibitor-free HIV-1 protease derived from NMR spectroscopy in a weakly oriented solution

Flexibility of the glycine-rich flaps is known to be essential for catalytic activity of the HIV-1 protease, but their exact conformations at the different stages of the enzymatic pathway remain subject to much debate. Although hundreds of crystal structures of protease-inhibitor complexes have been solved, only about a dozen inhibitor-free protease structures have been reported. These latter structures reveal a large diversity of flap conformations, ranging from closed to semi-open to wide open. To evaluate the average structure in solution, we measured residual dipolar couplings (RDCs) and compared these to values calculated for crystal structures representative of the closed, semi-open, and wide-open states. The RDC data clearly indicate that the inhibitor-free protease, on average, adopts a closed conformation in solution that is very similar to the inhibitor-bound state. By contrast, a highly drug-resistant protease mutant, PR20, adopts the wide-open flap conformation.

Dimerization of HIV-1 protease occurs through two steps relating to the mechanism of protease dimerization inhibition by darunavir

Dimerization of HIV-1 protease (PR) subunits is an essential process for PR's acquisition of proteolytic activity, which plays a critical role in the maturation of HIV-1. Recombinant wild-type PR (PR(WT)) proved to dimerize, as examined with electrospray ionization mass spectrometry; however, two active site interface PR mutants (PR(T26A) and PR(R87K)) remained monomeric. On the other hand, two termini interface PR mutants (PR(1-C95A) and PR(97/99)) took both monomeric and dimeric forms. Differential scanning fluorimetry indicated that PR(1-C95A) and PR(97/99) dimers were substantially less stable than PR(WT) dimers. These data indicate that intermolecular interactions of two monomers occur first at the active site interface, generating unstable or transient dimers, and interactions at the termini interface subsequently occur, generating stable dimers. Darunavir (DRV), an HIV-1 protease inhibitor, inhibits not only proteolytic activity but also PR dimerization. DRV bound to protease monomers in a one-to-one molar ratio, inhibiting the first step of PR dimerization, whereas conventional protease inhibitors (such as saquinavir) that inhibit enzymatic activity but not dimerization failed to bind to monomers. DRV also bound to mutant PRs containing the transframe region-added PR (TFR-PR(D25N) and TFR-PR(D25N-7AA)), whereas saquinavir did not bind to TFR-PR(D25N) or TFR-PR(D25N-7AA). Notably, DRV failed to bind to mutant PR containing four amino acid substitutions (V32I, L33F, I54M, and I84V) that confer resistance to DRV on HIV-1. To our knowledge, the present report represents the first demonstration of the two-step PR dimerization dynamics and the mechanism of dimerization inhibition by DRV, which should help design further, more potent novel PIs.

The role of select subtype polymorphisms on HIV-1 protease conformational sampling and dynamics

HIV-1 protease is an essential enzyme for viral particle maturation and is a target in the fight against HIV-1 infection worldwide. Several natural polymorphisms are also associated with drug resistance. Here, we utilized both pulsed electron double resonance, also called double electron-electron resonance, and NMR (15)N relaxation measurements to characterize equilibrium conformational sampling and backbone dynamics of an HIV-1 protease construct containing four specific natural polymorphisms commonly found in subtypes A, F, and CRF_01 A/E. Results show enhanced backbone dynamics, particularly in the flap region, and the persistence of a novel conformational ensemble that we hypothesize is an alternative flap orientation of a curled open state or an asymmetric configuration when interacting with inhibitors.

Investigation on the mechanism for the binding and drug resistance of wild type and mutations of G86 residue in HIV-1 protease complexed with Darunavir by molecular dynamic simulation and free energy calculation

Residue Gly86 is considered as the highly conversed residue in the HIV-1 protease. In our work, the detailed binding free energies for the wild-type (WT) and mutated proteases binding to the TMC-114 are estimated to investigate the protein-inhibitor binding and drug resistance mechanism by molecule dynamic simulations and molecular mechanics Poisson Boltzmann surface area (MM-PBSA) method. The binding affinities between the mutants and inhibitor are different than that in the wild-type complex and the major resistance to Darunavir (DRV) of G86A and G86S originate from the electrostatic energy and entropy, respectively. Furthermore, free energy decomposition analysis for the WT and mutated complexes on the basis of per-residue indicates that the mutagenesis influences the energy contribution of the residue located at three regions: active site region (residue 24-32), the flap region, and the region around the mutated residue G86 (residue 79-88), especially the flap region. Finally, further hydrogen bonds and structure analysis are carried out to detect the relationship between the energy and conformation. In all, the G86 mutations change the flap region's conformation. The experimental results are in good agreement with available results.

Understanding HIV-1 protease autoprocessing for novel therapeutic development

In the infected cell, HIV-1 protease (PR) is initially synthesized as part of the GagPol polyprotein. PR autoprocessing is a virus-specific process by which the PR domain embedded in the precursor catalyzes proteolytic reactions responsible for liberation of free mature PRs, which then recognize and cleave at least ten different peptide sequences in the Gag and GagPol polyproteins. Despite extensive structure and function studies of the mature PRs as well as the successful development of ten US FDA-approved catalytic-site inhibitors, the precursor autoprocessing mechanism remains an intriguing yet-to-be-solved puzzle. This article discusses current understanding of the autoprocessing mechanism, in an effort to prompt the development of novel anti-HIV drugs that selectively target precursor autoprocessing.

The maturation of HIV-1 protease precursor studied by discrete molecular dynamics

The equilibrium properties of a HIV-1-protease precursor are studied by means of an efficient molecular dynamics scheme, which allows for the simulation of the folding of the protein monomers and their dimerization into an active form and compare them with those of the mature protein. The results of the model provide, with atomic detail, an overall account of several experimental findings, including the NMR conformation of the mature dimer, the calorimetric properties of the system, the effects of the precursor tail on the dimerization constant, the secondary chemical shifts of the monomer, and the paramagnetic relaxation enhancement data associated with the conformations of the precursor. It is found that although the mature protein can dimerize in a unique, single way, the precursor populates several dimeric conformations in which monomers are always native-like, but their binding can be non-native.

Enhanced stability of monomer fold correlates with extreme drug resistance of HIV-1 protease

During treatment, mutations in HIV-1 protease (PR) are selected rapidly that confer resistance by decreasing affinity to clinical protease inhibitors (PIs). As these unique drug resistance mutations can compromise the fitness of the virus to replicate, mutations that restore conformational stability and activity while retaining drug resistance are selected on further evolution. Here we identify several compensating mechanisms by which an extreme drug-resistant mutant bearing 20 mutations (PR20) with >5-fold increased Kd and >4000-fold decreased affinity to the PI darunavir functions. (1) PR20 cleaves, albeit poorly, Gag polyprotein substrates essential for viral maturation. (2) PR20 dimer, which exhibits distinctly enhanced thermal stability, has highly attenuated autoproteolysis, thus likely prolonging its lifetime in vivo. (3) The enhanced stability of PR20 results from stabilization of the monomer fold. Both monomeric PR20(T26A) and dimeric PR20 exhibit Tm values 6-7.5 °C higher than those for their PR counterparts. Two specific mutations in PR20, L33F and L63P at sites of autoproteolysis, increase the Tm of monomeric PR(T26A) by ~8 °C, similar to PR20(T26A). However, without other compensatory mutations as seen in PR20, L33F and L63P substitutions, together, neither restrict autoproteolysis nor significantly reduce binding affinity to darunavir. To determine whether dimer stability contributes to binding affinity for inhibitors, we examined single-chain dimers of PR and PR(D25N) in which the corresponding identical monomer units were covalently linked by GGSSG sequence. Linking of the subunits did not appreciably change the ΔTm on inhibitor binding; thus stabilization by tethering appears to have little direct effect on enhancing inhibitor affinity.

Whiskers-less HIV-protease: a possible way for HIV-1 deactivation

BACKGROUND

Among viral enzymes, the human HIV-1 protease comprises the most interesting target for drug discovery. There are increasing efforts focused on designing more effective inhibitors for HIV-1 protease in order to prevent viral replication in AIDS patients. The frequent and continuous mutation of HIV-1 protease gene creates a formidable obstacle for enzyme inhibition which could not be overcome by the traditional single drug therapy. Nowadays, in vitro and in silico studies of protease inhibition constitute an advanced field in biological researches. In this article, we tried to simulate protease-substrate complexes in different states; a native state and states with whiskers deleted from one and two subunits. Molecular dynamic simulations were carried out in a cubic box filled with explicit water at 37°C and in 1atomsphere of pressure.

RESULTS

Our results showed that whisker truncation of protease subunits causes the dimer structure to decrease in compactness, disrupts substrate-binding site interactions and changes in flap status simultaneously.

CONCLUSIONS

Based on our findings we claim that whisker truncation even when applied to a single subunit, threats dimer association which probably leads to enzyme inactivation. We may postulate that inserting a gene to express truncated protease inside infected cells can interfere with protease dimerization. The resulted proteases would presumably have a combination of native and truncated subunits in their structures which exert no enzyme activities as evidenced by the present work. Our finding may create a new field of research in HIV gene therapy for protease inhibition, circumventing problems of drug resistance.

Elucidating a relationship between conformational sampling and drug resistance in HIV-1 protease

Enzyme targets in rapidly replicating systems, such as retroviruses, commonly respond to drug-selective pressure with mutations arising in the active site pocket that limit inhibitor effectiveness by introducing steric hindrance or by eliminating essential molecular interactions. However, these primary mutations are disposed to compromising pathogenic fitness. Emerging secondary mutations, which are often found outside of the binding cavity, may or can restore fitness while maintaining drug resistance. The accumulated drug pressure selected mutations could have an indirect effect in the development of resistance, such as altering protein flexibility or the dynamics of protein-ligand interactions. Here, we show that accumulation of mutations in a drug-resistant HIV-1 protease (HIV-1 PR) variant, D30N/M36I/A71V, changes the fractional occupancy of the equilibrium conformational sampling ensemble. Correlations are made among populations of the conformational states, namely, closed-like, semiopen, and open-like, with inhibition constants, as well as kinetic parameters. Mutations that stabilize a closed-like conformation correlate with enzymes of lowered activity and with higher affinity for inhibitors, which is corroborated by a further increase in the fractional occupancy of the closed state upon addition of inhibitor or substrate-mimic. Cross-resistance is found to correlate with combinations of mutations that increase the population of the open-like conformations at the expense of the closed-like state while retaining native-like occupancy of the semiopen population. These correlations suggest that at least three states are required in the conformational sampling model to establish the emergence of drug resistance in HIV-1 PR. More importantly, these results shed light on a possible mechanism whereby mutations combine to impart drug resistance while maintaining catalytic activity.

Differential Flap Dynamics in Wild-type and a Drug Resistant Variant of HIV-1 Protease Revealed by Molecular Dynamics and NMR Relaxation

In the rapidly evolving disease of HIV drug resistance readily emerges, nullifying the effectiveness of therapy. Drug resistance has been extensively studied in HIV-1 protease where resistance occurs when the balance between enzyme inhibition and substrate recognition and turn-over is perturbed to favor catalytic activity. Mutations which confer drug resistance can impact the dynamics and structure of both the bound and unbound forms of the enzyme. Flap+ is a multi-drug-resistant variant of HIV-1 protease with a combination of mutations at the edge of the active site, within the active site, and in the flaps (L10I, G48V, I54V, V82A). The impact of these mutations on the dynamics in the unliganded form in comparison with the wild-type protease was elucidated with Molecular Dynamic simulations and NMR relaxation experiments. The comparative analyses from both methods concur in showing that the enzyme's dynamics are impacted by the drug resistance mutations in Flap+ protease. These alterations in the enzyme dynamics, particularly within the flaps, likely modulate the balance between substrate turn-over and drug binding, thereby conferring drug resistance.

Terminal interface conformations modulate dimer stability prior to amino terminal autoprocessing of HIV-1 protease

The HIV-1 protease (PR) mediates its own release (autoprocessing) from the polyprotein precursor, Gag-Pol, flanked by the transframe region (TFR) and reverse transcriptase at its N- and C-termini, respectively. Autoprocessing at the N-terminus of PR mediates stable dimer formation essential for catalytic activity, leading to the formation of infectious virus. An antiparallel β-sheet interface formed by the four N- and C-terminal residues of each subunit is important for dimer stability. Here, we present the first high-resolution crystal structures of model protease precursor-clinical inhibitor (PI darunavir or saquinavir) complexes, revealing varying conformations of the N-terminal flanking (S(-4)FNF(-1)) and interface residues (P(1)QIT(4)). A 180° rotation of the T(4)-L(5) peptide bond is accompanied by a new Q(2)-L(5) hydrogen bond and complete disengagement of PQIT from the β-sheet dimer interface, which may be a feature for intramolecular autoprocessing. This result is consistent with drastically lower thermal stability by 14-20 °C of PI complexes of precursors and the mature PR lacking its PQIT residues (by 18.3 °C). Similar to the TFR-PR precursor, this deletion also results in a darunavir dissociation constant (2 × 10(4))-fold higher and a markedly increased dimer dissociation constant relative to the mature PR. The terminal β-sheet perturbations of the dimeric structure likely account for the drastically poorer inhibition of autoprocessing of TFR-PR relative to the mature PR, even though significant differences in active site-PI interactions in these structures were not observed. The novel conformations of the dimer interface may be exploited to target selectively the protease precursor prior to its N-terminal cleavage.

Inhibition of autoprocessing of natural variants and multidrug resistant mutant precursors of HIV-1 protease by clinical inhibitors

Self-cleavage at the N terminus of HIV-1 protease from the Gag-Pol precursor (autoprocessing) is crucial for stabilizing the protease dimer required for onset of mature-like catalytic activity, viral maturation, and propagation. Among nine clinical protease inhibitors (PIs), darunavir and saquinavir were the most effective in inhibiting wild-type HIV-1 group M precursor autoprocessing, with an IC(50) value of 1-2 μM, 3-5 orders of magnitude higher than their binding affinities to the corresponding mature protease. Accordingly, both group M and N precursor-PI complexes exhibit T(m)s 17-21 °C lower than those of the corresponding mature protease-PI complexes suggestive of markedly reduced stabilities of the precursor dimer-PI ensembles. Autoprocessing of group N (natural variant) and three group M precursors bearing 11-20 mutations associated with multidrug resistance was either weakly responsive or fully unresponsive to inhibitors at concentrations up to a practical limit of approximately 150 μM PI. This observation parallels decreases of up to 8 × 10(3)-fold (e.g., 5 pM to 40 nM) in the binding affinity of darunavir and saquinavir to mature multidrug resistant proteases relative to wild type, suggesting that inhibition of some of these mutant precursors will occur only in the high μM to mM range in extreme PI-resistance, which is an effect arising from coordinated multiple mutations. An extremely darunavir-resistant mutant precursor is more responsive to inhibition by saquinavir. These findings raise the questions whether clinical failure of PI therapy is related to lack of inhibition of autoprocessing and whether specific inhibitors can be designed with low-nM affinity to target autoprocessing.

The L76V drug resistance mutation decreases the dimer stability and rate of autoprocessing of HIV-1 protease by reducing internal hydrophobic contacts

The mature HIV-1 protease (PR) bearing the L76V drug resistance mutation (PR(L76V)) is significantly less stable, with a >7-fold higher dimer dissociation constant (K(d)) of 71 ± 24 nM and twice the sensitivity to urea denaturation (UC(50) = 0.85 M) relative to those of PR. Differential scanning calorimetry showed decreases in T(m) of 12 °C for PR(L76V) in the absence of inhibitors and 5-7 °C in the presence of inhibitors darunavir (DRV), saquinavir (SQV), and lopinavir (LPV), relative to that of PR. Isothermal titration calorimetry gave a ligand dissociation constant of 0.8 nM for DRV, ∼160-fold higher than that of PR, consistent with DRV resistance. Crystal structures of PR(L76V) in complexes with DRV and SQV were determined at resolutions of 1.45-1.46 Å. Compared to the corresponding PR complexes, the mutated Val76 lacks hydrophobic interactions with Asp30, Lys45, Ile47, and Thr74 and exhibits closer interactions with Val32 and Val56. The bound DRV lacks one hydrogen bond with the main chain of Asp30 in PR(L76V) relative to PR, possibly accounting for the resistance to DRV. SQV shows slightly improved polar interactions with PR(L76V) compared to those with PR. Although the L76V mutation significantly slows the N-terminal autoprocessing of the precursor TFR-PR(L76V) to give rise to the mature PR(L76V), the coselected M46I mutation counteracts the effect by enhancing this rate but renders the TFR-PR(M46I/L76V) precursor less responsive to inhibition by 6 μM LPV while preserving inhibition by SQV and DRV. The correlation of lowered stability, higher K(d), and impaired autoprocessing with reduced internal hydrophobic contacts suggests a novel molecular mechanism for drug resistance.

Autocatalytic maturation, physical/chemical properties, and crystal structure of group N HIV-1 protease: relevance to drug resistance

The mature protease from Group N human immunodeficiency virus Type 1 (HIV-1) (PR1(N)) differs in 20 amino acids from the extensively studied Group M protease (PR1(M)) at positions corresponding to minor drug-resistance mutations (DRMs). The first crystal structure (1.09 Å resolution) of PR1(N) with the clinical inhibitor darunavir (DRV) reveals the same overall structure as PR1(M), but with a slightly larger inhibitor-binding cavity. Changes in the 10s loop and the flap hinge propagate to shift one flap away from the inhibitor, whereas L89F and substitutions in the 60s loop perturb inhibitor-binding residues 29-32. However, kinetic parameters of PR1(N) closely resemble those of PR1(M), and calorimetric results are consistent with similar binding affinities for DRV and two other clinical PIs, suggesting that minor DRMs coevolve to compensate for the detrimental effects of drug-specific major DRMs. A miniprecursor (TFR 1-61-PR1(N)) comprising the transframe region (TFR) fused to the N-terminus of PR1(N) undergoes autocatalytic cleavage at the TFR/PR1(N) site concomitant with the appearance of catalytic activity characteristic of the dimeric, mature enzyme. This cleavage is inhibited at an equimolar ratio of precursor to DRV (∼6 μM), which partially stabilizes the precursor dimer from a monomer. However, cleavage at L34/W35 within the TFR, which precedes the TFR 1-61/PR1(N) cleavage at pH ≤ 5, is only partially inhibited. Favorable properties of PR1(N) relative to PR1(M) include its suitability for column fractionation by size under native conditions and >10-fold higher dimer dissociation constant (150 nM). Exploiting these properties may facilitate testing of potential dimerization inhibitors that perturb early precursor processing steps.

Autoprocessing of human immunodeficiency virus type 1 protease miniprecursor fusions in mammalian cells

Background

HIV protease (PR) is a virus-encoded aspartic protease that is essential for viral replication and infectivity. The fully active and mature dimeric protease is released from the Gag-Pol polyprotein as a result of precursor autoprocessing.

Results

We here describe a simple model system to directly examine HIV protease autoprocessing in transfected mammalian cells. A fusion precursor was engineered encoding GST fused to a well-characterized miniprecursor, consisting of the mature protease along with its upstream transframe region (TFR), and small peptide epitopes to facilitate detection of the precursor substrate and autoprocessing products. In HEK 293T cells, the resulting chimeric precursor undergoes effective autoprocessing, producing mature protease that is rapidly degraded likely via autoproteolysis. The known protease inhibitors Darunavir and Indinavir suppressed both precursor autoprocessing and autoproteolysis in a dose-dependent manner. Protease mutations that inhibit Gag processing as characterized using proviruses also reduced autoprocessing efficiency when they were introduced to the fusion precursor. Interestingly, autoprocessing of the fusion precursor requires neither the full proteolytic activity nor the majority of the N-terminal TFR region.

Conclusions

We suggest that the fusion precursors provide a useful system to study protease autoprocessing in mammalian cells, and may be further developed for screening of new drugs targeting HIV protease autoprocessing.

An assay to monitor HIV-1 protease activity for the identification of novel inhibitors in T-cells

The emergence of resistant HIV strains, together with the severe side-effects of existing drugs and lack of development of effective anti-HIV vaccines highlight the need for novel antivirals, as well as innovative methods to facilitate their discovery. Here, we have developed an assay in T-cells to monitor the proteolytic activity of the HIV-1 protease (PR). The assay is based on the inducible expression of HIV-1 PR fused within the Gal4 DNA-binding and transactivation domains. The fusion protein binds to the Gal4 responsive element and activates the downstream reporter, enhanced green fluorescent protein (eGFP) gene only in the presence of an effective PR Inhibitor (PI). Thus, in this assay, eGFP acts as a biosensor of PR activity, making it ideal for flow cytometry based screening. Furthermore, the assay was developed using retroviral technology in T-cells, thus providing an ideal environment for the screening of potential novel PIs in a cell-type that represents the natural milieu of HIV infection. Clones with the highest sensitivity, and robust, reliable and reproducible reporter activity, were selected. The assay is easily adaptable to other PR variants, a multiplex platform, as well as to high-throughput plate reader based assays and will greatly facilitate the search for novel peptide and chemical compound based PIs in T-cells.

Formation of transient dimers by a retroviral protease

Retroviral proteases have been shown previously to be only active as homodimers. They are essential to form the separate and active proteins from the viral precursors. Spumaretroviruses produce separate precursors for Gag and Pol, rather than a Gag and a Gag-Pol precursor. Nevertheless, processing of Pol into a PR (protease)-RT (reverse transcriptase) and integrase is essential in order to obtain infectious viral particles. We showed recently that the PR-RT from a simian foamy virus, as well as the separate PRshort (protease) domain, exhibit proteolytic activities, although only monomeric forms could be detected. In the present study, we demonstrate that PRshort and PR-RT can be inhibited by the putative dimerization inhibitor cholic acid. Various other inhibitors, including darunavir and tipranavir, known to prevent HIV-1 PR dimerization in cells, had no effect on foamy virus protease in vitro. 1H-15N HSQC (heteronuclear single quantum coherence) NMR analysis of PRshort indicates that cholic acid binds in the proposed PRshort dimerization interface and appears to impair formation of the correct dimer. NMR analysis by paramagnetic relaxation enhancement resulted in elevated transverse relaxation rates of those amino acids predicted to participate in dimer formation. Our results suggest transient PRshort homodimers are formed under native conditions but are only present as a minor transient species, which is not detectable by traditional methods.

Highly conserved glycine 86 and arginine 87 residues contribute differently to the structure and activity of the mature HIV-1 protease

The structural and functional role of conserved residue G86 in HIV-1 protease (PR) was investigated by NMR and crystallographic analyses of substitution mutations of glycine to alanine and serine (PR(G86A) and PR(G86S)). While PR(G86S) had undetectable catalytic activity, PR(G86A) exhibited approximately 6000-fold lower catalytic activity than PR. (1)H-(15)N NMR correlation spectra revealed that PR(G86A) and PR(G86S) are dimeric, exhibiting dimer dissociation constants (K(d)) of approximately 0.5 and approximately 3.2 muM, respectively, which are significantly lower than that seen for PR with R87K mutation (K(d) > 1 mM). Thus, the G86 mutants, despite being partially dimeric under the assay conditions, are defective in catalyzing substrate hydrolysis. NMR spectra revealed no changes in the chemical shifts even in the presence of excess substrate, indicating very poor binding of the substrate. Both NMR chemical shift data and crystal structures of PR(G86A) and PR(G86S) in the presence of active-site inhibitors indicated high structural similarity to previously described PR/inhibitor complexes, except for specific perturbations within the active site loop and around the mutation site. The crystal structures in the presence of the inhibitor showed that the region around residue 86 was connected to the active site by a conserved network of hydrogen bonds, and the two regions moved further apart in the mutants. Overall, in contrast to the role of R87 in contributing significantly to the dimer stability of PR, G86 is likely to play an important role in maintaining the correct geometry of the active site loop in the PR dimer for substrate binding and hydrolysis. Proteins 2010. (c) 2009 Wiley-Liss, Inc.

Revealing the dimer dissociation and existence of a folded monomer of the mature HIV-2 protease

Purification and in vitro protein-folding schemes were developed to produce monodisperse samples of the mature wild-type HIV-2 protease (PR2), enabling a comprehensive set of biochemical and biophysical studies to assess the dissociation of the dimeric protease. An E37K substitution in PR2 significantly retards autoproteolytic cleavage during expression. Furthermore, it permits convenient measurement of the dimer dissociation of PR2(E37K) (elevated K(d) approximately 20 nM) by enzyme kinetics. Differential scanning calorimetry reveals a T(m) of 60.5 for PR2 as compared with 65.7 degrees C for HIV-1 protease (PR1). Consistent with weaker binding of the clinical inhibitor darunavir (DRV) to PR2, the T(m) of PR2 increases by 14.8 degrees C in the presence of DRV as compared with 22.4 degrees C for PR1. Dimer interface mutations, such as a T26A substitution in the active site (PR2(T26A)) or a deletion of the C-terminal residues 96-99 (PR2(1-95)), drastically increase the K(d) (>10(5)-fold). PR2(T26A) and PR2(1-95) consist predominantly of folded monomers, as determined by nuclear magnetic resonance (NMR) and size-exclusion chromatography coupled with multiangle light scattering and refractive index measurements (SMR), whereas wild-type PR2 and its active-site mutant PR2(D25N) are folded dimers. Addition of twofold excess active-site inhibitor promotes dimerization of PR2(T26A) but not of PR2(1-95), indicating that subunit interactions involving the C-terminal residues are crucial for dimer formation. Use of SMR and NMR with PR2 facilitates probing for potential inhibitors that restrict protein folding and/or dimerization and, thus, may provide insights for the future design of inhibitors to circumvent drug resistance.

Approaches to the design of HIV protease inhibitors with improved resistance profiles

PURPOSE OF REVIEW

This review describes current approaches to HIV protease inhibitor design, with a focus on improving their profile against drug-resistant mutants. Potential explanations for the flat resistance profile of some potent protease inhibitors and discrepancies between the apparent fold change of potency at the enzyme level and in cell-based assays are discussed.

RECENT FINDINGS

Despite new ideas and a clear rationale for designing inhibitors that bind outside the enzyme active site, all current protease inhibitors with potent antiviral activity target this site. Several bis-tetrahydrofuran-containing compounds including darunavir, brecanavir, GS-8374, and Sequoia protease inhibitors exhibit excellent potency against mutant HIV strains that are resistant to clinically used protease inhibitors. The apparently flat resistance profiles of these and some other protease inhibitors may, at least in part, be explained by their high potency against wild-type enzyme. The substrate envelope and solvent-anchoring hypotheses have been used to design and/or rationalize improved resistance profiles. Traditional approaches yielded a lysine sulfonamide PL-100 with a unique resistance profile.

SUMMARY

Several theories on how to design HIV protease inhibitors with improved resistance profiles have been proposed during the review period. The general concepts that are incorporated into most design strategies include maximizing the interactions with the backbone and conserved side chains of the enzyme while minimizing inhibitor size and maintaining conformational flexibility to allow for modified binding modes.

Interactions of different inhibitors with active-site aspartyl residues of HIV-1 protease and possible relevance to pepsin

The importance of the active site region aspartyl residues 25 and 29 of the mature HIV-1 protease (PR) for the binding of five clinical and three experimental protease inhibitors [symmetric cyclic urea inhibitor DMP323, nonhydrolyzable substrate analog (RPB) and the generic aspartic protease inhibitor acetyl-pepstatin (Ac-PEP)] was assessed by differential scanning calorimetry. DeltaT(m) values, defined as the difference in T(m) for a given protein in the presence and absence of inhibitor, for PR with DRV, ATV, SQV, RTV, APV, DMP323, RPB, and Ac-PEP are 22.4, 20.8, 19.3, 15.6, 14.3, 14.7, 8.7, and 6.5 degrees C, respectively. Binding of APV and Ac-PEP is most sensitive to the D25N mutation, as shown by DeltaT(m) ratios [DeltaT(m)(PR)/DeltaT(m)(PR(D25N))] of 35.8 and 16.3, respectively, whereas binding of DMP323 and RPB (DeltaT(m) ratios of 1-2) is least affected. Binding of the substrate-like inhibitors RPB and Ac-PEP is nearly abolished (DeltaT(m)(PR)/DeltaT(m)(PR(D29N)) > or = 44) by the D29N mutation, whereas this mutation only moderately affects binding of the smaller inhibitors (DeltaT(m) ratios of 1.4-2.2). Of the nine FDA-approved clinical HIV-1 protease inhibitors screened, APV, RTV, and DRV competitively inhibit porcine pepsin with K(i) values of 0.3, 0.6, and 2.14 microM, respectively. DSC results were consistent with this relatively weak binding of APV (DeltaT(m) 2.7 degrees C) compared with the tight binding of Ac-PEP (DeltaT(m) > or = 17 degrees C). Comparison of superimposed structures of the PR/APV complex with those of PR/Ac-PEP and pepsin/pepstatin A complexes suggests a role for Asp215, Asp32, and Ser219 in pepsin, equivalent to Asp25, Asp25', and Asp29 in PR in the binding and stabilization of the pepsin/APV complex.

Uncoupling human immunodeficiency virus type 1 Gag and Pol reading frames: role of the transframe protein p6* in viral replication

Apart from its regulatory role in protease (PR) activation, little is known about the function of the human immunodeficiency virus type 1 transframe protein p6* in the virus life cycle. p6* is located between the nucleocapsid and PR domains in the Gag-Pol polyprotein precursor and is cleaved by PR during viral maturation. We have recently reported that the central region of p6* can be extensively mutated without abolishing viral infectivity and replication in vitro. However, mutagenesis of the entire p6*-coding sequence in the proviral context is not feasible without affecting the superimposed frameshift signal or the overlapping p1-p6(gag) sequences. To overcome these limitations, we created a novel NL4-3-derived provirus by displacing the original frameshift signal to the 3' end of the gag gene, thereby uncoupling the p6* gene sequence from the p1-p6(gag) reading frame. The resulting virus (AL) proved to be replication competent in different cell cultures and thus represents an elegant tool for detailed analysis of p6* function. Hence, extensive deletions or substitutions were introduced into the p6* gene sequence of the AL provirus, and effects on particle release, protein processing, and viral infectivity were evaluated. Interestingly, neither the deletion of 63% of all p6* residues nor the partial substitution by a heterologous sequence affected virus growth and infectivity, suggesting that p6* is widely dispensable for viral in vitro replication. However, the insertion of a larger reporter sequence interfered with virus production and maturation, implying that the length or conformation of this spacer region might be critical for p6* function.

Analysis and characterization of dimerization inhibition of a multi-drug-resistant human immunodeficiency virus type 1 protease using a novel size-exclusion chromatographic approach

Active-site inhibitors of HIV-1 PR (protease) block viral replication by preventing viral maturation. However, HIV-1 often develops resistance to active-site inhibitors through multiple mutations in PR and therefore recent efforts have focused on inhibiting PR dimerization as an alternative approach. Dimerization inhibitors have been identified using kinetic analysis, but additional characterization of the effect of these inhibitors on PR by physical methods has been difficult. In the present study, we identified a PR(MDR) (multi-drug-resistant HIV-1 PR) that was highly resistant to autoproteolysis. Using this PR and a novel size-exclusion chromatographic approach that incorporated fluorescence and MS detection, we were able to demonstrate inhibition of dimerization using P27 (peptide 27), a peptide dimerization inhibitor of PR previously identified on the basis of kinetic analysis. Incubation of PR(MDR) with P27, or other dimerization inhibitors, led to a dose- and time-dependent formation of PR monomers based on the change in elution time by size exclusion and its similar elution time to engineered forms of monomeric PR, namely PR(T26A) and glutathionylated PR. In contrast, incubation of PR(MDR) with a potent active-site inhibitor did not change the elution time for the PR(MDR) dimer. The monomeric PR induced by P27 had fluorescent characteristics which were consistent with unfolded PR. Structure-activity studies identified the active regions of P27 and experiments were performed to examine the effect of other dimerization inhibitors on PR. The present study is the first characterization of dimerization inhibition of PR(MDR), a prime target for these inhibitors, using a novel size-exclusion chromatographic approach.

The folding free-energy surface of HIV-1 protease: insights into the thermodynamic basis for resistance to inhibitors

Spontaneous mutations at numerous sites distant from the active site of human immunodeficiency virus type 1 protease enable resistance to inhibitors while retaining enzymatic activity. As a benchmark for probing the effects of these mutations on the conformational adaptability of this dimeric beta-barrel protein, the folding free-energy surface of a pseudo-wild-type variant, HIV-PR(*), was determined by a combination of equilibrium and kinetic experiments on the urea-induced unfolding/refolding reactions. The equilibrium unfolding reaction was well described by a two-state model involving only the native dimeric form and the unfolded monomer. The global analysis of the kinetic folding mechanism reveals the presence of a fully folded monomeric intermediate that associates to form the native dimeric structure. Independent analysis of a stable monomeric version of the protease demonstrated that a small-amplitude fluorescence phase in refolding and unfolding, not included in the global analysis of the dimeric protein, reflects the presence of a transient intermediate in the monomer folding reaction. The partially folded and fully folded monomers are only marginally stable with respect to the unfolded state, and the dimerization reaction provides a modest driving force at micromolar concentrations of protein. The thermodynamic properties of this system are such that mutations can readily shift the equilibrium from the dimeric native state towards weakly folded states that have a lower affinity for inhibitors but that could be induced to bind to their target proteolytic sites. Presumably, subsequent secondary mutations increase the stability of the native dimeric state in these variants and, thereby, optimize the catalytic properties of the resistant human immunodeficiency virus type 1 protease.

Visualizing transient events in amino-terminal autoprocessing of HIV-1 protease

HIV-1 protease processes the Gag and Gag-Pol polyproteins into mature structural and functional proteins, including itself, and is therefore indispensable for viral maturation. The mature protease is active only as a dimer with each subunit contributing catalytic residues. The full-length transframe region protease precursor appears to be monomeric yet undergoes maturation via intramolecular cleavage of a putative precursor dimer, concomitant with the appearance of mature-like catalytic activity. How such intramolecular cleavage can occur when the amino and carboxy termini of the mature protease are part of an intersubunit beta-sheet located distal from the active site is unclear. Here we visualize the early events in N-terminal autoprocessing using an inactive mini-precursor with a four-residue N-terminal extension that mimics the transframe region protease precursor. Using paramagnetic relaxation enhancement, a technique that is exquisitely sensitive to the presence of minor species, we show that the mini-precursor forms highly transient, lowly populated (3-5%) dimeric encounter complexes that involve the mature dimer interface but occupy a wide range of subunit orientations relative to the mature dimer. Furthermore, the occupancy of the mature dimer configuration constitutes a very small fraction of the self-associated species (accounting for the very low enzymatic activity of the protease precursor), and the N-terminal extension makes transient intra- and intersubunit contacts with the substrate binding site and is therefore available for autocleavage when the correct dimer orientation is sampled within the encounter complex ensemble.

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

Novel macromolecular inhibitors of human immunodeficiency virus-1 protease

An intracellularly expressed defective human immunodeficiency virus type-1 (HIV-1) protease (PR) monomer could function as a dominant-negative inhibitor of the enzyme that requires dimerization for activity. Based on in silico studies, two mutant PRs harboring hydrophilic mutations, capable of forming favorable inter- and intra-subunit interactions, were selected: PR(RE) containing Asp25Arg and Gly49Glu mutations, and PR(RER) containing an additional Ile50Arg mutation. The mutants were expressed and tested by PR assays, nuclear magnetic resonance (NMR) and cell culture experiments. The mutant PRs showed dose-dependent inhibition of the wild-type PR in a fluorescent microtiter plate PR assay. Furthermore, both mutants were retained by hexahistidine-tagged wild-type HIV-1 PR immobilized on nickel-chelate affinity resin. For the first time, heterodimerization between wild-type and dominant-negative mutant PRs were also demonstrated by NMR spectroscopy. (1)H-(15)N Heteronuclear Single Quantum Coherence NMR spectra showed that although PR(RE) has a high tendency to aggregate, PR(RER) exists mainly as a folded monomer at 25-35 microM concentration, but in the presence of wild-type PR in a ratio of 1:1, heterodimerization occurs with both mutants. While the recombinant virus containing the PR(RE) sequence showed only very low level of expression, expression of the viral proteins of the virus with the PR(RER) sequence was comparable with that of the wild-type. In cell culture experiments, infectivity of viral particles containing PR(RER) protein was reduced by 82%, at mutant to wild-type infective DNA ratio of 2:1.

Atomistic simulations of the HIV-1 protease folding inhibition

Biochemical experiments have recently revealed that the p-S8 peptide, with an amino-acid sequence identical to the conserved fragment 83-93 (S8) of the HIV-1 protease, can inhibit catalytic activity of the enzyme by interfering with protease folding and dimerization. In this study, we introduce a hierarchical modeling approach for understanding the molecular basis of the HIV-1 protease folding inhibition. Coarse-grained molecular docking simulations of the flexible p-S8 peptide with the ensembles of HIV-1 protease monomers have revealed structurally different complexes of the p-S8 peptide, which can be formed by targeting the conserved segment 24-34 (S2) of the folding nucleus (folding inhibition) and by interacting with the antiparallel termini beta-sheet region (dimerization inhibition). All-atom molecular dynamics simulations of the inhibitor complexes with the HIV-1 PR monomer have been independently carried out for the predicted folding and dimerization binding modes of the p-S8 peptide, confirming the thermodynamic stability of these complexes. Binding free-energy calculations of the p-S8 peptide and its active analogs are then performed using molecular dynamics trajectories of the peptide complexes with the HIV-1 PR monomers. The results of this study have provided a plausible molecular model for the inhibitor intervention with the HIV-1 PR folding and dimerization and have accurately reproduced the experimental inhibition profiles of the active folding inhibitors.

Effect of the active site D25N mutation on the structure, stability, and ligand binding of the mature HIV-1 protease

All aspartic proteases, including retroviral proteases, share the triplet DTG critical for the active site geometry and catalytic function. These residues interact closely in the active, dimeric structure of HIV-1 protease (PR). We have systematically assessed the effect of the D25N mutation on the structure and stability of the mature PR monomer and dimer. The D25N mutation (PR(D25N)) increases the equilibrium dimer dissociation constant by a factor >100-fold (1.3 +/- 0.09 microm) relative to PR. In the absence of inhibitor, NMR studies reveal clear structural differences between PR and PR(D25N) in the relatively mobile P1 loop (residues 79-83) and flap regions, and differential scanning calorimetric analyses show that the mutation lowers the stabilities of both the monomer and dimer folds by 5 and 7.3 degrees C, respectively. Only minimal differences are observed in high resolution crystal structures of PR(D25N) complexed to darunavir (DRV), a potent clinical inhibitor, or a non-hydrolyzable substrate analogue, Ac-Thr-Ile-Nle-r-Nle-Gln-Arg-NH(2) (RPB), as compared with PR.DRV and PR.RPB complexes. Although complexation with RPB stabilizes both dimers, the effect on their T(m) is smaller for PR(D25N) (6.2 degrees C) than for PR (8.7 degrees C). The T(m) of PR(D25N).DRV increases by only 3 degrees C relative to free PR(D25N), as compared with a 22 degrees C increase for PR.DRV, and the mutation increases the ligand dissociation constant of PR(D25N).DRV by a factor of approximately 10(6) relative to PR.DRV. These results suggest that interactions mediated by the catalytic Asp residues make a major contribution to the tight binding of DRV to PR.