Pentacycloundecane-based inhibitors of wild-type C-South African HIV-protease

In this study, we present the first account of pentacycloundecane (PCU) peptide based HIV-protease inhibitors. The inhibitor exhibiting the highest activity made use of a natural HIV-protease substrate peptide sequence, that is, attached to the cage (PCU-EAIS). This compound showed nanomolar IC₅₀ activity against the resistance-prone wild type C-South African HIV-protease (C-SA) catalytic site via a norstatine type functional group of the PCU hydroxy lactam. NMR was employed to determine a logical correlation between the inhibitory concentration (IC₅₀) results and the 3D structure of the corresponding inhibitors in solution. NMR investigations indicated that the activity is related to the chirality of the PCU moiety and its ability to induce conformations of the coupled peptide side chain. The results from docking experiments coincided with the experimental observed activities. These findings open up useful applications for this family of cage peptide inhibitors, considering the vast number of alternative disease related proteases that exist.     

Characterization of an active single polypeptide form of the human immunodeficiency virus type 1 protease

The pepsin-like aspartyl proteases consist of a single polypeptide chain with topologically similar amino- and carboxyl-terminal domains, each of which contributes 1 aspartic acid residue to the active site. This structure has been proposed to have evolved by gene duplication and fusion from a dimeric enzyme composed of two identical polypeptide chains, such as the aspartyl protease (PRT) of human immunodeficiency virus type 1 (HIV-1). To determine if a single polypeptide form of the HIV-1 protease would be enzymatically active, two protease coding regions were linked to form a dimeric gene (pFGGP). Expression of this gene in Escherichia coli yielded a protein with the expected molecular mass of 22 kDa. The in vitro kinetic parameters of PRT and FGGP (where FGGP is the single polypeptide form of the HIV-1 protease with 2 glycine residues connecting the two subunits) for three peptide substrates are similar. Construction and analysis of a CheY-GAG-FGGP fusion protein demonstrated that FGGP is capable of precursor processing in vivo. Mutation of one or both of the active site aspartates to either asparagine or glutamate rendered the enzyme inactive, demonstrating that both active site aspartate residues are required for enzymatic activity.     

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.     

F99 is Critical for Dimerization and Activation of South African HIV-1 Subtype C Protease

HIV-1 protease (PR) is an obligate homodimer which plays a pivotal role in the maturation and hence propagation of HIV. Although successful developments on PR active site inhibitors have been achieved, the major limiting factor has been the emergence of HIV drug-resistant strains. Disruption of the dimer interface serves as an alternative mechanism to inactivate the enzyme. The terminal residue, F99, was mutated to an alanine to investigate its contribution to dimer stability in the South African HIV-1 subtype C (C-SA) PR. The F99A PR and wild-type C-SA PR were overexpressed and purified. The activities of the PRs and their ability to bind an active site inhibitor, acetyl-pepstatin, were determined in vitro. The F99A PR showed no activity and the inability to bind to the inhibitor. Secondary and quaternary structure analysis were performed and revealed that the F99A PR is monomeric with reduced β-sheet content. The mutation of F99 to alanine disrupted the presumed ‘lock-and-key’ motif at the terminal dimer interface, in turn creating a cavity at the N- and C-terminal antiparallel β-sheet. These findings support the design of inhibitors targeting the C-terminus of the C-SA PR, centered on interactions with the bulky F99.     

Pentacycloundecane lactam vs lactone norstatine type protease HIV inhibitors: binding energy calculations and DFT study

Background

Novel pentacycloundecane (PCU)-lactone-CO-EAIS peptide inhibitors were designed, synthesized, and evaluated against wild-type C-South African (C-SA) HIV-1 protease. Three compounds are reported herein, two of which displayed IC₅₀ values of less than 1.00 μM. A comparative MM-PB(GB)SA binding free energy of solvation values of PCU-lactam and lactone models and their enantiomers as well as the PCU-lactam-NH-EAIS and lactone-CO-EAIS peptide inhibitors and their corresponding diastereomers complexed with South African HIV protease (C-SA) was performed. This will enable us to rationalize the considerable difference between inhibitory concentration (IC₅₀) of PCU-lactam-NH-EAIS and PCU-lactone-CO-EAIS peptides.

Results

The PCU-lactam model exhibited more negative calculated binding free energies of solvation than the PCU-lactone model. The same trend was observed for the PCU-peptide inhibitors, which correspond to the experimental activities for the PCU-lactam-NH-EAIS peptide (IC₅₀ = 0.076 μM) and the PCU-lactone-CO-EAIS peptide inhibitors (IC₅₀ = 0.850 μM). Furthermore, a density functional theory (DFT) study on the natural atomic charges of the nitrogen and oxygen atoms of the three PCU-lactam, PCU-lactim and PCU-lactone models were performed using natural bond orbital (NBO) analysis. Electrostatic potential maps were also used to visualize the electron density around electron-rich regions. The asymmetry parameter (η) and quadrupole coupling constant (χ) values of the nitrogen and oxygen nuclei of the model compounds were calculated at the same level of theory. Electronic molecular properties including polarizability and electric dipole moments were also calculated and compared. The Gibbs theoretical free solvation energies of solvation (∆G_solv) were also considered.

Conclusions

A general trend is observed that the lactam species appears to have a larger negative charge distribution around the heteroatoms, larger quadrupole constant, dipole moment and better solvation energy, in comparison to the PCU-lactone model. It can be argued that these characteristics will ensure better eletronic interaction between the lactam and the receptor, corresponding to the observed HIV protease activities in terms of experimental IC₅₀ data.

Electronic supplementary material

The online version of this article (doi:10.1186/s12929-015-0115-5) contains supplementary material, which is available to authorized users.     

Studies on flexibility and binding affinity of Asp25 of HIV-1 protease mutants

We have investigated and highlighted the behavior of binding residue, Asp25 by computational analysis, which play an important role in understanding docking process with drug molecule, Ritonavir (Norvir^®) and the flexibility nature of the Human Immunodeficiency Virus-1 (HIV-1) protease enzyme. It is well known that Ritonavir is a potent and a selective HIV-1 protease inhibitor. Molecular dockings were performed in order to gain insights regarding the binding mode of this inhibitor. In our analysis, we observed Ritonavir had different rank orders of scores against different mutant of this enzyme. Asp25 of the enzyme was found to be the active site for all the mutants. The results clearly suggest that Ritonavir is not able to appropriately bind at the active site of each HIV-1 protease mutant due to RMSD difference of the amino acid (Asp) at the position 25 of all mutants. These findings support the concept that 3D space of active site is a qualitative assessment for binding affinity of inhibitor with an enzyme. The investigation on the flexibility nature of Asp25 by normal mode analysis, show that binding residue posses less flexibility due to its solvation potential. The overall analysis of our study brings clarity to the binding behavior with respect to the different mutants with Ritonavir on the basis RMSD and also on the flexible nature of HIV-1 protease enzyme with respect to Asp25 position.     

Kinetic studies of human immunodeficiency virus type 1 protease and its active-site hydrogen bond mutant A28S.

Human immunodeficiency virus type 1 (HIV-1) protease optimally catalyzes in the pH range of 4-6 in contrast to nearly all of the other eukaryotic aspartic proteases, which catalyze best in the pH range of 2-4. A possible structural reason for the higher optimal pH of HIV-1 protease is the absence of a hydrogen bond to the carboxyl group of active-site Asp25, which is nearly universally present in others. To investigate this hypothesis, we have mutated residue 28 in HIV-1 protease from alanine to serine. Both the wild-type and the mutant A28S enzymes have been overexpressed in Escherichia coli using a chemically synthesized gene and purified for a comparative study in enzyme kinetics. The kcat and Km values were determined by a radiometric assay for the wild-type enzyme from pH 3.2 to 7.0, and for the mutant enzyme from pH 3.2 to 6.0. The low pK values of the active site of the free enzyme, pKe1, are 3.3 and 3.4 for the wild-type and mutant enzymes, respectively. The low pK values of the active site of the enzyme bound to substrate, pKes1, are 5.1 and 4.3 for the wild-type and mutant enzymes, respectively. The high pK values of the free enzyme, pKe2, are 6.8 and 5.6, and the corresponding ones for the substrate-bound enzyme, pKes2, are 6.9 and 6.0 for the wild-type and mutant enzymes, respectively. The lowering of pK values in mutant HIV-1 protease indicates that the hydroxyl group of Ser28 forms a new hydrogen bond to active-site Asp25 to increase its acidity.     

Nine Crystal Structures Determine the Substrate Envelope of the MDR HIV-1 Protease

Under drug selection pressure, emerging mutations render HIV-1 protease drug resistant, leading to the therapy failure in anti-HIV treatment. It is known that nine substrate cleavage site peptides bind to wild type (WT) HIV-1 protease in a conserved pattern. However, how the multidrug-resistant (MDR) HIV-1 protease binds to the substrate cleavage site peptides is yet to be determined. MDR769 HIV-1 protease (resistant mutations at residues 10, 36, 46, 54, 62, 63, 71, 82, 84, and 90) was selected for present study to understand the binding to its natural substrates. MDR769 HIV-1 protease was co-crystallized with nine substrate cleavage site hepta-peptides. Crystallographic studies show that MDR769 HIV-1 protease has an expanded substrate envelope with wide open flaps. Furthermore, ligand binding energy calculations indicate weaker binding in MDR769 HIV-1 protease-substrate complexes. These results help in designing the next generation of HIV-1 protease inhibitors by targeting the MDR HIV-1 protease.     

Interaction of ganoderic acid on HIV related target: molecular docking studies

Finding the ultimate HIV cure remain a challenging tasks for decades. Various active compounds have been tested against various components of the virus in the effort to halt the virus development in infected host. The idea of finding cure from known pharmacologically active natural occurring compounds is intriguing and practical. Ganoderma lucidum (Ling-Zhi or Reishi) is one of the most productive and pharmacologically active compounds found in Asian countries. It has been used traditionally for many years throughout different cultures. More than a decade ago, el-Mekkawy and co-workers (1998) have tested several active compounds found in this plant. They have successfully identified several active compounds with reasonable inhibitory activity against HIV protease however; no further studies were done on these compounds. This study aimed to elucidate interactions for one of the active compounds of Ganoderma lucidum namely ganoderic acid with HIV-1 protease using molecular docking simulation. This study revealed four hydrogen bonds formed between model34 of ganoderic acid B and 1HVR. Hydrogen bonds in 1HVR-Model34 complex were formed through ILE50, ILE50', ASP29 and ASP30 residues. Interestingly similar interactions were also observed in the native ligand in 1HVR. Furthermore, interactions involving ILE50 and ILE50' residues have been previously identified to play central roles in HIV-1 protease-ligand interactions.These observed interactions not only suggested HIV-1 protease in general is a suitable target for ganoderic acid B, they also indicated a huge potential for HIV drug discovery based on this compound.     

Defining the level of human immunodeficiency virus type 1 (HIV-1) protease activity required for HIV-1 particle maturation and infectivity.

The human immunodeficiency virus type 1 (HIV-1) protease is the enzyme required for processing of the Gag and Gag-Pol polyproteins to yield mature, infectious virions. Although the complete absence of proteolytic activity prevents maturation, the level of activity sufficient for maturation and subsequent infectivity has not been determined. Amino acid substitutions that reduce catalytic activity without affecting substrate recognition have been engineered into the active site of the HIV-1 protease. The catalytic efficiency (kcat) of the HIV-1 protease is decreased 4-fold when threonine 26 is replaced by serine (T26S) and approximately 50-fold when alanine 28 is replaced by serine (A28S). Genes containing these mutations were cloned into a proviral vector for analysis of their effects on virion maturation and infectivity. The results show that virions containing the T26S protease variant, in which only 25% of the protease is active, are very similar to wild-type virions, although slight reductions in infectivity are observed. Virions containing the A28S protease variant are not infectious, even though a limited amount of polyprotein processing does occur. There appears to be a linear correlation between the level of protease activity and particle infectivity. Our observations suggest that a threshold of protease activity exists between a 4-fold and 50-fold reduction, below which processing is insufficient to yield infectious particles. Our data also suggest that a reduction of protease activity by 50-fold or greater is sufficient to prevent the formation of infectious particles.     

Structural insights into the South African HIV-1 subtype C protease: impact of hinge region dynamics and flap flexibility in drug resistance

The HIV protease plays a major role in the life cycle of the virus and has long been a target in antiviral therapy. Resistance of HIV protease to protease inhibitors (PIs) is problematic for the effective treatment of HIV infection. The South African HIV-1 subtype C protease (C-SA PR), which contains eight polymorphisms relative to the consensus HIV-1 subtype B protease, was expressed in Escherichia coli, purified, and crystallized. The crystal structure of the C-SA PR was resolved at 2.7?Å, which is the first crystal structure of a HIV-1 subtype C protease that predominates in Africa. Structural analyses of the C-SA PR in comparison to HIV-1 subtype B proteases indicated that polymorphisms at position 36 of the homodimeric HIV-1 protease may impact on the stability of the hinge region of the protease, and hence the dynamics of the flap region. Molecular dynamics simulations showed that the flap region of the C-SA PR displays a wider range of movements over time as compared to the subtype B proteases. Reduced stability in the hinge region resulting from the absent E35-R57 salt bridge in the C-SA PR, most likely contributes to the increased flexibility of the flaps which may be associated with reduced susceptibility to PIs.

An animated interactive 3D complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:36     

Inhibition of HIV-1 protease by short peptides derived from the terminal segments of the protease.

The active HIV-1 protease is a homodimeric enzyme. A beta-sheet consisting of N- and C-terminal segments provides the main driving force for dimerization of the inactive protomers. Several short peptides with sequences derived from the N- and C-termini of the protease were tested for inhibition of protease activity and for inhibition of HIV-1 replication in lymphocytes. Medium inhibitory activity was found with each of the peptides in the enzyme test and no inhibition of the lymphocytes was found up to 200 micrograms/ml. The enzyme tests indicate that HIV-1 protease is the target of the inhibitory action. Synergistic action could not be found with pairs of the peptides derived from the two different termini. Prolonged incubation with one of the peptides increased inhibition indicating a slow dissociation of the protease dimers. No cytotoxic effect of the inhibitors could be found below 200 micrograms/ml.     

Crystal structures of multidrug-resistant HIV-1 protease in complex with two potent anti-malarial compounds

Two potent inhibitors (compounds 1 and 2) of malarial aspartyl protease, plasmepsin-II, were evaluated against wild type (NL4-3) and multidrug-resistant clinical isolate 769 (MDR) variants of human immunodeficiency virus type-1 (HIV-1) aspartyl protease. Enzyme inhibition assays showed that both 1 and 2 have better potency against NL4-3 than against MDR protease. Crystal structures of MDR protease in complex with 1 and 2 were solved and analyzed. Crystallographic analysis revealed that the MDR protease exhibits a typical wide-open conformation of the flaps (Gly48 to Gly52) causing an overall expansion in the active site cavity, which, in turn caused unstable binding of the inhibitors. Due to the expansion of the active site cavity, both compounds showed loss of direct contacts with the MDR protease compared to the docking models of NL4-3. Multiple water molecules showed a rich network of hydrogen bonds contributing to the stability of the ligand binding in the distorted binding pockets of the MDR protease in both crystal structures. Docking analysis of 1 and 2 showed a decrease in the binding affinity for both compounds against MDR supporting our structure-function studies. Thus, compounds 1 and 2 show promising inhibitory activity against HIV-1 protease variants and hence are good candidates for further development to enhance their potency against NL4-3 as well as MDR HIV-1 protease variants.     

In Silico Prediction of Mutant HIV-1 Proteases Cleaving a Target Sequence

HIV-1 protease represents an appealing system for directed enzyme re-design, since it has various different endogenous targets, a relatively simple structure and it is well studied. Recently Chaudhury and Gray (Structure (2009) 17: 1636–1648) published a computational algorithm to discern the specificity determining residues of HIV-1 protease. In this paper we present two computational tools aimed at re-designing HIV-1 protease, derived from the algorithm of Chaudhuri and Gray. First, we present an energy-only based methodology to discriminate cleavable and non cleavable peptides for HIV-1 proteases, both wild type and mutant. Secondly, we show an algorithm we developed to predict mutant HIV-1 proteases capable of cleaving a new target substrate peptide, different from the natural targets of HIV-1 protease. The obtained in silico mutant enzymes were analyzed in terms of cleavability and specificity towards the target peptide using the energy-only methodology. We found two mutant proteases as best candidates for specificity and cleavability towards the target sequence.     

Comparative study of some energetic and steric parameters of the wild type and mutants HIV-1 protease: a way to explain the viral resistance

Because, in vivo , the HIV-1 PR (HIV-1 protease) present a high mutation rate we performed a comparative study of the energetic behaviors of the wild type HIV-1 PR and four type of mutants: Val82/Asn; Val82/Asp; Gln7/Lys, Leu33/Ile, Leu63/Ile; Ala71/Thr, Val82/Ala. We suggest that the energetic fluctuation (electrostatic, van der Waals and torsion energy) of the mutants and the solvent accessible surface (SAS) values can be useful to explain the viral resistance process developed by HIV-1 PR. The number and localization of enzyme mutations induce important modifications of the van der Waals and torsional energy, while the electrostatic energy has an insignificant fluctuation. We showed that the viral resistance can be explored if the solvent accessible surfaces of the active site for the mutant structures are calculated. In this paper we have obtained the solvent accessible surface for a group of 15 mutants (11 mutants obtained by Protein Data Bank (PDB) file, 4 mutants modeled by CHARMM software) and for the wild type HIV-1 PR). Our study try to show that the number and localization of the mutations are factors which induce the HIV-1 PR viral resistance. The larger solvent accessible surface could be recorded for the point mutant Val 82/Phe.     

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)     

Probing the S1/S1' substrate binding pocket geometry of HIV-1 protease with modified aspartic acid analogues

Aspartates 25 and 125, the active site residues of HIV-1 protease, participate functionally in proteolysis by what is believed to be a general acid-general base mechanism. However, the structural role that these residues may play in the formation and maintenance of the neighboring S1/S1' substrate binding pockets remains largely unstudied. Because the active site aspartic acids are essential for catalysis, alteration of these residues to any other naturally occurring amino acid by conventional site-directed mutagenesis renders the protease inactive, and hence impossible to characterize functionally. To investigate whether Asp-25 and Asp-125 may also play a structural role that influences substrate processing, a series of active site protease mutants has been produced in a cell-free protein synthesizing system via readthrough of mRNA nonsense (UAG) codons by chemically misacylated suppressor tRNAs. The suppressor tRNAs were activated with the unnatural aspartic acid analogues erythro-beta-methylaspartic acid, threo-beta-methylaspartic acid, or beta,beta-dimethylaspartic acid. On the basis of the specific activity measurements of the mutants that were produced, the introduction of the beta-methyl moiety was found to alter protease function to varying extents depending upon its orientation. While a beta-methyl group in the erythro orientation was the least deleterious to the specific activity of the protease, a beta-methyl group in the threo orientation, present in the modified proteins containing threo-beta-methylaspartate and beta,beta-dimethylaspartate, resulted in specific activities between 0 and 45% of that of the wild type depending upon the substrate and the substituted active site position. Titration studies of pH versus specific activity and inactivation studies, using an aspartyl protease specific suicide inhibitor, demonstrated that the mutant proteases maintained bell-shaped pH profiles, as well as suicide-inhibitor susceptibilities that are characteristic of aspartyl proteases. A molecular dynamics simulation of the beta-substituted aspartates in position 25 of HIV-1 protease indicated that the threo-beta-methyl moiety may partially obstruct the adjacent S1' binding pocket, and also cause reorganization within the pocket, especially with regard to residues Val-82 and Ile-84. This finding, in conjunction with the biochemical studies, suggests that the active site aspartate residues are in proximity to the S1/S1' binding pocket and may be spatially influenced by the residues presented in these pockets upon substrate binding. It thus seems possible that the catalytic residues cooperatively interact with the residues that constitute the S1/S1' binding pockets and can be repositioned during substrate binding to orient the active site carboxylates with respect to the scissile amide bond, a process that likely affects the facility of proteolysis.     

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.