We developed a mathematical model to predict the outcomes of different possible COVID-19 treatments.

A model of reservoir activation and viral replication is introduced accounting for the production of 2-LTR HIV-1 DNA circles following antiviral intensification with the HIV integrase inhibitor raltegravir, considering contributions of de novo infection events and exogenous sources of infected cells, including quiescent infected cell activation. The model shows that a monotonic increase in measured 2-LTR concentration post intensification is consistent with limited de novo infection primarily maintained by sources of infected cells unaffected by raltegravir, such as quiescent cell activation, while a transient increase in measured 2-LTR concentration is consistent with significant levels of efficient ( R ₀ > 1) de novo infection. The model is validated against patient data from the INTEGRAL study and is shown to have a statistically significant fit relative to the null hypothesis of random measurement variation about a mean. We obtain estimates and confidence intervals for the model parameters, including 2-LTR half-life. Seven of the 13 patients with detectable 2-LTR concentrations from the INTEGRAL study have measured 2-LTR dynamics consistent with significant levels of efficient replication of the virus prior to treatment intensification.

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.

The Antibody Mediated Prevention (AMP) trials (NCT02716675 and NCT02568215) demonstrated that passive administration of the broadly neutralizing monoclonal antibody VRC01 could prevent some HIV-1 acquisition events. Here, we use mathematical modeling in a post hoc analysis to demonstrate that VRC01 influenced viral loads in AMP participants who acquired HIV. Instantaneous inhibitory potential (IIP), which integrates VRC01 serum concentration and VRC01 sensitivity of acquired viruses in terms of both IC50 and IC80, follows a dose-response relationship with first positive viral load (p = 0.03), which is particularly strong above a threshold of IIP = 1.6 (r = -0.6, p = 2e-4). Mathematical modeling reveals that VRC01 activity predicted from in vitro IC80s and serum VRC01 concentrations overestimates in vivo neutralization by 600-fold (95% CI: 300–1200). The trained model projects that even if future therapeutic HIV trials of combination monoclonal antibodies do not always prevent acquisition, reductions in viremia and reservoir size could be expected.