Two-way Bayesian hierarchical phylogenetic models: An application to the co-evolution of gp120 and gp41 during and after enfuvirtide treatment

Assessment of Antibody Interference of Enfuvirtide (T20) Function Shows Assay Dependent Variability

Background:

During HIV infection, fusion of the viral and cellular membranes is dependent on folding of the gp41 trimer into a six-helix bundle. Fusion inhibitors, such as the antiretroviral Enfuvirtide (T20), interfere with the formation of the gp41 six-helix bundle. Recent in vitro studies reveal that the gp41 immunodominant region one targeting antibody 3D6 can block T20 interference, but the clinical and pathophysiologic significance of this finding is unclear.

Objective/Method:

We have previously characterized a number of antibodies that target conformational epitopes on gp41and herein characterized their ability to interfere with T20 in multiple assays and assess their prevalence in HIV infected subjects.

Results:

The T20 interference by antibody 3D6 was confirmed in a CHO-HXB2 envelope/ HeLaT4+ cell culture assay. Antibodies that target an immunodominant region one epitope, as well as a gp41 discontinuous epitope, also interfered in this assay, however, not all antibodies that targeted these epitopes showed T20 interference. This response was not due to the direct binding of T20 by the antibodies and could not be replicated utilizing TZM-bl and HL2/3 cells. Notably, serum competition studies on a panel of HIV subjects demonstrate that these conformational targeting antibodies are common in the HIV population.

Conclusion:

The relatively common nature of antibodies targeting these epitopes, the disparate in vitro results, and lack of reported clinical failures ascribed to such antibodies leads us to conclude that antibody interference of T20 is likely not clinically relevant. However, this warrants continued consideration with the advancement of other fusion inhibitors.

Vicriviroc resistance decay and relative replicative fitness in HIV-1 clinical isolates under sequential drug selection pressures

We previously described an HIV-1-infected individual who developed resistance to vicriviroc (VCV), an investigational CCR5 antagonist, during 28 weeks of therapy (Tsibris AM et al., J. Virol. 82:8210-8214, 2008). To investigate the decay of VCV resistance mutations, a standard clonal analysis of full-length env (gp160) was performed on plasma HIV-1 samples obtained at week 28 (the time of VCV discontinuation) and at three subsequent time points (weeks 30, 42, and 48). During 132 days, VCV-resistant HIV-1 was replaced by VCV-sensitive viruses whose V3 loop sequences differed from the dominant pretreatment forms. A deep-sequencing analysis showed that the week 48 VCV-sensitive V3 loop form emerged from a preexisting viral variant. Enfuvirtide was added to the antiretroviral regimen at week 30; by week 48, enfuvirtide treatment selected for either the G36D or N43D HR-1 mutation. Growth competition experiments demonstrated that viruses incorporating the dominant week 28 VCV-resistant env were less fit than week 0 viruses in the absence of VCV but more fit than week 48 viruses. This week 48 fitness deficit persisted when G36D was corrected by either site-directed mutagenesis or week 48 gp41 domain swapping. The correction of N43D, in contrast, restored fitness relative to that of week 28, but not week 0, viruses. Virus entry kinetics correlated with observed fitness differences; the slower entry of enfuvirtide-resistant viruses corrected to wild-type rates in the presence of enfuvirtide. These findings suggest that while VCV and enfuvirtide select for resistance mutations in only one env subunit, gp120 and gp41 coevolve to maximize viral fitness under sequential drug selection pressures.

Impact of CCR5delta32 host genetic background and disease progression on HIV-1 intrahost evolutionary processes: efficient hypothesis testing through hierarchical phylogenetic models

The interplay between C-C chemokine receptor type 5 (CCR5) host genetic background, disease progression, and intrahost HIV-1 evolutionary dynamics remains unclear because differences in viral evolution between hosts limit the ability to draw conclusions across hosts stratified into clinically relevant populations. Similar inference problems are proliferating across many measurably evolving pathogens for which intrahost sequence samples are readily available. To this end, we propose novel hierarchical phylogenetic models (HPMs) that incorporate fixed effects to test for differences in dynamics across host populations in a formal statistical framework employing stochastic search variable selection and model averaging. To clarify the role of CCR5 host genetic background and disease progression on viral evolutionary patterns, we obtain gp120 envelope sequences from clonal HIV-1 variants isolated at multiple time points in the course of infection from populations of HIV-1-infected individuals who only harbored CCR5-using HIV-1 variants at all time points. Presence or absence of a CCR5 wt/Δ32 genotype and progressive or long-term nonprogressive course of infection stratify the clinical populations in a two-way design. As compared with the standard approach of analyzing sequences from each patient independently, the HPM provides more efficient estimation of evolutionary parameters such as nucleotide substitution rates and d(N)/d(S) rate ratios, as shown by significant shrinkage of the estimator variance. The fixed effects also correct for nonindependence of data between populations and results in even further shrinkage of individual patient estimates. Model selection suggests an association between nucleotide substitution rate and disease progression, but a role for CCR5 genotype remains elusive. Given the absence of clear d(N)/d(S) differences between patient groups, delayed onset of AIDS symptoms appears to be solely associated with lower viral replication rates rather than with differences in selection on amino acid fixation.