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 withmore »
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 more »
- Award ID(s):
- 1720625
- Publication Date:
- NSF-PAR ID:
- 10359477
- Journal Name:
- Viruses
- Volume:
- 14
- Issue:
- 3
- Page Range or eLocation-ID:
- 605
- ISSN:
- 1999-4915
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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 advantagesmore »
-
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 effectsmore »
-
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 cannotmore »
-
Optimization-based models have been used to predict cellular behavior for over 25 years. The constraints in these models are derived from genome annotations, measured macromolecular composition of cells, and by measuring the cell's growth rate and metabolism in different conditions. The cellular goal (the optimization problem that the cell is trying to solve) can be challenging to derive experimentally for many organisms, including human or mammalian cells, which have complex metabolic capabilities and are not well understood. Existing approaches to learning goals from data include (a) estimating a linear objective function, or (b) estimating linear constraints that model complex biochemical reactions and constrain the cell's operation. The latter approach is important because often the known reactions are not enough to explain observations; therefore, there is a need to extend automatically the model complexity by learning new reactions. However, this leads to nonconvex optimization problems, and existing tools cannot scale to realistically large metabolic models. Hence, constraint estimation is still used sparingly despite its benefits for modeling cell metabolism, which is important for developing novel antimicrobials against pathogens, discovering cancer drug targets, and producing value-added chemicals. Here, we develop the first approach to estimating constraint reactions from data that can scalemore »