Crystal structure of a tethered dimer of HIV-1 proteinase complexed with an inhibitor

A Single Amino Acid Substitution at the HIV-1 Protease Termini Dimer Interface Significantly Reduces Viral Particles Processing Efficiency

The dimeric form of HIV-1 protease (PR) is required for its full proteolytic activity. The stability of the dimer primarily depends on the termini interface, with N-terminal residues 1-4 of one monomer encountering C-terminal residues 96-99 of another. We made an alanine substitution for valine 3 (V3) or leucine 97 (L97) at the termini dimer interface and tested their proteolytic activity. We found that an alanine substitution for L97 (PRL97A) completely inhibited the proteolytic activity of the PR. However, an alanine substitution for V3 (PRV3A) partially impaired the proteolytic activity. We then introduced two forced-dimerization systems involving nucleocapsid (NC) replacement or the addition of 1-2 leucine zippers to determine whether the proteolytic activity of dimer-defective PRs could be restored. We found that two forced-dimerization systems compensated for the defect in PRV3A, but not in PRL97A. This implies that PRV3A and PRL97A potentially impair the PR via different mechanisms or cause defects in PR activity to different extents. These novel findings will likely serve as a foundation for developing new PR inhibitors for treating drug-resistant HIV-1 infections in the future.

Fusion-protein-assisted protein crystallization

Fusion proteins can be used directly in protein crystallization to assist crystallization in at least two different ways. In one approach, the `heterologous fusion-protein approach', the fusion partner can provide additional surface area to promote crystal contact formation. In another approach, the `fusion of interacting proteins approach', protein assemblies can be stabilized by covalently linking the interacting partners. The linker connecting the proteins plays different roles in the two applications: in the first approach a rigid linker is required to reduce conformational heterogeneity; in the second, conversely, a flexible linker is required that allows the native interaction between the fused proteins. The two approaches can also be combined. The recent applications of fusion-protein technology in protein crystallization from the work of our own and other laboratories are briefly reviewed.

A Direct Interaction with RNA Dramatically Enhances the Catalytic Activity of the HIV-1 Protease In Vitro

Though the steps of human immunodeficiency virus type 1 (HIV-1) virion maturation are well documented, the mechanisms regulating the proteolysis of the Gag and Gag-Pro-Pol polyproteins by the HIV-1 protease (PR) remain obscure. One proposed mechanism argues that the maturation intermediate p15NC must interact with RNA for efficient cleavage by the PR. We investigated this phenomenon and found that processing of multiple substrates by the HIV-1 PR was enhanced in the presence of RNA. The acceleration of proteolysis occurred independently from the substrate's ability to interact with nucleic acid, indicating that a direct interaction between substrate and RNA is not necessary for enhancement. Gel-shift assays demonstrated the HIV-1 PR is capable of interacting with nucleic acids, suggesting that RNA accelerates processing reactions by interacting with the PR rather than the substrate. All HIV-1 PRs examined have this ability; however, the HIV-2 PR does not interact with RNA and does not exhibit enhanced catalytic activity in the presence of RNA. No specific sequence or structure was required in the RNA for a productive interaction with the HIV-1 PR, which appears to be principally, though not exclusively, driven by electrostatic forces. For a peptide substrate, RNA increased the kinetic efficiency of the HIV-1 PR by an order of magnitude, affecting both turnover rate (k(cat)) and substrate affinity (K(m)). These results suggest that an allosteric binding site exists on the HIV-1 PR and that HIV-1 PR activity during maturation could be regulated in part by the juxtaposition of the enzyme with virion-packaged RNA.

Ligand binding analysis and screening by chemical denaturation shift

The identification of small molecule ligands is an important first step in drug development, especially drugs that target proteins with no intrinsic activity. Toward this goal, it is important to have access to technologies that are able to measure binding affinities for a large number of potential ligands in a fast and accurate way. Because ligand binding stabilizes the protein structure in a manner dependent on concentration and binding affinity, the magnitude of the protein stabilization effect elicited by binding can be used to identify and characterize ligands. For example, the shift in protein denaturation temperature (Tm shift) has become a popular approach to identify potential ligands. However, Tm shifts cannot be readily transformed into binding affinities, and the ligand rank order obtained at denaturation temperatures (≥60°C) does not necessarily coincide with the rank order at physiological temperature. An alternative approach is the use of chemical denaturation, which can be implemented at any temperature. Chemical denaturation shifts allow accurate determination of binding affinities with a surprisingly wide dynamic range (high micromolar to sub nanomolar) and in situations where binding changes the cooperativity of the unfolding transition. In this article, we develop the basic analytical equations and provide several experimental examples.

A single mutation in the core domain of the lac repressor reduces leakiness

BACKGROUND

The lac operon provides cells with the ability to switch from glucose to lactose metabolism precisely when necessary. This metabolic switch is mediated by the lac repressor (LacI), which in the absence of lactose binds to the operator DNA sequence to inhibit transcription. Allosteric rearrangements triggered by binding of the lactose isomer allolactose to the core domain of the repressor impede DNA binding and lift repression. In Nature, the ability to detect and respond to environmental conditions comes at the cost of the encoded enzymes being constitutively expressed at low levels. The readily-switched regulation provided by LacI has resulted in its widespread use for protein overexpression, and its applications in molecular biology represent early examples of synthetic biology. However, the leakiness of LacI that is essential for the natural function of the lac operon leads to an increased energetic burden, and potentially toxicity, in heterologous protein production.

RESULTS

Analysis of the features that confer promiscuity to the inducer-binding site of LacI identified tryptophan 220 as a target for saturation mutagenesis. We found that phenylalanine (similarly to tryptophan) affords a functional repressor that is still responsive to IPTG. Characterisation of the W220F mutant, LacIWF, by measuring the time dependence of GFP production at different IPTG concentrations and at various incubation temperatures showed a 10-fold reduction in leakiness and no decrease in GFP production. Cells harbouring a cytotoxic protein under regulatory control of LacIWF showed no decrease in viability in the early phases of cell growth. Changes in responsiveness to IPTG observed in vivo are supported by the thermal shift assay behaviour of purified LacIWF with IPTG and operator DNA.

CONCLUSIONS

In LacI, long-range communications are responsible for the transmission of the signal from the inducer binding site to the DNA binding domain and our results are consistent with the involvement of position 220 in modulating these. The mutation of this single tryptophan residue to phenylalanine generated an enhanced repressor with a 10-fold decrease in leakiness. By minimising the energetic burden and cytotoxicity caused by leakiness, LacIWF constitutes a useful switch for protein overproduction and synthetic biology.

Structural, kinetic, and thermodynamic studies of specificity designed HIV-1 protease

HIV-1 protease recognizes and cleaves more than 12 different substrates leading to viral maturation. While these substrates share no conserved motif, they are specifically selected for and cleaved by protease during viral life cycle. Drug resistant mutations evolve within the protease that compromise inhibitor binding but allow the continued recognition of all these substrates. While the substrate envelope defines a general shape for substrate recognition, successfully predicting the determinants of substrate binding specificity would provide additional insights into the mechanism of altered molecular recognition in resistant proteases. We designed a variant of HIV protease with altered specificity using positive computational design methods and validated the design using X-ray crystallography and enzyme biochemistry. The engineered variant, Pr3 (A28S/D30F/G48R), was designed to preferentially bind to one out of three of HIV protease's natural substrates; RT-RH over p2-NC and CA-p2. In kinetic assays, RT-RH binding specificity for Pr3 increased threefold compared to the wild-type (WT), which was further confirmed by isothermal titration calorimetry. Crystal structures of WT protease and the designed variant in complex with RT-RH, CA-p2, and p2-NC were determined. Structural analysis of the designed complexes revealed that one of the engineered substitutions (G48R) potentially stabilized heterogeneous flap conformations, thereby facilitating alternate modes of substrate binding. Our results demonstrate that while substrate specificity could be engineered in HIV protease, the structural pliability of protease restricted the propagation of interactions as predicted. These results offer new insights into the plasticity and structural determinants of substrate binding specificity of the HIV-1 protease.

Comparison of structure- and ligand-based virtual screening protocols considering hit list complementarity and enrichment factors

Structure- and ligand-based virtual-screening methods (docking, 2D- and 3D-similarity searching) were analysed for their effectiveness in virtual screening against four different targets: angiotensin-converting enzyme (ACE), cyclooxygenase 2 (COX-2), thrombin and human immunodeficiency virus 1 (HIV-1) protease. The relative performance of the tools was compared by examining their ability to recognise known active compounds from a set of actives and nonactives. Furthermore, we investigated whether the application of different virtual-screening methods in parallel provides complementary or redundant hit lists. Docking was performed with GOLD, Glide, FlexX and Surflex. The obtained docking poses were rescored by using nine different scoring functions in addition to the scoring functions implemented as objective functions in the docking algorithms. Ligand-based virtual screening was done with ROCS (3D-similarity searching), Feature Trees and Scitegic Functional Fingerprints (2D-similarity searching). The results show that structure- and ligand-based virtual-screening methods provide comparable enrichments in detecting active compounds. Interestingly, the hit lists that are obtained from different virtual-screening methods are generally highly complementary. These results suggest that a parallel application of different structure- and ligand-based virtual-screening methods increases the chance of identifying more (and more diverse) active compounds from a virtual-screening campaign.

HIV-1 protease function and structure studies with the simplicial neighborhood analysis of protein packing method

The Simplicial Neighborhood Analysis of Protein Packing (SNAPP) method was used to predict the effect of mutagenesis on the enzymatic activity of the HIV-1 protease (HIVP). SNAPP relies on a four-body statistical scoring function derived from the analysis of spatially nearest neighbor residue compositional preferences in a diverse and representative subset of protein structures from the Protein Data Bank. The method was applied to the analysis of HIVP mutants with residue substitutions in the hydrophobic core as well as at the interface between the two protease monomers. Both wild-type and tethered structures were employed in the calculations. We obtained a strong correlation, with R(2) as high as 0.96, between DeltaSNAPP score (i.e., the difference in SNAPP scores between wild-type and mutant proteins) and the protease catalytic activity for tethered structures. However, a weaker but significant correlation was obtained for nontethered structures. Our analysis identified residues both in the hydrophobic core and at the dimeric interface that are very important for the protease function. This study demonstrates a potential utility of the SNAPP method for rational design of mutagenesis studies and protein engineering.

Mutagenesis of the dimer interface residues of tethered and untethered HIV-1 protease result in differential activity and suggest multiple mechanisms of compensation

As is the case for all retroviruses, the protease of HIV-1 is only functional as a homodimer; dimerization of two protease monomers results in the formation of the enzyme active site. This dimer structure is supported primarily by interactions between the first four amino-terminal and the last four carboxy-terminal amino acids. These eight amino acids form a beta-sheet in which hydrophobic residues are oriented towards the core of the molecule and polar residues are directed towards the solvent. Although the structure of the dimer interface has been determined, the forces that support dimerization have not been fully characterized. Here, we describe a tethered construct in which two protease monomers are joined by a 5 amino acid linker. We evaluate the relative role of each dimer interface residue in functional homo- and heterodimers. Our studies indicate that the hydrophobic residues of the dimer interface are particularly important in maintaining enzyme activity and that enzyme activity is more sensitive to substitutions of the C-terminal amino acids. Further, we demonstrate that the presence of the tether is able to compensate for mutations within the dimer interface that inactivate the enzyme.

Real time NMR monitoring of local unfolding of HIV-1 protease tethered dimer driven by autolysis

Structural studies in proteases have been hampered because of their inherent autolytic function. However, since autolysis is known to be mediated via protein unfolding, careful monitoring of the autolytic reaction has the potential to throw light on the folding-unfolding equilibria. In this paper we describe real time nuclear magnetic resonance investigations on the tethered dimer construct of the human immunodeficiency virus-1 protease, which have yielded insights into the relative stabilities of several residues in the protein. The residues lying along the active site (bottom, side and top of the active site) and those in helix have lower unfolding free energy values than the other parts of the protein. The residue level stability differences suggest that the protein is well suited to adjust itself in almost all the regions of its structure, as and when perturbations occur, either due to ligand binding or due to mutations.

How does a symmetric dimer recognize an asymmetric substrate? A substrate complex of HIV-1 protease

The crystal structure of an actual HIV-1 protease-substrate complex is presented at 2.0 A resolution (R-value of 19.7 % (R(free) 23.3 %)) between an inactive variant (D25N) of HIV-1 protease and a long substrate peptide, Lys-Ala-Arg-Val-Leu-Ala-Glu-Ala-Met-Ser, which covers a full binding epitope of capsid(CA)-p2, cleavage site. The substrate peptide is asymmetric in both size and charge distribution. To accommodate this asymmetry the two protease monomers adopt different conformations burying a total of 1038 A(2) of surface area at the protease-substrate interface. The specificity for the CA-p2 substrate peptide is mainly hydrophobic, as most of the hydrogen bonds are made with the backbone of the peptide substrate. Two water molecules bridge the two monomers through the loops Gly49-Gly52 (Gly49'-Gly52') and Pro79'-Val82' (Pro79-Val82). When other complexes are compared, the mobility of these loops is correlated with the content of the P1 and P1' sites. Interdependence of the conformational changes allows the protease to exhibit its wide range of substrate specificity.

Small dipeptide-based HIV protease inhibitors containing the hydroxymethylcarbonyl isostere as an ideal transition-state mimic

The human immunodeficiency virus (HIV) codes for an aspartic protease known to be essential for retroviral maturation and replication. HIV protease is formed from two identical 99 amino acid peptides. We synthesized [(NHCH2CH2-S-CH2CO)51-52, Ala67,95]HIV-1 protease using the thioether chemical ligation method, and then prepared the [(NHCH2CH2-S-CH2CO)51-52, Ala67,95, Cys98]HIV-1 protease dimer analogue covalently linked by a disulfide bridge. These HIV-1 protease analogues effectively cleaved the Tyr-Phe-type substrate, but had weak affinity to the Tyr-Pro-type substrate. Consequently, the molecular recognition of the protease analogues differs from that of the wild-type enzyme. Based on the substrate transition state, we designed and synthesized a novel class of HIV protease inhibitors containing an unnatural amino acid, (2S, 3S)-3-amino-2-hydroxy-4-phenylbutyric acid, named allophenylnorstatine, with a hydroxymethylcarbonyl (HMC) isostere. The stereochemistry of the hydroxyl group was significant for the enzyme inhibition and the HMC group interacted excellently with the aspartic acid carboxyl groups of HIV protease active site in the essentially same hydrogen-bonding mode as the transition state. Small dipeptide-based HIV protease inhibitors containing the HMC isostere were studied as advantageous compounds. Among them, a dipeptide-based HIV protease inhibitor, KNI-577, exhibited potent antiviral activities, low cytotoxicity, and good pharmacokinetic properties.

Unfolding kinetics of tryptophan side chains in the dimerization and hinge regions of HIV-I protease tethered dimer by real time NMR spectroscopy

HIV I protease has been the target of extensive and variety of investigations in recent years because of its importance in the AIDS viral life cycle. We describe here real time NMR studies on the unfolding kinetics of two tryptophans, W6 and W42, which are located in the dimerization and hinge domains of the protein, respectively. Unfolding seems to get initiated in the dimerization domain. The kinetic data at two temperatures, 32 and 42 degrees C, can both be described by two-state models for both the tryptophans, and the final state reached at 42 degrees C does not depend on the path of unfolding. Unfolding free energy changes derived from the kinetic fitting parameters are less than 3 kJ/mol, indicating that the energy landscape is very shallow. The free energy values and the rates for the two tryptophans are different at 32 degrees C, but are nearly the same at 42 degrees C. These are interpreted in the light of the "new view" of protein folding and the relative behaviors of the two tryptophans suggest the existence of cooperative pathways in the unfolding reaction of the protein. These observations would provide valuable insights into protein function, stability, and effects of nonactive site mutations conferring drug resistance.

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.

DREAM++: flexible docking program for virtual combinatorial libraries

We present a set of programs, DREAM+2 (Docking and Reaction programs using Efficient seArch Methods written in C++), for docking computationally generated ligands into macromolecular binding sites. DREAM++ is composed of three programs: ORIENT++, REACT++ and SEARCH++. The program ORIENT++ positions molecules in a binding site with the DOCK algorithm. Its output can be used as input to REACT++ and SEARCH+2. The program REACT++ performs user-specific chemical reactions on a docked molecule, so that reaction products can be evaluated for three dimensional complementarity with the macromolecular site. The program SEARCH++ performs an efficient conformation search on the reaction products using a hybrid backtrack and incremental construction algorithm. We have applied the programs to HIV protease-inhibitor complexes as test systems. We found that we can differentiate high-affinity ligands based on several measures: interaction energies, occupancy of protein subsites and the number of successfully docked conformations for each product. Encouraged by the results in the test case, we applied the programs to propose novel inhibitors of HIV protease. These inhibitors can be generated by organic reactions using commercially available reagents. They are alternatives to the inhibitors synthesized by Glaxo.

Inhibitors of HIV-1 protease: a major success of structure-assisted drug design

Retroviral protease (PR) from the human immunodeficiency virus type 1 (HIV-1) was identified over a decade ago as a potential target for structure-based drug design. This effort was very successful. Four drugs are already approved, and others are undergoing clinical trials. The techniques utilized in this remarkable example of structure-assisted drug design included crystallography, NMR, computational studies, and advanced chemical synthesis. The development of these drugs is discussed in detail. Other approaches to designing HIV-1 PR inhibitors, based on the concepts of symmetry and on the replacement of a water molecule that had been found tetrahedrally coordinated between the enzyme and the inhibitors, are also discussed. The emergence of drug-induced mutations of HIV-1 PR leads to rapid loss of potency of the existing drugs and to the need to continue the development process. The structural basis of drug resistance and the ways of overcoming this phenomenon are mentioned.

Crystallographic analysis of human immunodeficiency virus 1 protease with an analog of the conserved CA-p2 substrate -- interactions with frequently occurring glutamic acid residue at P2' position of substrates

Human immunodeficiency virus type 1 (HIV-1) protease hydrolysis of the Gag CA-p2 cleavage site is crucial for virion maturation and is optimal at acidic pH. To understand the processing of the CA-p2 site, we have determined the structure of HIV-1 protease complexed with an analog of the CA-p2 site, the reduced peptide inhibitor Arg-Val-Leu-r-Phe-Glu-Ala-Ahx-NH2 [r denotes the reduced peptide bond and Ahx 2-aminohexanoic acid (norleucine), respectively]. The crystal structure was refined to an R-factor of 0.17 at 0.21-nm resolution. The crystals have nearly the same lattice as related complexes in P2(1)2(1)2(1) which have twofold disordered inhibitor, but are in space group P2(1). and the asymmetric unit contains two dimers of HIV-1 protease related by 180 degrees rotation. An approximate non-crystallographic symmetry has replaced the exact crystal symmetry resulting in well-ordered inhibitor structure. Each protease dimer binds one ordered inhibitor molecule, but in opposite orientations. The interactions of the inhibitor with the two dimers are very similar for the central P2 Val to P2' Glu residues, but show more variation for the distal P3 Arg and P4' Ahx residues. Importantly, the carboxylate oxygens of Glu at P2' in the inhibitor are within hydrogen-bonding distance of a carboxylate oxygen of Asp30 of the protease suggesting that the two side chains share a proton. This interaction suggests that the enzyme-substrate complex is additionally stabilized at lower pH. The importance of this interaction is emphasized by the absence of polymorphisms of Asp30 in the protease and variants of P2' Glu in the critical CA-p2 cleavage site.

The 80's loop (residues 78 to 85) is important for the differential activity of retroviral proteases

The abundance of structural data available for retroviral proteases affords a unique opportunity to investigate structure activity relationships. Our approach attempts to genetically engineer an HIV (human immunodeficiency virus)-1 protease that is functionally equivalent to the HIV-2 and the SIV (simian immunodeficiency virus) enzymes and conversely to engineer an HIV-2 protease that is functionally equivalent to the HIV-1 enzyme. For this purpose, the HIV-2 and SIV proteases were cloned and characterized in an Escherichia coli (E. coli) assay system along with 33 engineered HIV-1 and HIV-2 enzymes. The results of these experiments show that a relatively large S1 or S1' subsite volume, which is likely determined by the conformation of the 80's loop (residues 78 to 85), is necessary to fully accommodate the HIV-1 protease specificity site AETF*YCDG (the asterisk indicates the location scissile bond) during productive binding.

Activity of tethered human immunodeficiency virus 1 protease containing mutations in the flap region of one subunit

The tethered-dimer protease of human immunodeficiency virus 1 (HIV-1) [Cheng Y.-S. E., Yin, F.H., Foundling, S., Blomstrom, D. & Kettner, C. A. (1990) Proc. Natl Acad. Sci. USA 87, 9660-9664] and its mutants containing amino acid substitutions or deletions or both in only one flap region were expressed in Escherichia coli. These mutant enzymes showed various degrees of self-processing and significantly reduced catalytic activity toward oligopeptide substrates compared with the wild type. Kinetic parameters determined for one of the oligopeptide substrates showed a dramatic increase in K(m) and decrease in Kcat values. Unexpectedly, the substrate cleavage was more efficient in low salt concentration for a mutant containing a shortened hydrophilic flap. Assays with oligopeptides representing naturally occurring cleavage sites or oligopeptides containing single amino acid substitutions at the P2 and P2' substrate positions showed only moderate changes in the substrate specificity of the mutant proteases. Predicted structures for the mutants were constructed by molecular modeling and used to interpret the results of kinetic measurements. In general, the data suggest that the mutated part of the flaps does not have a major role in determining substrate specificity; rather, it provides the hydrophobic environment and hydrogen-bond interactions with the conserved water that are necessary for efficient substrate binding and catalysis.

Structure of HIV-1 protease with KNI-272, a tight-binding transition-state analog containing allophenylnorstatine

BACKGROUND

HIV-1 protease (HIV PR), an aspartic protease, cleaves Phe-Pro bonds in the Gag and Gag-Pol viral polyproteins. Substrate-based peptide mimics constitute a major class of inhibitors of HIV PR presently being developed for AIDS treatment. One such compound, KNI-272, which incorporates allophenylnorstatine (Apns)-thioproline (Thp) in place of Phe-Pro, has potent antiviral activity and is undergoing clinical trials. The structure of the enzyme-inhibitor complex should lead to an understanding of the structural basis for its tight binding properties and provide a framework for interpreting the emerging resistance to this drug.

RESULTS

The three-dimensional crystal structure of KNI-272 bound to HIV PR has been determined to 2.0 A resolution and used to analyze structure-activity data and drug resistance for the Arg8-->Gln and ILe84-->Val mutations in HIV PR. The conformationally constrained Apns-Thp linkage is favorably recognized in its low energy trans conformation, which results in a symmetric mode of binding to the active-site aspartic acids and also explains the unusual preference of HIV PR for the S, or syn, hydroxyl group of the Apns residue. The inhibitor recognizes the enzyme via hydrogen bonds to three bridging water molecules, including one that is coordinated directly to the catalytic Asp125 residue.

CONCLUSIONS

The structure of the HIV PR/KNI-272 complex illustrates the importance of limiting the conformational degrees of freedom and of using protein-bound water molecules for building potent inhibitors. The binding mode of HIV PR inhibitors can be predicted from the stereochemical relationship between adjacent hydroxyl-bearing and side chain bearing carbon atoms of the P1 substituent. Our structure also provides a framework for designing analogs targeted to drug-resistant mutant enzymes.