Modeling HIV infection dynamics: the role of recombination in the development of drug resistance

Viral kinetics suggests a reconciliation of the disparate observations of the modulation of claudin-1 expression on cells exposed to hepatitis C virus

The tight junction protein claudin-1 (CLDN1) is necessary for hepatitis C virus (HCV) entry into target cells. Recent studies have made disparate observations of the modulation of the expression of CLDN1 on cells following infection by HCV. In one study, the mean CLDN1 expression on cells exposed to HCV declined, whereas in another study HCV infected cells showed increased CLDN1 expression compared to uninfected cells. Consequently, the role of HCV in modulating CLDN1 expression, and hence the frequency of cellular superinfection, remains unclear. Here, we present a possible reconciliation of these disparate observations. We hypothesized that viral kinetics and not necessarily HCV-induced receptor modulation underlies these disparate observations. To test this hypothesis, we constructed a mathematical model of viral kinetics in vitro that mimicked the above experiments. Model predictions provided good fits to the observed evolution of the distribution of CLDN1 expression on cells following exposure to HCV. Cells with higher CLDN1 expression were preferentially infected and outgrown by cells with lower CLDN1 expression, resulting in a decline of the mean CLDN1 expression with time. At the same time, because the susceptibility of cells to infection increased with CLDN1 expression, infected cells tended to have higher CLDN1 expression on average than uninfected cells. Our study thus presents an explanation of the disparate observations of CLDN1 expression following HCV infection and points to the importance of considering viral kinetics in future studies of receptor expression on cells exposed to HCV.

Taking multiple infections of cells and recombination into account leads to small within-host effective-population-size estimates of HIV-1

Whether HIV-1 evolution in infected individuals is dominated by deterministic or stochastic effects remains unclear because current estimates of the effective population size of HIV-1 in vivo, N(e), are widely varying. Models assuming HIV-1 evolution to be neutral estimate N(e)~10²-10⁴, smaller than the inverse mutation rate of HIV-1 (~10⁵), implying the predominance of stochastic forces. In contrast, a model that includes selection estimates N(e)>10⁵, suggesting that deterministic forces would hold sway. The consequent uncertainty in the nature of HIV-1 evolution compromises our ability to describe disease progression and outcomes of therapy. We perform detailed bit-string simulations of viral evolution that consider large genome lengths and incorporate the key evolutionary processes underlying the genomic diversification of HIV-1 in infected individuals, namely, mutation, multiple infections of cells, recombination, selection, and epistatic interactions between multiple loci. Our simulations describe quantitatively the evolution of HIV-1 diversity and divergence in patients. From comparisons of our simulations with patient data, we estimate N(e)~10³-10⁴, implying predominantly stochastic evolution. Interestingly, we find that N(e) and the viral generation time are correlated with the disease progression time, presenting a route to a priori prediction of disease progression in patients. Further, we show that the previous estimate of N(e)>10⁵ reduces as the frequencies of multiple infections of cells and recombination assumed increase. Our simulations with N(e)~10³-10⁴ may be employed to estimate markers of disease progression and outcomes of therapy that depend on the evolution of viral diversity and divergence.

Estimating frequencies of minority nevirapine-resistant strains in chronically HIV-1-infected individuals naive to nevirapine by using stochastic simulations and a mathematical model

Nevirapine forms the mainstay of our efforts to curtail the pediatric AIDS epidemic through prevention of mother-to-child transmission of HIV-1. A key limitation, however, is the rapid selection of HIV-1 strains resistant to nevirapine following the administration of a single dose. This rapid selection of resistance suggests that nevirapine-resistant strains preexist in HIV-1 patients and may adversely affect outcomes of treatment. The frequencies of nevirapine-resistant strains in vivo, however, remain poorly estimated, possibly because they exist as a minority below current assay detection limits. Here, we employ stochastic simulations and a mathematical model to estimate the frequencies of strains carrying different combinations of the common nevirapine resistance mutations K103N, V106A, Y181C, Y188C, and G190A in chronically infected HIV-1 patients naïve to nevirapine. We estimate the relative fitness of mutant strains from an independent analysis of previous competitive growth assays. We predict that single mutants are likely to preexist in patients at frequencies ( approximately 0.01% to 0.001%) near or below current assay detection limits (>0.01%), emphasizing the need for more-sensitive assays. The existence of double mutants is subject to large stochastic variations. Triple and higher mutants are predicted not to exist. Our estimates are robust to variations in the recombination rate, cellular superinfection frequency, and the effective population size. Thus, with 10(7) to 10(8) infected cells in HIV-1 patients, even when undetected, nevirapine-resistant genomes may exist in substantial numbers and compromise efforts to prevent mother-to-child transmission of HIV-1, accelerate the failure of subsequent antiretroviral treatments, and facilitate the transmission of drug resistance.

Timing the emergence of resistance to anti-HIV drugs with large genetic barriers

New antiretroviral drugs that offer large genetic barriers to resistance, such as the recently approved inhibitors of HIV-1 protease, tipranavir and darunavir, present promising weapons to avert the failure of current therapies for HIV infection. Optimal treatment strategies with the new drugs, however, are yet to be established. A key limitation is the poor understanding of the process by which HIV surmounts large genetic barriers to resistance. Extant models of HIV dynamics are predicated on the predominance of deterministic forces underlying the emergence of resistant genomes. In contrast, stochastic forces may dominate, especially when the genetic barrier is large, and delay the emergence of resistant genomes. We develop a mathematical model of HIV dynamics under the influence of an antiretroviral drug to predict the waiting time for the emergence of genomes that carry the requisite mutations to overcome the genetic barrier of the drug. We apply our model to describe the development of resistance to tipranavir in in vitro serial passage experiments. Model predictions of the times of emergence of different mutant genomes with increasing resistance to tipranavir are in quantitative agreement with experiments, indicating that our model captures the dynamics of the development of resistance to antiretroviral drugs accurately. Further, model predictions provide insights into the influence of underlying evolutionary processes such as recombination on the development of resistance, and suggest guidelines for drug design: drugs that offer large genetic barriers to resistance with resistance sites tightly localized on the viral genome and exhibiting positive epistatic interactions maximally inhibit the emergence of resistant genomes.