Computational proteomics analysis of HIV-1 protease interactome

HIV-1 protease is a small homodimeric enzyme that ensures maturation of HIV virions by cleaving the viral precursor Gag and Gag-Pol polyproteins into structural and functional elements. The cleavage sites in the viral polyproteins share neither sequence homology nor binding motif and the specificity of the HIV-1 protease is therefore only partially understood. Using an extensive data set collected from 16 years of HIV proteome research we have here created a general and predictive rule-based model for HIV-1 protease specificity based on rough sets. We demonstrate that HIV-1 protease specificity is much more complex than previously anticipated, which cannot be defined based solely on the amino acids at the substrate's scissile bond or by any other single substrate amino acid position only. Our results show that the combination of at least three particular amino acids is needed in the substrate for a cleavage event to occur. Only by combining and analyzing massive amounts of HIV proteome data it was possible to discover these novel and general patterns of physico-chemical substrate cleavage determinants. Our study is an example how computational biology methods can advance the understanding of the viral interactomes.     

The dimer interfaces of protease and extra-protease domains influence the activation of protease and the specificity of GagPol cleavage

Activation of the human immunodeficiency virus type 1 (HIV-1) protease is an essential step in viral replication. As is the case for all retroviral proteases, enzyme activation requires the formation of protease homodimers. However, little is known about the mechanisms by which retroviral proteases become active within their precursors. Using an in vitro expression system, we have examined the determinants of activation efficiency and the order of cleavage site processing for the protease of HIV-1 within the full-length GagPol precursor. Following activation, initial cleavage occurs between the viral p2 and nucleocapsid proteins. This is followed by cleavage of a novel site located in the transframe domain. Mutational analysis of the dimer interface of the protease produced differential effects on activation and specificity. A subset of mutations produced enhanced cleavage at the amino terminus of the protease, suggesting that, in the wild-type precursor, cleavages that liberate the protease are a relatively late event. Replacement of the proline residue at position 1 of the protease dimer interface resulted in altered cleavage of distal sites and suggests that this residue functions as a cis-directed specificity determinant. In summary, our studies indicate that interactions within the protease dimer interface help determine the order of precursor cleavage and contribute to the formation of extended-protease intermediates. Assembly domains within GagPol outside the protease domain also influence enzyme activation.     

Enzymatic activity of a synthetic 99 residue protein corresponding to the putative HIV-1 protease

A protein corresponding to the putative protease of the human immunodeficiency virus 1 (HIV-1) has been prepared by total chemical synthesis. This 99 residue synthetic enzyme showed specific proteolytic activity on fragments of the natural gag precursor and on synthetic peptide substrates, two of which released fragments corresponding to the N terminus and C terminus of the protease molecule itself. The observed substrate specificity was not restricted to cleavage at Phe/Tyr-Pro bonds. Inhibition studies provided direct evidence that the HIV-1 protease belongs to the family of aspartic proteases. The availability of the HIV-1 protease as a defined molecular species has important implications for the design of specific inhibitors that do not interfere with the host cell metabolism as a possible route to antiviral agents against acquired immunodeficiency syndrome (AIDS).     

Comprehensive bioinformatic analysis of the specificity of human immunodeficiency virus type 1 protease

Rapidly developing viral resistance to licensed human immunodeficiency virus type 1 (HIV-1) protease inhibitors is an increasing problem in the treatment of HIV-infected individuals and AIDS patients. A rational design of more effective protease inhibitors and discovery of potential biological substrates for the HIV-1 protease require accurate models for protease cleavage specificity. In this study, several popular bioinformatic machine learning methods, including support vector machines and artificial neural networks, were used to analyze the specificity of the HIV-1 protease. A new, extensive data set (746 peptides that have been experimentally tested for cleavage by the HIV-1 protease) was compiled, and the data were used to construct different classifiers that predicted whether the protease would cleave a given peptide substrate or not. The best predictor was a nonlinear predictor using two physicochemical parameters (hydrophobicity, or alternatively polarity, and size) for the amino acids, indicating that these properties are the key features recognized by the HIV-1 protease. The present in silico study provides new and important insights into the workings of the HIV-1 protease at the molecular level, supporting the recent hypothesis that the protease primarily recognizes a conformation rather than a specific amino acid sequence. Furthermore, we demonstrate that the presence of 1 to 2 lysine residues near the cleavage site of octameric peptide substrates seems to prevent cleavage efficiently, suggesting that this positively charged amino acid plays an important role in hindering the activity of the HIV-1 protease.     

Autoprocessing of HIV-1 protease is tightly coupled to protein folding.

In the Gag-Pol polyprotein of HIV-1, the 99-amino acid protease is flanked at its N-terminus by a transframe region (TFR) composed of the transframe octapeptide (TFP) and 48 amino acids of the p6pol, separated by a protease cleavage site. The intact precursor (TFP-p6pol-PR) has very low dimer stability relative to that of the mature enzyme and exhibits negligible levels of stable tertiary structure. Thus, the TFR functions by destabilizing the native structure, unlike proregions found in zymogen forms of monomeric aspartic proteases. Cleavage at the p6pol-PR site to release a free N-terminus of protease is concomitant with the appearance of enzymatic activity and formation of a stable tertiary structure that is characteristic of the mature protease as demonstrated by nuclear magnetic resonance. The release of the mature protease from the precursor can either occur in two steps at pH values of 4 to 6 or in a single step above pH 6. The mature protease forms a dimer through a four-stranded beta-sheet at the interface. Residues 1-4 of the mature protease from each subunit constitute the outer strands of the beta-sheet, and are essential for maintaining the stability of the free protease but are not a prerequisite for the formation of tertiary structure and catalytic activity. Our experimental results provide the basis for the model proposed here for the regulation of the HIV-1 protease in the viral replication cycle.     

Analysis of retroviral protease cleavage sites reveals two types of cleavage sites and the structural requirements of the P1 amino acid

Retroviruses encode a protease which cleaves the viral Gag and Gag/Pol protein precursors into mature products. To understand the target sequence specificity of the viral protease, the amino acid sequences from 46 known processing sites from 10 diverse retroviruses were compared. Sequence preference was evident in positions P4 through P3' when compared to flanking sequences. Approximately 80% of all cleavage site sequences could be grouped into two classes based on the sequence composition flanking the scissile bond. The sequences at the amino-terminal cleavage site of the major capsid protein of Gag is always a member of one of the two classes while the carboxyl-terminal cleavage site is of the other class, suggesting a biological role for the two classes. Known processing site sequences proved useful in a motif searching strategy to identify processing sites in retroviral protein sequences, particularly in Gag. In all known cleavage sites, the P1 amino acid is hydrophobic and unbranched at the beta-carbon. The sequence requirements of the P1 position were tested by site-directed mutagenesis of the P1 Phe codon in an HIV-1 Pol cleavage site. Mutations were tested for protease-mediated cleavage of the Pol precursor expressed in Escherichia coli.     

Studies on the specificity of HIV protease: An application of Markov chain theory

A sequence-coupled (Markov chain) model is proposed to predict the cleavage sites in proteins by proteases with extended specificity subsites. In addition to the probability of an amino acid occurring at each of these subsites as observed from a training set of oligopeptides known cleavable by HIV protease, the conditional probabilities as reflected by the neighbor-coupled effect along the subsite sequence are also taken into account. These conditional probabilities are derived from an expanded training set consisting of sufficiently large peptide sequences generated by the Monte Carlo sampling process. Very high accuracy was obtained in predicting protein cleavage sites by both HIV-1 and HIV-2 proteases. The new method provides a rapid and accurate means for analyzing the specificity of HIV protease, and hence can be used to help find effective inhibitors of HIV protease as potential drugs against AIDS. The principle of this method can also be used to study the specificity of any multisubsite enzyme.     

Substrate shape determines specificity of recognition for HIV-1 protease: analysis of crystal structures of six substrate complexes

The homodimeric HIV-1 protease is the target of some of the most effective antiviral AIDS therapy, as it facilitates viral maturation by cleaving ten asymmetric and nonhomologous sequences in the Gag and Pol polyproteins. Since the specificity of this enzyme is not easily determined from the sequences of these cleavage sites alone, we solved the crystal structures of complexes of an inactive variant (D25N) of HIV-1 protease with six peptides that correspond to the natural substrate cleavage sites. When the protease binds to its substrate and buries nearly 1000 A2 of surface area, the symmetry of the protease is broken, yet most internal hydrogen bonds and waters are conserved. However, no substrate side chain hydrogen bond is conserved. Specificity of HIV-1 protease appears to be determined by an asymmetric shape rather than a particular amino acid sequence.     

The three-dimensional structure of the aspartyl protease from the HIV-1 isolate BRU.

The crystal structure of the aspartyl protease encoded by the gene pol of the human immunodeficiency virus (HIV-1, isolate BRU) has been determined to 2.7 A resolution. The enzyme, expressed as an insoluble denatured polypeptide in inclusion bodies of Escherichia coli has been renatured and crystallized. It differs by several amino acid replacements from the homologous enzymes of other HIV-1 isolates. A superposition of the C alpha-backbone of the BRU protease with that of the SF2 protease gives a roots mean square positional difference of 0.45 A. Thus, neither the denaturation/renaturation process nor the amino acid replacements have a noticeable effect on the three-dimensional structure of the BRU protease or on the detailed conformation of the catalytic site, which is very similar to that of other aspartyl proteases.     

3C-like protease of rabbit hemorrhagic disease virus: identification of cleavage sites in the ORF1 polyprotein and analysis of cleavage specificity.

Rabbit hemorrhagic disease virus, a positive-stranded RNA virus of the family Caliciviridae, encodes a trypsin-like cysteine protease as part of a large polyprotein. Upon expression in Escherichia coli, the protease releases itself from larger precursors by proteolytic cleavages at its N and C termini. Both cleavage sites were determined by N-terminal sequence analysis of the cleavage products. Cleavage at the N terminus of the protease occurred with high efficiency at an EG dipeptide at positions 1108 and 1109. Cleavage at the C terminus of the protease occurred with low efficiency at an ET dipeptide at positions 1251 and 1252. To study the cleavage specificity of the protease, amino acid substitutions were introduced at the P2, P1, and P1' positions at the cleavage site at the N-terminal boundary of the protease. This analysis showed that the amino acid at the P1 position is the most important determinant for substrate recognition. Only glutamic acid, glutamine, and aspartic acid were tolerated at this position. At the P1' position, glycine, serine, and alanine were the preferred substrates of the protease, but a number of amino acids with larger side chains were also tolerated. Substitutions at the P2 position had only little effect on the cleavage efficiency. Cell-free expression of the C-terminal half of the ORF1 polyprotein showed that the protease catalyzes cleavage at the junction of the RNA polymerase and the capsid protein. An EG dipeptide at positions 1767 and 1768 was identified as the putative cleavage site. Our data show that rabbit hemorrhagic disease virus encodes a trypsin-like cysteine protease that is similar to 3C proteases with regard to function and specificity but is more similar to 2A proteases with regard to size.     

Chemical synthesis of a biotinylated derivative of the simian immunodeficiency virus protease. Purification by avidin affinity chromatography and autocatalytic activation.

The protease from simian immunodeficiency virus (SIV) was chemically synthesized by automated solid-phase technology as an NH2-terminally extended derivative, capped with biotin. Biotin-linker-(SIV protease (1-99)): the linker segment, Gly-Gly-Asp-Arg-Gly-Phe-Ala-Ala, corresponds to the amino acid sequence preceding that of the protease in the SIV gag/pol precursor polyprotein. Accordingly, the Ala-Pro bond joining the octapeptide linker to the protease constitutes a site naturally cleaved by the protease during viral maturation. This strategy for synthesis was designed to facilitate purification of the biotinylated protein derivative from a complex mixture of reaction products by avidin/agarose-affinity chromatography and to provide the means for autocatalytic removal of the biotin-linker segment. As anticipated, folding of the full-length construct leads to activation of the enzyme and excision of the desired 99-residue SIV protease (overall yield, approximately). The specificity of the synthetic SIV protease toward a number of well characterized protein substrates was the same as observed for the nearly identical enzyme from human immunodeficiency virus type 2 (HIV-2 protease) and distinct from that of the more disparate HIV-1 protease. The same functional ordering with respect to the human retroviral proteases was reflected in Ki values observed with a number of protease inhibitors. Thus, the folded synthetic SIV protease shows patterns of specificity and susceptibility to inhibition that are in accord with what would be expected based upon its degree of structural similarity to proteases from HIV-1 and HIV-2.     

HIV-1 protease specificity of peptide cleavage is sufficient for processing of gag and pol polyproteins

The mature proteins of retroviruses originate as a result of proteolytic cleavages of polyprotein precursors. Retroviruses encode proteases responsible for several of these processing events, making them potential antiviral drug targets. A 99-amino acid HIV-1 protease, produced by chemical synthesis or by expression in bacteria, is shown here to hydrolyze peptides corresponding to all of the known cleavage sites in the HIV-1 gag and pol polyproteins. It does not hydrolyze peptides corresponding to an env cleavage site or a distantly related retroviral gag cleavage site.     

Molecular and enzymatic characterization of the porcine endogenous retrovirus protease

The protease of the porcine endogenous retrovirus (PERV) subtypes A/B and C was recombinantly expressed in Escherichia coli as proteolytically active enzyme and characterized. The PERV Gag precursor was also recombinantly produced and used as the substrate in an in vitro enzyme assay in parallel with synthetic nonapeptide substrates designed according to cleavage site sequences identified in the PERV Gag precursor. The proteases of all PERV subtypes consist of 127 amino acid residues with an M(r) of 14,000 as revealed by determining the protease N and C termini. The PERV proteases have a high specificity for PERV substrates and do not cleave human immunodeficiency virus (HIV)-specific substrates, nor are they inhibited by specific HIV protease inhibitors. Among the known retroviral proteases, the PERV proteases resemble most closely the protease of the murine leukemia retrovirus.     

Amino acid substitutions in Gag protein at non-cleavage sites are indispensable for the development of a high multitude of HIV-1 resistance against protease inhibitors.

Amino acid substitutions in human immunodeficiency virus type 1 (HIV-1) Gag cleavage sites have been identified in HIV-1 isolated from patients with AIDS failing chemotherapy containing protease inhibitors (PIs). However, a number of highly PI-resistant HIV-1 variants lack cleavage site amino acid substitutions. In this study we identified multiple novel amino acid substitutions including L75R, H219Q, V390D/V390A, R409K, and E468K in the Gag protein at non-cleavage sites in common among HIV-1 variants selected against the following four PIs: amprenavir, JE-2147, KNI-272, and UIC-94003. Analyses of replication profiles of various mutant clones including competitive HIV-1 replication assays demonstrated that these mutations were indispensable for HIV-1 replication in the presence of PIs. When some of these mutations were reverted to wild type amino acids, such HIV-1 clones failed to replicate. However, virtually the same Gag cleavage pattern was seen, indicating that the mutations affected Gag protein functions but not their cleavage sensitivity to protease. These data strongly suggest that non-cleavage site amino acid substitutions in the Gag protein recover the reduced replicative fitness of HIV-1 caused by mutations in the viral protease and may open a new avenue for designing PIs that resist the emergence of PI-resistant HIV-1.     

Amino acid preferences for a critical substrate binding subsite of retroviral proteases in type 1 cleavage sites

The specificities of the proteases of 11 retroviruses representing each of the seven genera of the family Retroviridae were studied using a series of oligopeptides with amino acid substitutions in the P2 position of a naturally occurring type 1 cleavage site (Val-Ser-Gln-Asn-Tyr Pro-Ile-Val-Gln; the arrow indicates the site of cleavage) in human immunodeficiency virus type 1 (HIV-1). This position was previously found to be one of the most critical in determining the substrate specificity differences of retroviral proteases. Specificities at this position were compared for HIV-1, HIV-2, equine infectious anemia virus, avian myeloblastosis virus, Mason-Pfizer monkey virus, mouse mammary tumor virus, Moloney murine leukemia virus, human T-cell leukemia virus type 1, bovine leukemia virus, human foamy virus, and walleye dermal sarcoma virus proteases. Three types of P2 preferences were observed: a subgroup of proteases preferred small hydrophobic side chains (Ala and Cys), and another subgroup preferred large hydrophobic residues (Ile and Leu), while the protease of HIV-1 preferred an Asn residue. The specificity distinctions among the proteases correlated well with the phylogenetic tree of retroviruses prepared solely based on the protease sequences. Molecular models for all of the proteases studied were built, and they were used to interpret the results. While size complementarities appear to be the main specificity-determining features of the S2 subsite of retroviral proteases, electrostatic contributions may play a role only in the case of HIV proteases. In most cases the P2 residues of naturally occurring type 1 cleavage site sequences of the studied proteases agreed well with the observed P2 preferences.     

Mutational analysis of human immunodeficiency virus type 1 protease suggests functional homology with aspartic proteinases.

Processing of the retroviral gag and pol gene products is mediated by a viral protease. Bacterial expression systems have been developed which permit genetic analysis of the human immunodeficiency virus type 1 protease as measured by cleavage of the pol protein precursor. Deletion analysis of the pol reading frame locates the sequences required to encode a protein with appropriate proteolytic activity near the left end of the pol reading frame but largely outside the gag-pol overlap region, which is at the extreme left end of pol. Most missense mutations within an 11-amino-acid domain highly conserved among retroviral proteases and with sequence similarity to the active site of aspartic proteinases abolish appropriate processing, suggesting that the retrovirus proteases share a catalytic mechanism with aspartic proteinases. Substitution of the amino acids flanking the scissile bond at three of the processing sites encoded by pol demonstrates distinct sequence requirements for cleavage at these different sites. The inclusion of a charged amino acid at the processing site blocks cleavage. A subset of these substitutions also inhibits processing at the nonmutated sites.     

A radiometric assay for HIV-1 protease

A rapid, high-throughput radiometric assay for HIV-1 protease has been developed using ion-exchange chromatography performed in 96-well filtration plates. The assay monitors the activity of the HIV-1 protease on the radiolabeled form of a heptapeptide substrate, [tyrosyl-3,5-3H]Ac-Ser-Gln-Asn-Tyr-Pro-Val-Val-NH2, which is based on the p17-p24 cleavage site found in the viral polyprotein substrate Pr55gag. Specific cleavage of this uncharged heptapeptide substrate by HIV-1 protease releases the anionic product [tyrosyl-3,5-3H]Ac-Ser-Gln-Asn-Tyr, which is retained upon minicolumns of the anion-exchange resin AG1-X8. Protease activity is determined from the recovery of this radiolabeled product following elution with formic acid. This facile and highly sensitive assay may be utilized for steady-state kinetic analysis of the protease, for measurements of enzyme activity during its purification, and as a routine assay for the evaluation of protease inhibitors from natural product or synthetic sources.     

Peptide substrates and inhibitors of the HIV-1 protease

Oligopeptides containing the consensus retroviral protease cleavage sequence Ser/Thr-X-Y-Tyr/Phe-Pro are substrates for purified recombinant HIV-1 protease with Km's in the millimolar range. The minimum sequence containing the consensus pentapeptide which serves as a good substrate is a heptapeptide spanning the P4-P3' residues. Substitution of reduced Phe-Pro or Tyr-Pro dipeptide isosteres or the statine analog 3-hydroxy-4-amino-5-phenylpentanoic acid for the scissile dipeptide afforded inhibitors of HIV-1 protease with Ki values in the micromolar range, three orders of magnitude better in affinity than the corresponding substrates. Inhibitors of HIV-1 protease may provide a novel and potentially useful therapeutic approach to the treatment of acquired immune deficiency syndrome (AIDS).     

Comparison of inhibitor binding in HIV-1 protease and in non-viral aspartic proteases: the role of the flap

The crystal structure of HIV-1 protease with an inhibitor has been compared with the structures of non-viral aspartic proteases complexed with inhibitors. In the dimeric HIV-1 protease, two 4-stranded beta-sheets are formed by half of the inhibitor, residues 27-29, and the flap from each monomer. In the monomeric non-viral enzyme the single flap does not form a beta-sheet with an inhibitor. The HIV-1 protease shows more interactions with a longer peptide inhibitor than are observed in non-viral aspartic protease-inhibitor complexes. This, and the large movement of the flaps, restricts the conformation of the protease cleavage sites in the retroviral polyprotein precursor.