Prediction Of Coreceptor Usage Research Articles

Adapting the geno2pheno[coreceptor] tool to HIV-1 subtype CRF01_AE by phenotypic validation using clinical isolates from South-East Asia.

Geno2pheno[coreceptor] is a widely used tool for the prediction of coreceptor usage (viral tropism) of HIV-1 samples. For HIV-1 CRF01_AE, a significant overcalling of X4-tropism is observed when using the standard settings of Geno2pheno[coreceptor]. The aim of this study was to provide the experimental backing for adaptations to the geno2pheno[coreceptor] algorithm in order to improve coreceptor usage predictions of clinical HIV-1 CRF01_AE isolates STUDY DESIGN: V3-sequences of 20 clinical HIV-1 subtype CRF01_AE samples were sequenced and analyzed by geno2pheno[coreceptor]. In parallel, coreceptor usage was determined for these samples by replicative phenotyping in human cells in the presence of specific X4- or R5-inhibitors. The sole introduction of the CRF01_AE V3 region into a full-length otherwise subtype B provirus failed to produce replication-competent viral progeny. A successive genome-replacement strategy revealed that also CRF01_AE derived gag and pol sequences are necessary to generate HIV genomes with sufficient replication competence. Subsequent phenotypic analysis confirmed overcalling of X4-tropism for CRF01_AE viruses using the current version and the standard cut-off at 10% false positive rate (FPR) of geno2pheno[coreceptor]. Lowering the FPR cut-off to 2.5% reduced the X4-overcalling in our sample collection, while still allowing a safe administration of Maraviroc (MCV). This study demonstrates the successful adjustment of geno2pheno[coreceptor] rules for subtype CRF01_AE. It also supports the unique strength of combining complementing methods, namely phenotyping and genotyping, for validating new bioinformatics tools prior to application in diagnostics.

Identification of novel molecular determinants of co-receptor usage in HIV-1 subtype F V3 envelope sequences

HIV-1 determinants of coreceptor usage within the gp120 V3 loop have been broadly studied over the past years. This information has led to the development of state-of the-art bioinformatic tools that are useful to predict co-receptor usage based on the V3 loop sequence mainly of subtypes B, C and A. However, these methods show a poor performance for subtype F V3 loops, which are found in an increasing number of HIV-1 strains worldwide. In the present work we investigated determinants of viral tropisms in the understudied subtype F by looking at genotypic and structural information of coreceptor:V3 loop interactions in a novel group of 40 subtype F V3 loops obtained from HIV-1 strains phenotypically characterized either as syncytium inducing or non-syncytium inducing by the MT-2 assay. We provide novel information about estimated interactions energies between a set of V3 loops with known tropism in subtype F, that allowed us to improve predictions of the coreceptor usage for this subtype. Understanding genetic and structural features underlying HIV coreceptor usage across different subtypes is relevant for the rational design of preventive and therapeutic strategies aimed at limiting the HIV-1 epidemic worldwide.

Molecular and geographic characterization of hiv-1 bf recombinant viruses

The Human Immunodeficiency Virus Type 1 (HIV-1) presents a wide genetic variability, which is represented by four groups, nine subtypes of group M and several recombinant forms. Among these, the BF recombinants have been distinguished by a high global dispersion and an increase in number and diversity. To date, 15 BF Circulating Recombinant Forms (CRFs) and diverse BF Unique Recombinant Forms (URFs) have been described. In Brazil, nine CRF_BF have been identified. The aim of this work was to perform molecular and geographic characterization of HIV-1 BF recombinant strains. Near full-length genomes of 265 BF recombinant viruses were collected from public databases and molecular analyses were performed. These sequences were originally retrieved between 1993–2006 and isolated from 16 countries (51.3% from Brazil). Diagnostic's year analysis showed that BF recombinants circulate in Brazil since at least 1985. Most sequences displayed recombination in the pol (84.9%), gag (69.3%) and env (51.4%) regions. The subtype B predominated in all accessory and regulatory genes, except in vif, in which the F subtype was predominant (40.4%). Twelve regions with a recombination rate higher than 10% were identified, especially one region inside p24 gene (1359–1397) whose recombination was present in more than 30% of the sequences. Coreceptor usage prediction during viral entry showed that BF recombinants preferentially use CCR5 (67.2%) and the most frequent tetrapeptides found in the V3 loop were GPGR (47.9%) and GPGQ (21.1%). The frequency of X4/dual viruses was lower amongst F subtype (25.8%) V3 sequences, compared with B subtype (43%). In addition, mutations associated with intermediate or high resistance levels to PI (10.6%), NRTI (15.0%), NNRTI (14.0%) and INSTI (2.6%) were identified. The great diversity of the recombination patterns evidences that the recombination between the subtypes B and F is frequent, reflecting a probable high rate of dual infection and the acquisition of advantageous characteristics for viral fitness.

Performance of Geno2Pheno[coreceptor] to infer coreceptor use in human immunodeficiency virus type 1 (HIV-1) subtype A

BackgroundAssessment of human immunodeficiency virus type 1 (HIV-1) coreceptor usage is required prior to treatment with the CCR5 antagonist maraviroc to exclude the presence of CXCR4-using (X4) strains. Genotype-based interpretation systems are mostly designed on subtype B and have been reported to be less accurate for subtype A/CRF02_AG. ObjectivesTo evaluate the performance of the widely used Geno2Pheno[coreceptor] (G2P[c]) algorithm for prediction of coreceptor usage with subtype A/CRF02_AG vs. subtype B. Study designCo-receptor tropism of 24 subtype A/CRF02_AG and 24 subtype B viruses was measured phenotypically by a homebrew single-cycle assay and genotypically by using G2P[c]. Samples with discrepant genotype-phenotype results were analyzed by next generation sequencing (NGS) and interpreted by the NGS Geno2Pheno algorithm (G2P[454]). ResultsAt 10% false positive rate (FPR), the G2P[c]/phenotype discordance rate was 12.5% (n = 3) for subtype A/CRF02_AG and 8.3% (n = 2) for subtype B. Minority X4 species escaping detection by bulk sequencing but documented by NGS explained the two subtype B and possibly one subtype A/CRF02_AG discordant case. The other two subtype A/CRF02_AG miscalled by G2P[c] could be explained by X4 overcalling at borderline FPR and/or by algorithm failure. DiscussionOur study did not demonstrate relevantly higher G2P[c] inaccuracy with subtype A/CRF02_AG with respect to subtype B. Genotype/phenotype discordances can be due to different reasons, including but not limited to, algorithm inaccuracy. Very large genotype/phenotype correlation panels are required to detect and explain the reason for any consistent difference in genotypic tropism prediction for subtype A/CRF02_AG vs. subtype B.

Individuals with HIV-1 Subtype C Infection and Cryptococcal Meningitis Exhibit Viral Genetic Intermixing of HIV-1 Between Plasma and Cerebrospinal Fluid and a High Prevalence of CXCR4-Using Variants.

The genotypic properties of human immunodeficiency virus type 1 (HIV-1) subtype C in individuals presenting with cryptococcal meningitis (CM) are not well established. Employing single-genome amplification as well as bulk PCR, cloning and sequencing strategies, we evaluated the genetic properties of HIV-1 subtype C env in 16 antiretroviral therapy-naive study participants with CM. Eleven of the 16 participants had matched blood plasma and cerebrospinal fluid (CSF) evaluated, with the rest having either a plasma or CSF sample evaluated. Before antiretroviral therapy initiation, matched plasma and CSF-derived env sequences of all 11 participants displayed genetic intermixing between the two compartments. Overall, 7 of the 16 (∼43.8%) participants harbored CXCR4-using variants in plasma and/or CSF, according to coreceptor usage prediction algorithms. This study suggests that HIV-1 subtype C genetic intermixing between peripheral blood and the central nervous system is common in individuals presenting with CM, and that CXCR4 usage is present in one or both compartments in approximately 44% of individuals.

After 18 months of antiretroviral therapy, total HIV DNA decreases more pronouncedly in patients infected by CRF01_AE than in those infected by subtype B and CRF07_BC.

Whether the amount of HIV DNA is associated with the subtype of HIV-1 after antiretroviral therapy (ART) has not been reported. In the present study, the amount of HIV DNA and RNA and CD4+T counts in blood and semen prior to and after 18 months of ART were compared in 48 patients infected by CRF01_AE, subtype B or CRF07_BC of HIV-1. Viral RNA was suppressed and CD4 cell count recovery achieved in all patients. The level of HIV DNA were similar before ART; however, patients with CRF01_AE had less HIV DNA after ART than those with subtype B and CRF07_BC infection. According to prediction of co-receptor usage by Geno2Pheno and PSSM in combination, more than 35.6% of clones for CRF01_AE were predicted as CXCR4-using before ART, whereas less than 6% of those for subtype B and CRF07_BC were predicted as CXCR4-using. After 18 months of ART, no CXCR4-using clones were predicted in any of the subtypes. Despite more HIV RNA and fewer CD4 + T cells in patients with CRF01_AE before therapy, no significant differences (P > 0.05) in viral RNA or CD4 cell counts were observed between the subtypes after 18 months of ART. Thus, 18 months of antiretroviral therapy was more efficient in patients with CRF01_AE. Considering that successful ART dramatically reduces the viral load in both blood and semen, risks of sexual transmission of HIV were reduced, contributing to prevention of rapid spread of HIV among men who have sex with men in the region.

Prediction of coreceptor usage by five bioinformatics tools in a large Ethiopian HIV-1 subtype C cohort.

BackgroundGenotypic tropism testing (GTT) has been developed largely on HIV-1 subtype B. Although a few reports have analysed the utility of GTT in other subtypes, more studies using HIV-1 subtype C (HIV-1C) are needed, considering the huge contribution of HIV-1C to the global epidemic.MethodsPlasma was obtained from 420 treatment-naïve HIV-1C infected Ethiopians recruited 2009–2011. The V3 region was sequenced and the coreceptor usage was predicted by five tools: Geno2Pheno clinical–and clonal–models, PhenoSeq-C, C-PSSM and Raymond’s algorithm. The impact of baseline tropism on antiretroviral treatment (ART) outcome was evaluated.ResultsOf 352 patients with successful baseline V3 sequences, the proportion of predicted R5 virus varied between the methods by 12.5% (78.1%-90.6%). However, only 58.2% of the predictions were concordant and only 1.7% were predicted to be X4-tropic across the five methods. Compared pairwise, the highest concordance was between C-PSSM and Geno2Pheno clonal (86.4%). In bivariate intention to treat (ITT) analysis, R5 infected patients achieved treatment success more frequently than X4 infected at month six as predicted by Geno2Pheno clinical (77.8% vs 58.7%, P = 0.004) and at month 12 by C-PSSM (61.9% vs 46.6%, P = 0.038). However, in the multivariable analysis adjusted for age, gender, baseline CD4 and viral load, only tropism as predicted by C-PSSM showed an impact on month 12 (P = 0.04, OR 2.47, 95% CI 1.06–5.79).ConclusionEach of the bioinformatics models predicted R5 tropism with comparable frequency but there was a large discordance between the methods. Baseline tropism had an impact on outcome of first line ART at month 12 in multivariable ITT analysis but only based on prediction by C-PSSM which thus possibly could be used for predicting outcome of ART in HIV-1C infected Ethiopians.

Predicted coreceptor usage at end-stage HIV disease in tissues derived from subjects on antiretroviral therapy with an undetectable plasma viral load

HIV cure research is increasingly focused on anatomical tissues as sites for residual HIV replication during combined antiretroviral therapy (cART). Tissue-based HIV could contribute to low-level immune activation and viral rebound over the course of infection and could also influence the development of diseases, such as atherosclerosis, neurological disorders and cancers. cART-treated subjects have a decreased and irregular presence of HIV among tissues, which has resulted in a paucity of actual evidence concerning how or if HIV persists, replicates and evolves in various anatomical sites during therapy. In this study, we pooled 1806 HIV envelope V3 loop sequences from twenty-six tissue types (seventy-one total tissues) of six pre-cART subjects, four subjects with an unknown cART history who died with profound AIDS, and five subjects who died while on cART with an undetectable plasma viral load. A computational approach was used to assess sequences for their ability to utilize specific cellular coreceptors (R5, R5 and X4, or X4). We found that autopsied tissues obtained from virally suppressed cART+ subjects harbored both integrated and expressed viruses with similar coreceptor usage profiles to subjects with no or ineffective cART therapy (i.e., significant plasma viral load at death). The study suggests that tissue microenvironments provide a sanctuary for the continued evolution of HIV despite cART.

A genotypic method for determining HIV-2 coreceptor usage enables epidemiological studies and clinical decision support.

BackgroundCCR5-coreceptor antagonists can be used for treating HIV-2 infected individuals. Before initiating treatment with coreceptor antagonists, viral coreceptor usage should be determined to ensure that the virus can use only the CCR5 coreceptor (R5) and cannot evade the drug by using the CXCR4 coreceptor (X4-capable). However, until now, no online tool for the genotypic identification of HIV-2 coreceptor usage had been available. Furthermore, there is a lack of knowledge on the determinants of HIV-2 coreceptor usage. Therefore, we developed a data-driven web service for the prediction of HIV-2 coreceptor usage from the V3 loop of the HIV-2 glycoprotein and used the tool to identify novel discriminatory features of X4-capable variants.ResultsUsing 10 runs of tenfold cross validation, we selected a linear support vector machine (SVM) as the model for geno2pheno[coreceptor-hiv2], because it outperformed the other SVMs with an area under the ROC curve (AUC) of 0.95. We found that SVMs were highly accurate in identifying HIV-2 coreceptor usage, attaining sensitivities of 73.5% and specificities of 96% during tenfold nested cross validation. The predictive performance of SVMs was not significantly different (p value 0.37) from an existing rules-based approach. Moreover, geno2pheno[coreceptor-hiv2] achieved a predictive accuracy of 100% and outperformed the existing approach on an independent data set containing nine new isolates with corresponding phenotypic measurements of coreceptor usage. geno2pheno[coreceptor-hiv2] could not only reproduce the established markers of CXCR4-usage, but also revealed novel markers: the substitutions 27K, 15G, and 8S were significantly predictive of CXCR4 usage. Furthermore, SVMs trained on the amino-acid sequences of the V1 and V2 loops were also quite accurate in predicting coreceptor usage (AUCs of 0.84 and 0.65, respectively).ConclusionsIn this study, we developed geno2pheno[coreceptor-hiv2], the first online tool for the prediction of HIV-2 coreceptor usage from the V3 loop. Using our method, we identified novel amino-acid markers of X4-capable variants in the V3 loop and found that HIV-2 coreceptor usage is also influenced by the V1/V2 region. The tool can aid clinicians in deciding whether coreceptor antagonists such as maraviroc are a treatment option and enables epidemiological studies investigating HIV-2 coreceptor usage. geno2pheno[coreceptor-hiv2] is freely available at http://coreceptor-hiv2.geno2pheno.org.Electronic supplementary materialThe online version of this article (doi:10.1186/s12977-016-0320-7) contains supplementary material, which is available to authorized users.

Maraviroc in addition to cART during primary HIV infection: Results from MAIN randomized clinical trial and 96-weeks follow-up

BackgroundMulti-targeted treatment strategies including maraviroc (MVC) during Primary HIV Infection (PHI) may benefit from the immune-modulatory properties of this CCR5-inhibitor. ObjectivesWe conducted a proof-of-concept clinical trial aimed at assessing whether maraviroc in addition of a combination antiretroviral therapy (cART) initiated during PHI would improve immunological and virological parameters. Study designThe MAIN (Maraviroc in HIV Acute INfection) study was a randomized open-label clinical trial (EUDRACT number: 2008-007004-29) which enrolled 29 patients with PHI. Subjects were randomly assigned to receive cART-only (cART), cART+8 weeks of MVC (ST-MVC) or cART+48 weeks of MVC (LT-MVC), regardless of predicted co-receptor usage. After 48 weeks patients in ST-MVC and LT-MVC groups discontinued MVC. Patients were evaluated at week 48 and at week 96 of follow-up to assess differences in CD4 T-cell gain and plasma HIV-RNA. ResultsTwenty-nine patients were enrolled. Seven patients (24%) had a predicted CXCR4 co-receptor usage. At week 48, 27 patients (93.1%) reached HIV-RNA<50cps/mL. Median CD4 T-cell count increase was 313 cells/μL (p<0.001, Wilcoxon signed-rank test). At multivariate linear regression analysis, LT-MVC arm had the greatest CD4 T-cell increase, while patients in ST-MVC arm had the least gain in CD4 T-cells (p=0.007). At week 96, multivariate analysis showed no associations between former treatment arm and CD4 T-cell gain. ConclusionsThe MAIN study showed that MVC for 48 weeks in addition to cART during PHI was able to enhance CD4 T-cell gain, regardless of co-receptor usage. After MVC discontinuation, the difference between treatment arms was lost.

On the Physicochemical and Structural Modifications Associated with HIV-1 Subtype B Tropism Transition.

HIV-1 enters immune cells via binding the viral envelope to a host cell CD4 receptor, and then a secondary co-receptor, usually CCR5 (R5) or CXCR4 (X4), and some HIV can utilize both co-receptors (R5X4). Although a small set of amino-acid properties such as charge and sequence length applied to HIV-1 envelope V3 loop sequence data can be used to predict co-receptor usage, we sought to expand the fundamental understanding of the physiochemical basis of tropism by analyzing many, perhaps less obvious, amino-acid properties over a diverse array of HIV sequences. We examined 74 amino-acid physicochemical scales over 1,559 V3 loop sequences with biologically tested tropisms downloaded from the Los Alamos HIV sequence database. Linear regressions were then calculated for each feature relative to three tropism transitions (R5→X4; R5→R5X4; R5X4→X4). Independent correlations were rank ordered to determine informative features. A structural analysis of the V3 loop was performed to better interpret these findings relative to HIV tropism states. Similar structural changes are required for R5 and R5X4 to transition to X4, thus suggesting that R5 and R5X4 types are more similar than either phenotype is to X4. Overall, the analysis suggests a continuum of viral tropism that is only partially related to charge; in fact, the analysis suggests that charge modification may be primarily attributed to decreased R5 usage, and further structural changes, particularly those associated with β-sheet structure, are likely required for full X4 usage.

A putative HIV-1 subtype C/CRF11_cpx unique recombinant from South Africa.

The HIV epidemic in South Africa is overwhelmingly driven by HIV-1 subtype C viruses. The HIV gag, pol, env (C2-V5) and nef sequences of virus 08MB26ZA, obtained from a 47 year old woman, were studied by phylogenetic analysis, REGA and the jumping Profile Hidden Markov Model (jPHMM) tools. The pol gene was further analyzed for recombination by Simplot. The pol and env sequences were examined for genetic drug resistance mutations and predicted co-receptor usage respectively. There was agreement in the assignment of the gag sequence as pure HIV-1 subtype C by phylogenetic, REGA and jPHMM analyses. The pol sequence clustered with CRF11_cpx and in the J-clade by phylogenetic analysis; and to a CRF11_cpx/subtype C recombinant by REGA. The assignment of pol to CRF11_cpx and subtype C was confirmed by Simplot. The recombinant was of the R5 biotype, with no important drug resistance mutations in the pol region. The epidemiologic and biologic significance of the virus are unknown. The finding suggests that complex viruses are being introduced into South Africa with potential implications for diagnosis. This is apparently the first report from South Africa of a putative unique recombinant involving CRF11_cpx and subtype C genomes.

Highly Accurate Structure-Based Prediction of HIV-1 Coreceptor Usage Suggests Intermolecular Interactions Driving Tropism.

HIV-1 entry into host cells is mediated by interactions between the V3-loop of viral glycoprotein gp120 and chemokine receptor CCR5 or CXCR4, collectively known as HIV-1 coreceptors. Accurate genotypic prediction of coreceptor usage is of significant clinical interest and determination of the factors driving tropism has been the focus of extensive study. We have developed a method based on nonlinear support vector machines to elucidate the interacting residue pairs driving coreceptor usage and provide highly accurate coreceptor usage predictions. Our models utilize centroid-centroid interaction energies from computationally derived structures of the V3-loop:coreceptor complexes as primary features, while additional features based on established rules regarding V3-loop sequences are also investigated. We tested our method on 2455 V3-loop sequences of various lengths and subtypes, and produce a median area under the receiver operator curve of 0.977 based on 500 runs of 10-fold cross validation. Our study is the first to elucidate a small set of specific interacting residue pairs between the V3-loop and coreceptors capable of predicting coreceptor usage with high accuracy across major HIV-1 subtypes. The developed method has been implemented as a web tool named CRUSH, CoReceptor USage prediction for HIV-1, which is available at http://ares.tamu.edu/CRUSH/.

Genetic Analyses of HIV-1 Strains Transmitted from Mother to Child in Northern Vietnam.

We previously reported mother-to-child transmission of HIV-1 in nine (6.7%) of 135 children on nevirapine prophylaxis in Vietnam. In the current study, we investigated the appearance and profile of antiretroviral drug (ARV) resistance mutations, the predicted coreceptor usage, and the genetic diversity of HIV-1 strains isolated from the eight pairs of HIV-1-infected mothers and their children, who were followed up to 12 months after birth. Portions of the pol and env C2V3 regions of the HIV-1 strains were analyzed genetically. HIV-1 CRF01_AE RNA was detected in four (50%) children at delivery. Y181C, a nevirapine resistance mutation, appeared in two (25%) children 1 and 3 months after birth, respectively. No ARV resistance mutation was detected in the mothers, though three mothers were on ARV prophylaxis. Five mothers and their children harbored CCR5-tropic (R5) viruses. Two mothers harbored both R5 and CXCR4-tropic (X4) viruses, but their children harbored only R5 viruses even though the X4 viruses were dominant in the mothers. In the remaining one mother, HIV-1 RNA was not amplified and her child harbored both R5 and X4 viruses at birth, but only X4 virus 12 months after delivery. The infants' viruses were more homogeneous than their mothers' viruses (mean distance: 0.5% vs. 1.1%, respectively). This is the first molecular epidemiological study of vertical HIV-1 infections in Vietnam. These findings may provide useful knowledge for the prevention of mother-to-child transmission of HIV-1 and the antiretroviral treatment of children in Vietnam.

Effect of lysine to arginine mutagenesis in the V3 loop of HIV-1 gp120 on viral entry efficiency and neutralization.

HIV-1 infection is characterized by an ongoing replication leading to T-lymphocyte decline which is paralleled by the switch from CCR5 to CXCR4 coreceptor usage. To predict coreceptor usage, several computer algorithms using gp120 V3 loop sequence data have been developed. In these algorithms an occupation of the V3 positions 11 and 25, by one of the amino acids lysine (K) or arginine (R), is an indicator for CXCR4 usage. Amino acids R and K dominate at these two positions, but can also be identified at positions 9 and 10. Generally, CXCR4-viruses possess V3 sequences, with an overall positive charge higher than the V3 sequences of R5-viruses. The net charge is calculated by subtracting the number of negatively charged amino acids (D, aspartic acid and E, glutamic acid) from the number of positively charged ones (K and R). In contrast to D and E, which are very similar in their polar and acidic properties, the characteristics of the R guanidinium group differ significantly from the K ammonium group. However, in coreceptor predictive computer algorithms R and K are both equally rated. The study was conducted to analyze differences in infectivity and coreceptor usage because of R-to-K mutations at the V3 positions 9, 10 and 11. V3 loop mutants with all possible RRR-to-KKK triplets were constructed and analyzed for coreceptor usage, infectivity and neutralization by SDF-1α and RANTES. Virus mutants R9R10R11 showed the highest infectivity rates, and were inhibited more efficiently in contrast to the K9K10K11 viruses. They also showed higher efficiency in a virus-gp120 paired infection assay. Especially V3 loop position 9 was relevant for a switch to higher infectivity when occupied by R. Thus, K-to-R exchanges play a role for enhanced viral entry efficiency and should therefore be considered when the viral phenotype is predicted based on V3 sequence data.

Covariance of charged amino acids at positions 322 and 440 of HIV-1 Env contributes to coreceptor specificity of subtype B viruses, and can be used to improve the performance of V3 sequence-based coreceptor usage prediction algorithms.

The ability to determine coreceptor usage of patient-derived human immunodeficiency virus type 1 (HIV-1) strains is clinically important, particularly for the administration of the CCR5 antagonist maraviroc. The envelope glycoprotein (Env) determinants of coreceptor specificity lie primarily within the gp120 V3 loop region, although other Env determinants have been shown to influence gp120-coreceptor interactions. Here, we determined whether conserved amino acid alterations outside the V3 loop that contribute to coreceptor usage exist, and whether these alterations improve the performance of V3 sequence-based coreceptor usage prediction algorithms. We demonstrate a significant covariant association between charged amino acids at position 322 in V3 and position 440 in the C4 Env region that contributes to the specificity of HIV-1 subtype B strains for CCR5 or CXCR4. Specifically, positively charged Lys/Arg at position 322 and negatively charged Asp/Glu at position 440 occurred more frequently in CXCR4-using viruses, whereas negatively charged Asp/Glu at position 322 and positively charged Arg at position 440 occurred more frequently in R5 strains. In the context of CD4-bound gp120, structural models suggest that covariation of amino acids at Env positions 322 and 440 has the potential to alter electrostatic interactions that are formed between gp120 and charged amino acids in the CCR5 N-terminus. We further demonstrate that inclusion of a “440 rule” can improve the sensitivity of several V3 sequence-based genotypic algorithms for predicting coreceptor usage of subtype B HIV-1 strains, without compromising specificity, and significantly improves the AUROC of the geno2pheno algorithm when set to its recommended false positive rate of 5.75%. Together, our results provide further mechanistic insights into the intra-molecular interactions within Env that contribute to coreceptor specificity of subtype B HIV-1 strains, and demonstrate that incorporation of Env determinants outside V3 can improve the reliability of coreceptor usage prediction algorithms.

GCUP: rapid GPU-based HIV-1 co-receptor usage prediction for next-generation sequencing.

Next-generation sequencing (NGS) has a large potential in HIV diagnostics, and genotypic prediction models have been developed and successfully tested in the recent years. However, albeit being highly accurate, these computational models lack computational efficiency to reach their full potential. In this study, we demonstrate the use of graphics processing units (GPUs) in combination with a computational prediction model for HIV tropism. Our new model named gCUP, parallelized and optimized for GPU, is highly accurate and can classify >175 000 sequences per second on an NVIDIA GeForce GTX 460. The computational efficiency of our new model is the next step to enable NGS technologies to reach clinical significance in HIV diagnostics. Moreover, our approach is not limited to HIV tropism prediction, but can also be easily adapted to other settings, e.g. drug resistance prediction. The source code can be downloaded at http://www.heiderlab.de d.heider@wz-straubing.de.

Bioinformatic analysis of HIV-1 entry and pathogenesis.

The evolution of human immunodeficiency virus type 1 (HIV-1) with respect to co-receptor utilization has been shown to be relevant to HIV-1 pathogenesis and disease. The CCR5-utilizing (R5) virus has been shown to be important in the very early stages of transmission and highly prevalent during asymptomatic infection and chronic disease. In addition, the R5 virus has been proposed to be involved in neuroinvasion and central nervous system (CNS) disease. In contrast, the CXCR4-utilizing (X4) virus is more prevalent during the course of disease progression and concurrent with the loss of CD4(+) T cells. The dual-tropic virus is able to utilize both co-receptors (CXCR4 and CCR5) and has been thought to represent an intermediate transitional virus that possesses properties of both X4 and R5 viruses that can be encountered at many stages of disease. The use of computational tools and bioinformatic approaches in the prediction of HIV-1 co-receptor usage has been growing in importance with respect to understanding HIV-1 pathogenesis and disease, developing diagnostic tools, and improving the efficacy of therapeutic strategies focused on blocking viral entry. Current strategies have enhanced the sensitivity, specificity, and reproducibility relative to the prediction of co-receptor use; however, these technologies need to be improved with respect to their efficient and accurate use across the HIV-1 subtypes. The most effective approach may center on the combined use of different algorithms involving sequences within and outside of the env-V3 loop. This review focuses on the HIV-1 entry process and on co-receptor utilization, including bioinformatic tools utilized in the prediction of co-receptor usage. It also provides novel preliminary analyses for enabling identification of linkages between amino acids in V3 with other components of the HIV-1 genome and demonstrates that these linkages are different between X4 and R5 viruses.

Determination of viral tropism by genotyping and phenotyping assays in Brazilian HIV-1-infected patients.

The clinical application of CCR5 antagonists involves first determining the coreceptor usage by the infecting viral strain. Bioinformatics programs that predict coreceptor usage could provide an alternative method to screen candidates for treatment with CCR5 antagonists, particularly in countries with limited financial resources. Thus, the present study aims to identify the best approach using bioinformatics tools for determining HIV-1 coreceptor usage in clinical practice. Proviral DNA sequences and Trofile results from 99 HIV-1-infected subjects under clinical monitoring were analyzed in this study. Based on the Trofile results, the viral variants present were 81.1% R5, 21.4% R5X4 and 1.8% X4. Determination of tropism using a Geno2pheno[coreceptor] analysis with a false positive rate of 10% gave the most suitable performance in this sampling: the R5 and X4 strains were found at frequencies of 78.5% and 28.4%, respectively, and there was 78.6% concordance between the phenotypic and genotypic results. Further studies are needed to clarify how genetic diversity amongst virus strains affects bioinformatics-driven approaches for determining tropism. Although this strategy could be useful for screening patients in developing countries, some limitations remain that restrict the wider application of coreceptor usage tests in clinical practice.

Profile of HIV Type 1 Coreceptor Tropism Among Kenyan Patients from 2009 to 2010

A switch of HIV coreceptor usage from CCR5 to CXCR4 occurs in AIDS pathogenesis and may play a critical role in the use of entry inhibitors. To determine the potential usefulness of maraviroc and other CCR5 antagonists among drug-naive and experienced patients in Kenya, the env-C2-V3 gene was successfully sequenced in samples from 176 (98 men and 78 female) consenting subjects between January 2009 and December 2012. In silico CPSSM, webPSSM/, and (ds) Kernel tools were used in predicting coreceptor usage. On the basis of the env V3 loop sequences, 84.1% (148) were reported with R-5 tropism, 4.5% (5) were dual tropic, while 13.4% (23) were of X4 tropism. However, similar to previous studies conducted in Kenya on genetic diversity, HIV-1 subtype A1 (73.9%; 130/176) still remains the most dominant subtype. The high levels of R5 tropism among the studied Kenyan infected populations suggested the potential use of CCR5 antagonists as new therapeutic options in Kenya.