Genotypic and phenotypic characterization of human immunodeficiency virus type 1 variants isolated from patients treated with the protease inhibitor nelfinavir

Molecular basis for reduced cleavage activity and drug resistance in D30N HIV-1 protease

Nelfinavir is one of the FDA-approved HIV-1 protease inhibitors and a part of highly active anti-retroviral therapy (HAART) for the treatment of HIV-AIDS. Nelfinavir was the first HIV-1 protease inhibitor to be approved as a paediatric formulation. The application of HAART had resulted in significant improvement in the lives of AIDS patients. However, the emergence of drug resistance in HIV-1 protease has limited the use of many of these drugs including nelfinavir. A unique mutation observed frequently in patients treated with nelfinavir is D30N as it is selected exclusively by nelfinavir. The D30N mutation imparts very high resistance to nelfinavir but unlike other primary mutations does not give cross-resistance to the majority of other drugs. D30N mutation also significantly reduces cleavage activity of HIV-1 protease and affects viral fitness. Here, we have determined crystal structures of D30N HIV-1 protease in unliganded form and in complex with nelfinavir. These structures provide the rationale for reduced cleavage activity and the molecular basis of drug resistance induced by D30N mutation. The loss of coulombic interaction part of a crucial hydrogen bond between the drug and the protease is likely to play a major role in reduced affinity and resistance towards nelfinavir. The decreased catalytic activity of D30N HIV-1 protease due to altered interaction with the substrates and reduced stability of folding core may be the reason for the reduced replicative capacity of the virus harboring mutant HIV-1 protease.Communicated by Ramaswamy H. Sarma.

Viral proteases as therapeutic targets

Some medically important viruses-including retroviruses, flaviviruses, coronaviruses, and herpesviruses-code for a protease, which is indispensable for viral maturation and pathogenesis. Viral protease inhibitors have become an important class of antiviral drugs. Development of the first-in-class viral protease inhibitor saquinavir, which targets HIV protease, started a new era in the treatment of chronic viral diseases. Combining several drugs that target different steps of the viral life cycle enables use of lower doses of individual drugs (and thereby reduction of potential side effects, which frequently occur during long term therapy) and reduces drug-resistance development. Currently, several HIV and HCV protease inhibitors are routinely used in clinical practice. In addition, a drug including an inhibitor of SARS-CoV-2 main protease, nirmatrelvir (co-administered with a pharmacokinetic booster ritonavir as Paxlovid®), was recently authorized for emergency use. This review summarizes the basic features of the proteases of human immunodeficiency virus (HIV), hepatitis C virus (HCV), and SARS-CoV-2 and discusses the properties of their inhibitors in clinical use, as well as development of compounds in the pipeline.

Protease Inhibitors for the Treatment of HIV/AIDS: Recent Advances and Future Challenges

Acquired Immunodeficiency Syndrome (AIDS) is a chronic disease characterized by multiple life-threatening illnesses caused by a retro-virus, Human Immunodeficiency Virus (HIV). HIV infection slowly destroys the immune system and increases the risk of various other infections and diseases. Although, there is no immediate cure for HIV infection/AIDS, several drugs targeting various cruxes of HIV infection are used to slow down the progress of the disease and to boost the immune system. One of the key therapeutic strategies is Highly Active Antiretroviral Therapy (HAART) or ' AIDS cocktail' in a general sense, which is a customized combination of anti-retroviral drugs designed to combat the HIV infection. Since HAART's inception in 1995, this treatment was found to be effective in improving the life expectancy of HIV patients over two decades. Among various classes of HAART treatment regimen, Protease Inhibitors (PIs) are known to be widely used as a major component and found to be effective in treating HIV infection/AIDS. For the past several years, a variety of protease inhibitors have been reported. This review outlines the drug design strategies of PIs, chemical and pharmacological characteristics of some mechanism-based inhibitors, summarizes the recent developments in small molecule based drug discovery with HIV protease as a drug target. Further discussed are the pharmacology, PI drug resistance on HIV PR, adverse effects of HIV PIs and challenges/impediments in the successful application of HIV PIs as an important class of drugs in HAART regimen for the effective treatment of AIDS.

Molecular Determinants of Epistasis in HIV-1 Protease: Elucidating the Interdependence of L89V and L90M Mutations in Resistance

Protease inhibitors have the highest potency among antiviral therapies against HIV-1 infections, yet the virus can evolve resistance. Darunavir (DRV), currently the most potent Food and Drug Administration-approved protease inhibitor, retains potency against single-site mutations. However, complex combinations of mutations can confer resistance to DRV. While the interdependence between mutations within HIV-1 protease is key for inhibitor potency, the molecular mechanisms that underlie this control remain largely unknown. In this study, we investigated the interdependence between the L89V and L90M mutations and their effects on DRV binding. These two mutations have been reported to be positively correlated with one another in HIV-1 patient-derived protease isolates, with the presence of one mutation making the probability of the occurrence of the second mutation more likely. The focus of our investigation is a patient-derived isolate, with 24 mutations that we call "KY"; this variant includes the L89V and L90M mutations. Three additional KY variants with back-mutations, KY(V89L), KY(M90L), and the KY(V89L/M90L) double mutation, were used to experimentally assess the individual and combined effects of these mutations on DRV inhibition and substrate processing. The enzymatic assays revealed that the KY(V89L) variant, with methionine at residue 90, is highly resistant, but its catalytic function is compromised. When a leucine to valine mutation at residue 89 is present simultaneously with the L90M mutation, a rescue of catalytic efficiency is observed. Molecular dynamics simulations of these DRV-bound protease variants reveal how the L90M mutation induces structural changes throughout the enzyme that undermine the binding interactions.

Tipranavir/Ritonavir (500/200 mg and 500/100 mg) Was Virologically Non-Inferior to Lopinavir/Ritonavir (400/100 mg) at Week 48 in Treatment-Naïve HIV-1-Infected Patients: A Randomized, Multinational, Multicenter Trial

Ritonavir-boosted tipranavir (TPV/r) was evaluated as initial therapy in treatment-naïve HIV-1-infected patients because of its potency, unique resistance profile, and high genetic barrier. Trial 1182.33, an open-label, randomized trial, compared two TPV/r dose combinations versus ritonavir-boosted lopinavir (LPV/r). Eligible adults, who had no prior antiretroviral therapy were randomized to twice daily (BID) 500/100 mg TPV/r, 500/200 mg TPV/r, or 400/100 mg LPV/r. Each treatment group also received Tenofovir 300 mg + Lamivudine 300 mg QD. The primary endpoint was a confirmed viral load (VL) <50 copies/mL at week 48 without prior antiretroviral regimen changes. Primary analyses examined CD4-adjusted response rates for non-inferiority, using a 15% non-inferiority margin. At week 48, VL<50 copies/mL was 68.4%, 69.9%, and 72.4% in TPV/r100, TPV/r200, and LPV/r groups, respectively, and TPV/r groups showed non-inferiority to LPV/r. Discontinuation due to adverse events was higher in TPV/r100 (10.3%) and TPV/r200 (15.3%) recipients versus LPV/r (3.2%) recipients. The frequency of grade ≥3 transaminase elevations was higher in the TPV/r200 group than the other groups, leading to closure of this group. However, upon continued treatment or following re-introduction after treatment interruption, transaminase elevations returned to grade ≤2 in >65% of patients receiving either TPV/r200 or TPV/r100. The trial was subsequently discontinued; primary objectives were achieved and continuing TPV/r100 was less tolerable than standard of care for initial highly active antiretroviral therapy. All treatment groups had similar 48-week treatment responses. TPV/r100 and TPV/r200 regimens resulted in sustained treatment responses, which were non-inferior to LPV/r at 48 weeks. When compared with the LPV/r regimen and examined in the light of more current regimens, these TPV/r regimens do not appear to be the best options for treatment-naïve patients based on their safety profiles.

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.

HIV-1 protease-substrate coevolution in nelfinavir resistance

Resistance to various human immunodeficiency virus type 1 (HIV-1) protease inhibitors (PIs) challenges the effectiveness of therapies in treating HIV-1-infected individuals and AIDS patients. The virus accumulates mutations within the protease (PR) that render the PIs less potent. Occasionally, Gag sequences also coevolve with mutations at PR cleavage sites contributing to drug resistance. In this study, we investigated the structural basis of coevolution of the p1-p6 cleavage site with the nelfinavir (NFV) resistance D30N/N88D protease mutations by determining crystal structures of wild-type and NFV-resistant HIV-1 protease in complex with p1-p6 substrate peptide variants with L449F and/or S451N. Alterations of residue 30's interaction with the substrate are compensated by the coevolving L449F and S451N cleavage site mutations. This interdependency in the PR-p1-p6 interactions enhances intermolecular contacts and reinforces the overall fit of the substrate within the substrate envelope, likely enabling coevolution to sustain substrate recognition and cleavage in the presence of PR resistance mutations.

IMPORTANCE

Resistance to human immunodeficiency virus type 1 (HIV-1) protease inhibitors challenges the effectiveness of therapies in treating HIV-1-infected individuals and AIDS patients. Mutations in HIV-1 protease selected under the pressure of protease inhibitors render the inhibitors less potent. Occasionally, Gag sequences also mutate and coevolve with protease, contributing to maintenance of viral fitness and to drug resistance. In this study, we investigated the structural basis of coevolution at the Gag p1-p6 cleavage site with the nelfinavir (NFV) resistance D30N/N88D protease mutations. Our structural analysis reveals the interdependency of protease-substrate interactions and how coevolution may restore substrate recognition and cleavage in the presence of protease drug resistance mutations.

HIV-1 Protease and Substrate Coevolution Validates the Substrate Envelope As the Substrate Recognition Pattern

Drug resistance of HIV-1 protease alters the balance in the molecular recognition events in favor of substrate processing versus inhibitor binding. To develop robust inhibitors targeting ensembles of drug-resistant variants, the code of this balance needs to be cracked. For this purpose, the principles governing the substrate recognition are required to be revealed. Previous crystallographic studies on the WT protease-substrate complexes showed that the substrates have a conserved consensus volume in the protease active site despite their low sequence homology. This consensus volume is termed as the substrate envelope. The substrate envelope was recently reevaluated by taking the substrate dynamics into account, and the dynamic substrate envelope was reported to better define the substrate specificity for HIV-1 protease. Drug resistance occurs mostly through mutations in the protease, occasionally accompanied by cleavage site mutations. In this study, three coevolved protease-substrate complexes (^AP2VNC-p1_V82A, ^LP1'Fp1-p6_D30N/N88D, and ^SP3'Np1-p6_D30N/N88D) were investigated for structural and dynamic properties by molecular modeling and dynamics simulations. The results show the substrate envelope is preserved by these cleavage site mutations in the presence of drug-resistance mutations in the protease, if not enhanced. This study on the conformational and mutational ensembles of protease-substrate complexes validates the substrate envelope as the substrate recognition motif for HIV-1 protease. The substrate envelope hypothesis allows for the elucidation of possible drug resistance mutation patterns in the polyprotein cleavage sites.

A multifaceted analysis of HIV-1 protease multidrug resistance phenotypes

Background

Great strides have been made in the effective treatment of HIV-1 with the development of second-generation protease inhibitors (PIs) that are effective against historically multi-PI-resistant HIV-1 variants. Nevertheless, mutation patterns that confer decreasing susceptibility to available PIs continue to arise within the population. Understanding the phenotypic and genotypic patterns responsible for multi-PI resistance is necessary for developing PIs that are active against clinically-relevant PI-resistant HIV-1 variants.

Results

In this work, we use globally optimal integer programming-based clustering techniques to elucidate multi-PI phenotypic resistance patterns using a data set of 398 HIV-1 protease sequences that have each been phenotyped for susceptibility toward the nine clinically-approved HIV-1 PIs. We validate the information content of the clusters by evaluating their ability to predict the level of decreased susceptibility to each of the available PIs using a cross validation procedure. We demonstrate the finding that as a result of phenotypic cross resistance, the considered clinical HIV-1 protease isolates are confined to ~6% or less of the clinically-relevant phenotypic space. Clustering and feature selection methods are used to find representative sequences and mutations for major resistance phenotypes to elucidate their genotypic signatures. We show that phenotypic similarity does not imply genotypic similarity, that different PI-resistance mutation patterns can give rise to HIV-1 isolates with similar phenotypic profiles.

Conclusion

Rather than characterizing HIV-1 susceptibility toward each PI individually, our study offers a unique perspective on the phenomenon of PI class resistance by uncovering major multidrug-resistant phenotypic patterns and their often diverse genotypic determinants, providing a methodology that can be applied to understand clinically-relevant phenotypic patterns to aid in the design of novel inhibitors that target other rapidly evolving molecular targets as well.

Interplay between single resistance-associated mutations in the HIV-1 protease and viral infectivity, protease activity, and inhibitor sensitivity

Resistance-associated mutations in the HIV-1 protease modify viral fitness through changes in the catalytic activity and altered binding affinity for substrates and inhibitors. In this report, we examine the effects of 31 mutations at 26 amino acid positions in protease to determine their impact on infectivity and protease inhibitor sensitivity. We found that primary resistance mutations individually decrease fitness and generally increase sensitivity to protease inhibitors, indicating that reduced virion-associated protease activity reduces virion infectivity and the reduced level of per virion protease activity is then more easily titrated by a protease inhibitor. Conversely, mutations at more variable positions (compensatory mutations) confer low-level decreases in sensitivity to all protease inhibitors with little effect on infectivity. We found significant differences in the observed effect on infectivity with a pseudotype virus assay that requires the protease to cleave the cytoplasmic tail of the amphotropic murine leukemia virus (MuLV) Env protein. Additionally, we were able to mimic the fitness loss associated with resistance mutations by directly reducing the level of virion-associated protease activity. Virions containing 50% of a D25A mutant protease were 3- to 5-fold more sensitive to protease inhibitors. This level of reduction in protease activity also resulted in a 2-fold increase in sensitivity to nonnucleoside inhibitors of reverse transcriptase and a similar increase in sensitivity to zidovudine (AZT), indicating a pleiotropic effect associated with reduced protease activity. These results highlight the interplay between enzyme activity, viral fitness, and inhibitor mechanism and sensitivity in the closed system of the viral replication complex.

The impact of the nelfinavir resistance-conferring mutation D30N on the susceptibility of HIV-1 subtype B to other protease inhibitors

The human immunodeficiency virus type 1 (HIV-1) protease mutation D30N is exclusively selected by the protease inhibitor (PI) nelfinavir and confers resistance to this drug. We demonstrate that D30N increases the susceptibility to saquinavir (SQV) and amprenavir in HIV-1 subtype B isolates and that the N88D mutation in a D30N background neutralizes this effect. D30N also suppresses indinavir (IDV) resistance caused by the M46I mutation. Interestingly, in patients with viruses originally containing the D30N mutation who were treated with IDV or SQV, the virus either reversed this mutation or acquired N88D, suggesting an antagonistic effect of D30N upon exposure to these PIs. These findings can improve direct salvage drug treatment in resource limited countries where subtype B is epidemiologically important and extend the value of first and second line PIs in these populations.

Efficient identification of human immunodeficiency virus type 1 mutants resistant to a protease inhibitor by using a random mutant library

Emergence of drug-resistant mutant viruses during the course of antiretroviral therapy is a major hurdle that limits the success of chemotherapeutic treatment to suppress human immunodeficiency virus type 1 (HIV-1) replication and AIDS progression. Development of new drugs and careful patient management based on resistance genotyping data are important for enhancing therapeutic efficacy. However, identifying changes leading to drug resistance can take years of clinical studies, and conventional in vitro assays are limited in generating reliable drug resistance data. Here we present an efficient in vitro screening assay for selecting drug-resistant variants from a library of randomly mutated HIV-1 strains generated by transposon-directed base-exchange mutagenesis. As a test of principle, we screened a library of mutant HIV-1 strains containing random mutations in the protease gene by using a reporter T-cell line in the presence of the protease inhibitor (PI) nelfinavir (NFV). Analysis of replicating viruses from a single round of infection identified 50 amino acid substitutions at 35 HIV-1 protease residue positions. The selected mutant viruses showed specific resistance to NFV and included most of the known NFV resistance mutations. Therefore, the new assay is efficient for identifying changes leading to drug resistance. The data also provide insights into the molecular mechanisms underlying the development of drug resistance.

Comparison of drug resistance scores for tipranavir in protease inhibitor-naive patients infected with HIV-1 B and non-B subtypes

Genotypes of samples from protease inhibitor-naïve patients in Frankfurt's HIV Cohort were analyzed with five tipranavir resistance prediction algorithms. Mean scores were higher in non-B than in B subtypes. The proportion of non-B subtypes increased with increasing scores, except in weighted algorithms. Virtual and in vitro phenotype analyses of samples with increased scores showed no reduced tipranavir susceptibility. Current algorithms appear suboptimal for interpretation of resistance to tipranavir in non-B subtypes; increased scores might reflect algorithm bias rather than "natural resistance."

Machine learning on normalized protein sequences

BACKGROUND

Machine learning techniques have been widely applied to biological sequences, e.g. to predict drug resistance in HIV-1 from sequences of drug target proteins and protein functional classes. As deletions and insertions are frequent in biological sequences, a major limitation of current methods is the inability to handle varying sequence lengths.

FINDINGS

We propose to normalize sequences to uniform length. To this end, we tested one linear and four different non-linear interpolation methods for the normalization of sequence lengths of 19 classification datasets. Classification tasks included prediction of HIV-1 drug resistance from drug target sequences and sequence-based prediction of protein function. We applied random forests to the classification of sequences into "positive" and "negative" samples. Statistical tests showed that the linear interpolation outperforms the non-linear interpolation methods in most of the analyzed datasets, while in a few cases non-linear methods had a small but significant advantage. Compared to other published methods, our prediction scheme leads to an improvement in prediction accuracy by up to 14%.

CONCLUSIONS

We found that machine learning on sequences normalized by simple linear interpolation gave better or at least competitive results compared to state-of-the-art procedures, and thus, is a promising alternative to existing methods, especially for protein sequences of variable length.

In Vivo validation of a bioinformatics based tool to identify reduced replication capacity in HIV-1

Although antiretroviral drug resistance is common in treated HIV infected individuals, it is not a consistent indicator of HIV morbidity and mortality. To the contrary, HIV resistance-associated mutations may lead to changes in viral fitness that are beneficial to infected individuals. Using a bioinformatics-based model to assess the effects of numerous drug resistance mutations, we determined that the D30N mutation in HIV-1 protease had the largest decrease in replication capacity among known protease resistance mutations. To test this in silico result in an in vivo environment, we constructed several drug-resistant mutant HIV-1 strains and compared their relative fitness utilizing the SCID-hu mouse model. We found HIV-1 containing the D30N mutation had a significant defect in vivo, showing impaired replication kinetics and a decreased ability to deplete CD4+ thymocytes, compared to the wild-type or virus without the D30N mutation. In comparison, virus containing the M184V mutation in reverse transcriptase, which shows decreased replication capacity in vitro, did not have an effect on viral fitness in vivo. Thus, in this study we have verified an in silico bioinformatics result with a biological assessment to identify a unique mutation in HIV-1 that has a significant fitness defect in vivo.

Molecular Basis for Drug Resistance in HIV-1 Protease

HIV-1 protease is one of the major antiviral targets in the treatment of patients infected with HIV-1. The nine FDA approved HIV-1 protease inhibitors were developed with extensive use of structure-based drug design, thus the atomic details of how the inhibitors bind are well characterized. From this structural understanding the molecular basis for drug resistance in HIV-1 protease can be elucidated. Selected mutations in response to therapy and diversity between clades in HIV-1 protease have altered the shape of the active site, potentially altered the dynamics and even altered the sequence of the cleavage sites in the Gag polyprotein. All of these interdependent changes act in synergy to confer drug resistance while simultaneously maintaining the fitness of the virus. New strategies, such as incorporation of the substrate envelope constraint to design robust inhibitors that incorporate details of HIV-1 protease’s function and decrease the probability of drug resistance, are necessary to continue to effectively target this key protein in HIV-1 life cycle.

Understanding the HIV-1 protease nelfinavir resistance mutation D30N in subtypes B and C through molecular dynamics simulations

A major concern in the antiretroviral (ARV) treatment of HIV infections with protease inhibitors (PI) is the emergence of resistance, which results from the selection of distinct mutations within the viral protease (PR) gene. Among patients who do not respond to treatment with the PI nelfinavir (NFV), the D30N mutation is often observed. However, several reports have shown that D30N emerges with different frequencies in distinct HIV-1 genetic forms or subtypes. In the present work, we analyzed the binding of NFV and the Gag substrate CA/p2 to PR from HIV-1 subtypes B and C through molecular dynamics (MD) simulations. The wild-type and drug-resistant D30N mutants were investigated in both subtypes. The compensatory mutations N83T and N88D, observed in vitro and in vivo when subtype C acquires D30N, were also studied. D30N appears to facilitate conformational changes in subtype B PR, but not in that from subtype C, and this could be associated with disestablishment of an alpha-helical region of the PR. Furthermore, the total contact areas of NFV or the CA/p2 substrate with the mutant PR correlated with changes in the resistance patterns and replicative capacity. Finally, we observed in our MD simulations that mutant PR proteins show different patterns for hydrophobic/van der Waals contact. These findings suggest that different molecular mechanisms contribute to resistance, and we propose that a single mutation has distinct impacts on different HIV-1 subtypes.

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.65A and 1.8A, 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.

Pairwise and higher-order correlations among drug-resistance mutations in HIV-1 subtype B protease

Background

The reaction of HIV protease to inhibitor therapy is characterized by the emergence of complex mutational patterns which confer drug resistance. The response of HIV protease to drugs often involves both primary mutations that directly inhibit the action of the drug, and a host of accessory resistance mutations that may occur far from the active site but may contribute to restoring the fitness or stability of the enzyme. Here we develop a probabilistic approach based on connected information that allows us to study residue, pair level and higher-order correlations within the same framework.

Results

We apply our methodology to a database of approximately 13,000 sequences which have been annotated by the treatment history of the patients from which the samples were obtained. We show that including pair interactions is essential for agreement with the mutational data, since neglect of these interactions results in order-of-magnitude errors in the probabilities of the simultaneous occurence of many mutations. The magnitude of these pair correlations changes dramatically between sequences obtained from patients that were or were not exposed to drugs. Higher-order effects make a contribution of as much as 10% for residues taken three at a time, but increase to more than twice that for 10 to 15-residue groups. The sequence data is insufficient to determine the higher-order effects for larger groups. We find that higher-order interactions have a significant effect on the predicted frequencies of sequences with large numbers of mutations. While relatively rare, such sequences are more prevalent after multi-drug therapy. The relative importance of these higher-order interactions increases with the number of drugs the patient had been exposed to.

Conclusion

Correlations are critical for the understanding of mutation patterns in HIV protease. Pair interactions have substantial qualitative effects, while higher-order interactions are individually smaller but may have a collective effect. Together they lead to correlations which could have an important impact on the dynamics of the evolution of cross-resistance, by allowing the virus to pass through otherwise unlikely mutational states. These findings also indicate that pairwise and possibly higher-order effects should be included in the models of protein evolution, instead of assuming that all residues mutate independently of one another.

Human immunodeficiency virus type 1 protease-correlated cleavage site mutations enhance inhibitor resistance

Drug resistance is an important cause of antiretroviral therapy failure in human immunodeficiency virus (HIV)-infected patients. Mutations in the protease render the virus resistant to protease inhibitors (PIs). Gag cleavage sites also mutate, sometimes correlating with resistance mutations in the protease, but their contribution to resistance has not been systematically analyzed. The present study examines mutations in Gag cleavage sites that associate with protease mutations and the impact of these associations on drug susceptibilities. Significant associations were observed between mutations in the nucleocapsid-p1 (NC-p1) and p1-p6 cleavage sites and various PI resistance-associated mutations in the protease. Several patterns were frequently observed, including mutations in the NC-p1 cleavage site in combination with I50L, V82A, and I84V within the protease and mutations within the p1-p6 cleavage site in combination with D30N, I50V, and I84V within the protease. For most patterns, viruses with mutations both in the protease and in either cleavage site were significantly less susceptible to specific PIs than viruses with mutations in the protease alone. Altered PI resistance in HIV-1 was found to be associated with the presence of Gag cleavage site mutations. These studies suggest that associated cleavage site mutations may contribute to PI susceptibility in highly specific ways depending on the particular combinations of mutations and inhibitors. Thus, cleavage site mutations should be considered when assessing the level of PI resistance.

Molecular basis of human immunodeficiency virus drug resistance: an update

Antiretroviral therapy has led to a significant decrease in human immunodeficiency virus (HIV)-related mortality. Approved antiretroviral drugs target different steps of the viral life cycle including viral entry (coreceptor antagonists and fusion inhibitors), reverse transcription (nucleoside and non-nucleoside inhibitors of the viral reverse transcriptase), integration (integrase inhibitors) and viral maturation (protease inhibitors). Despite the success of combination therapies, the emergence of drug resistance is still a major factor contributing to therapy failure. Viral resistance is caused by mutations in the HIV genome coding for structural changes in the target proteins that can affect the binding or activity of the antiretroviral drugs. This review provides an overview of the molecular mechanisms involved in the acquisition of resistance to currently used and promising investigational drugs, emphasizing the structural role of drug resistance mutations. The optimization of current antiretroviral drug regimens and the development of new drugs are still challenging issues in HIV chemotherapy. This article forms part of a special issue of Antiviral Research marking the 25th anniversary of antiretroviral drug discovery and development, Vol 85, issue 1, 2010.

Molecular characterization of clinical isolates of human immunodeficiency virus resistant to the protease inhibitor darunavir

Darunavir is the most recently approved human immunodeficiency virus (HIV) protease (PR) inhibitor (PI) and is active against many HIV type 1 PR variants resistant to earlier-generation PIs. Darunavir shows a high genetic barrier to resistance development, and virus strains with lower sensitivity to darunavir have a higher number of PI resistance-associated mutations than viruses resistant to other PIs. In this work, we have enzymologically and structurally characterized a number of highly mutated clinically derived PRs with high levels of phenotypic resistance to darunavir. With 18 to 21 amino acid residue changes, the PR variants studied in this work are the most highly mutated HIV PR species ever studied by means of enzyme kinetics and X-ray crystallography. The recombinant proteins showed major defects in substrate binding, while the substrate turnover was less affected. Remarkably, the overall catalytic efficiency of the recombinant PRs (5% that of the wild-type enzyme) is still sufficient to support polyprotein processing and particle maturation in the corresponding viruses. The X-ray structures of drug-resistant PRs complexed with darunavir suggest that the impaired inhibitor binding could be explained by change in the PR-inhibitor hydrogen bond pattern in the P2' binding pocket due to a substantial shift of the aminophenyl moiety of the inhibitor. Recombinant virus phenotypic characterization, enzyme kinetics, and X-ray structural analysis thus help to explain darunavir resistance development in HIV-positive patients.

PL-100, a novel HIV-1 protease inhibitor displaying a high genetic barrier to resistance: an in vitro selection study

The development of new HIV inhibitors with distinct resistance profiles is essential in order to combat the development of multi-resistant viral strains. A drug discovery program based on the identification of compounds that are active against drug-resistant viruses has produced PL-100, a novel potent protease inhibitor (PI) that incorporates a lysine-based scaffold. A selection for resistance against PL-100 in cord blood mononuclear cells was performed, using the laboratory-adapted IIIb strain of HIV-1, and it was shown that resistance appears to develop slower against this compound than against amprenavir, which was studied as a control. Four mutations in protease (PR) were selected after 25 weeks: two flap mutations (K45R and M46I) and two novel active site mutations (T80I and P81S). Site-directed mutagenesis revealed that all four mutations were required to develop low-level resistance to PL-100, which is indicative of the high genetic barrier of the compound. Importantly, these mutations did not cause cross-resistance to currently marketed PIs. In contrast, the P81S mutation alone caused hypersensitivity to two other PIs, saquinavir (SQV) and nelfinavir (NFV). Analysis of p55Gag processing showed that a marked defect in protease activity caused by mutation P81S could only be compensated when K45R and M46I were present. These data correlated well with the replication capacity (RC) of the mutant viruses as measured by a standard viral growth assay, since only viruses containing all four mutations approached the RC of wild type virus. X-ray crystallography provided insight on the structural basis of the resistance conferred by the identified mutations.

Short communication: Different mutation patterns in subtype CRF06_cpx after mother-to-child transmission

Abstract Development of drug resistance mutation patterns (DRMP) in HIV after treatment failure depends on the drugs used in the failing regimen. However, selected patterns may not be unique; there is evidence that selection of DRMP for nelfinavir is dependent on subtype and/or background polymorphisms. Here we describe the selection of DRMP in a mother and son infected with subtype CRF06_cpx by mother-to-child transmission. Four years after delivery the mother received stavudine/lamivudine/nelfinavir as first-line therapy. Genotypic resistance tests (GRT) during follow-up showed selection of M184V/L283I in reverse transcriptase (RT) and H63Q/A71V/L90M in protease (PR). The child started treatment 8 months after birth with stavudine/didanosine/nelfinavir followed by an intensification period with efavirenz. Due to toxicity, efavirenz was removed from the regimen again. GRT during follow-up showed selection of L74V/K103N/M184V/M230L in RT and M46I/H63Q/N88S in PR. The viral load (VL) of the mother was initially undetectable followed by intermediate replication (1000-21,000 copies/ml), whereas the child had both periods of undetectable VL and low-level replication. Although both patients were infected with the same virus and treated with the same protease inhibitor, different DRMPs were selected. Whether the nucleoside backbone, course of antiretroviral therapy, or different host environment is responsible for this variability must be determined in larger studies.

GRL-02031, a novel nonpeptidic protease inhibitor (PI) containing a stereochemically defined fused cyclopentanyltetrahydrofuran potent against multi-PI-resistant human immunodeficiency virus type 1 in vitro

We generated a novel nonpeptidic protease inhibitor (PI), GRL-02031, by incorporating a stereochemically defined fused cyclopentanyltetrahydrofuran (Cp-THF) which exerted potent activity against a wide spectrum of human immunodeficiency virus type 1 (HIV-1) isolates, including multidrug-resistant HIV-1 variants. GRL-02031 was highly potent against laboratory HIV-1 strains and primary clinical isolates, including subtypes A, B, C, and E (50% effective concentration [EC(50)] range, 0.015 to 0.038 microM), with minimal cytotoxicity (50% cytotoxic concentration, >100 microM in CD4(+) MT-2 cells), although it was less active against two HIV-2 strains (HIV-2(EHO) and HIV-2(ROD)) (EC(50), approximately 0.60 microM) than against HIV-1 strains. GRL-02031 at relatively low concentrations blocked the infection and replication of each of the HIV-1(NL4-3) variants exposed to and selected by up to 5 microM of saquinavir, amprenavir, indinavir, nelfinavir, or ritonavir and 1 microM of lopinavir or atazanavir (EC(50) range, 0.036 to 0.14 microM). GRL-02031 was also potent against multi-PI-resistant clinical HIV-1 variants isolated from patients who had no response to the conventional antiretroviral regimens that then existed, with EC(50)s ranging from 0.014 to 0.042 microM (changes in the EC(50)s were less than twofold the EC(50) for wild-type HIV-1). Upon selection of HIV-1(NL4-3) in the presence of GRL-02031, mutants carrying L10F, L33F, M46I, I47V, Q58E, V82I, I84V, and I85V in the protease-encoding region and G62R (within p17), L363M (p24-p2 cleavage site), R409K (within p7), and I437T (p7-p1 cleavage site) in the gag-encoding region emerged. GRL-02031 was potent against a variety of HIV-1(NL4-3)-based molecular infectious clones containing a single primary mutation reported previously or a combination of such mutations, although it was slightly less active against HIV-1 variants containing consecutive amino acid substitutions: M46I and I47V or I84V and I85V. Structural modeling analysis demonstrated a distinct bimodal binding of GRL-02031 to protease, which may provide advantages to GRL-02031 in blocking the replication of a wide spectrum of HIV-1 variants resistant to PIs and in delaying the development of resistance of HIV-1 to GRL-02031. The present data warrant the further development of GRL-02031 as a potential therapeutic agent for the treatment of infections with primary and multidrug-resistant HIV-1 variants.

Enzymatic and structural analysis of the I47A mutation contributing to the reduced susceptibility to HIV protease inhibitor lopinavir

Lopinavir (LPV) is a second-generation HIV protease inhibitor (PI) designed to overcome resistance development in patients undergoing long-term antiviral therapy. The mutation of isoleucine at position 47 of the HIV protease (PR) to alanine is associated with a high level of resistance to LPV. In this study, we show that recombinant PR containing a single I47A substitution has the inhibition constant (K(i) ) value for lopinavir by two orders of magnitude higher than for the wild-type PR. The addition of the I47A substitution to the background of a multiply mutated PR species from an AIDS patient showed a three-order-of-magnitude increase in K(i) in vitro relative to the patient PR without the I47A mutation. The crystal structure of I47A PR in complex with LPV showed the loss of van der Waals interactions in the S2/S2' subsites. This is caused by the loss of three side-chain methyl groups due to the I47A substitution and by structural changes in the A47 main chain that lead to structural changes in the flap antiparallel beta-strand. Furthermore, we analyzed possible interaction of the I47A mutation with secondary mutations V32I and I54V. We show that both mutations in combination with I47A synergistically increase the relative resistance to LPV in vitro. The crystal structure of the I47A/I54V PR double mutant in complex with LPV shows that the I54V mutation leads to a compaction of the flap, and molecular modeling suggests that the introduction of the I54V mutation indirectly affects the strain of the bound inhibitor in the PR binding cleft.

Optimal drug cocktail design: methods for targeting molecular ensembles and insights from theoretical model systems

Drug resistance is a significant obstacle in the effective treatment of diseases with rapidly mutating targets, such as AIDS, malaria, and certain forms of cancer. Such targets are remarkably efficient at exploring the space of functional mutants and at evolving to evade drug binding while still maintaining their biological role. To overcome this challenge, drug regimens must be active against potential target variants. Such a goal may be accomplished by one drug molecule that recognizes multiple variants or by a drug "cocktail"--a small collection of drug molecules that collectively binds all desired variants. Ideally, one wants the smallest cocktail possible due to the potential for increased toxicity with each additional drug. Therefore, the task of designing a regimen for multiple target variants can be framed as an optimization problem--find the smallest collection of molecules that together "covers" the relevant target variants. In this work, we formulate and apply this optimization framework to theoretical model target ensembles. These results are analyzed to develop an understanding of how the physical properties of a target ensemble relate to the properties of the optimal cocktail. We focus on electrostatic variation within target ensembles, as it is one important mechanism by which drug resistance is achieved. Using integer programming, we systematically designed optimal cocktails to cover model target ensembles. We found that certain drug molecules covered much larger regions of target space than others, a phenomenon explained by theory grounded in continuum electrostatics. Molecules within optimal cocktails were often dissimilar, such that each drug was responsible for binding variants with a certain electrostatic property in common. On average, the number of molecules in the optimal cocktails correlated with the number of variants, the differences in the variants' electrostatic properties at the binding interface, and the level of binding affinity required. We also treated cases in which a subset of target variants was to be avoided, modeling the common challenge of closely related host molecules that may be implicated in drug toxicity. Such decoys generally increased the size of the required cocktail and more often resulted in infeasible optimizations. Taken together, this work provides practical optimization methods for the design of drug cocktails and a theoretical, physics-based framework through which useful insights can be achieved.

Target cell APOBEC3C can induce limited G-to-A mutation in HIV-1

The evolutionary success of primate lentiviruses reflects their high capacity to mutate and adapt to new host species, immune responses within individual hosts, and, in recent years, antiviral drugs. APOBEC3G (A3G) and APOBEC3F (A3F) are host cell DNA-editing enzymes that induce extensive HIV-1 mutation that severely attenuates viral replication. The HIV-1 virion infectivity factor (Vif), expressed in vivo, counteracts the antiviral activity of A3G and A3F by inducing their degradation. Other APOBECs may contribute more to viral diversity by inducing less extensive mutations allowing viral replication to persist. Here we show that in APOBEC3C (A3C)-expressing cells infected with the patient-derived HIV-1 molecular clones 210WW, 210WM, 210MW, and 210MM, and the lab-adapted molecular clone LAI, viral G-to-A mutations were detected in the presence of Vif expression. Mutations occurred primarily in the GA context and were relatively infrequent, thereby allowing for spreading infection. The mutations were absent in cells lacking A3C but were induced after transient expression of A3C in the infected target cell. Inhibiting endogenous A3C by RNA interference in Magi cells prevented the viral mutations. Thus, A3C is necessary and sufficient for G-to-A mutations in some HIV-1 strains. A3C-induced mutations occur at levels that allow replication to persist and may therefore contribute to viral diversity. Developing drugs that inhibit A3C may be a novel strategy for delaying viral escape from immune or antiretroviral inhibition.

HIV has shown a chameleon-like nature, always changing to adapt to its environment. Defining the factors that drive and regulate genetic changes in HIV over time is key to understanding how HIV causes disease and escapes from the body's immune responses and drug treatment. The diversity of HIV has implications for the development of effective drugs and for fostering better immune responses. In this study, we showed that a human protein, called APOBEC3C (A3C), could induce certain mutations in some HIV-1 strains, including those derived from a patient who had developed rapid drug resistance. Eliminating the expression of this protein prevented the mutations in two different types of cells. Furthermore, short-term A3C expression was sufficient to cause mutation. We conclude that A3C is necessary and sufficient to induce signature mutations in HIV-1. A3C-induced mutations may provide potential benefit for the virus if the mutation rate is low enough such that the majority of viruses are able to replicate, while accumulating a limited number of novel mutations that may allow the virus to survive in the face of antiviral drugs or immune responses.

The contribution of naturally occurring polymorphisms in altering the biochemical and structural characteristics of HIV-1 subtype C protease

Fourteen subtype B and C protease variants have been engineered in an effort to study whether the preexistent baseline polymorphisms, by themselves or in combination with drug resistance mutations, differentially alter the biochemical and structural features of the subtype C protease when compared with those of subtype B protease. The kinetic studies performed in this work showed that the preexistent polymorphisms in subtype C protease, by themselves, do not provide for a greater level of resistance. Inhibition analysis with eight clinically used protease inhibitors revealed that the natural polymorphisms found in subtype C protease, in combination with drug resistance mutations, can influence enzymatic catalytic efficiency and inhibitor resistance. Structural analyses of the subtype C protease bound to nelfinavir and indinavir showed that these inhibitors form similar interactions with the residues in the active site of subtype B and C proteases. It also revealed that the naturally occurring polymorphisms could alter the position of the outer loops of the subtype C protease, especially the 60's loop.

Human immunodeficiency virus type 1 escape is restricted when conserved genome sequences are targeted by RNA interference

RNA interference (RNAi) is a cellular mechanism in which small interfering RNAs (siRNAs) mediate sequence-specific gene silencing by cleaving the targeted mRNA. RNAi can be used as an antiviral approach to silence the human immunodeficiency virus type 1 (HIV-1) through stable expression of short-hairpin RNAs (shRNAs). We previously reported efficient HIV-1 inhibition by an shRNA against the nonessential nef gene but also described viral escape by mutation or deletion of the nef target sequence. The objective of this study was to obtain insight in the viral escape routes when essential and highly conserved sequences are targeted in the Gag, protease, integrase, and Tat-Rev regions of HIV-1. Target sequences were analyzed of more than 500 escape viruses that were selected in T cells expressing individual shRNAs. Viruses acquired single point mutations, occasionally secondary mutations, but-in contrast to what is observed with nef-no deletions were detected. Mutations occurred predominantly at target positions 6, 8, 9, 14, and 15, whereas none were selected at positions 1, 2, 5, 18, and 19. We also analyzed the type of mismatch in the siRNA-target RNA duplex, and G-U base pairs were frequently selected. These results provide insight into the sequence requirements for optimal RNAi inhibition. This knowledge on RNAi escape may guide the design and selection of shRNAs for the development of an effective RNAi therapy for HIV-1 infections.

HIV-1 protease inhibitors: effects on HIV-2 replication and resistance

Novel antiretroviral drugs include protease (PR) inhibitors (e.g. atazanavir, tipranavir and darunavir) that block HIV-1 maturation and show remarkable antiviral potency on drug-resistant isolates. However, the strains used as prototypes in the design of the novel drugs belong to a specific clade (i.e. HIV-1 group M subtype B), which is the most prevalent in developed countries. At the same time, there is an increasing concern about the expansion of other HIV-1 clades as well as other related retroviruses, such as HIV-2. The HIV-2 PR is weakly inhibited by some PR inhibitors (e.g. amprenavir), and little is known of the mutational pathways leading to drug resistance in this virus. The design of specific PR inhibitors targeting HIV-2, or potent drugs showing broad specificity on HIV-1 and HIV-2 clades, remains a major challenge for the future.

Factors associated with the selection of mutations conferring resistance to protease inhibitors (PIs) in PI-experienced patients displaying treatment failure on darunavir

The objective of this study was to characterize the mutations selected by darunavir (DRV) use in protease inhibitor (PI)-experienced patients and the associated factors. We analyzed treatment failure in 54 PI-experienced human immunodeficiency virus (HIV)-infected patients on a DRV- and ritonavir-containing regimen. Viral genotyping was carried out at the baseline, at between 1 and 3 months of treatment, and at between 3 and 6 months of treatment to search for the selection of mutations conferring resistance to PIs. The median baseline HIV RNA level was 4.9 log(10) copies/ml, and the median CD4 count was 87 cells/mm(3). At the baseline, the median numbers of resistance mutations were as follows: 3 DRV resistance mutations, 4 major PI resistance mutations, and 10 minor PI resistance mutations. The most common mutations that emerged at rebound included V32I (44%), I54M/L (24%), L33F (25%), I84V (21%), and L89V (12%). Multivariate analysis showed that higher baseline HIV RNA levels and smaller numbers of nucleoside reverse transcriptase inhibitor simultaneously used with DRV were associated with a higher risk of DRV resistance mutation selection. By contrast, L76V, a known DRV resistance mutation, was found to decrease the risk of selection of another DRV resistance mutation. The occurrence of virological failure while a patient was on DRV was associated with the selection of mutations that increased the level of DRV resistance without affecting susceptibility to tipranavir (TPV). In these PI-treated patients who displayed treatment failure while they were on a DRV-containing regimen, we confirmed the set of emerging mutations associated with DRV failure and identified the factors associated with the selection of these mutations. TPV susceptibility does not seem to be affected by the selection of a DRV resistance mutation.

An evolutionary-network model reveals stratified interactions in the V3 loop of the HIV-1 envelope

The third variable loop (V3) of the human immunodeficiency virus type 1 (HIV-1) envelope is a principal determinant of antibody neutralization and progression to AIDS. Although it is undoubtedly an important target for vaccine research, extensive genetic variation in V3 remains an obstacle to the development of an effective vaccine. Comparative methods that exploit the abundance of sequence data can detect interactions between residues of rapidly evolving proteins such as the HIV-1 envelope, revealing biological constraints on their variability. However, previous studies have relied implicitly on two biologically unrealistic assumptions: (1) that founder effects in the evolutionary history of the sequences can be ignored, and; (2) that statistical associations between residues occur exclusively in pairs. We show that comparative methods that neglect the evolutionary history of extant sequences are susceptible to a high rate of false positives (20%-40%). Therefore, we propose a new method to detect interactions that relaxes both of these assumptions. First, we reconstruct the evolutionary history of extant sequences by maximum likelihood, shifting focus from extant sequence variation to the underlying substitution events. Second, we analyze the joint distribution of substitution events among positions in the sequence as a Bayesian graphical model, in which each branch in the phylogeny is a unit of observation. We perform extensive validation of our models using both simulations and a control case of known interactions in HIV-1 protease, and apply this method to detect interactions within V3 from a sample of 1,154 HIV-1 envelope sequences. Our method greatly reduces the number of false positives due to founder effects, while capturing several higher-order interactions among V3 residues. By mapping these interactions to a structural model of the V3 loop, we find that the loop is stratified into distinct evolutionary clusters. We extend our model to detect interactions between the V3 and C4 domains of the HIV-1 envelope, and account for the uncertainty in mapping substitutions to the tree with a parametric bootstrap.

Mutational patterns and correlated amino acid substitutions in the HIV-1 protease after virological failure to nelfinavir- and lopinavir/ritonavir-based treatments

Human immunodeficiency virus type 1 (HIV-1) antiviral drug resistance is a major consequence of therapy failure and compromises future therapeutic options. Nelfinavir and lopinavir/ritonavir-based therapies have been widely used in the treatment of HIV-infected patients, in combination with reverse transcriptase inhibitors. The aim of this observational study was the identification and characterization of mutations or combinations of mutations associated with resistance to nelfinavir and lopinavir/ritonavir in treated patients. Nucleotide sequences of 1,515 subtype B HIV-1 isolates from 1,313 persons with different treatment histories (including naïve and treated patients) were collected in 31 Spanish hospitals over the years 2002-2005. Chi-square contingency tests were performed to detect mutations associated with failure to protease inhibitor-based therapies, and correlated mutations were identified using statistical methods. Virological failure to nelfinavir was associated with two different mutational pathways. D30N and N88D appeared mostly in patients without previous exposure to protease inhibitors, while K20T was identified as a secondary resistance mutation in those patients. On the other hand, L90M together with L10I, I54V, A71V, G73S, and V82A were selected in protease inhibitor-experienced patients. A series of correlated mutations including L10I, M46I, I54V, A71V, G73S, and L90M appeared as a common cluster of amino acid substitutions, associated with failure to lopinavir/ritonavir-based treatments. Despite the relatively high genetic barrier of some protease inhibitors, a relatively small cluster of mutations, previously selected under drug pressure, can seriously compromise the efficiency of nelfinavir- and lopinavir/ritonavir-based therapies.

Molecular analysis of the HIV-1 resistance development: enzymatic activities, crystal structures, and thermodynamics of nelfinavir-resistant HIV protease mutants

Human immunodeficiency virus (HIV) encodes an aspartic protease (PR) that cleaves viral polyproteins into mature proteins, thus leading to the formation of infectious particles. Protease inhibitors (PIs) are successful virostatics. However, their efficiency is compromised by antiviral resistance. In the PR sequence of viral variants resistant to the PI nelfinavir, the mutations D30N and L90M appear frequently. However, these two mutations are seldom found together in vivo, suggesting that there are two alternative evolutionary pathways leading to nelfinavir resistance. Here we analyze the proteolytic activities, X-ray structures, and thermodynamics of inhibitor binding to HIV-1 PRs harboring the D30N and L90M mutations alone and in combination with other compensatory mutations. Vitality values obtained for recombinant mutant proteases and selected PR inhibitors confirm the crucial role of mutations in positions 30 and 90 for nelfinavir resistance. The combination of the D30N and L90M mutations significantly increases the enzyme vitality in the presence of nelfinavir, without a dramatic decrease in the catalytic efficiency of the recombinant enzyme. Crystal structures, molecular dynamics simulations, and calorimetric data for four mutants (D30N, D30N/A71V, D30N/N88D, and D30N/L90M) were used to augment our kinetic data. Calorimetric analysis revealed that the entropic contribution to the mutant PR/nelfinavir interaction is less favorable than the entropic contribution to the binding of nelfinavir by wild-type PR. This finding is supported by the structural data and simulations; nelfinavir binds most strongly to the wild-type protease, which has the lowest number of protein-ligand hydrogen bonds and whose structure exhibits the greatest degree of fluctuation upon inhibitor binding.

Clinical significance of human immunodeficiency virus type 1 replication fitness

The relative fitness of a variant, according to population genetics theory, is that variant's relative contribution to successive generations. Most drug-resistant human immunodeficiency virus type 1 (HIV-1) variants have reduced replication fitness, but at least some of these deficits can be compensated for by the accumulation of second-site mutations. HIV-1 replication fitness also appears to influence the likelihood of a drug-resistant mutant emerging during treatment failure and is postulated to influence clinical outcomes. A variety of assays are available to measure HIV-1 replication fitness in cell culture; however, there is no agreement regarding which assays best correlate with clinical outcomes. A major limitation is that there is no high-throughput assay that incorporates an internal reference strain as a control and utilizes intact virus isolates. Some retrospective studies have demonstrated statistically significant correlations between HIV-1 replication fitness and clinical outcomes in some patient populations. However, different studies disagree as to which clinical outcomes are most closely associated with fitness. This may be in part due to assay design, sample size limitations, and differences in patient populations. In addition, the strength of the correlations between fitness and clinical outcomes is modest, suggesting that, at present, it would be difficult to utilize these assays for clinical management.

Computational simulations of HIV-1 proteases--multi-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.

In vitro translation to study HIV protease activity

HIV-1 is an etiological agent of AIDS. One of the targets of the current anti-HIV-1 combination chemotherapy, called highly active antiretroviral therapy (HAART), is HIV-1 protease (PR), which is responsible for the processing of viral structural proteins and, therefore, essential for virus replication. Here, we describe an in vitro transcription/translation-based method of phenotyping HIV-1 PR. In this system, both substrate and PR for the assay can be prepared by in vitro transcription/translation. Protease activity is estimated by the cleavage of a substrate, as measured by enzyme-linked immunosorbent assay (ELISA). This assay is safe, rapid, and requires no special facility to be carried out. Our rapid phenotyping method of HIV-1 PR may help evaluate drug resistance, useful when choosing an appropriate therapeutic regiment, and could potentially facilitate the discovery of new drugs effective against HIV-1 PR.

In vitro antiviral activity and cross-resistance profile of PL-100, a novel protease inhibitor of human immunodeficiency virus type 1

Despite the success of highly active antiretroviral therapy, the current emergence and spread of drug-resistant variants of human immunodeficiency virus (HIV) stress the need for new inhibitors with distinct properties. We designed, produced, and screened a library of compounds based on an original l-lysine scaffold for their potentials as HIV type 1 (HIV-1) protease inhibitors (PI). One candidate compound, PL-100, emerged as a specific and noncytotoxic PI that exhibited potent inhibition of HIV-1 protease and viral replication in vitro (K(i), approximately 36 pM, and 50% effective concentration [EC(50)], approximately 16 nM, respectively). To confirm that PL-100 possessed a favorable resistance profile, we performed a cross-resistance study using a panel of 63 viral strains from PI-experienced patients selected for the presence of primary PI mutations known to confer resistance to multiple PIs now in clinical use. The results showed that PL-100 retained excellent antiviral activity against almost all of these PI-resistant viruses and that its performance in this regard was superior to those of atazanavir, amprenavir, indinavir, lopinavir, nelfinavir, and saquinavir. In almost every case, the increase in the EC(50) for PL-100 observed with viruses containing multiple mutations in protease was far less than that obtained with the other drugs tested. These data underscore the potential for PL-100 to be used in the treatment of drug-resistant HIV disease and argue for its further development.

ART suppresses plasma HIV-1 RNA to a stable set point predicted by pretherapy viremia

Current antiretroviral therapy is effective in suppressing but not eliminating HIV-1 infection. Understanding the source of viral persistence is essential for developing strategies to eradicate HIV-1 infection. We therefore investigated the level of plasma HIV-1 RNA in patients with viremia suppressed to less than 50-75 copies/ml on standard protease inhibitor- or non-nucleoside reverse transcriptase inhibitor-containing antiretroviral therapy using a new, real-time PCR-based assay for HIV-1 RNA with a limit of detection of one copy of HIV-1 RNA. Single copy assay results revealed that >80% of patients on initial antiretroviral therapy for 60 wk had persistent viremia of one copy/ml or more with an overall median of 3.1 copies/ml. The level of viremia correlated with pretherapy plasma HIV-1 RNA but not with the specific treatment regimen. Longitudinal studies revealed no significant decline in the level of viremia between 60 and 110 wk of suppressive antiretroviral therapy. These data suggest that the persistent viremia on current antiretroviral therapy is derived, at least in part, from long-lived cells that are infected prior to initiation of therapy.

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.

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.

HIV-1 protease catalytic efficiency effects caused by random single amino acid substitutions

Protein evolution has occurred by successive fixation of individual mutations. The probability of fixation depends on the fitness of the mutation, and the arising variant can be deleterious, neutral, or beneficial. Despite its relevance, only few studies have estimated the distribution of fitness effects caused by random single mutations on protein function. The human immunodeficiency virus type 1 (HIV-1) protease was chosen as a model protein to quantify protein's tolerability to random single mutations. After determining the enzymatic activity of 107 single random mutants, we found that 86% of single mutations were deleterious for the enzyme catalytic efficiency and 54% lethal. Only 2% of the mutations significantly increased the catalytic efficiency of the enzyme. These data demonstrate the vulnerability of HIV-1 protease to single random mutations. When a second random mutagenesis library was constructed from an HIV-1 protease carrying a highly deleterious single mutation (D30N), a higher proportion of mutations with neutral or beneficial effect were found, 26% and 9%, respectively. Importantly, antagonist epistasis was observed between deleterious mutations. In particular, the mutation N88D, lethal for the wild-type protease, restored the wild-type catalytic efficiency when combined with the highly deleterious mutation D30N. The low tolerability to single random substitutions shown here for the wild-type HIV-1 protease contrasts with its in vivo ability to generate an adaptive variation. Thus, the antagonist epistasis between deleterious or lethal mutations may be responsible for increasing the protein mutational robustness and evolvability.

Analysis of HIV-1 sequences before and after co-infecting syphilis

Increasing syphilis incidence among men who have sex with men (MSM) has been reported. The index case was a human immunodeficiency virus type 1 (HIV-1)-positive MSM who presented coincidentally with the secondary syphilis and a rebound of plasma viral load after complete suppression of HIV-1 (below 50 copies/ml) for 13 months with potent antiretroviral therapy (PART), suggesting a possibility of HIV-1 superinfection. We analyzed HIV-1 sequences before and after syphilis in four HIV-1-positive patients including the index case to explore drug resistance mutations (DRMs) and a possibility of HIV-1 superinfection. There were patients who obtained DRMs around syphilis infection but no evidence of HIV-1 superinfection was obtained. Our results underline the importance of strict adherence to PART.

Silencing of HIV-1 with RNA interference: a multiple shRNA approach

Double-stranded RNA can induce gene silencing via a process known as RNA interference (RNAi). Previously, we have shown that stable expression of a single shRNA targeting the HIV-1 Nef gene strongly inhibits HIV-1 replication. However, this was not sufficient to maintain inhibition. One of the hallmarks of RNAi, its sequence specificity, presented a way out for the virus, as single nucleotide substitutions in the target region abolished inhibition. For the development of a durable gene therapy that prevents viral escape, we proposed to combine multiple shRNAs against conserved HIV-1 regions. Therefore, we screened 86 different shRNAs targeting highly conserved regions. We identified multiple shRNAs that act as potent inhibitors of virus replication. We show, for the first time, that expression of three different shRNAs from a single lentiviral vector results in similar levels of inhibition per shRNA compared to single shRNA vectors. Thus, their combined expression results in a much stronger inhibition of virus production. Moreover, when we infected cells transduced with a double shRNA viral vector, virus escape was delayed. These results confirm that RNAi has great potential as an antiviral gene therapy approach and support our efforts to develop this strategy for treatment of HIV-1-infected individuals.

Evaluation and management of early virological failure

PURPOSE OF REVIEW

To describe the causes, evaluation, and management of early virological failure in patients treated with their first antiretroviral regimen.

RECENT FINDINGS

Combination antiretroviral therapy predictably induces a rapid virological response, with the majority of patients achieving an undetectable HIV-RNA load by week 24. In clinical trials and cohorts, rates of virological suppression have improved over time. Poor adherence to therapy remains the most common cause of virological failure, and genotype resistance testing is a critical step in evaluating the optimal subsequent approach. Most studies suggest that transient HIV-RNA elevations do not warrant changing treatment, and may be a consequence of laboratory variation. For those who experience virological failure, resistance to individual components of the antiretroviral regimen is dependent on the initial choice of treatment. Once virological failure is confirmed and adherence issues addressed, a prompt change in treatment is warranted to limit the selection of further drug resistance.

SUMMARY

Resistance that occurs with early virological failure follows typical patterns, with limited resistance to most antiretroviral agents in the nucleoside reverse transcriptase and protease inhibitor classes. The likelihood of achieving virological suppression with subsequent regimens should be high so long as adherence can be assured.