Exploring the drug resistance mechanism of active site,non-active site mutations and their cooperative effects in CRF01_AE HIV-1 protease: molecular dynamics simulations and free energy calculations

Human immunodeficiency virus type 1 protease is essential for virus replication and maturation and has been considered as one of the important drug target for the antiretroviral treatment of HIV infection. The majority of HIV infections are caused due to non-B subtypes in developing countries. Subtype AE is spreading rapidly and infecting huge population worldwide. Understanding the interdependence of active and non-active site mutations in conferring drug resistance is crucial for the development effective inhibitors in subtype AE protease. In this work, we have investigated the mechanism of resistance against indinavir (IDV) due to therapy selected active site mutation V82F, non-active site mutations PF82V and their cooperative effects PV82F in subtype AE-protease using molecular dynamics simulations and binding free energy calculations. The simulations suggested all the three complexes lead to decrease in binding affinity of IDV, whereas the PF82V complex resulted in an enhanced binding affinity compared to V82F and PV82F complexes. Large positional deviation of IDV was observed in V82F complex. The preservation of hydrogen bonds of IDV with active site Asp25/Asp25′ and flap residue Ile50/50′ via a water molecule is crucial for effective binding. Owing to the close contact of 80s loop with Ile50′ and Asp25, the alteration between residues Thr80 and Val82, further induces conformational change thereby resulting in loss of interactions between IDV and the residues in the active site cavity, leading to drug resistance. Our present study shed light on the effect of active, non-active site mutations and their cooperative effects in AE protease.

Communicated by Ramaswamy H. Sarma     

Accessory mutations balance the marginal stability of the HIV-1 protease in drug resistance

The HIV-1 protease is a major target of inhibitor drugs in AIDS therapies. The therapies are impaired by mutations of the HIV-1 protease that can lead to resistance to protease inhibitors. These mutations are classified into major mutations, which usually occur first and clearly reduce the susceptibility to protease inhibitors, and minor, accessory mutations that occur later and individually do not substantially affect the susceptibility to inhibitors. Major mutations are predominantly located in the active site of the HIV-1 protease and can directly interfere with inhibitor binding. Minor mutations, in contrast, are typically located distal to the active site. A central question is how these distal mutations contribute to resistance development. In this article, we present a systematic computational investigation of stability changes caused by major and minor mutations of the HIV-1 protease. As most small single-domain proteins, the HIV-1 protease is only marginally stable. Mutations that destabilize the folded, active state of the protease therefore can shift the conformational equilibrium towards the unfolded, inactive state. We find that the most frequent major mutations destabilize the HIV-1 protease, whereas roughly half of the frequent minor mutations are stabilizing. An analysis of protease sequences from patients in treatment indicates that the stabilizing minor mutations are frequently correlated with destabilizing major mutations, and that highly resistant HIV-1 proteases exhibit significant fractions of stabilizing mutations. Our results thus indicate a central role of minor mutations in balancing the marginal stability of the protease against the destabilization induced by the most frequent major mutations.     

Combining mutations in HIV-1 protease to understand mechanisms of resistance

HIV-1 develops resistance to protease inhibitors predominantly by selecting mutations in the protease gene. Studies of resistant mutants of HIV-1 protease with single amino acid substitutions have shown a range of independent effects on specificity, inhibition, and stability. Four double mutants, K45I/L90M, K45I/V82S, D30N/V82S, and N88D/L90M were selected for analysis on the basis of observations of increased or decreased stability or enzymatic activity for the respective single mutants. The double mutants were assayed for catalysis, inhibition, and stability. Crystal structures were analyzed for the double mutants at resolutions of 2.2-1.2 A to determine the associated molecular changes. Sequence-dependent changes in protease-inhibitor interactions were observed in the crystal structures. Mutations D30N, K45I, and V82S showed altered interactions with inhibitor residues at P2/P2', P3/P3'/P4/P4', and P1/P1', respectively. One of the conformations of Met90 in K45I/L90M has an unfavorably close contact with the carbonyl oxygen of Asp25, as observed previously in the L90M single mutant. The observed catalytic efficiency and inhibition for the double mutants depended on the specific substrate or inhibitor. In particular, large variation in cleavage of p6(pol)-PR substrate was observed, which is likely to result in defects in the maturation of the protease from the Gag-Pol precursor and hence viral replication. Three of the double mutants showed values for stability that were intermediate between the values observed for the respective single mutants. D30N/V82S mutant showed lower stability than either of the two individual mutations, which is possibly due to concerted changes in the central P2-P2' and S2-S2' sites. The complex effects of combining mutations are discussed.     

Modeling the full length HIV-1 Gag polyprotein reveals the role of its p6 subunit in viral maturation and the effect of non-cleavage site mutations in protease drug resistance

HIV polyprotein Gag is increasingly found to contribute to protease inhibitor resistance. Despite its role in viral maturation and in developing drug resistance, there remain gaps in the knowledge of the role of certain Gag subunits (e.g. p6), and that of non-cleavage mutations in drug resistance. As p6 is flexible, it poses a problem for structural experiments, and is hence often omitted in experimental Gag structural studies. Nonetheless, as p6 is an indispensable component for viral assembly and maturation, we have modeled the full length Gag structure based on several experimentally determined constraints and studied its structural dynamics. Our findings suggest that p6 can mechanistically modulate Gag conformations. In addition, the full length Gag model reveals that allosteric communication between the non-cleavage site mutations and the first Gag cleavage site could possibly result in protease drug resistance, particularly in the absence of mutations in Gag cleavage sites. Our study provides a mechanistic understanding to the structural dynamics of HIV-1 Gag, and also proposes p6 as a possible drug target in anti-HIV therapy.     

Amplification of the effects of drug resistance mutations by background polymorphisms in HIV-1 protease from African subtypes

The vast majority of HIV-1 infections worldwide are caused by the C and A viral subtypes rather than the B subtype prevalent in the United States and Western Europe. Genomic differences between subtypes give rise to sequence variations in the encoded proteins, including those identified as targets for antiretroviral therapies. In the case of the HIV-1 protease, we reported earlier [Velazquez-Campoy et al. (2001) Proc. Natl. Acad. Sci. U.S.A. 98, 6062-6067] that proteases from the C and A subtypes exhibit a higher biochemical fitness in the presence of widely prescribed protease inhibitors. In this paper we present a complete thermodynamic dissection of the differences between proteases from different subtypes and the effects of the V82F/I84V drug-resistant mutation within the framework of the B, C, and A subtypes. These studies involved four inhibitors in clinical use (indinavir, saquinavir, ritonavir, and nelfinavir) and a second-generation protease inhibitor (KNI-764). Naturally occurring amino acid polymorphisms found in proteases from the C and A subtypes lower the binding affinities of existing clinical inhibitors by factors ranging between 2 and 7.5 which by themselves are not enough to cause drug resistance. The preexisting lower affinity in the C and A subtypes, however, significantly amplifies the effects of the drug-resistant mutation. Relative to the wild-type B subtype protease, the V82F/I84V drug-resistant mutation within the C and A subtypes lowers the binding affinity of inhibitors by factors ranging between 40 and 3000. When the enzyme kinetic properties (k(cat) and K(m)) are included in the analysis, the biochemical fitness of the C and A subtype drug-resistant mutants can be up to 1000-fold higher than that of the wild-type B subtype protease in the presence of the studied inhibitors. From a thermodynamic standpoint, the combined effects of the drug-resistant mutations and the natural amino acid polymorphisms on the Gibbs energy are additive and involve significant alterations in the enthalpy and entropy changes associated with inhibitor binding. At the biochemical level, the combined effects of naturally existing polymorphisms and drug-resistant mutations might have important consequences on the long-term viability of current HIV-1 protease inhibitors.     

Drug resistance in HIV-1 protease: Flexibility-assisted mechanism of compensatory mutations

The emergence of drug-resistant variants is a serious side effect associated with acquired immune deficiency syndrome therapies based on inhibition of human immunodeficiency virus type 1 protease (HIV-1 PR). In these variants, compensatory mutations, usually located far from the active site, are able to affect the enzymatic activity via molecular mechanisms that have been related to differences in the conformational flexibility, although the detailed mechanistic aspects have not been clarified so far. Here, we perform multinanosecond molecular dynamics simulations on L63P HIV-1 PR, corresponding to the wild type, and one of its most frequently occurring compensatory mutations, M46I, complexed with the substrate and an enzymatic intermediate. The quality of the calculations is established by comparison with the available nuclear magnetic resonance data. Our calculations indicate that the dynamical fluctuations of the mutated enzyme differ from those in the wild type. These differences in the dynamic properties of the adducts with the substrate and with the gem-diol intermediate might be directly related to variations in the enzymatic activity and therefore offer an explanation of the observed changes in catalytic rate between wild type and mutated enzyme. We anticipate that this "flexibility-assisted" mechanism might be effective in the vast majority of compensatory mutations, which do not change the electrostatic properties of the enzyme.     

Insights into a mutation-assisted lateral drug escape mechanism from the HIV-1 protease active site

We provide insight into the first stages of a kinetic mechanism of lateral drug expulsion from the active site of HIV-1 protease, by conducting all atom molecular dynamics simulations with explicit solvent over a time scale of 24 ns for saquinavir bound to the wildtype, G48V, L90M and G48V/L90M mutant proteases. We find a consistent escape mechanism associated with the G48V mutation. First, increased hydrophilic and hydrophobic flap coupling and water mediated disruption of catalytic dyad hydrogen bonding induce drug motion away from the dyad and promote protease flap transition to the semi-open form. Conversely, flap-inhibitor motion is decoupled in the wildtype. Second, the decrease of total interactions causes unidirectional lateral inhibitor translation by up to 4 A toward the P3 subsite exit of the active site, increased P3 subsite exposure to solvent and a complete loss of hydrophobic interactions with the opposite end of the active site. The P1 subsite moves beyond the hydrophobic active site side pocket, the only remaining steric barrier to complete expulsion being the "breathable" residue, P81. Significant inhibitor deviation is reported over 24 ns, and subsequent complete expulsion, implemented using steered molecular dynamics simulations, is shown to occur most easily for the G48V-containing mutants. Our simulations thus provide compelling support for lateral drug escape from a protease in a semi-open flap conformation. It is likely that some mutations take advantage of this escape mechanism to increase the rate of inhibitor dissociation from the protease. Finally, unidirectional translation may be countered by designing inhibitors with terminal subsites that provide sufficient anchoring to the flaps, thus increasing the steric barrier for translation in either direction.     

Predicting drug resistance of the HIV-1 protease using molecular interaction energy components

Drug resistance significantly impairs the efficacy of AIDS therapy. Therefore, precise prediction of resistant viral mutants is particularly useful for developing effective drugs and designing therapeutic regimen. In this study, we applied a structure-based computational approach to predict mutants of the HIV-1 protease resistant to the seven FDA approved drugs. We analyzed the energetic pattern of the protease-drug interaction by calculating the molecular interaction energy components (MIECs) between the drug and the protease residues. Support vector machines (SVMs) were trained on MIECs to classify protease mutants into resistant and nonresistant categories. The high prediction accuracies for the test sets of cross-validations suggested that the MIECs successfully characterized the interaction interface between drugs and the HIV-1 protease. We conducted a proof-of-concept study on a newly approved drug, darunavir (TMC114), on which no drug resistance data were available in the public domain. Compared with amprenavir, our analysis suggested that darunavir might be more potent to combat drug resistance. To quantitatively estimate binding affinities of drugs and study the contributions of protease residues to causing resistance, linear regression models were trained on MIECs using partial least squares (PLS). The MIEC-PLS models also achieved satisfactory prediction accuracy. Analysis of the fitting coefficients of MIECs in the regression model revealed the important resistance mutations and shed light into understanding the mechanisms of these mutations to cause resistance. Our study demonstrated the advantages of characterizing the protease-drug interaction using MIECs. We believe that MIEC-SVM and MIEC-PLS can help design new agents or combination of therapeutic regimens to counter HIV-1 protease resistant strains.     

A cleavage enzyme-cytometric bead array provides biochemical profiling of resistance mutations in HIV-1 Gag and protease

Most protease-substrate assays rely on short, synthetic peptide substrates consisting of native or modified cleavage sequences. These assays are inadequate for interrogating the contribution of native substrate structure distal to a cleavage site that influences enzymatic cleavage or for inhibitor screening of native substrates. Recent evidence from HIV-1 isolates obtained from individuals resistant to protease inhibitors has demonstrated that mutations distal to or surrounding the protease cleavage sites in the Gag substrate contribute to inhibitor resistance. We have developed a protease-substrate cleavage assay, termed the cleavage enzyme- cytometric bead array (CE-CBA), which relies on native domains of the Gag substrate containing embedded cleavage sites. The Gag substrate is expressed as a fluorescent reporter fusion protein, and substrate cleavage can be followed through the loss of fluorescence utilizing flow cytometry. The CE-CBA allows precise determination of alterations in protease catalytic efficiency (k(cat)/K(M)) imparted by protease inhibitor resistance mutations in protease and/or gag in cleavage or noncleavage site locations in the Gag substrate. We show that the CE-CBA platform can identify HIV-1 protease present in cellular extractions and facilitates the identification of small molecule inhibitors of protease or its substrate Gag. Moreover, the CE-CBA can be readily adapted to any enzyme-substrate pair and can be utilized to rapidly provide assessment of catalytic efficiency as well as systematically screen for inhibitors of enzymatic processing of substrate.     

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.     

Molecular dynamics studies on HIV-1 protease drug resistance and folding pathways

Drug resistance to HIV-1 protease involves the accumulation of multiple mutations in the protein. We investigate the role of these mutations by using molecular dynamics simulations that exploit the influence of the native-state topology in the folding process. Our calculations show that sites contributing to phenotypic resistance of FDA-approved drugs are among the most sensitive positions for the stability of partially folded states and should play a relevant role in the folding process. Furthermore, associations between amino acid sites mutating under drug treatment are shown to be statistically correlated. The striking correlation between clinical data and our calculations suggest a novel approach to the design of drugs tailored to bind regions crucial not only for protein function, but for folding as well.     

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.     

Curling of flap tips in HIV-1 protease as a mechanism for substrate entry and tolerance of drug resistance

BACKGROUND: The human immunodeficiency virus type 1 (HIV-1) protease is an essential viral protein that is a major drug target in the fight against Acquired Immune Deficiency Syndrome (AIDS). Access to the active site of this homodimeric enzyme is gained when two large flaps, one from each monomer, open. The flap movements are therefore central to the function of the enzyme, yet determining how these flaps move at an atomic level has not been experimentally possible. RESULTS: In the present study, we observe the flaps of HIV-1 protease completely opening during a 10 ns solvated molecular dynamics simulation starting from the unliganded crystal structure. This movement is on the time scale observed by Nuclear Magnetic Resonance (NMR) relaxation data. The highly flexible tips of the flaps, with the sequence Gly-Gly-Ile-Gly-Gly, are seen curling back into the protein and thereby burying many hydrophobic residues. CONCLUSIONS: This curled-in conformational change has never been previously described. Previous models of this movement, with the flaps as rigid levers, are not consistent with the experimental data. The residues that participate in this hydrophobic cluster as a result of the conformational change are highly sensitive to mutation and often contribute to drug resistance when they do change. However, several of these residues are not part of the active site cavity, and their essential role in causing drug resistance could possibly be rationalized if this conformational change actually occurs. Trapping HIV-1 protease in this inactive conformation would provide a unique opportunity for future drug design.     

Structural insights into the mechanisms of drug resistance in HIV-1 protease NL4-3

The development of resistance to anti-retroviral drugs targeted against HIV is an increasing clinical problem in the treatment of HIV-1-infected individuals. Many patients develop drug-resistant strains of the virus after treatment with inhibitor cocktails (HAART therapy), which include multiple protease inhibitors. Therefore, it is imperative that we understand the mechanisms by which the viral proteins, in particular HIV-1 protease, develop resistance. We have determined the three-dimensional structure of HIV-1 protease NL4-3 in complex with the potent protease inhibitor TL-3 at 2.0 A resolution. We have also obtained the crystal structures of three mutant forms of NL4-3 protease containing one (V82A), three (V82A, M46I, F53L) and six (V82A, M46I, F53L, V77I, L24I, L63P) point mutations in complex with TL-3. The three protease mutants arose sequentially under ex vivo selective pressure in the presence of TL-3, and exhibit fourfold, 11-fold, and 30-fold resistance to TL-3, respectively. This series of protease crystal structures offers insights into the biochemical and structural mechanisms by which the enzyme can overcome inhibition by TL-3 while recovering some of its native catalytic activity.     

Prevalence and impact of minority variant drug resistance mutations in primary HIV-1 infection

Objective

To evaluate minority variant drug resistance mutations detected by the oligonucleotide ligation assay (OLA) but not consensus sequencing among subjects with primary HIV-1 infection.

Design/Methods

Observational, longitudinal cohort study. Consensus sequencing and OLA were performed on the first available specimens from 99 subjects enrolled after 1996. Survival analyses, adjusted for HIV-1 RNA levels at the start of antiretroviral (ARV) therapy, evaluated the time to virologic suppression (HIV-1 RNA<50 copies/mL) among subjects with minority variants conferring intermediate or high-level resistance.

Results

Consensus sequencing and OLA detected resistance mutations in 5% and 27% of subjects, respectively, in specimens obtained a median of 30 days after infection. Median time to virologic suppression was 110 (IQR 62–147) days for 63 treated subjects without detectable mutations, 84 (IQR 56–109) days for ten subjects with minority variant mutations treated with ≥3 active ARVs, and 104 (IQR 60–162) days for nine subjects with minority variant mutations treated with <3 active ARVs (p = .9). Compared to subjects without mutations, time to virologic suppression was similar for subjects with minority variant mutations treated with ≥3 active ARVs (aHR 1.2, 95% CI 0.6–2.4, p = .6) and subjects with minority variant mutations treated with <3 active ARVs (aHR 1.0, 95% CI 0.4–2.4, p = .9). Two subjects with drug resistance and two subjects without detectable resistance experienced virologic failure.

Conclusions

Consensus sequencing significantly underestimated the prevalence of drug resistance mutations in ARV-naïve subjects with primary HIV-1 infection. Minority variants were not associated with impaired ARV response, possibly due to the small sample size. It is also possible that, with highly-potent ARVs, minority variant mutations may be relevant only at certain critical codons.     

Resistance mechanism revealed by crystal structures of unliganded nelfinavir-resistant HIV-1 protease non-active site mutants N88D and N88S

Nelfinavir is an inhibitor of HIV-1 protease, and is used for treatment of patients suffering from HIV/AIDS. However, treatment results in drug resistant mutations in HIV-1 protease. N88D and N88S are two such mutations which occur in the non-active site region of the enzyme. We have determined crystal structures of unliganded N88D and N88S mutants of HIV-1 protease to resolution of 1.65 Å and 1.8 Å, respectively. These structures refined against synchrotron data lead to R-factors of 0.1859 and 0.1780, respectively. While structural effects of N88D are very subtle, the mutation N88S has caused a significant conformational change in D30, an active site residue crucial for substrate and inhibitor binding.     

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.