Structural and kinetic analysis of drug resistant mutants of HIV-1 protease.

Mutants of HIV-1 protease that are commonly selected on exposure to different drugs, V82S, G48V, N88D and L90M, showed reduced catalytic activity compared to the wild-type protease on cleavage site peptides, CA-p2, p6pol-PR and PR-RT, critical for viral maturation. Mutant V82S is the least active (2-20% of wild-type protease), mutants N88D, R8Q, and L90M exhibit activities ranging from 20 to 40% and G48V from 50 to 80% of the wild-type activity. In contrast, D30N is variable in its activity on different substrates (10-110% of wild-type), with the PR-RT site being the most affected. Mutants K45I and M46L, usually selected in combination with other mutations, showed activities that are similar to (60-110%) or greater than (110-530%) wild-type, respectively. No direct relationship was observed between catalytic activity, inhibition, and structural stability. The mutants D30N and V82S were similar to wild-type protease in their stability toward urea denaturation, while R8Q, G48V, and L90M showed 1.5 to 2.7-fold decreased stability, and N88D and K45I showed 1.6 to 1.7-fold increased stability. The crystal structures of R8Q, K45I and L90M mutants complexed with a CA-p2 analog inhibitor were determined at 2.0, 1.55 and 1.88 A resolution, respectively, and compared to the wild-type structure. The intersubunit hydrophobic contacts observed in the crystal structures are in good agreement with the relative structural stability of the mutant proteases. All these results suggest that viral resistance does not arise by a single mechanism.     

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.     

Replicative fitness of protease inhibitor-resistant mutants of human immunodeficiency virus type 1

The relative replicative fitness of human immunodeficiency virus type 1 (HIV-1) mutants selected by different protease inhibitors (PIs) in vivo was determined. Each mutant was compared to wild type (WT), NL4-3, in the absence of drugs by several methods, including clonal genotyping of cultures infected with two competing viral variants, kinetics of viral antigen production, and viral infectivity/virion particle ratios. A nelfinavir-selected protease D30N substitution substantially decreased replicative capacity relative to WT, while a saquinavir-selected L90M substitution moderately decreased fitness. The D30N mutant virus was also outcompeted by the L90M mutant in the absence of drugs. A major natural polymorphism of the HIV-1 protease, L63P, compensated well for the impairment of fitness caused by L90M but only slightly improved the fitness of D30N. Multiply substituted indinavir-selected mutants M46I/L63P/V82T/I84V and L10R/M46I/L63P/V82T/I84V were just as fit as WT. These results indicate that the mutations which are usually initially selected by nelfinavir and saquinavir, D30N and L90M, respectively, impair fitness. However, additional mutations may improve the replicative capacity of these and other drug-resistant mutants. Hypotheses based on the greater fitness impairment of the nelfinavir-selected D30N mutant are suggested to explain observations that prolonged responses to delayed salvage regimens, including alternate PIs, may be relatively common after nelfinavir failure.     

Drug resistance mutations can effect dimer stability of HIV-1 protease at neutral pH.

The monomer-dimer equilibrium for the human immunodeficiency virus type 1 (HIV-1) protease has been investigated under physiological conditions. Dimer dissociation at pH 7.0 was correlated with a loss in beta-sheet structure and a lower degree of ANS binding. An autolysis-resistant mutant, Q7K/L33I/L63I, was used to facilitate sedimentation equilibrium studies at neutral pH where the wild-type enzyme is typically unstable in the absence of bound inhibitor. The dimer dissociation constant (KD) of the triple mutant was 5.8 microM at pH 7.0 and was below the limit of measurement (approximately 100 nM) at pH 4.5. Similar studies using the catalytically inactive D25N mutant yielded a KD value of 1.0 microM at pH 7.0. These values differ significantly from a previously reported value of 23 nM obtained indirectly from inhibitor binding measurements (Darke et al., 1994). We show that the discrepancy may result from the thermodynamic linkage between the monomer-dimer and inhibitor binding equilibria. Under conditions where a significant degree of monomer is present, both substrates and competitive inhibitors will shift the equilibrium toward the dimer, resulting in apparent increases in dimer stability and decreases in ligand binding affinity. Sedimentation equilibrium studies were also carried out on several drug-resistant HIV-1 protease mutants: V82F, V82F/I84V, V82T/I84V, and L90M. All four mutants exhibited reduced dimer stability relative to the autolysis-resistant mutant at pH 7.0. Our results indicate that reductions in drug affinity may be due to the combined effects of mutations on both dimer stability and inhibitor binding.     

Unique thermodynamic response of tipranavir to human immunodeficiency virus type 1 protease drug resistance mutations

Drug resistance is a major problem affecting the clinical efficacy of antiretroviral agents, including protease inhibitors, in the treatment of infection with human immunodeficiency virus type 1 (HIV-1)/AIDS. Consequently, the elucidation of the mechanisms by which HIV-1 protease inhibitors maintain antiviral activity in the presence of mutations is critical to the development of superior inhibitors. Tipranavir, a nonpeptidic HIV-1 protease inhibitor, has been recently approved for the treatment of HIV infection. Tipranavir inhibits wild-type protease with high potency (K(i) = 19 pM) and demonstrates durable efficacy in the treatment of patients infected with HIV-1 strains containing multiple common mutations associated with resistance. The high potency of tipranavir results from a very large favorable entropy change (-TDeltaS = -14.6 kcal/mol) combined with a favorable, albeit small, enthalpy change (DeltaH = -0.7 kcal/mol, 25 degrees C). Characterization of tipranavir binding to wild-type protease, active site mutants I50V and V82F/I84V, the multidrug-resistant mutant L10I/L33I/M46I/I54V/L63I/V82A/I84V/L90M, and the tipranavir in vitro-selected mutant I13V/V32L/L33F/K45I/V82L/I84V was performed by isothermal titration calorimetry and crystallography. Thermodynamically, the good response of tipranavir arises from a unique behavior: it compensates for entropic losses by actual enthalpic gains or by sustaining minimal enthalpic losses when facing the mutants. The net result is a small loss in binding affinity. Structurally, tipranavir establishes a very strong hydrogen bond network with invariant regions of the protease, which is maintained with the mutants, including catalytic Asp25 and the backbone of Asp29, Asp30, Gly48 and Ile50. Moreover, tipranavir forms hydrogen bonds directly to Ile50, while all other inhibitors do so by being mediated by a water molecule.     

In vivo and in vitro characterization of TEV protease mutants

Tobacco etch virus protease (TEVp) is frequently applied in the cleavage of fusion protein. However, production of TEV protease in Escherichia coli is hampered by low yield and poor solubility, and auto-cleavage of wild type TEVp gives rise to the loss-of-function. Previously it was reported that TEVp S219V displayed more stability, and TEVp variant containing T17S/N68D/I77V and double mutant L56V/S135G resulted in the enhanced production and solubility, respectively. Here, we introduced T17S/N68D/I77V in TEVp S219V to generate TEVpM1 and combined five amino acid mutations (T17S/L56V/N68D/I77V/S135G) in TEVp S219V to create TEVpM2. Among TEVp S219V, and two constructed variants, TEVpM2 displayed highest solubility and catalytic activity in vivo, using EmGFP as the solubility reporter, and the designed fusion protein as in vivo substrate containing an N-terminal hexahistidine tagged GST, a peptide sequence for thrombin and TEV cut and E. coli diaminopropionate ammonia-lyase. The purified TEVp mutants fused with double hexahistidine-tag at N and C terminus showed highest yield, solubility and cleavage efficiency. Mutations of five amino acid residues in TEVpM2 slightly altered protein secondary structure conformed by circular dichroism assay.     

Secondary mutations M36I and A71V in the human immunodeficiency virus type 1 protease can provide an advantage for the emergence of the primary mutation D30N

Development of resistance mutations in enzymatic targets of human immunodeficiency virus 1 (HIV-1) hampers the ability to provide adequate therapy. Of special interest is the effect mutations outside the active site of HIV-1 protease have on inhibitor binding and virus viability. We engineered protease mutants containing the active site mutation D30N alone and with the nonactive site polymorphisms M36I and/or A71V. We determined the K(i) values for the inhibitors nelfinavir, ritonavir, indinavir, KNI272, and AG1776 as well as the catalytic efficiency of the mutants. Single and double mutation combinations exhibited a decrease in catalytic efficiency, while the triple mutant displayed catalytic efficiency greater than that of the wild type. Variants containing M36I or A71V alone did not display a significant change in binding affinities to the inhibitors tested. The variant containing mutation D30N displayed a 2-6-fold increase in K(i) for all inhibitors tested, with nelfinavir showing the greatest increase. The double mutants containing a combination of mutations D30N, M36I, and A71V displayed -0.5-fold to +6-fold changes in the K(i) of all inhibitors tested, with ritonavir and nelfinavir most affected. Only the triple mutant showed a significant increase (>10-fold) in K(i) for inhibitor nelfinavir, ritonavir, or AG-1776 displaying 22-, 19-, or 15-fold increases, respectively. Our study shows that the M36I and A71V mutations provide a greater level of inhibitor cross-resistance combined with active site mutation D30N. M36I and A71V, when present as natural polymorphisms, could aid the virus in developing active site mutations to escape inhibitor binding while maintaining catalytic efficiency.     

Functional correlation between a novel amino acid insertion at codon 19 in the protease of human immunodeficiency virus type 1 and polymorphism in the p1/p6 Gag cleavage site in drug resistance and replication fitness

Population-based sequence analysis revealed the presence of a variant of human immunodeficiency virus type 1 (HIV-1) containing an insertion of amino acid Ile in the protease gene at codon 19 (19I) and amino acid substitutions in the protease at codons 21 (E21D) and 22 (A22V) along with multiple mutations associated with drug resistance, M46I/P63L/A71V/I84V/I93L, in a patient who had failed protease inhibitor (PI) therapy. Longitudinal analysis revealed that the P63L/A71V/I93L changes were present prior to PI therapy. Polymorphisms in the Gag sequence were only seen in the p1/p6 cleavage site at the P1' position (Leu to Pro) and the P5' position (Pro to Leu). To characterize the role of these mutations in drug susceptibility and replication capacity, a chimeric HIV-1 strain containing the 19I/E21D/A22V mutations with the M46I/P63L/A71V/I84V/I93L and p1/p6 mutations was constructed. The chimera displayed high-level resistance to multiple PIs, but not to lopinavir, and grew to 30% of that of the wild type. To determine the relative contribution of each mutation to the phenotypic characteristic of the virus, a series of mutants was constructed using site-directed mutagenesis. A high level of resistance was only seen in mutants containing the 19I/A22V and p1/p6 mutations. The E21D mutation enhanced viral replication. These results suggest that the combination of the 19I/E21D/A22V mutations may emerge and lead to high-level resistance to multiple PIs. The combination of the 19I/A22V mutations may be associated with PI resistance; however, the drug resistance may be caused by the presence of a unique set of mutations in the p1/p6 mutations. The E21D mutation contributes to replication fitness rather than drug resistance.     

Identification of Genotypic Changes in Human Immunodeficiency Virus Protease That Correlate with Reduced Susceptibility to the Protease Inhibitor Lopinavir among Viral Isolates from Protease Inhibitor-Experienced Patients

The association of genotypic changes in human immunodeficiency virus (HIV) protease with reduced in vitro susceptibility to the new protease inhibitor lopinavir (previously ABT-378) was explored using a panel of viral isolates from subjects failing therapy with other protease inhibitors. Two statistical tests showed that specific mutations at 11 amino acid positions in protease (L10F/I/R/V, K20M/R, L24I, M46I/L, F53L, I54L/T/V, L63P, A71I/L/T/V, V82A/F/T, I84V, and L90M) were associated with reduced susceptibility. Mutations at positions 82, 54, 10, 63, 71, and 84 were most closely associated with relatively modest (4- and 10-fold) changes in phenotype, while the K20M/R and F53L mutations, in conjunction with multiple other mutations, were associated with >20- and >40-fold-reduced susceptibility, respectively. The median 50% inhibitory concentrations (IC(50)) of lopinavir against isolates with 0 to 3, 4 or 5, 6 or 7, and 8 to 10 of the above 11 mutations were 0.8-, 2.7-, 13.5-, and 44.0-fold higher, respectively, than the IC(50) against wild-type HIV. On average, the IC(50) of lopinavir increased by 1.74-fold per mutation in isolates containing three or more mutations. Each of the 16 viruses that displayed a >20-fold change in susceptibility contained mutations at residues 10, 54, 63, and 82 and/or 84, along with a median of three mutations at residues 20, 24, 46, 53, 71, and 90. The number of protease mutations from the 11 identified in these analyses (the lopinavir mutation score) may be useful for the interpretation of HIV genotypic resistance testing with respect to lopinavir-ritonavir (Kaletra) regimens and may provide insight into the genetic barrier to resistance to lopinavir-ritonavir in both antiretroviral therapy-naive and protease inhibitor-experienced patients.     

A major role for a set of non-active site mutations in the development of HIV-1 protease drug resistance

A major problem in the chemotherapy of HIV-1 infection is the appearance of drug resistance. In the case of HIV-1 protease inhibitors, resistance originates from mutations in the protease molecule that lower the affinity of inhibitors while still maintaining a viable enzymatic profile. Drug resistance mutations can be classified as active site or non-active site mutations depending on their location within the protease molecule. Active site mutations directly affect drug/target interactions, and their action can be readily understood in structural terms. Non-active site mutations influence binding from distal locations, and their mechanism of action is not immediately apparent. In this paper, we have characterized a mutant form of the HIV-1 protease, ANAM-11, identified in clinical isolates from HIV-1 infected patients treated with protease inhibitors. This mutant protease contains 11 mutations, 10 of which are located outside the active site (L10I/M36I/S37D/M46I/R57K/L63P/A71V/G73S/L90M/I93L) and 1 within the active site (I84V). ANAM-11 lowers the binding affinity of indinavir, nelfinavir, saquinavir, and ritonavir by factors of 4000, 3300, 5800, and 80000, respectively. Surprisingly, most of the loss in inhibitor affinity is due to the non-active site mutations as demonstrated by additional experiments performed with a protease containing only the 10 non-active site mutations (NAM-10) and another containing only the active site mutation (A-1). Kinetic analysis with two different substrates yielded comparable catalytic efficiencies for A-1, ANAM-11, NAM-10, and the wild-type protease. These studies demonstrate that non-active site mutations can be the primary source of resistance and that their role is not necessarily limited to compensate deleterious effects of active site mutations. Analysis of the structural stability of the proteases by differential scanning calorimetry reveals that ANAM-11 and NAM-10 are structurally more stable than the wild-type protease while A-1 is less stable. Together, the binding and structural thermodynamic results suggest that the non-active site mutants affect inhibitor binding by altering the geometry of the binding site cavity through the accumulation of mutations within the core of the protease molecule.     

The complete amino acid sequence of the low molecular weight cytosolic acid phosphatase

This paper presents the complete amino acid sequence of the low molecular weight acid phosphatase from bovine liver. This isoenzyme of the acid phosphatase family is located in the cytosol, is not inhibited by L-(+)-tartrate and fluoride ions, but is inhibited by sulfhydryl reagents. The enzyme consists of 157 amino acid residues, has an acetylated NH2 terminus, and has arginine as the COOH-terminal residue. All 8 half-cystine residues are in the free thiol form. The molecular weight calculated from the sequence is 17,953. The sequence was determined by characterizing the peptides purified by reverse-phase high performance liquid chromatography from tryptic, thermolytic, peptic, Staphylococcus aureus protease, and chymotryptic digests of the carboxymethylated protein. No sequence homologies were found with the two known acylphosphatase isoenzymes or the metalloproteins porcine uteroferrin and purple acid phosphatase from bovine spleen (both of which have acid phosphatase activity). Two half-cystines at or near the active site were identified through the reaction of the enzyme with [14C] iodoacetate in the presence or in the absence of a competitive inhibitor (i.e. inorganic phosphate). Ac-A E Q V T K S V L F V C L G N I C R S P I A E A V F R K L V T D Q N I S D N W V I D S G A V S D W N V G R S P N P R A V S C L R N H G I N T A H K A R Q V T K E D F V T F D Y I L C M D E S N L R D L N R K S N Q V K N C R A K I E L L G S Y D P Q K Q L I I E D P Y Y G N D A D F E T V Y Q Q C V R C C R A F L E K V R-OH.     

Multidrug resistance to HIV-1 protease inhibition requires cooperative coupling between distal mutations

The appearance of viral strains that are resistant to protease inhibitors is one of the most serious problems in the chemotherapy of HIV-1/AIDS. The most pervasive drug-resistant mutants are those that affect all inhibitors in clinical use. In this paper, we have characterized a multiple-drug-resistant mutant of the HIV-1 protease that affects indinavir, nelfinavir, saquinavir, ritonavir, amprenavir, and lopinavir. This mutant (MDR-HM) contains six amino acid mutations (L10I/M46I/I54V/V82A/I84V/L90M) located within and outside the active site of the enzyme. Microcalorimetric and enzyme kinetic measurements indicate that this mutant lowers the affinity of all inhibitors by 2-3 orders of magnitude. By comparison, the multiiple-drug-resistant mutant only increased the K(m) of the substrate by a factor of 2, indicating that the substrate is able to adapt to the changes caused by the mutations and maintain its binding affinity. To understand the origin of resistance, three submutants containing mutations in specific regions were also studied, i.e., the active site (V82A/I84V), flap region (M46I/I54V), and dimerization region (L10I/L90M). None of these sets of mutations by themselves lowered the affinity of inhibitors by more than 1 order of magnitude, and additionally, the sum of the effects of each set of mutations did not add up to the overall effect, indicating the presence of cooperative effects. A mutant containing only the four active site mutations (V82A/I84V/M46I/I54V) only showed a small cooperative effect, suggesting that the mutations at the dimer interface (L10I/L90M) play a major role in eliciting a cooperative response. These studies demonstrate that cooperative interactions contribute an average of 1.2 +/- 0.7 kcal/mol to the overall resistance, most of the cooperative effect (0.8 +/- 0.7 kcal/mol) being mediated by the mutations at the dimerization interface. Not all inhibitors in clinical use are affected the same by long-range cooperative interactions between mutations. These interactions can amplify the effects of individual mutations by factors ranging between 2 and 40 depending on the inhibitor. Dissection of the energetics of drug resistance into enthalpic and entropic components provides a quantitative account of the inhibitor response and a set of thermodynamic guidelines for the design of inhibitors with a lower susceptibility to this type of mutations.     

Mutational analysis of the uracil DNA glycosylase inhibitor protein and its interaction with Escherichia coli uracil DNA glycosylase

Uracil DNA glycosylase inhibitor (Ugi), a protein of 9.4 kDa consists of a five-stranded antiparallel beta sheet flanked on either side by single alpha helices, forms an exclusive complex with uracil DNA glycosylases (UDGs) that is stable in 8M urea. We report on the mutational analysis of various structural elements in Ugi, two of which (hydrophobic pocket and the beta1 edge) establish key interactions with Escherichia coli UDG. The point mutations in helix alpha1 (amino acid residues 3-14) do not affect the stability of the UDG-Ugi complexes in urea. And, while the complex of the deltaN13 mutant with UDG is stable in only approximately 4M urea, its overall structure and thermostability are maintained. The identity of P37, stacked between P26 and W68, was not important for the maintenance of the hydrophobic pocket or for the stability of the complex. However, the M24K mutation at the rim of the hydrophobic pocket lowered the stability of the complex in 6M urea. On the other hand, non-conservative mutations E49G, D61G (cancels the only ionic interaction with UDG) and N76K, in three of the loops connecting the beta strands, conferred no such phenotype. The L23R and S21P mutations (beta1 edge) at the UDG-Ugi interface, and the N35D mutation far from the interface resulted in poor stability of the complex. However, the stability of the complexes was restored in the L23A, S21T and N35A mutations. These analyses and the studies on the exchange of Ugi mutants in preformed complexes with the substrate or the native Ugi have provided insights into the two-step mechanism of UDG-Ugi complex formation. Finally, we discuss the application of the Ugi isolates in overproduction of UDG mutants, toxic to cells.     

Gated binding of ligands to HIV-1 protease: Brownian dynamics simulations in a coarse-grained model

The internal motions of proteins may serve as a "gate" in some systems, which controls ligand-protein association. This study applies Brownian dynamics simulations in a coarse-grained model to study the gated association rate constants of HIV-1 proteases and drugs. The computed gated association rate constants of three protease mutants, G48V/V82A/I84V/L90M, G48V, and L90M with three drugs, amprenavir, indinavir, and saquinavir, yield good agreements with experiments. The work shows that the flap dynamics leads to "slow gating". The simulations suggest that the flap flexibility and the opening frequency of the wild-type, the G48V and L90M mutants are similar, but the flaps of the variant G48V/V82A/I84V/L90M open less frequently, resulting in a lower gated rate constant. The developed methodology is fast and provides an efficient way to predict the gated association rate constants for various protease mutants and ligands.