Modeling and Estimation of Kinetic Parameters and Replicative Fitness of HIV-1 from Flow-Cytometry-Based Growth Competition Experiments

Growth competition assays have been developed to quantify the relative fitness of HIV-1 mutants. In this article, we develop mathematical models to describe viral/cellular dynamic interactions in the assay system from which the competitive fitness indices or parameters are defined. In our previous HIV-viral fitness experiments, the concentration of uninfected target cells was assumed to be constant (Wu et al. 2006). But this may not be true in some experiments. In addition, dual infection may frequently occur in viral fitness experiments and may not be ignorable. Here, we relax these two assumptions and extend our earlier viral fitness model (Wu et al. 2006). The resulting models then become nonlinear ODE systems for which closed-form solutions are not achievable. In the new model, the viral relative fitness is a function of time since it depends on the target cell concentration. First, we studied the structure identifiability of the nonlinear ODE models. The identifiability analysis showed that all parameters in the proposed models are identifiable from the flow-cytometry-based experimental data that we collected. We then employed a global optimization approach (the differential evolution algorithm) to directly estimate the kinetic parameters as well as the relative fitness index in the nonlinear ODE models using nonlinear least square regression based on the experimental data. Practical identifiability was investigated via Monte Carlo simulations.     

Pairwise Growth Competition Assay for Determining the Replication Fitness of Human Immunodeficiency Viruses

In vitro fitness assays are essential tools for determining viral replication fitness for viruses such as HIV-1. Various measurements have been used to extrapolate viral replication fitness, ranging from the number of viral particles per infectious unit, growth rate in cell culture, and relative fitness derived from multiple-cycle growth competition assays. Growth competition assays provide a particularly sensitive measurement of fitness since the viruses are competing for cellular targets under identical growth conditions. There are several experimental factors to consider when conducting growth competition assays, including the multiplicity of infection (MOI), sampling times, and viral detection and fitness calculation methods. Each factor can affect the end result and hence must be considered carefully during the experimental design. The protocol presented here includes steps from constructing a new recombinant HIV-1 clone to performing growth competition assays and analyzing the experimental results. This protocol utilizes experimental parameter values previously shown to yield consistent and robust results. Alternatives are discussed, as some parameters need to be adjusted according to the cell type and viruses being studied. The protocol contains two alternative viral detection methods to provide flexibility as the availability of instruments, reagents and expertise varies between laboratories.     

Evolution of human immunodeficiency virus type 1 cytotoxic T-lymphocyte epitopes: fitness-balanced escape

CD8(+) cytotoxic T lymphocytes (CTL) are strong mediators of human immunodeficiency virus type 1 (HIV-1) control, yet HIV-1 frequently mutates to escape CTL recognition. In an analysis of sequences in the Los Alamos HIV-1 database, we show that emerging CTL escape mutations were more often present at lower frequencies than the amino acid(s) that they replaced. Furthermore, epitopes that underwent escape contained amino acid sites of high variability, whereas epitopes persisting at high frequencies lacked highly variable sites. We therefore infer that escape mutations are likely to be associated with weak functional constraints on the viral protein. This was supported by an extensive analysis of one subject for whom all escape mutations within defined CTL epitopes were studied and by an analysis of all reported escape mutations of defined CTL epitopes in the HIV Immunology Database. In one of these defined epitopes, escape mutations involving the substitution of amino acids with lower database frequencies occurred, and the epitope soon reverted back to the sensitive form. We further show that this escape mutation substantially diminished viral fitness in in vitro competition assays. Coincident with the reversion in vivo, we observed the fixation of a mutation 3 amino acids C terminal to the epitope, coincident with the ablation of the corresponding CTL response. The C-terminal mutation did not restore replication fitness reduced by the escape mutation in the epitope and by itself had little effect on replication fitness. Therefore, this C-terminal mutation presumably impaired the processing and presentation of the epitope. Finally, for one persistent epitope, CTL cross-reactivity to a mutant form may have suppressed the mutant to undetected levels, whereas for two other persistent epitopes, each of two mutants showed poor cross-reactivity and appeared in the subject at later time points. Thus, a viral dynamic exists between the advantage of immune escape, peptide cross-reactivity, and the disadvantage of lost replication fitness, with the balance playing an important role in determining whether a CTL epitope will persist or decline during infection.     

Comparative Fitness of Multi-Dideoxynucleoside-Resistant Human Immunodeficiency Virus Type 1 (HIV-1) in an In Vitro Competitive HIV-1 Replication Assay

We examined whether human immunodeficiency virus type 1 (HIV-1) fitness was altered upon the acquisition of a set or subset of five mutations (A62V, V75I, F77L, F116Y, and Q151M) in the pol gene, which confers resistance to multiple dideoxynucleosides (MDR), as well as the zidovudine resistance-associated mutation T215Y, using a competitive HIV-1 replication assay in a setting of an HXB2D genetic background. Target H9 cells were exposed to a 50:50 mixture of paired infectious molecular clones, and HIV-1 in the culture supernatant was transmitted to new cultures every 7 to 10 days. The polymerase-encoding region of the virus was sequenced at various time points, and the relative proportion of the two viral populations was determined. In the absence of drugs, the comparative order for replicative fitness was HIV-162/75/77/116/151 > HIV-177/116/151 > HIV-1151 > wild-type HIV-1 (HIV-1wt) > HIV-175/77/116/151 > HIV-1151/215 > HIV-1215. In the presence of zidovudine or didanosine, the order was HIV-162/75/77/116/151 > HIV-177/116/151 > HIV-175/77/116/151 > HIV-1151 > HIV-1215. HIV-1215S(TCC), a putative intermediate infectious clone for HIV-1215, replicated comparably to HIV-1wt, while two putative intermediates for HIV-1151 [HIV-1151L(CTG) and HIV-1151K(AAG)] replicated much less efficiently than HIV-1wt and HIV-1151, suggesting that for HIV-1151 to develop, two base substitutions are likely to occur concurrently or within a short interval. These data may illustrate the molecular basis by which HIV-1151 emerges much less frequently than HIV-1215. The present data also demonstrate that several MDR HIV-1 variants are more fit than HIV-1wt in the absence of drugs and that resistance-associated mutations and drug pressure are critical variates for HIV-1 fitness.     

Changes in human immunodeficiency virus type 1 fitness and genetic diversity during disease progression

This study examined the relationship between ex vivo human immunodeficiency virus type 1 (HIV-1) fitness and viral genetic diversity during the course of HIV-1 disease. Primary HIV-1 isolates from 10 patients at different time points were competed against control HIV-1 strains in peripheral blood mononuclear cell (PBMC) cultures to determine relative fitness values. Patient HIV-1 isolates sequentially gained fitness during disease at a significant rate that directly correlated with viral load and HIV-1 env C2V3 diversity. A loss in both fitness and viral diversity was observed upon the initiation of antiretroviral therapy. A possible relationship between genotype and phenotype (virus replication efficiency) is supported by the parallel increases in ex vivo fitness and viral diversity during disease, of which the correlation is largely based on specific V3 sequences. Syncytium-inducing, CXCR4-tropic HIV-1 isolates did have higher relative fitness values than non-syncytium-inducing, CCR5-tropic HIV-1 isolates, as determined by dual virus competitions in PBMC, but increases in fitness during disease were not solely powered by a gradual switch in coreceptor usage. These data provide in vivo evidence that increasing HIV-1 replication efficiency may be related to a concomitant increase in HIV-1 diversity, which in turn may be a determining factor in disease progression.     

A dual infection/competition assay shows a correlation between ex vivo human immunodeficiency virus type 1 fitness and disease progression

This study was designed to examine the impact of human immunodeficiency virus type 1 (HIV-1) fitness on disease progression through the use of a dual competition/heteroduplex tracking assay (HTA). Despite numerous studies on the impact of HIV-1 diversity and HIV-specific immune response on disease progression, we still do not have a firm understanding of the long-term pathogenesis of this virus. Strong and early CD8-positive cytotoxic T-cell and CD4-positive T-helper cell responses directed toward HIV-infected cells appear to curb HIV pathogenesis. However, the rate at which the virus infects the CD4(+) T-cell population and possibly destroys the HIV-specific immune response may also alter the rate of disease progression. For HIV-1 fitness studies, we established conditions for dual HIV-1 infections of peripheral blood mononuclear cells (PBMC) and a sensitive HTA to measure relative virus production. A pairwise comparison was then performed to estimate the relative fitness of various non-syncytium-inducing/CCR5-tropic (NSI/R5) and syncytium-inducing/CXCR4-tropic (SI/X4) HIV-1 isolates. Four HIV-1 strains (two NSI/R5 and two SI/X4) with moderate ex vivo fitness were then selected as controls and competed against primary HIV-1 isolates from an HIV-infected Belgian cohort. HIV-1 isolates from long-term survivors (LTS) were outcompeted by control strains and were significantly less fit than HIV-1 isolates from patients with accelerated progression to AIDS (PRO). In addition, NSI/R5 HIV-1 isolates from PRO overgrew control SI/X4 strains, suggesting that not all SI/X4 HIV-1 isolates replicate more efficiently than all NSI/R5 isolates. Finally, there were strong, independent correlations between viral load and the total relative fitness values of HIV-1 isolates from PRO (r = 0.84, P = 0.033) and LTS (r = 0.86, P = 0.028). Separation of the PRO and LTS plots suggest that HIV-1 fitness together with viral load may be a strong predictor for the rate of disease progression.     

Amino acid substitutions in Gag protein at non-cleavage sites are indispensable for the development of a high multitude of HIV-1 resistance against protease inhibitors.

Amino acid substitutions in human immunodeficiency virus type 1 (HIV-1) Gag cleavage sites have been identified in HIV-1 isolated from patients with AIDS failing chemotherapy containing protease inhibitors (PIs). However, a number of highly PI-resistant HIV-1 variants lack cleavage site amino acid substitutions. In this study we identified multiple novel amino acid substitutions including L75R, H219Q, V390D/V390A, R409K, and E468K in the Gag protein at non-cleavage sites in common among HIV-1 variants selected against the following four PIs: amprenavir, JE-2147, KNI-272, and UIC-94003. Analyses of replication profiles of various mutant clones including competitive HIV-1 replication assays demonstrated that these mutations were indispensable for HIV-1 replication in the presence of PIs. When some of these mutations were reverted to wild type amino acids, such HIV-1 clones failed to replicate. However, virtually the same Gag cleavage pattern was seen, indicating that the mutations affected Gag protein functions but not their cleavage sensitivity to protease. These data strongly suggest that non-cleavage site amino acid substitutions in the Gag protein recover the reduced replicative fitness of HIV-1 caused by mutations in the viral protease and may open a new avenue for designing PIs that resist the emergence of PI-resistant HIV-1.     

Evaluation of the NucliSens EasyQ v2.0 Assay in Comparison with the Roche Amplicor v1.5 and the Roche CAP/CTM HIV-1 Test v2.0 in Quantification of C-Clade HIV-1 in Plasma

Human immunodeficiency virus type 1 (HIV-1) genetic diversity poses a challenge to reliable viral load monitoring. Discrepancies between different testing platforms have been observed, especially for non-clade-B virus. Therefore we compare, in antiretroviral therapy (ART)-naïve South African subjects predominantly infected with HIV-1 clade-C, three commercially available assays: the COBAS AmpliPrep/COBAS TaqMan HIV-1 Test version 2.0 by Roche (CAP/CTM v2.0), the BioMérieux NucliSens Version 2.0 Easy Q/Easy Mag (NucliSens v2.0) and the Roche COBAS Amplicor HIV-1 Monitor Test Version 1.5 (Amplicor v1.5). Strong linear correlation was observed and Bland-Altman analyses showed overall good agreement between the assays with mean viral load differences of 0.078 log cp/ml (NucliSens v2.0 – Amplicor v1.5), 0.260 log cp/ml (CAP/CTM v2.0 – Amplicor v1.5) and 0.164 log cp/ml (CAP/CTM v2.0 – NucliSens v2.0), indicating lower mean viral load results for the Amplicor v1.5 and higher mean readings for the CAP/CTM v2.0. Consistent with observations following previous comparisons of CAP/CTM v2.0 versus Amplicor v1.5, the CAP/CTM v2.0 assay detected low-level viremia (median 65 cp/ml) in more than one-third of those in whom viremia had been undetectable (<20 cp/ml) in assays using the NucliSens platform. These levels of viremia are of uncertain clinical significance but may be of importance in early detection of ART resistance in those on treatment. Overall the three assays showed good comparability of results but with consistent, albeit relatively small, discrepancies for HIV-1 clade-C samples, especially in the low-viremic range that should be taken into account when interpreting viral load data.     

Estimating relative fitness in viral competition experiments

The relative fitness of viral variants has previously been defined as the slope of the logarithmic ratio of the genotype or phenotype frequencies in time plots of pairwise competition experiments. Developing mathematical models for such experiments by employing the conventional coefficient of selection s, we demonstrate that this logarithmic ratio gives the fitness difference, rather than the relative fitness. This fitness difference remains proportional to the actual replication rate realized in the particular experimental setup and hence cannot be extrapolated to other situations. Conversely, the conventional relative fitness (1 + s) should be more generic. We develop an approach to compute the generic relative fitness in conventional competition experiments. This involves an estimation of the total viral replication during the experiment and requires an estimate of the average lifetime of productively infected cells. The novel approach is illustrated by estimating the relative fitness, i.e., the relative replication rate, of a set of zidovudine-resistant human immunodeficiency virus type 1 variants. A tool for calculating the relative fitness from observed changes in viral load and genotype (or phenotype) frequencies is publically available on the website at http://www-binf.bio.uu.nl/( approximately )rdb/fitness.html.     

Decreased infectivity of a neutralization-resistant equine infectious anemia virus variant can be overcome by efficient cell-to-cell spread

Two variants of equine infectious anemia virus (EIAV) that differed in sensitivity to broadly neutralizing antibody were tested in direct competition assays. No differences were observed in the growth curves and relative fitness scores of EIAVs of principal neutralizing domain variants of groups 1 (EIAV(PND-1)) and 5 (EIAV(PND-5)), respectively; however, the neutralization-resistant EIAV(PND-5) variant was less infectious in single-round replication assays. Infectious center assays indicated similar rates of cell-to-cell spread, which was approximately 1,000-fold more efficient than cell-free infectivity. These data indicate that efficient cell-to-cell spread can overcome the decreased infectivity that may accompany immune escape and should be considered in studies assessing the relative levels of fitness among lentivirus variants, including HIV-1.     

High Viral Fitness during Acute HIV-1 Infection

Several clinical studies have shown that, relative to disease progression, HIV-1 isolates that are less fit are also less pathogenic. The aim of the present study was to investigate the relationship between viral fitness and control of viral load (VL) in acute and early HIV-1 infection. Samples were obtained from subjects participating in two clinical studies. In the PULSE study, antiretroviral therapy (ART) was initiated before, or no later than six months following seroconversion. Subjects then underwent multiple structured treatment interruptions (STIs). The PHAEDRA study enrolled and monitored a cohort of individuals with documented evidence of primary infection. The subset chosen were individuals identified no later than 12 months following seroconversion to HIV-1, who were not receiving ART. The relative fitness of primary isolates obtained from study participants was investigated ex vivo. Viral DNA production was quantified using a novel real time PCR assay. Following intermittent ART, the fitness of isolates obtained from 5 of 6 PULSE subjects decreased over time. In contrast, in the absence of ART the fitness of paired isolates obtained from 7 of 9 PHAEDRA subjects increased over time. However, viral fitness did not correlate with plasma VL. Most unexpected was the high relative fitness of isolates obtained at Baseline from PULSE subjects, before initiating ART. It is widely thought that the fitness of strains present during the acute phase is low relative to strains present during chronic HIV-1 infection, due to the bottleneck imposed upon transmission. The results of this study provide evidence that the relative fitness of strains present during acute HIV-1 infection may be higher than previously thought. Furthermore, that viral fitness may represent an important clinical parameter to be considered when deciding whether to initiate ART during early HIV-1 infection.     

Infection with multiple human immunodeficiency virus type 1 variants is associated with faster disease progression

Human immunodeficiency virus type 1 (HIV-1)-infected individuals develop a genetically diverse virus population over time, but often only a limited number of viral variants are transmitted from a chronic carrier to a newly infected person. Interestingly, many women but few men are infected by multiple HIV-1 variants from a single partner. To determine whether the complexity of the infecting virus population influences clinical outcome, we examined viral diversity in the HIV-1 envelope sequences present at primary infection in 156 women from Kenya for whom we had follow-up data on viral RNA levels and CD4 T-cell counts. Eighty-nine women had multiple viral genotypes, while 67 women had a single genotype at primary infection. Women who acquired multiple viral genotypes had a significantly higher viral load (median, 4.84 versus 4.64 log(10) copies/ml, P = 0.04) and a significantly lower CD4(+)-T-cell count (median, 416 versus 617 cells/mm(3), P = 0.01) 4 to 24 months after infection compared to women who were infected with a single viral genotype. These studies suggest that early HIV-1 genetic diversity is linked to faster disease progression.     

Fitness ranking of individual mutants drives patterns of epistatic interactions in HIV-1

Fitness interactions between mutations, referred to as epistasis, can strongly impact evolution. For RNA viruses and retroviruses with their high mutation rates, epistasis may be particularly important to overcome fitness losses due to the accumulation of deleterious mutations and thus could influence the frequency of mutants in a viral population. As human immunodeficiency virus type 1 (HIV-1) resistance to azidothymidine (AZT) requires selection of sequential mutations, it is a good system to study the impact of epistasis. Here we present a thorough analysis of a classical AZT-resistance pathway (the 41-215 cluster) of HIV-1 variants by fitness measurements in single round infection assays covering physiological drug concentrations ex vivo. The sign and value of epistasis varied and did not predict the epistatic effect on the mutant frequency. This complex behavior is explained by the fitness ranking of the variants that strongly depends on environmental factors, i.e., the presence and absence of drugs and the host cells used. Although some interactions compensate fitness losses, the observed small effect on the relative mutant frequencies suggests that epistasis might be inefficient as a buffering mechanism for fitness losses in vivo. While the use of epistasis-based hypotheses to make general assumptions on the evolutionary dynamics of viral populations is appealing, our data caution their interpretation without further knowledge on the characteristics of the viral mutant spectrum under different environmental conditions.     

Role of the human immunodeficiency virus type 1 envelope gene in viral fitness

A human host offers a variety of microenvironments to the infecting human immunodeficiency virus type 1 (HIV-1), resulting in various selective pressures, most of them directed against the envelope (env) gene. Therefore, it seems evident that the replicative capacity of the virus is largely related to viral entry. In this study we have used growth competition experiments and TaqMan real-time PCR detection to measure the fitness of subtype B HIV-1 primary isolates and autologous env-recombinant viruses in order to analyze the contribution of wild-type env sequences to overall HIV-1 fitness. A significant correlation was observed between fitness values obtained for wild-type HIV-1 isolates and those for the corresponding env-recombinant viruses (r = 0.93; P = 0.002). Our results suggest that the env gene, which is linked to a myriad of viral characteristics (e.g., entry into the host cell, transmission, coreceptor usage, and tropism), plays a major role in fitness of wild-type HIV-1. In addition, this new recombinant assay may be useful for measuring the contribution of HIV-1 env to fitness in viruses resistant to novel antiretroviral entry inhibitors.     

Potent activity of the HIV-1 maturation inhibitor bevirimat in SCID-hu Thy/Liv mice

Background

The HIV-1 maturation inhibitor, 3-O-(3′,3′-dimethylsuccinyl) betulinic acid (bevirimat, PA-457) is a promising drug candidate with 10 nM in vitro antiviral activity against multiple wild-type (WT) and drug-resistant HIV-1 isolates. Bevirimat has a novel mechanism of action, specifically inhibiting cleavage of spacer peptide 1 (SP1) from the C-terminus of capsid which results in defective core condensation.

Methods and Findings

Oral administration of bevirimat to HIV-1-infected SCID-hu Thy/Liv mice reduced viral RNA by >2 log₁₀ and protected immature and mature T cells from virus-mediated depletion. This activity was observed at plasma concentrations that are achievable in humans after oral dosing, and bevirimat was active up to 3 days after inoculation with both WT HIV-1 and an AZT-resistant HIV-1 clinical isolate. Consistent with its mechanism of action, bevirimat caused a dose-dependent inhibition of capsid-SP1 cleavage in HIV-1-infected human thymocytes obtained from these mice. HIV-1 NL4-3 with an alanine-to-valine substitution at the N-terminus of SP1 (SP1/A1V), which is resistant to bevirimat in vitro, was also resistant to bevirimat treatment in the mice, and SP1/AIV had replication and thymocyte kinetics similar to that of WT NL4-3 with no evidence of fitness impairment in in vivo competition assays. Interestingly, protease inhibitor-resistant HIV-1 with impaired capsid-SP1 cleavage was hypersensitive to bevirimat in vitro with a 50% inhibitory concentration 140 times lower than for WT HIV-1.

Conclusions

These results support further clinical development of this first-in-class maturation inhibitor and confirm the usefulness of the SCID-hu Thy/Liv model for evaluation of in vivo antiretroviral efficacy, drug resistance, and viral fitness.     

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.