Comparison of the substrate specificity of the human T-cell leukemia virus and human immunodeficiency virus proteinases.

Human T-cell leukemia virus type-1 (HTLV-1) is associated with a number of human diseases. Based on the therapeutic success of human immunodeficiency virus type 1 (HIV-1) PR inhibitors, the proteinase (PR) of HTLV-1 is a potential target for chemotherapy. To facilitate the design of potent inhibitors, the subsite specificity of HTLV-1 PR was characterized and compared to that of HIV-1 PR. Two sets of substrates were used that contained single amino-acid substitutions in peptides representing naturally occurring cleavage sites in HIV-1 and HTLV-1. The original HIV-1 matrix/capsid cleavage site substrate and most of its substituted peptides were not hydrolyzed by the HTLV-1 enzyme, except for those with hydrophobic residues at the P4 and P2 positions. On the other hand, most of the peptides representing the HTLV-1 capsid/nucleocapsid cleavage site were substrates of both enzymes. A large difference in the specificity of HTLV-1 and HIV-1 proteinases was demonstrated by kinetic measurements, particularly with regard to the S4 and S2 subsites, whereas the S1 subsite appeared to be more conserved. A molecular model of the HTLV-1 PR in complex with this substrate was built, based on the crystal structure of the S9 mutant of Rous sarcoma virus PR, in order to understand the molecular basis of the enzyme specificity. Based on the kinetics of shortened analogs of the HTLV-1 substrate and on analysis of the modeled complex of HTLV-1 PR with substrate, the substrate binding site of the HTLV-1 PR appeared to be more extended than that of HIV-1 PR. Kinetic results also suggested that the cleavage site between the capsid and nucleocapsid protein of HTLV-1 is evolutionarily optimized for rapid hydrolysis.     

Solid phase synthesis of the proteinase of bovine leukemia virus. Comparison of its specificity to that of HIV-2 proteinase.

The 126-residue proteinase (PR) of bovine leukemia virus (BLV) was synthesized by solid-phase peptide synthesis and its activity was shown using various oligopeptide substrates representing cleavage sites in BLV, human T-cell leukemia virus type 1 (HTLV-1), murine leukemia virus (MuLV) and human immunodeficiency virus type 1 (HIV-1). The specificity of the BLV PR was also compared to that of chemically synthesized human immunodeficiency virus type 2 (HIV-2) PR. Many of the peptides were cleaved at the expected site, however, 6 out of 15 were hydrolyzed only by one of the PRs. Furthermore, one BLV peptide was processed differently by the two enzymes. These results, together with the relative activities and the lack of inhibition of BLV PR by two HIV-1 PR inhibitors, suggest that the BLV PR specificity is substantially different from that of HIV PRs.     

Narrow substrate specificity and sensitivity toward ligand-binding site mutations of human T-cell Leukemia virus type 1 protease

Human T-cell leukemia virus type 1 (HTLV-1) is associated with a number of human diseases; therefore, its protease is a potential target for chemotherapy. To compare the specificity of HTLV-1 protease with that of human immunodeficiency virus type 1 (HIV-1) protease, oligopeptides representing naturally occurring cleavage sites in various retroviruses were tested. The number of hydrolyzed peptides as well as the specificity constants suggested a substantially broader specificity of the HIV protease. Amino acid residues of HTLV-1 protease substrate-binding sites were replaced by equivalent ones of HIV-1 protease. Most of the single and multiple mutants had altered specificity and a dramatically reduced folding and catalytic capability, suggesting that mutations are not well tolerated in HTLV-1 protease. The catalytically most efficient mutant was that with the flap residues of HIV-1 protease. The inhibition profile of the mutants was also determined for five inhibitors used in clinical practice and inhibitor analogs of HTLV-1 cleavage sites. Except for indinavir, the HIV-1 protease inhibitors did not inhibit wild type and most of the mutant HTLV-1 proteases. The wild type HTLV-1 protease was inhibited by the reduced peptide bond-containing substrate analogs, whereas the mutants showed various degrees of weakened binding capability. Most interesting, the enzyme with HIV-1-like residues in the flap region was the most sensitive to the HIV-1 protease inhibitors and least sensitive to the HTLV-1 protease inhibitors, indicating that the flap plays an important role in defining the specificity differences of retroviral proteases.     

Effect of serine and tyrosine phosphorylation on retroviral proteinase substrates.

Vimentin, a cellular substrate of HIV type 1 (HIV-1) proteinase, contains a protein kinase C (PKC) phosphorylation site at one of its cleavage sites. Peptides representing this site were synthesized in P2 Ser-phosphorylated and nonphosphorylated forms. While the nonphosphorylated peptide was a fairly good substrate of the enzyme, phosphorylation prevented hydrolysis. Phosphorylation of human recombinant vimentin by PKC prevented its processing within the head domain, where the phosphorylation occurred. Oligopeptides representing naturally occurring cleavage sites at the C-terminus of the Rous sarcoma virus integrase were assayed as substrates of the avian proteinase. Unlike the nonphosphorylated peptides, a Ser-phosphorylated peptide was not hydrolyzed by the enzyme at the Ser-Pro bond, suggesting the role of previously established phosphorylation in processing at this site. Ser-phosphorylated and Tyr-phosphorylated forms of model substrates were also tested as substrates of the HIV-1 and the avian retroviral proteinases. In contrast to the moderate effect of P4 Ser phosphorylation, phosphorylation of P1 Tyr prevented substrate hydrolysis by HIV-1 proteinase. Substrate phosphorylation had substantially smaller effects on the hydrolysis by the avian retroviral proteinase. As the active retroviral proteinase as well as various protein kinases are incorporated into mature virions, substrate phosphorylation resulting in attenuation or prevention of proteolytic processing may have important consequences in the regulation of the retroviral life cycle as well as in virus-host cell interactions.     

Expression and characterization of human foamy virus proteinase

The human foamy virus proteinase was expressed in fusion with maltose binding protein in Escherichia coli and purified. The specific activity of the fusion protein was similar to that of the processed enzyme. The kinetic constants on foamy virus cleavage site substrates were very low but comparable to those obtained with the gag-encoded avian proteinase on its own substrates. The proteinase showed preference for high ionic strength and a pH optimum of 6.6. None of the tested retroviral cleavage site peptides were substrates, however, some peptides representing cleavage sites in retrotransposons were properly processed by the enzyme.     

Specificity studies on retroviral proteinase from myeloblastosis-associated virus.

The specificity of the p15 proteinase of myeloblastosis-associated virus (MAV) was tested with nonviral high molecular weight substrates and with synthetic peptides. Peptides with sequences spanning known cleavage sites in viral polyproteins of Rous sarcoma virus (RSV) and avian leukemia viruses, as well as in BSA and HSA, were synthesized, and the rate of their cleavage by the MAV proteinase was compared. Synthetic peptides require for successful cleavage at least 4 residues at the N-terminal side and 3 residues at the C-terminal side. The proteinase shows a preference for hydrophobic residues with bulky side chains (Met, Tyr, Phe) in P3, although Arg and Gln can also be accepted. Small hydrophobic residues are required in P2 and P2', and large hydrophobic residues (Tyr, Met, Phe/p-nitro-Phe) are preferred in both P1 and P1'. The difference between the specificity of the p15 proteinase and that of the HIV-1 proteinase mostly pertains to position P2' of the substrate, where bulkier side chains are accepted by the HIV-1 proteinase (Richards et al., 1990). A good chromogenic substrate for the MAV and RSV proteinases was developed and used to further characterize the MAV proteinase activity with respect to ionic strength and pH. The activity of the proteinase is strongly dependent on ionic strength and pH. Both the kcat and Km values contribute to a higher cleavage efficiency at higher salt concentrations and show a bell-shaped pH dependence curve with a sharp maximum at pH 5.5 (kcat) and 6.5 (Km).     

Cloning of the bovine leukemia virus proteinase in Escherichia coli and comparison of its specificity to that of human T-cell leukemia virus proteinase

The proteinase of bovine leukemia virus (BLV) was cloned into pMal-c2 vector with N-terminal or with N- as well as C-terminal flanking sequences, and expressed in fusion with maltose binding protein. The proteinase self-processed itself from the fusion protein during expression and formed inclusion bodies. The enzyme was purified from inclusion bodies by cation-exchange chromatography followed by gel filtration. Specificity of the enzyme was compared to that of human T-cell leukemia proteinase type 1. Although the two viruses belong to the same subfamily of retroviruses, the differences in their proteinase specificity, based on kinetics with oligopeptide substrates representing naturally occurring cleavage sites as well as on inhibition pattern, appear to be pronounced.     

Substrate requirements of human rhinovirus 3C protease for peptide cleavage in vitro

A series of synthetic peptides representing authentic proteolytic cleavage sites of human rhinovirus type 14 were assayed as substrates for purified 3C protease. Competition cleavage assays were employed to determine the relative specificity constants (Kcat/Km) for substrates with sequences related to the viral 2C-3A cleavage site. Variable length peptides representing the 2C-3A cleavage site were cleaved with comparable efficiency. These studies defined a minimum substrate of 6 amino acids (TLFQ/GP), although retention of the residue at position P5 (ETLFQ/GP) resulted in a better substrate by an order of magnitude. Amino acid substitutions at position P5, P4, P1', or P2' indicated that the identity of the residue at position P5 was not critical, whereas substitutions at position P4, P1' or P2' resulted in substrates with Kcat/Km values varying over 2 orders of magnitude. In contrast to the 2C-3A cleavage site, small peptide derivatives representative of the 3A-3B cleavage site were relatively poor substrates, which suggested that residues flanking the minimum core sequence may influence susceptibility to cleavage. The 3C protease of rhinovirus type 14 was also capable of cleaving peptides representing comparable cleavage sites predicted for coxsackie B virus and poliovirus.     

Functional identification of nucleotides conferring substrate specificity to retroviral integrase reactions.

The long terminal repeats (LTRs) that flank the retroviral DNA genome play a distinct role in the integration process by acting as specific substrates for the integrase (IN). The role of LTR sequences in providing substrate recognition and specificity to integration reactions was investigated for INs from human immunodeficiency virus type 1 (HIV-1), Moloney murine leukemia virus (M-MuLV), human T-cell leukemia virus type 1 (HTLV-1), and human T-cell leukemia virus type 2 (HTLV-2). Overall, these INs required specific LTR sequences for optimal catalysis of 3'-processing reactions, as opposed to strand transfer and disintegration reactions. It is of particular note that in strand transfer reactions the sites of integration were similar among the four INs. In the 3'-processing reaction, sequence specificity for each IN was traced to the three nucleotides proximal to the conserved CA. Reactions catalyzed by M-MuLV IN were additionally influenced by upstream regions. The nucleotide requirements for optimal catalysis differed for each IN. HIV-1 IN showed a broad range of substrate specificities, while HTLV-1 IN and HTLV-2 IN had more defined sequence requirements. M-MuLV IN exhibited greater activity with the heterologous LTR substrates than with its own wild-type substrate. This finding was further substantiated by the high levels of activity catalyzed by the IN on modified M-MuLV LTRs. This work suggests that unlike the other INs examined, M-MuLV IN has evolved with an IN-LTR interaction that is suboptimal.     

Effect of sequence polymorphism and drug resistance on two HIV-1 Gag processing sites.

The HIV-1 proteinase (PR) has proved to be a good target for antiretroviral therapy of AIDS, and various PR inhibitors are now in clinical use. However, there is a rapid selection of viral variants bearing mutations in the proteinase that are resistant to clinical inhibitors. Drug resistance also involves mutations of the nucleocapsid/p1 and p1/p6 cleavage sites of Gag, both in vitro and in vivo. Cleavages at these sites have been shown to be rate limiting steps for polyprotein processing and viral maturation. Furthermore, these sites show significant sequence polymorphism, which also may have an impact on virion infectivity. We have studied the hydrolysis of oligopeptides representing these cleavage sites with representative mutations found as natural variations or that arise as resistant mutations. Wild-type and five drug resistant PRs with mutations within or outside the substrate binding site were tested. While the natural variations showed either increased or decreased susceptibility of peptides toward the proteinases, the resistant mutations always had a beneficial effect on catalytic efficiency. Comparison of the specificity changes obtained for the various substrates suggested that the maximization of the van der Waals contacts between substrate and PR is the major determinant of specificity: the same effect is crucial for inhibitor potency. The natural nucleocapsid/p1 and p1/p6 sites do not appear to be optimized for rapid hydrolysis. Hence, mutation of these rate limiting cleavage sites can partly compensate for the reduced catalytic activity of drug resistant mutant HIV-1 proteinases.     

Effect of substrate residues on the P2' preference of retroviral proteinases.

The substrate sequence requirements for preference toward P2' Glu residue by human immunodeficiency virus type 1 (HIV-1) proteinase were studied in both the matrix protein/ capsid protein (MA/CA) and CA/p2 cleavage site sequence contexts. These sequences represent typical type 1 (-aromatic*Pro-) and type 2 (-hydrophobic* hydrophobic-) cleavage site sequences, respectively. While in the type 1 sequence context, the preference for P2' Glu over Ile or Gln was found to be strongly dependent on the ionic strength and the residues being outside the P2-P2' region of the substrate, it remained preferable in the type 2 substrates when typical type 1 substrate sequence residues were substituted into the outside regions. The pH profile of the specificity constants suggested a lower pH optimum for substrates having P2' Glu in contrast to those having uncharged residues, in both sequence contexts. The very low frequency of P2' Glu in naturally occurring retroviral cleavage sites of various retroviruses including equine infectious anemia virus (EIAV) and murine leukemia virus (MuLV) suggests that such a residue may not have a general regulatory role in the retroviral life cycle. In fact, unlike HIV-1 and HIV-2, EIAV and MuLV proteinases do not favor P2' Glu in either the MA/CA or CA/p2 sequence contexts.     

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.     

In vitro processing of HIV-1 nucleocapsid protein by the viral proteinase: effects of amino acid substitutions at the scissile bond in the proximal zinc finger sequence

The human immunodeficiency virus type 1 (HIV-1) nucleocapsid protein flanked by Gag sequences (r-preNC) was expressed in Escherichia coli and purified. HIV-1 proteinase cleaved r-preNC to the "mature" NCp7 form, which is comprised of 55 residues. Further incubation resulted in cleavages of NCp7 itself between Phe16 and Asn17 of the proximal zinc finger domain and between Cys49 and Thr50 in the C-terminal part. Kinetic parameters determined for the cleavage of oligopeptides corresponding to the cleavage sites in r-preNC correlated well with the sequential processing of r-preNC. Mutations of Asn17 were introduced to alter the susceptibility of NC protein to HIV-1 proteinase. While mutating Asn17 to Ala resulted in a protein which was processed in a manner similar to that of the wild type, mutating it to Phe or Leu resulted in proteins which were processed at a substantially higher rate at this site than the wild type. Mutation of Asn17 to Lys or Gly resulted in proteins which were very poorly cleaved at this site. Oligopeptides containing the same amino acid substitutions at the cleavage site of the proximal zinc finger domain were also tested as substrates of the proteinase, and the kinetic parameters agreed well with the semiquantitative results obtained with the protein substrates.     

Specificity and inhibition of proteases from human immunodeficiency viruses 1 and 2

Highly purified, recombinant preparations of the virally encoded proteases from human immunodeficiency viruses (HIV) 1 and 2 have been compared relative to 1) their specificities toward non-viral protein and synthetic peptide substrates, and 2) their inhibition by several P1-P1' pseudodipeptidyl-modified substrate analogs. Hydrolysis of the Leu-Leu and Leu-Ala bonds in the Pseudomonas exotoxin derivative, Lys-PE40, is qualitatively the same for HIV-2 protease as published earlier for the HIV-1 enzyme (Tomasselli, A. G., Hui, J. O., Sawyer, T. K., Staples, D. J., FitzGerald, D. J., Chaudhary, V. K., Pastan, I., and Heinrikson, R. L. (1990) J. Biol. Chem. 265, 408-413). However, the rates of cleavage at these two sites are reversed for the HIV-2 protease which prefers the Leu-Ala bond. The kinetics of hydrolysis of this protein substrate by both enzymes are mirrored by those obtained from cleavage of model peptides. Hydrolysis by the two proteases of other synthetic peptides modeled after processing sites in HIV-1 and HIV-2 gag polyproteins and selected analogs thereof demonstrated differences, as well as similarities, in selectivity. For example, while the two proteases were nearly identical in their rates of cleavage of the Tyr-Pro bond in the HIV-1 gag fragment, Val-Ser-Gln-Asn-Tyr-Pro-Ile-Val, the HIV-1 protease showed a 64-fold enhancement over the HIV-2 enzyme in hydrolysis of a Tyr-Val bond in the same template. Accordingly, the HIV-2 protease appears to have a different specificity than the HIV-1 enzyme; it is better able to hydrolyze substrates with small amino acids in P1 and P1', but is variable in its rate of hydrolysis of peptides with bulky substituents in these positions. In addition to these comparisons of the two proteases with respect to substrate specificity, we present inhibitor structure-activity data for the HIV-2 protease. Relative to P1-P1' statine or Phe psi [CH2N]Pro-modified pseudopeptidyl inhibitors, compounds having Xaa psi[CH(OH)CH2]Yaa inserts were found to show significantly higher affinities to both enzymes, generally binding from 10 to 100 times stronger to HIV-1 protease than to the HIV-2 enzyme. Molecular modeling comparisons based upon the sequence homology of the two enzymes and x-ray crystal structures of HIV-1 protease suggest that most of the nonconservative amino acid replacements occur in regions well outside the catalytic cleft, while only subtle structural differences exist within the active site.(ABSTRACT TRUNCATED AT 400 WORDS)     

A catalogue of putative HIV-1 protease host cell substrates

Abstract Processing of human immunodeficiency virus (HIV) proteins by the HIV-1 protease is essential for HIV infectivity. In addition, several studies have revealed cleavage of human proteins by this viral protease during infection; however, no large-scale HIV-1 protease degradomics study has yet been performed. To identify putative host substrates in an unbiased manner and on a proteome-wide scale, we used positional proteomics to identify peptides reporting protein processing by the HIV-1 protease, and a catalogue of over 120 cellular HIV-1 protease substrates processed in vitro was generated. This catalogue includes previously reported substrates as well as recently described interaction partners of HIV-1 proteins. Cleavage site alignments revealed a specificity profile in good correlation with previous studies, even though the ELLE consensus motif was not cleaved efficiently when incorporated into peptide substrates due to subsite cooperativity. Our results are further discussed in the context of HIV-1 infection and the complex substrate recognition by the viral protease.     

Molecular basis for the relative substrate specificity of human immunodeficiency virus type 1 and feline immunodeficiency virus proteases

We have used a random hexamer phage library to delineate similarities and differences between the substrate specificities of human immunodeficiency virus type 1 (HIV-1) and feline immunodeficiency virus (FIV) proteases (PRs). Peptide sequences were identified that were specifically cleaved by each protease, as well as sequences cleaved equally well by both enzymes. Based on amino acid distinctions within the P3-P3' region of substrates that appeared to correlate with these cleavage specificities, we prepared a series of synthetic peptides within the framework of a peptide sequence cleaved with essentially the same efficiency by both HIV-1 and FIV PRs, Ac-KSGVF/VVNGLVK-NH(2) (arrow denotes cleavage site). We used the resultant peptide set to assess the influence of specific amino acid substitutions on the cleavage characteristics of the two proteases. The findings show that when Asn is substituted for Val at the P2 position, HIV-1 PR cleaves the substrate at a much greater rate than does FIV PR. Likewise, Glu or Gln substituted for Val at the P2' position also yields peptides specifically susceptible to HIV-1 PR. In contrast, when Ser is substituted for Val at P1', FIV PR cleaves the substrate at a much higher rate than does HIV-1 PR. In addition, Asn or Gln at the P1 position, in combination with an appropriate P3 amino acid, Arg, also strongly favors cleavage by FIV PR over HIV PR. Structural analysis identified several protease residues likely to dictate the observed specificity differences. Interestingly, HIV PR Asp30 (Ile-35 in FIV PR), which influences specificity at the S2 and S2' subsites, and HIV-1 PR Pro-81 and Val-82 (Ile-98 and Gln-99 in FIV PR), which influence specificity at the S1 and S1' subsites, are residues which are often involved in development of drug resistance in HIV-1 protease. The peptide substrate KSGVF/VVNGK, cleaved by both PRs, was used as a template for the design of a reduced amide inhibitor, Ac-GSGVF Psi(CH(2)NH)VVNGL-NH(2.) This compound inhibited both FIV and HIV-1 PRs with approximately equal efficiency. These findings establish a molecular basis for distinctions in substrate specificity between human and feline lentivirus PRs and offer a framework for development of efficient broad-based inhibitors.     

Specificity assay of serine proteinases by reverse-phase high-performance liquid chromatography analysis of competing oligopeptide substrate library

In this paper we present an HPLC method developed for quick activity and specificity analysis of serine proteinases. The method applies a carefully designed peptide library in which the individual components differ only at the potential cleavage site for enzymes. The library has seven members representing seven different cleavage sites and it offers substrates for both trypsin and chymotrypsin-like enzymes. The individual peptide substrates compete for the proteinase during the enzymatic reaction. The reaction is monitored by RP-HPLC separation of the components. We describe the systematic design of the competitive peptide substrate library and the test of the system with eight different serine proteinases. The specificity profiles of the investigated enzymes as determined by the new method were essentially identical to the ones reported in the literature, verifying the ability of the system to characterize substrate specificity. The tests also demonstrated that the system could detect even subtle specificity differences of two isoforms of an enzyme. In addition to recording qualitative specificity profiles, data provided by the system can be analyzed quantitatively, yielding specificity constant values. This method can be a useful tool for quick analysis of uncharacterized gene products as well as new forms of enzymes generated by protein engineering.     

Re-mapping the molecular features of the human immunodeficiency virus type 1 and human T-cell lymphotropic virus type 1 Brazilian sequences using a bioinformatics unit established in Salvador, Bahia, Brazil, to give support to the viral epidemiology studies

The analysis of genetic data for human immunodeficiency virus type 1 (HIV-1) and human T-cell lymphotropic virus type 1 (HTLV-1) is essential to improve treatment and public health strategies as well as to select strains for vaccine programs. However, the analysis of large quantities of genetic data requires collaborative efforts in bioinformatics, computer biology, molecular biology, evolution, and medical science. The objective of this study was to review and improve the molecular epidemiology of HIV-1 and HTLV-1 viruses isolated in Brazil using bioinformatic tools available in the Laboratório Avan?ado de Sáude Pública (Lasp) bioinformatics unit. The analysis of HIV-1 isolates confirmed a heterogeneous distribution of the viral genotypes circulating in the country. The Brazilian HIV-1 epidemic is characterized by the presence of multiple subtypes (B, F1, C) and B/F1 recombinant virus while, on the other hand, most of the HTLV-1 sequences were classified as Transcontinental subgroup of the Cosmopolitan subtype. Despite the high variation among HIV-1 subtypes, protein glycosylation and phosphorylation domains were conserved in the pol, gag, and env genes of the Brazilian HIV-1 strains suggesting constraints in the HIV-1 evolution process. As expected, the functional protein sites were highly conservative in the HTLV-1 env gene sequences. Furthermore, the presence of these functional sites in HIV-1 and HTLV-1 strains could help in the development of vaccines that pre-empt the viral escape process.     

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.