skip to main content


Title: Multiscale Model of Antiviral Timing, Potency, and Heterogeneity Effects on an Epithelial Tissue Patch Infected by SARS-CoV-2
We extend our established agent-based multiscale computational model of infection of lung tissue by SARS-CoV-2 to include pharmacokinetic and pharmacodynamic models of remdesivir. We model remdesivir treatment for COVID-19; however, our methods are general to other viral infections and antiviral therapies. We investigate the effects of drug potency, drug dosing frequency, treatment initiation delay, antiviral half-life, and variability in cellular uptake and metabolism of remdesivir and its active metabolite on treatment outcomes in a simulated patch of infected epithelial tissue. Non-spatial deterministic population models which treat all cells of a given class as identical can clarify how treatment dosage and timing influence treatment efficacy. However, they do not reveal how cell-to-cell variability affects treatment outcomes. Our simulations suggest that for a given treatment regime, including cell-to-cell variation in drug uptake, permeability and metabolism increase the likelihood of uncontrolled infection as the cells with the lowest internal levels of antiviral act as super-spreaders within the tissue. The model predicts substantial variability in infection outcomes between similar tissue patches for different treatment options. In models with cellular metabolic variability, antiviral doses have to be increased significantly (>50% depending on simulation parameters) to achieve the same treatment results as with the homogeneous cellular metabolism.  more » « less
Award ID(s):
1720625
NSF-PAR ID:
10359477
Author(s) / Creator(s):
; ; ; ; ; ; ;
Publisher / Repository:
MDPI
Date Published:
Journal Name:
Viruses
Volume:
14
Issue:
3
ISSN:
1999-4915
Page Range / eLocation ID:
605
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The antiviral remdesivir has been approved by regulatory bodies such as the European Medicines Agency (EMA) and the US Food and Drug administration (FDA) for the treatment of COVID-19. However, its efficacy is debated and toxicity concerns might limit the therapeutic range of this drug. Computational models that aid in balancing efficacy and toxicity would be of great help. Parametrizing models is difficult because the prodrug remdesivir is metabolized to its active form (RDV-TP) upon cell entry, which complicates dose–activity relationships. Here, we employ a computational model that allows drug efficacy predictions based on the binding affinity of RDV-TP for its target polymerase in SARS-CoV-2. We identify an optimal infusion rate to maximize remdesivir efficacy. We also assess drug efficacy in suppressing both wild-type and resistant strains, and thereby describe a drug regimen that may select for resistance. Our results differ from predictions using prodrug dose–response curves (pseudo-EC50s). We expect that reaching 90% inhibition (EC90) is insufficient to suppress SARS-CoV-2 in the lungs. While standard dosing mildly inhibits viral polymerase and therefore likely reduces morbidity, we also expect selection for resistant mutants for most realistic parameter ranges. To increase efficacy and safeguard against resistance, we recommend more clinical trials with dosing regimens that substantially increase the levels of RDV-TP and/or pair remdesivir with companion antivirals. 
    more » « less
  2. Abstract

    The rise of the novel virus SARS‐CoV2 which causes the disease known as COVID‐19 has led to a global pandemic claiming millions of lives. With no clinically approved treatment for COVID‐19, physicians initially struggled to treat the disease, and a need remains for improved antiviral therapies in this area. It is conceived early in the pandemic that an inhalable formulation of the drug remdesivir which directly targets the virus at the site of infection could improve therapeutic outcomes in COVID‐19. A set of requirements are developed that would be conducive to rapid drug approval: 1) try to use GRAS reagents 2) minimize excipient concentration and 3) achieve a working concentration of 5 mg/mL remdesivir to obtain a deliverable dose which is 5–10% of the IV dose. In this work, it is discovered that Poly(2‐oxazoline) block copolymers can stabilize drug nanocrystal suspensions and provide suitable formulation characteristics for aerosol delivery while maintaining antiviral efficacy. The authors believe POx block copolymers can be used as a semi‐ubiquitous stabilizer for the rapid development of nanocrystal formulations for new and existing diseases.

     
    more » « less
  3. Tumor microenvironment is a complex niche consisting of cancer cells and stromal cells in a network of extracellular matrix proteins and various soluble factors. Dynamic interactions among cellular and non-cellular components of the tumor microenvironment regulate tumor initiation and progression. Fibroblasts are the most abundant stromal cell type and dynamically interact with cancer cells both in primary tumors and in metastases. Cancer cells activate resident fibroblasts to produce and secrete soluble signaling molecules that support proliferation, migration, matrix invasion, and drug resistance of cancer cell and tumor angiogenesis. In recent years, various forms of three-dimensional tumor models have been developed to study tumor–stromal interactions and to identify anti-cancer drugs that block these interactions. There is currently a technological gap in development of tumor models that are physiologically relevant, scalable, and allow convenient, on-demand addition of desired components of the tumor microenvironment. In this review, we discuss three studies from our group that focus on developing bioengineered models to study tumor-stromal signaling. We will present these studies chronologically and based on their increasing complexity. We will discuss the validation of the models using a CXCL12-CXCR4 chemokine-receptor signaling present among activated fibroblasts and breast cancer cells in solid tumors, highlight the advantages and shortcomings of the models, and conclude with our perspectives on their applications. Impact statement Tumor stroma plays an important role in progression of cancers to a fatal metastatic disease. Modern treatment strategies are considering targeting tumor stroma to improve outcomes for cancer patients. A current challenge to develop stroma-targeting therapeutics is the lack of preclinical physiologic tumor models. Animal models widely used in cancer research lack human stroma and are not amenable to screening of chemical compounds for cancer drug discovery. In this review, we outline in vitro three-dimensional tumor models that we have developed to study the interactions among cancer cells and stromal cells. We describe development of the tumor models in a modular fashion, from a spheroid model to a sophisticated organotypic model, and discuss the importance of using correct physiologic models to recapitulate tumor-stromal signaling. These biomimetic tumor models will facilitate understanding of tumor-stromal signaling biology and provide a scalable approach for testing and discovery of cancer drugs. 
    more » « less
  4. Faeder, James R. (Ed.)
    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. 
    more » « less
  5. null (Ed.)
    Mitochondria are dynamic organelles that undergo fission and fusion. While they are essential for cellular metabolism, the effect of dysregulated mitochondrial dynamics on cellular metabolism is not fully understood. We previously found that transmembrane protein 135 ( Tmem135) plays a role in the regulation of mitochondrial dynamics in mice. Mice homozygous for a Tmem135 mutation ( Tmem135 FUN025/FUN025 ) display accelerated aging and age-related disease pathologies in the retina including the retinal pigment epithelium (RPE). We also generated a transgenic mouse line globally overexpressing the Tmem135 gene ( Tmem135 TG). In several tissues and cells that we studied such as the retina, heart, and fibroblast cells, we observed that the Tmem135 mutation causes elongated mitochondria, while overexpression of Tmem135 results in fragmented mitochondria. To investigate how abnormal mitochondrial dynamics affect metabolic signatures of tissues and cells, we identified metabolic changes in primary RPE cell cultures as well as heart, cerebellum, and hippocampus isolated from Tmem135 FUN025/FUN025 mice (fusion > fission) and Tmem135 TG mice (fusion < fission) using nuclear magnetic resonance spectroscopy. Metabolomics analysis revealed a tissue-dependent response to Tmem135 alterations, whereby significant metabolic changes were observed in the heart of both Tmem135 mutant and TG mice as compared to wild-type, while negligible effects were observed in the cerebellum and hippocampus. We also observed changes in Tmem135 FUN025/FUN025 and Tmem135 TG RPE cells associated with osmosis and glucose and phospholipid metabolism. We observed depletion of NAD + in both Tmem135 FUN025/FUN025 and Tmem135 TG RPE cells, indicating that imbalance in mitochondrial dynamics to both directions lowers the cellular NAD + level. Metabolic changes identified in this study might be associated with imbalanced mitochondrial dynamics in heart tissue and RPE cells which can likely lead to functional abnormalities. Impact statement Mitochondria are dynamic organelles undergoing fission and fusion. Proper regulation of this process is important for healthy aging process, as aberrant mitochondrial dynamics are associated with several age-related diseases/pathologies. However, it is not well understood how imbalanced mitochondrial dynamics may lead to those diseases and pathologies. Here, we aimed to determine metabolic alterations in tissues and cells from mouse models with over-fused (fusion > fission) and over-fragmented (fusion < fission) mitochondria that display age-related disease pathologies. Our results indicated tissue-dependent sensitivity to these mitochondrial changes, and metabolic pathways likely affected by aberrant mitochondrial dynamics. This study provides new insights into how dysregulated mitochondrial dynamics could lead to functional abnormalities of tissues and cells. 
    more » « less