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.

Modeling and Analysis of HIV-1 Pol Polyprotein as a Case Study for Predicting Large Polyprotein Structures.

Acquired immunodeficiency syndrome (AIDS) is caused by human immunodeficiency virus (HIV). HIV protease, reverse transcriptase, and integrase are targets of current drugs to treat the disease. However, anti-viral drug-resistant strains have emerged quickly due to the high mutation rate of the virus, leading to the demand for the development of new drugs. One attractive target is Gag-Pol polyprotein, which plays a key role in the life cycle of HIV. Recently, we found that a combination of M50I and V151I mutations in HIV-1 integrase can suppress virus release and inhibit the initiation of Gag-Pol autoprocessing and maturation without interfering with the dimerization of Gag-Pol. Additional mutations in integrase or RNase H domain in reverse transcriptase can compensate for the defect. However, the molecular mechanism is unknown. There is no tertiary structure of the full-length HIV-1 Pol protein available for further study. Therefore, we developed a workflow to predict the tertiary structure of HIV-1 NL4.3 Pol polyprotein. The modeled structure has comparable quality compared with the recently published partial HIV-1 Pol structure (PDB ID: 7SJX). Our HIV-1 NL4.3 Pol dimer model is the first full-length Pol tertiary structure. It can provide a structural platform for studying the autoprocessing mechanism of HIV-1 Pol and for developing new potent drugs. Moreover, the workflow can be used to predict other large protein structures that cannot be resolved via conventional experimental methods.

Preclinical characterization of a non-peptidomimetic HIV protease inhibitor with improved metabolic stability.

Protease inhibitors (PIs) remain an important component of antiretroviral therapy for the treatment of HIV-1 infection due to their high genetic barrier to resistance development. Nevertheless, the two most commonly prescribed HIV PIs, atazanavir and darunavir, still require co-administration with a pharmacokinetic boosting agent to maintain sufficient drug plasma levels which can lead to undesirable drug-drug interactions. Herein, we describe GS-9770, a novel investigational non-peptidomimetic HIV PI with unboosted once-daily oral dosing potential due to improvements in its metabolic stability and its pharmacokinetic properties in preclinical animal species. This compound demonstrates potent inhibitory activity and high on-target selectivity for recombinant HIV-1 protease versus other aspartic proteases tested. In cell culture, GS-9770 inhibits Gag polyprotein cleavage and shows nanomolar anti-HIV-1 potency in primary human cells permissive to HIV-1 infection and against a broad range of HIV subtypes. GS-9770 demonstrates an improved resistance profile against a panel of patient-derived HIV-1 isolates with resistance to atazanavir and darunavir. In resistance selection experiments, GS-9770 prevented the emergence of breakthrough HIV-1 variants at all fixed drug concentrations tested and required multiple protease substitutions to enable outgrowth of virus exposed to escalating concentrations of GS-9770. This compound also remained fully active against viruses resistant to drugs from other antiviral classes and showed no in vitro antagonism when combined pairwise with drugs from other antiretroviral classes. Collectively, these preclinical data identify GS-9770 as a potent, non-peptidomimetic once-daily oral HIV PI with potential to overcome the persistent requirement for pharmacological boosting with this class of antiretroviral agents.

Cellular Targets of HIV-1 Protease: Just the Tip of the Iceberg?

Human immunodeficiency virus 1 (HIV-1) viral protease (PR) is one of the most studied viral enzymes and a crucial antiviral target. Despite its well-characterized role in virion maturation, an increasing body of research is starting to focus on its ability to cleave host cell proteins. Such findings are apparently in contrast with the dogma of HIV-1 PR activity being restricted to the interior of nascent virions and suggest catalytic activity within the host cell environment. Given the limited amount of PR present in the virion at the time of infection, such events mainly occur during late viral gene expression, mediated by newly synthesized Gag-Pol polyprotein precursors, rather than before proviral integration. HIV-1 PR mainly targets proteins involved in three different processes: those involved in translation, those controlling cell survival, and restriction factors responsible for innate/intrinsic antiviral responses. Indeed, by cleaving host cell translation initiation factors, HIV-1 PR can impair cap-dependent translation, thus promoting IRES-mediated translation of late viral transcripts and viral production. By targeting several apoptotic factors, it modulates cell survival, thus promoting immune evasion and viral dissemination. Additionally, HIV-1 PR counteracts restriction factors incorporated in the virion that would otherwise interfere with nascent virus vitality. Thus, HIV-1 PR appears to modulate host cell function at different times and locations during its life cycle, thereby ensuring efficient viral persistency and propagation. However, we are far from having a complete picture of PR-mediated host cell modulation, which is emerging as a field that needs further investigation.

Different Mutation Tolerance of Lentiviral (HIV-1) and Deltaretroviral (BLV and HTLV) Protease Precursors.

The bovine leukemia virus (BLV) and the human T-lymphothropic viruses (HTLVs) are members of the deltaretrovirus genus of Retroviridae family. An essential event of the retroviral life cycle is the processing of the polyproteins by the viral protease (PR); consequently, these enzymes became important therapeutic targets of the anti-retroviral drugs. As compared to human immunodeficiency viruses (HIVs), the deltaretroviruses have a different replication strategy, as they replicate predominantly in the DNA form, by forcing the infected cell to divide, unlike HIV-1, which replicates mainly by producing a vast number of progeny virions and by reinfection. Due to bypassing the error-prone reverse transcription step of replication, the PRs of deltaretroviruses did not undergo such extensive evolution as HIV PRs and remained more highly conserved. In this work, we studied the abilities of wild-type and modified BLV, HTLV (type 1, 2 and 3), and HIV-1 PRs (fused to an N-terminal MBP tag) for self-processing. We designed a cleavage site mutant MBP-fused BLV PR precursor as well, this recombinant enzyme was unable for self-proteolysis, the MBP fusion tag decreased its catalytic efficiency but showed an unusually low Ki for the IB-268 protease inhibitor. Our results show that the HTLV and BLV deltaretrovirus PRs exhibit lower mutation tolerance as compared to HIV-1 PR, and are less likely to retain their activity upon point mutations at various positions, indicating a higher flexibility of HIV-1 PR in tolerating mutations under selective pressure.

Protease and gag diversity and drug resistance mutations among treatment-naive Mexican people living with HIV.

INTRODUCTION: In Mexico, HIV genotyping is performed in people living with HIV (PLWH) failing their first-line antiretroviral (ARV) regimen; it is not routinely done for all treatment-naive PLWH before ARV initiation. The first nationally representative survey published in 2016 reported that the prevalence of pretreatment drug mutations in treatment-naive Mexican PLWH was 15.5% to any antiretroviral drug and 10.6% to non-nucleoside reverse transcriptase inhibitors (NNRTIs) using conventional Sanger sequencing. Most reports in Mexico focus on HIV pol gene and nucleoside and non-nucleoside reverse transcriptase inhibitor (NRTI and NNRTI) drug resistance mutations (DRMs) prevalence, using Sanger sequencing, next-generation sequencing (NGS) or both. To our knowledge, NGS has not be used to detect pretreatment drug resistance mutations (DRMs) in the HIV protease (PR) gene and its substrate the Gag polyprotein. METHODS: Treatment-naive adult Mexican PLWH were recruited between 2016 and 2019. HIV Gag and protease sequences were obtained by NGS and DRMs were identified using the WHO surveillance drug resistance mutation (SDRM) list. RESULTS: One hundred PLWH attending a public national reference hospital were included. The median age was 28 years-old, and most were male. The median HIV viral load was 4.99 [4.39-5.40] log copies/mL and median CD4 cell count was 150 [68.0-355.78] cells/mm3. As expected, most sequences clustered with HIV-1 subtype B (97.9%). Major PI resistance mutations were detected: 8 (8.3%) of 96 patients at a detection threshold of 1% and 3 (3.1%) at a detection threshold of 20%. A total of 1184 mutations in Gag were detected, of which 51 have been associated with resistance to PI, most of them were detected at a threshold of 20%. Follow-up clinical data was available for 79 PLWH at 6 months post-ART initiation, seven PLWH failed their first ART regimen; however no major PI mutations were identified in these individuals at baseline. CONCLUSIONS: The frequency of DRM in the HIV protease was 7.3% at a detection threshold of 1% and 3.1% at a detection threshold of 20%. NGS-based HIV drug resistance genotyping provide improved detection of DRMs. Viral load was used to monitor ARV response and treatment failure was 8.9%.

Co-evolution of drug resistance and broadened substrate recognition in HIV protease variants isolated from an Escherichia coli genetic selection system.

A genetic selection system for activity of HIV protease is described that is based on a synthetic substrate constructed as a modified AraC regulatory protein that when cleaved stimulate l-arabinose metabolism in an Escherichia coli araC strain. Growth stimulation on selective plates was shown to depend on active HIV protease and the scissile bond in the substrate. In addition, the growth of cells correlated well with the established cleavage efficiency of the sites in the viral polyprotein, Gag, when these sites were individually introduced into the synthetic substrate of the selection system. Plasmids encoding protease variants selected based on stimulation of cell growth in the presence of saquinavir or cleavage of a site not cleaved by wild-type protease, were indistinguishable with respect to both phenotypes. Also, both groups of selected plasmids encoded side chain substitutions known from clinical isolates or displayed different side chain substitutions but at identical positions. One highly frequent side chain substitution, E34V, not regarded as a major drug resistance substitution was found in variants obtained under both selective conditions and is suggested to improve protease processing of the synthetic substrate. This substitution is away from the substrate-binding cavity and together with other substitutions in the selected reading frames supports the previous suggestion of a substrate-binding site extended from the active site binding pocket itself.

Amino acid substitutions at the HIV-1 transframe region significantly impair virus infectivity.

A transframe region within HIV-1 Gag-Pol (referred to as p6* or p6pol), directly linked to the protease (PR) N-terminus, plays a pivotal role in modulating PR activation. To identify specific p6* residues involved in PR activation, we created a series of p6* mutants by making substitutions for conserved p6* residues. Our results indicate that some p6* mutants were defective in terms of virus infectivity, despite displaying a wild-type virus particle processing pattern. Mutations at p6* F8 reduced virus infectivity associated with insufficient virus processing, due in part to impaired PR maturation and RT packaging. Our data strongly suggest that conserved Phe (F) residues at position 8 of p6* are involved in the PR maturation process.

Multiple Molecular Dynamics Simulations and Free-Energy Predictions Uncover the Susceptibility of Variants of HIV-1 Protease against Inhibitors Darunavir and KNI-1657.

HIV-1 protease (PR) is considered to be the main targets of anti-AIDS drug design because of its role in the proteolytic processing of viral polyproteins. However, the emergence of drug-resistant HIV has become a major problem in the therapy of HIV-1-infected patients. Focused on the complexes of wild type (WT) PR and two mutant PRs (V32I/L33F/I54M/V82I and V32I/L33F/I54M/I84 V) with inhibitors Darunavir (DRV) and KNI-1657 (KNI), respectively, we have conducted research on the conformational dynamics and the resistance mechanism caused by residue mutations through multiple molecular dynamics (MD) simulations combined with an energy (MM-PBSA and solvated interaction energy (SIE)) prediction. The results indicate that mutated residues of PR alter the distance between flap regions and catalytic sites, the volume of the inner catalytic site, and the curling degree of the flap tips, thereby affecting DRV and KNI inhibitor binding to PR. These mutated residues reduced the binding affinity of the two mutant PRs to DRV, resulting in drug resistance, whereas the two mutant PRs increase the binding affinity with KNI, indicating they enhance the sensitivity to KNI. Compared with the WT PR, the changes in van der Waals interaction and electrostatic interaction in the two variant PRs play a vital part in the binding of PR with DRV and KNI. These results may supply valuable guidance for the design of anti-AIDS drugs targeting PR.

Darunavir-Resistant HIV-1 Protease Constructs Uphold a Conformational Selection Hypothesis for Drug Resistance.

Multidrug resistance continues to be a barrier to the effectiveness of highly active antiretroviral therapy in the treatment of human immunodeficiency virus 1 (HIV-1) infection. Darunavir (DRV) is a highly potent protease inhibitor (PI) that is oftentimes effective when drug resistance has emerged against first-generation inhibitors. Resistance to darunavir does evolve and requires 10-20 amino acid substitutions. The conformational landscapes of six highly characterized HIV-1 protease (PR) constructs that harbor up to 19 DRV-associated mutations were characterized by distance measurements with pulsed electron double resonance (PELDOR) paramagnetic resonance spectroscopy, namely double electron-electron resonance (DEER). The results show that the accumulated substitutions alter the conformational landscape compared to PI-naïve protease where the semi-open conformation is destabilized as the dominant population with open-like states becoming prevalent in many cases. A linear correlation is found between values of the DRV inhibition parameter Ki and the open-like to closed-state population ratio determined from DEER. The nearly 50% decrease in occupancy of the semi-open conformation is associated with reduced enzymatic activity, characterized previously in the literature.

Structure-Based Design of Highly Potent HIV-1 Protease Inhibitors Containing New Tricyclic Ring P2-Ligands: Design, Synthesis, Biological, and X-ray Structural Studies.

We describe here design, synthesis, and biological evaluation of a series of highly potent HIV-1 protease inhibitors containing stereochemically defined and unprecedented tricyclic furanofuran derivatives as P2 ligands in combination with a variety of sulfonamide derivatives as P2' ligands. These inhibitors were designed to enhance the ligand-backbone binding and van der Waals interactions in the protease active site. A number of inhibitors containing the new P2 ligand, an aminobenzothiazole as the P2' ligand and a difluorophenylmethyl as the P1 ligand, displayed very potent enzyme inhibitory potency and also showed excellent antiviral activity against a panel of highly multidrug-resistant HIV-1 variants. The tricyclic P2 ligand has been synthesized efficiently in an optically active form using enzymatic desymmetrization of meso-1,2-(dihydroxymethyl)cyclohex-4-ene as the key step. We determined high-resolution X-ray structures of inhibitor-bound HIV-1 protease. These structures revealed extensive interactions with the backbone atoms of HIV-1 protease and provided molecular insights into the binding properties of these new inhibitors.

HIV-1 Mutant Assembly, Processing and Infectivity Expresses Pol Independent of Gag.

The pol retrovirus gene encodes required enzymes for virus replication and maturation. Unlike HIV-1 Pol (expressed as a Gag-Pol fusion protein), foamy virus (described as an ancient retrovirus) expresses Pol without forming Gag-Pol polyproteins. We placed a "self-cleaving" 2A peptide between HIV-1 Gag and Pol. This construct, designated G2AP, is capable of producing virions with the same density as a wild-type (wt) HIV-1 particle. The 2A peptide allows for Pol to be packaged into virions independently from Gag following co-translationally cleaved from Gag. We found that G2AP exhibited only one-third the virus infectivity of the wt, likely due, at least in part, to defects in Pol packaging. Attenuated protease (PR) activity, or a reduction in Pol expression due to the placement of 2A-mediated Pol in a normal Gag-Pol frameshift context, resulted in significant increases in virus yields and/or titers. This suggests that reduced G2AP virus yields were largely due to increased PR activity associated with overexpressed Pol. Our data suggest that HIV-1 adopts a gag/pol ribosomal frameshifting mechanism to support virus assembly via the efficient modulation of Gag-Pol/Gag expression, as well as to promote viral enzyme packaging. Our results help clarify the molecular basis of HIV-1 gene expression and assembly.

Mexican HIV-1 Protease Sequence Diversity.

Protease is one of three enzymes encoded within HIV's pol gene, responsible for the cleavage of viral Gag-Pol polypeptide into mature viral proteins and a target of current anti-retroviral therapy. Protease diversity analysis in Latin America has been lacking in spite of extensive studies of protease-inhibitor resistance mutations. We studied the diversity of 777 Mexican protease sequences and found that all were subtype B except one (CRF02_AG). Phylogenetic analysis suggested the existence of six different clades with geospecific contributions. Thirty-three percent of sites were conserved, 25% had conservative substitutions, and 41% exhibited physicochemical changes. The most conserved regions surrounded the active site, most of the flap domain, and a region between the 60's loop and C-terminal triad. A single sequence exhibited an active site mutation (T26S). Variable sites were mapped to a crystallographic structure, providing further insight into the distribution and functional relevance of variable sites among Mexican isolates.

Energetics of Flap Opening in HIV-1 Protease: String Method Calculations.

HIV-1 protease (PR) is the viral protein responsible for virion maturation, and its mechanisms of action remain incompletely understood. PR is dimeric and contains two flexible, symmetry-related flaps, which act as a gate to inhibit access to the binding pocket and hold the polypeptide substrate in the binding pocket once bound. Wide flap opening, a conformational change assumed to be necessary for substrate binding, is a rare event in the closed and bound form. In this study, we use molecular dynamics (MD) simulations and advanced MD techniques including temperature acceleration and string method in collective variables to study the conformational changes associated with substrate unbinding of both wild-type and F99Y mutant PR. The F99Y mutation is shown via MD to decouple the closing of previously unrecognized distal pockets from substrate unbinding. To determine whether or not the F99Y mutation affects the energetic cost of wide flap opening, we use string method in collective variables to determine the minimum free-energy mechanism for wide flap opening in concert with distal pocket closing. The results indicate that the major energetic cost in flap opening is disengagement of the two flap-tip Ile50 residues from each other and is not affected by the F99Y mutation.

Picomolar to Micromolar: Elucidating the Role of Distal Mutations in HIV-1 Protease in Conferring Drug Resistance.

Drug resistance continues to be a growing global problem. The efficacy of small molecule inhibitors is threatened by pools of genetic diversity in all systems, including antibacterials, antifungals, cancer therapeutics, and antivirals. Resistant variants often include combinations of active site mutations and distal "secondary" mutations, which are thought to compensate for losses in enzymatic activity. HIV-1 protease is the ideal model system to investigate these combinations and underlying molecular mechanisms of resistance. Darunavir (DRV) binds wild-type (WT) HIV-1 protease with a potency of <5 pM, but we have identified a protease variant that loses potency to DRV 150 000-fold, with 11 mutations in and outside the active site. To elucidate the roles of these mutations in DRV resistance, we used a multidisciplinary approach, combining enzymatic assays, crystallography, and molecular dynamics simulations. Analysis of protease variants with 1, 2, 4, 8, 9, 10, and 11 mutations showed that the primary active site mutations caused ∼50-fold loss in potency (2 mutations), while distal mutations outside the active site further decreased DRV potency from 13 nM (8 mutations) to 0.76 µM (11 mutations). Crystal structures and simulations revealed that distal mutations induce subtle changes that are dynamically propagated through the protease. Our results reveal that changes remote from the active site directly and dramatically impact the potency of the inhibitor. Moreover, we find interdependent effects of mutations in conferring high levels of resistance. These mechanisms of resistance are likely applicable to many other quickly evolving drug targets, and the insights may have implications for the design of more robust inhibitors.

PERSIA for Direct Fluorescence Measurements of Transcription, Translation, and Enzyme Activity in Cell-Free Systems.

Quantification of biology's central dogma (transcription and translation) is pursued by a variety of methods. Direct, immediate, and ongoing quantification of these events is difficult to achieve. Common practice is to use fluorescent or luminescent proteins to report indirectly on prior cellular events, such as turning on a gene in a genetic circuit. We present an alternative approach, PURExpress-ReAsH-Spinach In-vitro Analysis (PERSIA). PERSIA provides information on the production of RNA and protein during cell-free reactions by employing short RNA and peptide tags. Upon synthesis, these tags yield quantifiable fluorescent signal without interfering with other biochemical events. We demonstrate the applicability of PERSIA in measuring cell-free transcription, translation, and other enzymatic activity in a variety of applications: from sequence-structure-function studies, to genetic code engineering, to testing antiviral drug resistance.

New heteroaryl carbamates: Synthesis and biological screening in vitro and in mammalian cells of wild-type and mutant HIV-protease inhibitors.

New heteroaryl HIV-protease inhibitors bearing a carbamoyl spacer were synthesized in few steps and high yield, from commercially available homochiral epoxides. Different substitution patterns were introduced onto a given isopropanoyl-sulfonamide core that can have either H or benzyl group. The in vitro inhibition activity against recombinant protease showed a general beneficial effect of both carbamoyl moiety and the benzyl group, ranging the IC50 values between 11 and 0.6 nM. In particular, benzofuryl and indolyl derivatives showed IC50 values among the best for such structurally simple inhibitors. Docking analysis allowed to identify the favorable situation of such derivatives in terms of number of interactions in the active site, supporting the experimental results. The inhibition activity was also confirmed in HEK293 mammalian cells and was maintained against protease mutants. Furthermore, the metabolic stability was comparable with that of the commercially available inhibitors.

HIV transmission networks among transgender women in Los Angeles County, CA, USA: a phylogenetic analysis of surveillance data.

BACKGROUND: Transgender women are among the groups at highest risk for HIV infection, with a prevalence of 27·7% in the USA; and despite this known high risk, undiagnosed infection is common in this population. We set out to identify transgender women and their partners in a molecular transmission network to prioritise public health activities. METHODS: Since 2006, HIV protease and reverse transcriptase gene (pol) sequences from drug resistance testing have been reported to the Los Angeles County Department of Public Health and linked to demographic data, gender, and HIV transmission risk factor data for each case in the enhanced HIV/AIDS Reporting System. We reconstructed a molecular transmission network by use of HIV-TRAnsmission Cluster Engine (with a pairwise genetic distance threshold of 0·015 substitutions per site) from the earliest pol sequences from 22 398 unique individuals, including 412 (2%) self-identified transgender women. We examined the possible predictors of clustering with multivariate logistic regression. We characterised the genetically linked partners of transgender women and calculated assortativity (the tendency for people to link to other people with the same attributes) for each transmission risk group. FINDINGS: 8133 (36·3%) of 22 398 individuals clustered in the network across 1722 molecular transmission clusters. Transgender women who indicated a sexual risk factor clustered at the highest frequency in the network, with 147 (43%) of 345 being linked to at least one other person (adjusted odds ratio [aOR] 2·0, p=0·0002). Transgender women were assortative in the network (assortativity 0·06, p<0·001), indicating that they tended to link to other transgender women. Transgender women were more likely than expected to link to other transgender women (OR 4·65, p<0·001) and cisgender men who did not identify as men who have sex with men (MSM; OR 1·53, p<0·001). Transgender women were less likely than expected to link to MSM (OR 0·75, p<0·001), despite the high prevalence of HIV among MSM. Transgender women were distributed across 126 clusters, and cisgender individuals linked to one transgender woman were 9·2 times more likely to link to a second transgender woman than other individuals in the surveillance database. Reconstruction of the transmission network is limited by sample availability, but sequences were available for more than 40% of diagnoses. INTERPRETATION: Clustering of transgender women and the observed tendency for linkage with cisgender men who did not identify as MSM, shows the potential to use molecular epidemiology both to identify clusters that are likely to include undiagnosed transgender women with HIV and to improve the targeting of public health prevention and treatment services to transgender women. FUNDING: California HIV and AIDS Research Program and National Institutes of Health-National Institute of Allergy and Infectious Diseases.

The effect of primary drug resistance on CD4+ cell decline and the viral load set-point in HIV-positive individuals before the start of antiretroviral therapy.

OBJECTIVE: To evaluate the effect of primary resistance and selected polymorphic amino-acid substitutions in HIV reverse transcriptase and protease on the CD4 cell count and viral load set point before the start of antiretroviral treatment. DESIGN: Prospective cohort study. METHODS: A total of 6180 individuals with a resistance test prior to starting antiretroviral treatment accessing care in HIV clinics across Europe who had at least one viral load and one CD4+ test available were included in the analysis. The impact of amino-acid substitutions variants on viral load and CD4+ trends was investigated using linear mixed models. Clusters of mutations were studied using principal component analysis. RESULTS: Overall, the detection of any primary resistance was not associated with either the speed of CD4+ cell decline or the viral load set point. However, transmitted nucleoside reverse transcriptase inhibitor and protease inhibitor resistance appeared to be weakly associated with lower viral load set points, as were the polymorphic G16E or Q92K protease mutations. There was some evidence suggesting that these effects varied according to HIV subtype, with the effects of transmitted nucleoside reverse transcriptase inhibitor and protease resistance being particularly marked among individuals with a subtype B virus. A cluster of five polymorphic protease substitutions at position 20, 13, 36, 69 and 89 was associated with less steep CD4+ cell declines and lower viral load set points. CONCLUSION: Although we found little evidence for an association between primary resistance and CD4+ speed of decline and viral load set point, the potential role of polymorphic protease (alone or in clusters) and their interplay with HIV subtype needs to be further evaluated.