RESUMO
Interactions between cancer cells and immune cells in the tumor microenvironment influence tumor growth and can contribute to the response to cancer immunotherapies. It is difficult to gain mechanistic insights into the effects of cell-cell interactions in tumors using a purely experimental approach. However, computational modeling enables quantitative investigation of the tumor microenvironment, and agent-based modeling, in particular, provides relevant biological insights into the spatial and temporal evolution of tumors. Here, we develop a novel agent-based model (ABM) to predict the consequences of intercellular interactions. Furthermore, we leverage our prior work that predicts the transitions of CD8+ T cells from a naïve state to a terminally differentiated state using Boolean modeling. Given the details incorporated to predict T cell state, we apply the integrated Boolean-ABM framework to study how the properties of CD8+ T cells influence the composition and spatial organization of tumors and the efficacy of an immune checkpoint blockade. Overall, we present a mechanistic understanding of tumor evolution that can be leveraged to study targeted immunotherapeutic strategies.
RESUMO
Mathematical models of biomolecular networks are commonly used to study cellular processes; however, their usefulness to explain and predict dynamic behaviors is often questioned due to the unclear relationship between parameter uncertainty and network dynamics. In this work, we introduce PyDyNo (Python dynamic analysis of biochemical networks), a non-equilibrium reaction-flux based analysis to identify dominant reaction paths within a biochemical reaction network calibrated to experimental data. We first show, in a simplified apoptosis execution model, that despite the thousands of parameter vectors with equally good fits to experimental data, our framework identifies the dynamic differences between these parameter sets and outputs three dominant execution modes, which exhibit varying sensitivity to perturbations. We then apply our methodology to JAK2/STAT5 network in colony-forming unit-erythroid (CFU-E) cells and provide previously unrecognized mechanistic explanation for the survival responses of CFU-E cell population that would have been impossible to deduce with traditional protein-concentration based analyses.
RESUMO
T cells play a key role in a variety of immune responses, including infection and cancer. Upon stimulation, naïve CD8+ T cells proliferate and differentiate into a variety of memory and effector cell types; however, failure to clear antigens causes prolonged stimulation of CD8+ T cells, ultimately leading to T cell exhaustion (TCE). The functional and phenotypic changes that occur during CD8+ T cell differentiation are well characterized, but the underlying gene expression state changes are not completely understood. Here, we utilize a previously published data-driven Boolean model of gene regulatory interactions shown to mediate TCE. Our network analysis and modeling reveal the final gene expression states that correspond to TCE, along with the sequence of gene expression patterns that give rise to those final states. With a model that predicts the changes in gene expression that lead to TCE, we could evaluate strategies to inhibit the exhausted state. Overall, we demonstrate that a common pathway model of CD8+ T cell gene regulatory interactions can provide insights into the transcriptional changes underlying the evolution of cell states in TCE.
Assuntos
Neoplasias , Exaustão das Células T , Humanos , Linfócitos T CD8-Positivos/metabolismo , Diferenciação Celular/genética , Ativação Linfocitária , Neoplasias/genética , Neoplasias/metabolismoRESUMO
Necroptosis is a form of regulated cell death associated with degenerative disorders, autoimmune and inflammatory diseases, and cancer. To better understand the biochemical mechanisms regulating necroptosis, we constructed a detailed computational model of tumor necrosis factor-induced necroptosis based on known molecular interactions from the literature. Intracellular protein levels, used as model inputs, were quantified using label-free mass spectrometry, and the model was calibrated using Bayesian parameter inference to experimental protein time course data from a well-established necroptosis-executing cell line. The calibrated model reproduced the dynamics of phosphorylated mixed lineage kinase domain-like protein, an established necroptosis reporter. A subsequent dynamical systems analysis identified four distinct modes of necroptosis signal execution, distinguished by rate constant values and the roles of the RIP1 deubiquitinating enzymes A20 and CYLD. In one case, A20 and CYLD both contribute to RIP1 deubiquitination, in another RIP1 deubiquitination is driven exclusively by CYLD, and in two modes either A20 or CYLD acts as the driver with the other enzyme, counterintuitively, inhibiting necroptosis. We also performed sensitivity analyses of initial protein concentrations and rate constants to identify potential targets for modulating necroptosis sensitivity within each mode. We conclude by associating numerous contrasting and, in some cases, counterintuitive experimental results reported in the literature with one or more of the model-predicted modes of necroptosis execution. In all, we demonstrate that a consensus pathway model of tumor necrosis factor-induced necroptosis can provide insights into unresolved controversies regarding the molecular mechanisms driving necroptosis execution in numerous cell types under different experimental conditions.
Assuntos
Sinais (Psicologia) , Necroptose , Humanos , Necrose/metabolismo , Necrose/patologia , Teorema de Bayes , Fator de Necrose Tumoral alfa/farmacologia , ApoptoseRESUMO
Small cell lung cancer (SCLC) tumors comprise heterogeneous mixtures of cell states, categorized into neuroendocrine (NE) and non-neuroendocrine (non-NE) transcriptional subtypes. NE to non-NE state transitions, fueled by plasticity, likely underlie adaptability to treatment and dismal survival rates. Here, we apply an archetypal analysis to model plasticity by recasting SCLC phenotypic heterogeneity through multi-task evolutionary theory. Cell line and tumor transcriptomics data fit well in a five-dimensional convex polytope whose vertices optimize tasks reminiscent of pulmonary NE cells, the SCLC normal counterparts. These tasks, supported by knowledge and experimental data, include proliferation, slithering, metabolism, secretion, and injury repair, reflecting cancer hallmarks. SCLC subtypes, either at the population or single-cell level, can be positioned in archetypal space by bulk or single-cell transcriptomics, respectively, and characterized as task specialists or multi-task generalists by the distance from archetype vertex signatures. In the archetype space, modeling single-cell plasticity as a Markovian process along an underlying state manifold indicates that task trade-offs, in response to microenvironmental perturbations or treatment, may drive cell plasticity. Stifling phenotypic transitions and plasticity may provide new targets for much-needed translational advances in SCLC. A record of this paper's Transparent Peer Review process is included in the supplemental information.