A Genotypic Test for HIV-1 Tropism Combining Sanger Sequencing with Ultradeep Sequencing Predicts Virologic Response in Treatment-Experienced Patients

A tropism test is required prior to initiation of CCR5 antagonist therapy in HIV-1 infected individuals, as these agents are not effective in patients harboring CXCR4 (X4) coreceptor-using viral variants. We developed a clinical laboratory-based genotypic tropism test for detection of CCR5-using (R5) or X4 variants that utilizes triplicate population sequencing (TPS) followed by ultradeep sequencing (UDS) for samples classified as R5. Tropism was inferred using the bioinformatic algorithms geno2pheno_[coreceptor] and PSSM_x4r5. Virologic response as a function of tropism readout was retrospectively assessed using blinded samples from treatment-experienced subjects who received maraviroc (N = 327) in the MOTIVATE and A4001029 clinical trials. MOTIVATE patients were classified as R5 and A4001029 patients were classified as non-R5 by the original Trofile test. Virologic response was compared between the R5 and non-R5 groups determined by TPS, UDS alone, the reflex strategy and the Trofile Enhanced Sensitivity (TF-ES) test. UDS had greater sensitivity than TPS to detect minority non-R5 variants. The median log₁₀ viral load change at week 8 was −2.4 for R5 subjects, regardless of the method used for classification; for subjects with non-R5 virus, median changes were −1.2 for TF-ES or the Reflex Test and −1.0 for UDS. The differences between R5 and non-R5 groups were highly significant in all 3 cases (p<0.0001). At week 8, the positive predictive value was 66% for TF-ES and 65% for both the Reflex test and UDS. Negative predictive values were 59% for TF-ES, 58% for the Reflex Test and 61% for UDS. In conclusion, genotypic tropism testing using UDS alone or a reflex strategy separated maraviroc responders and non-responders as well as a sensitive phenotypic test, and both assays showed improved performance compared to TPS alone. Genotypic tropism tests may provide an alternative to phenotypic testing with similar discriminating ability.     

Genotypic Prediction of Tropism of Highly Diverse HIV-1 Strains from Cameroon

Background

The use of CCR5 antagonists involves determination of HIV-1 tropism prior to initiation of treatment. HIV-1 tropism can be assessed either by phenotypic or genotypic methods. Genotypic methods are extensively used for tropism prediction. However, their validation in predicting tropism of viral isolates belonging to group M non-B subtypes remains challenging. In Cameroon, the genetic diversity of HIV-1 strains is the broadest reported worldwide. To facilitate the integration of CCR5 antagonists into clinical practice in this region, there is a need to evaluate the performance of genotypic methods for predicting tropism of highly diverse group M HIV-1 strains.

Methods

Tropism of diverse HIV-1 strains isolated from PBMCs from Cameroon was determined using the GHOST cell assay. Prediction, based on V3 sequences from matched plasma samples, was determined using bioinformatics algorithms and rules based on position 11/25 and net charge applied independently or combined according to Delobel''s and Garrido''s rules. Performance of genotypic methods was evaluated by comparing prediction generated with tropism assigned by the phenotypic assay.

Results

Specificity for predicting R5-tropic virus was high, ranging from 83.7% to 97.7% depending on the genotypic methods used. Sensitivity for X4-tropic viruses was fairly low, ranging from 33.3% to 50%. In our study, overall, genotypic methods were less able to accurately predict X4-tropic virus belonging to subtype CRF02_AG. In addition, it was found that of the methods we used the Garrido rule has the highest sensitivity rate of over 50% with a specificity of 93%.

Conclusion

Our study demonstrated that overall, genotypic methods were less sensitive for accurate prediction of HIV-1 tropism in settings where diverse HIV-1 strains co-circulate. Our data suggest that further optimization of genotypic methods is needed and that larger studies to determine their utility for tropism prediction of diverse HIV-1 strains may be warranted.     

Use of Four Next-Generation Sequencing Platforms to Determine HIV-1 Coreceptor Tropism

HIV-1 coreceptor tropism assays are required to rule out the presence of CXCR4-tropic (non-R5) viruses prior treatment with CCR5 antagonists. Phenotypic (e.g., Trofile™, Monogram Biosciences) and genotypic (e.g., population sequencing linked to bioinformatic algorithms) assays are the most widely used. Although several next-generation sequencing (NGS) platforms are available, to date all published deep sequencing HIV-1 tropism studies have used the 454™ Life Sciences/Roche platform. In this study, HIV-1 co-receptor usage was predicted for twelve patients scheduled to start a maraviroc-based antiretroviral regimen. The V3 region of the HIV-1 env gene was sequenced using four NGS platforms: 454™, PacBio® RS (Pacific Biosciences), Illumina®, and Ion Torrent™ (Life Technologies). Cross-platform variation was evaluated, including number of reads, read length and error rates. HIV-1 tropism was inferred using Geno2Pheno, Web PSSM, and the 11/24/25 rule and compared with Trofile™ and virologic response to antiretroviral therapy. Error rates related to insertions/deletions (indels) and nucleotide substitutions introduced by the four NGS platforms were low compared to the actual HIV-1 sequence variation. Each platform detected all major virus variants within the HIV-1 population with similar frequencies. Identification of non-R5 viruses was comparable among the four platforms, with minor differences attributable to the algorithms used to infer HIV-1 tropism. All NGS platforms showed similar concordance with virologic response to the maraviroc-based regimen (75% to 80% range depending on the algorithm used), compared to Trofile (80%) and population sequencing (70%). In conclusion, all four NGS platforms were able to detect minority non-R5 variants at comparable levels suggesting that any NGS-based method can be used to predict HIV-1 coreceptor usage.     

Prediction of HIV-1 Coreceptor Usage (Tropism) by Sequence Analysis using a Genotypic Approach

Maraviroc (MVC) is the first licensed antiretroviral drug from the class of coreceptor antagonists. It binds to the host coreceptor CCR5, which is used by the majority of HIV strains in order to infect the human immune cells (Fig. 1). Other HIV isolates use a different coreceptor, the CXCR4. Which receptor is used, is determined in the virus by the Env protein (Fig. 2). Depending on the coreceptor used, the viruses are classified as R5 or X4, respectively. MVC binds to the CCR5 receptor inhibiting the entry of R5 viruses into the target cell. During the course of disease, X4 viruses may emerge and outgrow the R5 viruses. Determination of coreceptor usage (also called tropism) is therefore mandatory prior to administration of MVC, as demanded by EMA and FDA. The studies for MVC efficiency MOTIVATE, MERIT and 1029 have been performed with the Trofile assay from Monogram, San Francisco, U.S.A. This is a high quality assay based on sophisticated recombinant tests. The acceptance for this test for daily routine is rather low outside of the U.S.A., since the European physicians rather tend to work with decentralized expert laboratories, which also provide concomitant resistance testing. These laboratories have undergone several quality assurance evaluations, the last one being presented in 2011¹.For several years now, we have performed tropism determinations based on sequence analysis from the HIV env-V3 gene region (V3)². This region carries enough information to perform a reliable prediction. The genotypic determination of coreceptor usage presents advantages such as: shorter turnover time (equivalent to resistance testing), lower costs, possibility to adapt the results to the patients'' needs and possibility of analysing clinical samples with very low or even undetectable viral load (VL), particularly since the number of samples analysed with VL<1000 copies/μl roughly increased in the last years (Fig. 3).The main steps for tropism testing (Fig. 4) demonstrated in this video: 1. Collection of a blood sample 2. Isolation of the HIV RNA from the plasma and/or HIV proviral DNA from blood mononuclear cells 3. Amplification of the env region 4. Amplification of the V3 region 5. Sequence reaction of the V3 amplicon 6. Purification of the sequencing samples 7. Sequencing the purified samples 8. Sequence editing 9. Sequencing data interpretation and tropism prediction     

Parameters Influencing Baseline HIV-1 Genotypic Tropism Testing Related to Clinical Outcome in Patients on Maraviroc

ObjectivesWe analysed the impact of different parameters on genotypic tropism testing related to clinical outcome prediction in 108 patients on maraviroc (MVC) treatment.Methods87 RNA and 60 DNA samples were used. The viral tropism was predicted using the geno2pheno_[coreceptor] and T-CUP tools with FPR cut-offs ranging from 1%-20%. Additionally, 27 RNA and 28 DNA samples were analysed in triplicate, 43 samples with the ESTA assay and 45 with next-generation sequencing. The influence of the genotypic susceptibility score (GSS) and 16 MVC-resistance mutations on clinical outcome was also studied.ResultsConcordance between single-amplification testing compared to ESTA and to NGS was in the order of 80%. Concordance with NGS was higher at lower FPR cut-offs. Detection of baseline R5 viruses in RNA and DNA samples by all methods significantly correlated with treatment success, even with FPR cut-offs of 3.75%-7.5%. Triple amplification did not improve the prediction value but reduced the number of patients eligible for MVC. No influence of the GSS or MVC-resistance mutations but adherence to treatment, on the clinical outcome was detected.ConclusionsProviral DNA is valid to select candidates for MVC treatment. FPR cut-offs of 5%-7.5% and single amplification from RNA or DNA would assure a safe administration of MVC without excluding many patients who could benefit from this drug. In addition, the new prediction system T-CUP produced reliable results.     

Phylodynamic Inference of Bacterial Outbreak Parameters Using Nanopore Sequencing

Nanopore sequencing and phylodynamic modeling have been used to reconstruct the transmission dynamics of viral epidemics, but their application to bacterial pathogens has remained challenging. Cost-effective bacterial genome sequencing and variant calling on nanopore platforms would greatly enhance surveillance and outbreak response in communities without access to sequencing infrastructure. Here, we adapt random forest models for single nucleotide polymorphism (SNP) polishing developed by Sanderson and colleagues (2020. High precision Neisseria gonorrhoeae variant and antimicrobial resistance calling from metagenomic nanopore sequencing. Genome Res. 30(9):1354–1363) to estimate divergence and effective reproduction numbers (R_e) of two methicillin-resistant Staphylococcus aureus (MRSA) outbreaks from remote communities in Far North Queensland and Papua New Guinea (PNG; n = 159). Successive barcoded panels of S. aureus isolates (2 × 12 per MinION) sequenced at low coverage (>5× to 10×) provided sufficient data to accurately infer genotypes with high recall when compared with Illumina references. Random forest models achieved high resolution on ST93 outbreak sequence types (>90% accuracy and precision) and enabled phylodynamic inference of epidemiological parameters using birth–death skyline models. Our method reproduced phylogenetic topology, origin of the outbreaks, and indications of epidemic growth (R_e > 1). Nextflow pipelines implement SNP polisher training, evaluation, and outbreak alignments, enabling reconstruction of within-lineage transmission dynamics for infection control of bacterial disease outbreaks on portable nanopore platforms. Our study shows that nanopore technology can be used for bacterial outbreak reconstruction at competitive costs, providing opportunities for infection control in hospitals and communities without access to sequencing infrastructure, such as in remote northern Australia and PNG.     

HIV-1 Tropism Dynamics and Phylogenetic Analysis from Longitudinal Ultra-Deep Sequencing Data of CCR5- and CXCR4-Using Variants

Objective

Coreceptor switch from CCR5 to CXCR4 is associated with HIV disease progression. The molecular and evolutionary mechanisms underlying the CCR5 to CXCR4 switch are the focus of intense recent research. We studied the HIV-1 tropism dynamics in relation to coreceptor usage, the nature of quasispecies from ultra deep sequencing (UDPS) data and their phylogenetic relationships.

Methods

Here, we characterized C2-V3-C3 sequences of HIV obtained from 19 patients followed up for 54 to 114 months using UDPS, with further genotyping and phylogenetic analysis for coreceptor usage. HIV quasispecies diversity and variability as well as HIV plasma viral load were measured longitudinally and their relationship with the HIV coreceptor usage was analyzed. The longitudinal UDPS data were submitted to phylogenetic analysis and sampling times and coreceptor usage were mapped onto the trees obtained.

Results

Although a temporal viral genetic structuring was evident, the persistence of several viral lineages evolving independently along the infection was statistically supported, indicating a complex scenario for the evolution of viral quasispecies. HIV X4-using variants were present in most of our patients, exhibiting a dissimilar inter- and intra-patient predominance as the component of quasispecies even on antiretroviral therapy. The viral populations from some of the patients studied displayed evidences of the evolution of X4 variants through fitness valleys, whereas for other patients the data favored a gradual mode of emergence.

Conclusions

CXCR4 usage can emerge independently, in multiple lineages, along the course of HIV infection. The mode of emergence, i.e. gradual or through fitness valleys seems to depend on both virus and patient factors. Furthermore, our analyses suggest that, besides becoming dominant after population-level switches, minor proportions of X4 viruses might exist along the infection, perhaps even at early stages of it. The fate of these minor variants might depend on both viral and host factors.     

HIV-1 GP120 V3 conformational and informational entropies

In an attempt to analyze structure, function and evolution of HIV-1 GP120 V3, interactions among the Hartree–Fock energy, the conformational entropy and the Shannon entropy were determined for the 1NJ0 set of antibody-bound V3 loop conformers. The Hartree–Fock energy of each conformer was determined at the MINI level with GAMESS. The conformational entropy was determined per conformer and per residue from the mass-weighted covariance matrices. The Shannon entropy per residue was determined from sequence-substitution frequencies. Correlations were determined by linear regression analysis. There was a negative correlation between the Hartree–Fock energy and the conformational entropy (R=−0.4840, p=0.0078, df =28) that enhanced the negative Helmholtz-free-energy change for the binding of the GP120 ligand to target CD4. The Shannon entropy of V3 was a function of the conformational entropy variance (R=0.7225, p=0.00157, df=15) and of the V3 Hartree–Fock energy. Biological implications of this work are that (1) conformational entropy interacts with V3 Hartree–Fock energy to enhance GP120 binding to CD4 cell receptors and that (2) the Hartree–Fock energy of V3 interacts with the evolutionary system to participate in the regulation of V3 diversity.     

Biophysical characterization of V3-lipopeptide liposomes influencing HIV-1 infectivity

The V3-loop of the HIV-1 gp120 alters host cell immune function and modulates infectivity. We investigated biophysical parameters of liposome constructs with embedded lipopeptides from the principle neutralizing domain of the V3-loop and their influence on viral infectivity. Dynamic light scattering measurements showed liposome supramolecular structures with hydrodynamic radius of the order of 900 and 1300nm for plain and V3-lipopeptide liposomes. Electron paramagnetic resonance measurements showed almost identical local microenvironment. The difference in liposome hydrodynamic radius was attributed to the fluctuating ionic environment of the V3-lipopeptide liposomes. In vitro HIV-1 infectivity assays showed that plain liposomes reduced virus production in all cell cultures, probably due to the hydrophobic nature of the aggregates. Liposomes carrying V3-lipopeptides with different cationic potentials restored and even enhanced infectivity (p<0.05). These results highlight the need for elucidation of the involvement of lipid bilayers as dynamic components in supramolecular structures and in HIV-1 fusion mechanisms.     

HIV-1 Tropism Testing in Subjects Achieving Undetectable HIV-1 RNA: Diagnostic Accuracy,Viral Evolution and Compartmentalization

Background

Technically, HIV-1 tropism can be evaluated in plasma or peripheral blood mononuclear cells (PBMCs). However, only tropism testing of plasma HIV-1 has been validated as a tool to predict virological response to CCR5 antagonists in clinical trials. The preferable tropism testing strategy in subjects with undetectable HIV-1 viremia, in whom plasma tropism testing is not feasible, remains uncertain.

Methods & Results

We designed a proof-of-concept study including 30 chronically HIV-1-infected individuals who achieved HIV-1 RNA <50 copies/mL during at least 2 years after first-line ART initiation. First, we determined the diagnostic accuracy of 454 and population sequencing of gp120 V3-loops in plasma and PBMCs, as well as of MT-2 assays before ART initiation. The Enhanced Sensitivity Trofile Assay (ESTA) was used as the technical reference standard. 454 sequencing of plasma viruses provided the highest agreement with ESTA. The accuracy of 454 sequencing decreased in PBMCs due to reduced specificity. Population sequencing in plasma and PBMCs was slightly less accurate than plasma 454 sequencing, being less sensitive but more specific. MT-2 assays had low sensitivity but 100% specificity. Then, we used optimized 454 sequence data to investigate viral evolution in PBMCs during viremia suppression and only found evolution of R5 viruses in one subject. No de novo CXCR4-using HIV-1 production was observed over time. Finally, Slatkin-Maddison tests suggested that plasma and cell-associated V3 forms were sometimes compartmentalized.

Conclusions

The absence of tropism shifts during viremia suppression suggests that, when available, testing of stored plasma samples is generally safe and informative, provided that HIV-1 suppression is maintained. Tropism testing in PBMCs may not necessarily produce equivalent biological results to plasma, because the structure of viral populations and the diagnostic performance of tropism assays may sometimes vary between compartments. Thereby, proviral DNA tropism testing should be specifically validated in clinical trials before it can be applied to routine clinical decision-making.     

Membranotropic effects of peptides from the V3 loop of HIV-1

Summary The V3 loop from HIV-1 envelope glycoprotein gp120 is involved in viral entry and determines the cellular tropism and HIV-1-induced cell-cell fusion. Earlier we have shown that V3 loop peptides representing the sequences of syncytia-inducing HIV strains have high membranotropic activity. These peptides caused the lysis of liposomes of various lipid compositions, could fuse negatively charged liposomes and induced hemolysis of erythrocytes. In contrast, peptides mimicking the sequences of non-syncytia-inducing viruses showed no lytic or fusion activities at the same concentrations. Now we have found that the V3 loop synthetic peptides containing the conserved GPGR region, derived from T-lymphotropic strains (BRU and MN), as opposed to peptides containing the GPGQ region, are able to cause a pronounced membrane permeabilization (dissipation of the ΔpH and the Δψ) of human peripheral blood lymphocytes, erythrocytes and plasma membrane vesicles at micromolar concentrations with a dose-dependent kinetics. Analysis of the secondary structures of the peptides by circular dichroism revealed conformational changes in V3 loop peptides depending on solvent hydrophobicity: from random coil in water to an α-helix/β-sheet conformation in trifluoroethanol. Such structural changes of the V3 loop together with the membrane insertion of the gp41 N-terminal fusion peptide may promote the formation of the fusion pore during virus-cell fusion.     

Membranotropic effects of peptides from the V3 loop of HIV-1

The V3 loop from HIV-1 envelope glycoprotein gp120 is involved in viral entry and determines the cellular tropism and HIV-1-induced cell–cell fusion. Earlier we have shown that V3 loop peptides representing the sequences of syncytia-inducing HIV strains have high membranotropic activity. These peptides caused the lysis of liposomes of various lipid compositions, could fuse negatively charged liposomes and induced hemolysis of erythrocytes. In contrast, peptides mimicking the sequences of non-syncytia-inducing viruses showed no lytic or fusion activities at the same concentrations. Now we have found that the V3 loop synthetic peptides containing the conserved GPGR region, derived from T-lymphotropic strains (BRU and MN), as opposed to peptides containing the GPGQ region, are able to cause a pronounced membrane permeabilization (dissipation of the pH and the of human peripheral blood lymphocytes, erythrocytes and plasma membrane vesicles at micromolar concentrations with a dose-dependent kinetics. Analysis of the secondary structures of the peptides by circular dichroism revealed conformational changes in V3 loop peptides depending on solvent hydrophobicity: from random coil in water to an -helix/-sheet conformation in trifluoroethanol. Such structural changes of the V3 loop together with the membrane insertion of the gp41 N-terminal fusion peptide may promote the formation of the fusion pore during virus–cell fusion.     

Isolation of TAK-779-resistant HIV-1 from an R5 HIV-1 GP120 V3 loop library

The human immunodeficiency virus (HIV-1) envelope glycoprotein (GP) 120 interacts with CD4 and the CCR5 coreceptor for viral entry. The V3 loop in GP120 is a crucial region for determining coreceptor usage during viral entry, and a variety of amino acid substitutions has been observed in clinical isolates. To construct an HIV-1 V3 loop library, we chose 10 amino acid positions in the V3 loop and incorporated random combinations (27,648 possibilities) of the amino acid substitutions derived from 31 R5 viruses into the V3 loop of HIV-1(JR-FL) proviral DNA. The constructed HIV-1 library contained 6.6 x 10(6) independent clones containing a set of 0-10 amino acid substitutions in the V3 loop. To address whether restricted steric alteration in the V3 loop could confer resistance to an entry inhibitor, TAK-779, we selected entry inhibitor-resistant HIV-1 by increasing the concentration of TAK-779 from 0.10 to 0.30 microM in PM1-CCR5 cells with high expression of CCR5. The selected viruses at passage 8 contained five amino acid substitutions in the V3 loop without any other mutations in GP120 and showed 15-fold resistance compared with the parental virus. These results indicated that a certain structure of the V3 loop containing amino acid substitutions derived from 31 R5 viruses can contribute to the acquisition of resistance to entry inhibitors binding to CCR5. Taken together, this type of HIV-1 V3 loop library is useful for isolating and analyzing the specific biological features of HIV-1 with respect to alterations of the V3 loop structure.     

Evolution of the HIV-1 V3 region in the Italian epidemic

The epidemic of Human Immunodeficiency Virus type 1 infection in Italy is mostly ascribed to the B subtype, which represents the prevalent subtype in Western Countries. The virus isolates of the B subtype, moreover, show an increasing nucleotide heterogeneity over time, indicating a continuous intra-subtype dynamic evolution, typical of long-lasting epidemics. In recent years, however, the progressive decrease in the transmission rate among the historically defined risk groups (i.e. homosexuals and IDUs) and the parallel increase in heterosexual transmission are slowly introducing variants of non-B subtype into the Italian HIV-1 epidemic. This appears to be strictly linked to the growing number of immigrants from non-Western Countries, where non-B clades and CRFs are prevalent, and consequent inter-racial blending. The distribution of these novel genetic forms needs to be evaluated by continuous molecular monitoring nationwide to verify whether they will overcome the pre-existing B-clade epidemic, which could have significant implications for diagnosis, treatment and vaccine development. Here we review the genetic evolution of HIV-1 spreading within the Italian epidemic.     

HIV-1 Quasispecies Delineation by Tag Linkage Deep Sequencing

Trade-offs between throughput, read length, and error rates in high-throughput sequencing limit certain applications such as monitoring viral quasispecies. Here, we describe a molecular-based tag linkage method that allows assemblage of short sequence reads into long DNA fragments. It enables haplotype phasing with high accuracy and sensitivity to interrogate individual viral sequences in a quasispecies. This approach is demonstrated to deduce ∼2000 unique 1.3 kb viral sequences from HIV-1 quasispecies in vivo and after passaging ex vivo with a detection limit of ∼0.005% to ∼0.001%. Reproducibility of the method is validated quantitatively and qualitatively by a technical replicate. This approach can improve monitoring of the genetic architecture and evolution dynamics in any quasispecies population.     

Immunophilins and HIV-1 V3 Loop For Structure-Based Anti-AIDS Drug Design

Abstract

The model of the structural complex of cyclophilin A (CycA) belonging to the immunophilins family with the HIV-MN gpl20 V3 loop was generated, and the computer-aided design of the immunophilin-derived peptide able to mask the biologically crucial V3 segments was implemented.

To this end, the following problems were solved: (i) the NMR-based conformational analysis of the HIV-MN V3 loop was put into effect, and its low energy structure fitting the input experimental observations was determined; (ii) molecular docking of this V3 structure with the X-ray conformation of CycA was carried out, and the energy refining the simulated structural complex was performed; (iii) the matrix of inter-atomic distances for the amino acids of the molecules forming part of the built over-molecular ensemble was computed, the types of interactions responsible for its stabilization were analyzed, and the CycA stretch, which accounts for the binding to V3, was identified; (iv) the most probable 3D structure for this stretch in the unbound state was predicted, and its collation with the X-ray structure for the corresponding site of CycA was performed; (v) the potential energy function and its constituents were studied for the structural complex generated by molecular docking of the V3 loop with the CycA peptide offering the virtual molecule that imitates the CycA segment, making a key contribution to the interactions of the native protein with the HIV-1 principal neutralizing determinant; (vi) as a result of the studies above, the designed molecule was shown to be capable of the efficacious blockading the functionally crucial V3 sites; and (vii) based on the joint analysis of the evidence obtained previously and in the present study, the composition of the peptide cocktail presenting the promising anti-AIDS pharmacological substance was developed.

The molecules simulated here by molecular modeling methods may become the first representatives of a new class of the chemical compounds (immunophilin-derived peptides) offering the looking-forward basic structures for the design of efficacious and safe antiviral agents.     

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

Local structural properties of the V3 loop of Thailand HIV-1 isolate

The model of locally accurate conformation for the HIV-Thailand principal neutralizing determinant (PND) located within the V3 loop of the virus envelope protein gp120 was built in terms of NMR spectroscopy data. To this end, the NMR-based conformational analysis of synthetic molecule representing the peptide copy of the fragment under study was carried out using the published sequential d connectivity data and values of spin-spin coupling constants. As a result, (i) the local structure for the V3 loop from Thailand isolate was determined, (ii) the conformations of its irregular segments were analyzed, and the secondary structure elements identified, (iii) the ensemble of conformers matching the experimental and theoretical data was derived for the stretch forming the neutralizing epitope of the HIV-Thailand PND, (iv) to estimate the probability of realizing each of these conformers in solution, the results obtained were collated with the X-ray data for corresponding segments in synthetic molecules imitating the central region of the HIV-MN PND as well as for homologous segments 39-44 in Bence-Jonce REI protein (BJRP), 41-46 in immunoglobulin lambda (Ig lambda), and 50-55 in beta-chain of horse hemoglobin (HH), (v) to find the conserved structural motifs inside diverse HIV-1 isolates, the structure determined was compared with the one derived earlier for the HIV-MN PND from NMR spectroscopy data, (vi) on the basis of all data obtained, the 3D structure model describing the set of biologically relevant conformations, which may present different antigenic determinants to the immune system in various HIV-1 isolates, was proposed for the immunogenic crown of the V3 loop. The results obtained are discussed in conjunction with the data on the structure for the HIV-1 PND reported in literature.