While most individuals suffer progressive disease following HIV infection, a small fraction spontaneously controls the infection. Although CD8 T-cells have been implicated in this natural control, their mechanistic roles are yet to be established. Here, we combined mathematical modeling and analysis of previously published data from 16 SIV-infected macaques, of which 12 were natural controllers, to elucidate the role of CD8 T-cells in natural control. For each macaque, we considered, in addition to the canonical in vivo plasma viral load and SIV DNA data, longitudinal ex vivo measurements of the virus suppressive capacity of CD8 T-cells. Available mathematical models do not allow analysis of such combined in vivo-ex vivo datasets. We explicitly modeled the ex vivo assay, derived analytical approximations that link the ex vivo measurements with the in vivo effector function of CD8-T cells, and integrated them with an in vivo model of virus dynamics, thus developing a new learning framework that enabled the analysis. Our model fit the data well and estimated the recruitment rate and/or maximal killing rate of CD8 T-cells to be up to 2-fold higher in controllers than non-controllers (p = 0.013). Importantly, the cumulative suppressive capacity of CD8 T-cells over the first 4-6 weeks of infection was associated with virus control (Spearman's ρ = -0.51; p = 0.05). Thus, our analysis identified the early cumulative suppressive capacity of CD8 T-cells as a predictor of natural control. Furthermore, simulating a large virtual population, our model quantified the minimum capacity of this early CD8 T-cell response necessary for long-term control. Our study presents new, quantitative insights into the role of CD8 T-cells in the natural control of HIV infection and has implications for remission strategies.

Households are a major driver of transmission of acute respiratory viruses, such as SARS-CoV-2 or Influenza. Until now antiviral treatments have mostly been used as a curative treatment in symptomatic individuals. During an outbreak, more aggressive strategies involving pre- or post-exposure prophylaxis (PrEP or PEP) could be employed to further reduce the risk of severe disease but also prevent transmission to household contacts. In order to understand the effectiveness of such strategies and the factors that may modulate them, we developed a multi-scale model that follows the infection at both the individual-level (viral dynamics) and the population-level (transmission dynamics) in households. Using a simulation study we explored different antiviral treatment strategies, evaluating their effectiveness on reducing the transmission risk and the virological burden in households for a range of virus characteristics (e.g., secondary attack rate-SAR, or time to peak viral load). We found that when the index case can be identified and treated before symptom onset, both transmission and virological burden are reduced by > 75% for most SAR values and time to peak viral load, with minimal benefit to treat additionally household contacts. While treatment initiated after index symptom onset does not reduce the risk of transmission, it can still reduce the virological burden in the household, a proxy for severe disease and subsequent transmission risk outside the household. In that case optimal strategies involve treatment of both index case and household contacts as PEP, with efficacy > 50% when peak viral load occurs after symptom onset, and 30-50% otherwise. In all the considered cases, antiviral treatment strategies were optimal for SAR ranging 20-60%, and for larger household sizes. This study highlights the opportunity of antiviral drug-based interventions in households during an outbreak to minimize viral transmission and disease burden.

International audience

The impact of variants of concern (VoC) on SARS-CoV-2 viral dynamics remains poorly understood and essentially relies on observational studies subject to various sorts of biases. In contrast, experimental models of infection constitute a powerful model to perform controlled comparisons of the viral dynamics observed with VoC and better quantify how VoC escape from the immune response. Here we used molecular and infectious viral load of 78 cynomolgus macaques to characterize in detail the effects of VoC on viral dynamics. We first developed a mathematical model that recapitulate the observed dynamics, and we found that the best model describing the data assumed a rapid antigen-dependent stimulation of the immune response leading to a rapid reduction of viral infectivity. When compared with the historical variant, all VoC except beta were associated with an escape from this immune response, and this effect was particularly sensitive for delta and omicron variant (p<10-6 for both). Interestingly, delta variant was associated with a 1.8-fold increased viral production rate (p = 0.046), while conversely omicron variant was associated with a 14-fold reduction in viral production rate (p<10-6). During a natural infection, our models predict that delta variant is associated with a higher peak viral RNA than omicron variant (7.6 log10 copies/mL 95% CI 6.8-8 for delta; 5.6 log10 copies/mL 95% CI 4.8-6.3 for omicron) while having similar peak infectious titers (3.7 log10 PFU/mL 95% CI 2.4-4.6 for delta; 2.8 log10 PFU/mL 95% CI 1.9-3.8 for omicron). These results provide a detailed picture of the effects of VoC on total and infectious viral load and may help understand some differences observed in the patterns of viral transmission of these viruses.

Experimental Zika virus infection in non-human primates results in acute viral load dynamics that can be well-described by mathematical models. The inoculum dose that would be received in a natural infection setting is likely lower than the experimental infections and how this difference affects the viral dynamics and immune response is unclear. Here we study a dataset of experimental infection of non-human primates with a range of doses of Zika virus. We develop new models of infection incorporating both an innate immune response and viral interference with that response. We find that such a model explains the data better than models with no interaction between virus and the immune response. We also find that larger inoculum doses lead to faster dynamics of infection, but approximately the same total amount of viral production.

Repurposed drugs that are safe and immediately available constitute a first line of defense against new viral infections. Despite limited antiviral activity against SARS-CoV-2, several drugs are being tested as medication or as prophylaxis to prevent infection. Using a stochastic model of early phase infection, we evaluate the success of prophylactic treatment with different drug types to prevent viral infection. We find that there exists a critical efficacy that a treatment must reach in order to block viral establishment. Treatment by a combination of drugs reduces the critical efficacy, most effectively by the combination of a drug blocking viral entry into cells and a drug increasing viral clearance. Below the critical efficacy, the risk of infection can nonetheless be reduced. Drugs blocking viral entry into cells or enhancing viral clearance reduce the risk of infection more than drugs that reduce viral production in infected cells. The larger the initial inoculum of infectious virus, the less likely is prevention of an infection. In our model, we find that as long as the viral inoculum is smaller than 10 infectious virus particles, viral infection can be prevented almost certainly with drugs of 90% efficacy (or more). Even when a viral infection cannot be prevented, antivirals delay the time to detectable viral loads. The largest delay of viral infection is achieved by drugs reducing viral production in infected cells. A delay of virus infection flattens the within-host viral dynamic curve, possibly reducing transmission and symptom severity. Thus, antiviral prophylaxis, even with reduced efficacy, could be efficiently used to prevent or alleviate infection in people at high risk.

Non-human primates infected with SARS-CoV-2 exhibit mild clinical signs. Here we used a mathematical model to characterize in detail the viral dynamics in 31 cynomolgus macaques for which nasopharyngeal and tracheal viral load were frequently assessed. We identified that infected cells had a large burst size (>104 virus) and a within-host reproductive basic number of approximately 6 and 4 in nasopharyngeal and tracheal compartment, respectively. After peak viral load, infected cells were rapidly lost with a half-life of 9 hours, with no significant association between cytokine elevation and clearance, leading to a median time to viral clearance of 10 days, consistent with observations in mild human infections. Given these parameter estimates, we predict that a prophylactic treatment blocking 90% of viral production or viral infection could prevent viral growth. In conclusion, our results provide estimates of SARS-CoV-2 viral kinetic parameters in an experimental model of mild infection and they provide means to assess the efficacy of future antiviral treatments.

Lassa fever is an haemorrhagic fever caused by Lassa virus (LASV). There is no vaccine approved against LASV and the only recommended antiviral treatment relies on ribavirin, despite limited evidence of efficacy. Recently, the nucleotide analogue favipiravir showed a high antiviral efficacy, with 100% survival obtained in an otherwise fully lethal non-human primate (NHP) model of Lassa fever. However the mechanism of action of the drug is not known and the absence of pharmacokinetic data limits the translation of these results to the human setting. Here we aimed to better understand the antiviral effect of favipiravir by developping the first mathematical model recapitulating Lassa viral dynamics and treatment. We analyzed the viral dynamics in 24 NHPs left untreated or treated with ribavirin or favipiravir, and we put the results in perspective with those obtained with the same drugs in the context of Ebola infection. Our model estimates favipiravir EC50 in vivo to 2.89 μg.mL-1, which is much lower than what was found against Ebola virus. The main mechanism of action of favipiravir was to decrease virus infectivity, with an efficacy of 91% at the highest dose. Based on our knowledge acquired on the drug pharmacokinetics in humans, our model predicts that favipiravir doses larger than 1200 mg twice a day should have the capability to strongly reduce the production infectious virus and provide a milestone towards a future use in humans.

Fecal excretion of antibiotics and resistant bacteria in the environment are major public health threats associated with extensive farming and modern medical care. Innovative strategies that can reduce the intestinal antibiotic concentrations during treatments are in development. However, the effect of lower exposure on the amount of resistant enterobacteria excreted has not been quantified, making it difficult to anticipate the impact of these strategies. Here, we introduce a bacterial kinetic model to capture the complex relationships between drug exposure, loss of susceptible enterobacteria and growth of resistant strains in the feces of piglets receiving placebo, 1.5 or 15 mg/kg/day ciprofloxacin, a fluoroquinolone, for 5 days. The model could well describe the kinetics of drug susceptible and resistant enterobacteria observed during treatment, and up to 22 days after treatment cessation. Next, the model was used to predict the expected amount of resistant enterobacteria excreted over an average piglet's lifetime (150 days) when varying drug exposure and treatment duration. For the clinically relevant dose of 15 mg/kg/day for 5 days, the total amount of resistant enterobacteria excreted was predicted to be reduced by 75% and 98% when reducing treatment duration to 3 and 1 day treatment, respectively. Alternatively, for a fixed 5-days treatment, the level of resistance excreted could be reduced by 18%, 33%, 57.5% and 97% if 3, 5, 10 and 30 times lower levels of colonic drug concentrations were achieved, respectively. This characterization on in vivo data of the dynamics of resistance to antibiotics in the colonic flora could provide new insights into the mechanism of dissemination of resistance and can be used to design strategies aiming to reduce it.

The current paradigm for studying hepatitis C virus (HCV) dynamics in patients utilizes a standard viral dynamic model that keeps track of uninfected (target) cells, infected cells, and virus. The model does not account for the dynamics of intracellular viral replication, which is the major target of direct-acting antiviral agents (DAAs). Here we describe and study a recently developed multiscale age-structured model that explicitly considers the potential effects of DAAs on intracellular viral RNA production, degradation, and secretion as virus into the circulation. We show that when therapy significantly blocks both intracellular viral RNA production and virus secretion, the serum viral load decline has three phases, with slopes reflecting the rate of serum viral clearance, the rate of loss of intracellular viral RNA, and the rate of loss of intracellular replication templates and infected cells, respectively. We also derive analytical approximations of the multiscale model and use one of them to analyze data from patients treated for 14 days with the HCV protease inhibitor danoprevir. Analysis suggests that danoprevir significantly blocks intracellular viral production (with mean effectiveness 99.2%), enhances intracellular viral RNA degradation about 5-fold, and moderately inhibits viral secretion (with mean effectiveness 56%). The multiscale model can be used to study viral dynamics in patients treated with other DAAs and explore their mechanisms of action in treatment of hepatitis C.