Modeling suggests SARS-CoV-2 rebound after nirmatrelvir-ritonavir treatment is driven by target cell preservation coupled with incomplete viral clearance.

In a subset of SARS-CoV-2 infected individuals treated with the oral antiviral nirmatrelvir-ritonavir, the virus rebounds following treatment. The mechanisms driving this rebound are not well understood. We used a mathematical model to describe the longitudinal viral load dynamics of 51 individuals treated with nirmatrelvir-ritonavir, 20 of whom rebounded. Target cell preservation, either by a robust innate immune response or initiation of nirmatrelvir-ritonavir near the time of symptom onset, coupled with incomplete viral clearance, appear to be the main factors leading to viral rebound. Moreover, the occurrence of viral rebound is likely influenced by time of treatment initiation relative to the progression of the infection, with earlier treatments leading to a higher chance of rebound. Finally, our model demonstrates that extending the course of nirmatrelvir-ritonavir treatment, in particular to a 10-day regimen, may greatly diminish the risk for rebound in people with mild-to-moderate COVID-19 and who are at high risk of progression to severe disease. Altogether, our results suggest that in some individuals, a standard 5-day course of nirmatrelvir-ritonavir starting around the time of symptom onset may not completely eliminate the virus. Thus, after treatment ends, the virus can rebound if an effective adaptive immune response has not fully developed. These findings on the role of target cell preservation and incomplete viral clearance also offer a possible explanation for viral rebounds following other antiviral treatments for SARS-CoV-2.

Understanding early HIV-1 rebound dynamics following antiretroviral therapy interruption: The importance of effector cell expansion.

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.

A multiscale model of the action of a capsid assembly modulator for the treatment of chronic hepatitis B.

Chronic hepatitis B virus (HBV) infection is strongly associated with increased risk of liver cancer and cirrhosis. While existing treatments effectively inhibit the HBV life cycle, viral rebound occurs rapidly following treatment interruption. Consequently, functional cure rates of chronic HBV infection remain low and there is increased interest in a novel treatment modality, capsid assembly modulators (CAMs). Here, we develop a multiscale mathematical model of CAM treatment in chronic HBV infection. By fitting the model to participant data from a phase I trial of the first-generation CAM vebicorvir, we estimate the drug's dose-dependent effectiveness and identify the physiological mechanisms that drive the observed biphasic decline in HBV DNA and RNA, and mechanistic differences between HBeAg-positive and negative infection. Finally, we demonstrate analytically and numerically that HBV RNA is more sensitive than HBV DNA to increases in CAM effectiveness.

Understanding early HIV-1 rebound dynamics following antiretroviral therapy interruption: The importance of effector cell expansion.

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.

Modeling the emergence of viral resistance for SARS-CoV-2 during treatment with an anti-spike monoclonal antibody.

To mitigate the loss of lives during the COVID-19 pandemic, emergency use authorization was given to several anti-SARS-CoV-2 monoclonal antibody (mAb) therapies for the treatment of mild-to-moderate COVID-19 in patients with a high risk of progressing to severe disease. Monoclonal antibodies used to treat SARS-CoV-2 target the spike protein of the virus and block its ability to enter and infect target cells. Monoclonal antibody therapy can thus accelerate the decline in viral load and lower hospitalization rates among high-risk patients with variants susceptible to mAb therapy. However, viral resistance has been observed, in some cases leading to a transient viral rebound that can be as large as 3-4 orders of magnitude. As mAbs represent a proven treatment choice for SARS-CoV-2 and other viral infections, evaluation of treatment-emergent mAb resistance can help uncover underlying pathobiology of SARS-CoV-2 infection and may also help in the development of the next generation of mAb therapies. Although resistance can be expected, the large rebounds observed are much more difficult to explain. We hypothesize replenishment of target cells is necessary to generate the high transient viral rebound. Thus, we formulated two models with different mechanisms for target cell replenishment (homeostatic proliferation and return from an innate immune response antiviral state) and fit them to data from persons with SARS-CoV-2 treated with a mAb. We showed that both models can explain the emergence of resistant virus associated with high transient viral rebounds. We found that variations in the target cell supply rate and adaptive immunity parameters have a strong impact on the magnitude or observability of the viral rebound associated with the emergence of resistant virus. Both variations in target cell supply rate and adaptive immunity parameters may explain why only some individuals develop observable transient resistant viral rebound. Our study highlights the conditions that can lead to resistance and subsequent viral rebound in mAb treatments during acute infection.

How robust are estimates of key parameters in standard viral dynamic models?

Mathematical models of viral infection have been developed, fitted to data, and provide insight into disease pathogenesis for multiple agents that cause chronic infection, including HIV, hepatitis C, and B virus. However, for agents that cause acute infections or during the acute stage of agents that cause chronic infections, viral load data are often collected after symptoms develop, usually around or after the peak viral load. Consequently, we frequently lack data in the initial phase of viral growth, i.e., when pre-symptomatic transmission events occur. Missing data may make estimating the time of infection, the infectious period, and parameters in viral dynamic models, such as the cell infection rate, difficult. However, having extra information, such as the average time to peak viral load, may improve the robustness of the estimation. Here, we evaluated the robustness of estimates of key model parameters when viral load data prior to the viral load peak is missing, when we know the values of some parameters and/or the time from infection to peak viral load. Although estimates of the time of infection are sensitive to the quality and amount of available data, particularly pre-peak, other parameters important in understanding disease pathogenesis, such as the loss rate of infected cells, are less sensitive. Viral infectivity and the viral production rate are key parameters affecting the robustness of data fits. Fixing their values to literature values can help estimate the remaining model parameters when pre-peak data is missing or limited. We find a lack of data in the pre-peak growth phase underestimates the time to peak viral load by several days, leading to a shorter predicted growth phase. On the other hand, knowing the time of infection (e.g., from epidemiological data) and fixing it results in good estimates of dynamical parameters even in the absence of early data. While we provide ways to approximate model parameters in the absence of early viral load data, our results also suggest that these data, when available, are needed to estimate model parameters more precisely.

Is the exquisite specificity of lymphocytes generated by thymic selection or due to evolution?

We have previously argued that the antigen receptors of T and B lymphocytes evolved to be sufficiently specific to avoid massive deletion of clonotypes by negative selection. Their optimal 'specificity' level, i.e., probability of binding any particular epitope, was shown to be inversely related to the number of self-antigens that the cells have to be tolerant to. Experiments have demonstrated that T lymphocytes also become more specific during negative selection in the thymus, because cells expressing the most crossreactive receptors have the highest likelihood of binding a self-antigen, and hence to be tolerized (i.e., deleted, anergized, or diverted into a regulatory T cell phenotype). Thus, there are two -not mutually exclusive- explanations for the exquisite specificity of T cells, one involving evolution and the other thymic selection. To better understand the impact of both, we extend a previously developed mathematical model by allowing for T cells with very different binding probabilities in the pre-selection repertoire. We confirm that negative selection tends to tolerize the most crossreactive clonotypes. As a result, the average level of specificity in the functional post-selection repertoire depends on the number of self-antigens, even if there is no evolutionary optimization of binding probabilities. However, the evolutionary optimal range of binding probabilities in the pre-selection repertoire also depends on the number of self-antigens. Species with more self antigens need more specific pre-selection repertoires to avoid excessive loss of T cells during thymic selection, and hence mount protective immune responses. We conclude that both evolution and negative selection are responsible for the high level of specificity of lymphocytes.

Multiscale modeling of HBV infection integrating intra- and intercellular viral propagation to analyze extracellular viral markers.

Chronic infection with hepatitis B virus (HBV) is caused by the persistence of closed circular DNA (cccDNA) in the nucleus of infected hepatocytes. Despite available therapeutic anti-HBV agents, eliminating the cccDNA remains challenging. Thus, quantifying and understanding the dynamics of cccDNA are essential for developing effective treatment strategies and new drugs. However, such study requires repeated liver biopsy to measure the intrahepatic cccDNA, which is basically not accepted because liver biopsy is potentially morbid and not common during hepatitis B treatment. We here aimed to develop a noninvasive method for quantifying cccDNA in the liver using surrogate markers in peripheral blood. We constructed a multiscale mathematical model that explicitly incorporates both intracellular and intercellular HBV infection processes. The model, based on age-structured partial differential equations, integrates experimental data from in vitro and in vivo investigations. By applying this model, we roughly predicted the amount and dynamics of intrahepatic cccDNA within a certain range using specific viral markers in serum samples, including HBV DNA, HBsAg, HBeAg, and HBcrAg. Our study represents a significant step towards advancing the understanding of chronic HBV infection. The noninvasive quantification of cccDNA using our proposed method holds promise for improving clinical analyses and treatment strategies. By comprehensively describing the interactions of all components involved in HBV infection, our multiscale mathematical model provides a valuable framework for further research and the development of targeted interventions.

Modeling resistance to the broadly neutralizing antibody PGT121 in people living with HIV-1.

PGT121 is a broadly neutralizing antibody in clinical development for the treatment and prevention of HIV-1 infection via passive administration. PGT121 targets the HIV-1 V3-glycan and demonstrated potent antiviral activity in a phase I clinical trial. Resistance to PGT121 monotherapy rapidly occurred in the majority of participants in this trial with the sampled rebound viruses being entirely resistant to PGT121 mediated neutralization. However, two individuals experienced long-term ART-free viral suppression following antibody infusion and retained sensitivity to PGT121 upon viral rebound. Here, we develop mathematical models of the HIV-1 dynamics during this phase I clinical trial. We utilize these models to understand the dynamics leading to PGT121 resistance and to identify the mechanisms driving the observed long-term viral control. Our modeling highlights the importance of the relative fitness difference between PGT121 sensitive and resistant subpopulations prior to treatment. Specifically, by fitting our models to data, we identify the treatment-induced competitive advantage of previously existing or newly generated resistant population as a primary driver of resistance. Finally, our modeling emphasizes the high neutralization ability of PGT121 in both participants who exhibited long-term viral control.

Induction of durable remission by dual immunotherapy in SHIV-infected ART-suppressed macaques.

The eradication of the viral reservoir represents the major obstacle to the development of a clinical cure for established HIV-1 infection. Here, we demonstrate that the administration of N-803 (brand name Anktiva) and broadly neutralizing antibodies (bNAbs) results in sustained viral control after discontinuation of antiretroviral therapy (ART) in simian-human AD8 (SHIV-AD8)-infected, ART-suppressed rhesus macaques. N-803+bNAbs treatment induced immune activation and transient viremia but only limited reductions in the SHIV reservoir. Upon ART discontinuation, viral rebound occurred in all animals, which was followed by durable control in approximately 70% of all N-803+bNAb-treated macaques. Viral control was correlated with the reprogramming of CD8+ T cells by N-803+bNAb synergy. Thus, complete eradication of the replication-competent viral reservoir is likely not a prerequisite for the induction of sustained remission after discontinuation of ART.

The reproduction number and its probability distribution for stochastic viral dynamics.

We consider stochastic models of individual infected cells. The reproduction number, R, is understood as a random variable representing the number of new cells infected by one initial infected cell in an otherwise susceptible (target cell) population. Variability in R results partly from heterogeneity in the viral burst size (the number of viral progeny generated from an infected cell during its lifetime), which depends on the distribution of cellular lifetimes and on the mechanism of virion release. We analyse viral dynamics models with an eclipse phase: the period of time after a cell is infected but before it is capable of releasing virions. The duration of the eclipse, or the subsequent infectious, phase is non-exponential, but composed of stages. We derive the probability distribution of the reproduction number for these viral dynamics models, and show it is a negative binomial distribution in the case of constant viral release from infectious cells, and under the assumption of an excess of target cells. In a deterministic model, the ultimate in-host establishment or extinction of the viral infection depends entirely on whether the mean reproduction number is greater than, or less than, one, respectively. Here, the probability of extinction is determined by the probability distribution of R, not simply its mean value. In particular, we show that in some cases the probability of infection is not an increasing function of the mean reproduction number.

Last in first out: SIV proviruses seeded later in infection are harbored in short-lived CD4+ T cells.

HIV can persist in a latent form as integrated DNA (provirus) in resting CD4+ T cells of infected individuals and as such is unaffected by antiretroviral therapy (ART). Despite being a major obstacle for eradication efforts, the genetic variation and timing of formation of this latent reservoir remains poorly understood. Previous studies on when virus is deposited in the latent reservoir have come to contradictory conclusions. To reexamine the genetic variation of HIV in CD4+ T cells during ART, we determined the divergence in envelope sequences collected from 10 SIV infected rhesus macaques. We found that the macaques displayed a biphasic decline of the viral divergence over time, where the first phase lasted for an average of 11.6 weeks (range 4-28 weeks). Motivated by recent observations that the HIV-infected CD4+ T cell population is composed of short- and long-lived subsets, we developed a model to study the divergence dynamics. We found that SIV in short-lived cells was on average more diverged, while long-lived cells harbored less diverged virus. This suggests that the long-lived cells harbor virus deposited starting earlier in infection and continuing throughout infection, while short-lived cells predominantly harbor more recent virus. As these cell populations decayed, the overall proviral divergence decline matched that observed in the empirical data. This model explains previous seemingly contradictory results on the timing of virus deposition into the latent reservoir, and should provide guidance for future eradication efforts.

Biphasic decay of intact SHIV genomes following initiation of antiretroviral therapy complicates analysis of interventions targeting the reservoir.

The latent reservoir for HIV-1 in resting CD4+ T cells persists despite antiretroviral therapy (ART) and precludes cure. Reservoir-targeting interventions are evaluated in ART-treated macaques infected with simian immunodeficiency virus (SIV) or simian-human immunodeficiency virus (SHIV). Efficacy is determined by reservoir measurements before and after the intervention. However, most proviruses persisting in the setting of ART are defective. In addition, intact HIV-1 and SIV genomes undergo complex, multiphasic decay observable when new infection events are blocked by ART. Intervention-induced elimination of latently infected cells must be distinguished from natural decay. Here, we address these issues for SHIV. We describe an intact proviral DNA assay that allows digital counting of SHIV genomes lacking common fatal defects. We show that intact SHIV genomes in circulating CD4+ T cells undergo biphasic decay during the first year of ART, with a rapid first phase (t1/2 = 30.1 d) and a slower second phase (t1/2 = 8.1 mo) that is still more rapid that the slow decay observed in people with HIV-1 on long-term ART (t1/2 = 3.7 y). In SHIV models, most interventions are tested during 2nd phase decay. Natural 2nd phase decay must be considered in evaluating interventions as most infected cells present at this time do not become part of the stable reservoir. In addition, for interventions tested during 2nd phase decay, a caveat is that the intervention may not be equally effective in people with HIV on long-term ART whose reservoirs are dominated by latently infected cells with a slower decay rate.

CD8+ T cells control SIV infection using both cytolytic effects and non-cytolytic suppression of virus production.

Whether CD8+ T lymphocytes control human immunodeficiency virus infection by cytopathic or non-cytopathic mechanisms is not fully understood. Multiple studies highlighted non-cytopathic effects, but one hypothesis is that cytopathic effects of CD8+ T cells occur before viral production. Here, to examine the role of CD8+ T cells prior to virus production, we treated SIVmac251-infected macaques with an integrase inhibitor combined with a CD8-depleting antibody, or with either reagent alone. We analyzed the ensuing viral dynamics using a mathematical model that included infected cells pre- and post- viral DNA integration to compare different immune effector mechanisms. Macaques receiving the integrase inhibitor alone experienced greater viral load decays, reaching lower nadirs on treatment, than those treated also with the CD8-depleting antibody. Models including CD8+ cell-mediated reduction of viral production (non-cytolytic) were found to best explain the viral profiles across all macaques, in addition an effect in killing infected cells pre-integration (cytolytic) was supported in some of the best models. Our results suggest that CD8+ T cells have both a cytolytic effect on infected cells before viral integration, and a direct, non-cytolytic effect by suppressing viral production.

AZD5582 plus SIV-specific antibodies reduce lymph node viral reservoirs in antiretroviral therapy-suppressed macaques.

The main barrier to HIV cure is a persistent reservoir of latently infected CD4+ T cells harboring replication-competent provirus that fuels rebound viremia upon antiretroviral therapy (ART) interruption. A leading approach to target this reservoir involves agents that reactivate latent HIV proviruses followed by direct clearance of cells expressing induced viral antigens by immune effector cells and immunotherapeutics. We previously showed that AZD5582, an antagonist of inhibitor of apoptosis proteins and mimetic of the second mitochondrial-derived activator of caspases (IAPi/SMACm), induces systemic reversal of HIV/SIV latency but with no reduction in size of the viral reservoir. In this study, we investigated the effects of AZD5582 in combination with four SIV Env-specific Rhesus monoclonal antibodies (RhmAbs) ± N-803 (an IL-15 superagonist) in SIV-infected, ART-suppressed rhesus macaques. Here we confirm the efficacy of AZD5582 in inducing SIV reactivation, demonstrate enhancement of latency reversal when AZD5582 is used in combination with N-803 and show a reduction in total and replication-competent SIV-DNA in lymph-node-derived CD4+ T cells in macaques treated with AZD5582 + RhmAbs. Further exploration of this therapeutic approach may contribute to the goal of achieving an HIV cure.

Modeling the emergence of viral resistance for SARS-CoV-2 during treatment with an anti-spike monoclonal antibody.

The COVID-19 pandemic has led to over 760 million cases and 6.9 million deaths worldwide. To mitigate the loss of lives, emergency use authorization was given to several anti-SARS-CoV-2 monoclonal antibody (mAb) therapies for the treatment of mild-to-moderate COVID-19 in patients with a high risk of progressing to severe disease. Monoclonal antibodies used to treat SARS-CoV-2 target the spike protein of the virus and block its ability to enter and infect target cells. Monoclonal antibody therapy can thus accelerate the decline in viral load and lower hospitalization rates among high-risk patients with susceptible variants. However, viral resistance has been observed, in some cases leading to a transient viral rebound that can be as large as 3-4 orders of magnitude. As mAbs represent a proven treatment choice for SARS-CoV-2 and other viral infections, evaluation of treatment-emergent mAb resistance can help uncover underlying pathobiology of SARS-CoV-2 infection and may also help in the development of the next generation of mAb therapies. Although resistance can be expected, the large rebounds observed are much more difficult to explain. We hypothesize replenishment of target cells is necessary to generate the high transient viral rebound. Thus, we formulated two models with different mechanisms for target cell replenishment (homeostatic proliferation and return from an innate immune response anti-viral state) and fit them to data from persons with SARS-CoV-2 treated with a mAb. We showed that both models can explain the emergence of resistant virus associated with high transient viral rebounds. We found that variations in the target cell supply rate and adaptive immunity parameters have a strong impact on the magnitude or observability of the viral rebound associated with the emergence of resistant virus. Both variations in target cell supply rate and adaptive immunity parameters may explain why only some individuals develop observable transient resistant viral rebound. Our study highlights the conditions that can lead to resistance and subsequent viral rebound in mAb treatments during acute infection.

Variant-Specific Viral Kinetics in Acute COVID-19.

Understanding variant-specific differences in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) viral kinetics may explain differences in transmission efficiency and provide insights on pathogenesis and prevention. We evaluated SARS-CoV-2 kinetics from nasal swabs across multiple variants (Alpha, Delta, Epsilon, Gamma) in placebo recipients of the ACTIV-2/A5401 trial. Delta variant infection led to the highest maximum viral load and shortest time from symptom onset to viral load peak. There were no significant differences in time to viral clearance across the variants. Viral decline was biphasic with first- and second-phase decays having half-lives of 11 hours and 2.5 days, respectively, with differences among variants, especially in the second phase. These results suggest that while variant-specific differences in viral kinetics exist, post-peak viral load all variants appeared to be efficiently cleared by the host. Clinical Trials Registration. NCT04518410.

The latent reservoir of inducible, infectious HIV-1 does not decrease despite decades of antiretroviral therapy.

HIV-1 persists in a latent reservoir in resting CD4+ T cells despite antiretroviral therapy (ART). The reservoir decays slowly over the first 7 years of ART (t1/2 = 44 months). However, whether decay continues with long-term ART is unclear. Recent integration site studies indicate gradual selection against inducible, intact proviruses, raising speculation that decades of ART might allow treatment interruption without viral rebound. Therefore, we measured the reservoir in 42 people on long-term ART (mean 22 years) using a quantitative viral outgrowth assay. After 7 years of ART, there was no long-term decrease in the frequency of inducible, replication-competent proviruses but rather an increase with an estimated doubling time of 23 years. Another reservoir assay, the intact proviral DNA assay, confirmed that reservoir decay with t1/2 of 44 months did not continue with long-term ART. The lack of decay reflected proliferation of infected cells. Most inducible, replication-competent viruses (79.8%) had env sequences identical to those of other isolates from the same sample. Thus, although integration site analysis indicates changes in reservoir composition, the proliferation of CD4+ T cells counteracts decay, maintaining the frequency of inducible, replication-competent proviruses at roughly constant levels over the long term. These results reinforce the need for lifelong ART.

An explanation for SARS-CoV-2 rebound after Paxlovid treatment.

In a fraction of SARS-CoV-2 infected individuals treated with the oral antiviral Paxlovid, the virus rebounds following treatment. The mechanism driving rebound is not understood. Here, we show that viral dynamic models based on the hypothesis that Paxlovid treatment near the time of symptom onset halts the depletion of target cells, but may not fully eliminate the virus, which can lead to viral rebound. We also show that the occurrence of viral rebound is sensitive to model parameters, and the time treatment is initiated, which may explain why only a fraction of individuals develop viral rebound. Finally, the models are used to test the therapeutic effects of two alternative treatment schemes. These findings also provide a possible explanation for rebounds following other antiviral treatments for SARS-CoV-2.

Making waves: Integrating wastewater surveillance with dynamic modeling to track and predict viral outbreaks.

Wastewater surveillance has proved to be a valuable tool to track the COVID-19 pandemic. However, most studies using wastewater surveillance data revolve around establishing correlations and lead time relative to reported case data. In this perspective, we advocate for the integration of wastewater surveillance data with dynamic within-host and between-host models to better understand, monitor, and predict viral disease outbreaks. Dynamic models overcome emblematic difficulties of using wastewater surveillance data such as establishing the temporal viral shedding profile. Complementarily, wastewater surveillance data bypasses the issues of time lag and underreporting in clinical case report data, thus enhancing the utility and applicability of dynamic models. The integration of wastewater surveillance data with dynamic models can enhance real-time tracking and prevalence estimation, forecast viral transmission and intervention effectiveness, and most importantly, provide a mechanistic understanding of infectious disease dynamics and the driving factors. Dynamic modeling of wastewater surveillance data will advance the development of a predictive and responsive monitoring system to improve pandemic preparedness and population health.