Dataset Information

A modular framework for multiscale multicellular spatial modeling of viral infection, immune response and drug therapy timing and efficacy in epithelial tissues: A multiscale model of viral infection in epithelial tissues.

ABSTRACT: Development of predictive quantitative models of all aspects of COVID-19 is essential for rapidly understanding the causes of differing disease outcomes and vulnerabilities, suggesting drug and therapeutic targets, and designing optimized personalized interventions. Easy to implement, predictive multiscale modeling frameworks to integrate the wide variety of clinical and research datasets into actionable insights, which could inform therapeutic regime strategies are lacking. We present a multiscale, multicellular, spatiotemporal model of the infection of epithelial tissue by a generic virus, a simplified cellular immune response and viral and immune-induced tissue damage. Our initial model is built of modular components to allow it to be easily extended and adapted in a collaborative fashion to describe specific viral infections, tissue types and immune responses. The model allows us to define three parameter regimes: where viral infection coincides with a massive cytopathic effect, where the immune System rapidly controls the virus and where the immune System controls the virus but extensive tissue damage occurs. We use the model in a proof-of-concept application to evaluate a number of drug therapy concepts. Inhibition of viral internalization and faster immune-cell recruitment lead to containment of infection. Fast viral internalization and slower immune response lead to uncontrolled spread of infection. Simulation of a drug, whose mode of action is to reduce production of viral RNAs, shows that a relatively limited reduction of viral replication at the beginning of infection greatly decreases the total area of tissue damage and maximal viral load, while even a treatment that greatly reduces the rate of genomic replication rapidly loses efficacy as the infection progresses. A number of simulation conditions lead to stochastically variable outcomes, with some replicas clearing or controlling the virus, while others see virus-induced damage sweep the simulated lung patch. The model is open-source and modular, allowing rapid development and extension of its components by groups working in parallel.

SUBMITTER: Sego TJ

PROVIDER: S-EPMC7263495 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:We designed a strategy for microarray analysis that is based on the identification of transcriptional modules formed by genes coordinately expressed in multiple disease data sets. Mapping changes in gene expression at the module level generated disease-specific transcriptional fingerprints that provide a stable framework for the visualization and functional interpretation of microarray data. The first step of the module-construction process analyzes expression patterns of transcripts across samples for individual diseases: sets of coordinately expressed transcripts were identified with an unsupervised clustering algorithm; in this case, the GeneSpring Version 7.1 (Agilent) implementation of the K-Means algorithm (k = 30). All transcripts detected in at least one sample were used as input; no screening for differential expression was performed. The second step of the module-construction process analyzed the “clustering behavior” of transcripts across diseases, taking into account the possibility that genes may cocluster in some diseases but not in others. Also, in our example, the transcripts that clustered together across all eight diseases were grouped to form a set of modules (round 1 of selection), and the stringency of the analysis was then decreased gradually to identify transcripts that belong to a similar K-means cluster in only a subset of diseases (round 2: seven out of eight diseases; round 3: six out of eight diseases). It is important to note that the module-selection process is “data-driven” and does not involve manual selection of genes by the investigator. We implemented the module-construction strategy described above, using as input a total of 239 peripheral-blood mononuclear cell (PBMC) samples obtained from individuals with one of the following conditions: systemic juvenile idiopathic arthritis (n = 47), systemic lupus erythematosus (n = 40), type I diabetes (n = 20), metastatic melanoma (n = 39), acute infections (Escherichia coli [n = 22], Staphylococcus aureus [n = 18], Influenza A [n = 16]) or liver-transplant recipients undergoing immunosuppressive therapy (n = 37). Transcriptional profiles were generated with Affymetrix U133A and U133B GeneChips (> 44,000 probe sets). A total of 4742 transcripts, distributed among 28 sets, were selected after running of the module-construction algorithm described above. Each module is assigned a unique identifier indicating the round and order of selection (i.e., M3.1 is the first module identified in the third round of selection). The stringency of this algorithm was tested statistically by implementation of the same module-construction procedure after randomization of the original data set. This process was repeated 200 times, without a single module identified. Therefore, the analysis of gene-cluster membership across multiple diseases provided a stringent means to identify PBMC transcriptional modules.

Project description:Mathematical modeling provides the predictive ability to understand the metabolic reprogramming and complex pathways that mediate cancer cells' proliferation. We present a mathematical model using a multiscale, multicellular approach to simulate avascular tumor growth, applied to pancreatic cancer. The model spans three distinct spatial and temporal scales. At the extracellular level, reaction diffusion equations describe nutrient concentrations over a span of seconds. At the cellular level, a lattice-based energy driven stochastic approach describes cellular phenomena including adhesion, proliferation, viability and cell state transitions, occurring on the timescale of hours. At the sub-cellular level, we incorporate a detailed kinetic model of intracellular metabolite dynamics on the timescale of minutes, which enables the cells to uptake and excrete metabolites and use the metabolites to generate energy and building blocks for cell growth. This is a particularly novel aspect of the model. Certain defined criteria for the concentrations of intracellular metabolites lead to cancer cell growth, proliferation or death. Overall, we model the evolution of the tumor in both time and space. Starting with a cluster of tumor cells, the model produces an avascular tumor that quantitatively and qualitatively mimics experimental measurements of multicellular tumor spheroids. Through our model simulations, we can investigate the response of individual intracellular species under a metabolic perturbation and investigate how that response contributes to the response of the tumor as a whole. The predicted response of intracellular metabolites under various targeted strategies are difficult to resolve with experimental techniques. Thus, the model can give novel predictions as to the response of the tumor as a whole, identifies potential therapies to impede tumor growth, and predicts the effects of those therapeutic strategies. In particular, the model provides quantitative insight into the dynamic reprogramming of tumor cells at the intracellular level in response to specific metabolic perturbations. Overall, the model is a useful framework to study targeted metabolic strategies for inhibiting tumor growth.

Dataset Information

A modular framework for multiscale multicellular spatial modeling of viral infection, immune response and drug therapy timing and efficacy in epithelial tissues: A multiscale model of viral infection in epithelial tissues.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets