Characterization and structural analysis of novel mutations in human immunodeficiency virus type 1 reverse transcriptase involved in the regulation of resistance to nonnucleoside inhibitors

The Phenomenon of Antiretroviral Drug Resistance in the Context of Human Immunodeficiency Virus Treatment: Dynamic and Ever Evolving Subject Matter

Human immunodeficiency virus (HIV) is a significant global health issue that affects a substantial number of individuals across the globe, with a total of 39 million individuals living with HIV/AIDS. ART has resulted in a reduction in HIV-related mortality. Nevertheless, the issue of medication resistance is a significant obstacle in the management of HIV/AIDS. The unique genetic composition of HIV enables it to undergo rapid mutations and adapt, leading to the emergence of drug-resistant forms. The development of drug resistance can be attributed to various circumstances, including noncompliance with treatment regimens, insufficient dosage, interactions between drugs, viral mutations, preexposure prophylactics, and transmission from mother to child. It is therefore essential to comprehend the molecular components of HIV and the mechanisms of antiretroviral medications to devise efficacious treatment options for HIV/AIDS.

APOBEC3 selects V179I in HIV-1 reverse transcriptase to provide selective advantage for non-nucleoside reverse transcriptase inhibitor-resistant mutants

The V179I substitution in human immunodeficiency virus type 1 (HIV-1) reverse transcriptase (RT) is selected in humans or mouse models treated with certain nonnucleoside reverse transcriptase inhibitors (NNRTIs). While it is often observed together with other NNRTI resistance mutations, V179I does not confer drug resistance. To understand how V179I arises during NNRTI treatment, we characterized it in HIV-1 molecular clones with or without the NNRTI resistance mutations Y181C or Y181V. While V179I alone did not confer resistance to any NNRTIs tested, when present with Y181C/V it enhanced drug resistance to some NNRTIs by 3- to 8-fold. In replication competition experiments in the presence of the NNRTI rilpivirine (RPV), V179I modestly enhanced Y181C HIV-1 or Y181V HIV-1 replication compared to viruses without V179I. As V179I arises from a G to A mutation, we evaluated whether it could arise due to host APOBEC3 deaminase activity and be maintained in the presence of a NNRTI to provide a selective advantage for the virus. V179I was detected in some humanized mice treated with RPV and was associated with G to A mutations characteristic of APOBEC3 activity. In RPV selection experiments, the frequency of V179I in HIV-1 was accelerated in CD4+ T cells expressing higher APOBEC3F and APOBEC3G levels. Our results provide evidence that V179I in HIV-1 RT can arise due to APOBEC-mediated G to A hypermutation and can confer a selective advantage to drug-resistant HIV-1 isolates in the presence of some NNRTIs.

Evaluation of HIV-1 integrase variability by combining computational and probabilistic approaches

This study aimed at updating previous data on HIV-1 integrase variability, by using effective bioinformatics methods combining different statistical instruments from simple entropy and mutation rate to more specific approaches such as Hellinger distance. A total of 2133 HIV-1 integrase sequences were analyzed in: i) 1460 samples from drug-naïve [DN] individuals; ii) 386 samples from drug-experienced but INI-naïve [IN] individuals; iii) 287 samples from INI-experienced [IE] individuals. Within the three groups, 76 amino acid positions were highly conserved (≤0.2% variation, Hellinger distance: <0.25%), with 35 fully invariant positions; while, 80 positions were conserved (>0.2% to <1% variation, Hellinger distance: <1%). The H12-H16-C40-C43 and D64-D116-E152 motifs were all well conserved. Some residues were affected by dramatic changes in their mutation distributions, especially between DN and IE samples (Hellinger distance ≥1%). In particular, 15 positions (D6, S24, V31, S39, L74, A91, S119, T122, T124, T125, V126, K160, N222, S230, C280) showed a significant decrease of mutation rate in IN and/or IE samples compared to DN samples. Conversely, 8 positions showed significantly higher mutation rate in samples from treated individuals (IN and/or IE) compared to DN. Some of these positions, such as E92, T97, G140, Y143, Q148 and N155, were already known to be associated with resistance to integrase inhibitors; other positions including S24, M154, V165 and D270 are not yet documented to be associated with resistance. Our study confirms the high conservation of HIV-1 integrase and identified highly invariant positions using robust and innovative methods. The role of novel mutations located in the critical region of HIV-1 integrase deserves further investigation.

HIV-1 drug resistance profiling using amino acid sequence space cartography

MOTIVATION

Human immunodeficiency virus (HIV) drug resistance is a global healthcare issue. The emergence of drug resistance influenced the efficacy of treatment regimens, thus stressing the importance of treatment adaptation. Computational methods predicting the drug resistance profile from genomic data of HIV isolates are advantageous for monitoring drug resistance in patients. However, existing computational methods for drug resistance prediction are either not suitable for emerging HIV strains with complex mutational patterns or lack interpretability, which is of paramount importance in clinical practice. The approach reported here overcomes these limitations and combines high accuracy of predictions and interpretability of the models.

RESULTS

In this work, a new methodology based on generative topographic mapping (GTM) for biological sequence space representation and quantitative genotype-phenotype relationships prediction purposes was introduced. The GTM-based resistance landscapes allowed us to predict the resistance of HIV strains based on sequencing and drug resistance data for three viral proteins [integrase (IN), protease (PR) and reverse transcriptase (RT)] from Stanford HIV drug resistance database. The average balanced accuracy for PR inhibitors was 0.89 ± 0.01, for IN inhibitors 0.85 ± 0.01, for non-nucleoside RT inhibitors 0.73 ± 0.01 and for nucleoside RT inhibitors 0.84 ± 0.01. We have demonstrated in several case studies that GTM-based resistance landscapes are useful for visualization and analysis of sequence space as well as for treatment optimization purposes. Here, GTMs were applied for the in-depth analysis of the relationships between mutation pattern and drug resistance using mutation landscapes. This allowed us to predict retrospectively the importance of the presence of particular mutations (e.g. V32I, L10F and L33F in HIV PR) for the resistance development. This study highlights some perspectives of GTM applications in clinical informatics and particularly in the field of sequence space exploration.

AVAILABILITY AND IMPLEMENTATION

https://github.com/karinapikalyova/ISIDASeq.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Detection of Gag C-terminal mutations among HIV-1 non-B subtypes in a subset of Cameroonian patients

Response to ritonavir-boosted-protease inhibitors (PI/r)-based regimen is associated with some Gag mutations among HIV-1 B-clade. There is limited data on Gag mutations and their covariation with mutations in protease among HIV-1 non-B-clades at PI/r-based treatment failure. Thus, we characterized Gag mutations present in isolates from HIV-1 infected individuals treated with a PI/r-regimen (n = 143) and compared them with those obtained from individuals not treated with PI/r (ART-naïve [n = 101] or reverse transcriptase inhibitors (RTI) treated [n = 118]). The most frequent HIV-1 subtypes were CRF02_AG (54.69%), A (13.53%), D (6.35%) and G (4.69%). Eighteen Gag mutations showed a significantly higher prevalence in PI/r-treated isolates compared to ART-naïve (p < 0.05): Group 1 (prevalence < 1% in drug-naïve): L449F, D480N, L483Q, Y484P, T487V; group 2 (prevalence 1–5% in drug-naïve): S462L, I479G, I479K, D480E; group 3 (prevalence ≥ 5% in drug-naïve): P453L, E460A, R464G, S465F, V467E, Q474P, I479R, E482G, T487A. Five Gag mutations (L449F, P453L, D480E, S465F, Y484P) positively correlated (Phi ≥ 0.2, p < 0.05) with protease-resistance mutations. At PI/r-failure, no significant difference was observed between patients with and without these associated Gag mutations in term of viremia or CD4 count. This analysis suggests that some Gag mutations show an increased frequency in patients failing PIs among HIV-1 non-B clades.

Long-Acting Rilpivirine (RPV) Preexposure Prophylaxis Does Not Inhibit Vaginal Transmission of RPV-Resistant HIV-1 or Select for High-Frequency Drug Resistance in Humanized Mice

The antiretroviral drug rilpivirine was developed into a long-acting formulation (RPV LA) to improve adherence for preexposure prophylaxis (PrEP) to prevent HIV-1 transmission. A concern is that RPV LA will not inhibit transmission of drug-resistant HIV-1 and may select for drug-resistant virus. In female humanized mice, we found that RPV LA inhibited vaginal transmission of WT or 3-fold RPV-resistant HIV-1 but not virus with 30-fold RPV resistance. In animals that became infected despite RPV LA PrEP, WT HIV-1 dissemination was delayed until genital and plasma RPV concentrations waned. RPV resistance was detected at similar low frequencies in untreated and PrEP-treated mice that became infected. These results indicate the importance of maintaining RPV at a sustained threshold after virus exposure to prevent dissemination of HIV-1 after vaginal infection and low-frequency resistance mutations conferred low-level resistance, suggesting that RPV resistance is difficult to develop after HIV-1 infection during RPV LA PrEP.

As a long-acting formulation of the nonnucleoside reverse transcriptase inhibitor rilpivirine (RPV LA) has been proposed for use as preexposure prophylaxis (PrEP) and the prevalence of transmitted RPV-resistant viruses can be relatively high, we evaluated the efficacy of RPV LA to inhibit vaginal transmission of RPV-resistant HIV-1 in humanized mice. Vaginal challenges of wild-type (WT), Y181C, and Y181V HIV-1 were performed in mice left untreated or after RPV PrEP. Plasma viremia was measured for 7 to 10 weeks, and single-genome sequencing was performed on plasma HIV-1 RNA in mice infected during PrEP. RPV LA significantly prevented vaginal transmission of WT HIV-1 and Y181C HIV-1, which is 3-fold resistant to RPV. However, it did not prevent transmission of Y181V HIV-1, which has 30-fold RPV resistance in the viruses used for this study. RPV LA did delay WT HIV-1 dissemination in infected animals until genital and plasma RPV concentrations waned. Animals that became infected despite RPV LA PrEP did not acquire new RPV-resistant mutations above frequencies in untreated mice or untreated people living with HIV-1, and the mutations detected conferred low-level resistance. These data suggest that high, sustained concentrations of RPV were required to inhibit vaginal transmission of HIV-1 with little or no resistance to RPV but could not inhibit virus with high resistance. HIV-1 did not develop high-level or high-frequency RPV resistance in the majority of mice infected after RPV LA treatment. However, the impact of low-frequency RPV resistance on virologic outcome during subsequent antiretroviral therapy still is unclear.

IMPORTANCE The antiretroviral drug rilpivirine was developed into a long-acting formulation (RPV LA) to improve adherence for preexposure prophylaxis (PrEP) to prevent HIV-1 transmission. A concern is that RPV LA will not inhibit transmission of drug-resistant HIV-1 and may select for drug-resistant virus. In female humanized mice, we found that RPV LA inhibited vaginal transmission of WT or 3-fold RPV-resistant HIV-1 but not virus with 30-fold RPV resistance. In animals that became infected despite RPV LA PrEP, WT HIV-1 dissemination was delayed until genital and plasma RPV concentrations waned. RPV resistance was detected at similar low frequencies in untreated and PrEP-treated mice that became infected. These results indicate the importance of maintaining RPV at a sustained threshold after virus exposure to prevent dissemination of HIV-1 after vaginal infection and low-frequency resistance mutations conferred low-level resistance, suggesting that RPV resistance is difficult to develop after HIV-1 infection during RPV LA PrEP.

Identification of novel bifunctional HIV-1 reverse transcriptase inhibitors

Objectives

The increasing prevalence of mutations in HIV-1 reverse transcriptase (RT) that confer resistance to existing NRTIs and NNRTIs underscores the need to develop RT inhibitors with novel mode-of-inhibition and distinct resistance profiles.

Methods

Biochemical assays were employed to identify inhibitors of RT activity and characterize their mode of inhibition. The antiviral activity of the inhibitors was assessed by cell-based assays using laboratory HIV-1 isolates and MT4 cells. RT variants were purified via avidin affinity columns.

Results

Compound A displayed equal or greater potency against many common NNRTI-resistant RTs (K103N and Y181C RTs) relative to WT RT. Despite possessing certain NNRTI-like properties, such as being unable to inhibit an engineered variant of RT lacking an NNRTI-binding pocket, we found that compound A was dependent on Mg2+ for binding to RT. Optimization of compound A led to more potent analogues, which retained similar activities against WT and K103N mutant viruses with submicromolar potency in a cell-based assay. One of the analogues, compound G, was crystallized in complex with RT and the structure was determined at 2.6 Å resolution. The structure indicated that compound G simultaneously interacts with the active site (Asp186), the highly conserved primer grip region (Leu234 and Trp229) and the NNRTI-binding pocket (Tyr188).

Conclusions

These findings reveal a novel class of RT bifunctional inhibitors that are not sensitive to the most common RT mutations, which can be further developed to address the deficiency of current RT inhibitors.

HIV Reverse Transcriptase and Protease Genes Variability Can Be a Biomarker Associated with HIV and Hepatitis B or C Coinfection

Variability of the HIV reverse transcriptase (RT) and protease (PR) genes has been used as indicators of drug resistance and as a mean to evaluate phylogenetic relationships among circulating virus. However, these studies have been carried in HIV mono-infected populations. The goal of this study was to evaluate, for the first time, the HIV PR and RT sequences from HIV/HBV and HIV/HCV co-infected patients. HIV PR and RT genes were amplificated and sequenced to resistance analysis. The bioinformatics analysis was performed to infer about sequences clustering and molecular evolution. The results showed that the most frequent amino acid substitutions in RT were L214F (67.6%), I135T (55.9%), and in PR was V15I (41.2%). The molecular clock analysis showed that the HIV circulating in co-infected patients were separated in two clusters in the years 1999-2000. Some patients included as HIV mono-infected according patients' medical records and inside the co-infected cluster were, in fact, co-infected by PCR analysis. Analysis of the decision trees showed susceptibility to lamivudine and emtricitabine were important attribute to characterize co-infected patients. In conclusion, the results obtained in this study suggest, for the first time, that HIV RT and PR genes variability could be a genetic biomarker to coinfection.

Polymorphisms and Mutational Covariation Associated with Death in a Prospective Cohort of HIV/AIDS Patients Receiving Long-Term ART in China

Background

HIV drug resistance is associated with faster clinical progression of AIDS. However, the effect of significant polymorphisms and mutational covariation on mortality among HIV/AIDS patients receiving long-term antiretroviral therapy (ART), have rarely been studied.

Methods

In this prospective cohort study from December 2003 to December 2014, we present a new computational modelling approach based on bioinformatics-based models and several statistical methods to elucidate the molecular mechanisms involved in the acquisition of polymorphisms and mutations on death in HIV/AIDS patients receiving long-term ART in China.

Results

This study involved 654 ART-treated patients, who had been followed for 5473.4 person-years, a median of 9.8 years, and 178 died (25.2%, 3.3/100 person-years). The first regimens included AZT/d4T + NVP+ ddI (78.9%) or AZT/d4T + NVP+ 3TC (20.0%). We calculated an individual Ka/Ks value for each specific amino acid mutation. Result showed that 20 polymorphisms (E6D, Q18H, E35D, S37N, T39A, K43E, S68N, L74I, I93L, K103N, V106A, E169D, Y181C, G190A, Q197K, T200V, T200E, T215I, E224D and P225H) were strongly associated with AIDS related deaths. Among them, 7 polymorphisms (L74I, K103N, V106A, Y181C, G190A, T215I and P225H) were known to be drug resistance mutations, 7 polymorphisms (E6D, E35D, S37N, I93L, E169D, T200V and T200E were considered to be potential drug resistance mutations, and 6 polymorphisms (T39A, K43E, S68N, Q197K, T200V and E224D) were newly found to have an association with drug resistance mutations, which formed a complex network of relationships.

Conclusions

Some polymorphisms and mutational covariation may be the important influencing factors in the failure of treatment. Understanding these mechanisms is essential for the development of new therapies, designing optimal drug combinations, and determining effective clinical management of individual patients.

Novel Mutations L228I and Y232H Cause Nonnucleoside Reverse Transcriptase Inhibitor Resistance in Combinational Pattern

The emergence of drug resistance mutations is increasing after the implementation of highly active antiretroviral therapy. To characterize two novel mutations L228I and Y232H in the primer grip of reverse transcriptase (RT) of HIV-1 circulating recombination form 08_BC (CRF08_BC) subtype, both mutant clones were constructed to determine their impacts on viral phenotypic susceptibility and replication capacity (RC). Results showed that the novel mutation, L228I, conferred a low-level resistance to etravirine by itself. L228I in combination with Y188C displayed a high level of cross-resistance to both nevirapine (NVP) and efavirenz (EFV). The copresence of A139V and Y232H induced a moderate level of resistance to NVP and EFV. Mutations Y188C/L228I, A139V, Y232H, and A139V/Y232H reduced more than 55% of viral RC compared with that of the wild-type (WT) reference virus. Modeling study suggested that the copresence of Y188C/L228I or A139V/Y232H might induce conformational changes to RT, which might result in reduced drug susceptibility and viral RC due to abolished hydrogen bonding or complex interaction with vicinal residues. Our results demonstrated that L228I and Y232H were novel accessory nonnucleoside reverse transcriptase inhibitor resistance-related mutations and provided valuable information for clinicians to design more effective treatment to patients infected with HIV-1 subtype CRF08_BC.

Rilpivirine resistance mutation E138K in HIV-1 reverse transcriptase predisposed by prevalent polymorphic mutations

BACKGROUND

Rilpivirine is listed as a recommended or alternative key drug in the current ART guidelines. E138K in HIV-1 reverse transcriptase (RT) is a primary mutation in resistance to rilpivirine, although in vitro experiments showed it confers only <3-fold resistance. An unidentified mechanism could amplify resistance to rilpivirine conferred by E138K.

OBJECTIVES

The objective of this study was to reveal the mechanism amplifying rilpivirine resistance conferred by E138K.

PATIENTS AND METHODS

HIV-1 RT sequences were compared in patients who failed rilpivirine-containing ART virologically. The effects of mutations commonly identified with E138K on rilpivirine susceptibility were analysed by using recombinant HIV-1 variants.

RESULTS

Rilpivirine-containing ART was introduced in 162 HIV-1-infected patients at the outpatient clinic of the AIDS Clinical Center (National Center for Global Health and Medicine, Tokyo, Japan) between May 2012 and June 2015. Virological treatment failure occurred in six of these patients. E138K emerged in three patients while other rilpivirine resistance mutations emerged in the other three patients. I135T/L were identified in only three patients with E138K and existed before the introduction of rilpivirine-containing ART. Analysis of recombinant HIV-1 variants indicated that E138K conferred low-level rilpivirine resistance and that coexistence of I135T/L with E138K amplified the resistance.

CONCLUSIONS

I135T/L, escape mutations from HLA-B*51/52-restricted cytotoxic T lymphocytes, which are prevalent in Japan, may predispose HIV-1 to harbour E138K upon failure of rilpivirine-containing ART. The mutation patterns of drug resistance may vary due to baseline polymorphic mutations.

Characterization of two HIV-1 infectors during initial antiretroviral treatment, and the emergence of phenotypic resistance in reverse transcriptase-associated mutation patterns

BACKGROUND

Highly active antiretroviral therapy (HAART) is recommended to control the infection of HIV-1. HIV-1 drug resistance becomes an obstacle to HAART due to the accumulation of specific mutations in the RT coding region. The development of resistance mutations may be more complex than previously thought.

METHODS

We followed two HIV-1 infectors from a HIV-1 drug resistance surveillance cohort in Henan province and evaluated CD4+ T-cell number and viral load thereafter at ten time-periods and characterized their reverse transcriptase-associated mutation patterns at each time point. Then we constructed the recombinant virus strains with these mutation patterns to mimick the viruses and test the phenotypic resistance caused by the mutation patterns on TZM-b1 cells.

RESULTS

CD4+ T-cell number initially increased and then decreased rapidly, while viral load decreased and then dropped sharply during initial antiretroviral treatment. The number of mutations and the combination patterns of mutations increased over time. According to the phenotypic resistance performed by recombinant virus strains, VirusT215Y/V179E/Y181C/H221Y exhibited high levels of resistance to EFV (5.57-fold), and T215Y/V179E-containing virus increased 20.20-fold in AZT resistance (p < 0.01). VirusT215Y/V179E/Y181C increased markedly in EFV resistance (p < 0.01). The IC50 for VirusT215Y/V179E/H221Y was similar to that for VirusT215Y/V179E/Y181C. VirusT215Y/K103N/Y181C/H221Y induced a dramatic IC50 increase of all the four agents (Efavirenz EFV, Zidovudine AZT, Lamivudine 3TC, and Stavudine d4T) (p < 0.01). As for VirusT215Y/K103N/Y181C, only the IC50 of EFV was significantly increased. T215Y/K103N resulted in a 26.36-fold increase in EFV (p < 0.01). T215Y/K103N/H221Y significantly increased the resistance to AZT and 3TC. The IC50 of EFV with T215Y/V179E was lower than with T215Y/K103N (F = 93.10, P < 0.0001). With T215Y/V179E, Y181C significantly increase in EFV resistance, while the interaction between 181 and 221 in EFV was not statistically significant (F = 1.20, P = 0.3052). With T215Y/K103N, neither H221Y nor Y181C showed a significant increase in EFV resistance, but the interaction between 181 and 221 was statistically significant (F = 38.12, P = 0.0003).

CONCLUSIONS

Data in this study suggests that pathways of viral evolution toward drug resistance appear to proceed through distinct steps and at different rates. Phenotypic resistance using recombinant virus strains with different combination of mutation patterns reveals that interactions among mutations may provide information on the impact of these mutations on drug resistance. All the result provides reference to optimize clinical treatment schedule.

A novel mutation, D404N, in the connection subdomain of reverse transcriptase of HIV-1 CRF08_BC subtype confers cross-resistance to NNRTIs

OBJECTIVES

Growing evidence suggests that mutations in the connection domain of the HIV-1 reverse transcriptase (RT) can contribute to viral resistance to RT inhibitors. This work was designed to determine the effects of a novel mutation, D404N, in the connection subdomain of RT of HIV-1 CRF08_BC subtype on drug resistance, viral replication capacity (RC) and RT activity.

METHODS

Mutation D404N, alone or together with the other reported mutations, was introduced into an HIV-1 CRF08_BC subtype infectious clone by site-directed mutagenesis. Viral susceptibility to nine RT inhibitors, viral RC and the DNA polymerase activity of viral RT of the constructed virus mutants were investigated. A modelling study using the server SWISS-MODEL was conducted to explore the possible structure-related drug resistance mechanism of the mutation D404N.

RESULTS

Single mutations D404N and H221Y conferred low-level resistance to nevirapine, efavirenz, rilpivirine and zidovudine. Double mutations Y181C/D404N and Y181C/H221Y significantly reduced susceptibility to NNRTIs. The most pronounced resistance to NNRTIs was observed with the triple mutation Y181C/D404N/H221Y. Virus containing D404N as the only mutation displayed ∼50% RC compared with the WT virus. The modelling study suggested that the D404N mutation might abolish the hydrogen bonds between residues 404 and K30 in p51 or K431 in p66, leading to impaired RT subunit structure and enhanced drug resistance.

CONCLUSIONS

These results indicate that D404N is a novel NNRTI-associated mutation in the HIV-1 subtype CRF08_BC and provides information valuable for the monitoring of clinical RTI resistance.

HIV-1 transmitted drug resistance-associated mutations and mutation co-variation in HIV-1 treatment-naïve MSM from 2011 to 2013 in Beijing, China

Background

Transmitted drug resistance (TDR) is an important public health issue, because TDR-associated mutation may affect the outcome of antiretroviral treatment potentially or directly. Men who have sex with men (MSM) constitute a major risk group for HIV transmission. However, current reports are scarce on HIV TDR-associated mutations and their co-variation among MSM.

Methods

Blood samples from 262 newly diagnosed HIV-positive, antiretroviral therapy (ART)-naïve MSM, were collected from January 2011 and December 2013 in Beijing. The polymerase viral genes were sequenced to explore TDR-associated mutations and mutation co-variation.

Results

A total of 223 samples were sequenced and analyzed. Among them, HIV-1 CRF01_AE are accounted for 60.5%, followed by CRF07_BC (27.8%), subtype B (9.9%), and others. Fifty-seven samples had at least one TDR-associated mutation, mainly including L10I/V (6.3%), A71L/T/V (6.3%), V179D/E (5.4%), and V106I (2.7%), with different distributions of TDR-associated mutations by different HIV-1 subtypes and by each year. Moreover, eight significant co-variation pairs were found between TDR-associated mutations (V179D/E) and seven overlapping polymorphisms in subtype CRF01_AE.

Conclusions

To date, this work consists the most comprehensive genetic characterization of HIV-1 TDR-associated mutations prevalent among MSM. It provides important information for understanding TDR and viral evolution among Chinese MSM, a population currently at particularly high risk of HIV transmission.

Electronic supplementary material

The online version of this article (doi:10.1186/s12879-014-0689-7) contains supplementary material, which is available to authorized users.

The development of drug resistance mutations K103N Y181C and G190A in long term Nevirapine-containing antiviral therapy

Objective

We built a cohort study of HIV patients taking long-term first-line Antiretroviral Therapy in 2003. In this assay, we focused on the development of primary drug resistance mutations against Non-Nucleoside Reverse Transcriptase Inhibitor (NNRTI), K103N, Y181C and G190A.

Method

The cohort study was built in Henan province, China. We used Single Genome Amplification (SGA) to analyze the frequency of K103N, Y181C and G190A in serial plasma samples of three individual patients. We also performed standard genotype HIV drug resistance assay in 204 patients of this cohort study to analyze the frequency of these mutations.

Result

In the SGA sequences, the K103N decreased and vanished, while the frequency of Y181C and G190A increased in individual patient receiving long-term Antiretroviral Therapy (ART). In the sequences of standard genotype HIV drug resistance assay, the frequency of K103N, Y181C and G190A had the similar pattern with that in SGA sequences. Among these patients, the viral suppression were still sufficient after receiving ART for 72 months, and 78.6% (160/204) patients could have their CD4 count over than 200cells/ul.

Conclusion

In some patients, first-line ART had the possibility to provide sufficient treatment effect for over than 72 months, but in long-term treatment, the dominant NNRTI drug resistance mutation K103N could reduced, while the proportion of variants with mutation Y181C or G190A may increased. This result was not similar with that in vitro study, which state that variant with K103N or Y181C had an equal viral fitness with wild type.

Prediction of HIV drug resistance from genotype with encoded three-dimensional protein structure

BACKGROUND

Drug resistance has become a severe challenge for treatment of HIV infections. Mutations accumulate in the HIV genome and make certain drugs ineffective. Prediction of resistance from genotype data is a valuable guide in choice of drugs for effective therapy.

RESULTS

In order to improve the computational prediction of resistance from genotype data we have developed a unified encoding of the protein sequence and three-dimensional protein structure of the drug target for classification and regression analysis. The method was tested on genotype-resistance data for mutants of HIV protease and reverse transcriptase. Our graph based sequence-structure approach gives high accuracy with a new sparse dictionary classification method, as well as support vector machine and artificial neural networks classifiers. Cross-validated regression analysis with the sparse dictionary gave excellent correlation between predicted and observed resistance.

CONCLUSION

The approach of encoding the protein structure and sequence as a 210-dimensional vector, based on Delaunay triangulation, has promise as an accurate method for predicting resistance from sequence for drugs inhibiting HIV protease and reverse transcriptase.

Impact of Y181C and/or H221Y mutation patterns of HIV-1 reverse transcriptase on phenotypic resistance to available non-nucleoside and nucleoside inhibitors in China

BACKGROUND

The aim of this study was to investigate the role of K101Q, Y181C and H221Y emerging in HIV-1 reverse transcriptase with different mutations patterns in phenotypic susceptibility to currently available NNRTIs (nevirapine NVP, efavirenz EFV) and NRTIs (zidovudine AZT, lamivudine 3TC, stavudine d4T) in China.

METHODS

Phenotype testing of currently available NNRTIs (NVP, EFV) and NRTIs (AZT, 3TC, d4T) was performed on TZM-b1 cells using recombined virus strains. P ≤ 0.05 was defined significant considering the change of 50% inhibitory drug concentration (IC50) compared with the reference, while P ≤ 0.01 was considered to be statistically significant considering multiple comparisons.

RESULTS

Triple-mutation K101Q/Y181C/H221Y and double-mutation K101Q/Y181C resulted in significant increase in NVP resistance (1253.9-fold and 986.4-fold), while only K101Q/Y181C/H221Y brought a 5.00-fold significant increase in EFV resistance. Remarkably, K101Q/H221Y was hypersusceptible to EFV (FC = 0.04), but was significantly resistant to the three NRTIs. Then, the interaction analysis suggested the interaction was not significant to NVP (F = 0.77, P = 0.4061) but significant to EFV and other three NRTIs.

CONCLUSION

Copresence of mutations reported to be associated with NNRTIs confers significant increase to NVP resistance. Interestingly, some may increase the susceptibility to EFV. Certainly, the double mutation (K101Q/H221Y) also changes the susceptibility of viruses to NRTIs. Interaction between two different sites makes resistance more complex.

Identification of the critical sites of NNRTI-resistance in reverse transcriptase of HIV-1 CRF_BC strains

Background

The polymorphisms involved in drug resistance to non-nucleoside reverse transcriptase inhibitors (NNRTIs) in HIV-1 CRF_BC, the most prevalent HIV-1 strain in China, have been poorly characterized.

Results

To reveal the drug resistance mutations, we compared the gene sequences of pol region of HIV-1 CRF_BC from 631 treatment-naïve and 363 treatment-experienced patients using the selection pressure-based method. We calculated an individual Ka/Ks value for each specific amino acid mutation. Result showed that eight polymorphic mutations (W88C, K101Q, I132L, R135L, T139K/R, H221Y and L228R) in RT for treatment-experienced patients were identified, while they, except for R135L, were completely absent in those from treatment-naïve patients. The I132L and T139K/R mutants exhibited high-level resistance to DLV and NVP and moderate resistance to TMC-125 and EFV, while the K101Q and H221Y mutants exhibited an increased resistance to all four NNRTIs tested. The W88C, R135L, and L228R may be RTI-induced adaptive mutations. Y181C+K101Q mutant showed a 2.5-, 4.4-, and 4.7-fold higher resistance to TMC-125, NVP and EFV, respectively, than Y181C alone mutant, while Y181C+H221Y or K103N+H221Y mutants had significantly higher resistance to all four NNRTIs than Y181C or K103N mutants. K103N+T139K and G190A+T139K mutant induce higher resistance (2.0∼14.2-fold and 1.5∼7.2-fold, respectively) to all four NNRTIs than K103N or G190A alone mutation.

Conclusions

I132L and T139K/R are rare but critical mutations associated with NNRTI-resistance for some NNRTIs. K101Q, H221Y and T139K can enhance K103N/Y181C/G190A-assocated NNRTI-resistance. Monitoring these mutations will provide useful information for rational design of the NNRTI-based antiretroviral regimen for HIV-1 CRF_BC-infected patients.

A polymorphism at position 400 in the connection subdomain of HIV-1 reverse transcriptase affects sensitivity to NNRTIs and RNaseH activity

Reverse transcriptase (RT) plays an essential role in HIV-1 replication, and inhibition of this enzyme is a key component of HIV-treatment. However, the use of RT inhibitors can lead to the emergence of drug-resistant variants. Until recently, most clinically relevant resistance mutations were found in the polymerase domain of RT. Lately, an increasing number of resistance mutations has been identified in the connection and RNaseH domain. To further explore the role of these domains we analyzed the complete RT sequence of HIV-1 subtype B patients failing therapy. Position A/T400 in the connection subdomain is polymorphic, but the proportion of T400 increases from 41% in naïve patients to 72% in patients failing therapy. Previous studies suggested a role for threonine in conferring resistance to nucleoside RT inhibitors. Here we report that T400 also mediates resistance to non-nucleoside RT inhibitors. The susceptibility to NVP and EFV was reduced 5-fold and 2-fold, respectively, in the wild-type subtype B NL4.3 background. We show that substitution A400T reduces the RNaseH activity. The changes in enzyme activity are remarkable given the distance to both the polymerase and RNaseH active sites. Molecular dynamics simulations were performed, which provide a novel atomistic mechanism for the reduction in RNaseH activity induced by T400. Substitution A400T was found to change the conformation of the RNaseH primer grip region. Formation of an additional hydrogen bond between residue T400 and E396 may play a role in this structural change. The slower degradation of the viral RNA genome may provide more time for dissociation of the bound NNRTI from the stalled RT-template/primer complex, after which reverse transcription can resume.

State-of-the-art and dissemination of computational tools for drug-design purposes: a survey among Italian academics and industrial institutions

During the first edition of the Computationally Driven Drug Discovery meeting, held in November 2011 at Dompé Pharma (L'Aquila, Italy), a questionnaire regarding the diffusion and the use of computational tools for drug-design purposes in both academia and industry was distributed among all participants. This is a follow-up of a previously reported investigation carried out among a few companies in 2007. The new questionnaire implemented five sections dedicated to: research group identification and classification; 18 different computational techniques; software information; hardware data; and economical business considerations. In this article, together with a detailed history of the different computational methods, a statistical analysis of the survey results that enabled the identification of the prevalent computational techniques adopted in drug-design projects is reported and a profile of the computational medicinal chemist currently working in academia and pharmaceutical companies in Italy is highlighted.

What do molecules do when we are not looking? State sequence analysis for stochastic chemical systems

Many biomolecular systems depend on orderly sequences of chemical transformations or reactions. Yet, the dynamics of single molecules or small-copy-number molecular systems are significantly stochastic. Here, we propose state sequence analysis--a new approach for predicting or visualizing the behaviour of stochastic molecular systems by computing maximum probability state sequences, based on initial conditions or boundary conditions. We demonstrate this approach by analysing the acquisition of drug-resistance mutations in the human immunodeficiency virus genome, which depends on rare events occurring on the time scale of years, and the stochastic opening and closing behaviour of a single sodium ion channel, which occurs on the time scale of milliseconds. In both cases, we find that our approach yields novel insights into the stochastic dynamical behaviour of these systems, including insights that are not correctly reproduced in standard time-discretization approaches to trajectory analysis.

Identification of drug resistant mutations in HIV-1 CRF07_BC variants selected by nevirapine in vitro

Since the antiretroviral therapy (ART) was introduced to patients infected by human immunodeficiency virus (HIV), the HIV related mortality and morbidity have been significantly reduced. The major obstacle for long-term successful anti-HIV treatment is the emergence of drug resistant mutants. Current data of drug resistance was mainly obtained on HIV-1 subtype B but rarely on non-B virus, even more rare with newly emerged circulating recombinant forms (CRFs). The lack of such data limits the rational management of ART for the increasing number of patients infected by non-subtype B virus. In this study, a HIV-1 CRF07_BC strain CNGZD was isolated from a HIV patient and its genome was sequenced and deposited in GenBank (JQ423923). Potential drug resistant mutants of this CRF07_BC virus strain were selected in PBMCs cultures in the presence of Nevirapine (NVP), which is the most frequently used antiretroviral drug in China. Four combination profiles of mutations were identified in the NVP-selected mutants, which were initiated with A98G, V108I, Y181C and I135T/I382L and followed by more than two other mutations at the end of the selections, respectively. A total of seven previously reported mutations (A98G, V106M, V108I, I135T, Y181C, V189I, K238N) and seven novel mutations (P4H, T48I, I178M, V314A, I382L/V, T386A) in the reverse transcriptase gene were found in these NVP-selected mutants. Phenotypic analysis in the NVP-selected mutants showed that all the mutations, except P4H, contribute to NVP resistance. Among them, V106M and Y181C reduce NVP susceptibility for more than 20-fold, while the other mutations cause less than 20 folds drug resistance. Although the information obtained in this in vitro selection study may not fully cover resistant mutations which will actually occur in patients, it has still provided useful information for rational management of ART in patients infected with HIV CRF_BC subtype.

Prediction of mutational tolerance in HIV-1 protease and reverse transcriptase using flexible backbone protein design

Predicting which mutations proteins tolerate while maintaining their structure and function has important applications for modeling fundamental properties of proteins and their evolution; it also drives progress in protein design. Here we develop a computational model to predict the tolerated sequence space of HIV-1 protease reachable by single mutations. We assess the model by comparison to the observed variability in more than 50,000 HIV-1 protease sequences, one of the most comprehensive datasets on tolerated sequence space. We then extend the model to a second protein, reverse transcriptase. The model integrates multiple structural and functional constraints acting on a protein and uses ensembles of protein conformations. We find the model correctly captures a considerable fraction of protease and reverse-transcriptase mutational tolerance and shows comparable accuracy using either experimentally determined or computationally generated structural ensembles. Predictions of tolerated sequence space afforded by the model provide insights into stability-function tradeoffs in the emergence of resistance mutations and into strengths and limitations of the computational model.

Many related protein sequences can be consistent with the structure and function of a given protein, suggesting that proteins may be quite robust to mutations. This tolerance to mutations is frequently exploited by pathogens. In particular, pathogens can rapidly evolve mutated proteins that have a new function - resistance against a therapeutic inhibitor - without abandoning other functions essential for the pathogen. This principle may also hold more generally: Proteins tolerant to mutational changes can more easily acquire new functions while maintaining their existing properties. The ability to predict the tolerance of proteins to mutation could thus help both to analyze the emergence of resistance mutations in pathogens and to engineer proteins with new functions. Here we develop a computational model to predict protein mutational tolerance towards point mutations accessible by single nucleotide changes, and validate it using two important pathogenic proteins and therapeutic targets: the protease and reverse transcriptase from HIV-1. The model provides insights into how resistance emerges and makes testable predictions on mutations that have not been seen yet. Similar models of mutational tolerance should be useful for characterizing and reengineering the functions of other proteins for which a three-dimensional structure is available.

Clinical, virological and biochemical evidence supporting the association of HIV-1 reverse transcriptase polymorphism R284K and thymidine analogue resistance mutations M41L, L210W and T215Y in patients failing tenofovir/emtricitabine therapy

Background

Thymidine analogue resistance mutations (TAMs) selected under treatment with nucleoside analogues generate two distinct genotypic profiles in the HIV-1 reverse transcriptase (RT): (i) TAM1: M41L, L210W and T215Y, and (ii) TAM2: D67N, K70R and K219E/Q, and sometimes T215F. Secondary mutations, including thumb subdomain polymorphisms (e.g. R284K) have been identified in association with TAMs. We have identified mutational clusters associated with virological failure during salvage therapy with tenofovir/emtricitabine-based regimens. In this context, we have studied the role of R284K as a secondary mutation associated with mutations of the TAM1 complex.

Results

The cross-sectional study carried out with >200 HIV-1 genotypes showed that virological failure to tenofovir/emtricitabine was strongly associated with the presence of M184V (P < 10^-10) and TAMs (P < 10^-3), while K65R was relatively uncommon in previously-treated patients failing antiretroviral therapy. Clusters of mutations were identified, and among them, the TAM1 complex showed the highest correlation coefficients. Covariation of TAM1 mutations and V118I, V179I, M184V and R284K was observed. Virological studies showed that the combination of R284K with TAM1 mutations confers a fitness advantage in the presence of zidovudine or tenofovir. Studies with recombinant HIV-1 RTs showed that when associated with TAM1 mutations, R284K had a minimal impact on zidovudine or tenofovir inhibition, and in their ability to excise the inhibitors from blocked DNA primers. However, the mutant RT M41L/L210W/T215Y/R284K showed an increased catalytic rate for nucleotide incorporation and a higher RNase H activity in comparison with WT and mutant M41L/L210W/T215Y RTs. These effects were consistent with its enhanced chain-terminated primer rescue on DNA/DNA template-primers, but not on RNA/DNA complexes, and can explain the higher fitness of HIV-1 having TAM1/R284K mutations.

Conclusions

Our study shows the association of R284K and TAM1 mutations in individuals failing therapy with tenofovir/emtricitabine, and unveils a novel mechanism by which secondary mutations are selected in the context of drug-resistance mutations.

Characterization of drug resistance mutations in naïve and ART-treated patients infected with HIV-1 in Yaounde, Cameroon

Currently the prevalence of HIV-1 infection in Cameroon is 5.1%, CRF02_AG subtype is responsible for about 50% of infections. Since an HIV-1 drug resistance test is not yet available widely, accurate data on the prevalence of resistant viral strains are missing. The objective of this study was to determine HIV-1 genetic diversity and to characterize HIV-1 mutations conferring drug resistance among antiretroviral therapy (ART)-naïve and ART-treated patients. A cohort of 239 patients infected with HIV were followed-up between January 2007 and July 2010 in Cameroon. Two hundred and sixteen plasma samples were sequenced for phylogenetic analysis and identification of drug resistance mutations in the HIV-1 pol region. A significant genetic diversity was found: Seven pure subtypes (A1, A3, D, F1, F2, G, H), nine circulating recombinant forms (CRFs: 01_AE, 02_AG, 06cpx, 09cpx, 11cpx, 13cpx, 16cpx, 18cpx, 37cpx) and one new unique recombinant form (URF) (G/F2). The rate of transmitted drug resistance (TDR) in naïve patients was 8.2% (4/49). Around 80% of patients failing a first-line ART harbored a virus with at least one resistance mutation to two antiretroviral (ARV) classes, and 36% of those failing a second-line regimen carried a virus with at least one resistant mutation to three ARV classes. The high level of drug resistance observed in the cohort is alarming because this occurred as a result of only few years of treatment. Adherence to therapy, adequate education of physicians, and the appropriate use of genotypic resistance assay are critical points of intervention for the improvement of patient care.

Phenotypic analysis of HIV-1 genotypic drug-resistant isolates from China, using a single-cycle system

BACKGROUND AND OBJECTIVES

Drug resistance in HIV-1 is one of the main causes of failure of antiretroviral therapy. Phenotypic detection of drug-resistant HIV-1 can provide guidance in selecting the optimal treatment regimen. Traditional phenotype assays are labor intensive and time consuming. Thus, a rapid and convenient phenotype assay with a single cycle of replication was developed and used in this study.

METHODS

Two restriction endonuclease sites, ApaI and AgeI, were inserted into the plasmid pSG3Δenv(,) using site-directed mutagenesis. The reverse transcriptase and protease genes of HIV-1 were amplified from patients and cloned into the modified pSG3Δenv. Sixteen original recombinant pseudoviruses were generated. The phenotypic susceptibility of these 16 recombinant pseudoviruses to 12 antiretroviral drugs was determined using a luciferase reporter system, and the phenotype and genotype results were compared.

RESULTS

A modified phenotype assay with a single-cycle system was established, and its reproducibility and feasibility were validated. Approximately 89% of the phenotype results were in agreement with the genotype results; this slight disagreement may have been due to complex and multiple resistance mutations. The phenotype results showed that individual pseudoviruses with four thymidine analog mutations (TAMs).[M41L, T67N, L210W, and T215Y] in combination with various other mutations had different levels of resistance to nucleoside reverse transcriptase inhibitors (NRTIs). Mutations E44A, T69D, and V118I influenced the pattern of resistance of TAMs. The level of resistance to non-NRTIs (NNRTIs) was also variable when different NNRTI-resistance mutations were combined.

CONCLUSION

The single-cycle pseudovirus phenotypic susceptibility detection system reflects HIV-1 drug resistance, especially for complex resistance mutants, and could be used to screen new antiretroviral candidates.

Impact of Novel Resistance Profiles in HIV-1 Reverse Transcriptase on Phenotypic Resistance to NVP

Objective. To clarify the impact of H221Y mutation on drug resistance to NVP. Methods. 646 bp HIV-1 pol gene fragments (from 592 to 1237 nucleotide) with different NNRTIs mutation profiles from AIDS patients receiving antiretroviral therapy containing NVP regimens were introduced into pNL4-3 backbone plasmid. H221Y and (or) Y181C mutations were reverted to wild type amino acids by site-directed mutagenesis, then strains containing various mutation patterns were packaged. Phenotypic drug resistance was analyzed on TZM-bl cells. Results. 12 strains containing different drug-resistant mutation profiles were constructed, including the K101Q series (K101Q/Y181C/H221Y, K101Q/Y181C, K101Q/H221Y, and K101Q), the V179D series (V179D/Y181C/H221Y, V179D/Y181C, V179D/H221Y, and V179D), and the K103N series (K103N/Y181C/H221Y, K103N/Y181C, K103N/H221Y, K103N). For strains containing the mutation profiles (K101Q/Y181C, K101Q, V179D/Y181C, V179D, K103N/Y181C, and K103N), the presence of H221Y reduced NVP susceptibility by 2.1 ± 0.5 to 3.6 ± 0.5 fold. To the mutation profiles K101Q/H221Y, K101Q, V179D/H221Y, V179D, K103N/H221Y, and K103N, the presence of Y181C reduced NVP susceptibility by 41.9 ± 8.4 to 1297.0 ± 289.1 fold. For the strains containing K101Q, V179D, and K103N, the presence of Y181C/H221Y combination decreased NVP susceptibility by 100.6 ± 32.5 to 3444.6 ± 834.5 fold. Conclusion. On the bases of various NNRTIs mutation profiles, Y181C remarkably improved the IC₅₀ to NVP, although H221Ymutation alone just increases 2.1 ∼ 3.6-fold resistance to NVP, the mutation could improve 100.6 ∼ 3444.6-fold resistance to NVP when it copresent with Y181C, the phenotypic drug resistance fold was improved extremely. For strains containing the mutation profiles (K101Q/Y181C, K101Q, V179D/Y181C, V179D, K103N/Y181C, and K103N), the presence of H221Y reduced NVP susceptibility by 2.1 ± 0.5 to 3.6 ± 0.5 fold.

Polymorphic mutations associated with the emergence of the multinucleoside/tide resistance mutations 69 insertion and Q151M

BACKGROUND

We hypothesized that polymorphic mutations exist that are associated with the emergence of the multinucleoside resistance mutations (MNR), 69 insertion and Q151M.

METHODS

The Swiss HIV Cohort Study was screened, and the frequencies of polymorphic mutations in HIV-1 (subtype B) were compared between patients detected with the 69 insertion (n = 17), Q151M (n = 29), ≥2 thymidine analogue mutations (TAM) 1 (n = 400) or ≥2 TAM 2 (n = 249). Logistic regressions adjusted for the antiretroviral treatment history were performed to analyze the association of the polymorphic mutations with MNR.

RESULTS

The 69 insertion and TAM 1 were strongly associated and occurred in 94.1% (16 of 17) together. The 69 insertion seemed to emerge as a consequence of the TAM 1 pathway (median years until detection: 6.8 compared with 4.4 for ≥2 TAM 1, P Wilcoxon = 0.009). Frequencies of 8 polymorphic mutations (K43E, V60I, S68G, S162C, T165I, I202V, R211K, F214L) were significantly different between groups. Logistic regression showed that F214L and V60I were associated with the emergence of Q151M/TAM 2 opposed to 69 insertion/TAM 1. S68G, T165I, and I202V were associated with Q151M instead of TAM 2.

CONCLUSIONS

Besides antiretroviral therapy, polymorphic mutations may contribute to the emergence of specific MNR mutations.

Docking analysis and resistance evaluation of clinically relevant mutations associated with the HIV-1 non-nucleoside reverse transcriptase inhibitors nevirapine, efavirenz and etravirine

An integrated computational and statistical approach was used to determine the association of non-nucleoside reverse transcriptase inhibitors (NNRTIs) nevirapine, efavirenz and etravirine with resistance mutations that cause therapeutic failure and their impact on NNRTI resistance. Mutations detected for nevirapine virological failure with a prevalence greater than 10% in the used patient set were: K103N, Y181C, G190A, and K101E. A support vector regression model, based on matched genotypic/phenotypic data (n=850), showed that among 6365 analyzed mutations, K103N, Y181C and G190A have the first, third, and sixth greatest significance for nevirapine resistance, respectively. The most common indicator of treatment failure for efavirenz was K103N mutation present in 56.7% of the patients where the drug failed, followed by V108I, L100I, and G190A. For efavirenz resistance, K103N, G190, and L100I have the first, fourth, and eighth greatest significance, respectively, as determined in support vector regression model. No positive interactions were observed among nevirapine resistance mutations, while a more complex situation was observed with treatment failure of efavirenz and etravirine, characterized by the accumulation of multiple mutations. Docking simulations and free energy analysis based on docking scores of mutated human immunodeficiency virus (HIV) RT complexes were used to evaluate the influence of selected mutations on drug recognition. Results from support vector regression were confirmed by docking analysis. In particular, for nevirapine and efavirenz, a single mutation K103N was associated with the most unfavorable energetic profile compared to the wild-type sequence. This is in line with recent clinical data reporting that diarylpyrimidine etravirine, a very potent third generation drug effective against a wide range of drug-resistant HIV-1 variants, shows increased affinity towards K103N/S mutants due to its high conformational flexibility.

HIV-1 mutational pathways under multidrug therapy

BACKGROUND

Genotype-derived drug resistance profiles are a valuable asset in HIV-1 therapy decisions. Therapy decisions could be further improved, both in terms of predicting length of current therapy success and in preserving followup therapy options, through better knowledge of mutational pathways- here defined as specific locations on the viral genome which, when mutant, alter the risk that additional specific mutations arise. We limit the search to locations in the reverse transcriptase region of the HIV-1 genome which host resistance mutations to nucleoside (NRTI) and non-nucleoside (NNRTI) reverse transcriptase inhibitors (as listed in the 2008 International AIDS Society report), or which were mutant at therapy start in 5% or more of the therapies studied.

METHODS

A Cox proportional hazards model was fit to each location with the hazard of a mutation at that location during therapy proportional to the presence/absence of mutations at the remaining locations at therapy start. A pathway from preexisting to occurring mutation was indicated if the covariate was both selected as important via smoothly clipped absolute deviation (a form of regularized regression) and had a small p-value. The Cox model also allowed controlling for non-genetic parameters and potential nuisance factors such as viral resistance and number of previous therapies. Results were based on 1981 therapies given to 1495 distinct patients drawn from the EuResist database.

RESULTS

The strongest influence on the hazard of developing NRTI resistance was having more than four previous therapies, not any one existing resistance mutation. Known NRTI resistance pathways were shown, and previously speculated inhibition between the thymidine analog pathways was evidenced. Evidence was found for a number of specific pathways between NRTI and NNRTI resistance sites. A number of common mutations were shown to increase the hazard of developing both NRTI and NNRTI resistance. Viral resistance to the therapy compounds did not materially effect the hazard of mutation in our model.

CONCLUSIONS

The accuracy of therapy outcome prediction tools may be increased by including the number of previous treatments, and by considering locations in the HIV genome which increase the hazard of developing resistance mutations.

Selective inhibition of Human Immunodeficiency Virus type 1 (HIV-1) by a novel family of tricyclic nucleosides

Nucleoside 1, with an unusual tricyclic carbohydrate moiety, specifically inhibits HIV-1 replication while being inactive against HIV-2 or other (retro) viruses. In an attempt to increase the inhibitory efficacy against HIV-1, and to further explore the structural features required for anti-HIV-1 activity, different types of modifications have been carried out on this prototype compound. These include substitution of the ethoxy group at the C-4″ position by alkoxy groups of different length, branching, conformational freedom or functionalization. In addition, the 4″-ethoxy group has been removed or substituted by other functional groups. The role of the tert-butyldimethylsilyl (TBDMS) group at the 2' position has also been studied by preparing the corresponding 2'-deprotected derivative or by replacing it by other silyl (tert-hexyldimethylsilyl) or acyl (acetyl) moieties. Finally, the thymine of the prototype compound has been replaced by N-3-methylthymine, uracil or thiophenyl. Some of these compounds were endowed with a 6- to 7-fold higher selectivity than the prototype 1. The tricyclic nucleosides here described represent a novel type of selective anti HIV-1 inhibitors, targeted at the HIV-1-encoded reverse transcriptase.

Indolylarylsulfones as HIV-1 non-nucleoside reverse transcriptase inhibitors: new cyclic substituents at indole-2-carboxamide

New indolylarylsulfone derivatives bearing cyclic substituents at indole-2-carboxamide linked through a methylene/ethylene spacer were potent inhibitors of the WT HIV-1 replication in CEM and PBMC cells with inhibitory concentrations in the low nanomolar range. Against the mutant L100I and K103N RT HIV-1 strains in MT-4 cells, compounds 20, 24-26, 36, and 40 showed antiviral potency superior to that of NVP and EFV. Against these mutant strains, derivatives 20, 24-26, and 40 were equipotent to ETV. Molecular docking experiments on this novel series of IAS analogues have also suggested that the H-bond interaction between the nitrogen atom in the carboxamide chain of IAS and Glu138:B is important in the binding of these compounds. These results are in accordance with the experimental data obtained on the WT and on the mutant HIV-1 strains tested.

Molecular and structural aspects of clinically relevant mutations related to the approved non-nucleoside inhibitors of HIV-1 reverse transcriptase

In recent years relevant progress has been made in the treatment of HIV-1 with a consequent decrease in mortality. The availability of potent antiretroviral drugs and the ability of viral load assays that accurately evaluate the true level of viral replication, have led to a better understanding of pathogenesis of the disease and how to obtain improved therapeutic profiles. The highly active antiretroviral therapy (HAART), based on a combination of three or more antiretroviral drugs, has radically changed the clinical outcome of HIV. In particular, reverse transcriptase non-nucleoside inhibitors (NNRTIs) play an essential role in most protocols and are often used in first line treatment. The high specificity of these inhibitors towards HIV-1 has increased the number of structural and molecular modeling studies of enzyme complexes and that have led to chemical syntheses of more selective second and third-generation NNRTIs. However, a considerable percentage of new HIV-1 infections are caused by the emergence of drug-resistant mutant viruses that complicate treatment strategies. In this review we discuss relevant clinical and structural aspects for the management of antiretroviral drug resistance, with detailed explanations of mechanisms and mutation patterns useful to better understand the relation between drug resistance and therapy failure.

Mechanisms involved in the selection of HIV-1 reverse transcriptase thumb subdomain polymorphisms associated with nucleoside analogue therapy failure

Previous studies showed an increased prevalence of human immunodeficiency virus type 1 (HIV-1) reverse transcriptase (RT) thumb subdomain polymorphisms Pro272, Arg277, and Thr286 in patients failing therapy with nucleoside analogue combinations. Interestingly, wild-type HIV-1(BH10) RT contains Pro272, Arg277, and Thr286. Here, we demonstrate that in the presence of zidovudine, HIV-1(BH10) RT mutations P272A/R277K/T286A produce a significant reduction of the viral replication capacity in peripheral blood mononuclear cells in both the absence and presence of M41L/T215Y. In studies carried out with recombinant enzymes, we show that RT thumb subdomain mutations decrease primer-unblocking activity on RNA/DNA complexes, but not on DNA/DNA template-primers. These effects were observed with primers terminated with thymidine analogues (i.e., zidovudine and stavudine) and carbovir (the relevant derivative of abacavir) and were more pronounced when mutations were introduced in the wild-type HIV-1(BH10) RT sequence context. RT thumb subdomain mutations increased by 2-fold the apparent dissociation equilibrium constant (K(d)) for RNA/DNA without affecting the K(d) for DNA/DNA substrates. RNase H assays carried out with RNA/DNA complexes did not reveal an increase in the reaction rate or in secondary cleavage events that could account for the decreased excision activity. The interaction of Arg277 with the phosphate backbone of the RNA template in HIV-1 RT bound to RNA/DNA and the location of Thr286 close to the RNA strand are consistent with thumb polymorphisms playing a role in decreasing nucleoside RT inhibitor excision activity on RNA/DNA template-primers by affecting interactions with the template-primer duplex without involvement of the RNase H activity of the enzyme.

Analysis of selection pressure and mutational pattern of HIV type 1 reverse transcriptase region among treated and nontreated patients

Variation of the HIV-1 subtype C reverse transcriptase region (RT) resulting in response to the selective pressures of drug therapy remains poorly characterized. Here, we compared the genetic variation resulting in the presence and absence of antiretroviral drug selective pressures on HIV-1 subtype C RT among nontreated and treated patients. The nucleotide variability, nonsynonymous and synonymous ratio, and the positively selected mutations were determined by comparing the RT sequences isolated at two time points among nontreated (baseline and follow-up) and treated patients (baseline and treatment failure). Compared to the nontreated patients, the intrapatient nucleotide variability, the number of nonsynonymous and synonymous substitutions was significantly higher among the treated patients. Among the mutations positively selected, the frequency of D121Y, I135R, and Q207E increased and the frequency of mutation S48T decreased significantly during treatment failure. Further studies are essential to discover the role of these mutations during treatment in HIV-1 subtype C.

Constrained patterns of covariation and clustering of HIV-1 non-nucleoside reverse transcriptase inhibitor resistance mutations

Objectives

We characterized pairwise and higher order patterns of non-nucleoside reverse transcriptase inhibitor (NNRTI)-selected mutations because multiple mutations are usually required for clinically significant resistance to second-generation NNRTIs.

Patients and methods

We analysed viruses from 13 039 individuals with sequences containing at least one of 52 published NNRTI-selected mutations, including 1133 viruses from individuals who received efavirenz but no other NNRTI and 1510 viruses from individuals who received nevirapine but no other NNRTI. Of the 17 reported etravirine resistance-associated mutations (RAMs), Y181C/I/V, L100I, K101P and M230L were considered major based on published in vitro susceptibility data.

Results

Efavirenz preferentially selected for 16 mutations, including L100I (14% versus 0.1%, P < 0.001), K101P (3.3% versus 0.4%, P < 0.001) and M230L (2.8% versus 1.3%, P = 0.004), whereas nevirapine preferentially selected for 12 mutations, including Y181C/I/V (48% versus 6.9%, P < 0.001). Twenty-nine pairs of NNRTI-selected mutations covaried significantly, including Y181C with seven other mutations (A98G, K101E/H, V108I, G190A/S and H221Y), L100I with K103N, and K101P with K103S. Two pairs (Y181C + V179F and Y181C + G190S) were predicted to confer >10-fold decreased etravirine susceptibility. Seventeen percent of sequences had three or more NNRTI-selected mutations, mostly in clusters of covarying mutations. Many clusters had Y181C plus a non-major etravirine RAM; few had more than one major etravirine RAM.

Conclusions

Although major etravirine RAMs rarely occur in combination, 2 of 29 pairs of covarying mutations were associated with >10-fold decreased etravirine susceptibility. Viruses with three or more NNRTI-selected mutations often contained Y181C in combination with one or more minor etravirine RAMs; however, phenotypic and clinical correlates for most of these higher order combinations have not been published.

Resistance pathways of human immunodeficiency virus type 1 against the combination of zidovudine and lamivudine

A better understanding of human immunodeficiency virus type 1 drug-resistance evolution under the selective pressure of combination treatment is important for the design of long-term effective treatment strategies. We applied Bayesian network learning to sequences from patients treated with the reverse transcriptase inhibitor combination of zidovudine (AZT) and lamivudine (3TC) to identify the role of many treatment-selected mutations in the development of resistance. Based on the Bayesian network structure, an in vivo fitness landscape was built, reflecting the necessary selective pressure under treatment, to evolve naive sequences to sequences obtained from patients treated with the combination. This landscape, combined with an evolutionary model, was used to predict resistance evolution in longitudinal sequence pairs. In our analysis, mutations 41L, 70R, 184V and 215F/Y were identified as major resistance mutations to the combination of AZT and 3TC, as they were associated directly with treatment experience. The network also suggested a possible role in resistance development for a number of novel mutations. Estimated fitness, using the landscape, correlated significantly with in vitro resistance phenotype in genotype-phenotype pairs (R(2)=0.70). Variation in predicted evolution under selective pressure correlated significantly with observed in vivo evolution during AZT plus 3CT treatment. In conclusion, we confirmed current knowledge on resistance development to the combination of AZT and 3CT, but additional novel mutations were identified. Moreover, a model to predict resistance evolution during AZT and 3CT treatment has been built and validated.

Impact of human leukocyte antigen-B*51-restricted cytotoxic T-lymphocyte pressure on mutation patterns of nonnucleoside reverse transcriptase inhibitor resistance

OBJECTIVE

The objective of this study is to determine the impact of human leukocyte antigen (HLA)-B*51-restricted cytotoxic T-lymphocyte (CTL) pressure on the development of nonnucleoside reverse transcriptase inhibitor (NNRTI) resistance.

DESIGN

The prevalence of HIV-1 harboring an escape mutation, I135X, in a major epitope of HLA-B*51-restricted CTL located in reverse transcriptase is increasing worldwide. We analyzed the effects of escape mutations on the emerging mutation patterns of NNRTI resistance.

METHODS

Monoclonal HIV-1 sequences harboring each of the escape mutations, including I135L (HIV-1I135L), I135V (HIV-1I135V), I135T (HIV-1I135T), and I135R (HIV-1I135R) in reverse transcriptase, and a wild-type monoclonal HIV-1 (HIV-1WT) were cultured in the presence of increasing concentrations of efavirenz. Induced mutations during culture passages of the culture were analyzed.

RESULTS

E138K emerged during the cultural passages of HIV-1I135V, HIV-1I135T, and HIV-1I135R, but not during the passages of HIV-1WT. The combination of I135T, the most frequent escape mutation, and E138K (HIV-1I135T/E138K) conferred significant resistance to efavirenz, nevirapine, and etravirine. The HIV-1I135L/E138K and HIV-1I135R/E138K were significantly resistant to nevirapine and etravirine, respectively, though each solo of escape mutations and E138K did not confer significant resistance to NNRTI. Computational analysis indicated that I135T and E138K cooperatively extend the gap between the binding site of reverse transcriptase and NNRTI.

CONCLUSION

HLA-B*51-restricted CTL can induce novel mutation patterns of NNRTI resistance by selecting escape mutations. The spread of CTL escape variants may alter the mutation patterns of drug resistance.

Five-year follow up of genotypic resistance patterns in HIV-1 subtype C infected patients in Botswana after failure of thymidine analogue-based regimens

Objective

Our objective was to establish genotypic resistance profiles among the 4% of Batswana patients who experienced virologic failure while being followed within Botswana's National Antiretroviral Treatment Program between 2002 and 2007.

Methods

At the beginning of the national program in 2002, almost all patients received stavudine (d4T), together with didanosine (ddI), as part of their first nucleoside reverse transcriptase inhibitor (NRTI)-based regimen (Group 1). In contrast, the standard of care for all patients subsequently enrolled (2002-2007) included zidovudine/lamivudine (ZDV/3TC) (Group 2). Genotypes were analyzed in 26 patients from Group 1 and 37 patients from Group 2. Associations between mutations were determined using Pearson's correlation coefficient and Jaccard's coefficient of similarity.

Results

Seventy-eight percent of genotyped patients possessed mutations associated with protease inhibitor (PI) resistance while 87% and 90%, respectively, exhibited mutations associated with NRTIs and non-nucleoside reverse transcriptase inhibitors (NNRTIs). The most frequent PI mutations involving resistance to NFV were L90M (25.2%) and D30N (16.2%), but mutations at positions K45Q and D30N were often observed in tandem (P = 60.5, J = 50; p = 0.002; Group 2) alongside Q61E in 42.8% of patients who received ZDV/3TC. Both major patterns of thymidine analogue mutations, TAM 1 (48%) and TAM 2 (59%), were represented in patients from Group 1 and 2, although M184V was higher among individuals who had initially received ddI (61% versus 40.5%). In contrast, L74V was more frequent among individuals from Group 2 (16.2% versus 7.7%). Differences in regard to NNRTI mutations were also observed between Group 1 and Group 2 patients.

Conclusion

Despite a low rate of therapeutic failure (4%) among these patients, those who failed possessed high numbers of resistance mutations as well as novel resistance mutations and/or polymorphisms at sites within reverse transcriptase and protease.

HIV-I reverse transcriptase variation in plasma and genital secretion of antiretroviral-naive females

The reverse transcriptase (RT) enzyme of HIV type 1 (HIV-1) is largely targeted by the host immune selection pressure and would differ in the anatomical compartments, thereby having a drastic impact on viral quasi-species evolution. The HIV-1 RT region sequenced from plasma and genital secretions of 8 antiretroviral treatment (ART)-naive females was analyzed for the pattern of amino acid mutations and the ratio of synonymous and nonsynonymous substitutions to determine whether it is under different selection pressure in both the compartments. Phylogenetic and mutational analysis of the HIV-1 RT in plasma and genital secretions of HIV-1-infected ART-naive females showed limited variation likely reflecting the absence of differential selection pressure and therefore genetic variation in these compartments.

The profile of mutational clusters associated with lamivudine resistance can be constrained by HBV genotypes

BACKGROUND/AIMS

To investigate the different clusters of mutations associated with lamivudine resistance in HBV genotypes D and A.

METHODS

HBV reverse transcriptase sequences of 89 HBV-infected patients failing lamivudine treatment were analyzed. The association of mutations with HBV genotypes was assessed by Chi-Squared test and multivariate logistic regression analysis. Covariate analysis was based on hierarchical clustering.

RESULTS

In genotype A, the rtM204V (prevalence: 68.2%) was the main sign of lamivudine failure. Multivariate analysis confirmed that genotype A is the only predictor for rtM204V emergence (OR: 14.5 [95% CI: 1.3-158], P=0.02). Covariate analysis showed that rtM204V clusters with rtL180M, rtL229V (corresponding to sF220L in the HBsAg), and, interestingly, with HBsAg mutation sS207N (bootstrap=0.95). Both sF220L and sS207N co-localized in the fourth transmembrane HBsAg domain. In contrast, in genotype D the primary mutations rtM204V and rtM204I occurred with similar prevalence (39.1% versus 45.3%, P=0.47), and showed a distinct pattern of compensatory mutations. rtM204V clusters with mutations localized in the RT-B domain (rtV173L, rtL180M, and rtT184A/S) (bootstrap=0.94), while rtM204I clusters with mutations localized in the RT-A domain (rtS53N, rtT54Y, and rtL80I/V) (bootstrap=0.96) (without associations with HBsAg specific mutations).

CONCLUSIONS

HBV genotype plays an important role in driving RT evolution under lamivudine treatment, and thus can be relevant for therapeutic sequencing, immunological response and disease progression.