Mathematical modeling to reveal breakthrough mechanisms in the HIV Antibody Mediated Prevention (AMP) trials

Modelling Role of Protective and Nonprotective HLA Allele Inducing Different HIV Infection Outcomes

Human immunodeficiency virus (HIV) infects CD4+ cells and causes progressive immune function failure, and CD8+ cells lyse infected CD4+ cell via recognising peptide presented by human leukocyte antigens (HLA). Variations in HLA allele lead to observed different HIV infection outcomes. Within-host HIV dynamics involves virus replication within infected cells and lysing of infected cells by CD8+ cells, but how variations in HLA alleles determine different infection outcomes was far from clear. Here, we used mathematical modelling and parameter inference with a new analysis of published virus inhibition assay data to estimate CD8+ cell lysing efficiency, and found that lysing efficiency fall in the gap between low bound (0.1-0.2 day-1 (Elemans et al. in PLoS Comput Biol 8(2):e1002381, 2012)) and upper boundary (6.5-8.4 day-1 (Wick et al. in J Virol 79(21):13579-13586, 2005)). Our outcomes indicate that both lysing efficiency and viral inoculum size jointly determine observed different infection outcomes. Low lysing rate associated with non-protective HLA alleles leads to monostable viral kinetic to high viral titre and oscillatory viral kinetics. High lysing rate associated with protective HLA alleles leads monostable viral kinetic to low viral titre and bistable viral kinetics; at a specific interval of CD8+ cell counts, small viral inoculum sizes are inhibited but not large viral inoculum sizes remain infectious. Further, with CD8+ cell recruitment, HIV kinetics always exhibit oscillatory kinetics, but lysing rate is negatively correlated with range of CD8+ cell count. Our finding highlights role of HLA allele determining different infection outcomes, thereby providing a potential mechanistic explanation for observed good and bad HIV infection outcomes induced by protective HLA allele.

A unifying model to explain frequent SARS-CoV-2 rebound after nirmatrelvir treatment and limited prophylactic efficacy

In a pivotal trial (EPIC-HR), a 5-day course of oral ritonavir-boosted nirmatrelvir, given early during symptomatic SARS-CoV-2 infection (within three days of symptoms onset), decreased hospitalization and death by 89.1% and nasal viral load by 0.87 log relative to placebo in high-risk individuals. Yet, nirmatrelvir/ritonavir failed as post-exposure prophylaxis in a trial, and frequent viral rebound has been observed in subsequent cohorts. We develop a mathematical model capturing viral-immune dynamics and nirmatrelvir pharmacokinetics that recapitulates viral loads from this and another clinical trial (PLATCOV). Our results suggest that nirmatrelvir's in vivo potency is significantly lower than in vitro assays predict. According to our model, a maximally potent agent would reduce the viral load by approximately 3.5 logs relative to placebo at 5 days. The model identifies that earlier initiation and shorter treatment duration are key predictors of post-treatment rebound. Extension of treatment to 10 days for Omicron variant infection in vaccinated individuals, rather than increasing dose or dosing frequency, is predicted to lower the incidence of viral rebound significantly.

A unifying model to explain high nirmatrelvir therapeutic efficacy against SARS-CoV-2, despite low post-exposure prophylaxis efficacy and frequent viral rebound

In a pivotal trial (EPIC-HR), a 5-day course of oral ritonavir-boosted nirmatrelvir, given early during symptomatic SARS-CoV-2 infection (within three days of symptoms onset), decreased hospitalization and death by 89.1% and nasal viral load by 0.87 log relative to placebo in high-risk individuals. Yet, nirmatrelvir/ritonavir failed as post-exposure prophylaxis in a trial, and frequent viral rebound has been observed in subsequent cohorts. We developed a mathematical model capturing viral-immune dynamics and nirmatrelvir pharmacokinetics that recapitulated viral loads from this and another clinical trial (PLATCOV). Our results suggest that nirmatrelvir's in vivo potency is significantly lower than in vitro assays predict. According to our model, a maximally potent agent would reduce the viral load by approximately 3.5 logs relative to placebo at 5 days. The model identifies that earlier initiation and shorter treatment duration are key predictors of post-treatment rebound. Extension of treatment to 10 days for Omicron variant infection in vaccinated individuals, rather than increasing dose or dosing frequency, is predicted to lower the incidence of viral rebound significantly.

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.

High monoclonal neutralization titers reduced breakthrough HIV-1 viral loads in the Antibody Mediated Prevention trials

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.

Application of the SLAPNAP statistical learning tool to broadly neutralizing antibody HIV prevention research

Combination monoclonal broadly neutralizing antibody (bnAb) regimens are in clinical development for HIV prevention, necessitating additional knowledge of bnAb neutralization potency/breadth against circulating viruses. Williamson et al. (2021) described a software tool, Super LeArner Prediction of NAb Panels (SLAPNAP), with application to any HIV bnAb regimen with sufficient neutralization data against a set of viruses in the Los Alamos National Laboratory's Compile, Neutralize, and Tally Nab Panels repository. SLAPNAP produces a proteomic antibody resistance (PAR) score for Env sequences based on predicted neutralization resistance and estimates variable importance of Env amino acid features. We apply SLAPNAP to compare HIV bnAb regimens undergoing clinical testing, finding improved power for downstream sieve analyses and increased precision for comparing neutralization potency/breadth of bnAb regimens due to the inclusion of PAR scores of Env sequences with much larger sample sizes available than for neutralization outcomes. SLAPNAP substantially improves bnAb regimen characterization, ranking, and down-selection.

Extended safety and tolerability of subcutaneous CAP256V2LS and VRC07-523LS in HIV-negative women: study protocol for the randomised, placebo-controlled double-blinded, phase 2 CAPRISA 012C trial

INTRODUCTION

Women-controlled HIV prevention technologies that overcome adherence challenges of available daily oral pre-exposure prophylaxis and give women a choice of options are urgently needed. Broadly neutralising monoclonal antibodies (bnAbs) administered passively may offer a valuable non-antiretroviral biological intervention for HIV prevention. Animal and human studies have demonstrated that bnAbs which neutralise HIV can prevent infection. The optimal plasma antibody concentrations to confer protection against HIV infection in humans is under intense study. The Centre for the AIDS Programme of Research in South Africa (CAPRISA) 012C trial will evaluate extended safety and pharmacokinetics of CAP256V2LS and VRC07-523LS among young HIV-negative South African and Zambian women. The study design also allows for an evaluation of a signal of HIV prevention efficacy.

METHODS AND ANALYSIS

CAPRISA 012 is a series of trials with three distinct protocols. The completed CAPRISA 012A and 012B phase 1 trials provided critical data for the CAPRISA 012C trial, which is divided into parts A and B. In part A, 90 participants were randomised to receive both CAP256V2LS and VRC07-523LS at 20 mg/kg or placebo, subcutaneously every 16 or 24 weeks. Part B will enrol 900 participants in South Africa and Zambia who will be randomised in a 1:1 ratio and receive an initial loading dose of 1.2 g of CAP256V2LS and VRC07-523LS or placebo followed by 600 mg of CAP256V2LS and 1.2 g of VRC07-523LS or placebo subcutaneously every 6 months. Safety will be assessed by frequency and severity of reactogenicity and other related adverse events. Pharmacokinetics of both antibodies will be measured in systemic and mucosal compartments over time, while participants will be monitored for breakthrough HIV infections.

ETHICS AND DISSEMINATION OF STUDY FINDINGS

The University of KwaZulu-Natal Biomedical Research Ethics Committee and South African Health Products Regulatory Authority have approved the trial (BREC/00002492/2021, SAHPRA20210317). Results will be disseminated through conference presentations, peer-reviewed publications and the clinical trial registry.

TRIAL REGISTRATION NUMBER

PACTR202112683307570.

Pharmacokinetic serum concentrations of VRC01 correlate with prevention of HIV-1 acquisition

BACKGROUND

The phase 2b proof-of-concept Antibody Mediated Prevention (AMP) trials showed that VRC01, an anti-HIV-1 broadly neutralising antibody (bnAb), prevented acquisition of HIV-1 sensitive to VRC01. To inform future study design and dosing regimen selection of candidate bnAbs, we investigated the association of VRC01 serum concentration with HIV-1 acquisition using AMP trial data.

METHODS

The case-control sample included 107 VRC01 recipients who acquired HIV-1 and 82 VRC01 recipients who remained without HIV-1 during the study. We measured VRC01 serum concentrations with a qualified pharmacokinetic (PK) Binding Antibody Multiplex Assay. We employed nonlinear mixed effects PK modelling to estimate daily-grid VRC01 concentrations. Cox regression models were used to assess the association of VRC01 concentration at exposure and baseline body weight, with the hazard of HIV-1 acquisition and prevention efficacy as a function of VRC01 concentration. We also compared fixed dosing vs. body weight-based dosing via simulations.

FINDINGS

Estimated VRC01 concentrations in VRC01 recipients without HIV-1 were higher than those in VRC01 recipients who acquired HIV-1. Body weight was inversely associated with HIV-1 acquisition among both placebo and VRC01 recipients but did not modify the prevention efficacy of VRC01. VRC01 concentration was inversely correlated with HIV-1 acquisition, and positively correlated with prevention efficacy of VRC01. Simulation studies suggest that fixed dosing may be comparable to weight-based dosing in overall predicted prevention efficacy.

INTERPRETATION

These findings suggest that bnAb serum concentration may be a useful marker for dosing regimen selection, and operationally efficient fixed dosing regimens could be considered for future trials of HIV-1 bnAbs.

FUNDING

Was provided by the National Institutes of Health, National Institute of Allergy and Infectious Diseases (NIAID) (UM1 AI068614, to the HIV Vaccine Trials Network [HVTN]; UM1 AI068635, to the HVTN Statistical Data and Management Center [SDMC], Fred Hutchinson Cancer Center [FHCC]; 2R37 054165 to the FHCC; UM1 AI068618, to HVTN Laboratory Center, FHCC; UM1 AI068619, to the HPTN Leadership and Operations Center; UM1 AI068613, to the HIV Prevention Trials Network [HPTN] Laboratory Center; UM1 AI068617, to the HPTN SDMC; and P30 AI027757, to the Center for AIDS Research, Duke University (AI P30 AI064518) and University of Washington (P30 AI027757) Centers for AIDS Research; R37AI054165 from NIAID to the FHCC; and OPP1032144 CA-VIMC Bill & Melinda Gates Foundation.

Neutralization titer biomarker for antibody-mediated prevention of HIV-1 acquisition

The Antibody Mediated Prevention trials showed that the broadly neutralizing antibody (bnAb) VRC01 prevented acquisition of human immunodeficiency virus-1 (HIV-1) sensitive to VRC01. Using AMP trial data, here we show that the predicted serum neutralization 80% inhibitory dilution titer (PT80) biomarker-which quantifies the neutralization potency of antibodies in an individual's serum against an HIV-1 isolate-can be used to predict HIV-1 prevention efficacy. Similar to the results of nonhuman primate studies, an average PT80 of 200 (meaning a bnAb concentration 200-fold higher than that required to reduce infection by 80% in vitro) against a population of probable exposing viruses was estimated to be required for 90% prevention efficacy against acquisition of these viruses. Based on this result, we suggest that the goal of sustained PT80 <200 against 90% of circulating viruses can be achieved by promising bnAb regimens engineered for long half-lives. We propose the PT80 biomarker as a surrogate endpoint for evaluatinon of bnAb regimens, and as a tool for benchmarking candidate bnAb-inducing vaccines.

Design of an optimal combination therapy with broadly neutralizing antibodies to suppress HIV-1

Infusion of broadly neutralizing antibodies (bNAbs) has shown promise as an alternative to anti-retroviral therapy against HIV. A key challenge is to suppress viral escape, which is more effectively achieved with a combination of bNAbs. Here, we propose a computational approach to predict the efficacy of a bNAb therapy based on the population genetics of HIV escape, which we parametrize using high-throughput HIV sequence data from bNAb-naive patients. By quantifying the mutational target size and the fitness cost of HIV-1 escape from bNAbs, we predict the distribution of rebound times in three clinical trials. We show that a cocktail of three bNAbs is necessary to effectively suppress viral escape, and predict the optimal composition of such bNAb cocktail. Our results offer a rational therapy design for HIV, and show how genetic data can be used to predict treatment outcomes and design new approaches to pathogenic control.

Reappraising the Value of HIV-1 Vaccine Correlates of Protection Analyses

With the much-debated exception of the modestly reduced acquisition reported for the RV144 efficacy trial, HIV-1 vaccines have not protected humans against infection, and a vaccine of similar design to that tested in RV144 was not protective in a later trial, HVTN 702. Similar vaccine regimens have also not consistently protected nonhuman primates (NHPs) against viral acquisition. Conversely, experimental vaccines of different designs have protected macaques from viral challenges but then failed to protect humans, while many other HIV-1 vaccine candidates have not protected NHPs. While efficacy varies more in NHPs than humans, vaccines have failed to protect in the most stringent NHP model. Intense investigations have aimed to identify correlates of protection (CoPs), even in the absence of net protection. Unvaccinated animals and humans vary vastly in their susceptibility to infection and in their innate and adaptive responses to the vaccines; hence, merely statistical associations with factors that do not protect are easily found. Systems biological analyses, including artificial intelligence, have identified numerous candidate CoPs but with no clear consistency within or between species. Proposed CoPs sometimes have only tenuous mechanistic connections to immune protection. In contrast, neutralizing antibodies (NAbs) are a central mechanistic CoP for vaccines that succeed against other viruses, including SARS-CoV-2. No HIV-1 vaccine candidate has yet elicited potent and broadly active NAbs in NHPs or humans, but narrow-specificity NAbs against the HIV-1 isolate corresponding to the immunogen do protect against infection by the autologous virus. Here, we analyze why so many HIV-1 vaccines have failed, summarize the outcomes of vaccination in NHPs and humans, and discuss the value and pitfalls of hunting for CoPs other than NAbs. We contrast the failure to find a consistent CoP for HIV-1 vaccines with the identification of NAbs as the principal CoP for SARS-CoV-2.

Optimizing clinical dosing of combination broadly neutralizing antibodies for HIV prevention

Broadly neutralizing antibodies (bNAbs) are promising agents to prevent HIV infection and achieve HIV remission without antiretroviral therapy (ART). As with ART, bNAb combinations are likely needed to cover HIV's extensive diversity. Not all bNAbs are identical in terms of their breadth, potency, and in vivo longevity (half-life). Given these differences, it is important to optimally select the composition, or dose ratio, of combination bNAb therapies for future clinical studies. We developed a model that synthesizes 1) pharmacokinetics, 2) potency against a wide HIV diversity, 3) interaction models for how drugs work together, and 4) correlates that translate in vitro potency to clinical protection. We found optimization requires drug-specific balances between potency, longevity, and interaction type. As an example, tradeoffs between longevity and potency are shown by comparing a combination therapy to a bi-specific antibody (a single protein merging both bNAbs) that takes the better potency but the worse longevity of the two components. Then, we illustrate a realistic dose ratio optimization of a triple combination of VRC07, 3BNC117, and 10-1074 bNAbs. We apply protection estimates derived from both a non-human primate (NHP) challenge study meta-analysis and the human antibody mediated prevention (AMP) trials. In both cases, we find a 2:1:1 dose emphasizing VRC07 is nearly optimal. Our approach can be immediately applied to optimize the next generation of combination antibody prevention and cure studies.

Estimation of the in vivo neutralization potency of eCD4Ig and conditions for AAV-mediated production for SHIV long-term remission

The engineered protein eCD4Ig has emerged as a promising approach to achieve HIV remission in the absence of antiviral therapy. eCD4Ig neutralizes nearly all HIV-1 isolates and induces antibody-dependent cell-mediated cytotoxicity (ADCC) in vitro. To characterize the in vivo antiviral neutralization and possible ADCC effects of eCD4Ig, we fit mathematical models to eCD4Ig, anti–eCD4Ig-drug antibody (ADA), and viral load kinetics from healthy and simian-human immunodeficiency virus AD8 (SHIV-AD8) infected nonhuman primates that were treated with single or sequentially dosed eCD4Ig passive administrations. Our model predicts that eCD4Ig transiently decreases SHIV viral loads due to neutralization only with an in vivo IC₅₀ of ~25 μg/ml but with limited effect due to ADA. Simulations suggest that endogenous, continuous expression of eCD4Ig at levels greater than 10⁵ μg/day, as is possible with Adeno-associated virus (AAV) vector-based production, could overcome the diminishing effects of ADA and allow for long-term remission of SHIV viremia in nonhuman primates.

Stochastic investigation of HIV infection and the emergence of drug resistance

Drug-resistant HIV-1 has caused a growing concern in clinic and public health. Although combination antiretroviral therapy can contribute massively to the suppression of viral loads in patients with HIV-1, it cannot lead to viral eradication. Continuing viral replication during sub-optimal therapy (due to poor adherence or other reasons) may lead to the accumulation of drug resistance mutations, resulting in an increased risk of disease progression. Many studies also suggest that events occurring during the early stage of HIV-1 infection (i.e., the first few hours to days following HIV exposure) may determine whether the infection can be successfully established. However, the numbers of infected cells and viruses during the early stage are extremely low and stochasticity may play a critical role in dictating the fate of infection. In this paper, we use stochastic models to investigate viral infection and the emergence of drug resistance of HIV-1. The stochastic model is formulated by a continuous-time Markov chain (CTMC), which is derived based on an ordinary differential equation model proposed by Kitayimbwa et al. that includes both forward and backward mutations. An analytic estimate of the probability of the clearance of HIV infection of the CTMC model near the infection-free equilibrium is obtained by a multitype branching process approximation. The analytical predictions are validated by numerical simulations. Unlike the deterministic dynamics where the basic reproduction number R₀ serves as a sharp threshold parameter (i.e., the disease dies out if R₀<1 and persists if R₀>1), the stochastic models indicate that there is always a positive probability for HIV infection to be eradicated in patients. In the presence of antiretroviral therapy, our results show that the chance of clearance of the infection tends to increase although drug resistance is likely to emerge.

Mathematical Modeling of Within-Host, Untreated, Cytomegalovirus Infection Dynamics after Allogeneic Transplantation

Cytomegalovirus (CMV) causes significant morbidity and mortality in recipients of allogeneic hematopoietic cell transplantation (HCT). Whereas insights gained from mathematical modeling of other chronic viral infections such as HIV, hepatitis C, and herpes simplex virus-2 have aided in optimizing therapy, previous CMV modeling has been hindered by a lack of comprehensive quantitative PCR viral load data from untreated episodes of viremia in HCT recipients. We performed quantitative CMV DNA PCR on stored, frozen serum samples from the placebo group of participants in a historic randomized controlled trial of ganciclovir for the early treatment of CMV infection in bone marrow transplant recipients. We developed four main ordinary differential Equation mathematical models and used model selection theory to choose between 38 competing versions of these models. Models were fit using a population, nonlinear, mixed-effects approach. We found that CMV kinetics from untreated HCT recipients are highly variable. The models that recapitulated the observed patterns most parsimoniously included explicit, dynamic immune cell compartments and did not include dynamic target cell compartments, consistent with the large number of tissue and cell types that CMV infects. In addition, in our best-fitting models, viral clearance was extremely slow, suggesting severe impairment of the immune response after HCT. Parameters from our best model correlated well with participants’ clinical risk factors and outcomes from the trial, further validating our model. Our models suggest that CMV dynamics in HCT recipients are determined by host immune response rather than target cell limitation in the absence of antiviral treatment.

Predicting in vivo escape dynamics of HIV-1 from a broadly neutralizing antibody

Broadly neutralizing antibodies are promising candidates for treatment and prevention of HIV-1 infections. Such antibodies can temporarily suppress viral load in infected individuals; however, the virus often rebounds by escape mutants that have evolved resistance. In this paper, we map a fitness model of HIV-1 interacting with broadly neutralizing antibodies using in vivo data from a recent clinical trial. We identify two fitness factors, antibody dosage and viral load, that determine viral reproduction rates reproducibly across different hosts. The model successfully predicts the escape dynamics of HIV-1 in the course of an antibody treatment, including a characteristic frequency turnover between sensitive and resistant strains. This turnover is governed by a dosage-dependent fitness ranking, resulting from an evolutionary trade-off between antibody resistance and its collateral cost in drug-free growth. Our analysis suggests resistance-cost trade-off curves as a measure of antibody performance in the presence of resistance evolution.

Timing HIV infection with a simple and accurate population viral dynamics model

Clinical trials for HIV prevention can require knowledge of infection times to subsequently determine protective drug levels. Yet, infection timing is difficult when study visits are sparse. Using population nonlinear mixed-effects (pNLME) statistical inference and viral loads from 46 RV217 study participants, we developed a relatively simple HIV primary infection model that achieved an excellent fit to all data. We also discovered that Aptima assay values from the study strongly correlated with viral loads, enabling imputation of very early viral loads for 28/46 participants. Estimated times between infecting exposures and first positives were generally longer than prior estimates (average of two weeks) and were robust to missing viral upslope data. On simulated data, we found that tighter sampling before diagnosis improved estimation more than tighter sampling after diagnosis. Sampling weekly before and monthly after diagnosis was a pragmatic design for good timing accuracy. Our pNLME timing approach is widely applicable to other infections with existing mathematical models. The present model could be used to simulate future HIV trials and may help estimate protective thresholds from the recently completed antibody-mediated prevention trials.

Neutralization diversity of HIV-1 Indian subtype C envelopes obtained from cross sectional and followed up individuals against broadly neutralizing monoclonal antibodies having distinct gp120 specificities

BACKGROUND

The potential use of the broadly neutralizing monoclonal antibodies (bnAbs) towards prophylaxis and treatment to HIV-1 is currently being explored. While a number of promising bnAbs have been discovered and a few of them have progressed towards clinical development, their extent of neutralization coverage with respect to global HIV-1 variants given the existence of genetically distinct subtypes and recombinants circulating globally is not clearly known. In the present study, we examined the variation in the neutralization susceptibility of pseudoviruses expressing 71 full length primary HIV-1 subtype C envs obtained from limited cross-sectional individuals over different time points against four bnAbs that target gp120 with distinct specificities: VRC01, CAP256-VRC26.25, PGDM1400 and PGT121.

RESULTS

We found significant variations in the susceptibility of Indian clade C to these four bnAbs. These variations were found to be distinct to that observed in African subtype C based on the existing datasets and concordant with their sequence diversity. Trend analysis indicated an increasing neutralization resistance observed over time with CAP25-VRC26.25, PGDM1400 and PGT121 when tested on pseudoviruses expressing envs obtained from 1999 to 2016. However, inconsistent trend in neutralization susceptibility was observed, when pseudoviruses expressing envs obtained from three followed up individuals were examined. Finally, through predictive analysis of the 98 Indian subtype C including those assessed in the present study by employing additive model implemented in CombiNAber ( http://www.hiv.lanl.gov ), we observed two possibilities where combinations of three bnAbs (VRC01/CAP56-VRC26.25/PGT121 and PGDM1400/CAP256-VRC26.25/PGT121) could achieve near 100% neutralization coverage.

CONCLUSIONS

Our findings not only indicate disparate intra-clade C genetic vis-à-vis neutralization diversities but also warrant the need for more comprehensive study using additional isolates towards comparing inter and intra-clade neutralization diversities which will be necessary for selecting the bnAb combinations suitable for optimal coverage of the region-specific HIV-1 circulating subtypes. Expanding these efforts is imperative for designing efficacious bnAb based intervention strategies for India as well as subtype C in general.

Success of prophylactic antiviral therapy for SARS-CoV-2: Predicted critical efficacies and impact of different drug-specific mechanisms of action

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.

Pre-existing resistance in the latent reservoir can compromise VRC01 therapy during chronic HIV-1 infection

Passive immunization with broadly neutralizing antibodies (bNAbs) of HIV-1 appears a promising strategy for eliciting long-term HIV-1 remission. When administered concomitantly with the cessation of antiretroviral therapy (ART) to patients with established viremic control, bNAb therapy is expected to prolong remission. Surprisingly, in clinical trials on chronic HIV-1 patients, the bNAb VRC01 failed to prolong remission substantially. Identifying the cause of this failure is important for improving VRC01-based therapies and unraveling potential vulnerabilities of other bNAbs. In the trials, viremia resurged rapidly in most patients despite suppressive VRC01 concentrations in circulation, suggesting that VRC01 resistance was the likely cause of failure. ART swiftly halts viral replication, precluding the development of resistance during ART. If resistance were to emerge post ART, virological breakthrough would have taken longer than without VRC01 therapy. We hypothesized therefore that VRC01-resistant strains must have been formed before ART initiation, survived ART in latently infected cells, and been activated during VRC01 therapy, causing treatment failure. Current assays preclude testing this hypothesis experimentally. We developed a mathematical model based on the hypothesis and challenged it with available clinical data. The model integrated within-host HIV-1 evolution, stochastic latency reactivation, and viral dynamics with multiple-dose VRC01 pharmacokinetics. The model predicted that single but not higher VRC01-resistant mutants would pre-exist in the latent reservoir. We constructed a virtual patient population that parsimoniously recapitulated inter-patient variations. Model predictions with this population quantitatively captured data of VRC01 failure from clinical trials, presenting strong evidence supporting the hypothesis. We attributed VRC01 failure to single-mutant VRC01-resistant proviruses in the latent reservoir triggering viral recrudescence, particularly when VRC01 was at trough levels. Pre-existing resistant proviruses in the latent reservoir may similarly compromise other bNAbs. Our study provides a framework for designing bNAb-based therapeutic protocols that would avert such failure and maximize HIV-1 remission.

Antiretroviral therapy (ART) can control but not eradicate HIV-1. Stopping ART leads to rapid viral resurgence and progressive disease. ART is therefore administered lifelong. Tremendous efforts are ongoing to devise strategies that will enable stopping ART and yet prevent viral resurgence. One such strategy involves the administration of broadly neutralizing antibodies (bNAbs) of HIV-1 at the time of stopping ART. This strategy is expected to delay if not prevent viral resurgence. Surprisingly, treatment with VRC01, a potent bNAb, resulted in hardly any improvement in viral remission. In this study, we elucidate the cause of this failure. We hypothesized that VRC01-resistant strains may pre-exist in latently infected cells, which are unaffected by ART. They can thus outlast ART and get reactivated, triggering VRC01 failure. We built a detailed mathematical model based on this hypothesis and showed that it quantitatively captured observations of VRC01 failure in clinical trials on chronic HIV-1 patients. Our study thus identifies a potential vulnerability of bNAbs, namely, bNAb-resistant strains pre-existing in latently infected cells. Our model offers a framework for predicting bNAb-based treatment protocols that would preclude failure due to pre-existing resistance and maximally prolong remission.

Potency and timing of antiviral therapy as determinants of duration of SARS-CoV-2 shedding and intensity of inflammatory response

To affect the COVID-19 pandemic, lifesaving antiviral therapies must be identified. The number of clinical trials that can be performed is limited. We developed mathematical models to project multiple therapeutic approaches. Our models recapitulate off-treatment viral dynamics and predict a three-phase immune response. Simulated treatment with remdesivir, selinexor, neutralizing antibodies, or cellular immunotherapy demonstrates that rapid viral elimination is possible if in vivo potency is sufficiently high. Therapies dosed soon after peak viral load when symptoms develop may decrease shedding duration and immune response intensity but have little effect on viral area under the curve (AUC), which is driven by high early viral loads. Potent therapy dosed before viral peak during presymptomatic infection could lower AUC. Drug resistance may emerge with a moderately potent agent dosed before viral peak. Our results support early treatment for COVID-19 if shedding duration, not AUC, is most predictive of clinical severity.

Potency and timing of antiviral therapy as determinants of duration of SARS-CoV-2 shedding and intensity of inflammatory response

To affect the COVID-19 pandemic, lifesaving antiviral therapies must be identified. The number of clinical trials that can be performed is limited. We developed mathematical models to project multiple therapeutic approaches. Our models recapitulate off-treatment viral dynamics and predict a three-phase immune response. Simulated treatment with remdesivir, selinexor, neutralizing antibodies, or cellular immunotherapy demonstrates that rapid viral elimination is possible if in vivo potency is sufficiently high. Therapies dosed soon after peak viral load when symptoms develop may decrease shedding duration and immune response intensity but have little effect on viral area under the curve (AUC), which is driven by high early viral loads. Potent therapy dosed before viral peak during presymptomatic infection could lower AUC. Drug resistance may emerge with a moderately potent agent dosed before viral peak. Our results support early treatment for COVID-19 if shedding duration, not AUC, is most predictive of clinical severity.