Non-active Site Mutations Research Articles

The Effect of Clade-Specific Sequence Polymorphisms on HIV-1 Protease Activity and Inhibitor Resistance Pathways

The majority of HIV-1 infections around the world result from non-B clade HIV-1 strains. The CRF01_AE (AE) strain is seen principally in Southeast Asia. AE protease differs by approximately 10% in amino acid sequence from clade B protease and carries several naturally occurring polymorphisms that are associated with drug resistance in clade B. AE protease has been observed to develop resistance through a nonactive-site N88S mutation in response to nelfinavir (NFV) therapy, whereas clade B protease develops both the active-site mutation D30N and the nonactive-site mutation N88D. Structural and biochemical studies were carried out with wild-type and NFV-resistant clade B and AE protease variants. The relationship between clade-specific sequence variations and pathways to inhibitor resistance was also assessed. AE protease has a lower catalytic turnover rate than clade B protease, and it also has weaker affinity for both NFV and darunavir (DRV). This weaker affinity may lead to the nonactive-site N88S variant in AE, which exhibits significantly decreased affinity for both NFV and DRV. The D30N/N88D mutations in clade B resulted in a significant loss of affinity for NFV and, to a lesser extent, for DRV. A comparison of crystal structures of AE protease shows significant structural rearrangement in the flap hinge region compared with those of clade B protease and suggests insights into the alternative pathways to NFV resistance. In combination, our studies show that sequence polymorphisms within clades can alter protease activity and inhibitor binding and are capable of altering the pathway to inhibitor resistance.

Computational Mutation Scanning and Drug Resistance Mechanisms of HIV-1 Protease Inhibitors

The drug resistance of various clinically available HIV-1 protease inhibitors has been studied using a new computational protocol, that is, computational mutation scanning (CMS), leading to valuable insights into the resistance mechanisms and structure-resistance correction of the HIV-1 protease inhibitors associated with a variety of active site and nonactive site mutations. By using the CMS method, the calculated mutation-caused shifts of the binding free energies linearly correlate very well with those derived from the corresponding experimental data, suggesting that the CMS protocol may be used as a generalized approach to predict drug resistance associated with amino acid mutations. Because it is essentially important for understanding the structure-resistance correlation and for structure-based drug design to develop an effective computational protocol for drug resistance prediction, the reasonable and computationally efficient CMS protocol for drug resistance prediction should be valuable for future structure-based design and discovery of antiresistance drugs in various therapeutic areas.

Computational Simulations of the Interactions between Acetyl-Coenzyme-A Carboxylase and Clodinafop: Resistance Mechanism Due to Active and Nonactive Site Mutations

Grass weed populations resistant to acetyl-CoA carboxylase-inhibiting (ACCase; EC 6.4.1.2) herbicides represent a major problem for the sustainable development of modern agriculture. In the present study, extensive computational simulations, including homology modeling, molecular dynamics (MD) simulations, and molecular mechanics-Poisson-Boltzmann surface area (MM/PBSA) calculations, have been carried out to uncover the detailed molecular mechanism of Alopecurus myosuroides resistance to clodinafop, a commercial herbicide targeting ACCase. All the computational model and energetic results indicated that W374C, I388N, D425G, and G443A mutations have great effects on the conformational change of the binding pocket and the hydrogen-bonding interactions. The pi-pi interaction between ligand and the residue of Phe377 and Tyr161', playing an important contribution to the binding affinity, were decreased after mutations. In addition, the hydrogen-bonding interactions between clodinafop and the residues (Ile158' and Ala54') disappeared or decreased significantly upon mutation. As a result, the mutant-type ACCase has a lower affinity for the inhibitor binding than the wild-type enzyme, which accounts for the molecular basis of herbicidal resistance. The structural and mechanistic insights obtained from the present study will provide a valuable clue for future designing of a promising inhibitor to reduce drug resistance associated with both active and nonactive site mutations.

Computational Characterization of Structural Role of the Non-active Site Mutation M36I of Human Immunodeficiency Virus Type 1 Protease

A prominent characteristic of human immunodeficiency virus type 1 (HIV-1) is its high genetic variability, which generates diversity of the virus and often causes a serious problem of the emergence of drug-resistant mutants. Subtype B HIV-1 is dominant in advanced countries, and the mortality rate due to subtype B HIV-1 has been decreased during the past decade. In contrast, the number of patients with non-subtype B viruses is still increasing in developing countries. One of the reasons for the prevalence of non-subtype B viruses is lack of information about the biological and therapeutic differences between subtype B and non-subtype B viruses. M36I is the most frequently observed polymorphism in non-subtype B HIV-1 proteases. However, since the 36th residue is located at a non-active site of the protease and has no direct interaction with any ligands, the structural role of M36I remains unclear. Here, we performed molecular dynamics (MD) simulations of M36I protease in complex with nelfinavir and revealed the influence of the M36I mutation. The results show that M36I regulates the size of the binding cavity of the protease. The reason for the rare emergence of D30N variants in non-subtype B HIV-1 proteases was also clarified from our computational analysis.

Mechanism of Drug Resistance Due to N88S in CRF01_AE HIV-1 Protease, Analyzed by Molecular Dynamics Simulations

Nelfinavir (NFV) is a currently available HIV-1 protease (PR) inhibitor. Patients in whom NFV treatment has failed predominantly carry D30N mutants of HIV-1 PRs if they have been infected with the subtype B virus. In contrast, N88S mutants of HIV-1 PRs predominantly emerge in patients in whom NFV treatment has failed and who carry the CRF01_AE virus. Both D30N and N88S confer resistance against NFV. However, it remains unclear why the nonactive site mutation N88S confers resistance against NFV. In this study, we examined the resistance mechanism through computational simulations. The simulations suggested that despite the nonactive site mutation, N88S causes NFV resistance by reducing interactions between PR and NFV. We also investigated why the emergence rates of D30N and N88S differ between subtype B and CRF01_AE HIV-1. The simulations suggested that polymorphisms of CRF01_AE PR are involved in the emergence rate of the drug-resistant mutants.

Computational Simulations of HIV-1 ProteasesMulti-drug Resistance Due to Nonactive Site Mutation L90M

Human immunodeficiency virus type 1 protease (HIV-1 PR) is one of the proteins that currently available anti-HIV-1 drugs target. Inhibitors of HIV-1 PR have become available, and they have lowered the rate of mortality from acquired immune deficiency syndrome (AIDS) in advanced countries. However, the rate of emergence of drug-resistant HIV-1 variants is quite high because of their short retroviral life cycle and their high mutation rate. Serious drug-resistant mutations against HIV-1 PR inhibitors (PIs) frequently appear at the active site of PR. Exceptionally, some other mutations such as L90M cause drug resistance, although these appear at nonactive sites. The mechanism of resistance due to nonactive site mutations is difficult to explain. In this study, we carried out computational simulations of L90M PR in complex with each of three kinds of inhibitors and one typical substrate, and we clarified the mechanism of resistance. The L90M mutation causes changes in interaction between the side chain atoms of the 90th residue and the main chain atoms of the 25th residue, and a slight dislocation of the 25th residue causes rotation of the side chain at the 84th residue. The rotation of the 84th residue leads to displacement of the inhibitor from the appropriate binding location, resulting in a collision with the flap or loop region. The difference in levels of resistance to the three inhibitors has been explained from energetic and structural viewpoints, which provides the suggestion for promising drugs keeping its efficacy even for the L90M mutant.

Analysis of HIV-1 CRF_01 A/E protease inhibitor resistance: structural determinants for maintaining sensitivity and developing resistance to atazanavir.

A series of HIV-1 protease mutants has been designed in an effort to analyze the contribution to drug resistance provided by natural polymorphisms as well as therapy-selective (active and non-active site) mutations in the HIV-1 CRF_01 A/E (AE) protease when compared to that of the subtype B (B) protease. Kinetic analysis of these variants using chromogenic substrates showed differences in substrate specificity between pretherapy B and AE proteases. Inhibition analysis with ritonavir, indinavir, nelfinavir, amprenavir, saquinavir, lopinavir, and atazanavir revealed that the natural polymorphisms found in A/E can influence inhibitor resistance. It was also apparent that a high level of resistance in the A/E protease, as with B protease, is due to it aquiring a combination of active site and non-active site mutations. Structural analysis of atazanavir bound to a pretherapy B protease showed that the ability of atazanavir to maintain its binding affinity for variants containing some resistance mutations is due to its unique interactions with flap residues. This structure also explains why the I50L and I84V mutations are important in decreasing the binding affinity of atazanavir.

Molecular Dynamics Simulations of HIV-1 Protease Suggest Different Mechanisms Contributing to Drug Resistance.

A major problem in the antiretroviral treatment of HIV-infections with protease-inhibitors is the emergence of resistance, resulting from the occurrence of distinct mutations within the protease molecule. In the present work molecular dynamics simulations of an active-site mutation (D30N) and a nonactive-site mutation (N88S) of HIV-1 protease that both directly confer resistance to the protease inhibitor Nelfinavir but not to Amprenavir were performed and compared to wild-type HIV-protease. A decreased interaction energy between protease and Nelfinavir was observed for the D30N mutant giving a plausible explanation for resistance, while the N88S mutation did not significantly affect the interaction energies in the bound form. Structural analysis including both ligand-bound and unliganded HIV-1 proteases revealed that the free N88S mutant protease shows significant differences in its hydrogen bonding pattern compared to free or Nelfinavir-bound wild-type protease. In particular, Asp30 forms more frequently a hydrogen bond with Ser88 in the unbound N88S mutant thus interfering with the Asp30-Nelfinavir interaction. These findings suggest that different molecular mechanisms contribute to resistance in active-site and nonactive-site mutants and propose a mechanism for the N88S mutant that is based on a shift of the conformational equilibrium of the unbound protease.

Resistant Mechanism against Nelfinavir of Human Immunodeficiency Virus Type 1 Proteases

Inhibitors against human immunodeficiency virus type-1 (HIV-1) proteases are finely effective for anti-HIV-1 treatments. However, the therapeutic efficacy is reduced by the rapid emergence of inhibitor-resistant variants of the protease. Among patients who failed in the inhibitor nelfinavir (NFV) treatment, D30N, N88D, and L90M mutations of HIV-1 protease are often observed. Despite the serious clinical problem, it is not clear how these mutations, especially nonactive site mutations N88D and L90M, affect the affinity of NFV or why they cause the resistance to NFV. In this study, we executed molecular dynamics simulations of the NFV-bound proteases in the wild-type and D30N, N88D, D30N/N88D, and L90M mutants. Our simulations clarified the conformational change at the active site of the protease and the change of the affinity with NFV for all of these mutations, even though the 88th and 90th residues are not located in the NFV-bound cavity and not able to directly interact with NFV. D30N mutation causes the disappearance of the hydrogen bond between the m-phenol group of NFV and the 30th residue. N88D mutation alters the active site conformation slightly and induces a favorable hydrophobic contact. L90M mutation dramatically changes the conformation at the flap region and leads to an unfavorable distortion of the binding pocket of the protease, although 90M is largely far apart from the flap region. Furthermore, the changes of binding energies of the mutants from the wild-type protease are shown to be correlated with the mutant resistivity previously reported by the phenotypic experiments.

Comparing the accumulation of active- and nonactive-site mutations in the HIV-1 protease.

Protease inhibitor resistance still poses one of the greatest challenges in treating HIV. To better design inhibitors able to target resistant proteases, a deeper understanding is needed of the effects of accumulating mutations and the contributions of active- and nonactive-site mutations to the resistance. We have engineered a series of variants containing the nonactive-site mutations M46I and I54V and the active-site mutation I84V. These mutations were added to a protease clone (V6) isolated from a pediatric patient on ritonavir therapy. This variant possessed the ritonavir-resistance-associated mutations in the active-site (V32I and V82A) and nonactive-site mutations (K20R, L33F, M36I, L63P, A71V, and L90M). The I84V mutation had the greatest effect on decreasing catalytic efficiency, 10-fold when compared to the pretherapy clone LAI. The decrease in catalytic efficiency was partially recovered by the addition of mutations M46I and I54V. The M46I and I54V were just as effective at decreasing inhibitor binding as the I84V mutation when compared to V6 and LAI. The V6(54/84) variant showed over 1000-fold decrease in inhibitor-binding strength to ritonavir, indinavir, and nelfinavir when compared to LAI and V6. Crystal-structure analysis of the V6(54/84) variant bound to ritonavir and indinavir shows structural changes in the 80's loops and active site, which lead to an enlarged binding cavity when compared to pretherapy structures in the Protein Data Bank. Structural changes are also seen in the 10's and 30's loops, which suggest possible changes in the dynamics of flap opening and closing.

A non-active site mutation in human hypoxanthine guanine phosphoribosyltransferase expands substrate specificity

Human hypoxanthine guanine phosphoribosyltransferase (HGPRT) lacks the ability to phosphoribosylate xanthine, a property exhibited by HGPRTs from many parasitic protozoa. Using random mutagenesis we have obtained a mutant, F36L, of human HGPRT that phosphoribosylates xanthine. Examination of the structure indicates that F36 does not make direct contact with the purine, but long-range modulation via loop IV, a segment contacting purine at C2 position, could influence substrate specificity. Expanded substrate specificity to include xanthine probably arises from increased flexibility of loop IV as a consequence of mutation at F36. Mutation of the corresponding residue, L44 in Plasmodium falciparum HGPRT, also results in alteration of K m and k cat for xanthine, substantiating its role in affecting purine base affinity. Our studies show that mutation of this residue in the core of the protein also affects the stability of both enzymes.

A non-active site mutation in human hypoxanthine guanine phosphoribosyltransferase expands substrate specificity

A non-active site mutation in human hypoxanthine guanine phosphoribosyltransferase expands substrate specificity

Crystal structures of HIV protease V82A and L90M mutants reveal changes in the indinavir-binding site.

The crystal structures of the wild-type HIV-1 protease (PR) and the two resistant variants, PR(V82A) and PR(L90M), have been determined in complex with the antiviral drug, indinavir, to gain insight into the molecular basis of drug resistance. V82A and L90M correspond to an active site mutation and nonactive site mutation, respectively. The inhibition (K(i)) of PR(V82A) and PR(L90M) was 3.3- and 0.16-fold, respectively, relative to the value for PR. They showed only a modest decrease, of 10-15%, in their k(cat)/K(m) values relative to PR. The crystal structures were refined to resolutions of 1.25-1.4 A to reveal critical features associated with inhibitor resistance. PR(V82A) showed local changes in residues 81-82 at the site of the mutation, while PR(L90M) showed local changes near Met90 and an additional interaction with indinavir. These structural differences concur with the kinetic data.

Structural and Thermodynamic Basis of Resistance to HIV-1 Protease Inhibition: Implications for Inhibitor Design

One of the most serious side effects associated with the therapy of HIV-1 infection is the appearance of viral strains that exhibit resistance to protease inhibitors. At the molecular level, resistance to protease inhibition predominantly takes the form of mutations within the protease molecule that preferentially lower the affinity of protease inhibitors with respect to protease substrates, while still maintaining a viable catalytic activity. Mutations associated with drug resistance occur within the active site cavity as well as distal sites. Active site mutations affect directly inhibitor/protease interactions while non-active site mutations affect inhibitor binding through long range cooperative perturbations. The effects of mutations associated with drug resistance are compounded by the presence of naturally occurring polymorphisms, especially those observed in non-B subtypes of HIV-1. The binding thermodynamics of all clinical inhibitors against the wild type protease, drug resistant mutations and non-B subtype HIV-1 proteases has been determined by high sensitivity isothermal titration calorimetry. In conjunction with structural information, these data have provided a precise characterization of the binding mechanism of different inhibitors and their response to mutations. Inhibitors that exhibit extremely high affinity and low susceptibility to the effects of mutations share common features and binding determinants even if they belong to different chemical scaffolds. These binding determinants define a set of rules and constraints for the design of better HIV-1 protease inhibitors.

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.

Resistance testing in children changing human immunodeficiency virus type 1 protease inhibitor.

To assess genotypic and phenotypic resistance testing in HIV-1-infected children failing a first protease inhibitor (PI) therapy. In a multicenter observational study 21 children, ages 3 to 16 years, were given two reverse transcriptase inhibitors and one PI (mainly ritonavir, n = 18). They were subsequently treated with single or dual PI-based therapy (predominantly nelfinavir, n = 10, or ritonavir-saquinavir, n = 7). Resistance testing was performed at the time of therapy switch via direct sequencing and a recombinant virus susceptibility assay. A total of 21 genotypic and 15 phenotypic resistance profiles were obtained. Most viruses displayed several reverse transcriptase mutations; however, 7 isolates maintained a wild-type protease. Ritonavir targeted the well-known pathway containing 82, 54, 46 and other secondary (nonactive site) mutations including T74A. No in vitro cross-resistance, i.e. > or = 8-fold resistance to saquinavir or amprenavir, was encountered. Secondary mutations enhanced the prediction of ritonavir resistance (i.e. L10I) and in vitro nelfinavir cross-resistance (i.e. K20R/I) conferred by primary (active site) resistance mutations. Either the 82, 54, 46 mutational genotype or the phenotype showing > or = 8-fold nelfinavir cross-resistance predicted a poorer virologic response to nelfinavir salvage therapy. In a small cohort of heavily pretreated pediatric patients, resistance testing appears to predict the response to nelfinavir as salvage for a ritonavir-based therapy. This is further supported by the correlation between ritonavir-selected mutations and in vitro nelfinavir cross-resistance. Prospective studies should assess clinical outcome in children undergoing regimen changes based on resistance testing.

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°C, can both be described by two-state models for both the tryptophans, and the final state reached at 42°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°C, but are nearly the same at 42°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.