Search

Motivated by viral persistence in HIV+ patients on long-term anti-retroviral treatment (ART), we present a stochastic model of HIV viral dynamics in the blood stream. We consider the hypothesis that the residual viremia in patients on ART can be explained principally by the activation of cells latently infected by HIV before the initiation of ART and that viral blips (clinically-observed short periods of detectable viral load) represent large deviations from the mean. We model the system as a continuous-time, multi-type branching process. Deriving equations for the probability generating function we use a novel numerical approach to extract the probability distributions for latent reservoir sizes and viral loads. We find that latent reservoir extinction-time distributions underscore the importance of considering reservoir dynamics beyond simply the half-life. We calculate blip amplitudes and frequencies by computing complete viral load probability distributions, and study the duration of viral blips via direct numerical simulation. We find that our model qualitatively reproduces short small-amplitude blips detected in clinical studies of treated HIV infection. Stochastic models of this type provide insight into treatment-outcome variability that cannot be found from deterministic models.

Most people living with HIV-1 experience rapid viral rebound once antiretroviral therapy is interrupted; however, a small fraction remain in viral remission for an extended duration. Understanding the factors that determine whether viral rebound is likely after treatment interruption can enable the development of optimal treatment regimens and therapeutic interventions to potentially achieve a functional cure for HIV-1. We built upon the theoretical framework proposed by Conway and Perelson to construct dynamic models of virus-immune interactions to study factors that influence viral rebound dynamics. We evaluated these models using viral load data from 24 individuals following antiretroviral therapy interruption. The best-performing model accurately captures the heterogeneity of viral dynamics and highlights the importance of the effector cell expansion rate. Our results show that post-treatment controllers and non-controllers can be distinguished based on the effector cell expansion rate in our models. Furthermore, these results demonstrate the potential of using dynamic models incorporating an effector cell response to understand early viral rebound dynamics post-antiretroviral therapy interruption. Author summary: Most people living with HIV-1 experience rapid viral rebound once antiretroviral therapy is interrupted; however, a small fraction remain in viral remission for an extended duration. The factors that determine viral rebound dynamics after treatment interruption are not well understood. In this study, we built upon a previous theoretical framework to construct dynamic models of virus-immune interactions to study factors that influence viral rebound dynamics. We evaluated these models using viral load data from 24 individuals following antiretroviral therapy interruption. The best-performing model accurately captures the heterogeneity of viral dynamics and robustly shows that people who remain in viral remission for an extended duration after treatment interruption have a higher effector cell expansion rate compared to those who rapidly rebound. These results demonstrate the potential of using viral dynamic models incorporating an effector cell response to understand early viral rebound dynamics post-antiretroviral therapy interruption. [ABSTRACT FROM AUTHOR]

Antiretroviral therapy (ART) effectively controls HIV infection, suppressing HIV viral loads. Suspension of therapy is followed by rebound of viral loads to high, pre-therapy levels. However, there is significant heterogeneity in speed of rebound, with some rebounds occurring within days, weeks, or sometimes years. We present a stochastic mathematical model to gain insight into these post-treatment dynamics, specifically characterizing the dynamics of short term viral rebounds (≤ 60 days). Li et al. (2016) report that the size of the expressed HIV reservoir, i.e., cell-associated HIV RNA levels, and drug regimen correlate with the time between ART suspension and viral rebound to detectable levels. We incorporate this information and viral rebound times to parametrize our model. We then investigate insights offered by our model into the underlying dynamics of the latent reservoir. In particular, we refine previous estimates of viral recrudescence after ART interruption by accounting for heterogeneity in infection rebound dynamics, and determine a recrudescence rate of once every 2-4 days. Our parametrized model can be used to aid in design of clinical trials to study viral dynamics following analytic treatment interruption. We show how to derive informative personalized testing frequencies from our model and offer a proof-of-concept example. Our results represent first steps towards a model that can make predictions on a person living with HIV (PLWH)’s rebound time distribution based on biomarkers, and help identify PLWH with long viral rebound delays. [ABSTRACT FROM AUTHOR]

In order to prevent and/or control infections it is necessary to understand their early-time dynamics. However, this is precisely the phase of HIV about which the least is known. To investigate the initial stages of HIV infection within a host we have developed a multitype continuoustime branching process model. This model is a stochastic extension of the standard viral dynamics model, under the assumption that the number of cell targets for viral infection is constant. We use our model to investigate three important clinical characteristics of early HIV infection following intravenous challenge: risk of infection, time to infection clearance (assuming failed infection), and time to infection detection. Our focus is on the impact of errors in viral replication that result in noninfectious virus production on these characteristics. Only a small fraction of circulating virus in any chronically infected individual is capable of infecting susceptible cells: estimates range from 1/104 to 1/103. Characterization and quantification of the processes by which virus becomes defective remains incomplete. We consider two mechanisms that result in defective virus: (1) Copying errors, i.e., lethal errors in reverse transcription, which introduce mutations into the HIV-1 proviral genome, some of which may cripple the viral genome produced, and (2) Packaging errors, i.e., errors during viral packaging, at the end of the viral replication cycle, which cause a defective virus by packaging new virions without, for example, viral RNA or key proteins required for infectivity. We show that assumptions on mechanisms of defective virus production can significantly impact early HIV infection model predictions. For example, the risk of infection is orders of magnitude higher if all defective virus is associated with packaging errors, but infection is predicted to be detectable sooner following HIV exposure if all defective virus is associated with copying errors. Thus, in order to make reliable predictions of risk, clearance time, and detection time, better characterization of viral replication is required. [ABSTRACT FROM AUTHOR]

HIV infection cannot be detected immediately after exposure because plasma viral loads are too small initially. The duration of this phase of infection (the “eclipse period”) is difficult to estimate because precise dates of exposure are rarely known. Therefore, the reliability of clinical HIV testing during the first few weeks of infection is unknown, creating anxiety among HIV-exposed individuals and their physicians. We address this by fitting stochastic models of early HIV infection to detailed viral load records for 78 plasma donors, taken during the period of exposure and infection. We first show that the classic stochastic birth-death model does not satisfactorily describe early infection. We therefore apply a different stochastic model that includes infected cells and virions separately. Since every plasma donor in our data eventually becomes infected, we must condition the model to reflect this bias, before fitting to the data. Applying our best estimates of unknown parameter values, we estimate the mean eclipse period to be 8–10 days. We further estimate the reliability of a negative test t days after potential exposure. [ABSTRACT FROM AUTHOR]

HIV-1 is subject to immune pressure exerted by the host, giving variants that escape the immune response an advantage. Virus released from activated latent cells competes against variants that have continually evolved and adapted to host immune pressure. Nevertheless, there is increasing evidence that virus displaying a signal of latency survives in patient plasma despite having reduced fitness due to long-term immune memory. We investigated the survival of virus with latent envelope genomic fragments by simulating within-host HIV-1 sequence evolution and the cycling of viral lineages in and out of the latent reservoir. Our model incorporates a detailed mutation process including nucleotide substitution, recombination, latent reservoir dynamics, diversifying selection pressure driven by the immune response, and purifying selection pressure asserted by deleterious mutations. We evaluated the ability of our model to capture sequence evolution in vivo by comparing our simulated sequences to HIV-1 envelope sequence data from 16 HIV-infected untreated patients. Empirical sequence divergence and diversity measures were qualitatively and quantitatively similar to those of our simulated HIV-1 populations, suggesting that our model invokes realistic trends of HIV-1 genetic evolution. Moreover, reconstructed phylogenies of simulated and patient HIV-1 populations showed similar topological structures. Our simulation results suggest that recombination is a key mechanism facilitating the persistence of virus with latent envelope genomic fragments in the productively infected cell population. Recombination increased the survival probability of latent virus forms approximately 13-fold. Prevalence of virus with latent fragments in productively infected cells was observed in only 2% of simulations when we ignored recombination, while the proportion increased to 27% of simulations when we allowed recombination. We also found that the selection pressures exerted by different fitness landscapes influenced the shape of phylogenies, diversity trends, and survival of virus with latent genomic fragments. Our model predicts that the persistence of latent genomic fragments from multiple different ancestral origins increases sequence diversity in plasma for reasonable fitness landscapes. [ABSTRACT FROM AUTHOR]