Timing of Antiviral Treatment Initiation is Critical to Reduce SARS-CoV-2 Viral Load

Modeling of SARS-CoV-2 Treatment Effects for Informed Drug Repurposing

Several repurposed drugs are currently under investigation in the fight against coronavirus disease 2019 (COVID-19). Candidates are often selected solely by their effective concentrations in vitro, an approach that has largely not lived up to expectations in COVID-19. Cell lines used in in vitro experiments are not necessarily representative of lung tissue. Yet, even if the proposed mode of action is indeed true, viral dynamics in vivo, host response, and concentration-time profiles must also be considered. Here we address the latter issue and describe a model of human SARS-CoV-2 viral kinetics with acquired immune response to investigate the dynamic impact of timing and dosing regimens of hydroxychloroquine, lopinavir/ritonavir, ivermectin, artemisinin, and nitazoxanide. We observed greatest benefits when treatments were given immediately at the time of diagnosis. Even interventions with minor antiviral effect may reduce host exposure if timed correctly. Ivermectin seems to be at least partially effective: given on positivity, peak viral load dropped by 0.3-0.6 log units and exposure by 8.8-22.3%. The other drugs had little to no appreciable effect. Given how well previous clinical trial results for hydroxychloroquine and lopinavir/ritonavir are explained by the models presented here, similar strategies should be considered in future drug candidate prioritization efforts.

Mechanistic Modeling of SARS-CoV-2 and Other Infectious Diseases and the Effects of Therapeutics

Modern viral kinetic modeling and its application to therapeutics is a field that attracted the attention of the medical, pharmaceutical, and modeling communities during the early days of the AIDS epidemic. Its successes led to applications of modeling methods not only to HIV but a plethora of other viruses, such as hepatitis C virus (HCV), hepatitis B virus and cytomegalovirus, which along with HIV cause chronic diseases, and viruses such as influenza, respiratory syncytial virus, West Nile virus, Zika virus, and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which generally cause acute infections. Here we first review the historical development of mathematical models to understand HIV and HCV infections and the effects of treatment by fitting the models to clinical data. We then focus on recent efforts and contributions of applying these models towards understanding SARS-CoV-2 infection and highlight outstanding questions where modeling can provide crucial insights and help to optimize nonpharmaceutical and pharmaceutical interventions of the coronavirus disease 2019 (COVID-19) pandemic. The review is written from our personal perspective emphasizing the power of simple target cell limited models that provided important insights and then their evolution into more complex models that captured more of the virology and immunology. To quote Albert Einstein, "Everything should be made as simple as possible, but not simpler," and this idea underlies the modeling we describe below.

Combination of Hydroxychloroquine Plus Azithromycin As Potential Treatment for COVID-19 Patients: Safety Profile, Drug Interactions, and Management of Toxicity

The coronavirus disease 2019 (COVID-19), caused by infection with severe acute respiratory syndrome coronavirus 2, has recently emerged worldwide. In this context, there is an urgent need to identify safe and effective therapeutic strategies for treatment of such highly contagious disease. We recently reported promising results of combining hydroxychloroquine and azithromycin as an early treatment option. Although ongoing clinical trials are challenging the efficacy of this combination, many clinicians claim the authorization to or have already begun to use it to treat COVID-19 patients worldwide. The aim of this article is to share pharmacology considerations contributing to the rationale of this combination, and to provide safety information to prevent toxicity and drug-drug interactions, based on available evidence.

A quantitative model used to compare within-host SARS-CoV-2, MERS-CoV, and SARS-CoV dynamics provides insights into the pathogenesis and treatment of SARS-CoV-2

The scientific community is focused on developing antiviral therapies to mitigate the impacts of the ongoing novel coronavirus disease 2019 (COVID-19) outbreak. This will be facilitated by improved understanding of viral dynamics within infected hosts. Here, using a mathematical model in combination with published viral load data, we compare within-host viral dynamics of SARS-CoV-2 with analogous dynamics of MERS-CoV and SARS-CoV. Our quantitative analyses using a mathematical model revealed that the within-host reproduction number at symptom onset of SARS-CoV-2 was statistically significantly larger than that of MERS-CoV and similar to that of SARS-CoV. In addition, the time from symptom onset to the viral load peak for SARS-CoV-2 infection was shorter than those of MERS-CoV and SARS-CoV. These findings suggest the difficulty of controlling SARS-CoV-2 infection by antivirals. We further used the viral dynamics model to predict the efficacy of potential antiviral drugs that have different modes of action. The efficacy was measured by the reduction in the viral load area under the curve (AUC). Our results indicate that therapies that block de novo infection or virus production are likely to be effective if and only if initiated before the viral load peak (which appears 2-3 days after symptom onset), but therapies that promote cytotoxicity of infected cells are likely to have effects with less sensitivity to the timing of treatment initiation. Furthermore, combining a therapy that promotes cytotoxicity and one that blocks de novo infection or virus production synergistically reduces the AUC with early treatment. Our unique modeling approach provides insights into the pathogenesis of SARS-CoV-2 and may be useful for development of antiviral therapies.

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.

Modeling SARS-CoV-2 viral kinetics and association with mortality in hospitalized patients from the French COVID cohort

The characterization of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) viral kinetics in hospitalized patients and its association with mortality is unknown. We analyzed death and nasopharyngeal viral kinetics in 655 hospitalized patients from the prospective French COVID cohort. The model predicted a median peak viral load that coincided with symptom onset. Patients with age ≥65 y had a smaller loss rate of infected cells, leading to a delayed median time to viral clearance occurring 16 d after symptom onset as compared to 13 d in younger patients (P < 10-4). In multivariate analysis, the risk factors associated with mortality were age ≥65 y, male gender, and presence of chronic pulmonary disease (hazard ratio [HR] > 2.0). Using a joint model, viral dynamics after hospital admission was an independent predictor of mortality (HR = 1.31, P < 10-3). Finally, we used our model to simulate the effects of effective pharmacological interventions on time to viral clearance and mortality. A treatment able to reduce viral production by 90% upon hospital admission would shorten the time to viral clearance by 2.0 and 2.9 d in patients of age <65 y and ≥65 y, respectively. Assuming that the association between viral dynamics and mortality would remain similar to that observed in our population, this could translate into a reduction of mortality from 19 to 14% in patients of age ≥65 y with risk factors. Our results show that viral dynamics is associated with mortality in hospitalized patients. Strategies aiming to reduce viral load could have an effect on mortality rate in this population.

Model-informed drug repurposing: A pharmacometric approach to novel pathogen preparedness, response and retrospection

During a pandemic caused by a novel pathogen (NP), drug repurposing offers the potential of a rapid treatment response via a repurposed drug (RD) while more targeted treatments are developed. Five steps of model‐informed drug repurposing (MIDR) are discussed: (i) utilize RD product label and in vitro NP data to determine initial proof of potential, (ii) optimize potential posology using clinical pharmacokinetics (PK) considering both efficacy and safety, (iii) link events in the viral life cycle to RD PK, (iv) link RD PK to clinical and virologic outcomes, and optimize clinical trial design, and (v) assess RD treatment effects from trials using model‐based meta‐analysis. Activities which fall under these five steps are categorized into three stages: what can be accomplished prior to an NP emergence (preparatory stage), during the NP pandemic (responsive stage) and once the crisis has subsided (retrospective stage). MIDR allows for extraction of a greater amount of information from emerging data and integration of disparate data into actionable insight.

In silico dynamics of COVID-19 phenotypes for optimizing clinical management

Understanding the underlying mechanisms of COVID-19 progression and the impact of various pharmaceutical interventions is crucial for the clinical management of the disease. We developed a comprehensive mathematical framework based on the known mechanisms of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, incorporating the renin-angiotensin system and ACE2, which the virus exploits for cellular entry, key elements of the innate and adaptive immune responses, the role of inflammatory cytokines, and the coagulation cascade for thrombus formation. The model predicts the evolution of viral load, immune cells, cytokines, thrombosis, and oxygen saturation based on patient baseline condition and the presence of comorbidities. Model predictions were validated with clinical data from healthy people and COVID-19 patients, and the results were used to gain insight into identified risk factors of disease progression including older age; comorbidities such as obesity, diabetes, and hypertension; and dysregulated immune response. We then simulated treatment with various drug classes to identify optimal therapeutic protocols. We found that the outcome of any treatment depends on the sustained response rate of activated CD8+ T cells and sufficient control of the innate immune response. Furthermore, the best treatment-or combination of treatments-depends on the preinfection health status of the patient. Our mathematical framework provides important insight into SARS-CoV-2 pathogenesis and could be used as the basis for personalized, optimal management of COVID-19.

Using in silico viral kinetic models to guide therapeutic strategies during a pandemic: An example in SARS-CoV-2

AIMS

We propose the use of in silico mathematical models to provide insights that optimize therapeutic interventions designed to effectively treat respiratory infection during a pandemic. A modelling and simulation framework is provided using SARS-CoV-2 as an example, considering applications for both treatment and prophylaxis.

METHODS

A target cell-limited model was used to quantify the viral infection dynamics of SARS-CoV-2 in a pooled population of 105 infected patients. Parameter estimates from the resulting model were used to simulate and compare the impact of various interventions against meaningful viral load endpoints.

RESULTS

Robust parameter estimates were obtained for the basic reproduction number, viral release rate and infected-cell mortality from the infection model. These estimates were informed by the largest dataset currently available for SARS-CoV-2 viral time course. The utility of this model was demonstrated using simulations, which hypothetically introduced inhibitory or stimulatory drug mechanisms at various target sites within the viral life-cycle. We show that early intervention is crucial to achieving therapeutic benefit when monotherapy is administered. In contrast, combination regimens of two or three drugs may provide improved outcomes if treatment is initiated late. The latter is relevant to SARS-CoV-2, where the period between infection and symptom onset is relatively long.

CONCLUSIONS

The use of in silico models can provide viral load predictions that can rationalize therapeutic strategies against an emerging viral pathogen.

Dynamical characterization of antiviral effects in COVID-19

Mathematical models describing SARS-CoV-2 dynamics and the corresponding immune responses in patients with COVID-19 can be critical to evaluate possible clinical outcomes of antiviral treatments. In this work, based on the concept of virus spreadability in the host, antiviral effectiveness thresholds are determined to establish whether or not a treatment will be able to clear the infection. In addition, the virus dynamic in the host - including the time-to-peak and the final monotonically decreasing behavior - is characterized as a function of the time to treatment initiation. Simulation results, based on nine patient data, show the potential clinical benefits of a treatment classification according to patient critical parameters. This study is aimed at paving the way for the different antivirals being developed to tackle SARS-CoV-2.

Randomized controlled open label trial on the use of favipiravir combined with inhaled interferon beta-1b in hospitalized patients with moderate to severe COVID-19 pneumonia

OBJECTIVE

To evaluate the therapeutic effectiveness of favipiravir combined with inhaled interferon beta-1b in adult patients hospitalized with moderate to severe COVID-19 pneumonia.

METHODS

A randomized, open-label controlled trial of oral favipiravir in adults hospitalized with moderate to severe COVID-19 pneumonia from June 22nd 2020 to August 13th 2020 was conducted. Patients were randomly assigned to receive either a combination of favipiravir with interferon beta-1b by inhalation aerosol or hydroxychloroquine (HCQ). The outcome endpoints included improvement in inflammatory markers, lower length of hospital stay (LOS), discharges and lower overall 14-day mortality.

RESULTS

A total of 89 patients underwent randomization with 49% (n = 44) assigned to favipiravir and 51% (n = 45) assigned HCQ. The overall mean age was 55 ± 14 years and 58% (n = 52) were males. There were no significant differences in the inflammatory biomarkers at hospital discharge between the two groups; C-reactive protein (p = 0.413), ferritin (p = 0.968), lactate dehydrogenase (p = 0.259) and interleukin 6 (p = 0.410). There were also no significant differences between the two groups with regards to the overall LOS (7 vs 7 days; p = 0.948), transfers to the ICU (18.2% vs 17.8%; p = 0.960), discharges (65.9% vs 68.9%; p = 0.764) and overall mortality (11.4% vs 13.3%; p = 0.778).

CONCLUSIONS

No differences in clinical outcomes were found between favipiravir plus inhaled interferon beta-1b and hydroxychloroquine in adults hospitalized with moderate to severe COVID-19 pneumonia.

Targeting TMPRSS2 and Cathepsin B/L together may be synergistic against SARS-CoV-2 infection

The entry of SARS-CoV-2 into target cells requires the activation of its surface spike protein, S, by host proteases. The host serine protease TMPRSS2 and cysteine proteases Cathepsin B/L can activate S, making two independent entry pathways accessible to SARS-CoV-2. Blocking the proteases prevents SARS-CoV-2 entry in vitro. This blockade may be achieved in vivo through 'repurposing' drugs, a potential treatment option for COVID-19 that is now in clinical trials. Here, we found, surprisingly, that drugs targeting the two pathways, although independent, could display strong synergy in blocking virus entry. We predicted this synergy first using a mathematical model of SARS-CoV-2 entry and dynamics in vitro. The model considered the two pathways explicitly, let the entry efficiency through a pathway depend on the corresponding protease expression level, which varied across cells, and let inhibitors compromise the efficiency in a dose-dependent manner. The synergy predicted was novel and arose from effects of the drugs at both the single cell and the cell population levels. Validating our predictions, available in vitro data on SARS-CoV-2 and SARS-CoV entry displayed this synergy. Further, analysing the data using our model, we estimated the relative usage of the two pathways and found it to vary widely across cell lines, suggesting that targeting both pathways in vivo may be important and synergistic given the broad tissue tropism of SARS-CoV-2. Our findings provide insights into SARS-CoV-2 entry into target cells and may help improve the deployability of drug combinations targeting host proteases required for the entry.

Antagonism of Type I Interferon by Severe Acute Respiratory Syndrome Coronavirus 2

The coronavirus disease 2019 (COVID-19) pandemic is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), warranting urgent study of the molecular mechanisms of SARS-CoV-2 infection and host immune response. Type I interferon (IFN-I) is a key component of host innate immune system responsible for eliminating the virus at the early stage of infection. In contrast, SARS-CoV-2 has evolved multiple strategies to evade innate immune response to facilitate viral replication, transmission, and pathogenesis. This review summarizes the recent progresses on SARS-CoV-2 proteins that antagonize host IFN-I production and/or signaling. These progresses have provided knowledge for new vaccine and antiviral development to prevent and control COVID-19.

A modular framework for multiscale, multicellular, spatiotemporal modeling of acute primary viral infection and immune response in epithelial tissues and its application to drug therapy timing and effectiveness

Simulations of tissue-specific effects of primary acute viral infections like COVID-19 are essential for understanding disease outcomes and optimizing therapies. Such simulations need to support continuous updating in response to rapid advances in understanding of infection mechanisms, and parallel development of components by multiple groups. We present an open-source platform for multiscale spatiotemporal simulation of an epithelial tissue, viral infection, cellular immune response and tissue damage, specifically designed to be modular and extensible to support continuous updating and parallel development. The base simulation of a simplified patch of epithelial tissue and immune response exhibits distinct patterns of infection dynamics from widespread infection, to recurrence, to clearance. Slower viral internalization and faster immune-cell recruitment slow infection and promote containment. Because antiviral drugs can have side effects and show reduced clinical effectiveness when given later during infection, we studied the effects on progression of treatment potency and time-of-first treatment after infection. In simulations, even a low potency therapy with a drug which reduces the replication rate of viral RNA greatly decreases the total tissue damage and virus burden when given near the beginning of infection. Many combinations of dosage and treatment time lead to stochastic outcomes, with some simulation replicas showing clearance or control (treatment success), while others show rapid infection of all epithelial cells (treatment failure). Thus, while a high potency therapy usually is less effective when given later, treatments at late times are occasionally effective. We illustrate how to extend the platform to model specific virus types (e.g., hepatitis C) and add additional cellular mechanisms (tissue recovery and variable cell susceptibility to infection), using our software modules and publicly-available software repository.

Disentangling the dynamical underpinnings of differences in SARS-CoV-2 pathology using within-host ecological models

Health outcomes following infection with Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) are remarkably variable. The way the virus spreads inside hosts, and how this spread interacts with host immunity and physiology, is likely to determine variation in health outcomes. Decades of data and dynamical analyses of how other viruses spread and interact with host cells could shed light on SARS-CoV-2 within-host trajectories. We review how common axes of variation in within-host dynamics and emergent pathology (such as age and sex) might be combined with ecological principles to understand the case of SARS-CoV-2. We highlight pitfalls in application of existing theoretical frameworks relevant to the complexity of the within-host context and frame the discussion in terms of growing knowledge of the biology of SARS-CoV-2. Viewing health outcomes for SARS-CoV-2 through the lens of ecological models underscores the value of repeated measures on individuals, especially since many lines of evidence suggest important contingence on trajectory.

Pharmacokinetic Basis of the Hydroxychloroquine Response in COVID-19: Implications for Therapy and Prevention

BACKGROUND AND OBJECTIVES

Chloroquine/hydroxychloroquine has recently been the subject of intense debate regarding its potential antiviral activity against SARS-Cov-2, the etiologic agent of COVID-19. Some report possible curative effects; others do not. Therefore, the objective of this study was to simulate possible scenarios of response to hydroxychloroquine in COVID-19 patients using mathematical modeling.

METHODS

To shed some light on this controversial topic, we simulated hydroxychloroquine-based interventions on virus/host cell dynamics using a basic system of previously published differential equations. Mathematical modeling was implemented using Python programming language v 3.7.

RESULTS

According to mathematical modeling, hydroxychloroquine may have an impact on the amplitude of the viral load peak and viral clearance if the drug is administered early enough (i.e., when the virus is still confined within the pharyngeal cavity). The effects of chloroquine/hydroxychloroquine may be fully explained only when also considering the capacity of this drug to increase the death rate of SARS-CoV-2-infected cells, in this case by enhancing the cell-mediated immune response.

CONCLUSIONS

These considerations may not only be applied to chloroquine/hydroxychloroquine but may have more general implications for development of anti-COVID-19 combination therapies and prevention strategies through an increased death rate of the infected cells.

Quantifying the effect of remdesivir in rhesus macaques infected with SARS-CoV-2

The world is in the midst of a pandemic caused by a novel coronavirus and is desperately searching for possible treatments. The antiviral remdesivir has shown some effectiveness against SARS-CoV-2 in vitro and in a recent animal study. We use data from a study of remdesivir in rhesus macaques to fit a viral kinetics model in an effort to determine the most appropriate mathematical descripton of the effect of remdesivir. We find statistically significant differences in the viral decay rate and use this to inform a possible mathematical formulation of the effect of remdesivir. Unfortunately, this model formulation suggests that the application of remdesivir will lengthen SARS-CoV-2 infections, putting into question its potential clinical benefit.

Drug repurposing and cytokine management in response to COVID-19: A review

Coronavirus disease 2019 (COVID-19), the infectious disease caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is an aggressive disease that attacks the respiratory tract and has a higher fatality rate than seasonal influenza. The COVID-19 pandemic is a global health crisis, and no specific therapy or drug has been formally recommended for use against SARS-CoV-2 infection. In this context, it is a rational strategy to investigate the repurposing of existing drugs to use in the treatment of COVID-19 patients. In the meantime, the medical community is trialing several therapies that target various antiviral and immunomodulating mechanisms to use against the infection. There is no doubt that antiviral and supportive treatments are important in the treatment of COVID-19 patients, but anti-inflammatory therapy also plays a pivotal role in the management COVID-19 patients due to its ability to prevent further injury and organ damage or failure. In this review, we identified drugs that could modulate cytokines levels and play a part in the management of COVID-19. Several drugs that possess an anti-inflammatory profile in others illnesses have been studied in respect of their potential utility in the treatment of the hyperinflammation induced by SAR-COV-2 infection. We highlight a number of antivirals, anti-rheumatic, anti-inflammatory, antineoplastic and antiparasitic drugs that have been found to mitigate cytokine production and consequently attenuate the "cytokine storm" induced by SARS-CoV-2. Reduced hyperinflammation can attenuate multiple organ failure, and even reduce the mortality associated with severe COVID-19. In this context, despite their current unproven clinical efficacy in relation to the current pandemic, the repurposing of drugs with anti-inflammatory activity to use in the treatment of COVID-19 has become a topic of great interest.

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.

Initial Inoculum and the Severity of COVID-19: A Mathematical Modeling Study of the Dose-Response of SARS-CoV-2 Infections

SARS-CoV-2 (Severe acute respiratory syndrome coronavirus 2) causes a variety of responses in those who contract the virus, ranging from asymptomatic infections to acute respiratory failure and death. While there are likely multiple mechanisms triggering severe disease, one potential cause of severe disease is the size of the initial inoculum. For other respiratory diseases, larger initial doses lead to more severe outcomes. We investigate whether there is a similar link for SARS-CoV-2 infections using the combination of an agent-based model (ABM) and a partial differential equation model (PDM). We use the model to examine the viral time course for different sizes of initial inocula, generating dose-response curves for peak viral load, time of viral peak, viral growth rate, infection duration, and area under the viral titer curve. We find that large initial inocula lead to short infections, but with higher viral titer peaks; and that smaller initial inocula lower the viral titer peak, but make the infection last longer.

A modular framework for multiscale, multicellular, spatiotemporal modeling of acute primary viral infection and immune response in epithelial tissues and its application to drug therapy timing and effectiveness: A multiscale model of viral infection in epithelial tissues

Simulations of tissue-specific effects of primary acute viral infections like COVID-19 are essential for understanding disease outcomes and optimizing therapies. Such simulations need to support continuous updating in response to rapid advances in understanding of infection mechanisms, and parallel development of components by multiple groups. We present an open-source platform for multiscale spatiotemporal simulation of an epithelial tissue, viral infection, cellular immune response and tissue damage, specifically designed to be modular and extensible to support continuous updating and parallel development. The base simulation of a simplified patch of epithelial tissue and immune response exhibits distinct patterns of infection dynamics from widespread infection, to recurrence, to clearance. Slower viral internalization and faster immune-cell recruitment slow infection and promote containment. Because antiviral drugs can have side effects and show reduced clinical effectiveness when given later during infection, we studied the effects on progression of treatment potency and time-of-first treatment after infection. In simulations, even a low potency therapy with a drug which reduces the replication rate of viral RNA greatly decreases the total tissue damage and virus burden when given near the beginning of infection. Many combinations of dosage and treatment time lead to stochastic outcomes, with some simulation replicas showing clearance or control (treatment success), while others show rapid infection of all epithelial cells (treatment failure). Thus, while a high potency therapy usually is less effective when given later, treatments at late times are occasionally effective. We illustrate how to extend the platform to model specific virus types (e.g., hepatitis C) and add additional cellular mechanisms (tissue recovery and variable cell susceptibility to infection), using our software modules and publicly-available software repository.

In silico dynamics of COVID-19 phenotypes for optimizing clinical management

Understanding the underlying mechanisms of COVID-19 progression and the impact of various pharmaceutical interventions is crucial for the clinical management of the disease. We developed a comprehensive mathematical framework based on the known mechanisms of the SARS-CoV-2 virus infection, incorporating the renin-angiotensin system and ACE2, which the virus exploits for cellular entry, key elements of the innate and adaptive immune responses, the role of inflammatory cytokines and the coagulation cascade for thrombus formation. The model predicts the evolution of viral load, immune cells, cytokines, thrombosis, and oxygen saturation based on patient baseline condition and the presence of co-morbidities. Model predictions were validated with clinical data from healthy people and COVID-19 patients, and the results were used to gain insight into identified risk factors of disease progression including older age, co-morbidities such as obesity, diabetes, and hypertension, and dysregulated immune response 1,2 . We then simulated treatment with various drug classes to identify optimal therapeutic protocols. We found that the outcome of any treatment depends on the sustained response rate of activated CD8 + T cells and sufficient control of the innate immune response. Furthermore, the best treatment -or combination of treatments - depends on the pre-infection health status of the patient. Our mathematical framework provides important insight into SARS-CoV-2 pathogenesis and could be used as the basis for personalized, optimal management of COVID-19.

A pharmacological perspective of chloroquine in SARS-CoV-2 infection: An old drug for the fight against a new coronavirus?

The coronavirus disease 2019 (COVID-19) pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is having serious consequences on health and the economy worldwide. All evidence-based treatment strategies need to be considered to combat this new virus. Drugs need to be considered on scientific grounds of efficacy, safety and cost. Chloroquine (CQ) and hydroxychloroquine (HCQ) are old drugs used in the treatment of malaria. Moreover, their antiviral properties have been previously studied, including against coronaviruses, where evidence of efficacy has been found. In the current race against time triggered by the COVID-19 pandemic, the search for new antivirals is very important. However, consideration should be given to old drugs with known anti-coronavirus activity, such as CQ and HCQ. These could be integrated into current treatment strategies while novel treatments are awaited, also in light of the fact that they display an anticoagulant effect that facilitates the activity of low-molecular-weight heparin, aimed at preventing acute respiratory distress syndrome (ARDS)-associated thrombotic events. The safety of CQ and HCQ has been studied for over 50 years, however recently published data raise concerns for cardiac toxicity of CQ/HCQ in patients with COVID-19. This review also re-examines the real information provided by some of the published alarming reports, although concluding that cardiac toxicity should in any case be stringently monitored in patients receiving CQ/HCQ.

Predictions of Systemic, Intracellular, and Lung Concentrations of Azithromycin With Different Dosing Regimens Used in COVID-19 Clinical Trials

Azithromycin (AZ), a broad-spectrum macrolide antibiotic, is being investigated in patients with coronavirus disease 2019 (COVID-19). A population pharmacokinetic model was implemented to predict lung, intracellular poly/mononuclear cell (peripheral blood monocyte (PBM)/polymorphonuclear leukocyte (PML)), and alveolar macrophage (AM) concentrations using published data and compared against preclinical effective concentration 90% (EC₉₀ ) for severe acute respiratory syndrome-coronavirus 2 (SARS-CoV-2). The final model described the data reported in eight publications adequately. Consistent with its known properties, concentrations were higher in AM and PBM/PML, followed by lung tissue, and lowest systemically. Simulated PBM/PML concentrations exceeded EC₉₀ following the first dose and for ~ 14 days following 500 mg q.d. for 3 days or 500 mg q.d. for 1 day/250 mg q.d. on days 2-5, 10 days following a single 1,000 mg dose, and for > 20 days with 500 mg q.d. for 10 days. AM concentrations exceeded the 90% inhibitory concentration for > 20 days for all regimens. These data will better inform optimization of dosing regimens for AZ clinical trials.

Time to Step Up: A Call to Action For the Clinical and Quantitative Pharmacology Community to Accelerate Therapeutics for COVID-19

The global response to finding therapeutics for coronavirus disease 2019 (COVID‐19) is chaotic even if well intentioned. There is an opportunity, but more importantly, an obligation for the global clinical and quantitative pharmacology community to come together and use our state‐of‐the‐art tools and expertise to help society accelerate therapeutics to fight COVID‐19. This brief commentary is a call to action and highlights how the global pharmacology community should contribute to the COVID‐19 pandemic and prepare for future pandemics.