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

Identification of a novel small-molecule inhibitor of the HIV-1 reverse transcriptase activity with a non-nucleoside mode of action

BACKGROUND

Human immunodeficiency virus-1 (HIV-1) is the causative agent of acquired immunodeficiency syndrome, which is a major global health problem. Although combination antiretroviral therapy (cART) successfully expands the lifespan of HIV-1-infected patients, long-term cART often increases drug resistance and adverse effects. Therefore, efforts are ongoing to develop novel anti-HIV-1 drugs.

METHODS

The anti-HIV-1 activities of compounds were investigated using TZM-bl reporter cell line, A3.01 T cell line, and peripheral blood mononuclear cells infected with several HIV-1 strains, including wild type and drug-resistance associated mutants. Next-generation sequencing analysis and in silico molecular docking studies were employed to determine the mode of action of the compound.

RESULTS

We identified a small-molecule inhibitor consisting of a thiadiazole core appended to two pyrazoles (BPPT), which exerted a highly potent inhibitory effect on HIV-1 infectivity, with a half-maximal effective concentration (EC50) of 60 nM, without causing cytotoxicity. In experiments with various HIV-1 strains and cell types, the potency of BPPT was found to be comparable to that of commercial antiretroviral agents (azidothymidine, nevirapine, and others). Further analysis of the mode of action demonstrated that BPPT is a novel type of HIV-1 non-nucleoside reverse transcriptase inhibitor (NNRTI). Analysis of viruses harboring drug-resistance-associated mutations showed that BPPT was potent against G190A (C or S) mutations in reverse transcriptase (RTase), exhibiting high-level resistance to other NNRTIs. Next-generation sequencing analysis of long-term treatment with BPPT displayed an RTase mutation profile different from that in the case of established NNRTIs. Given these data, in silico molecular docking studies demonstrated the molecular mechanism underlying the BPPT-mediated inhibition of RTase.

CONCLUSION

Our data suggest that BPPT is a novel small-molecule inhibitor of HIV-1 RTase and could serve as a promising chemical scaffold to complement or replace conventional treatments, particularly for overcoming resistance associated with the G190 mutation.

Exploring novel HIV-1 reverse transcriptase inhibitors with drug-resistant mutants: A double mutant surprise

HIV-1 reverse transcriptase (RT) remains a key target for HIV drug development. As successful management of the disease requires lifelong treatment, the emergence of resistance mutations is inevitable, making development of new RT inhibitors, which remain effective against resistant variants crucial. To this end, previous computationally guided drug design efforts have resulted in catechol diether compounds, which inhibit wildtype RT with picomolar affinities and appear to be promising preclinical candidates. To confirm that these compounds remain potent against Y181C, a widespread mutation conferring resistance to first generation inhibitors, they were screened against the HIV-1 N119 clinical isolate, reported as a Y181C single mutant. In comparison to a molecular clone with the same mutation, N119 appears less susceptible to inhibition by our preclinical candidate compounds. A more detailed sequencing effort determined that N119 was misidentified and carries V106A in combination with Y181C. While both indolizine and naphthalene substituted catechol diethers are potent against the classical Y181C single mutant, the addition of V106A confers more resistance against the indolizine derivatives than the naphthalene derivatives. Crystal structures presented in this study highlight key features of the naphthyl group, which allow these compounds to remain potent in the double mutant, including stronger interactions with F227 and less reliance on V106 for stabilization of the ethoxy-uracil ring, which makes critical hydrogen bonds with other residues in the binding pocket.

Low prevalence of doravirine-associated resistance mutations among polish human immunodeficiency-1 (HIV-1)–infected patients

Introduction

Doravirine (DOR) is a novel non-nucleoside reverse transcriptase inhibitor (NNRTI) that retains activity against common NNRTI resistance mutations. In this study, we aimed to investigate the prevalence of DOR resistance mutations compared with that of resistance mutations for other NNRTIs among HIV-1-infected treatment‐experienced and -naïve patients from Poland.

Methods

Resistance to DOR and other NNRTIs was assessed in two datasets: 1760 antiretroviral treatment-naïve HIV-1 patients and 200 treatment‐experienced patients. All 1960 sequences were derived from the patients using bulk sequencing. For resistance analyses, Stanford HIV drug resistance database scores were used.

Results

Overall, DOR resistance was present in 32 patients (1.62%), of whom 13 (0.74%) were naïve and 19 (9.50%) were treatment-experienced. The most common DOR resistance mutations observed among the naïve patients were A98G and K101E (0.2% each), and those among cART-experienced patients were L100I (2.0%), K101E, V108I, H221Y, and P225H (1.5% each). Furthermore, among the naïve patients, less common resistance to DOR (0.7%) compared with that to nevirapine (NVP) (2.1%; p = 0.0013) and rilpivirine (5.40%; p < 0.0001) was observed. For sequences obtained from treatment-experienced patients, the frequency of resistance to DOR (9.5%) was lower than that for efavirenz (25.5%; p < 0.0001) and NVP (26.0%; p < 0.0001).

Conclusions

The frequency of transmitted drug resistance to DOR is low, allowing for effective treatment of antiretroviral treatment-naïve patients and rapid treatment initiation. In cART-experienced patients, this agent remains an attractive NNRTI option with a higher genetic barrier to resistance.

Nationwide Study of Drug Resistance Mutations in HIV-1 Infected Individuals under Antiretroviral Therapy in Brazil

The success of antiretroviral treatment (ART) is threatened by the emergence of drug resistance mutations (DRM). Since Brazil presents the largest number of people living with HIV (PLWH) in South America we aimed at understanding the dynamics of DRM in this country. We analyzed a total of 20,226 HIV-1 sequences collected from PLWH undergoing ART between 2008–2017. Results show a mild decline of DRM over the years but an increase of the K65R reverse transcriptase mutation from 2.23% to 12.11%. This increase gradually occurred following alterations in the ART regimens replacing zidovudine (AZT) with tenofovir (TDF). PLWH harboring the K65R had significantly higher viral loads than those without this mutation (p < 0.001). Among the two most prevalent HIV-1 subtypes (B and C) there was a significant (p < 0.001) association of K65R with subtype C (11.26%) when compared with subtype B (9.27%). Nonetheless, evidence for K65R transmission in Brazil was found both for C and B subtypes. Additionally, artificial neural network-based immunoinformatic predictions suggest that K65R could enhance viral recognition by HLA-B27 that has relatively low prevalence in the Brazilian population. Overall, the results suggest that tenofovir-based regimens need to be carefully monitored particularly in settings with subtype C and specific HLA profiles.

Genetic diversity and drug resistance of HIV-1 among infected pregnant women newly diagnosed in Luanda, Angola

Monitoring genetic diversity and drug resistance mutations (DRMs) is critical for understanding HIV epidemiology. Here, we report HIV-1 genetic diversity and DRMs in blood samples from 42 HIV-positive pregnant women naive to antiretroviral therapy (ART), in Luanda. The samples were subjected to nested-PCR, followed by sequencing of HIV-1 pol gene, targeting the protease and reverse transcriptase fragments. HIV-1 diversity was analyzed using the REGA HIV-1 subtyping tool and DRMs were identified using the Calibrated Population Resistance tool. A total of 34 sequences were obtained. The data revealed wide HIV-1 subtypes heterogeneity, with subtype C (38%, 13/34) the most frequent, followed by the subtypes F1 (18%, 6/34), A1 (9%, 3/34), G (9%, 3/34), D (6%, 2/34) and H (3%, 1/34). In addition, recombinants strains were detected, with CRF02_AG (6%, 2/34) the most frequent, followed by CRF37_cpx, F1/C, A1/G and H/G, all with 3% (1/34). A total of 6/34 (18%) of the sequences presented DRMs. The non-nucleoside reverse transcriptase inhibitors presented 15% (5/34) of resistance. Moreover, 1/34 (3%) sequence presented resistance against both non-nucleoside reverse transcriptase inhibitors and nucleoside reverse transcriptase inhibitors, simultaneously. Despite the small sample size, our results suggest the need to update currently used ART regimens. Surveillance of HIV-1 subtypes and DRMs are necessary to understand HIV epidemiology and to guide modification of ART guidelines in Angola.

Monotypic low-level HIV viremias during antiretroviral therapy are associated with disproportionate production of X4 virions and systemic immune activation

OBJECTIVE

During effective antiretroviral therapy (ART), low-level plasma viremias (LLV) (HIV RNA >30-1000 copies/ml) can be detected intermittently. We hypothesized that systemic inflammation is associated with LLV either as the cause or result of the production of virions from clonally expanded cells.

METHODS

Prospective cohort study of HIV-infected ART-naive Peruvians enrolled prior to ART and followed for 2 years. Plasma HIV RNA and peripheral blood mononuclear cell (PBMC) HIV DNA concentrations were quantified pre-ART from individuals whose plasma HIV RNA was ART-suppressed. Inflammatory biomarker concentrations were measured pre and during ART. Single-genome amplification (SGA) derived HIV env and pol genotypes from pre-ART and LLV specimens. Antiretroviral levels during ART assessed adherence. Statistical associations and phylogenetic relationships were examined.

RESULTS

Among 82 participants with median plasma HIV RNA less than 30 copies/ml, LLV were detected in 33 of 82 (40%), with a LLV median HIV RNA of 73 copies/ml. Participants with vs. without LLV had significantly higher pre-ART plasma HIV RNA (P < 0.001) and PBMC HIV DNA (P < 0.007); but, during ART, their antiretroviral drug levels were similar. LLV env sequences were monotypic in 17 of 28 (61%) and diverse in 11 of 28 (39%) participants. Those with the monotypic vs. diverse LLV pattern had elevated hsCRP and sCD163 (P = 0.004) and LLV with more X4 variants (P = 0.02).

CONCLUSION

In individuals with monotypic LLV sequences, higher levels of pre-ART HIV DNA and RNA, systemic inflammation and X4 viruses suggest an interaction between inflammation and the production of virions from proliferating infected cells, and that naïve T cells may be a source of LLV.

KRDS: a web server for evaluating drug resistance mutations in kinases by molecular docking

Kinases are major targets of anti-cancer therapies owing to their importance in signaling processes that regulate cell growth and proliferation. However, drug resistance has emerged as a major obstacle to cancer therapy. Resistance to drugs has various underlying mechanisms, including the acquisition of mutations at drug binding sites and the resulting reduction in drug binding affinity. Therefore, the identification of mutations that are relevant to drug resistance may be useful to overcome this issue. We hypothesized that these mutations can be identified by combining recent advances in computational methods for protein structure modeling and ligand docking simulation. Hence, we developed a web-based tool named the Kinase Resistance Docking System (KRDS) that enables the assessment of the effects of mutations on kinase-ligand interactions. KRDS receives a list of mutations in kinases, generates structural models of the mutants, performs docking simulations, and reports the results to users. The changes in docking scores and docking conformations can be analyzed to infer the effects of mutations on drug binding and drug resistance. We expect our tool to improve our understanding of drug binding mechanisms and facilitate the development of effective new drugs to overcome resistance related to kinases; it may be particularly useful for biomedical researchers who are not familiar with computational environments. Our tool is available at http://bcbl.kaist.ac.kr/KRDS/.

The HIV-1 reverse transcriptase polymorphism A98S improves the response to tenofovir disoproxil fumarate+emtricitabine-containing HAART both in vivo and in vitro

The impact of baseline HIV-1 reverse transcriptase (RT) polymorphisms on response to first-line modern HAART containing tenofovir disoproxil fumarate (TDF) and emtricitabine (FTC) was evaluated. The impact of each RT polymorphism on virological success (VS) was evaluated in 604 HIV-1 subtype B-infected patients starting TDF+FTC-containing HAART. TDF and FTC antiviral activity was also tested in PBMCs infected by mutagenised HIV. Structural analysis based on docking simulations was performed. A98S was the only mutation significantly correlated with an increased proportion of patients achieving VS at 24 weeks (94.0% vs. 84.3%; P=0.03). Multivariate regression and Cox model analyses confirmed this result. At concentrations close to the minimal concentration achieved in patient plasma, TDF and FTC exhibited higher potency in the presence of A98S-mutated virus compared with wild-type (IC_90,TDF, 8.6±1.1 vs. 19.3±3.5nM; and IC_90,FTC, 12.4±7.7 vs. 16.8±9.8nM, respectively). The efficacy of FTC, abrogated by M184V, was partially restored by A98S (IC_90,FTC, 5169±5931nM for A98S+M184V vs. 18477±12478nM for M184V alone). Docking analysis showed the higher potency of TDF and FTC in the presence of A98S-mutated virus was mainly due to higher binding affinity between drugs and mutated RT compared with wild-type. In the presence of FTC, A98S also partially restored the RT binding affinity impaired by M184V alone. A98S polymorphism improves virological response to TDF+FTC-containing HAART. This may help clinicians in the choice of the optimal NRTI backbone aimed at achieving maximal virological inhibition.

In vitro resistance selection with doravirine (MK-1439), a novel nonnucleoside reverse transcriptase inhibitor with distinct mutation development pathways

Doravirine (DOR, formerly known as MK-1439) is a human immunodeficiency type 1 virus (HIV-1) nonnucleoside reverse transcriptase inhibitor (NNRTI) that is currently in phase 2b clinical trials. In vitro resistance selection of subtype B virus (MT4-green fluorescent protein [GFP] cells), as well as subtype A and C viruses (MT4-GFP/CCR5 cells) was conducted with DOR, rilpivirine (RPV), and efavirine (EFV) under low-multiplicity-of-infection conditions in a 96-well format. Resistance selection was performed with escalating concentrations of the NNRTIs ranging from the 95% effective concentration (1 × EC(95)) to 1,000 × EC(95) in the presence of 10% fetal bovine serum. In the resistance selection of subtype B virus with DOR, a V106A mutant virus led to two mutation pathways, followed by the emergence separately of either F227L or L234I. In the resistance selection of subtype A and C viruses, similar mutation development pathways were detected, in which a V106A or V106M mutant was also the starting virus in the pathways. Mutations that are commonly associated with RPV and EFV in clinical settings were also identified in subtype B viruses such as the E138K and K103N mutants, respectively, in this in vitro resistance selection study. The susceptibility of subtype B mutant viruses selected by DOR, RPV, and EFV to NNRTIs was evaluated. Results suggest that mutant viruses selected by DOR are susceptible to RPV and EFV and mutants selected by RPV and EFV are susceptible to DOR. When the replication capacity of the V106A mutant was compared with that of the wild-type (WT) virus, the mutant virus was 4-fold less fit than the WT virus.

Hepatitis C virus genetic variability and the presence of NS5B resistance-associated mutations as natural polymorphisms in selected genotypes could affect the response to NS5B inhibitors

Because of the extreme genetic variability of hepatitis C virus (HCV), we analyzed the NS5B polymerase genetic variability in circulating HCV genotypes/subtypes and its impact on the genetic barrier for the development of resistance to clinically relevant nucleoside inhibitors (NIs)/nonnucleoside inhibitors (NNIs). The study included 1,145 NS5B polymerase sequences retrieved from the Los Alamos HCV database and GenBank. The genetic barrier was calculated for drug resistance emergence. Prevalence and genetic barrier were calculated for 1 major NI and 32 NNI resistance variants (13 major and 19 minor) at 21 total NS5B positions. Docking calculations were used to analyze sofosbuvir affinity toward the diverse HCV genotypes. Overall, NS5B polymerase was moderately conserved among all HCV genotypes, with 313/591 amino acid residues (53.0%) showing ≤1% variability and 83/591 residues (14.0%) showing high variability (≥25.1%). Nine NNI resistance variants (2 major variants, 414L and 423I; 7 minor variants, 316N, 421V, 445F, 482L, 494A, 499A, and 556G) were found as natural polymorphisms in selected genotypes. In particular, 414L and 423I were found in HCV genotype 4 (HCV-4) (n = 14/38, 36.8%) and in all HCV-5 sequences (n = 17, 100%), respectively. Regardless of HCV genotype, the 282T major NI resistance variant and 10 major NNI resistance variants (316Y, 414L, 423I/T/V, 448H, 486V, 495L, 554D, and 559G) always required a single nucleotide substitution to be generated. Conversely, the other 3 major NNI resistance variants (414T, 419S, and 422K) were associated with a different genetic barrier score development among the six HCV genotypes. Sofosbuvir docking analysis highlighted a better ligand affinity toward HCV-2 than toward HCV-3, in agreement with the experimental observations. The genetic variability among HCV genotypes, particularly with the presence of polymorphisms at NNI resistance positions, could affect their responsiveness to NS5B inhibitors. A pretherapy HCV NS5B sequencing could help to provide patients with the full efficacy of NNI-containing regimens.

In vitro characterization of MK-1439, a novel HIV-1 nonnucleoside reverse transcriptase inhibitor

Nonnucleoside reverse transcriptase inhibitors (NNRTIs) are a mainstay of therapy for treating human immunodeficiency type 1 virus (HIV-1)-infected patients. MK-1439 is a novel NNRTI with a 50% inhibitory concentration (IC50) of 12, 9.7, and 9.7 nM against the wild type (WT) and K103N and Y181C reverse transcriptase (RT) mutants, respectively, in a biochemical assay. Selectivity and cytotoxicity studies confirmed that MK-1439 is a highly specific NNRTI with minimum off-target activities. In the presence of 50% normal human serum (NHS), MK-1439 showed excellent potency in suppressing the replication of WT virus, with a 95% effective concentration (EC95) of 20 nM, as well as K103N, Y181C, and K103N/Y181C mutant viruses with EC95 of 43, 27, and 55 nM, respectively. MK-1439 exhibited similar antiviral activities against 10 different HIV-1 subtype viruses (a total of 93 viruses). In addition, the susceptibility of a broader array of clinical NNRTI-associated mutant viruses (a total of 96 viruses) to MK-1439 and other benchmark NNRTIs was investigated. The results showed that the mutant profile of MK-1439 was superior overall to that of efavirenz (EFV) and comparable to that of etravirine (ETR) and rilpivirine (RPV). Furthermore, E138K, Y181C, and K101E mutant viruses that are associated with ETR and RPV were susceptible to MK-1439 with a fold change (FC) of <3. A two-drug in vitro combination study indicated that MK-1439 acts nonantagonistically in the antiviral activity with each of 18 FDA-licensed drugs for HIV infection. Taken together, these in vitro data suggest that MK-1439 possesses the desired properties for further development as a new antiviral agent.

Non-nucleoside reverse transcriptase inhibitor (NNRTI) cross-resistance: implications for preclinical evaluation of novel NNRTIs and clinical genotypic resistance testing

OBJECTIVES

The introduction of two new non-nucleoside reverse transcriptase inhibitors (NNRTIs) in the past 5 years and the identification of novel NNRTI-associated mutations have made it necessary to reassess the extent of phenotypic NNRTI cross-resistance.

METHODS

We analysed a dataset containing 1975, 1967, 519 and 187 genotype-phenotype correlations for nevirapine, efavirenz, etravirine and rilpivirine, respectively. We used linear regression to estimate the effects of RT mutations on susceptibility to each of these NNRTIs.

RESULTS

Sixteen mutations at 10 positions were significantly associated with the greatest contribution to reduced phenotypic susceptibility (≥10-fold) to one or more NNRTIs, including: 14 mutations at six positions for nevirapine (K101P, K103N/S, V106A/M, Y181C/I/V, Y188C/L and G190A/E/Q/S); 10 mutations at six positions for efavirenz (L100I, K101P, K103N, V106M, Y188C/L and G190A/E/Q/S); 5 mutations at four positions for etravirine (K101P, Y181I/V, G190E and F227C); and 6 mutations at five positions for rilpivirine (L100I, K101P, Y181I/V, G190E and F227C). G190E, a mutation that causes high-level nevirapine and efavirenz resistance, also markedly reduced susceptibility to etravirine and rilpivirine. K101H, E138G, V179F and M230L mutations, associated with reduced susceptibility to etravirine and rilpivirine, were also associated with reduced susceptibility to nevirapine and/or efavirenz.

CONCLUSIONS

The identification of novel cross-resistance patterns among approved NNRTIs illustrates the need for a systematic approach for testing novel NNRTIs against clinical virus isolates with major NNRTI-resistance mutations and for testing older NNRTIs against virus isolates with mutations identified during the evaluation of a novel NNRTI.

Molecular mechanism of antagonism between the Y181C and E138K mutations in HIV-1 reverse transcriptase

Etravirine (ETR) is an expanded-spectrum nonnucleoside reverse transcriptase inhibitor (NNRTI) approved for use as an antiretroviral agent in treatment-experienced patients. Y181C and E138K in HIV-1 RT are among 20 different drug resistance mutations associated with ETR. However, E138K can be consistently selected by ETR when wild-type viruses but not viruses containing Y181C are grown in tissue culture. This study was carried out to evaluate any possible mechanisms that might explain antagonism between the Y181C and E138K mutations. Accordingly, we performed tissue culture studies to investigate the evolutionary dynamics of E138K in both a wild-type (WT) and a Y181C background. We also generated recombinant enzymes containing Y181C and E138K alone or in combination in order to study enzyme processivity, rates of processive DNA synthesis, enzyme kinetics, and susceptibility to ETR. We now show that the presence of the Y181C mutation prevented the emergence of E138K in cell culture and that the simultaneous presence of E138K and Y181C impaired each of enzyme activity, processivity, rate of processive DNA synthesis, and deoxynucleoside triphosphate (dNTP) affinity. The addition of E138K to Y181C also decreased the level of resistance to ETR compared to that obtained with Y181C alone.