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.     

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.     

Prevalence and evolution of drug resistance HIV-1 variants in Henan, China

To understand the prevalence and evolution of drug resistant HIV strains in Henan China after the implementation of free antiretroviral therapy for AIDS patients. 45 drug naive AIDS patients, 118 AIDS patients who received three months antiretroviral therapy and 124 AIDS patients who received six months antiretroviral treatment were recruited in the southern part of Henan province. Information on general condition, antiretroviral medicines, adherence and clinical syndromes were collected by face to face interview. Meanwhile, 14ml EDTA anticoagulant blood was drawn. CD4/CD8 T cell count, viral load and genotypic drug resistance were tested. The rates of clinical improvement were 55.1% and 50.8% respectively three months and six months after antiretroviral therapy. The mean CD4 cell count after antiretroviral therapy was significantly higher than in drug naive patients. The prevalence rate of drug resistant HIV strains were 13.9%, 45.4% and 62.7% in drug naive patients, three month treatment patients and six month treatment patients, respectively. The number of resistance mutation codons and the frequency of mutations increased significantly with continued antiretroviral therapy. The mutation sites were primarily at the 103, 106 and 215 codons in the three-month treatment group and they increased to 15 codon mutations in the six-month treatment group. From this result, the evolution of drug resistant strains was inferred to begin with the high level NNRTI resistant strain, and then develop low level resistant strains to NRTIs. The HIV strains with high level resistance to NVP and low level resistance to AZT and DDI were highly prevalent because of the AZT＋DDI＋NVP combination therapy. These HIV strains were also cross resistant to DLV, EFV, DDC and D4T. Poor adherence to therapy was believed to be the main reason for the emergence and prevalence of drug resistant HIV strains. The prevalence of drug resistant HIV strains was increased with the continuation of antiretroviral therapy in the southern part of Henan province. Measures, that could promote high level adherence, provide new drugs and change ART regimens in failing patients, should be implemented as soon as possible.     

Prevalence and evolution of drug resistance HIV-1 variants in Henan, China

To understand the prevalence and evolution of drug resistant HIV strains in Henan China after the implementation of free antiretroviral therapy for AIDS patients. 45 drug naive AIDS patients, 118 AIDS patients who received three months antiretroviral therapy and 124 AIDS patients who received six months antiretroviral treatment were recruited in the southern part of Henan province. Information on general condition, antiretroviral medicines, adherence and clinical syndromes were collected by face to face interview. Meanwhile, 14ml EDTA anticoagulant blood was drawn. CD4/CD8 T cell count, viral load and genotypic drug resistance were tested. The rates of clinical improvement were 55.1% and 50.8% respectively three months and six months after antiretroviral therapy. The mean CD4 cell count after antiretroviral therapy was significantly higher than in drug naive patients. The prevalence rate of drug resistant HIV strains were 13.9%, 45.4% and 62.7% in drug naive patients, three month treatment patients and six month treatment patients, respectively.The number of resistance mutation codons and the frequency of mutations increased significantly with continued antiretroviral therapy. The mutation sites were primarily at the 103, 106 and 215 codons in the three-month treatment group and they increased to 15 codon mutations in the six-month treatment group. From this result, the evolution of drug resistant strains was inferred to begin with the high level NNRTI resistant strain, and then develop low level resistant strains to NRTIs. The HIV strains with high level resistance to NVP and low level resistance to AZT and DDI were highly prevalent because of the AZT DDI NVP combination therapy. These HIV strains were also cross resistant to DLV, EFV, DDC and D4T. Poor adherence to therapy was believed to be the main reason for the emergence and prevalence of drug resistant HIV strains. The prevalence of drug resistant HIV strains was increased with the continuation of antiretroviral therapy in the southern part of Henan province. Measures, that could promote high level adherence,provide new drugs and change ART regimens in failing patients, should be implemented as soon as possible.     

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     

Binding to the open conformation of HIV-1 protease

A recent crystal structure of HIV-1 protease (HIVp) was the first to experimentally observe a ligand targeting an open-flap conformation. Researchers studying a symmetric pyrrolidine inhibitor found that two ligands cocrystallized with the protease, forcing an unusual configuration and unique crystallographic contacts. One molecule is centered in the traditional binding site (α pose) and the other binds between the flaps (β pose). The ligands stack against each other in a region termed the "eye" site. Ligands bound to the eye site should prevent flap closure, but it is unclear if the pyrrolidine inhibitors or the crystal packing are causing the open state. Molecular dynamics simulations were used to examine the solution-state behavior of three possible binding modes: the ternary complex of HIVp+αβ and the binary complexes, HIVp+α and HIVp+β. We show that HIVp+α is the most stable of the three states. During conformational sampling, α takes an asymmetric binding pose, with one naphthyl ring occupying the eye site and the other reoriented down to occupy positions seen with traditional inhibitors. This finding supports previous studies that reveal a requirement for asymmetric binding at the eye site. In fact, if the α pose is modified to splay both naphthyl rings across the binding site like traditional inhibitors, one ring consistently flips to occupy the eye site. Our simulations reveal that interactions to the eye site encourage a conformationally restrained state, and understanding those contacts may aid the design of ligands to specifically target alternate conformations of the protease.     

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.     

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.     

DynaMut2: Assessing changes in stability and flexibility upon single and multiple point missense mutations

Predicting the effect of missense variations on protein stability and dynamics is important for understanding their role in diseases, and the link between protein structure and function. Approaches to estimate these changes have been proposed, but most only consider single‐point missense variants and a static state of the protein, with those that incorporate dynamics are computationally expensive. Here we present DynaMut2, a web server that combines Normal Mode Analysis (NMA) methods to capture protein motion and our graph‐based signatures to represent the wildtype environment to investigate the effects of single and multiple point mutations on protein stability and dynamics. DynaMut2 was able to accurately predict the effects of missense mutations on protein stability, achieving Pearson's correlation of up to 0.72 (RMSE: 1.02 kcal/mol) on a single point and 0.64 (RMSE: 1.80 kcal/mol) on multiple‐point missense mutations across 10‐fold cross‐validation and independent blind tests. For single‐point mutations, DynaMut2 achieved comparable performance with other methods when predicting variations in Gibbs Free Energy (ΔΔG) and in melting temperature (ΔT_m). We anticipate our tool to be a valuable suite for the study of protein flexibility analysis and the study of the role of variants in disease. DynaMut2 is freely available as a web server and API at http://biosig.unimelb.edu.au/dynamut2 .     

Crystal structure of an in vivo HIV-1 protease mutant in complex with saquinavir: insights into the mechanisms of drug resistance

Saquinavir is a widely used HIV-1 protease inhibitor drug for AIDS therapy. Its effectiveness, however, has been hindered by the emergence of resistant mutations, a common problem for inhibitor drugs that target HIV-1 viral enzymes. Three HIV-1 protease mutant species, G48V, L90M, and G48V/L90M double mutant, are associated in vivo with saquinavir resistance by the enzyme (Jacobsen et al., 1996). Kinetic studies on these mutants demonstrate a 13.5-, 3-, and 419-fold increase in Ki values, respectively, compared to the wild-type enzyme (Ermolieff J, Lin X, Tang J, 1997, Biochemistry 36:12364-12370). To gain an understanding of how these mutations modulate inhibitor binding, we have solved the HIV-1 protease crystal structure of the G48V/L90M double mutant in complex with saquinavir at 2.6 A resolution. This mutant complex is compared with that of the wild-type enzyme bound to the same inhibitor (Krohn A, Redshaw S, Richie JC, Graves BJ, Hatada MH, 1991, J Med Chem 34:3340-3342). Our analysis shows that to accommodate a valine side chain at position 48, the inhibitor moves away from the protease, resulting in the formation of larger gaps between the inhibitor P3 subsite and the flap region of the enzyme. Other subsites also demonstrate reduced inhibitor interaction due to an overall change of inhibitor conformation. The new methionine side chain at position 90 has van der Waals interactions with main-chain atoms of the active site residues resulting in a decrease in the volume and the structural flexibility of S1/S1' substrate binding pockets. Indirect interactions between the mutant methionine side chain and the substrate scissile bond or the isostere part of the inhibitor may differ from those of the wild-type enzyme and therefore may facilitate catalysis by the resistant mutant.     

Oculocutaneous albinism type 1: link between mutations,tyrosinase conformational stability,and enzymatic activity

Oculocutaneous albinism type 1 (OCA1) is an autosomal recessive disorder caused by mutations in the tyrosinase gene. Two subtypes of OCA1 have been described: severe OCA1A with complete absence of tyrosinase activity and less severe OCA1B with residual tyrosinase activity. Here, we characterize the recombinant human tyrosinase intramelanosomal domain and mutant variants, which mimic genetic changes in both subtypes of OCA1 patients. Proteins were prepared using site‐directed mutagenesis, expressed in insect larvae, purified by chromatography, and characterized by enzymatic activities, tryptophan fluorescence, and Gibbs free energy changes. The OCA1A mutants showed very low protein expression and protein yield and are enzymatically inactive. Mutants mimicking OCA1B were biochemically similar to the wild type, but exhibited lower specific activities and protein stabilities. The results are consistent with clinical data, which indicates that OCA1A mutations inactivate tyrosinase and result in severe phenotype, while OCA1B mutations partially inactivate tyrosinase and result in OCA1B albinism.     

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.     

Evolutionary analysis of HIV-1 protease inhibitors: Methods for design of inhibitors that evade resistance

Drug-resistant strains are rapidly selected during AIDS therapy because of the high rate of mutation in HIV. In this report, we present an evolutionary simulation method for analysis of viral mutation and its use for optimization of HIV-1 protease drugs to improve their robustness in the face of resistance mutation. We first present an analysis of the range of resistant mutants that produce viable viruses by using a volume-based viral fitness model. Then, we analyze how this range of mutant proteases allows development of resistance to an optimal inhibitor previously designed by computational coevolution techniques. Finally, we evaluate the resistance patterns of commercially available drugs, and we discuss how resistance might be overcome by optimizing the size of specific side-chains of these inhibitors.     

Study of HIV-1 drug resistance in patients receiving free antiretroviral therapy in China

To investigate the prevalence of drug-resistance mutations, resistance to antiretroviral drugs, and the subsequent virological response to therapy in treatment-naive and antiretroviral-treated patients infected with HIV/AIDS in Henan, China, a total of 431 plasma samples were collected in Queshan county between 2003 and 2004, from patients undergoing the antiretroviral regimen Zidovudine + Didanosine + Nevirapine (Azt+Ddi+Nvp). Personal information was collected by face to face interview. Viral load and genotypic drug resistance were tested. Drug resistance mutation data were obtained by analyzing patient-derived sequences through the HIVdb Program (http://hivdb.stanford.edu). Overall, 38.5% of treatment-naive patients had undetectable plasma viral load (VL), the rate significantly increased to 61.9% in 0 to 6 months treatment patients (mean 3 months) (P<0.005) but again significantly decrease to 38.6% in 6 to 12 months treatment patients (mean 9 months) (P<0.001) and 40.0% in patients receiving more than 12 months treatment (mean 16 months) (P<0.005). The prevalence of drug resistance in patients who had a detectable VL and available sequences were 7.0%, 48.6%, 70.8%, 72.3% in treatment-naïve, 0 to 6 months treatment, 6 to 12 months treatment, and treatment for greater than 12 months patients, respectively. No mutation associated with resistance to Protease inhibitor (PI) was detected in this study. Nucleoside RT inhibitor (NRTI) mutations always emerged after non-nucleoside RT inhibitor (NNRTI) mutations, and were only found in patients treated for more than 6 months, with a frequency less than 5%, with the exception of mutation T215Y (12.8%, 6/47) which occurred in patients treated for more than 12 months. NNRTI mutations emerged quickly after therapy begun, and increased significantly in patients treated for more than 6 months (P<0.005), and the most frequent mutations were K103N, V106A, Y181C, G190A. There had been optimal viral suppression in patients undergoing treatment for less than 6 months in Queshan, Henan. The drug resistance strains were highly prevalent in antiretroviral-treated patients, and increased with the continuation of therapy, with many patients encountering virological failure after 6 months therapy.     

Clarifying allosteric control of flap conformations in the 1TW7 crystal structure of HIV-1 protease

The 1TW7 crystal structure of HIV-1 protease shows the flaps placed wider and more open than what is seen in other examples of the semi-open, apo form. It has been proposed that this might be experimental evidence of allosteric control, because crystal packing creates contacts to the "elbow region" of the protease, which may cause deformation of the flaps. Recent dynamics simulations have shown that the conformation seen in 1TW7 relaxes into the typical semi-open conformation in the absence of the crystal contacts, definitively showing that the crystal contacts cause the deformation (Layten et al., J Am Chem Soc 2006;128:13360-13361). However, this does not prove or disprove allosteric modulation at the elbow. In this study, we have conducted additional simulations, supplemented with experimental testing, to further probe the possibility of 1TW7 providing an example of allosteric control of the flap region. We show that the contacts are unstable and do not restrict the conformational sampling of the flaps. The deformation seen in the 1TW7 crystal structure is simply opportunistic crystal packing and not allosteric control.     

Combining cross-sectional and prospective data methods to improve transition parameter estimation for characterizing the accumulation of HIV-1 drug resistance mutations

The order and rate of acquisition of HIV drug resistance mutations have been estimated from longitudinal and cross-sectional data using Markov models and branching trees, respectively. This article proposes methods that make use of both longitudinal and cross-sectional data simultaneously by employing link functions between the two parameter sets. Most functions that link the two parameter sets also depend on the time on treatment before the start of the study-information that may not be available. Nonetheless, under certain assumptions, some link functions eliminate the dependence on time. Using such functions, the two sources of information can be combined to improve the precision of parameter estimation. The method also accommodates error in the link functions from uncertainty in the assumptions required for the links or other reasons. These methods are applied to data from AIDS Clinical Trial Group protocol 398, a randomized comparison of mono- versus dual-protease inhibitor use in heavily treatment experienced HIV patients. Combining the two sources of information allows detection of differences between rates of transition that are not detectable using prospective data alone.     

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.