Toward a universal inhibitor of retroviral proteases: comparative analysis of the interactions of LP-130 complexed with proteases from HIV-1, FIV, and EIAV

Dimer Interface Organization is a Main Determinant of Intermonomeric Interactions and Correlates with Evolutionary Relationships of Retroviral and Retroviral-Like Ddi1 and Ddi2 Proteases

The life cycles of retroviruses rely on the limited proteolysis catalyzed by the viral protease. Numerous eukaryotic organisms also express endogenously such proteases, which originate from retrotransposons or retroviruses, including DNA damage-inducible 1 and 2 (Ddi1 and Ddi2, respectively) proteins. In this study, we performed a comparative analysis based on the structural data currently available in Protein Data Bank (PDB) and Structural summaries of PDB entries (PDBsum) databases, with a special emphasis on the regions involved in dimerization of retroviral and retroviral-like Ddi proteases. In addition to Ddi1 and Ddi2, at least one member of all seven genera of the Retroviridae family was included in this comparison. We found that the studied retroviral and non-viral proteases show differences in the mode of dimerization and density of intermonomeric contacts, and distribution of the structural characteristics is in agreement with their evolutionary relationships. Multiple sequence and structure alignments revealed that the interactions between the subunits depend mainly on the overall organization of the dimer interface. We think that better understanding of the general and specific features of proteases may support the characterization of retroviral-like proteases.

Comparison of a retroviral protease in monomeric and dimeric states

Retroviral proteases (RPs) are of high interest owing to their crucial role in the maturation process of retroviral particles. RPs are obligatory homodimers, with a pepsin-like active site built around two aspartates (in DTG triads) that activate a water molecule, as the nucleophile, under two flap loops. Mason-Pfizer monkey virus (M-PMV) is unique among retroviruses as its protease is also stable in the monomeric form, as confirmed by an existing crystal structure of a 13 kDa variant of the protein (M-PMV PR) and its previous biochemical characterization. In the present work, two mutants of M-PMV PR, D26N and C7A/D26N/C106A, were crystallized in complex with a peptidomimetic inhibitor and one mutant (D26N) was crystallized without the inhibitor. The crystal structures were solved at resolutions of 1.6, 1.9 and 2.0 Å, respectively. At variance with the previous study, all of the new structures have the canonical dimeric form of retroviral proteases. The protomers within a dimer differ mainly in the flap-loop region, with the most extreme case observed in the apo structure, in which one flap loop is well defined while the other flap loop is not defined by electron density. The presence of the inhibitor molecules in the complex structures was assessed using polder maps, but some details of their conformations remain ambiguous. In all of the presented structures the active site contains a water molecule buried deeply between the Asn26-Thr27-Gly28 triads of the protomers. Such a water molecule is completely unique not only in retropepsins but also in aspartic proteases in general. The C7A and C106A mutations do not influence the conformation of the protein. The Cys106 residue is properly placed at the homodimer interface area for a disulfide cross-link, but the reducing conditions of the crystallization experiment prevented S-S bond formation. An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:Acta_Cryst_D:S2059798319011355.

Dynamic and Electrostatic Effects on the Reaction Catalyzed by HIV-1 Protease

HIV-1 Protease (HIV-1 PR) is one of the three enzymes essential for the replication process of HIV-1 virus, which explains why it has been the main target for design of drugs against acquired immunodeficiency syndrome (AIDS). This work is focused on exploring the proteolysis reaction catalyzed by HIV-1 PR, with special attention to the dynamic and electrostatic effects governing its catalytic power. Free energy surfaces for all possible mechanisms have been computed in terms of potentials of mean force (PMFs) within hybrid QM/MM potentials, with the QM subset of atoms described at semiempirical (AM1) and DFT (M06-2X) level. The results suggest that the most favorable reaction mechanism involves formation of a gem-diol intermediate, whose decomposition into the product complex would correspond to the rate-limiting step. The agreement between the activation free energy of this step with experimental data, as well as kinetic isotope effects (KIEs), supports this prediction. The role of the protein dynamic was studied by protein isotope labeling in the framework of the Variational Transition State Theory. The predicted enzyme KIEs, also very close to the values measured experimentally, reveal a measurable but small dynamic effect. Our calculations show how the contribution of dynamic effects to the effective activation free energy appears to be below 1 kcal·mol^-1. On the contrary, the electric field created by the protein in the active site of the enzyme emerges as being critical for the electronic reorganization required during the reaction. These electrostatic properties of the active site could be used as a mold for future drug design.

Structure of RC1339/APRc from Rickettsia conorii, a retropepsin-like aspartic protease

The crystal structures of two constructs of RC1339/APRc from Rickettsia conorii, consisting of either residues 105-231 or 110-231 followed by a His tag, have been determined in three different crystal forms. As predicted, the fold of a monomer of APRc resembles one-half of the mandatory homodimer of retroviral pepsin-like aspartic proteases (retropepsins), but the quaternary structure of the dimer of APRc differs from that of the canonical retropepsins. The observed dimer is most likely an artifact of the expression and/or crystallization conditions since it cannot support the previously reported enzymatic activity of this bacterial aspartic protease. However, the fold of the core of each monomer is very closely related to the fold of retropepsins from a variety of retroviruses and to a single domain of pepsin-like eukaryotic enzymes, and may represent a putative common ancestor of monomeric and dimeric aspartic proteases.

Virtual drug screen schema based on multiview similarity integration and ranking aggregation

The current drug virtual screen (VS) methods mainly include two categories. i.e., ligand/target structure-based virtual screen and that, utilizing protein-ligand interaction fingerprint information based on the large number of complex structures. Since the former one focuses on the one-side information while the later one focuses on the whole complex structure, they are thus complementary and can be boosted by each other. However, a common problem faced here is how to present a comprehensive understanding and evaluation of the various virtual screen results derived from various VS methods. Furthermore, there is still an urgent need for developing an efficient approach to fully integrate various VS methods from a comprehensive multiview perspective. In this study, our virtual screen schema based on multiview similarity integration and ranking aggregation was tested comprehensively with statistical evaluations, providing several novel and useful clues on how to perform drug VS from multiple heterogeneous data sources. (1) 18 complex structures of HIV-1 protease with ligands from the PDB were curated as a test data set and the VS was performed with five different drug representations. Ritonavir ( 1HXW ) was selected as the query in VS and the weighted ranks of the query results were aggregated from multiple views through four similarity integration approaches. (2) Further, one of the ranking aggregation methods was used to integrate the similarity ranks calculated by gene ontology (GO) fingerprint and structural fingerprint on the data set from connectivity map, and two typical HDAC and HSP90 inhibitors were chosen as the queries. The results show that rank aggregation can enhance the result of similarity searching in VS when two or more descriptions are involved and provide a more reasonable similarity rank result. Our study shows that integrated VS based on multiple data fusion can achieve a remarkable better performance compared to that from individual ones and, thus, serves as a promising way for efficient drug screening, taking advantages of the rapidly accumulated molecule representations and heterogeneous data in the pharmacological area.

High-resolution structure of a retroviral protease folded as a monomer

Mason-Pfizer monkey virus (M-PMV), a D-type retrovirus assembling in the cytoplasm, causes simian acquired immunodeficiency syndrome (SAIDS) in rhesus monkeys. Its pepsin-like aspartic protease (retropepsin) is an integral part of the expressed retroviral polyproteins. As in all retroviral life cycles, release and dimerization of the protease (PR) is strictly required for polyprotein processing and virion maturation. Biophysical and NMR studies have indicated that in the absence of substrates or inhibitors M-PMV PR should fold into a stable monomer, but the crystal structure of this protein could not be solved by molecular replacement despite countless attempts. Ultimately, a solution was obtained in mr-rosetta using a model constructed by players of the online protein-folding game Foldit. The structure indeed shows a monomeric protein, with the N- and C-termini completely disordered. On the other hand, the flap loop, which normally gates access to the active site of homodimeric retropepsins, is clearly traceable in the electron density. The flap has an unusual curled shape and a different orientation from both the open and closed states known from dimeric retropepsins. The overall fold of the protein follows the retropepsin canon, but the C(α) deviations are large and the active-site 'DTG' loop (here NTG) deviates up to 2.7 Å from the standard conformation. This structure of a monomeric retropepsin determined at high resolution (1.6 Å) provides important extra information for the design of dimerization inhibitors that might be developed as drugs for the treatment of retroviral infections, including AIDS.

Structural basis for drug and substrate specificity exhibited by FIV encoding a chimeric FIV/HIV protease

A chimeric feline immunodeficiency virus (FIV) protease (PR) has been engineered that supports infectivity but confers sensitivity to the human immunodeficiency virus (HIV) PR inhibitors darunavir (DRV) and lopinavir (LPV). The 6s-98S PR has five replacements mimicking homologous residues in HIV PR and a sixth which mutated from Pro to Ser during selection. Crystal structures of the 6s-98S FIV PR chimera with DRV and LPV bound have been determined at 1.7 and 1.8 Å resolution, respectively. The structures reveal the role of a flexible 90s loop and residue 98 in supporting Gag processing and infectivity and the roles of residue 37 in the active site and residues 55, 57 and 59 in the flap in conferring the ability to specifically recognize HIV PR drugs. Specifically, Ile37Val preserves tertiary structure but prevents steric clashes with DRV and LPV. Asn55Met and Val59Ile induce a distinct kink in the flap and a new hydrogen bond to DRV. Ile98Pro→Ser and Pro100Asn increase 90s loop flexibility, Gln99Val contributes hydrophobic contacts to DRV and LPV, and Pro100Asn forms compensatory hydrogen bonds. The chimeric PR exhibits a comparable number of hydrogen bonds, electrostatic interactions and hydrophobic contacts with DRV and LPV as in the corresponding HIV PR complexes, consistent with IC(50) values in the nanomolar range.

Analysis of structural water and CH···π interactions in HIV-1 protease and PTP1B complexes using a hydrogen bond prediction tool, HBPredicT

A hydrogen bond prediction tool HBPredicT is developed for detecting structural water molecules and CH···π interactions in PDB files of protein-ligand complexes. The program adds the missing hydrogen atoms to the protein, ligands, and oxygen atoms of water molecules and subsequently all the hydrogen bonds in the complex are located using specific geometrical criteria. Hydrogen bonds are classified into various types based on (i) donor and acceptor atoms, and interactions such as (ii) protein-protein, (iii) protein-ligand, (iv) protein-water, (v) ligand-water, (vi) water-water, and (vii) protein-water-ligand. Using the information in category (vii), the water molecules which form hydrogen bonds with the ligand and the protein simultaneously-the structural water-is identified and retrieved along with the associated ligand and protein residues. For CH···π interactions, the relevant portions of the corresponding structures are also extracted in the output. The application potential of this program is tested using 19 HIV-1 protease and 11 PTP1B inhibitor complexes. All the systems showed presence of structural water molecules and in several cases, the CH···π interaction between ligand and protein are detected. A rare occurrence of CH···π interactions emanating from both faces of a phenyl ring of the inhibitor is identified in HIV-1 protease 1D4L.

FIV as a Model for HIV: An Overview

Animal models for human immunodeficiency virus (HIV) infection play a key role in understanding the pathogenesis of AIDS and the development of therapeutic agents and vaccines. As the only lentivirus that causes an immunodeficiency resembling that of HIV infection, in its natural host, feline immunodeficiency virus (FIV) has been a unique and powerful model for AIDS research. FIV was first described in 1987 by Niels Pedersen and co-workers as the causative agent for a fatal immunodeficiency syndrome observed in cats housed in a cattery in Petaluma, California. Since this landmark observation, multiple studies have shown that natural and experimental infection of cats with biological isolates of FIV produces an AIDS syndrome very similar in pathogenesis to that observed for human AIDS. FIV infection induces an acute viremia associated with Tcell alterations including depressed CD4 :CD8 T-cell ratios and CD4 T-cell depletion, peripheral lymphadenopathy, and neutropenia. In later stages of FIV infection, the host suffers from chronic persistent infections that are typically self-limiting in an immunocompetent host, as well as opportunistic infections, chronic diarrhea and wasting, blood dyscracias, significant CD4 T-cell depletion, neurologic disorders, and B-cell lymphomas. Importantly, chronic FIV infection induces a progressive lymphoid and CD4 T-cell depletion in the infected cat. The primary mode of natural FIV transmission appears to be blood-borne facilitated by fighting and biting. However, experimental infection through transmucosal routes (rectal and vaginal mucosa and perinatal) have been well documented for specific FIV isolates. Accordingly, FIV disease pathogenesis exhibits striking similarities to that described for HIV-1 infection.

Structural characterization of B and non-B subtypes of HIV-protease: insights into the natural susceptibility to drug resistance development

Although a majority of HIV-1 infections in Brazil are caused by the subtype B virus (also prevalent in the United States and Western Europe), viral subtypes F and C are also found very frequently. Genomic differences between the subtypes give rise to sequence variations in the encoded proteins, including the HIV-1 protease. The current anti-HIV drugs have been developed primarily against subtype B and the effects arising from the combination of drug-resistance mutations with the naturally existing polymorphisms in non-B HIV-1 subtypes are only beginning to be elucidated. To gain more insights into the structure and function of different variants of HIV proteases, we have determined a 2.1 A structure of the native subtype F HIV-1 protease (PR) in complex with the protease inhibitor TL-3. We have also solved crystal structures of two multi-drug resistant mutant HIV PRs in complex with TL-3, from subtype B (Bmut) carrying the primary mutations V82A and L90M, and from subtype F (Fmut) carrying the primary mutation V82A plus the secondary mutation M36I, at 1.75 A and 2.8 A resolution, respectively. The proteases Bmut, Fwt and Fmut exhibit sevenfold, threefold, and 54-fold resistance to TL-3, respectively. In addition, the structure of subtype B wild type HIV-PR in complex with TL-3 has been redetermined in space group P6(1), consistent with the other three structures. Our results show that the primary mutation V82A causes the known effect of collapsing the S1/S1' pockets that ultimately lead to the reduced inhibitory effect of TL-3. Our results further indicate that two naturally occurring polymorphic substitutions in subtype F and other non-B HIV proteases, M36I and L89M, may lead to early development of drug resistance in patients infected with non-B HIV subtypes.

Evolutionarily conserved functional mechanics across pepsin-like and retroviral aspartic proteases

The biological function of the aspartic protease from HIV-1 has recently been related to the conformational flexibility of its structural scaffold. Here, we use a multistep strategy to investigate whether the same mechanism affects the functionality in the pepsin-like fold. (i) We identify the set of conserved residues by using sequence-alignment techniques. These residues cluster in three distinct regions: near the cleavage-site cavity, in the four beta-sheets cross-linking the two lobes, and in a solvent-exposed region below the long beta-hairpin in the N-terminal lobe. (ii) We elucidate the role played by the conserved residues for the enzymatic functionality of one representative member of the fold family, the human beta-secretase, by means of classical molecular dynamics (MD). The conserved regions exhibit little overall mobility and yet are involved into the most important modes of structural fluctuations. These modes influence the substrate-catalytic aspartates distance through a relative rotation of the N- and C-terminal lobes. (iii) We investigate the effects of this modulation by estimating the reaction free energy at different representative substrate/enzyme conformations. The activation free energy is strongly affected by large-scale protein motions, similarly to what has been observed in the HIV-1 enzyme. (iv) We extend our findings to all other members of the two eukaryotic and retroviral fold families by recurring to a simple, topology-based, energy functional. This analysis reveals a sophisticated mechanism of enzymatic activity modulation common to all aspartic proteases. We suggest that aspartic proteases have been evolutionarily selected to possess similar functional motions despite the observed fold variations.

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.

Inhibitor design by wrapping packing defects in HIV-1 proteins

Two viral proteins, HIV-1 protease and HIV-1 integrase, have been targeted for inhibitor design to prevent assembly and maturation of HIV-1 virions. The enzymatic mechanism of these proteins involves side-chain groups that serve as general acids or bases. Furthermore, catalytic activity requires that water be removed from the microenvironment surrounding the chemical reaction site or be constrained to serve as an activated nucleophile. Here, we identify previously unrecognized structural features that promote water removal from polar catalytic regions. Packing defects in the form of hydrogen bonds that are insufficiently dehydrated intramolecularly, named "dehydrons," are strategically placed in the structure to induce an anhydrous enzymatic pathway. Dehydrons become electrostatically enhanced and stabilized upon further desolvation. Thus, packing defects act synergistically with the polar active groups to enhance the enzymatic electrostatics. However, because dehydrons are sticky, they constitute targets for inhibitor design. We noticed that inhibitors attach to polar surfaces by further desolvating dehydrons, thus blocking the active sites or the sites involved in harnessing the substrate. The dehydrons are thus required for functional reasons, making them suitable targets. The differences in success when targeting HIV-1 protease, feline immunodeficiency virus protease, and HIV-1 integrase are rationalized in terms of the dehydron distribution, revealing possible improvements in the targeting strategy. Principles of design optimization are proposed to create an inhibitor that can be neutralized only at the expense of the loss of catalytic function. The possibility of using drugs that wrap dehydrons to block protein-protein associations is also discussed.

Rabbit endogenous retrovirus-H encodes a functional protease

Recent studies have revealed that 'human retrovirus-5' sequences found in human samples belong to a rabbit endogenous retrovirus family named RERV-H. A part of the gag-pro region of the RERV-H genome was amplified by PCR from DNA in human samples and several forms of RERV-H protease were expressed in bacteria. The RERV-H protease was able to cleave itself from a precursor protein and was also able to cleave the RERV-H Gag polyprotein precursor in vitro whereas a form of the protease with a mutation engineered into the active site was inactive. Potential N- and C-terminal autocleavage sites were characterized. The RERV-H protease was sensitive to pepstatin A, showing it to be an aspartic protease. Moreover, it was strongly inhibited by PYVPheStaAMT, a pseudopeptide inhibitor specific for Mason-Pfizer monkey virus and avian myeloblastosis-associated virus. A structural model of the RERV-H protease was constructed that, together with the activity data, confirms that this is a retroviral aspartic protease.

Unusual binding mode of an HIV-1 protease inhibitor explains its potency against multi-drug-resistant virus strains

Protease inhibitors (PIs) are an important class of drugs for the treatment of HIV infection. However, in the course of treatment, resistant viral variants with reduced sensitivity to PIs often emerge and become a major obstacle to successful control of viral load. On the basis of a compound equipotently inhibiting HIV-1 and 2 proteases (PR), we have designed a pseudopeptide inhibitor, QF34, that efficiently inhibits a wide variety of PR variants. In order to analyze the potency of the inhibitor, we constructed PR species harboring the typical (signature) mutations that confer resistance to commercially available PIs. Kinetic analyses showed that these mutated PRs were inhibited up to 1,000-fold less efficiently by the clinically approved PIs. In contrast, all PR species were effectively inhibited by QF34. In a clinical study, we have monitored 30 HIV-positive patients in the Czech Republic undergoing highly active antiretroviral therapy, and have identified highly PI resistant variants. Kinetic analyses revealed that QF34 retained its subnanomolar potency against multi-drug resistant PR variants. X-ray crystallographic analysis and molecular modeling experiments explained the wide specificity of QF34: this inhibitor binds to the PR in an unusual manner, thus avoiding contact sites that are mutated upon resistance development, and the unusual binding mode and consequently the binding energy is therefore preserved in the complex with a resistant variant. These results suggest a promising route for the design of second-generation PIs that are active against a variety of resistant PR variants.

Combined effect of zidovudine (ZDV), lamivudine (3TC) and abacavir (ABC) antiretroviral therapy in suppressing in vitro FIV replication

In view of close similarities at the molecular and clinical levels, feline immunodeficiency virus (FIV) infection of the domestic cat is subject of increasing attention as an animal model for human immunodeficiency virus (HIV) infection. A range of reverse transcriptase inhibitors effective against HIV are also active against FIV, allowing successful use of the cat model to investigate drug interactions and resistance development. Nevertheless, while combined nucleoside analog and protease inhibitor usage has proven remarkably effective in treating HIV infection, combination antiretroviral therapy of FIV infection has been hampered by lack of protease inhibitors specific for FIV. In an attempt to circumvent this problem, we have examined the feasibility of applying in the FIV system combination protocols lacking a protease inhibitor. We now report that, as observed during HIV infection, the nucleoside analog abacavir (ABC or 1592U89) is able to effectively block in vitro FIV-replication. Furthermore, we demonstrate that combined usage of ABC with the nucleoside analogs zidovudine (ZDV or AZT) and lamivudine (3TC) also blocks in vitro FIV replication in a synergistic manner. However, in contrast to its effect on HIV replication, the ribonucleotide reductase inhibitor hydroxyurea (HU) is unable to effectively control in vitro FIV replication.

Structural implications of drug-resistant mutants of HIV-1 protease: high-resolution crystal structures of the mutant protease/substrate analogue complexes

Emergence of drug-resistant mutants of HIV-1 protease is an ongoing problem in the fight against AIDS. The mechanisms governing resistance are both complex and varied. We have determined crystal structures of HIV-1 protease mutants, D30N, K45I, N88D, and L90M complexed with peptide inhibitor analogues of CA-p2 and p2-NC cleavage sites in the Gag-pol precursor in order to study the structural mechanisms underlying resistance. The structures were determined at 1.55-1.9-A resolution and compared with the wild-type structure. The conformational disorder seen for most of the hydrophobic side-chains around the inhibitor binding site indicates flexibility of binding. Eight water molecules are conserved in all 9 structures; their location suggests that they are important for catalysis as well as structural stability. Structural differences among the mutants were analyzed in relation to the observed changes in protease activity and stability. Mutant L90M shows steric contacts with the catalytic Asp25 that could destabilize the catalytic loop at the dimer interface, leading to its observed decreased dimer stability and activity. Mutant K45I reduces the mobility of the flap and the inhibitor and contributes to an enhancement in structural stability and activity. The side-chain variations at residue 30 relative to wild-type are the largest in D30N and the changes are consistent with the altered activity observed with peptide substrates. Polar interactions in D30N are maintained, in agreement with the observed urea sensitivity. The side-chains of D30N and N88D are linked through a water molecule suggesting correlated changes at the two sites, as seen with clinical inhibitors. Structural changes seen in N88D are small; however, water molecules that mediate interactions between Asn88 and Thr74/Thr31/Asp30 in other complexes are missing in N88D.

1.9 A x-ray study shows closed flap conformation in crystals of tethered HIV-1 PR

Three-dimensional structure of an asymmetrically mutated (C95M) tethered human immunodeficiency virus type 1 protease enzyme (HIV-1 PR) has been determined in an unliganded form using X-ray diffraction data to 1.9 A resolution. The structure, refined using X-PLOR to an R factor of 19.5%, is unexpectedly similar to the ligand-bound native enzyme, rather than to the ligand-free native enzyme. In particular, the two flaps in the tethered dimer are in a closed configuration. The environments around M95 and C1095 are identical, showing no structural effect of this asymmetric mutation at position 95. Oxidation of Cys1095 has been observed for the first time. There is one well-defined water molecule that hydrogen bonds to both carboxyl groups of the essential aspartic acids in the active site. Proteins 2001;43:57-64.

Comparison of inhibitor binding to feline and human immunodeficiency virus proteases: structure-based drug design and the resistance problem

The design and synthesis of compounds targeted against human immunodeficiency virus 1 (HIV-1) protease have resulted in effective antiviral therapies. However, the rapid replication of the virus and the inherent mutability of the viral genome result in the outgrowth of resistant strains in the majority of patients. Thus, there is a continuing need to develop new antiprotease compounds that may bind more effectively to the resistant forms of protease. This contribution examines the binding of a single inhibitor to two different retroviral proteases, HIV-1 protease and feline immunodeficiency virus protease. Despite the overall similarity of the related retroviral enzymes, specific substitutions within the binding site cavity provide a distinctly different binding landscape that dramatically alters the affinity of compounds. Through this comparison, insights have been obtained into new strategies for drug design. New compounds based on these concepts have been tested against the two enzymes.

Structural and biochemical studies of retroviral proteases

Retroviral proteases form a unique subclass of the family of aspartic proteases. These homodimeric enzymes from a number of viral sources have by now been extensively characterized, both structurally and biochemically. The importance of such knowledge to the development of new drugs against AIDS has been, to a large extent, the driving force behind this progress. High-resolution structures are now available for enzymes from human immunodeficiency virus types 1 and 2, simian immunodeficiency virus, feline immunodeficiency virus, Rous sarcoma virus, and equine infectious anemia virus. In this review, structural and biochemical data for retroviral proteases are compared in order to analyze the similarities and differences between the enzymes from different sources and to enhance our understanding of their properties.

Structural studies of FIV and HIV-1 proteases complexed with an efficient inhibitor of FIV protease

Three forms of feline immunodeficiency virus protease (FIV PR), the wild type (wt) and two single point mutants, V59I and Q99V, as well as human immunodeficiency virus type 1 protease (HIV-1 PR), were cocrystallized with the C2-symmetric inhibitor, TL-3. The mutants of FIV PR were designed to replace residues involved in enzyme-ligand interactions by the corresponding HIV-1 PR residues at the structurally equivalent position. TL-3 shows decreased (improved) inhibition constants with these FIV PR mutants relative to wt FIV PR. Despite similar modes of binding of the inhibitor to all PRs (from P3 to P3'), small differences are evident in the conformation of the Phe side chains of TL-3 at the P1 and P1' positions in the complexes with the mutated FIV PRs. The differences mimick the observed binding of TL-3 in HIV-1 PR and correlate with a significant improvement in the inhibition constants of TL-3 with the two mutant FIV PRs. Large differences between the HIV-1 and FIV PR complexes are evident in the binding modes of the carboxybenzyl groups of TL-3 at P4 and P4'. In HIV-1 PR:TL-3, these groups bind over the flap region, whereas in the FIV PR complexes, the rings are located along the major axis of the active site. A significant difference in the location of the flaps in this region of the HIV-1 and FIV PRs correlates with the observed conformational changes in the binding mode of the peptidomimetic inhibitor at the P4 and P4' positions. These findings provide a structural explanation of the observed Ki values for TL-3 with the different PRs and will further assist in the development of improved inhibitors.

Methodology for protein-ligand binding studies: application to a model for drug resistance, the HIV/FIV protease system

A protocol for the rapid energetic analysis of protein-ligand complexes has been developed. This protocol involves the generation of protein-ligand complex ensembles followed by an analysis of the binding free energy components. We apply this methodology toward understanding the origin of binding specificity within the human immunodeficiency virus/feline immunodeficiency virus (HIV/FIV) protease system, a model system for drug resistance studies. A distinct difference in the internal strain of an inhibitor within each protein environment clearly favors the HIV protease complex, as observed experimentally. Our analysis also predicts that residues within the S2-S3 pockets of the FIV protease active site are responsible for this strain. Close examination of the active site residue contributions to interaction energy and desolvation energy identifies specific amino acids that may also play a role in determining the binding preferences of these two enzymes. Proteins 1999;36:318-331.