From Intercalation to External Binding: Ru(II) Complexes with a Spiro Ligand for TAR RNA Selective Binding and HIV-1 Reverse Transcriptase Inhibition.

As a typical RNA virus, the genetic information on HIV-1 is entirely stored in RNA. The reverse transcription activity of HIV-1 reverse transcriptase (RT) plays a crucial role in the replication and transmission of the virus. Non-nucleoside RT inhibitors (NNRTIs) block the function of RT by binding to the RNA binding site on RT, with very few targeting viral RNA. In this study, by transforming planar conjugated ligands into a spiro structure, we convert classical Ru(II) DNA intercalators into a nonintercalator. This enables selective binding to HIV-1 transactivation response (TAR) RNA on the outer side of nucleic acids through dual interactions involving hydrogen bonds and electrostatic attraction, effectively inhibiting HIV-1 RT and serving as a selective fluorescence probe for TAR RNA.

Determinants of the interaction between the 5'-leader of HIV-1 genome and human lysyl-tRNA synthetase in reverse transcription primer release process.

Reverse transcription of human immunodeficiency virus type 1 (HIV-1) initiates from the 3' end of human tRNALys3. The primer tRNALys3 is selectively packaged into the virus in the form of a complex with human lysyl-tRNA synthetase (LysRS). To facilitate reverse transcription initiation, part of the 5' leader (5'L) of HIV-1 genomic RNA (gRNA) evolves a tRNA anticodon-like element (TLE), which binds LysRS and releases tRNALys3 for primer annealing and reverse transcription initiation. Although TLE has been identified as a key element in 5'L responsible for LysRS binding, how the conformations and various hairpin structures of 5'L regulate 5'L-LysRS interaction is not fully understood. Here, these factors have been individually investigated using direct and competitive fluorescence anisotropy binding experiments. Our data showed that the conformation of 5'L significantly influences its binding affinity with LysRS. The 5'L conformation favoring gRNA dimerization and packaging exhibits much weaker binding affinity with LysRS compared to the alternative 5'L conformation that is not selected for packaging. Additionally, dimerization of 5'L impairs LysRS-5'L interaction. Furthermore, among various regions of 5'L, both the primer binding site/TLE domain and the stem-loop 3 are important for LysRS interaction, whereas the dimerization initiation site and the splicing donor plays a minor role. In contrast, the presence of the transacting responsive and the polyadenylation signal hairpins slightly inhibit LysRS binding. These findings reveal that the conformation and various regions of the 5'L of HIV-1 genome regulate its interaction with human LysRS and the reverse transcription primer release process.

Interplay between protease and reverse transcriptase dimerization in a model HIV-1 polyprotein.

The Gag-Pol polyprotein in human immunodeficiency virus type I (HIV-1) encodes enzymes that are essential for virus replication: protease (PR), reverse transcriptase (RT), and integrase (IN). The mature forms of PR, RT and IN are homodimer, heterodimer and tetramer, respectively. The precise mechanism underlying the formation of dimer or tetramer is not yet understood. Here, to gain insight into the dimerization of PR and RT in the precursor, we prepared a model precursor, PR-RT, incorporating an inactivating mutation at the PR active site, D25A, and including two residues in the p6* region, fused to a SUMO-tag, at the N-terminus of the PR region. We also prepared two mutants of PR-RT containing a dimer dissociation mutation either in the PR region, PR(T26A)-RT, or in the RT region, PR-RT(W401A). Size exclusion chromatography showed both monomer and dimer fractions in PR-RT and PR(T26A)-RT, but only monomer in PR-RT(W401A). SEC experiments of PR-RT in the presence of protease inhibitor, darunavir, significantly enhanced the dimerization. Additionally, SEC results suggest an estimated PR-RT dimer dissociation constant that is higher than that of the mature RT heterodimer, p66/p51, but slightly lower than the premature RT homodimer, p66/p66. Reverse transcriptase assays and RT maturation assays were performed as tools to assess the effects of the PR dimer-interface on these functions. Our results consistently indicate that the RT dimer-interface plays a crucial role in the dimerization in PR-RT, whereas the PR dimer-interface has a lesser role.

Premature Activation of the HIV-1 Protease Is Influenced by Polymorphisms in the Hinge Region.

HIV-1 protease inhibitors are an essential component of antiretroviral therapy. However, drug resistance is a pervasive issue motivating a persistent search for novel therapies. Recent reports found that when protease activates within the host cell's cytosol, it facilitates the pyroptotic killing of infected cells. This has led to speculation that promoting protease activation, rather than inhibiting it, could help to eradicate infected cells and potentially cure HIV-1 infection. Here, we used a nanoscale flow cytometry-based assay to characterize protease resistance mutations and polymorphisms. We quantified protease activity, viral concentration, and premature protease activation and confirmed previous findings that major resistance mutations generally destabilize the protease structure. Intriguingly, we found evidence that common polymorphisms in the hinge domain of protease can influence its susceptibility to premature activation. This suggests that viral heterogeneity could pose a considerable challenge for therapeutic strategies aimed at inducing premature protease activation in the future.

Performance of the Applied Biosystems HIV-1 Genotyping Kit with Integrase.

HIV genotyping is used to assess HIV susceptibility to antiretroviral drugs. The Applied Biosystems HIV-1 Genotyping Kit with Integrase (AB kit, Thermo Fisher Scientific) detects resistance-associated mutations (RAMs) in HIV protease (PR), reverse transcriptase (RT), and integrase (IN). We compared results from the AB kit with results obtained previously with the ViroSeq HIV-1 Genotyping System. DNA amplicons from the AB kit were also analyzed using next-generation sequencing (NGS). HIV RNA was extracted using the MagNA Pure 24 instrument (Roche Diagnostics; 96 plasma samples, HIV subtype B, viral load range: 530-737,741 copies/mL). FASTA files were generated from AB kit data using Exatype (Hyrax Biosciences). DNA amplicons from the AB kit were also analyzed by NGS using the Nextera XT kit (Illumina). Drug resistance was predicted using the Stanford HIV Drug Resistance Database. The mean genetic distance for sequences from ViroSeq and the AB kit was 0.02% for PR/RT and 0.04% for IN; 103 major RAMs were detected by both methods. Four additional major RAMs were detected by the AB kit only. These four major RAMs were also detected by NGS (detected in 18.1%-38.2% of NGS reads). NGS detected 27 major RAMs that were not detected with either of the Sanger sequencing-based kits. All major RAMs detected with ViroSeq were detected with the AB kit; additional RAMs were detected with the AB kit only. DNA amplicons from the AB kit can be used for NGS for more sensitive detection of RAMs.

Design, synthesis and biological evaluation of Thiazolo[3, 2-a]Pyrimidine derivatives as novel RNase H inhibitors.

Targeting Ribonuclease H (RNase H) has been considered a viable strategy for HIV therapy. In this study, a series of novel thiazolo[3, 2-a]pyrimidine derivatives were firstly designed and synthesized as potential inhibitors of HIV-1 RNase H. Among these compounds, A28 exhibited the most potent inhibition against HIV-1 RNase H with an IC50 value of 4.14 µM, which was about 5-fold increase in potency than the hit compound A1 (IC50 = 21.49 µM). To gain deeper insights into the structure-activity relationship (SAR), a CoMFA model was constructed to yield reasonable statistical results (q2 = 0.658 and R2 = 0.969). Results from magnesium ion chelation experiments and molecular docking studies revealed that these thiazolopyrimidine inhibitors may exert their inhibitory activity by binding to an allosteric site on RNase H at the interface between subunits p51 and p66. Furthermore, this analog demonstrated favorable physicochemical properties. Our findings provide valuable groundwork for further development of allosteric inhibitors targeting HIV-1 RNase H.

Comparison of feline and human immunodeficiency virus reverse transcriptase enzymes through chemical screening and computational analysis.

Feline immunodeficiency virus (FIV) is a common infection found in domesticated and wild cats worldwide. Despite the wealth of therapeutic understanding of the disease in humans, considerably less information exists regarding the treatment of the disease in felines. Current treatment relies on drugs developed for the related human immunodeficiency virus (HIV) and includes compounds of the popular non-nucleotide reverse transcriptase (NNRTI) class. This is despite FIV-RT being only 67% similar to HIV-1 RT at the enzyme level, increasing to 88% for the allosteric pocket targeted by NNRTIs. The goal of this project was to try to quantify how well the more extensive pharmacological knowledge available for human disease translates to felines. To this end we screened known NNRTIs and 10 diverse pyrimidine analogs identified virtually. We use this chemo-centric probe approach to (a) assess the similarity between the two related RT targets based on the observed experimental inhibition values, (b) try to identify more potent inhibitors at FIV, and (c) gain a better appreciation of the structure-activity relationships (SAR). We found the correlation between IC50s at the two targets to be strong (r2 = 0.87) and identified compound 1 as the most potent inhibitor of FIV with IC50 of 0.030 µM ± 0.009. This compared to FIV IC50 values of 0.22 ± 0.17 µM, 0.040 ± 0.010 µM and >160 µM for known anti HIV-1 RT drugs Efavirenz, Rilpivirine, and Nevirapine, respectively. This knowledge, along with an understanding of the structural origin that give rise to any differences could improve the way HIV drugs are repurposed for FIV.

Structure-Based Design of Novel Thiazolone[3,2-a]pyrimidine Derivatives as Potent RNase H Inhibitors for HIV Therapy.

Ribonuclease H (RNase H) was identified as an important target for HIV therapy. Currently, no RNase H inhibitors have reached clinical status. Herein, a series of novel thiazolone[3,2-a]pyrimidine-containing RNase H inhibitors were developed, based on the hit compound 10i, identified from screening our in-house compound library. Some of these derivatives exhibited low micromolar inhibitory activity. Among them, compound 12b was identified as the most potent inhibitor of RNase H (IC50 = 2.98 µM). The experiment of magnesium ion coordination was performed to verify that this ligand could coordinate with magnesium ions, indicating its binding ability to the catalytic site of RNase H. Docking studies revealed the main interactions of this ligand with RNase H. A quantitative structure activity relationship (QSAR) was also conducted to disclose several predictive mathematic models. A molecular dynamics simulation was also conducted to determine the stability of the complex. Taken together, thiazolone[3,2-a]pyrimidine can be regarded as a potential scaffold for the further development of RNase H inhibitors.

Structure-guided design of novel biphenyl-quinazoline derivatives as potent non-nucleoside reverse transcriptase inhibitors featuring improved anti-resistance, selectivity, and solubility.

In pursuit of enhancing the anti-resistance efficacy and solubility of our previously identified NNRTI 1, a series of biphenyl-quinazoline derivatives were synthesized employing a structure-based drug design strategy. Noteworthy advancements in anti-resistance efficacy were discerned among some of these analogs, prominently exemplified by compound 7ag, which exhibited a remarkable 1.37 to 602.41-fold increase in potency against mutant strains (Y181C, L100I, Y188L, F227L + V106A, and K103N + Y181C) in comparison to compound 1. Compound 7ag also demonstrated comparable anti-HIV activity against both WT HIV and K103N, albeit with a marginal reduction in activity against E138K. Of significance, this analog showed augmented selectivity index (SI > 5368) relative to compound 1 (SI > 37764), Nevirapine (SI > 158), Efavirenz (SI > 269), and Etravirine (SI > 1519). Moreover, it displayed a significant enhancement in water solubility, surpassing that of compound 1, Etravirine, and Rilpivirine. To elucidate the underlying molecular mechanisms, molecular docking studies were undertaken to probe the critical interactions between 7ag and both WT and mutant strains of HIV-1 RT. These findings furnish invaluable insights driving further advancements in the development of DAPYs for HIV therapy.

HIV-1 Integrase Assembles Multiple Species of Stable Synaptic Complex Intasomes That Are Active for Concerted DNA Integration In vitro.

Retroviral DNA integration is mediated by nucleoprotein complexes (intasomes) in which a pair of viral DNA ends are bridged by a multimer of integrase (IN). Most of the high-resolution structures of HIV-1 intasomes are based on an HIV-1 IN with an Sso7d protein domain fused to the N-terminus. Sso7d-IN aggregates much less than wild-type IN and has been critical for structural studies of HIV-1 intasomes. Unexpectedly, these structures revealed that the common core architecture that mediates catalysis could be assembled in various ways, giving rise to both tetrameric and dodecameric intasomes, together with other less well-characterized species. This differs from related retroviruses that assemble unique multimeric intasomes, although the number of protomers in the intasome varies between viruses. The question of whether the additional Sso7d domain contributes to the heterogeneity of HIV-1 intasomes is therefore raised. We have addressed this by biochemical and structural studies of intasomes assembled with wild-type HIV-1 IN. Negative stain and cryo-EM reveal a similar range of multimeric intasome species as with Sso7d-IN with the same common core architecture. Stacks of intasomes resulting from domain swapping are also seen with both wild-type and Sso7d-IN intasomes. The propensity to assemble multimeric intasome species is, therefore, an intrinsic property of HIV-1 IN and is not conferred by the presence of the Sso7d domain. The recently solved intasome structures of different retroviral species, which have been reported to be tetrameric, octameric, dodecameric, and hexadecameric, highlight how a common intasome core architecture can be assembled in different ways for catalysis.

Role of I182, R187, and K188 Amino Acid Residues in the Catalytic Domain of HIV-1 Integrase in the Processes of Reverse Transcription and Integration.

Structural organization of HIV-1 integrase is based on a tetramer formed by two protein dimers. Within this tetramer, the catalytic domain of one subunit of the first dimer interacts with the N-terminal domain of the second dimer subunit. It is the tetrameric structure that allows both ends of the viral DNA to be correctly positioned relative to the cellular DNA and to realize catalytic functions of integrase, namely 3'-processing and strand transfer. However, during the HIV-1 replicative cycle, integrase is responsible not only for the integration stage, it is also involved in reverse transcription and is necessary at the stage of capsid formation of the newly formed virions. It has been suggested that HIV-1 integrase is a structurally dynamic protein and its biological functions depend on its structure. Accordingly, studying interactions between the domains of integrase that provide its tetrameric structure is important for understanding its multiple functions. In this work, we investigated the role of three amino acids of the catalytic domain, I182, R187, and K188, located in the contact region of two integrase dimers in the tetramer structure, in reverse transcription and integration. It has been shown that the R187 residue is extremely important for formation of the correct integrase structure, which is necessary at all stages of its functional activity. The I182 residue is necessary for successful integration and is not important for reverse transcription, while the K188 residue, on the contrary, is involved in formation of the integrase structure, which is important for the effective reverse transcription.

Prevalence and Phenotypic Susceptibility to Doravirine of the HIV-1 Reverse Transcriptase V106I Polymorphism in B and Non-B Subtypes.

BACKGROUND: Limited data are available regarding the susceptibility of the reverse transcriptase V106 polymorphism to doravirine. METHODS: Doravirine susceptibility was measured in site-directed mutants (SDMs) containing V106I, V106A, V106M, and Y188L mutations in subtype B (NL4-3, HXB2) and CRF02_AG background and in recombinant viruses with RT harboring V106I alone derived from 50 people with HIV. RESULTS: HIV-1 B subtype was detected in 1523 of 2705 cases. Prevalence of V106I was 3.2% in B and 2.5% in non-B subtypes, and was higher in subtype F (8.1%) and D (14.3%). Fold-changes (FC) in susceptibility for SDMs were below doravirine biological cutoff (3.0) for V106I, but not for V106A, V106M, and Y188L. Clinically derived viruses tested included 22 B (median FC, 1.2; interquartile range [IQR], 0.9-1.6) and 28 non-B subtypes (median FC, 1.8; IQR, 0.9-3.0). Nine (18%) viruses showed FC values equal or higher than the doravirine biological FC cutoff. CONCLUSIONS: The prevalence of the HIV-1 RT V106I polymorphism in MeditRes HIV consortium remains low, but significantly more prevalent in subtypes D and F. V106I minimally decreased the susceptibility to doravirine in SDMs and most clinical isolates. Reduced susceptibility seems to occur at increased frequency in subtype F1; however, the clinical impact remains to be investigated. CLINICAL TRIALS REGISTRATION: NCT04894357.

HIV-1 reverse transcriptase stability correlates with Gag cleavage efficiency: reverse transcriptase interaction implications for modulating protease activation.

Proteolytic processing of human immunodeficiency virus type 1 particles mediated by viral protease (PR) is essential for acquiring virus infectivity. Activation of PR embedded in Gag-Pol is triggered by Gag-Pol dimerization during virus assembly. We previously reported that amino acid substitutions at the RT tryptophan repeat motif destabilize virus-associated RT and attenuate the ability of efavirenz (EFV, an RT dimerization enhancer) to increase PR-mediated Gag cleavage efficiency. Furthermore, a single amino acid change at RT significantly reduces virus yields due to enhanced Gag cleavage. These data raise the possibility of the RT domain contributing to PR activation by promoting Gag-Pol dimerization. To test this hypothesis, we investigated the putative involvement of a hydrophobic leucine repeat motif (LRM) spanning RT L282 to L310 in RT/RT interactions. We found that LRM amino acid substitutions led to RT instability and that RT is consequently susceptible to degradation by PR. The LRM mutants exhibited reduced Gag cleavage efficiencies while attenuating the EFV enhancement of Gag cleavage. In addition, an RT dimerization-defective mutant, W401A, reduced enhanced Gag cleavage via a leucine zipper (LZ) motif inserted at the deleted Gag-Pol region. Importantly, the presence of RT and integrase domains failed to counteract the LZ enhancement of Gag cleavage. A combination of the Gag cleavage enhancement factors EFV and W402A markedly impaired Gag cleavage, indicating a disruption of W402A Gag-Pol dimerization following EFV binding to W402A Gag-Pol. Our results support the idea that RT modulates PR activation by affecting Gag-Pol/Gag-Pol interaction. IMPORTANCE A stable reverse transcriptase (RT) p66/51 heterodimer is required for HIV-1 genome replication in host cells following virus entry. The activation of viral protease (PR) to mediate virus particle processing helps viruses acquire infectivity following cell release. RT and PR both appear to be major targets for inhibiting HIV-1 replication. We found a strong correlation between impaired p66/51RT stability and deficient PR-mediated Gag cleavage, suggesting that RT/RT interaction is critical for triggering PR activation via the promotion of adequate Gag-Pol dimerization. Accordingly, RT/RT interaction is a potentially advantageous method for anti-HIV/AIDS therapy if it is found to simultaneously block PR and RT enzymatic activity.

Prediction of HIV-1 protease cleavage site from octapeptide sequence information using selected classifiers and hybrid descriptors.

BACKGROUND: In most parts of the world, especially in underdeveloped countries, acquired immunodeficiency syndrome (AIDS) still remains a major cause of death, disability, and unfavorable economic outcomes. This has necessitated intensive research to develop effective therapeutic agents for the treatment of human immunodeficiency virus (HIV) infection, which is responsible for AIDS. Peptide cleavage by HIV-1 protease is an essential step in the replication of HIV-1. Thus, correct and timely prediction of the cleavage site of HIV-1 protease can significantly speed up and optimize the drug discovery process of novel HIV-1 protease inhibitors. In this work, we built and compared the performance of selected machine learning models for the prediction of HIV-1 protease cleavage site utilizing a hybrid of octapeptide sequence information comprising bond composition, amino acid binary profile (AABP), and physicochemical properties as numerical descriptors serving as input variables for some selected machine learning algorithms. Our work differs from antecedent studies exploring the same subject in the combination of octapeptide descriptors and method used. Instead of using various subsets of the dataset for training and testing the models, we combined the dataset, applied a 3-way data split, and then used a "stratified" 10-fold cross-validation technique alongside the testing set to evaluate the models. RESULTS: Among the 8 models evaluated in the "stratified" 10-fold CV experiment, logistic regression, multi-layer perceptron classifier, linear discriminant analysis, gradient boosting classifier, Naive Bayes classifier, and decision tree classifier with AUC, F-score, and B. Acc. scores in the ranges of 0.91-0.96, 0.81-0.88, and 80.1-86.4%, respectively, have the closest predictive performance to the state-of-the-art model (AUC 0.96, F-score 0.80 and B. Acc. ~ 80.0%). Whereas, the perceptron classifier and the K-nearest neighbors had statistically lower performance (AUC 0.77-0.82, F-score 0.53-0.69, and B. Acc. 60.0-68.5%) at p < 0.05. On the other hand, logistic regression, and multi-layer perceptron classifier (AUC of 0.97, F-score > 0.89, and B. Acc. > 90.0%) had the best performance on further evaluation on the testing set, though linear discriminant analysis, gradient boosting classifier, and Naive Bayes classifier equally performed well (AUC > 0.94, F-score > 0.87, and B. Acc. > 86.0%). CONCLUSIONS: Logistic regression and multi-layer perceptron classifiers have comparable predictive performances to the state-of-the-art model when octapeptide sequence descriptors consisting of AABP, bond composition and standard physicochemical properties are used as input variables. In our future work, we hope to develop a standalone software for HIV-1 protease cleavage site prediction utilizing the linear regression algorithm and the aforementioned octapeptide sequence descriptors.

Computational and Crystallographic Analysis of Binding Structures of Inhibitory Compounds for HIV-1 RNase H Activity.

Chemotherapy of human immunodeficiency virus type-1 (HIV-1) has significantly developed over the last three decades. The emergence of drug-resistant variants is, however, still a severe problem. The RNase H activity of HIV-1 reverse transcriptase is an attractive target for a new class of antiviral drugs because there is no approved inhibitor. The nitro-furan-carbonyl and nitro-thiophene-carbonyl groups are potent scaffolds for the HIV-1 RNase H inhibitor. In this work, the binding structures of six inhibitory compounds were obtained by X-ray crystal analysis in a complex with a recombinant protein of HIV-1 RNase H domain. Every inhibitory compound was found to be bound to the catalytic site with the furan- or thiophene-ring coordinated to two divalent metal ions at the binding pocket. All the atoms in nitro, furan, carbonyl, and two metals were aligned in the nitro-furan derivatives. The straight line connecting nitro and carboxyl groups was parallel to the plane made by two metal ions and a furan O atom. The binding modes of the nitro-thiophene derivatives were slightly different from those of the nitro-furan ones. The nitro and carbonyl groups deviated from the plane made by two metals and a thiophene S atom. Molecular dynamics simulations suggested that the furan O or thiophene S atom and carbonyl O atom were firmly coordinated to the metal ions. The simulations made the planar nitro-furan moiety well aligned to the line connecting the two metal ions. In contrast, the nitro-thiophene derivatives were displaced from the initial positions after the simulations. The computational findings will be a sound basis for developing potent inhibitors for HIV-1 RNase H activity.

Effectively predicting HIV-1 protease cleavage sites by using an ensemble learning approach.

BACKGROUND: The site information of substrates that can be cleaved by human immunodeficiency virus 1 proteases (HIV-1 PRs) is of great significance for designing effective inhibitors against HIV-1 viruses. A variety of machine learning-based algorithms have been developed to predict HIV-1 PR cleavage sites by extracting relevant features from substrate sequences. However, only relying on the sequence information is not sufficient to ensure a promising performance due to the uncertainty in the way of separating the datasets used for training and testing. Moreover, the existence of noisy data, i.e., false positive and false negative cleavage sites, could negatively influence the accuracy performance. RESULTS: In this work, an ensemble learning algorithm for predicting HIV-1 PR cleavage sites, namely EM-HIV, is proposed by training a set of weak learners, i.e., biased support vector machine classifiers, with the asymmetric bagging strategy. By doing so, the impact of data imbalance and noisy data can thus be alleviated. Besides, in order to make full use of substrate sequences, the features used by EM-HIV are collected from three different coding schemes, including amino acid identities, chemical properties and variable-length coevolutionary patterns, for the purpose of constructing more relevant feature vectors of octamers. Experiment results on three independent benchmark datasets demonstrate that EM-HIV outperforms state-of-the-art prediction algorithm in terms of several evaluation metrics. Hence, EM-HIV can be regarded as a useful tool to accurately predict HIV-1 PR cleavage sites.

Cryo-EM structures of wild-type and E138K/M184I mutant HIV-1 RT/DNA complexed with inhibitors doravirine and rilpivirine.

Structures trapping a variety of functional and conformational states of HIV-1 reverse transcriptase (RT) have been determined by X-ray crystallography. These structures have played important roles in explaining the mechanisms of catalysis, inhibition, and drug resistance and in driving drug design. However, structures of several desired complexes of RT could not be obtained even after many crystallization or crystal soaking experiments. The ternary complexes of doravirine and rilpivirine with RT/DNA are such examples. Structural study of HIV-1 RT by single-particle cryo-electron microscopy (cryo-EM) has been challenging due to the enzyme's relatively smaller size and higher flexibility. We optimized a protocol for rapid structure determination of RT complexes by cryo-EM and determined six structures of wild-type and E138K/M184I mutant RT/DNA in complexes with the nonnucleoside inhibitors rilpivirine, doravirine, and nevirapine. RT/DNA/rilpivirine and RT/DNA/doravirine complexes have structural differences between them and differ from the typical conformation of nonnucleoside RT inhibitor (NNRTI)-bound RT/double-stranded DNA (dsDNA), RT/RNA-DNA, and RT/dsRNA complexes; the primer grip in RT/DNA/doravirine and the YMDD motif in RT/DNA/rilpivirine have large shifts. The DNA primer 3'-end in the doravirine-bound structure is positioned at the active site, but the complex is in a nonproductive state. In the mutant RT/DNA/rilpivirine structure, I184 is stacked with the DNA such that their relative positioning can influence rilpivirine in the pocket. Simultaneously, E138K mutation opens the NNRTI-binding pocket entrance, potentially contributing to a faster rate of rilpivirine dissociation by E138K/M184I mutant RT, as reported by an earlier kinetic study. These structural differences have implications for understanding molecular mechanisms of drug resistance and for drug design.

The HIV-1 Viral Protease Is Activated during Assembly and Budding Prior to Particle Release.

HIV-1 encodes a viral protease that is essential for the maturation of infectious viral particles. While protease inhibitors are effective antiretroviral agents, recent studies have shown that prematurely activating, rather than inhibiting, protease function leads to the pyroptotic death of infected cells, with exciting implications for efforts to eradicate viral reservoirs. Despite 40 years of research into the kinetics of protease activation, it remains unclear exactly when protease becomes activated. Recent reports have estimated that protease activation occurs minutes to hours after viral release, suggesting that premature protease activation is challenging to induce efficiently. Here, monitoring viral protease activity with sensitive techniques, including nanoscale flow cytometry and instant structured illumination microscopy, we demonstrate that the viral protease is activated within cells prior to the release of free virions. Using genetic mutants that lock protease into a precursor conformation, we further show that both the precursor and mature protease have rapid activation kinetics and that the activity of the precursor protease is sufficient for viral fusion with target cells. Our finding that HIV-1 protease is activated within producer cells prior to release of free virions helps resolve a long-standing question of when protease is activated and suggests that only a modest acceleration of protease activation kinetics is required to induce potent and specific elimination of HIV-infected cells. IMPORTANCE HIV-1 protease inhibitors have been a mainstay of antiretroviral therapy for more than 2 decades. Although antiretroviral therapy is effective at controlling HIV-1 replication, persistent reservoirs of latently infected cells quickly reestablish replication if therapy is halted. A promising new strategy to eradicate the latent reservoir involves prematurely activating the viral protease, which leads to the pyroptotic killing of infected cells. Here, we use highly sensitive techniques to examine the kinetics of protease activation during and shortly after particle formation. We found that protease is fully activated before virus is released from the cell membrane, which is hours earlier than recent estimates. Our findings help resolve a long-standing debate as to when the viral protease is initially activated during viral assembly and confirm that prematurely activating HIV-1 protease is a viable strategy to eradicate infected cells following latency reversal.

Prediction and molecular field view of drug resistance in HIV-1 protease mutants.

Conquering the mutational drug resistance is a great challenge in anti-HIV drug development and therapy. Quantitatively predicting the mutational drug resistance in molecular level and elucidating the three dimensional structure-resistance relationships for anti-HIV drug targets will help to improve the understanding of the drug resistance mechanism and aid the design of resistance evading inhibitors. Here the MB-QSAR (Mutation-dependent Biomacromolecular Quantitative Structure Activity Relationship) method was employed to predict the molecular drug resistance of HIV-1 protease mutants towards six drugs, and to depict the structure resistance relationships in HIV-1 protease mutants. MB-QSAR models were constructed based on a published data set of Ki values for HIV-1 protease mutants against drugs. Reliable MB-QSAR models were achieved and these models display both well internal and external prediction abilities. Interpreting the MB-QSAR models supplied structural information related to the drug resistance as well as the guidance for the design of resistance evading drugs. This work showed that MB-QSAR method can be employed to predict the resistance of HIV-1 protease caused by polymorphic mutations, which offer a fast and accurate method for the prediction of other drug target within the context of 3D structures.

Drug Resistance in the HIV-1 Subtype C Protease Enzyme: A High Throughput Virtual Screening Approach in Search of New Ligands with Activity.

BACKGROUND: HIV-1 subtype C protease is a strategic target for antiretroviral treatment. However, resistance to protease inhibitors appears after months of treatment. Chromones and 2- biscoumarin derivatives show potential for inhibition of the HIV- subtype C protease. OBJECTIVE: Different heterocyclic structures from the ZINC database were docked against Human Immunodeficiency Virus-1 (HIV) subtype C protease crystal structure 2R5Q and 2R5P. The 5 best molecules were selected to be docked against 62 homology models based on HIV-protease sequences from infants failing antiretroviral protease treatment. This experimentation was performed with two molecular docking programs: Autodock and Autodock Vina. These molecules were modified by substituting protons with different moieties, and the derivatives were docked against the same targets. Ligand-protein interactions, physical/chemical proprieties of the molecules, and dynamics simulations were analyzed. METHODS: Docking of all of the molecules was performed to find out the binding sites of HIV-1 subtype C proteases. An in-house script was made to substitute protons of molecules with different moieties. According to the Lipinski rule of five, physical and chemical properties were determined. Complexes of certain ligands-protease were compared to the protein alone in molecular dynamics simulations. RESULTS: From the first docking results, the 5 best (lowest energy) ligands (dibenz[a,h]acridine, dibenz[a, i]acridine, NSC114903, dibenz[c,h]acridine, benzo[a]acridine) were selected. The binding energy of the modified ligands increased, including the poorest-performing molecules. A correlation between nature, the position, and the resulting binding energy was observed. According to the Lipinski rules, the physico-chemical characteristics of the five best-modified ligands are ideal for oral bioavailability. Molecular dynamics simulations show that some lead-protease complexes were stable. CONCLUSION: Dibenz[a,h]acridine, dibenz[a, i]acridine, NSC114903, dibenz[c,h]acridine, benzo[ a]acridine and their derivatives might be considered as promising HIV-1 subtype C protease inhibitors. This could be confirmed through synthesis and subsequent in vitro assays.