Search

Supplemental Video 2

Explanation for Supplementary Videos

Immunotherapies are an emerging strategy for treatment of solid tumors. Improved understanding of the mechanisms employed by cytotoxic T lymphocytes (CTL) to control tumors will aid in the development of immunotherapies. CTLs can directly kill tumor cells in a contact-dependent manner or may exert indirect effects on tumor cells via secretion of cytokines. Here, we aim to quantify the importance of these mechanisms in murine thymoma EL4/EG7 cells. We developed an agent-based model (ABM) and an ordinary differential equation model of tumor regression after adoptive transfer of a population of CTLs. Models were parameterized based on in vivo measurements of CTL infiltration and killing rates applied to EL4/EG7 tumors and OTI T cells. We quantified whether infiltrating CTLs are capable of controlling tumors through only direct, contact-dependent killing. Both models agreed that the low measured killing rate of CTLs in vivo was insufficient to cause tumor regression. In our ABM, we also simulated CTL production of the cytokine IFNγ in order to explore how an antiproliferative effect of IFNγ might aid CTLs in tumor control. In this model, IFNγ substantially reduced tumor growth compared with direct killing alone. Collectively, these data demonstrate that contact-dependent killing is insufficient for EL4 regression in vivo and highlight the potential importance of cytokine-induced antiproliferative effects in T-cell–mediated tumor control.Significance:Computational modeling highlights the importance of cytokine-induced antiproliferative effects in T-cell–mediated control of tumor progression.

Supplementary Methods

While human regulatory risk assessment (RA) still largely relies on animal studies, new approach methodologies (NAMs) based on in vitro, in silico or non-mammalian alternative models are increasingly used to evaluate chemical hazards. Moreover, human epidemiological studies with biomarkers of effect (BoE) also play an invaluable role in identifying health effects associated with chemical exposures. To move towards the next generation risk assessment (NGRA), it is therefore crucial to establish bridges between NAMs and standard approaches, and to establish processes for increasing mechanistically-based biological plausibility in human studies. The Adverse Outcome Pathway (AOP) framework constitutes an important tool to address these needs but, despite a significant increase in knowledge and awareness, the use of AOPs in chemical RA remains limited. The objective of this paper is to address issues related to using AOPs in a regulatory context from various perspectives as it was discussed in a workshop organized within the European Union partnerships HBM4EU and PARC in spring 2022. The paper presents examples where the AOP framework has been proven useful for the human RA process, particularly in hazard prioritization and characterization, in integrated approaches to testing and assessment (IATA), and in the identification and validation of BoE in epidemiological studies. Nevertheless, several limitations were identified that hinder the optimal usability and acceptance of AOPs by the regulatory community including the lack of quantitative information on response-response relationships and of efficient ways to map chemical data (exposure and toxicity) onto AOPs. The paper summarizes suggestions, ongoing initiatives and third-party tools that may help to overcome these obstacles and thus assure better implementation of AOPs in the NGRA., European Commission 733032 857560 101057014, Ministry of Education, Youth and Sports by the RECETOX Research Infrastructure LM2018121, OP RDE project CETOCOEN Excellence CZ.02.1.01/0.0/0.0/17_043/0009632, Japan Agency for Medical Research and Development (AMED) JP21mk0101216 JP22mk0101216, Ministry of Education, Culture, Sports, Science and Technology, Japan (MEXT), Japan Society for the Promotion of Science, Grants-in-Aid for Scientific Research (KAKENHI) 21K12133

Epithelial–mesenchymal transition (EMT) and immunoevasion through upregulation of programmed death-ligand 1 (PD-L1) are important drivers of cancer progression. While EMT has been proposed to facilitate PD-L1-mediated immunosuppression, molecular mechanisms of their interaction remain obscure. Here, we provide insight into these mechanisms by proposing a mathematical model that describes the crosstalk between EMT and interferon gamma (IFN γ )-induced PD-L1 expression. Our model shows that via interaction with microRNA-200 (miR-200), the multi-stability of the EMT regulatory circuit is mirrored in PD-L1 levels, which are further amplified by IFN γ stimulation. This IFN γ -mediated effect is most prominent for cells in a fully mesenchymal state and less strong for those in an epithelial or partially mesenchymal state. In addition, bidirectional crosstalk between miR-200 and PD-L1 implies that IFN γ stimulation allows cells to undergo EMT for lower amounts of inducing signal, and the presence of IFN γ accelerates EMT and decelerates mesenchymal–epithelial transition (MET). Overall, our model agrees with published findings and provides insight into possible mechanisms behind EMT-mediated immune evasion, and primary, adaptive, or acquired resistance to immunotherapy. Our model can be used as a starting point to explore additional crosstalk mechanisms, as an improved understanding of these mechanisms is indispensable for developing better diagnostic and therapeutic options for cancer patients.

Technological development has led to ever-increasing amounts of data in high-content screening. For utilizing such data in an efficient, thorough, and user-friendly manner we developed the R package H5CellProfiler. H5CellProfiler is based on R packages data.table, ggvis, ggplot2 and shiny. H5CellProfiler launches a browser that allows scientists to analyze large single-cell datasets and make statistical summaries and graphs on their local desktop in a fast and memory-efficient manner. In addition, single-cell track labels are calculated and broken tracks are re-connected based on user-defined thresholds resulting in unique sets of annotated tracks.

Toxicity is not only a function of damage mechanisms, but is also determined by cellular resilience factors. Glutathione has been reported as essential element to counteract negative influences. The present work hence pursued the question how intracellular glutathione can be elevated transiently to render cells more resistant toward harmful conditions. The antibiotic nitrofurantoin (NFT) was identified to stimulate de novo synthesis of glutathione in the human hepatoma cell line, HepG2, and in primary human hepatocytes. In intact cells, activation of NFT yielded a radical anion, which subsequently initiated nuclear-factor-erythroid 2-related-factor-2 (Nrf2)-dependent induction of glutamate cysteine ligase (GCL). Application of siRNA-based intervention approaches confirmed the involvement of the Nrf2-GCL axis in the observed elevation of intracellular glutathione levels. Quantitative activation of Nrf2 by NFT, and the subsequent rise in glutathione, were similar as observed with the potent experimental Nrf2 activator diethyl maleate. The elevation of glutathione levels, observed even 48 h after withdrawal of NFT, rendered cells resistant to different stressors such as the mitochondrial inhibitor rotenone, the redox cycler paraquat, the proteasome inhibitors MG-132 or bortezomib, or high concentrations of NFT. Repurpose of the antibiotic NFT as activator of Nrf2 could thus be a promising strategy for a transient and targeted activation of the endogenous antioxidant machinery.

Epithelial-Mesenchymal Transition (EMT) and immunoevasion through upregulation of Programmed Death-Ligand 1 (PD-L1) are important drivers of cancer progression. While EMT has been proposed to facilitate PD-L1-mediated immunosuppression, the molecular mechanisms of their interaction remain obscure. Here we provide insight into these mechanisms by proposing a mathematical model that describes the crosstalk between EMT and Interferon gamma (IFNγ)-induced PD-L1 expression. Our model shows that via interaction with microRNA-200 (miR-200), the multistability of the EMT regulatory circuit is mirrored in the PD-L1 levels, which are further amplified by IFNγ stimulation. This IFNγ-mediated effect is most prominent for cells in a fully mesenchymal state, and less strong for those in an epithelial or partially mesenchymal state. Additionally, bi-directional crosstalk between miR-200 and PD-L1 implies that IFNγ stimulation allows cells to undergo EMT for lower amounts of inducing signal, and that IFNγ presence accelerates EMT and decelerates Mesenchymal-Epithelial Transition (MET). Overall, our model agrees with published findings and provides insight into possible mechanisms behind EMT-mediated immune-evasion; and primary, adaptive, or acquired resistance to immunotherapy. Our model can be used as a starting point to explore additional crosstalk mechanisms, as an improved understanding of these mechanisms is indispensable for developing better diagnostic and therapeutic options for cancer patients.

Mitochondrial perturbation is a key event in chemical-induced organ toxicities that is incompletely understood. Here, we studied how electron transport chain (ETC) complex I, II, or III (CI, CII and CIII) inhibitors affect mitochondrial functionality, stress response activation, and cell viability using a combination of high-content imaging and TempO-Seq in HepG2 hepatocyte cells. CI and CIII inhibitors perturbed mitochondrial membrane potential (MMP) and mitochondrial and cellular ATP levels in a concentration- and time-dependent fashion and, under conditions preventing a switch to glycolysis attenuated cell viability, whereas CII inhibitors had no effect. TempO-Seq analysis of changes in mRNA expression pointed to a shared cellular response to CI and CIII inhibition. First, to define specific ETC inhibition responses, a gene set responsive toward ETC inhibition (and not to genotoxic, oxidative, or endoplasmic reticulum stress) was identified using targeted TempO-Seq in HepG2. Silencing of one of these genes, NOS3, exacerbated the impact of CI and CIII inhibitors on cell viability, indicating its functional implication in cellular responses to mitochondrial stress. Then by monitoring dynamic responses to ETC inhibition using a HepG2 GFP reporter panel for different classes of stress response pathways and applying pathway and gene network analysis to TempO-Seq data, we looked for downstream cellular events of ETC inhibition and identified the amino acid response (AAR) as being triggered in HepG2 by ETC inhibition. Through in silico approaches we provide evidence indicating that a similar AAR is associated with exposure to mitochondrial toxicants in primary human hepatocytes. Altogether, we (i) unravel quantitative, time- and concentration-resolved cellular responses to mitochondrial perturbation, (ii) identify a gene set associated with adaptation to exposure to active ETC inhibitors, and (iii) show that ER stress and an AAR accompany ETC inhibition in HepG2 and primary hepatocytes. Supplementary Information The online version contains supplementary material available at 10.1007/s00204-021-03160-7.

The ability of cancer cells to invade neighboring tissue from primary tumors is an important determinant of metastatic behavior. Quantification of cell migration characteristics such as migration speed and persistence helps to understand the requirements for such invasiveness. One factor that may influence invasion is how local tumor cell density shapes cell migration characteristics, which we here investigate with a combined experimental and computational modeling approach. First, we generated and analyzed time-lapse imaging data on two aggressive Triple-Negative Breast Cancer (TNBC) cell lines, HCC38 and Hs578T, during 2D migration assays at various cell densities. HCC38 cells exhibited a counter-intuitive increase in speed and persistence with increasing density, whereas Hs578T did not exhibit such an increase. Moreover, HCC38 cells exhibited strong cluster formation with active pseudopod-driven migration, especially at low densities, whereas Hs578T cells maintained a dispersed positioning. In order to obtain a mechanistic understanding of the density-dependent cell migration characteristics and cluster formation, we developed realistic spatial simulations using a Cellular Potts Model (CPM) with an explicit description of pseudopod dynamics. Model analysis demonstrated that strong coordination between pseudopods within single cells could explain the experimentally observed increase in speed and persistence with increasing density in HCC38 cells. Thus, the density-dependent migratory behavior could be an emergent property of single-cell characteristics without the need for additional mechanisms. This implies that coordination amongst pseudopods may play a role in the aggressive nature of cancers through mediating dispersal.

Cells can adjust their mitochondrial morphology by altering the balance between mitochondrial fission and fusion to adapt to stressful conditions. The connection between a chemical perturbation, changes in mitochondrial function, and altered mitochondrial morphology is not well understood. Here, we made use of high-throughput high-content confocal microscopy to assess the effects of distinct classes of oxidative phosphorylation (OXPHOS) complex inhibitors on mitochondrial parameters in a concentration and time resolved manner. Mitochondrial morphology phenotypes were clustered based on machine learning algorithms and mitochondrial integrity patterns were mapped. In parallel, changes in mitochondrial membrane potential (MMP), mitochondrial and cellular ATP levels, and viability were microscopically assessed. We found that inhibition of MMP, mitochondrial ATP production, and oxygen consumption rate (OCR) using sublethal concentrations of complex I and III inhibitors did not trigger mitochondrial fragmentation. Instead, complex V inhibitors that suppressed ATP and OCR but increased MMP provoked a more fragmented mitochondrial morphology. In agreement, complex V but not complex I or III inhibitors triggered proteolytic cleavage of the mitochondrial fusion protein, OPA1. The relation between increased MMP and fragmentation did not extend beyond OXPHOS complex inhibitors: increasing MMP by blocking the mPTP pore did not lead to OPA1 cleavage or mitochondrial fragmentation and the OXPHOS uncoupler FCCP was associated with OPA1 cleavage and MMP reduction. Altogether, our findings connect vital mitochondrial functions and phenotypes in a high-throughput high-content confocal microscopy approach that help understanding of chemical-induced toxicity caused by OXPHOS complex perturbing chemicals.

Mitochondria are the main bioenergetic organelles of cells. Exposure to chemicals targeting mitochondria therefore generally results in the development of toxicity. The cellular response to perturbations in cellular energy production is a balance between adaptation, by reorganisation and organelle biogenesis, and sacrifice, in the form of cell death. In homeostatic conditions, aerobic mitochondrial energy production requires the maintenance of a mitochondrial membrane potential (MMP). Chemicals can perturb this MMP, and the extent of this perturbation depends both on the pharmacokinetics of the chemicals and on downstream MMP dynamics. Here we obtain a quantitative understanding of mitochondrial adaptation upon exposure to various mitochondrial respiration inhibitors by applying mathematical modeling to partially published high-content imaging time-lapse confocal imaging data, focusing on MMP dynamics in HepG2 cells over a period of 24 h. The MMP was perturbed using a set of 24 compounds, either acting as uncoupler or as mitochondrial complex inhibitor targeting complex I, II, III or V. To characterize the effect of chemical exposure on MMP dynamics, we adapted an existing differential equation model and fitted this model to the observed MMP dynamics. Complex III inhibitor data were better described by the model than complex I data. Incorporation of pharmacokinetic decay into the model was required to obtain a proper fit for the uncoupler FCCP. Furthermore, oligomycin (complex V inhibitor) model fits were improved by either combining pharmacokinetic (PK) decay and ion leakage or a concentration-dependent decay. Subsequent mass spectrometry measurements showed that FCCP had a significant decay in its PK profile as predicted by the model. Moreover, the measured oligomycin PK profile exhibited only a limited decay at high concentration, whereas at low concentrations the compound remained below the detection limit within cells. This is consistent with the hypothesis that oligomycin exhibits a concentration-dependent decay, yet awaits further experimental verification with more sensitive detection methods. Overall, we show that there is a complex interplay between PK and MMP dynamics within mitochondria and that data-driven modeling is a powerful combination to unravel such complexity.

Simple Summary Cytotoxic T lymphocytes (CTLs) play an important role in controlling tumours, and an improved understanding of how they accomplish this will benefit immunotherapeutic cancer treatment strategies. Stimulation of CTLs by targeting their CD137 receptor is a strategy currently under investigation for enhancing responses against tumours, yet so far only limited quantitative knowledge regarding the effects of such stimulation upon CTLs has been obtained. Here, we develop mathematical models to describe dynamic in vivo two-photon imaging of tumour infiltrating CTLs, to characterise differences in their function either in the presence or absence of a CD137 agonist antibody. We showed that an increased antiproliferative effect and a more sustained presence of CTLs within the tumour were the most important effects associated with anti-CD137 treatment. Abstract Several immunotherapeutic strategies for the treatment of cancer are under development. Two prominent strategies are adoptive cell transfer (ACT) of CTLs and modulation of CTL function with immune checkpoint inhibitors or with costimulatory antibodies. Despite some success with these approaches, there remains a lack of detailed and quantitative descriptions of the events following CTL transfer and the impact of immunomodulation. Here, we have applied ordinary differential equation models to two photon imaging data derived from a B16F10 murine melanoma. Models were parameterised with data from two different treatment conditions: either ACT-only, or ACT with intratumoural costimulation using a CD137 targeted antibody. Model dynamics and best fitting parameters were compared, in order to assess the mode of action of the CTLs and examine how the CD137 antibody influenced their activities. We found that the cytolytic activity of the transferred CTLs was minimal without CD137 costimulation, and that the CD137 targeted antibody did not enhance the per-capita killing ability of the transferred CTLs. Instead, the results of our modelling study suggest that an antiproliferative effect of CTLs exerted upon the tumour likely accounted for the majority of the reduction in tumour growth after CTL transfer. Moreover, we found that CD137 most likely improved tumour control via enhancement of this antiproliferative effect, as well as prolonging the period in which CTLs were inside the tumour, leading to a sustained duration of their antitumour effects following CD137 stimulation.

Drug-induced liver injury (DILI) is the most prevalent adversity encountered in drug development and clinical settings leading to urgent needs to understand the underlying mechanisms. In this study, we have systematically investigated the dynamics of the activation of cellular stress response pathways and cell death outcomes upon exposure of a panel of liver toxicants using live cell imaging of fluorescent reporter cell lines. We established a comprehensive temporal dynamic response profile of a large set of BAC-GFP HepG2 cell lines representing the following components of stress signaling: i) unfolded protein response (UPR) [ATF4, XBP1, BIP and CHOP]; ii) oxidative stress [NRF2, SRXN1, HMOX1]; iii) DNA damage [P53, P21, BTG2, MDM2]; and iv) NF-κB pathway [A20, ICAM1]. We quantified the single cell GFP expression as a surrogate for endogenous protein expression using live cell imaging over > 60 h upon exposure to 14 DILI compounds at multiple concentrations. Using logic-based ordinary differential equation (Logic-ODE), we modelled the dynamic profiles of the different stress responses and extracted specific descriptors potentially predicting the progressive outcomes. We identified the activation of ATF4-CHOP axis of the UPR as the key pathway showing the highest correlation with cell death upon DILI compound perturbation. Knocking down main components of the UPR provided partial protection from compound-induced cytotoxicity, indicating a complex interplay among UPR components as well as other stress pathways. Our results suggest that a systematic analysis of the temporal dynamics of ATF4-CHOP axis activation can support the identification of DILI risk for new candidate drugs. published

Several immunotherapeutic strategies for the treatment of cancer are under development. Two prominent strategies are adoptive cell transfer (ACT) of cytotoxic T lymphocytes (CTLs) and modulation of CTL function with immune checkpoint inhibitors or with costimulatory antibodies. Despite some success with these approaches, there remains a lack of detailed and quantitative descriptions of the events following CTL transfer and the impact of immunomodulation. Here, we have applied ordinary differential equation models to two photon imaging data derived from a B16F10 murine melanoma. Models were parameterised with data from two different treatment conditions: either ACT-only, or ACT with intratumoural costimulation using a CD137 targeted antibody. Model dynamics and best fitting parameters were compared, in order to assess the mode of action of the CTLs and examine how the CD137 antibody influenced their activities. We found that the cytolytic activity of the transferred CTLs was minimal without CD137 costimulation, and that the CD137 targeted antibody did not enhance the per-capita killing ability of the transferred CTLs. Instead, the results of our modelling study suggest that an antiproliferative effect of CTLs exerted upon the tumour likely accounted for the majority of the reduction in tumour growth after CTL transfer. We found that CD137 most likely improved tumour control via enhancement of this antiproliferative effect, as well as prolonging the period in which CTLs were inside the tumour, leading to a sustained duration of their antitumour effects following CD137 stimulation.

Following local inflammation, CD8+ T cells develop into tissue-resident memory T cells (TRM). Through combined lineage tracing and single-cell transcriptome analysis, the authors show that the circulating T cell pool harbors TRM precursors that commit to the TRM fate before tissue entry., An increasing body of evidence emphasizes the role of tissue-resident memory T cells (TRM) in the defense against recurring pathogens and malignant neoplasms. However, little is known with regard to the origin of these cells and their kinship to other CD8+ T cell compartments. To address this issue, we followed the antigen-specific progeny of individual naive CD8+ T cells to the T effector (TEFF), T circulating memory (TCIRCM), and TRM pools by lineage-tracing and single-cell transcriptome analysis. We demonstrate that a subset of T cell clones possesses a heightened capacity to form TRM, and that enriched expression of TRM–fate-associated genes is already apparent in the circulating TEFF offspring of such clones. In addition, we demonstrate that the capacity to generate TRM is permanently imprinted at the clonal level, before skin entry. Collectively, these data provide compelling evidence for early stage TRM fate decisions and the existence of committed TRM precursor cells in the circulatory TEFF compartment., Graphical Abstract

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Adverse outcome pathways (AOPs) are a recent toxicological construct that connects, in a formalized, transparent and quality-controlled way, mechanistic information to apical endpoints for regulatory purposes. AOP links a molecular initiating event (MIE) to the adverse outcome (AO) via key events (KE), in a way specified by key event relationships (KER). Although this approach to formalize mechanistic toxicological information only started in 2010, over 200 AOPs have already been established. At this stage, new requirements arise, such as the need for harmonization and re-assessment, for continuous updating, as well as for alerting about pitfalls, misuses and limits of applicability. In this review, the history of the AOP concept and its most prominent strengths are discussed, including the advantages of a formalized approach, the systematic collection of weight of evidence, the linkage of mechanisms to apical end points, the examination of the plausibility of epidemiological data, the identification of critical knowledge gaps and the design of mechanistic test methods. To prepare the ground for a broadened and appropriate use of AOPs, some widespread misconceptions are explained. Moreover, potential weaknesses and shortcomings of the current AOP rule set are addressed (1) to facilitate the discussion on its further evolution and (2) to better define appropriate vs. less suitable application areas. Exemplary toxicological studies are presented to discuss the linearity assumptions of AOP, the management of event modifiers and compensatory mechanisms, and whether a separation of toxicodynamics from toxicokinetics including metabolism is possible in the framework of pathway plasticity. Suggestions on how to compromise between different needs of AOP stakeholders have been added. A clear definition of open questions and limitations is provided to encourage further progress in the field. © 2017, Springer-Verlag GmbH Germany.

Cazaux et al. use intravital imaging to dissect anti-CD19 CAR T cell activity. This study uncovers both anatomical and functional diversity in the outcome of anti-CD19 CAR T cell interactions with tumor cells impacting engraftment, killing dynamics, and tumor immunoediting., CAR T cells represent a potentially curative strategy for B cell malignancies. However, the outcome and dynamics of CAR T cell interactions in distinct anatomical sites are poorly understood. Using intravital imaging, we tracked interactions established by anti-CD19 CAR T cells in B cell lymphoma–bearing mice. Circulating targets trapped CAR T cells in the lungs, reducing their access to lymphoid organs. In the bone marrow, tumor apoptosis was largely due to CAR T cells that engaged, killed, and detached from their targets within 25 min. Notably, not all CAR T cell contacts elicited calcium signaling or killing while interacting with tumors, uncovering extensive functional heterogeneity. Mathematical modeling revealed that direct killing was sufficient for tumor regression. Finally, antigen-loss variants emerged in the bone marrow, but not in lymph nodes, where CAR T cell cytotoxic activity was reduced. Our results identify a previously unappreciated level of diversity in the outcomes of CAR T cell interactions in vivo, with important clinical implications., Graphical Abstract

Immunotherapies are an emerging strategy for treatment of solid tumors. Improved understanding of the mechanisms employed by cytotoxic T lymphocytes (CTL) to control tumors will aid in the development of immunotherapies. CTLs can directly kill tumor cells in a contact-dependent manner or may exert indirect effects on tumor cells via secretion of cytokines. Here, we aim to quantify the importance of these mechanisms in murine thymoma EL4/EG7 cells. We developed an agent-based model (ABM) and an ordinary differential equation model of tumor regression after adoptive transfer of a population of CTLs. Models were parameterized based on in vivo measurements of CTL infiltration and killing rates applied to EL4/EG7 tumors and OTI T cells. We quantified whether infiltrating CTLs are capable of controlling tumors through only direct, contact-dependent killing. Both models agreed that the low measured killing rate of CTLs in vivo was insufficient to cause tumor regression. In our ABM, we also simulated CTL production of the cytokine IFNγ in order to explore how an antiproliferative effect of IFNγ might aid CTLs in tumor control. In this model, IFNγ substantially reduced tumor growth compared with direct killing alone. Collectively, these data demonstrate that contact-dependent killing is insufficient for EL4 regression in vivo and highlight the potential importance of cytokine-induced antiproliferative effects in T-cell–mediated tumor control. Significance: Computational modeling highlights the importance of cytokine-induced antiproliferative effects in T-cell–mediated control of tumor progression.

While hazard assessment of chemicals can make direct use of descriptive adverse outcome pathways (AOPs), risk assessment requires quantitative relationships from exposure to effect timing and magnitude. To seamlessly integrate the data generated by alternative methods or in vivo testing, quantitative AOPs (qAOPs) providing dose-time-response predictions are more valuable than qualitative AOPs. Here, we compare three approaches to qAOP building: empirical dose-response modeling, Bayesian network (BN) calibration, and systems biology (SB) modeling. These methods were applied to the quantification of a simplified oxidative stress induced chronic kidney disease AOP, on the basis of in vitro data obtained on RPTEC/TERT1 cells exposed to potassium bromate. Effectopedia was used to store the experimental data and the developed models in a unified representation so they can be compared and further analyzed. We argue that despite the fact that dose-response models give adequate fits to the data they should be accompanied by mechanistic SB modeling to gain a proper perspective on the quantification. BNs can be both more precise than dose-response models and simpler than SB models, but more experience with their usage is needed.

T cell responses upon infection display a remarkably reproducible pattern of expansion, contraction, and memory formation. If the robustness of this pattern builds entirely on signals derived from other cell types or if activated T cells themselves contribute to the orchestration of these population dynamics-akin to bacterial quorum regulation-is unclear. Here, we examined this question using time-lapse microscopy, genetic perturbation, bioinformatic predictions, and mathematical modeling. We found that ICAM-1-mediated cell clustering enabled CD8+ T cells to collectively regulate the balance between proliferation and apoptosis. Mechanistically, T cell expressed CD80 and CD86 interacted with the receptors CD28 and CTLA-4 on neighboring T cells; these interactions fed two nested antagonistic feedback circuits that regulated interleukin 2 production in a manner dependent on T cell density as confirmed by in vivo modulation of this network. Thus, CD8+ T cell-population-intrinsic mechanisms regulate cellular behavior, thereby promoting robustness of population dynamics.

Emerging data show that tissue-resident memory T (T

Oxidative stress leads to the activation of the Nuclear factor-erythroid-2-related factor 2 (Nrf2) pathway. While most studies have focused on the activation of the Nrf2 pathway after single chemical treatment, little is known about the dynamic regulation of the Nrf2 pathway in the context of repeated exposure scenarios. Here we employed single cell live imaging to quantitatively monitor the dynamics of the Nrf2 pathway during repeated exposure, making advantage of two HepG2 fluorescent protein reporter cell lines, expressing GFP tagged Nrf2 or sulfiredoxin 1 (Srxn1), a direct downstream target of Nrf2. High throughput live confocal imaging was used to measure the temporal dynamics of these two components of the Nrf2 pathway after repeated exposure to an extensive concentration range of diethyl maleate (DEM) and tert-butylhydroquinone (tBHQ). Single treatment with DEM or tBHQ induced Nrf2 and Srxn1 over time in a concentration-dependent manner. The Nrf2 response to a second treatment was lower than the response to the first exposure with the same concentration, indicating that the response is adaptive. Moreover, a limited fraction of individual cells committed themselves into the Nrf2 response during the second treatment. Despite the suppression of the Nrf2 pathway, the second treatment resulted in a three-fold higher Srxn1-GFP response compared to the first treatment, with all cells participating in the response. While after the first treatment Srxn1-GFP response was linearly related to Nrf2-GFP nuclear translocation, such a linear relationship was less clear for the second exposure. siRNA-mediated knockdown demonstrated that the second response is dependent on the activity of Nrf2. Several other, clinically relevant, compounds (i.e., sulphorophane, nitrofurantoin and CDDO-Me) also enhanced the induction of Srxn1-GFP upon two consecutive repeated exposure. Together the data indicate that adaptation towards pro-oxidants lowers the Nrf2 activation capacity, but simultaneously primes cells for the enhancement of an antioxidant response which depends on factors other than just Nrf2. These data provide further insight in the overall dynamics of stress pathway activation after repeated exposure and underscore the complexity of responses that may govern repeated dose toxicity.

Introduction: Drug-induced liver injury (DILI) is a significant threat to human health and a major problem in drug development. It is hard to predict due to its idiosyncratic nature and which does not show up in animal trials. Hepatic adaptive stress response pathway activation is generally observed in drug-induced liver injury. Dynamical pathway modeling has the potential to foresee adverse effects of drugs before they go in trial. Ordinary differential equation modeling can offer mechanistic insight, and allows us to study the dynamical behavior of stress pathways involved in DILI. Areas covered: This review provides an overview on the progress of the dynamical modeling of stress and death pathways pertinent to DILI, i.e. pathways relevant for oxidative stress, inflammatory stress, DNA damage, unfolded proteins, heat shock and apoptosis. We also discuss the required steps for applying such modeling to the liver. Expert opinion: Despite the strong progress made since the turn of the century, models of stress pathways have only rarely been specifically applied to describe pathway dynamics for DILI. We argue that with minor changes, in some cases only to parameter values, many of these models can be repurposed for application in DILI research. Combining both dynamical models with in vitro testing might offer novel screening methods for the harmful side-effects of drugs. KEYWORDS: Drug induced liver injury (DILI), dynamical modeling, ordinary differential equations, stress pathway, death pathway, predictive toxicology, cell fate, adverse outcome pathway

In healthy individuals, immune-checkpoint molecules prevent autoimmune responses and limit immune cell-mediated tissue damage. Tumors frequently exploit these molecules to evade eradication by the immune system. Over the past years, immune-checkpoint blockade of cytotoxic T lymphocyte antigen-4 and programed death-1 emerged as promising strategies to activate antitumor cytotoxic T cell responses. Although complete regression and long-term survival is achieved in some patients, not all patients respond. This review describes promising, novel combination approaches involving immune-checkpoint blockade in the context of the cancer-immunity cycle, aimed at increasing response rates to the single treatments. Specifically, we discuss combinations that promote antigen release and presentation, that further amplify T cell activation, that inhibit trafficking of regulatory T cells or MSDCs, that stimulate intratumoral T cell infiltration, that increase cancer recognition by T cells, and that stimulate tumor killing.

There is strong evidence that both adoptive T cell transfer and T cell checkpoint blockade can lead to regression of human melanoma. However, little data are available on the effect of these cancer therapies on the tumor-reactive T cell compartment. To address this issue we have profiled therapy-induced T cell reactivity against a panel of 145 melanoma-associated CD8(+) T cell epitopes. Using this approach, we demonstrate that individual tumor-infiltrating lymphocyte cell products from melanoma patients contain unique patterns of reactivity against shared melanoma-associated antigens, and that the combined magnitude of these responses is surprisingly low. Importantly, TIL therapy increases the breadth of the tumor-reactive T cell compartment in vivo, and T cell reactivity observed post-therapy can almost in full be explained by the reactivity observed within the matched cell product. These results establish the value of high-throughput monitoring for the analysis of immuno-active therapeutics and suggest that the clinical efficacy of TIL therapy can be enhanced by the preparation of more defined tumor-reactive T cell products.

One of the main questions in the field of imaging immune cell migration in living tissues is whether cells fulfill their functionality via random or nonrandom migration processes. For some applications, this issue has remained controversial even after publication of various imaging studies. A prime example is B-cell migration in germinal centers (GCs) where somatic hypermutation and clonal selection of B cells are thought to occur within morphologically distinct regions termed dark zone (DZ) and light zone (LZ). Here, we reanalyze a previously published dataset on GC B-cell migration, applying a sensitive analysis technique to detect directed migration and using a procedure to correct for a number of artifacts that frequently occur in time-lapse imaging experiments. Although B cells roughly perform a persistent random walk, we present evidence that they have a small preference (of on average about 0.2−0.3 μm min −1 ) to migrate from DZ to LZ, which is consistent with classical views of the GC reaction. This preference is most pronounced among a large subset of almost half of the B-cell population migrating along relatively straight tracks. Using a computational model to generate long-lasting B-cell tracks based on the experimental motility data (including the small directional preference), we predict a time course to travel from DZ to LZ of a few hours. This is consistent with experimental observations, and we show that at the observed cellular motility such a time course cannot be explained without the small preferential migration from DZ to LZ.

The visualization of the dynamic behaviour of and interactions between immune cells using time-lapse video microscopy has an important role in modern immunology. To draw robust conclusions, quantification of such cell migration is required. However, imaging experiments are associated with various artefacts that can affect the estimated positions of the immune cells under analysis, which form the basis of any subsequent analysis. Here, we describe potential artefacts that could affect the interpretation of data sets on immune cell migration. We propose how these errors can be recognized and corrected, and suggest ways to prevent the data analysis itself leading to biased results.