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

Computer-Aided Prediction of the Interactions of Viral Proteases with Antiviral Drugs: Antiviral Potential of Broad-Spectrum Drugs

Human society is facing the threat of various viruses. Proteases are promising targets for the treatment of viral infections. In this study, we collected and profiled 170 protease sequences from 125 viruses that infect humans. Approximately 73 of them are viral 3-chymotrypsin-like proteases (3CLpro), and 11 are pepsin-like aspartic proteases (PAPs). Their sequences, structures, and substrate characteristics were carefully analyzed to identify their conserved nature for proposing a pan-3CLpro or pan-PAPs inhibitor design strategy. To achieve this, we used computational prediction and modeling methods to predict the binding complex structures for those 73 3CLpro with 4 protease inhibitors of SARS-CoV-2 and 11 protease inhibitors of HCV. Similarly, the complex structures for the 11 viral PAPs with 9 protease inhibitors of HIV were also obtained. The binding affinities between these compounds and proteins were also evaluated to assess their pan-protease inhibition via MM-GBSA. Based on the drugs targeting viral 3CLpro and PAPs, repositioning of the active compounds identified several potential uses for these drug molecules. As a result, Compounds 1-2, modified based on the structures of Ray1216 and Asunaprevir, indicate potential inhibition of DENV protease according to our computational simulation results. These studies offer ideas and insights for future research in the design of broad-spectrum antiviral drugs.

Genetic Diversity and Low Therapeutic Impact of Variant-Specific Markers in HIV-1 Pol Proteins

The emergence and spread of new HIV-1 variants pose a challenge for the effectiveness of antiretrovirals (ARV) targeting Pol proteins. During viral evolution, non-synonymous mutations have fixed along the viral genome, leading to amino acid (aa) changes that can be variant-specific (V-markers). Those V-markers fixed in positions associated with drug resistance mutations (DRM), or R-markers, can impact drug susceptibility and resistance pathways. All available HIV-1 Pol sequences from ARV-naïve subjects were downloaded from the United States Los Alamos HIV Sequence Database, selecting 59,733 protease (PR), 6,437 retrotranscriptase (RT), and 6,059 integrase (IN) complete sequences ascribed to the four HIV-1 groups and group M subtypes and circulating recombinant forms (CRFs). Using a bioinformatics tool developed in our laboratory (EpiMolBio), we inferred the consensus sequences for each Pol protein and HIV-1 variant to analyze the aa conservation in Pol. We analyzed the Wu–Kabat protein variability coefficient (WK) in PR, RT, and IN group M to study the susceptibility of each site to evolutionary replacements. We identified as V-markers the variant-specific aa changes present in >75% of the sequences in variants with >5 available sequences, considering R-markers those V-markers that corresponded to DRM according to the IAS-USA2019 and Stanford-Database 9.0. The mean aa conservation of HIV-1 and group M consensus was 82.60%/93.11% in PR, 88.81%/94.07% in RT, and 90.98%/96.02% in IN. The median group M WK was 10 in PR, 4 in RT, and 5 in IN. The residues involved in binding or catalytic sites showed a variability <0.5%. We identified 106 V-markers: 31 in PR, 28 in RT, and 47 in IN, present in 11, 12, and 13 variants, respectively. Among them, eight (7.5%) were R-markers, present in five variants, being minor DRM with little potential effect on ARV susceptibility. We present a thorough analysis of Pol variability among all HIV-1 variants circulating to date. The relatively high aa conservation observed in Pol proteins across HIV-1 variants highlights their critical role in the viral cycle. However, further studies are needed to understand the V-markers’ impact on the Pol proteins structure, viral cycle, or treatment strategies, and periodic variability surveillance studies are also required to understand PR, RT, and IN evolution.

Structural and functional insights into the DNA damage-inducible protein 1 (Ddi1) from protozoa

Ddi1 is a multidomain protein that belongs to the ubiquitin receptor family of proteins. The Ddi1 proteins contain a highly conserved retroviral protease (RVP)-like domain along with other domains. The severity of opportunistic infections, caused by parasitic protozoa in AIDS patients, was found to decline when HIV protease inhibitors were used in antiretroviral therapy. Parasite growth was shown to be suppressed by a few of the inhibitors targeting Ddi1 present in these parasites. In this study, the binding of HIV protease inhibitors to the RVP domain of Ddi1 from Toxoplasma gondii and Cryptosporidium hominis; and the binding of ubiquitin to the ubiquitin-associated domain of Ddi1 from these two parasites were established using Biolayer Interferometry. The crystal structures of the RVP domains of Ddi1 from T. gondii and C. hominis were determined; they form homodimers similar to those observed in HIV protease and the reported structures of the same domain from Saccharomyces cerevisiae, Leishmania major and humans. The native form of the domain showed an open dimeric structure and a normal mode analysis revealed that it can take up a closed conformation resulting from relative movements of the subunits. Based on the crystal structure of the RVP domain of Ddi1 from L. major, a seven residue peptide inhibitor was designed and it was shown to bind to the RVP domain of Ddi1 from L. major by Biolayer Interferometry. This peptide was modified using computational methods and was shown to have a better affinity than the initial peptide.

•

Crystal structures of the retroviral protease (RVP)-like domains of DNA damage inducible protein 1 (Ddi1) from Toxoplasma gondii and Cryptosporidium hominis have been determined.

•

ToxoDdi1-RVP and CrypDdi1-RVP crystal structures illustrate the conserved structural features with HIV protease.

•

HIV protease inhibitors were found to bind to ToxoDdi1-RVP and CrypDdi1-RVP in solution showing the potential to be developed as drugs for protozoal diseases.

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 K_i 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.

Identification of Highly Potent Human Immunodeficiency Virus Type-1 Protease Inhibitors against Lopinavir and Darunavir Resistant Viruses from Allophenylnorstatine-Based Peptidomimetics with P2 Tetrahydrofuranylglycine

The emergence of drug-resistant HIV from a widespread antiviral chemotherapy targeting HIV protease in the past decades is unavoidable and provides a challenge to develop alternative inhibitors. We synthesized a series of allophenylnorstatine-based peptidomimetics with various P₃, P₂, and P₂́ moieties. The derivatives with P₂ tetrahydrofuranylglycine (Thfg) were found to be potent against wild type HIV-1 protease and the virus, leading to a highly potent compound 21f (KNI-1657) against lopinavir/ritonavir- or darunavir-resistant strains. Co-crystal structures of 21f and the wild-type protease revealed numerous key hydrogen bonding interactions with Thfg. These results suggest that the strategy to design allophenylnorstatine-based peptidomimetics combined with Thfg residue would be promising for generating candidates to overcome multidrug resistance.

Effects of Hinge-region Natural Polymorphisms on Human Immunodeficiency Virus-Type 1 Protease Structure, Dynamics, and Drug Pressure Evolution

Multidrug resistance to current Food and Drug Administration-approved HIV-1 protease (PR) inhibitors drives the need to understand the fundamental mechanisms of how drug pressure-selected mutations, which are oftentimes natural polymorphisms, elicit their effect on enzyme function and resistance. Here, the impacts of the hinge-region natural polymorphism at residue 35, glutamate to aspartate (E35D), alone and in conjunction with residue 57, arginine to lysine (R57K), are characterized with the goal of understanding how altered salt bridge interactions between the hinge and flap regions are associated with changes in structure, motional dynamics, conformational sampling, kinetic parameters, and inhibitor affinity. The combined results reveal that the single E35D substitution leads to diminished salt bridge interactions between residues 35 and 57 and gives rise to the stabilization of open-like conformational states with overall increased backbone dynamics. In HIV-1 PR constructs where sites 35 and 57 are both mutated (e.g. E35D and R57K), x-ray structures reveal an altered network of interactions that replace the salt bridge thus stabilizing the structural integrity between the flap and hinge regions. Despite the altered conformational sampling and dynamics when the salt bridge is disrupted, enzyme kinetic parameters and inhibition constants are similar to those obtained for subtype B PR. Results demonstrate that these hinge-region natural polymorphisms, which may arise as drug pressure secondary mutations, alter protein dynamics and the conformational landscape, which are important thermodynamic parameters to consider for development of inhibitors that target for non-subtype B PR.

Pulsed EPR characterization of HIV-1 protease conformational sampling and inhibitor-induced population shifts

The conformational landscape of HIV-1 protease (PR) can be experimentally characterized by pulsed-EPR double electron-electron resonance (DEER). For this characterization, nitroxide spin labels are attached to an engineered cysteine residue in the flap region of HIV-1 PR. DEER distance measurements from spin-labels contained within each flap of the homodimer provide a detailed description of the conformational sampling of apo-enzyme as well as induced conformational shifts as a function of inhibitor binding. The distance distribution profiles are further interpreted in terms of a conformational ensemble scheme that consists of four unique states termed "curled/tucked", "closed", "semi-open" and "wide-open" conformations. Reported here are the DEER results for a drug-resistant variant clinical isolate sequence, V6, in the presence of FDA approved protease inhibitors (PIs) as well as a non-hydrolyzable substrate mimic, CaP2. Results are interpreted in the context of the current understanding of the relationship between conformational sampling, drug resistance, and kinetic efficiency of HIV-1PR as derived from previous DEER and kinetic data for a series of HIV-1PR constructs that contain drug-pressure selected mutations or natural polymorphisms. Specifically, these collective results support the notion that inhibitor-induced closure of the flaps correlates with inhibitor efficiency and drug resistance. This body of work also suggests DEER as a tool for studying conformational sampling in flexible enzymes as it relates to function.

Spin labeling and Double Electron-Electron Resonance (DEER) to Deconstruct Conformational Ensembles of HIV Protease

An understanding of macromolecular conformational equilibrium in biological systems is oftentimes essential to understand function, dysfunction, and disease. For the past few years, our lab has been utilizing site-directed spin labeling (SDSL), coupled with electron paramagnetic resonance (EPR) spectroscopy, to characterize the conformational ensemble and ligand-induced conformational shifts of HIV-1 protease (HIV-1PR). The biomedical importance of characterizing the fractional occupancy of states within the conformational ensemble critically impacts our hypothesis of a conformational selection mechanism of drug-resistance evolution in HIV-1PR. The purpose of the following chapter is to give a timeline perspective of our SDSL EPR approach to characterizing conformational sampling of HIV-1PR. We provide detailed instructions for the procedure utilized in analyzing distance profiles for HIV-1PR obtained from pulsed electron-electron double resonance (PELDOR). Specifically, we employ a version of PELDOR known as double electron-electron resonance (DEER). Data are processed with the software package "DeerAnalysis" (http://www.epr.ethz.ch/software), which implements Tikhonov regularization (TKR), to generate a distance profile from electron spin-echo amplitude modulations. We assign meaning to resultant distance profiles based upon a conformational sampling model, which is described herein. The TKR distance profiles are reconstructed with a linear combination of Gaussian functions, which is then statistically analyzed. In general, DEER has proven powerful for observing structural ensembles in proteins and, more recently, nucleic acids. Our goal is to present our advances in order to aid readers in similar applications.

The role of select subtype polymorphisms on HIV-1 protease conformational sampling and dynamics

HIV-1 protease is an essential enzyme for viral particle maturation and is a target in the fight against HIV-1 infection worldwide. Several natural polymorphisms are also associated with drug resistance. Here, we utilized both pulsed electron double resonance, also called double electron-electron resonance, and NMR (15)N relaxation measurements to characterize equilibrium conformational sampling and backbone dynamics of an HIV-1 protease construct containing four specific natural polymorphisms commonly found in subtypes A, F, and CRF_01 A/E. Results show enhanced backbone dynamics, particularly in the flap region, and the persistence of a novel conformational ensemble that we hypothesize is an alternative flap orientation of a curled open state or an asymmetric configuration when interacting with inhibitors.

Natural polymorphisms and unusual mutations in HIV-1 protease with potential antiretroviral resistance: a bioinformatic analysis

Background

The correlations of genotypic and phenotypic tests with treatment, clinical history and the significance of mutations in viruses of HIV-infected patients are used to establish resistance mutations to protease inhibitors (PIs). Emerging mutations in human immunodeficiency virus type 1 (HIV-1) protease confer resistance to PIs by inducing structural changes at the ligand interaction site. The aim of this study was to establish an in silico structural relationship between natural HIV-1 polymorphisms and unusual HIV-1 mutations that confer resistance to PIs.

Results

Protease sequences isolated from 151 Mexican HIV-1 patients that were naïve to, or subjected to antiretroviral therapy, were examined. We identified 41 unrelated resistance mutations with a prevalence greater than 1%. Among these mutations, nine exhibited positive selection, three were natural polymorphisms (L63S/V/H) in a codon associated with drug resistance, and six were unusual mutations (L5F, D29V, L63R/G, P79L and T91V). The D29V mutation, with a prevalence of 1.32% in the studied population, was only found in patients treated with antiretroviral drugs. Using in silico modelling, we observed that D29V formed unstable protease complexes when were docked with lopinavir, saquinavir, darunavir, tipranavir, indinavir and atazanavir.

Conclusions

The structural correlation of natural polymorphisms and unusual mutations with drug resistance is useful for the identification of HIV-1 variants with potential resistance to PIs. The D29V mutation likely confers a selection advantage in viruses; however, in silico, presence of this mutation results in unstable enzyme/PI complexes, that possibly induce resistance to PIs.

Consensus HIV-1 subtype A integrase and its raltegravir-resistant variants: design and characterization of the enzymatic properties

Model studies of the subtype B and non-subtype B integrases are still required to compare their susceptibility to antiretroviral drugs, evaluate the significance of resistance mutations and identify the impact of natural polymorphisms on the level of enzymatic reactivity. We have therefore designed the consensus integrase of the HIV-1 subtype A strain circulating in the former Soviet Union territory (FSU-A) and two of its variants with mutations of resistance to the strand transfer inhibitor raltegravir. Their genes were synthesized, and expressed in E coli; corresponding His-tagged proteins were purified using the affinity chromatography. The enzymatic properties of the consensus integrases and their sensitivity to raltegravir were examined in a series of standard in vitro reactions and compared to the properties of the integrase of HIV-1 subtype B strain HXB2. The consensus enzyme demonstrated similar DNA-binding properties, but was significantly more active than HXB-2 integrase in the reactions of DNA cleavage and integration. All integrases were equally susceptible to inhibition by raltegravir and elvitegravir, indicating that the sporadic polymorphisms inherent to the HXB-2 enzyme have little effect on its susceptibility to drugs. Insensitivity of the mutated enzymes to the inhibitors of strand transfer occurred at a cost of a 30-90% loss of the efficacies of both 3'-processing and strand transfer. This is the first study to describe the enzymatic properties of the consensus integrase of HIV-1 clade A and the effects of the resistance mutations when the complex actions of sporadic sequence polymorphisms are excluded.

Inhibitor-induced conformational shifts and ligand-exchange dynamics for HIV-1 protease measured by pulsed EPR and NMR spectroscopy

Double electron-electron resonance (DEER) spectroscopy was utilized to investigate shifts in conformational sampling induced by nine FDA-approved protease inhibitors (PIs) and a nonhydrolyzable substrate mimic for human immunodeficiency virus type 1 protease (HIV-1 PR) subtype B, subtype C, and CRF_01 A/E. The ligand-bound subtype C protease has broader DEER distance profiles, but trends for inhibitor-induced conformational shifts are comparable to those previously reported for subtype B. Ritonavir, one of the strong-binding inhibitors for subtypes B and C, induces less of the closed conformation in CRF_01 A/E. (1)H-(15)N heteronuclear single-quantum coherence (HSQC) spectra were acquired for each protease construct titrated with the same set of inhibitors. NMR (1)H-(15)N HSQC titration data show that inhibitor residence time in the protein binding pocket, inferred from resonance exchange broadening, shifting or splitting correlates with the degree of ligand-induced flap closure measured by DEER spectroscopy. These parallel results show that the ligand-induced conformational shifts resulting from protein-ligand interactions characterized by DEER spectroscopy of HIV-1 PR obtained at the cryogenic temperature are consistent with more physiological solution protein-ligand interactions observed by solution NMR spectroscopy.

Co-circulation HIV-1 subtypes B, C, and CRF31_BC in a drug-naïve population from Southernmost Brazil: analysis of primary resistance mutations

In Southernmost Brazil HIV-1 subtypes B, C, and CRF31_BC co-circulates and, since 1996 with the implementation of free access to highly active antiretroviral treatment (HAART), this epidemic is under a quite characteristic selective pressure. The profile of mutations and polymorphisms in the protease (PR) and reverse transcriptase (RT) genes of HIV-1 from untreated patients living in Porto Alegre, Southernmost Brazil were evaluated in order to identify the subtypes and circulating drug resistant genotypes. Blood samples from 99 HIV-1 positive drugs-naïve patients were collected from 2006 to 2007 in Porto Alegre, Brazil. HIV PR and RT genes were amplified, sequenced, and subtyped. The HIV-1 genotyping was performed by partial sequence analysis of the pol in the HIV Drug Resistance Database of Stanford University. Phylogenetic analyses allowed to classify the HIV samples according to their subtypes: B (26.2%), C (39.4%), F (1.1%), CRF31_CB (19.2%), and URF (14.1%). Eight (8.1%) samples showed primary resistance mutations according to the Calibrated Population Resistance tool based in the 2009 Surveillance Drug Resistance Mutation list. Two samples presented resistance mutations to PI, three NRTI and three NNRTI. There was no significant association between presence of resistant genotypes and subtypes, but resistance mutations seem to be less frequent in the subtype C. In addition, this study describes for the first time the mutational profile of CRF31_BC to PI, NRTI, and NNRTI. Genetic analyses of HIV-1 from naïve patients are a promising and important method for surveillance of HIV infection.

Structural and biochemical characterization of the inhibitor complexes of xenotropic murine leukemia virus-related virus protease

Interactions between the protease (PR) encoded by the xenotropic murine leukemia virus-related virus and a number of potential inhibitors have been investigated by biochemical and structural techniques. It was observed that several inhibitors used clinically against HIV PR exhibit nanomolar or even subnanomolar values of K(i) , depending on the exact experimental conditions. Both TL-3, a universal inhibitor of retroviral PRs, and some inhibitors originally shown to inhibit plasmepsins were also quite potent, whereas inhibition by pepstatin A was considerably weaker. Crystal structures of the complexes of xenotropic murine leukemia virus-related virus PR with TL-3, amprenavir and pepstatin A were solved at high resolution and compared with the structures of complexes of these inhibitors with other retropepsins. Whereas TL-3 and amprenavir bound in a predictable manner, spanning the substrate-binding site of the enzyme, two molecules of pepstatin A bound simultaneously in an unprecedented manner, leaving the catalytic water molecule in place.

Autocatalytic maturation, physical/chemical properties, and crystal structure of group N HIV-1 protease: relevance to drug resistance

The mature protease from Group N human immunodeficiency virus Type 1 (HIV-1) (PR1(N)) differs in 20 amino acids from the extensively studied Group M protease (PR1(M)) at positions corresponding to minor drug-resistance mutations (DRMs). The first crystal structure (1.09 Å resolution) of PR1(N) with the clinical inhibitor darunavir (DRV) reveals the same overall structure as PR1(M), but with a slightly larger inhibitor-binding cavity. Changes in the 10s loop and the flap hinge propagate to shift one flap away from the inhibitor, whereas L89F and substitutions in the 60s loop perturb inhibitor-binding residues 29-32. However, kinetic parameters of PR1(N) closely resemble those of PR1(M), and calorimetric results are consistent with similar binding affinities for DRV and two other clinical PIs, suggesting that minor DRMs coevolve to compensate for the detrimental effects of drug-specific major DRMs. A miniprecursor (TFR 1-61-PR1(N)) comprising the transframe region (TFR) fused to the N-terminus of PR1(N) undergoes autocatalytic cleavage at the TFR/PR1(N) site concomitant with the appearance of catalytic activity characteristic of the dimeric, mature enzyme. This cleavage is inhibited at an equimolar ratio of precursor to DRV (∼6 μM), which partially stabilizes the precursor dimer from a monomer. However, cleavage at L34/W35 within the TFR, which precedes the TFR 1-61/PR1(N) cleavage at pH ≤ 5, is only partially inhibited. Favorable properties of PR1(N) relative to PR1(M) include its suitability for column fractionation by size under native conditions and >10-fold higher dimer dissociation constant (150 nM). Exploiting these properties may facilitate testing of potential dimerization inhibitors that perturb early precursor processing steps.

Impact of HIV type 1 subtype on drug resistance mutations in Nigerian patients failing first-line therapy

A diverse array of non-subtype B HIV-1 viruses circulates in Africa and dominates the global pandemic. It is important to understand how drug resistance mutations in non-B subtypes may develop differently from the patterns described in subtype B. HIV-1 reverse transcriptase and protease sequences from 338 patients with treatment failure to first-line ART regimens were evaluated. Multivariate logistic regression was used to examine the effect of subtype on each mutation controlling for regimen, time on therapy, and total mutations. The distribution of HIV-1 subtypes included CRF02_AG (45.0%), G (37.9%), CRF06_cpx (4.4%), A (3.6%), and other subtypes or recombinant sequences (9.2%). The most common NRTI mutations were M184V (89.1%) and thymidine analog mutations (TAMs). The most common NNRTI mutations were Y181C (49.7%), K103N (36.4%), G190A (26.3%), and A98G (19.5%). Multivariate analysis showed that CRF02_AG was less likely to have the M41L mutation compared to other subtypes [adjusted odds ratio (AOR) = 0.35; p = 0.022]. Subtype A patients showed a 42.5-fold increased risk (AOR = 42.5, p = 0.001) for the L210W mutation. Among NNRTI mutations, subtype G patients had an increased risk for A98G (AOR = 2.40, p = 0.036) and V106I (AOR = 6.15, p = 0.010), whereas subtype CRF02_AG patients had an increased risk for V90I (AOR = 3.16; p = 0.003) and a decreased risk for A98G (AOR = 0.48, p = 0.019). Five RT mutations were found to vary significantly between different non-B West African subtypes. Further study to understand the clinical impact of subtype-specific diversity on drug resistance will be critically important to the continued success of ART scale-up in resource-limited settings.

Molecular Basis for Drug Resistance in HIV-1 Protease

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

Short communication: antiretroviral therapy resistance mutations present in the HIV type 1 subtype C pol and env regions from therapy-naive patients in Zambia

The prevalence of antiretroviral therapy (ART) resistance mutations present in HIV-1 subtype C pol and env regions of the proviral DNA was analyzed and compared from therapy-naive individuals before (Cohort A) and after (Cohort B) the availability of free ART in Zambia. Mutations present in sequences published in a previous study from Zambian ART-naive individuals infected with subtype C were analyzed using current parameters for the classification of ART drug resistance and compared with Cohorts A and B. No statistically significant differences were observed when comparing mutations present in the pol and env of these cohorts. However, an increase in the number of minor, borderline, or partial resistance mutations as well as the presence of major resistance mutations were observed in Cohort B. These results suggest there is an increasing trend of drug resistance-associated mutations that could be a result of the availability of free ART in Zambia. Moreover, the high prevalence of resistance mutations observed for maraviroc and vicriviroc in both cohorts may suggest a limited efficacy of entry inhibitors on HIV-1 subtype C.

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

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

Structure of the unbound form of HIV-1 subtype A protease: comparison with unbound forms of proteases from other HIV subtypes

The crystal structure of the unbound form of HIV-1 subtype A protease (PR) has been determined to 1.7 A resolution and refined as a homodimer in the hexagonal space group P6(1) to an R(cryst) of 20.5%. The structure is similar in overall shape and fold to the previously determined subtype B, C and F PRs. The major differences lie in the conformation of the flap region. The flaps in the crystal structures of the unbound subtype B and C PRs, which were crystallized in tetragonal space groups, are either semi-open or wide open. In the present structure of subtype A PR the flaps are found in the closed position, a conformation that would be more anticipated in the structure of HIV protease complexed with an inhibitor. The amino-acid differences between the subtypes and their respective crystal space groups are discussed in terms of the differences in the flap conformations.

Lead expansion and virtual screening of Indinavir derivate HIV-1 protease inhibitors using pharmacophoric - shape similarity scoring function

Indinavir (Crivaxan®) is a potent inhibitor of the HIV (human immunodeficiency virus) protease. This enzyme has an important role in viral replication and is considered to be very attractive target for new antiretroviral drugs. However, it becomes less effective due to highly resistant new viral strains of HIV, which have multiple mutations in their proteases. For this reason, we used a lead expansion method to create a new set of compounds with a new mode of action to protease binding site. 1300 compounds chemically diverse from the initial hit were generated and screened to determine their ability to interact with protease and establish their QSAR properties. Further computational analyses revealed one unique compound with different protease binding ability from the initial hit and its role for possible new class of protease inhibitors is discussed in this report.

Subtype polymorphisms among HIV-1 protease variants confer altered flap conformations and flexibility

Human immunodeficiency virus type 1 (HIV-1) protease plays a fundamental role in the maturation and life cycle of the retrovirus HIV-1, as it functions in regulating post-translational processing of the viral polyproteins gag and gag-pol; thus, it is a key target of AIDS antiviral therapy. Accessibility of substrate to the active site is mediated by two flaps, which must undergo a large conformational change from an open to a closed conformation during substrate binding and catalysis. The electron paramagnetic resonance (EPR) method of site-directed spin labeling (SDSL) with double electron-electron resonance (DEER) spectroscopy was utilized to monitor the conformations of the flaps in apo HIV-1 protease (HIV-1PR), subtypes B, C, and F, CRF01_A/E, and patient isolates V6 and MDR 769. The distance distribution profiles obtained from analysis of the dipolar modulated echo curves were reconstructed to yield a set of Gaussian-shaped populations, which provide an analysis of the flap conformations sampled. The relative percentages of each conformer population described as "tucked/curled", "closed", "semi-open", and "wide-open" were determined and compared for various constructs. The results and analyses show that sequence variations among subtypes, CRFs, and patient isolates of apo HIV-1PR alter the average flap conformation in a way that can be understood as inducing shifts in the relative populations, or conformational sampling, of the previously described four conformations for HIV-1PR.

Genetic determinants in HIV-1 Gag and Env V3 are related to viral response to combination antiretroviral therapy with a protease inhibitor

OBJECTIVE

To identify novel viral determinants in HIV-1 protease, Gag, and envelope V3 that relate to outcomes to initial protease inhibitor-based antiretroviral therapy.

DESIGN

A longitudinal cohort study of protease inhibitor-naive, HIV-infected individuals was designed to identify genetic variables in viral Gag and envelope sequences associated with response to antiretroviral therapy.

METHODS

Genetic and statistical models, including amino acid profiles, phylogenetic analyses, receiver operating characteristic analyses, and covariation analyses, were used to evaluate viral sequences and clinical variables from individuals who developed immune reconstitution with or without suppression of viral replication.

RESULTS

Pretherapy chemokine (C-X-C motif) receptor 4-using V3 regions had significant associations with viral failure (P = 0.04). Amino acid residues in protease covaried with Gag residues, particularly in p7(NC), independent of cleavage sites. Pretherapy V3 charge combined with p6(Pol) and p2/p7(NC) cleavage site genotypes produced the best three-variable model to predict viral suppression in 88% of individuals. Combinations of baseline CD4 cell percentage with genetic determinants in Gag-protease predicted viral fitness in 100% of individuals who failed to suppress viral replication.

CONCLUSION

Baseline genetic determinants in Gag p6(Pol) and p2/p7(NC), as well as envelope, provide novel combinations of biomarkers for predicting emergence of viral resistance to initial therapy regimens.

Analysis of binding parameters of HIV-1 integrase inhibitors: correlates of drug inhibition and resistance

This study undertook an exploratory data analysis of the binding parameters of HIV-1 integrase inhibitors. The study group involved inhibitors in preclinical development from the diketo acid, pyrroloquinoline and naphthyridine carboxamide families and the most advanced inhibitors Raltegravir and Elvitegravir. Distinct differences were observed in the energetics of binding between the studied classes of inhibitors that also correlated with drug resistant patterns. Quantitative-property-activity-relationships correlated experimental IC(50) values to the binding energy and the logarithm of the partition coefficient between n-octanol and water (clogP). The approach followed here serves as an improved basis for the development of 'second generation' integrase inhibitors.

Implications of HIV-1 M group polymorphisms on integrase inhibitor efficacy and resistance: genetic and structural in silico analyses

The extensive polymorphisms among HIV-1 subtypes have been implicated in drug resistance development. Integrase inhibitors represent the latest addition to the treatment of HIV-1, and their efficacy and resistance patterns among M group strains are currently under investigation. This study analyzed the intersubtype variation within 108 integrase sequences from seven subtypes. The residues associated with catalytic activity and primary resistance to raltegravir were highly conserved among all strains. Variations were observed in residues associated with secondary resistance. Molecular modeling studies indicated a two-way binding mode of raltegravir that explains the resistance pathways and the implication of nonconservative mutations in integrase-raltegravir interactions.

Active-site mutations in the South african human immunodeficiency virus type 1 subtype C protease have a significant impact on clinical inhibitor binding: kinetic and thermodynamic study

Human immunodeficiency virus (HIV) infections in sub-Saharan Africa represent about 56% of global infections. Study of active-site mutations (the V82A single mutation and the V82F I84V double mutation) in the less-studied South African HIV type 1 subtype C (C-SA) protease indicated that neither mutation had a significant impact on the proteolytic functioning of the protease. However, the binding affinities of, and inhibition by, saquinavir, ritonavir, indinavir, and nelfinavir were weaker for each variant than for the wild-type protease, with the double mutant exhibiting the most dramatic change. Therefore, our results show that the C-SA V82F I84V double mutation decreased the binding affinities of protease inhibitors to levels significantly lower than that required for effective inhibition.

Effect of flap mutations on structure of HIV-1 protease and inhibition by saquinavir and darunavir

HIV-1 (human immunodeficiency virus type 1) protease (PR) and its mutants are important antiviral drug targets. The PR flap region is critical for binding substrates or inhibitors and catalytic activity. Hence, mutations of flap residues frequently contribute to reduced susceptibility to PR inhibitors in drug-resistant HIV. Structural and kinetic analyses were used to investigate the role of flap residues Gly48, Ile50, and Ile54 in the development of drug resistance. The crystal structures of flap mutants PR(I50V) (PR with I50V mutation), PR(I54V) (PR with I54V mutation), and PR(I54M) (PR with I54M mutation) complexed with saquinavir (SQV) as well as PR(G48V) (PR with G48V mutation), PR(I54V), and PR(I54M) complexed with darunavir (DRV) were determined at resolutions of 1.05-1.40 A. The PR mutants showed changes in flap conformation, interactions with adjacent residues, inhibitor binding, and the conformation of the 80s loop relative to the wild-type PR. The PR contacts with DRV were closer in PR(G48V)-DRV than in the wild-type PR-DRV, whereas they were longer in PR(I54M)-DRV. The relative inhibition of PR(I54V) and that of PR(I54M) were similar for SQV and DRV. PR(G48V) was about twofold less susceptible to SQV than to DRV, whereas the opposite was observed for PR(I50V). The observed inhibition was in agreement with the association of G48V and I50V with clinical resistance to SQV and DRV, respectively. This analysis of structural and kinetic effects of the mutants will assist in the development of more effective inhibitors for drug-resistant HIV.