Immune-viral dynamics modeling for SARS-CoV-2 drug development

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 frequently occurs 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 the relative change of HBV RNA more accurately reflects the antiviral effectiveness of a CAM than the relative change in HBV DNA.

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 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 N-R near the time of symptom onset, coupled with incomplete viral clearance, appears to be the main factor leading to viral rebound. Moreover, the occurrence of viral rebound is likely influenced by the time of treatment initiation relative to the progression of the infection, with earlier treatments leading to a higher chance of rebound. A comparison with an untreated cohort suggests that early treatments with nirmatrelvir-ritonavir may be associated with a delay in the onset of an adaptive immune response. Nevertheless, our model demonstrates that extending the course of nirmatrelvir-ritonavir treatment to a 10-day regimen may greatly diminish the chance of 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.

IMPORTANCE

Nirmatrelvir-ritonavir is an effective treatment for SARS-CoV-2. In a subset of individuals treated with nirmatrelvir-ritonavir, the initial reduction in viral load is followed by viral rebound once treatment is stopped. We show that the timing of treatment initiation with nirmatrelvir-ritonavir may influence the risk of viral rebound. Nirmatrelvir-ritonavir stops viral growth and preserves target cells but may not lead to full clearance of the virus. Thus, once treatment ends, if an effective adaptive immune response has not adequately developed, the remaining virus can lead to rebound. Our results provide insights into the mechanisms of rebound and can help develop better treatment strategies to minimize this possibility.

Exploring the treatment of SARS-CoV-2 with modified vesicular stomatitis virus

SARS-CoV-2 caused a global pandemic and is now an endemic virus that will require continued antiviral and vaccine development. A possible new treatment modality was recently suggested that would use vesicular stomatitis virus (VSV) modified to express the ACE2 receptor. Since the modified VSV expresses the cell surface receptor that is used by the SARS-CoV-2 spike protein, the thought is that SARS-CoV-2 virions would bind to the modified VSV and thus be neutralized. Additionally, since SARS-CoV-2 infected cells also express the spike protein, the modified VSV could potentially infect these cells, allowing for its own replication, but also potentially interfering with replication of SARS-CoV-2. This idea has not yet been tested experimentally, but we can investigate the feasibility of this possible treatment theoretically. In this manuscript, we develop a mathematical model of this suggested treatment and explore conditions under which it might be effective. We find that treatment with modified VSV does little to change the SARS-CoV-2 time course except when the treatment is applied at the onset of the SARS-CoV-2 infection at very high doses. In this case, VSV reduces the peak SARS-CoV-2 viral load, but lengthens the duration of the SARS-CoV-2 infection. Thus, we find that modified VSV treatment is unlikely to be effective largely because it does not prevent infection of cells by SARS-CoV-2.

Development of a SARS-CoV-2 viral dynamic model for patients with COVID-19 based on the amount of viral RNA and viral titer

The target-cell limited model, which is one of the mathematical modeling approaches providing a quantitative understanding of viral dynamics, has been applied to describe viral RNA profiles of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in previous studies. However, these models have been developed mainly using patient data from the early phase of the pandemic. Furthermore, no reports focused on the profiles of the viral titer. In this study, the dynamics of both viral RNA and viral titer were characterized using data reflecting the current clinical situation in which the Omicron variant has become epidemic and vaccines for SARS-CoV-2 have become available. Consecutive data for 5212 viral RNA levels and 5216 viral titers were obtained from 720 patients with coronavirus disease 2019 (COVID-19) in a phase II/III study for ensitrelvir. Our model assumed that productively infected cells would produce only infectious viruses, which could be transformed into non-infectious viruses, and has been used to describe the dynamics of both viral RNA levels and viral titer. The time from infection to symptom onset (tinf) of unvaccinated patients was estimated to be 3.0 days, which was shorter than that of the vaccinated patients. The immune-related parameter as a power function for the vaccinated patients was 1.1 times stronger than that for the unvaccinated patients. Our model allows the prediction of the viral dynamics in patients with COVID-19 from the time of infection to symptom onset. Vaccination status was identified as a factor influencing tinf and the immune function.

Within-host evolution of SARS-CoV-2: how often are de novo mutations transmitted from symptomatic infections?

Despite a relatively low mutation rate, the large number of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections has allowed for substantial genetic change, leading to a multitude of emerging variants. Using a recently determined mutation rate (per site replication), as well as within-host parameter estimates for symptomatic SARS-CoV-2 infection, we apply a stochastic transmission-bottleneck model to describe the survival probability of de novo SARS-CoV-2 mutations as a function of bottleneck size and selection coefficient. For narrow bottlenecks, we find that mutations affecting per-target-cell attachment rate (with phenotypes associated with fusogenicity and ACE2 binding) have similar transmission probabilities to mutations affecting viral load clearance (with phenotypes associated with humoral evasion). We further find that mutations affecting the eclipse rate (with phenotypes associated with reorganization of cellular metabolic processes and synthesis of viral budding precursor material) are highly favoured relative to all other traits examined. We find that mutations leading to reduced removal rates of infected cells (with phenotypes associated with innate immune evasion) have limited transmission advantage relative to mutations leading to humoral evasion. Predicted transmission probabilities, however, for mutations affecting innate immune evasion are more consistent with the range of clinically estimated household transmission probabilities for de novo mutations. This result suggests that although mutations affecting humoral evasion are more easily transmitted when they occur, mutations affecting innate immune evasion may occur more readily. We examine our predictions in the context of a number of previously characterized mutations in circulating strains of SARS-CoV-2. Our work offers both a null model for SARS-CoV-2 mutation rates and predicts which aspects of viral life history are most likely to successfully evolve, despite low mutation rates and repeated transmission bottlenecks.

Viral burden rebound in hospitalised patients with COVID-19 receiving oral antivirals in Hong Kong: a population-wide retrospective cohort study

BACKGROUND

Viral rebound after nirmatrelvir-ritonavir treatment has implications for the clinical management and isolation of patients with COVID-19. We evaluated an unselected, population-wide cohort to identify the incidence of viral burden rebound and associated risk factors and clinical outcomes.

METHODS

We did a retrospective cohort study of hospitalised patients with a confirmed diagnosis of COVID-19 in Hong Kong, China, for an observation period from Feb 26 to July 3, 2022 (during the omicron BA.2.2 variant wave). Adult patients (age ≥18 years) admitted 3 days before or after a positive COVID-19 test were selected from medical records held by the Hospital Authority of Hong Kong. We included patients with non-oxygen-dependent COVID-19 at baseline receiving either molnupiravir (800 mg twice a day for 5 days), nirmatrelvir-ritonavir (nirmatrelvir 300 mg with ritonavir 100 mg twice a day for 5 days), or no oral antiviral treatment (control group). Viral burden rebound was defined as a reduction in cycle threshold (Ct) value (≥3) on quantitative RT-PCR test between two consecutive measurements, with such decrease sustained in an immediately subsequent Ct measurement (for those patients with ≥3 Ct measurements). Logistic regression models were used to identify prognostic factors for viral burden rebound, and to assess associations between viral burden rebound and a composite clinical outcome of mortality, intensive care unit admission, and invasive mechanical ventilation initiation, stratified by treatment group.

FINDINGS

We included 4592 hospitalised patients with non-oxygen-dependent COVID-19 (1998 [43·5%] women and 2594 [56·5%] men). During the omicron BA.2.2 wave, viral burden rebound occurred in 16 of 242 patients (6·6% [95% CI 4·1-10·5]) receiving nirmatrelvir-ritonavir, 27 of 563 (4·8% [3·3-6·9]) receiving molnupiravir, and 170 of 3787 (4·5% [3·9-5·2]) in the control group. The incidence of viral burden rebound did not differ significantly across the three groups. Immunocompromised status was associated with increased odds of viral burden rebound, regardless of antiviral treatment (nirmatrelvir-ritonavir: odds ratio [OR] 7·37 [95% CI 2·56-21·26], p=0·0002; molnupiravir: 3·05 [1·28-7·25], p=0·012; control: 2·21 [1·50-3·27], p<0·0001). Among patients receiving nirmatrelvir-ritonavir, the odds of viral burden rebound were higher in those aged 18-65 years (vs >65 years; 3·09 [1·00-9·53], p=0·050), those with high comorbidity burden (score >6 on the Charlson Comorbidity Index; 6·02 [2·09-17·38], p=0·0009), and those concomitantly taking corticosteroids (7·51 [1·67-33·82], p=0·0086); whereas the odds were lower in those who were not fully vaccinated (0·16 [0·04-0·67], p=0·012). In patients receiving molnupiravir, those aged 18-65 years (2·68 [1·09-6·58], p=0·032) or on concomitant corticosteroids (3·11 [1·23-7·82], p=0·016) had increased odds of viral burden rebound. We found no association between viral burden rebound and occurrence of the composite clinical outcome from day 5 of follow-up (nirmatrelvir-ritonavir: adjusted OR 1·90 [0·48-7·59], p=0·36; molnupiravir: 1·05 [0·39-2·84], p=0·92; control: 1·27 [0·89-1·80], p=0·18).

INTERPRETATION

Viral burden rebound rates are similar between patients with antiviral treatment and those without. Importantly, viral burden rebound was not associated with adverse clinical outcomes.

FUNDING

Health and Medical Research Fund, Health Bureau, The Government of the Hong Kong Special Administrative Region, China.

TRANSLATION

For the Chinese translation of the abstract see Supplementary Materials section.

A quantitative systems pharmacology model of the pathophysiology and treatment of COVID-19 predicts optimal timing of pharmacological interventions

A quantitative systems pharmacology (QSP) model of the pathogenesis and treatment of SARS-CoV-2 infection can streamline and accelerate the development of novel medicines to treat COVID-19. Simulation of clinical trials allows in silico exploration of the uncertainties of clinical trial design and can rapidly inform their protocols. We previously published a preliminary model of the immune response to SARS-CoV-2 infection. To further our understanding of COVID-19 and treatment, we significantly updated the model by matching a curated dataset spanning viral load and immune responses in plasma and lung. We identified a population of parameter sets to generate heterogeneity in pathophysiology and treatment and tested this model against published reports from interventional SARS-CoV-2 targeting mAb and antiviral trials. Upon generation and selection of a virtual population, we match both the placebo and treated responses in viral load in these trials. We extended the model to predict the rate of hospitalization or death within a population. Via comparison of the in silico predictions with clinical data, we hypothesize that the immune response to virus is log-linear over a wide range of viral load. To validate this approach, we show the model matches a published subgroup analysis, sorted by baseline viral load, of patients treated with neutralizing Abs. By simulating intervention at different time points post infection, the model predicts efficacy is not sensitive to interventions within five days of symptom onset, but efficacy is dramatically reduced if more than five days pass post symptom onset prior to treatment.

A simple model of COVID-19 explains disease severity and the effect of treatments

Considerable effort has been made to better understand why some people suffer from severe COVID-19 while others remain asymptomatic. This has led to important clinical findings; people with severe COVID-19 generally experience persistently high levels of inflammation, slower viral load decay, display a dysregulated type-I interferon response, have less active natural killer cells and increased levels of neutrophil extracellular traps. How these findings are connected to the pathogenesis of COVID-19 remains unclear. We propose a mathematical model that sheds light on this issue by focusing on cells that trigger inflammation through molecular patterns: infected cells carrying pathogen-associated molecular patterns (PAMPs) and damaged cells producing damage-associated molecular patterns (DAMPs). The former signals the presence of pathogens while the latter signals danger such as hypoxia or lack of nutrients. Analyses show that SARS-CoV-2 infections can lead to a self-perpetuating feedback loop between DAMP expressing cells and inflammation, identifying the inability to quickly clear PAMPs and DAMPs as the main contributor to hyperinflammation. The model explains clinical findings and reveal conditions that can increase the likelihood of desired clinical outcome from treatment administration. In particular, the analysis suggest that antivirals need to be administered early during infection to have an impact on disease severity. The simplicity of the model and its high level of consistency with clinical findings motivate its use for the formulation of new treatment strategies.

Modeling explains prolonged SARS-CoV-2 nasal shedding relative to lung shedding in remdesivir treated rhesus macaques

In clinical trials, remdesivir decreased recovery time in hospitalized patients with SARS- CoV-2 and prevented hospitalization when given early during infection, despite not reducing nasal viral loads. In rhesus macaques, early remdesivir prevented pneumonia and lowered lung viral loads, but viral loads increased in nasal passages after five days. We developed mathematical models to explain these results. Our model raises the hypotheses that: 1) in contrast to nasal passages viral load monotonically decreases in lungs during therapy because of infection-dependent generation of refractory cells, 2) slight reduction in lung viral loads with an imperfect agent may result in a substantial decrease in lung damage, and 3) increases in nasal viral load may occur due to a blunting of peak viral load which decreases the intensity of the innate immune response. We demonstrate that a higher potency drug could lower viral loads in nasal passages and lung.