Novel strategies for targeting the dimerization interface of HIV protease with cross-linked interfacial peptides

Computer-aided design and discovery of protein-protein interaction inhibitors as agents for anti-HIV therapy

Protein-protein interactions (PPI) are involved in most of the essential processes that occur in organisms. In recent years, PPI have become the object of increasing attention in drug discovery, particularly for anti-HIV drugs. Although the use of combinations of existing drugs, termed highly active antiretroviral therapy (HAART), has revolutionized the treatment of HIV/AIDS, problems with these agents, such as the rapid emergence of drug-resistant HIV-1 mutants and serious adverse effects, have highlighted the need for further discovery of new drugs and new targets. Numerous investigations have shown that PPI play a key role in the virus's life cycle and that blocking or modulating them has a significant therapeutic potential. Here we summarize the recent progress in computer-aided design of PPI inhibitors, mainly focusing on the selection of the drug targets (HIV enzymes and virus entry machinery) and the utilization of peptides and small molecules to prevent a variety of protein-protein interactions (viral-viral or viral-host) that play a vital role in the progression of HIV infection.

Mechanism of dissociative inhibition of HIV protease and its autoprocessing from a precursor

Dimerization is indispensible for release of the human immunodeficiency virus protease (PR) from its precursor (Gag-Pol) and ensuing mature-like catalytic activity that is crucial for virus maturation. We show that a single-chain Fv fragment (scFv) of a previously reported monoclonal antibody (mAb1696), which recognizes the N-terminus of PR, dissociates a dimeric mature D25N PR mutant with an enhanced dimer dissociation constant (K(d)) in the sub-micromolar range to form predominantly a monomer-scFv complex at a 1:1 ratio, along with small (5-10%) amounts of a dimer-scFv complex. Enzyme kinetics indicate a mixed mechanism of inhibition of the wild-type PR, which exhibits a K(d)<10nM, with effects both on K(m) and k(cat) at an scFv-to-PR ratio of 10:1. ScFv binds to the N-terminal peptide P(1)QITLW(6) of PR and to PR monomers with dissociation constants of ≤30 nM and ~100 nM, respectively. Consistent with an ~400-fold increase in the dissociation of the antibody (K(Ab)) on even addition of an acetyl group to P(1) of the peptide, the antibody fails to inhibit N-terminal autoprocessing of the PR from a model precursor (at ~5 μM). However, subsequent to this cleavage, it sequesters the PR, thus blocking autoprocessing at its C-terminus. A second monoclonal antibody [PRM1 (human monoclonal antibody to PR)], which recognizes part of the flap region (residues 41-47) of the mature PR and its precursor, does not inhibit autoprocessing and ensuing catalytic activity. However, its failure to recognize drug-resistant clinical mutants of PR may be beneficial to monitor the selection of mutations in this region under drug pressure.

Effects of HIV-1 protease on cellular functions and their potential applications in antiretroviral therapy

Human Immunodeficiency Virus Type 1 (HIV-1) protease inhibitors (PIs) are the most potent class of drugs in antiretroviral therapies. However, viral drug resistance to PIs could emerge rapidly thus reducing the effectiveness of those drugs. Of note, all current FDA-approved PIs are competitive inhibitors, i.e., inhibitors that compete with substrates for the active enzymatic site. This common inhibitory approach increases the likelihood of developing drug resistant HIV-1 strains that are resistant to many or all current PIs. Hence, new PIs that move away from the current target of the active enzymatic site are needed. Specifically, allosteric inhibitors, inhibitors that prohibit PR enzymatic activities through non-competitive binding to PR, should be sought. Another common feature of current PIs is they were all developed based on the structure-based design. Drugs derived from a structure-based strategy may generate target specific and potent inhibitors. However, this type of drug design can only target one site at a time and drugs discovered by this method are often associated with strong side effects such as cellular toxicity, limiting its number of target choices, efficacy, and applicability. In contrast, a cell-based system may provide a useful alternative strategy that can overcome many of the inherited shortcomings associated with structure-based drug designs. For example, allosteric PIs can be sought using a cell-based system without considering the site or mechanism of inhibition. In addition, a cell-based system can eliminate those PIs that have strong cytotoxic effect. Most importantly, a simple, economical, and easy-to-maintained eukaryotic cellular system such as yeast will allow us to search for potential PIs in a large-scaled high throughput screening (HTS) system, thus increasing the chances of success. Based on our many years of experience in using fission yeast as a model system to study HIV-1 Vpr, we propose the use of fission yeast as a possible surrogate system to study the effects of HIV-1 protease on cellular functions and to explore its utility as a HTS system to search for new PIs to battle HIV-1 resistant strains.

Multiple routes and milestones in the folding of HIV-1 protease monomer

Proteins fold on a time scale incompatible with a mechanism of random search in conformational space thus indicating that somehow they are guided to the native state through a funneled energetic landscape. At the same time the heterogeneous kinetics suggests the existence of several different folding routes. Here we propose a scenario for the folding mechanism of the monomer of HIV–1 protease in which multiple pathways and milestone events coexist. A variety of computational approaches supports this picture. These include very long all-atom molecular dynamics simulations in explicit solvent, an analysis of the network of clusters found in multiple high-temperature unfolding simulations and a complete characterization of free-energy surfaces carried out using a structure-based potential at atomistic resolution and a combination of metadynamics and parallel tempering. Our results confirm that the monomer in solution is stable toward unfolding and show that at least two unfolding pathways exist. In our scenario, the formation of a hydrophobic core is a milestone in the folding process which must occur along all the routes that lead this protein towards its native state. Furthermore, the ensemble of folding pathways proposed here substantiates a rational drug design strategy based on inhibiting the folding of HIV–1 protease.

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.

Chemogenomics with protein secondary-structure mimetics

During molecular recognition of proteins in biological systems, helices, reverse turns, and beta-sheets are dominant motifs. Often there are therapeutic reasons for blocking such recognition sites, and significant progress has been made by medicinal chemists in the design and synthesis of semirigid molecular scaffolds on which to display amino acid side chains. The basic premise is that preorganization of the competing ligand enhances the binding affinity and potential selectivity of the inhibitor. In this chapter, current progress in these efforts is reviewed.

Sidechain-linked inhibitors of HIV-1 protease dimerization

There is a great need for alternative modes of inhibition for the design of anti-HIV therapies, due to the increased resistance of HIV to currently approved drugs. A novel strategy for generating potent dimerization inhibitors of HIV-1 protease is described based on sidechain-linked interfacial peptides. In a number of cases the activity of these agents against HIV-1 protease was found to be among the most potent reported, with inhibitory constants in the low nM range.

Molecular dynamics simulations of HIV‐1 protease monomer: Assembly of N‐terminus and C‐terminus into β‐sheet in water solution

HIV-1 protease (HIV-PR) consists of two identical subunits that are united together through a four-stranded antiparallel beta-sheet formed of the peptide termini of each monomer. Since the active site exists only in the dimer, a strategy that is attracting more and more attention in inhibitor design and which may overcome the serious drug resistance caused by competitive inhibitors is to block the peptide termini of the monomer, thereby interfering with formation of the active dimer. In the present work, we performed several extensive molecular dynamics (MD) simulations of the HIV-PR monomer in water to illustrate its solvated conformation and dynamics behavior. We found that the peptide termini usually assembled into beta-sheet after several nanoseconds' simulation, and became much less flexible. This beta-sheet is stabilized by intramolecular interactions and is not easily disaggregated under the present MD simulation conditions. This transformation may be an important transition during the relaxing and equilibrating of the HIV-PR monomer in aqueous solution, and the terminal beta-sheet may be one of the major conformations of the solvated HIV-PR monomer termini in water. This work may provide new insights into the dynamics behavior and dimerization mechanism of HIV-PR, and more significantly, offer a more rational receptor model for the design and discovery of novel dimerization inhibitors than crystalline structures.

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

An experimental protocol for folding the mature human immunodeficiency virus-1 (HIV-1) protease is presented that facilitates NMR studies at a low protein concentration of approximately 20 micoM. Under these conditions, NMR spectra show that the mature protease lacking its terminal beta-sheet residues 1-4 and 96-99 (PR(5-95)) exhibits a stable monomer fold spanning the region 10-90 that is similar to that of the single subunit of the wild-type dimer and the dimer bearing a D25N mutation (PR(D25N)). Urea-induced unfolding monitored both by changes in (1)H-(15)N heteronuclear single quantum correlation spectra and by protein fluorescence indicates that although PR(5-95) monomer displays a transition profile similar to that of the PR(D25N) dimer (50% unfolded (U(50)) = approximately 1.9 M), extending the protease with 4 residues (SFNF) of its N-terminally flanking sequence in the Gag-Pol precursor ((SFNF)PR(D25N)) decreases the stability of the fold (U(50) = approximately 1.5 M). Assigned backbone chemical shifts were used to elucidate differences in the stability of the PR(T26A) (U(50) = 2.5 M) and (SFNF)PR(D25N) monomers and compared with PR(D25N/T26A) monomer. Discernible differences in the backbone chemical shifts were observed for N-terminal protease residues 3-6 of (SFNF)PR(D25N) that may relate to the increase in the equilibrium dissociation constant (K(d)) and the very low catalytic activity of the protease prior to its autoprocessing at its N terminus from the Gag-Pol precursor.

Differential susceptibility of Plasmodium falciparum versus yeast and mammalian enolases to dissociation into active monomers

In the past, several unsuccessful attempts have been made to dissociate homodimeric enolases into their active monomeric forms. The main objective of these studies had been to understand whether intersubunit interactions are essential for the catalytic and structural stability of enolases. Further motivation to investigate the properties of monomeric enolase has arisen from several recent reports on the involvement of enolase in diverse nonglycolytic (moonlighting) functions, where it may occur in monomeric form. Here, we report successful dissociation of dimeric enolases from Plasmodium falciparum, yeast and rabbit muscle into active and isolatable monomers. Dimeric enolases could be dissociated into monomers by high concentrations ( approximately 250 mm) of imidazole and/or hydrogen ions. Two forms were separated using Superdex-75 gel filtration chromatography. A detailed comparison of the kinetic and structural properties of monomeric and dimeric forms of recombinant P. falciparum enolase showed differences in specific activity, salt-induced inhibition and inactivation, thermal stability, etc. Furthermore, we found that enolases from the three species differ in their dimer dissociation profiles. Specifically, on challenge with imidazole, Mg(II) protected the enolases of yeast and rabbit muscle but not of P. falciparum from dissociation. The observed differential stability of the P. falciparum enolase dimer interface with respect to mammalian enolases could be exploited to selectively dissociate the dimeric parasite enzyme into its catalytically inefficient, thermally unstable monomeric form. Thus enolase could be a novel therapeutic target for malaria.

Ddi1, a eukaryotic protein with the retroviral protease fold

Retroviral aspartyl proteases are homodimeric, whereas eukaryotic aspartyl proteases tend to be large, monomeric enzymes with 2-fold internal symmetry. It has been proposed that contemporary monomeric aspartyl proteases evolved by gene duplication and fusion from a primordial homodimeric enzyme. Recent sequence analyses have suggested that such "fossil" dimeric aspartyl proteases are still encoded in the eukaryotic genome. We present evidence for retention of a dimeric aspartyl protease in eukaryotes. The X-ray crystal structure of a domain of the Saccharomyces cerevisiae protein Ddi1 shows that it is a dimer with a fold similar to that of the retroviral proteases. Furthermore, the double Asp-Thr-Gly-Ala amino acid sequence motif at the active site of HIV protease is found with identical geometry in the Ddi1 structure. However, the putative substrate binding groove is wider in Ddi1 than in the retroviral proteases, suggesting that Ddi1 accommodates bulkier substrates. Ddi1 belongs to a family of proteins known as the ubiquitin receptors, which have in common the ability to bind ubiquitinated substrates and the proteasome. Ubiquitin receptors contain an amino-terminal ubiquitin-like (UBL) domain and a carboxy-terminal ubiquitin-associated (UBA) domain, but Ddi1 is the only representative in which the UBL and UBA domains flank an aspartyl protease-like domain. The remarkable structural similarity between the central domain of Ddi1 and the retroviral proteases, in the global fold and in active-site detail, suggests that Ddi1 functions proteolytically during regulated protein turnover in the cell.

Computational mutagenesis studies of protein structure‐function correlations

Topological scores, measures of sequence-structure compatibility, are calculated for all 1,881 single point mutants of the human immunodeficiency virus (HIV)-1 protease using a four-body statistical potential function based on Delaunay tessellation of protein structure. Comparison of the mutant topological score data with experimental data from alanine scan studies specifically on the dimer interface residues supports previous findings that 1) L97 and F99 contribute greatly to the Gibbs energy of HIV-1 protease dimerization, 2) Q2 and T4 contribute the least toward the Gibbs energy, and 3) C-terminal residues are more sensitive to mutations than those at the N-terminus. For a more comprehensive treatment of the relationship between protease structure and function, mutant topological scores are compared with the activity levels for a set of 536 experimentally synthesized protease mutants, and a significant correlation is observed. Finally, this structure-function correlation is similarly identified by examining model systems consisting of 2,015 single point mutants of bacteriophage T4 lysozyme as well as 366 single point mutants of HIV-1 reverse transcriptase and is hypothesized to be a property generally applicable to all proteins.

Rapid Synthesis and In Situ Screening of Potent HIV-1 Protease Dimerization Inhibitors

A library of dimerization inhibitors of HIV-1 protease is described based on crosslinked interfacial peptides. The 54 component library was designed to contain two modifications to the starting structure, one each in the Northern and Southern fragments. A rapid synthesis and in situ screening method in microtiter plates was developed to facilitate the generation and evaluation of the library members. More than 90% of the doubly modified agents were more potent than their respective singly mutated parent compounds, and five of the most potent dimerization inhibitors of HIV-1 protease described to date were identified. The free energy of binding for the combined two modifications was generally found to be additive, demonstrating the predictive value of earlier libraries.

Dimerization of beta-site beta-amyloid precursor protein-cleaving enzyme

Cleavage of the beta-amyloid precursor protein (APP) by the aspartyl protease beta-site APP-cleaving enzyme (BACE) is the first step in the generation of the amyloid beta-peptide, which is deposited in the brain of Alzheimer's disease patients. Whereas the subsequent cleavage by gamma-secretase was shown to originate from the cooperation of a multicomponent complex, it is currently unknown whether in a cellular environment BACE is enzymatically active as a monomer or in concert with other proteins. Using blue native gel electrophoresis we found that endogenous and overexpressed BACE has a molecular mass of 140 kDa instead of the expected mass of 70 kDa under denaturing conditions. This suggests that under native conditions BACE exists as a homodimer. Homodimerization was confirmed by co-immunoprecipitation of full-length BACE carrying different epitope tags. In contrast, the soluble active BACE ectodomain was exclusively present as a monomer both under native and denaturing conditions. A domain analysis revealed that the BACE ectodomain dimerized as long as it was attached to the membrane, whereas the cytoplasmic domain and the transmembrane domain were dispensable for dimerization. By adding a KKXX-endoplasmic reticulum retention signal to BACE, we demonstrate that dimerization of BACE occurs already before full maturation and pro-peptide cleavage. Furthermore, kinetic analysis of the purified native BACE dimer revealed a higher affinity and turnover rate in comparison to the monomeric soluble BACE. Dimerization of BACE might, thus, facilitate binding and cleavage of physiological substrates.

A unidirectional crosslinking strategy for HIV-1 protease dimerization inhibitors

A novel strategy to identify potent HIV-1 protease dimerization inhibitors was developed using 12-aminododecanoic acid as a tether to crosslink interfacial peptides. The directionality of the southern peptide was changed from N-->C to C-->N as compared to previously reported inhibitors. The terminal amine of the southern peptide and side chains were further diversified to find essential functional groups for dimerization inhibition of HIV-1 protease.

Small-Molecule Dimerization Inhibitors of Wild-Type and Mutant HIV Protease: A Focused Library Approach

We demonstrate that a focused library based on truncated, cross-linked interfacial peptides of HIV-1 protease produces effective dimerization inhibitors of the enzyme. By combining individual changes of the library into a single compound, we obtained a significantly more potent agent and found that an additive increase in inhibitor efficacy was obtained. The good activity of library members against an active-site drug-resistant protease mutant bodes well for dimerization inhibition as a complementary method to targeting the active site.

The folding and dimerization of HIV-1 protease: evidence for a stable monomer from simulations

HIV-1 protease (PR) is a major drug target in combating AIDS, as it plays a key role in maturation and replication of the virus. Six FDA-approved drugs are currently in clinical use, all designed to inhibit enzyme activity by blocking the active site, which exists only in the dimer. An alternative inhibition mode would be required to overcome the emergence of drug-resistance through the accumulation of mutations. This might involve inhibiting the formation of the dimer itself. Here, the folding of HIV-1 PR dimer is studied with several simulation models appropriate for folding mechanism studies. Simulations with an off-lattice Gō-model, which corresponds to a perfectly funneled energy landscape, indicate that the enzyme is formed by association of structured monomers. All-atom molecular dynamics simulations strongly support the stability of an isolated monomer. The conjunction of results from a model that focuses on the protein topology and a detailed all-atom force-field model suggests, in contradiction to some reported equilibrium denaturation experiments, that monomer folding and dimerization are decoupled. The simulation result is, however, in agreement with the recent NMR detection of folded monomers of HIV-1 PR mutants with a destabilized interface. Accordingly, the design of dimerization inhibitors should not focus only on the flexible N and C termini that constitute most of the dimer interface, but also on other structured regions of the monomer. In particular, the relatively high phi values for residues 23-35 and 79-87 in both the folding and binding transition states, together with their proximity to the interface, highlight them as good targets for inhibitor design.

Interfacial peptide inhibitors of HIV-1 integrase activity and dimerization

Peptides derived from the interfacial region of dimeric HIV-1 integrase were evaluated as inhibitors of integrase's 3'-endonuclease activity. Three peptides were found to be moderately potent inhibitors with IC(50) values in the low micromolar range. The mode of inhibition was probed through protein crosslinking experiments. Active interfacial peptides were found to inhibit crosslinking of the dimeric form of integrase. Interfacial peptides that were poor inhibitors had no effect on integrase crosslinking.