|Home | About | Journals | Submit | Contact Us | Français|
GRL-02031 (1) is an HIV-1 protease (PR) inhibitor containing a novel P1′ (R)-aminomethyl-2-pyrrolidinone group. Crystal structures at resolutions of 1.25 to 1.55 Å were analyzed for complexes of 1 with the PR containing major drug resistant mutations, PRI47V, PRL76V, PRV82A and PRN88D. Mutations of I47V and V82A alter residues in the inhibitor-binding site, while L76V and N88D are distal mutations having no direct contact with the inhibitor. Substitution of a smaller amino acid in PRI47V and PRL76V, and the altered charge of PRN88D are associated with significant local structural changes compared to the wild-type PRWT, while substitution of alanine in PRV82A increases the size of the S1′ subsite. The P1′ pyrrolidinone group of 1 accommodates to these local changes by assuming two different conformations. Overall, the conformation and interactions of 1 with PR mutants resemble those of PRWT with similar inhibition constants in good agreement with the antiviral potency on multidrug resistant HIV-1.
Human immunodeficiency virus type 1 (HIV-1) protease (PR) is a retroviral aspartyl protease with an essential role in the final step of viral maturation.1 The proteolytic cleavage of Gag and Gag-Pol polyproteins by PR releases the individual viral proteins2 and leads to morphological rearrangements3, which are required for the viral infectivity. Thus PR serves as an excellent target for antiretroviral drugs. Since 1995, nine protease inhibitors (PIs) have been approved by the Food and Drug Administration (FDA) for clinical use together with reverse transcriptase inhibitors and other antiviral agents, as part of the highly active antiretroviral therapy (HAART) regimens.4-6 Despite the success of HAART regimens, the treatments suffer from complications such as drug cytotoxicity, side effects and poor oral bioavailability, as well as drug resistance.7 Probably the most critical obstacle for long-term success of therapy is the emergence of multidrug-resistant HIV-1 variants due to the high drug selection pressure and extremely error-prone viral reverse transcriptase that lacks the proofreading step.8, 9
Among the clinical PIs, darunavir (DRV) exerts high antiviral activity against a wide spectrum of HIV-1 variants10, 11 with the enzyme inhibitory potency in the low picomolar range (Ki = 16 pM).12 DRV is a second generation PI that utilized the “backbone binding” strategy to maximize the interaction between inhibitor and PR backbone atoms.13 The recently described nonpeptidic PI GRL-02031 (1) (Figure 1A) based on the DRV scaffold retains potent activity against laboratory and primary HIV-1 strains.14 Compound 1 offers additional benefits over other clinical PIs with its low dose cytotoxicity (CC50 >100 μM) and high selectivity index (CC50/EC50 >3,600) against HIV-1LAI. Furthermore, compound 1 inhibited a wide spectrum of multidrug-resistant HIV-1 variants with less than 2-fold change in the EC50 relative to the effect on the wild-type HIV-1ERS104pre.15
HIV-1 PR is a member of the aspartyl protease family comprising 99-residues in each monomer of the homodimer. Two conserved triplet residues (Asp25-Thr26-Gly27) interact closely in the center of the active site formed by two monomers.16 Important features of the PR dimer are: the catalytic pocket, the dimer interface and the mobile flaps. Inhibitors bind within a cavity between the catalytic triplets and the flaps (Figure 1B). In the presence of substrate or inhibitor, the flaps undergo conformation changes from an open state to a closed state.17, 18
HIV-1 develops resistance to PIs by selection of mutations in various regions of the PR. Resistance mutations have been observed in nearly forty different residues. Seventeen of the sites are considered major mutations providing high levels of resistance to one or more PIs.19 In this work, we have selected four of the major mutations that are involved in the cross-resistance of PIs. These mutations alter residues in distinct structural environments. Two of the mutations, I47V and V82A, are located inside the inhibitor binding site (Figure 1B). Residues Ile47 and Val82 contact the P2/P2′ and P1/P1′ groups of inhibitor 1, respectively.14, 15 These two residues are located in very different environments with Ile47 buried in a hydrophobic region of the catalytic pocket while Val82 is located on a surface loop of the protein. The other two residues, Leu76 and Asn88, are located in distal regions farther from the catalytic pocket and separated by greater than 5 Å from the inhibitor (Figure 1B). Residue 76 forms hydrophobic contacts with several residues in the hydrophobic core, as described previously.20 In contrast, residue 88 lies on the short helix near the protein surface.
Mutation of Ile47 shows a strong association with decreased susceptibility to amprenavir (APV), darunavir (DRV), lopinavir (LPV) and tipranavir (TPV).19, 21 Similarly, L76V mutation shows decreased susceptibility for APV, DRV and LPV. Interestingly, this mutation has an opposing effect on other drugs as it becomes hypersensitive to atazanavir (ATV), saquinavir (SQV) and TPV.22, 23 Mutations of V82A/T/F/S/L are very common in PI resistance and show reduced susceptibility to all the clinical PIs except DRV.19 Lastly, mutation of N88D alone is associated primarily with resistance to nelfinavir (NFV). However, N88D is common in association with D30N, which exerts a synergistic effect on resistance to NFV, and the double mutant shows 260-fold increased Ki value for NFV compared to wild-type.24-27 To analyze binding of compound 1 with various drug-resistant mutants, crystal structures were determined of its complexes with PRI47V, PRL76V, PRV82A and PRN88D. Analysis of these structures and the relative inhibition values will assist in the design of improved drugs for resistant HIV.
Crystal structures of PRI47V, PRL76V, PRV82A and PRN88D drug-resistant mutants were determined with inhibitor 1 bound in the active site. All four structures were solved in the same space group P21212 as the PRWT/1 complex (PDB ID 3H5B)14 with X-ray data in the high resolution range of 1.25 to 1.55 Å. The data collection and refinement statistics are listed in Table 1. The asymmetric unit contains one homodimer PR with residues labeled 1-99 and 1′-99′. Superposition of wild-type and mutant dimers shows almost identical backbone conformations with a root mean square deviation (RMSD) of 0.12-0.19 Å over 198 Cα atoms. Two opposite orientations of inhibitor 1 were seen in the PRI47V and PRV82A mutants with occupancies of 0.7/0.3 and 0.65/0.35, respectively, comparable to the values of 0.8/0.2 in the PRWT complex. These alternate orientations of inhibitor related by 180° commonly occur due to the pseudosymmetry in the protease dimer. The conformation of the inhibitor is almost identical in all structures with pairwise RMSD of 0.1-0.3 Å for all non-hydrogen atoms; the largest difference is seen for the major orientation in PRI47V. PRL76V and PRN88D mutants showed a single orientation of the inhibitor. Furthermore, the P1′ group of 1 in PRI47V and PRL76V also shows two alternate conformations resembling the major orientation of PRWT/1. For PRI47V mutant, the alternate conformations of the P1′ group have 0.5/0.2 relative occupancy for the major orientation of 1, and a single conformation of P1′ group was refined for the minor inhibitor orientation. In the structure of PRL76V, the alternate conformations of P1′ group have 0.6/0.4 relative occupancy. In all four mutant structures, the inhibitor shows similar numbers of hydrogen bonds and hydrophobic interactions with the enzyme. Both major and minor orientations of inhibitor in the PRI47V and PRL76V structures have comparable interactions with the protease.
The inhibition constants (Ki) for the wild-type and mutants, PRI47V, PRL76V, PRV82A and PRN88D, were determined using a p2/NC based fluorogenic hexapeptide substrate (Table 2). The compound was effective with all four mutants showing less than 2-fold changes in the inhibition constant. Furthermore, Koh et al. previously reported that 1 was a potent PI with comparable activity to DRV.15
The kinetic parameters of the mutants were similar to those of wild-type PR with less than 2-fold changes in the catalytic efficiency (kcat/Km). Interestingly, V82A and N88D mutations have opposite effects on the enzyme activity; PRV82A increases catalytic efficiency and PRN88D decreases catalytic efficiency. In our previous study, the catalytic efficiency of PRN88D had a 2.6 to 4.4-fold reduction over the wild-type value with the tested peptide substrates.28 In contrast, PRV82A showed differences in catalytic efficiency with two substrates based on the MA/CA or CA/p2 cleavage sites.29, 30 The contrasting observations are consistent with the location of the V82A mutation in the substrate binding site, while N88D is a distal mutation without direct effect on the substrate.
Asn88 lies in the N-terminal helix near the protein surface in PRWT. Interestingly, the polar side chain of this residue is directed inward toward the 7th β-strand and loop between the 3rd and 4th β-strand rather than facing the protein surface. The side chain of Asn88 forms a network of several direct and water mediated hydrogen bonds with Asp30, Thr31 and Thr74 connecting to the inhibitor (Figure 2A). The carboxylate group of Asn88 interacts with the hydroxyl group and main chain amide of Thr31 with interatomic distances of 2.8 and 3.2 Å, respectively. The side chain amino group of Asn88 forms a hydrogen bond with the main chain carbonyl oxygen of Thr74 at 2.9 Å and a longer interaction with a water molecule at 3.4 Å. This water molecule is centered between Asp30, Thr31, Thr74 and Asn88 with hydrogen bonds to the carbonyl oxygen of Thr31 and Thr74, main chain amide of Thr31, and the side chain amide of Asn88 (Figure 2A). The mutation in the PRN88D/1 complex introduces a negatively charged side chain showing different conformations in the two subunits. These changes are accompanied by significant alteration in the hydrogen bonding network seen for Asn88 in PRWT. In monomer 1′-99′ of PRN88D, the mutated Asp88′ side chain has rotated 79° outward with a 1.9 Å relative translation at the Cγ position, and in the other monomer, the Asp88 side chain has rotated 14° inward with a 0.4 Å shift at Cγ position (Figure 2B and 2C).
The outward movement of the Asp88′ side chain in the mutant creates a cavity between Thr31′ and Thr74′, and as a result the hydrogen bond interactions with the carbonyl oxygen and side chain hydroxyl group of Thr31′ are abolished. Interestingly, the enzymatic activity is not altered significantly by the hydrogen bond disruption as a water molecule apparently assumes the role of the Asn88′ side chain. The newly created cavity is occupied by a water molecule in a new position that maintains hydrogen bond interactions with Thr31′, Thr74′ and Asp88′ (Figure 2A). Thus the cavity between the 7th β-strand and 3rd/4th β-strand loop generates a more compact structure with a 0.6-0.7 Å backbone shift in the 6th and 7th β-strands of monomer 1′-99′. In addition, the Asp30′ side chain has moved closer to Thr74′, which increases the C-H…O contact distance between the carboxylate group of Asp30′ and methyl group of inhibitor 1 by 0.4 Å (Figure 2B). In contrast, the inward movement of the Asp88 side chain in the other monomer breaks the hydrogen bond interactions with the main chain of Thr74 and shifts the backbone of 6th and 7th β-strand by 0.5-0.6 Å away from residue 88 (Figure 2C). Furthermore, the side chain movement of Asp88 has created an 11° rotation at the carbonyl oxygen to form a shorter hydrogen bond with the hydroxyl group of Thr91 of 2.9 Å in PRN88D relative to 3.2 Å in the PRWT complex. This inward movement of the side chain of Asp88 was also observed in the mutant structures of PRN88D with substrate analog peptides CA/p2 and p2/NC (PDB ID 1FG8 and 1FG6).31 In the PRWT structure, two alternative conformations were modeled for the pyrrolidinone ring of inhibitor 1, while in the PRN88D structure only a single inhibitor conformation was observed. Despite these differences, the protease interaction with inhibitor 1 was not disturbed by the N88D mutation.
Leu76 is located in a hydrophobic environment toward the end of the 7th β-strand with no direct contact to the inhibitor (closest distance of >4.7 Å from P2/P2′ of inhibitor 1). Leu76 is positioned between Asp30 and Gln58 with its side chain directed toward Lys45 and Val56 of the flap. The Leu76 side chain forms multiple hydrophobic contacts with the side chains of Asp30, Val32, Lys45, Ile47, Val56, Gln58 and Thr74.20 The substitution of leucine to the shorter valine in the mutant eliminates several hydrophobic interactions and the side chain rotates by 180° to face Asp30. However, the mutant shows closer hydrophobic contacts between the side chains of the Val76 mutation and Val32 by up to 0.6 Å (Figure 3). Similar changes were also reported in the structures of PRL76V/DRV and PRL76V/SQV (PDB ID 3PWM and 3PWR).20 Aside from the changes in van der Waals interactions, the smaller size of the valine side chain allows movement of the side chains of Asp30/30′ and Lys45/45′ to form an ionic interaction with interatomic distances of 2.6-3.5 Å (Figure 3). A similar effect was seen in the structures of PRL76V/DRV and PRL76V/SQV, although there was a somewhat weaker ionic interaction (3.3-4.2 Å) between Asp30/30′ and Lys45/45′. Overall, the L76V and N88D mutations appear to disrupt the internal van der Waals contacts and hydrogen bonding network, which may alter the stability of the enzyme with minimal effect on the binding of inhibitors. In agreement with the loss of hydrophobic contacts, PRL76V exhibited increased sensitivity to urea and thermal denaturation as well as increased dimer dissociation relative to the PRWT.20
Residue Ile47 lies in the flap, which comprises two anti-parallel β-strands with most of the hydrophobic side chains extended into the catalytic pocket. Ile47 forms hydrophobic interactions with the side chains of Val32, Lys45, Ile54, Val56, Leu76 and Ile50′ from the other flap in the dimer. In addition, Ile47 forms van der Waals contacts with the p-methoxyphenyl and cyclopentanyltetrahydrofuran (Cp-THF) groups of inhibitor 1. Examination of various HIV protease complexes shows that the Ile47 side chain is maintained in a relatively fixed orientation with Cγ1 and Cγ2 positioned toward the side chain of Val56 and the Cδ methyl group extended toward the side chain of Val32 (Figure 4A).
The drug-resistant mutations of Ile47 show very conservative changes with valine as the most frequent substitution. The loss of a methyl group from the isoleucine to valine substitution alters the van der Waals contacts and side chain conformation. In the PRI47V structure, the Val47 side chain has adopted two conformations with equal occupancy. One of the Val47 side chain conformations resembles the wild-type isoleucine with both methyl groups directed toward Val56. Despite this similarity, the mutated Val47 side chain has fewer van der Waals contacts with neighboring residues and abolished contacts with Val32. In the PRWT/1 complex, Ile47 and Ile47′ side chains have a total of 23 van der Waals contacts with the neighboring residues (11 contacts in S2 subsite and 12 contacts in S2′ subsite), whereas Val47 in the mutant has only 12 van der Waals contacts to the neighboring residues. In the other conformation, Val47 shows a significant movement at the main chain with both phi and psi angles rotated by 8-13° and the Cα is shifted by 0.5-0.6 Å. The side chain is rotated 142° and shifted by 0.8-1.0 Å at Cβ and the two methyl groups are shifted closer to Val56 and P2/P2′ of inhibitor. As a result, Val47 in the mutant has lost contacts with Lys45 and shows only 4-6 internal hydrophobic contacts with neighboring residues. Despite the loss of internal interactions, the new side chain position of the alternate conformation of Val47/47′ has gained hydrophobic contacts with the p-methoxyphenyl and Cp-THF of inhibitor 1 (Figure 4B).
Superposition of PRWT and PRI47V inhibitor complexes showed no significant change in the protein backbone with a RMSD of less than 0.2 Å. The only notable variation is seen at residues 79′-82′ of the 80’s loop with 0.6 Å shift at Cα of Pro81′. In the mutant structure, Pro81′ moves away from the oxygen of the pyrrolidinone ring by 0.5 Å and Val82′ moves 0.5 Å closer to the oxygen of the pyrrolidinone ring (Figure 4C). The pyrrolidinone ring has two alternate conformations with a 0.5/0.2 occupancy that favors interaction with the 80’s loop. The two different conformations of pyrrolidinone are also observed in the PRWT complex with a 0.4/0.4 occupancy. Overall, the inhibitor conformations and protease interactions are not significantly altered by the I47V mutation consistent with essentially identical inhibition of 0.4 nM (Table 2).
Mutations of Val82 are common in drug-resistant strains. Superposition of PRWT and PRV82A structures showed minimal changes in the backbone, except for the two surface loops formed by residues 16′-18′ and 79′-82′, which show a Cα deviation of 0.5-0.7 Å (Figure 5A). Residues 16′-18′ form the turn between the 2nd and 3rd β-strand at a region remote from the inhibitor binding site. Residues 79′-82′ are located on the loop between 7th and 8th β-strand and contribute direct contacts to the bound inhibitor. Residues from the 80’s loops and the flaps interact with the P1/P1′ groups of inhibitor 1. Val82 lies perpendicular to the P1 phenylalanine of inhibitor forming extensive C-H…π interactions (Figure 5B). In addition, the phenylalanine forms multiple van der Waals contacts with side chains of Asp25, Pro81, Ile84, Gly27′ and main chain of Gly48′, Gly49′, Ile50′ of the flap. The valine to alanine substitution changes the type of interaction from C-H...π to a single van der Waals contact at 4.1 Å from the P1 group.
In the wild-type S1′ subsite, the major orientation with 80% occupancy of inhibitor 1 shows two distinct conformations of the pyrrolidinone ring; one favors van der Waals interaction, and the other one favors hydrogen bond interaction with PRWT (Figure 5C). One conformation of the pyrrolidinone ring forms van der Waals contacts with Gly49 and Ile50 from the flap, Thr80′, Pro81′, Val82′and Ile84′ from the 80’s loop, and Gly27, Leu23′ and Asp25′ from the catalytic site. The alternative conformation has fewer van der Waals interactions and no contact with Pro81′ and Ile50. The loss of hydrophobic interactions is compensated by the two hydrogen bonds between the pyrrolidinone N15 and the main chain carbonyl oxygen of Gly27, and a water-mediated interaction between the O46 oxygen and Arg8′. In the mutant structure, only a single pyrrolidinone conformation is observed forming hydrogen bond interactions with the Gly27 and water molecules (Figure 5C). The valine to alanine substitution increases the flexibility of the 80’s loop so that the mutated side chain of Ala82′ no longer forms van der Waals contacts with Leu23′ and Ile84′ thus enlarging the S1′ pocket. This shift decreases hydrophobic interactions of the protease and P1′ pyrrolidinone with interatomic distances increased up to 0.7 Å. A similar movement in the 80’s loop was also observed in the PRV82A complexes with inhibitors DRV, IDV and SQV (PDB ID 2IDW, 1SDV and 2NMZ).30, 32, 33 Despite the variation in the P1′ conformation, the inhibition is not affected by the mutation as the pyrrolidinone ring maintains interactions with the mutant independent of changes at Pro81′ and Val82′/V82A′.
Inhibitor 1 exerts potent activity against a wide spectrum of HIV-1 isolates and multidrug-resistant clinical strains.14, 15 Here, we present four high resolution crystal structures of 1 in complex with PRI47V, PRL76V, PRV82A and PRN88D bearing drug resistance mutations. All four single mutants shared comparable inhibition values as the wild-type enzyme with the largest difference of only 2-fold observed in the Ki of PRL76V and PRN88D. In both cases, the mutation increases the mobility of the Asp30 side chain, thus decreasing interactions with the inhibitor. The dynamic side chain position of Asp30 has been observed in various D30N27 single mutants and double mutants such as D30N/N88D and D30N/N88S34. This mutation has significantly reduced the effectiveness of the clinical PIs as observed from the kinetic studies.27 Moreover, D30N is one of the major mutations associated with NFV resistance.19, 35 Despite the variable position of the Asp30 side chain, the p-methoxysulfonamide in 1 maintains a hydrogen bond with the main chain amide of Asp30 and van der Waals contacts with side chain of Ala28 and Asp29. Although a mutation at residue 28 has been observed in the in vitro studies, mutations at Ala28 and Asp29 have not been identified in the clinical PI-resistance isolates.19, 36
The A28S mutation was identified in the UIC-94003 (also referred to as TMC-126)37 or GRL-139838 treated laboratory HIV-1 variants, suggesting an association between the presence of p-methoxysulfonamide in both compounds and the development of A28S resistance mutation. Interestingly, structural analysis and kinetic studies give a different result. Superposition of PRWT/1 and PRA28S/U-89360E (PDB ID 3H5B and 1AXA)39 suggests the A28S side chain may form multiple interactions with p-methoxysulfonamide and Cp-THF of inhibitor 1. Furthermore, kinetic studies showed the mutation of A28S significantly reduces the catalytic activity of the enzyme and viral infectivity, consistent with the absence of this mutation in clinical isolates.40, 41 A further study is needed to clarify the association of p-methoxysulfonamide and Ala28.
In the S2/S2′ subsite, Ile47 lies in the internal hydrophobic region with its side chain aligned parallel to the P2/P2′ group of inhibitor 1. In the mutant structure, I47V loses contacts with neighboring residues, which would potentially decrease the stability of the flaps as suggested by molecular dynamics simulations.42 In addition, the loss of the methyl group would likely affect the binding of substrates and other PIs. As observed in the PRWT/JE-2147 structure (PDB ID 1KZK)43, C5 of Ile47 forms multiple van der Waals contacts with the methylbenzamide group at P2′ position, and the isoleucine to valine mutation would reduce contacts with inhibitor JE-2147 as suggested by Yoshimaru.44 However, compound 1 remains a strong inhibitor of PRI47V. The Cp-THF in the P2 position forms multiple contacts with the flap residues Gly48 and Gly49, thus acting to stabilize the flap and S2 subsite. Notably, the PRI47V structure illustrates an important interaction between Val32 and Ile47.45 The Val32 side chain would potentially restrict the side chain conformations of Ile47 to maintain optimal van der Waals contacts among neighboring residues. This isoleucine-valine pair is retained in other lentivirus PR structures, although the residue types are reversed becoming Ile32 and Val47 in human immunodeficiency virus type 2 (HIV-2) PR, simian immunodeficiency virus (SIV) PR and many HIV-1 drug-resistant mutants.46
The S1/S1′ subsite is composed of the flap and 80s loop from each monomer. One unique feature is the flexibility of the 80s loop that accommodates different peptide or inhibitor groups at P1/P1′.47 A comparison of various inhibitor bound V82A mutant structures shows increased variability of the conformation of 80s loop, which suggested the Val82 side chain plays a role in stabilizing the loop. Val82 shows highly mobile Cγ1 and Cγ2 methyl groups in the electron density maps. The mobility of the Val82 side chain enables hydrophobic contacts with Arg8, Leu23 and/or Ile84, which cannot occur for the short alanine side chain in the PRV82A mutant. The diminished hydrophobic interactions of Ala82′ in the mutant shifts the 80’s loop farther away from the P1′ group of inhibitor 1. Interestingly, the rearrangement of the 80’s loop has a minimal effect on the inhibition. Instead of forming van der Waals contacts with Pro81′, the inhibitor 1 forms hydrogen bond interactions with Gly27 and Arg8′ side chain as described in the previous study of the wild-type complex14. The unusual ability of the P1′-pyrrolidinone ring to assume two alternate conformations has several advantages. First, the inhibitor is less dependent on the P1′ interaction with the 80s loop that is the location of several drug resistance mutations. Second, the inhibitor interactions with PR backbone atoms will be less affected by the high viral mutation rate. Lastly, Gly27 and Arg8′ provide a more stable inhibitor interaction as both residues are highly conserved and excluded from PI-resistant mutations seen in clinical isolates.
The new nonpeptidic PI 1 shows excellent nanomolar inhibition against four of the drug-resistant mutants with the Ki comparable to that of the wild-type enzyme. Moreover, it demonstrated potent antiviral activity against multidrug-resistant clinical variants of HIV-1.14, 15 The crystal structures of inhibitor 1 complexes with four of the major HIV-1 PR drug-resistant mutants, PRI47V, PRL76V, PRV82A and PRN88D, have provided new molecular insights for the interaction of the inhibitor and mutants. Importantly, this structural analysis demonstrates that the P1′ pyrrolidinone group is capable of adapting to fit the binding site of drug resistant mutants. Drug resistance mutations in PR often cause changes in the inhibitor binding site that weaken the PR-inhibitor interactions. One potential improvement to antiviral inhibitors might be the introduction of flexible groups, as observed for compound 1 where the alternate binding conformations of the P1′ pyrrolidinone ring within the S1′ subsite maintain strong inhibition of mutants. Moreover, flexible groups at the P2 or P2′ positions of the inhibitors would permit accommodation to mutations that alter the S2 and S2′ subsites. The flexible groups, however, might weaken the inhibition of wild-type PR. This information can improve the design of the next generation PIs and help to control the emerging multidrug-resistant HIV variants.
Inhibitor 1 has shown analytical purity >99% by HPLC analysis.14 HPLC conditions: linear gradient from 25% acetonitrile and 75% water to 90% acetonitrile and 10% water over 20 min; flow rate 1.5 mL/min, XDB-C18 5 μM, 4.6 × 150 mm column, retention time 5.42 min. HRMS (m/z) calcd for C30H40N3O8S [M+H]+ 602.2536, found 602.2536.
The HIV-1 PR (GenBank HIVHXB2CG) gene was cloned into pET11a expression vector (Novagen) and further optimized to restrict autoproteolysis (Q7K, L33I and L63I) and cysteine-thiol oxidation (C67A and C95A).1 HIV-1 PR point mutants (I47V, L76V, V82A and N88D) were constructed using optimized HIV-1 PR template with Quikchange site-directed mutagenesis kit (Stratagene). Plasmids were transformed into Escherichia coli BL21(DE3) for protein expression. The proteases were purified from inclusion bodies as described previously.48
The kinetic parameters of HIV-1 PR were determined by monitoring the hydrolysis of the fluorescence substrate, Abz-Thr-Ile-Nle-p-nitro-Phe-Gln-Arg-NH2 (where Abz is anthranilic acid, Nle is norleucine and p-nitro-Phe is p-Nitrophenylalanine), that is derived from the p2/NC cleavage site of the viral Gag-Pol polyprotein (Bachem Americas).49 The enzymatic assay was performed as follows; 10 μL protease to final concentration of 12-45 nM was pre-incubated with 98 μL of reaction buffer (0.1 M MES, pH 5.6, 0.4 M NaCl, 1 mM EDTA and 5% glycerol) and 2 μL of DMSO or inhibitor at 37 °C for 5 minutes. The reaction was initiated by addition of 90 μL anthranilyl to a final concentration of 6-96 μM and monitored over a period of 5 minutes using a BMG LABTECH POLARstar OPTIMA 96-well plate reader with excitation wavelength set at 340 nm and emission at 420 nm. The reactions were analyzed using FLUOstar OPTIMA 1.30-0 (BMG LABTECH). Values of Km and kcat were determined by fitting the data to the Michaelis-Menten equation using SigmaPlot 9.01 (Systat Software). Inhibitor Ki values were determined from the fluorescence assay of IC50 dose-response inhibition curve using the equation Ki = (IC50 − [E]/2)/(1 + [S]/Km), where [E] and [S] are protease and anthranilyl concentrations, respectively.48
PRI47V, PRL76V, PRV82A and PRN88D were co-crystallized with compound 1 by the hanging-drop vapor diffusion method at 23 °C. Protein (4.0-7.0 mg/mL) was pre-incubated with the inhibitor at 5-fold molar excess (1:5, protein-inhibitor ratio) for 30 minutes at 4 °C. Each crystallization drop contained equal volumes of protein and reservoir solution (0.1 M sodium acetate buffer, pH 4.2-5.4 and 0.6-1.5 M NaCl). The crystals grew within 1-3 days as bricks, plates or clusters. The cluster crystals were cut and separated into individual crystals before mounting. Prior to the diffraction experiment, crystals were soaked for 30 seconds in reservoir solution containing 30% (v/v) glycerol as cryoprotectant then flash-frozen in liquid nitrogen.
X-ray diffraction data were collected at the SER-CAT beamline of the Advanced Photon Source at Argonne National Laboratory (Argonne, IL). All the experiments were done using a 100 K stream of gaseous nitrogen for crystal cooling. Redundant data sets were collected and processed with the HKL-2000 package.50 Initial phases for PRI47V, PRL76V, PRV82A and PRN88D were determined by molecular replacement using the CCP4 program suite.51 PRWT/DRV (PDB ID 2IEN)32 was used as the initial model for PRI47V and PRL76V. PRWT/1 (PDB ID 3H5B) was used as the initial model for PRV82A and PRN88D. Structural refinement was performed using SHELX52 with the molecular graphics program O53 and Coot54. Alternate conformations were modeled for the protease residues, inhibitor and solvent molecules based on the observed electron density maps. The relative occupancy was calculated during refinement. Anisotropic B-factors were refined for all atoms. Hydrogen atoms were added in the last round of refinement using SHELXL.52 The models were validated with SFCHECK51 and MolProbity55. The atomic coordinates and structure factors were deposited in the Protein Data Bank (PDB).56 Data collection and refinement statistics are listed in Table 1. Structural figures were generated using PyMOL program suite.57
This work was supported in part by the Georgia State University Molecular Basis of Disease Program, the Georgia State University Research Program Enhancement award in Bioinformatics, the Georgia Research Alliance, and the National Institutes of Health grants GM062920 and GM053386. We thank Ting-Yi Chiu and Jing Wang for assisting with construction of PRI47V. We thank the staff at the SER6 CAT beamline at the Advanced Photon Source, Argonne National Laboratory, for assistance during X6 ray data collection. Use of the Advanced Photon Source was supported by the US Department of Energy, Basic Energy Sciences, Office of Science, under Contract No. W-31-109-Eng-38.
PROTEIN DATA BANK: Crystallographic data are available for inhibitor complexes with: PRI47V (PDB ID: 3VF5), PRL76V (PDB ID: 3VF7), PRV82A (PDB ID: 3VFA), and PRN88D (PDB ID: 3VFB).