Explaining in-vitro to in-vivo efficacy correlations in oncology pre-clinical development via a semi-mechanistic mathematical model.

In-vitro to in-vivo correlations (IVIVC), relating in-vitro parameters like IC50 to in-vivo drug exposure in plasma and tumour growth, are widely used in oncology for experimental design and dose decisions. However, they lack a deeper understanding of the underlying mechanisms. Our paper therefore focuses on linking empirical IVIVC relations for small-molecule kinase inhibitors with a semi-mechanistic tumour-growth model. We develop an approach incorporating parameters like the compound's peak-trough ratio (PTR), Hill coefficient of in-vitro dose-response curves, and xenograft-specific properties. This leads to formulas for determining efficacious doses for tumor stasis under linear pharmacokinetics equivalent to traditional empirical IVIVC relations, but enabling more systematic analysis. Our findings reveal that in-vivo xenograft-specific parameters, specifically the growth rate (g) and decay rate (d), along with the average exposure, are generally more significant determinants of tumor stasis and effective dose than the compound's peak-trough ratio. However, as the Hill coefficient increases, the dependency of tumor stasis on the PTR becomes more pronounced, indicating that the compound is more influenced by its maximum or trough values rather than the average exposure. Furthermore, we discuss the translation of our method to predict population dose ranges in clinical studies and propose a resistance mechanism that solely relies on specific in-vivo xenograft parameters instead of IC50 exposure coverage. In summary, our study aims to provide a more mechanistic understanding of IVIVC relations, emphasizing the importance of xenograft-specific parameters and PTR on tumor stasis.

Radius additivity score: a novel combination index for tumour growth inhibition in fixed-dose xenograft studies.

The effect of combination therapies in many cancers has often been shown to be superior to that of monotherapies. This success is commonly attributed to drug synergies. Combinations of two (or more) drugs in xenograft tumor growth inhibition (TGI) studies are typically designed at fixed doses for each compound. The available methods for assessing synergy in such study designs are based on combination indices (CI) and model-based analyses. The former methods are suitable for screening exercises but are difficult to verify in in vivo studies, while the latter incorporate drug synergy in semi-mechanistic frameworks describing disease progression and drug action but are unsuitable for screening. In the current study, we proposed the empirical radius additivity (Rad-add) score, a novel CI for synergy detection in fixed-dose xenograft TGI combination studies. The Rad-add score approximates model-based analysis performed using the semi-mechanistic constant-radius growth TGI model. The Rad-add score was compared with response additivity, defined as the addition of the two response values, and the bliss independence model in combination studies derived from the Novartis PDX dataset. The results showed that the bliss independence and response additivity models predicted synergistic interactions with high and low probabilities, respectively. The Rad-add score predicted synergistic probabilities that appeared to be between those predicted with response additivity and the Bliss model. We believe that the Rad-add score is particularly suitable for assessing synergy in the context of xenograft combination TGI studies, as it combines the advantages of CI approaches suitable for screening exercises with those of semi-mechanistic TGI models based on a mechanistic understanding of tumor growth.

Complementary CRISPR screen highlights the contrasting role of membrane-bound and soluble ICAM-1 in regulating antigen-specific tumor cell killing by cytotoxic T cells.

Cytotoxic CD8 +T lymphocytes (CTLs) are key players of adaptive anti-tumor immunity based on their ability to specifically recognize and destroy tumor cells. Many cancer immunotherapies rely on unleashing CTL function. However, tumors can evade killing through strategies which are not yet fully elucidated. To provide deeper insight into tumor evasion mechanisms in an antigen-dependent manner, we established a human co-culture system composed of tumor and primary immune cells. Using this system, we systematically investigated intrinsic regulators of tumor resistance by conducting a complementary CRISPR screen approach. By harnessing CRISPR activation (CRISPRa) and CRISPR knockout (KO) technology in parallel, we investigated gene gain-of-function as well as loss-of-function across genes with annotated function in a colon carcinoma cell line. CRISPRa and CRISPR KO screens uncovered 187 and 704 hits, respectively, with 60 gene hits overlapping between both. These data confirmed the role of interferon-γ (IFN-γ), tumor necrosis factor α (TNF-α) and autophagy pathways and uncovered novel genes implicated in tumor resistance to killing. Notably, we discovered that ILKAP encoding the integrin-linked kinase-associated serine/threonine phosphatase 2 C, a gene previously unknown to play a role in antigen specific CTL-mediated killing, mediate tumor resistance independently from regulating antigen presentation, IFN-γ or TNF-α responsiveness. Moreover, our work describes the contrasting role of soluble and membrane-bound ICAM-1 in regulating tumor cell killing. The deficiency of membrane-bound ICAM-1 (mICAM-1) or the overexpression of soluble ICAM-1 (sICAM-1) induced resistance to CTL killing, whereas PD-L1 overexpression had no impact. These results highlight the essential role of ICAM-1 at the immunological synapse between tumor and CTL and the antagonist function of sICAM-1.

Transcriptomic and Proteomic Changes Driving Pulmonary Fibrosis Resolution in Young and Old Mice.

Bleomycin-induced pulmonary fibrosis in mice mimics major hallmarks of idiopathic pulmonary fibrosis. Yet in this model, it spontaneously resolves over time. We studied molecular mechanisms of fibrosis resolution and lung repair, focusing on transcriptional and proteomic signatures and the effect of aging. Old mice showed incomplete and delayed lung function recovery 8 weeks after bleomycin instillation. This shift in structural and functional repair in old bleomycin-treated mice was reflected in a temporal shift in gene and protein expression. We reveal gene signatures and signaling pathways that underpin the lung repair process. Importantly, the downregulation of WNT, BMP, and TGFß antagonists Frzb, Sfrp1, Dkk2, Grem1, Fst, Fstl1, and Inhba correlated with lung function improvement. Those genes constitute a network with functions in stem cell pathways, wound, and pulmonary healing. We suggest that insufficient and delayed downregulation of those antagonists during fibrosis resolution in old mice explains the impaired regenerative outcome. Together, we identified signaling pathway molecules with relevance to lung regeneration that should be tested in-depth experimentally as potential therapeutic targets for pulmonary fibrosis.

Antifibrotic Drug Nintedanib Inhibits CSF1R to Promote IL-4-associated Tissue Repair Macrophages.

Profibrotic and prohomeostatic macrophage phenotypes remain ill-defined, both in vivo and in vitro, impeding the successful development of drugs that reprogram macrophages as an attractive therapeutic approach to manage fibrotic disease. The goal of this study was to reveal profibrotic and prohomeostatic macrophage phenotypes that could guide the design of new therapeutic approaches targeting macrophages to treat fibrotic disease. This study used nintedanib, a broad kinase inhibitor approved for idiopathic pulmonary fibrosis, to dissect lung macrophage phenotypes during fibrosis-linked inflammation by combining in vivo and in vitro bulk and single-cell RNA-sequencing approaches. In the bleomycin model, nintedanib drove the expression of IL-4/IL-13-associated genes important for tissue regeneration and repair at early and late time points in lung macrophages. These findings were replicated in vitro in mouse primary bone marrow-derived macrophages exposed to IL-4/IL-13 and nintedanib. In addition, nintedanib promoted the expression of IL-4/IL-13 pathway genes in human macrophages in vitro. The molecular mechanism was connected to inhibition of the colony stimulating factor 1 (CSF1) receptor in both human and mouse macrophages. Moreover, nintedanib counterbalanced the effects of TNF on IL-4/IL-13 in macrophages to promote expression of IL-4/IL-13-regulated tissue repair genes in fibrotic contexts in vivo and in vitro. This study demonstrates that one of nintedanib's antifibrotic mechanisms is to increase IL-4 signaling in macrophages through inhibition of the CSF1 receptor, resulting in the promotion of tissue repair phenotypes.

Graft-Versus-Host Disease Prevention by In Vitro-Generated Myeloid-Derived Suppressor Cells Is Exclusively Mediated by the CD11b+CD11c+ MDSC Subpopulation.

Myeloid-derived suppressor cells (MDSCs) are a heterogeneous population of myeloid progenitor cells that dampen overwhelming adaptive immune responses through multiple mechanisms and are recognized as an attractive novel immune intervention therapy for counteracting the destructive effects of graft- versus -host disease (GVHD) developing after allogeneic bone marrow transplantation (BMT). MDSCs can be produced in great numbers for cellular therapy, but they present a mixture of subsets whose functions in GVHD prevention are undefined. Here, we generated MDSCs in vitro from murine BM cells in the presence of GM-CSF and defined the integrin CD11c as a marker to subdivide MDSCs into two functional subgroups: CD11b+CD11c+ and CD11b+CD11c- MDSCs. Isolated CD11b+CD11c+ and CD11b+CD11c- MDSCs both inhibited alloantigen-stimulated T-cell proliferation in vitro, although CD11b+CD11c+ MDSCs were more efficient and expressed higher levels of different immunosuppressive molecules. Likewise, expression of surface markers such as MHC class II, CD80, CD86, or PD-L1 further delineated both subsets. Most importantly, only the adoptive transfer of CD11b+CD11c+ MDSCs into a single MHC class I-disparate allogeneic BMT model prevented GVHD development and strongly decreased disease-induced mortality, while CD11b+CD11c- MDSCs were totally ineffective. Surprisingly, allogeneic T-cell homing and expansion in lymphatic and GVHD target organs were not affected by cotransplanted CD11b+CD11c+ MDSCs indicating a clear contradiction between in vitro and in vivo functions of MDSCs. However, CD11b+CD11c+ MDSCs shifted immune responses towards type 2 immunity reflected by increased Th2-specific cytokine expression of allogeneic T cells. Induction of type 2 immunity was mandatory for GVHD prevention, since CD11b+CD11c+ MDSCs were ineffective if recipients were reconstituted with STAT6-deficient T cells unable to differentiate into Th2 cells. Most importantly, the beneficial graft- versus -tumor (GVT) effect was maintained in the presence of CD11b+CD11c+ MDSCs since syngeneic tumor cells were efficiently eradicated. Strong differences in the transcriptomic landscape of both subpopulations underlined their functional differences. Defining CD11b+CD11c+ MDSCs as the subset of in vitro-generated MDSCs able to inhibit GVHD development might help to increase efficiency of MDSC therapy and to further delineate relevant target molecules and signaling pathways responsible for GVHD prevention.

CD276 is an important player in macrophage recruitment into the tumor and an upstream regulator for PAI-1.

More than 70% of colorectal, prostate, ovarian, pancreatic and breast cancer specimens show expression of CD276 (B7-H3), a potential immune checkpoint family member. Several studies have shown that high CD276 expression in cancer cells correlates with a poor clinical prognosis. This has been associated with the presence of lower tumor infiltrating leukocytes. Among those, tumor-associated macrophages can comprise up to 50% of the tumor mass and are thought to support tumor growth through various mechanisms. However, a lack of information on CD276 function and interaction partner(s) impedes rigorous evaluation of CD276 as a therapeutic target in oncology. Therefore, we aimed to understand the relevance of CD276 in tumor-macrophage interaction by employing a 3D spheroid coculture system with human cells. Our data show a role for tumor-expressed CD276 on the macrophage recruitment into the tumor spheroid, and also in regulation of the extracellular matrix modulator PAI-1. Furthermore, our experiments focusing on macrophage-expressed CD276 suggest that the antibody-dependent CD276 engagement triggers predominantly inhibitory signaling networks in human macrophages.

Response to IL-6 trans- and IL-6 classic signalling is determined by the ratio of the IL-6 receptor α to gp130 expression: fusing experimental insights and dynamic modelling.

BACKGROUND: Interleukin-6 is a pleiotropic cytokine with high clinical relevance and an important mediator of cellular communication, orchestrating both pro- and anti-inflammatory processes. Interleukin-6-induced signalling is initiated by binding of IL-6 to the IL-6 receptor α and subsequent binding to the signal transducing receptor subunit gp130. This active receptor complex initiates signalling through the Janus kinase/signal transducer and activator of transcription pathway. Of note, IL-6 receptor α exists in a soluble and a transmembrane form. Binding of IL-6 to membrane-bound IL-6 receptor α induces anti-inflammatory classic signalling, whereas binding of IL-6 to soluble IL-6 receptor α induces pro-inflammatory trans-signalling. Trans-signalling has been described to be markedly stronger than classic signalling. Understanding the molecular mechanisms that drive differences between trans- and classic signalling is important for the design of trans-signalling-specific therapies. These differences will be addressed here using a combination of dynamic mathematical modelling and molecular biology. METHODS: We apply an iterative systems biology approach using set-based modelling and validation approaches combined with quantitative biochemical and cell biological analyses. RESULTS: The combination of experimental analyses and dynamic modelling allows to relate the observed differences between IL-6-induced trans- and classic signalling to cell-type specific differences in the expression and ratios of the individual subunits of the IL-6 receptor complex. Canonical intracellular Jak/STAT signalling is indifferent in IL-6-induced trans- and classic signalling. CONCLUSION: This study contributes to the understanding of molecular mechanisms of IL-6 signal transduction and underlines the power of combined dynamical modelling, model-based validation and biological experiments. The opposing pro- and anti-inflammatory responses initiated by IL-6 trans- and classic signalling depend solely on the expression ratios of the subunits of the entire receptor complex. By pointing out the importance of the receptor expression ratio for the strength of IL-6 signalling this study lays a foundation for future precision medicine approaches that aim to selectively block pro-inflammatory trans-signalling. Furthermore, the derived models can be used for future therapy design.

The balance of interleukin (IL)-6, IL-6·soluble IL-6 receptor (sIL-6R), and IL-6·sIL-6R·sgp130 complexes allows simultaneous classic and trans-signaling.

Interleukin (IL-)6 is the major pro-inflammatory cytokine within the IL-6 family. IL-6 signals via glycoprotein 130 (gp130) and the membrane-bound or soluble IL-6 receptor (IL-6R), referred to as classic or trans-signaling, respectively. Whereas inflammation triggers IL-6 expression, eventually rising to nanogram/ml serum levels, soluble IL-6R (sIL-6R) and soluble gp130 (sgp130) are constitutively present in the upper nanogram/ml range. Calculations based on intermolecular affinities have suggested that systemic IL-6 is immediately trapped in IL-6·sIL-6R and IL-6·sIL-6R·sgp130 complexes, indicating that sIL-6R and sgp130 constitute a buffer system that increases the serum half-life of IL-6 or restricts systemic IL-6 signaling. However, this scenario has not been experimentally validated. Here, we quantified IL-6·sIL-6R and IL-6·sIL-6R·sgp130 complexes over a wide concentration range. The amounts of IL-6 used in this study reflect concentrations found during active inflammatory events. Our results indicated that most IL-6 is free and not complexed with sIL-6R or sgp130, indicating that the level of endogenous sgp130 in the bloodstream is not sufficient to block IL-6 trans-signaling via sIL-6R. Importantly, addition of the single-domain antibody VHH6, which specifically stabilizes IL-6·sIL-6R complexes but did not bind to IL-6 or sIL-6R alone, drove free IL-6 into IL-6·sIL-6R complexes and boosted trans-signaling but not classic signaling, demonstrating that endogenous sIL-6R has at least the potential to form complexes with IL-6. Our findings indicate that even though high concentrations of sIL-6R and sgp130 are present in human serum, the relative ratio of free IL-6 to IL-6·sIL-6R allows for simultaneous classic and trans-signaling.

Doxorubicin-induced chronic dilated cardiomyopathy-the apoptosis hypothesis revisited.

The chemotherapeutic agent doxorubicin (DOX) has significantly increased survival rates of pediatric and adult cancer patients. However, 10% of pediatric cancer survivors will 10-20 years later develop severe dilated cardiomyopathy (DCM), whereby the exact molecular mechanisms of disease progression after this long latency time remain puzzling. We here revisit the hypothesis that elevated apoptosis signaling or its increased likelihood after DOX exposure can lead to an impairment of cardiac function and cause a cardiac dilation. Based on recent literature evidence, we first argue why a dilated phenotype can occur when little apoptosis is detected. We then review findings suggesting that mature cardiomyocytes are protected against DOX-induced apoptosis downstream, but not upstream of mitochondrial outer membrane permeabilisation (MOMP). This lack of MOMP induction is proposed to alter the metabolic phenotype, induce hypertrophic remodeling, and lead to functional cardiac impairment even in the absence of cardiomyocyte apoptosis. We discuss findings that DOX exposure can lead to increased sensitivity to further cardiomyocyte apoptosis, which may cause a gradual loss in cardiomyocytes over time and a compensatory hypertrophic remodeling after treatment, potentially explaining the long lag time in disease onset. We finally note similarities between DOX-exposed cardiomyocytes and apoptosis-primed cancer cells and propose computational system biology as a tool to predict patient individual DOX doses. In conclusion, combining recent findings in rodent hearts and cardiomyocytes exposed to DOX with insights from apoptosis signal transduction allowed us to obtain a molecularly deeper insight in this delayed and still enigmatic pathology of DCM.

Signal transduction analysis of the NLRP3-inflammasome pathway after cellular damage and its paracrine regulation.

Activation of the NLRP3-inflammasome pathway and production of the inflammatory cytokine IL-1B after cellular damage caused by infarct or infection is a key process in several diseases such as acute myocardial infarction and inflammatory bowel disease. However, while the molecular triggers of the NLRP3-pathway after cellular damage are well known, the mechanisms that sustain or confine its activity are currently under investigation. We present here an Ordinary Differential Equation-based model that investigates the mechanisms of inflammasome activation and regulation in monocytes to predict IL-1ß activation kinetics upon a two-step activation by Damage-Associate-Molecular-Particles (DAMP) and extracellular ATP. Assuming both activation signals to be concomitantly present or present with a delay of 12h, the model predicted a transient IL-1ß activation at different concentration levels dependent on signal synchronisation. Introducing a positive feedback loop mediated by active IL-1ß resulted in a sustained IL-1ß activation, hence arguing for a paracrine signalling between inflammatory cells to guarantee a temporally stable inflammatory response. We then investigate mechanisms that control termination of inflammation using two recently identified molecular intervention points in the inflammasome pathway. We found that a more upstream regulation, by attenuating production of the IL-1ß-proform, was more potent in attenuating active IL-1ß production than direct inhibition of the NLRP3-inflammasome. Interestingly, ablating this upstream negative feedback led to a high variability of IL-1ß production in monocytes from different subjects, consistent with a recent pre-clinical study. We finally discuss the relevance and implications of our findings in disease models of acute myocardial infarction and spontaneous colitis.

Modelling tumour cell proliferation from vascular structure using tissue decomposition into avascular elements.

Computer models allow the mechanistically detailed study of tumour proliferation and its dependency on nutrients. However, the computational study of large vascular tumours requires detailed information on the 3-dimensional vessel network and rather high computation times due to complex geometries. This study puts forward the idea of partitioning vascularised tissue into connected avascular elements that can exchange cells and nutrients between each other. Our method is able to rapidly calculate the evolution of proliferating as well as dead and quiescent cells, and hence a proliferative index, from a given amount and distribution of vascularisation of arbitrary complexity. Applying our model, we found that a heterogeneous vessel distribution provoked a higher proliferative index, suggesting increased malignancy, and increased the amount of dead cells compared to a more static tumour environment when a homogenous vessel distribution was assumed. We subsequently demonstrated that under certain amounts of vascularisation, cell proliferation may even increase when vessel density decreases, followed by a subsequent decrease of proliferation. This effect was due to a trade-off between an increase in compensatory proliferation for replacing dead cells and a decrease of cell population due to lack of oxygen supply in lowly vascularised tumours. Findings were illustrated by an ectopic colorectal cancer mouse xenograft model. Our presented approach can be in the future applied to study the effect of cytostatic, cytotoxic and anti-angiogenic chemotherapy and is ideally suited for translational systems biology, where rapid interaction between theory and experiment is essential.

Computational Analysis of AMPK-Mediated Neuroprotection Suggests Acute Excitotoxic Bioenergetics and Glucose Dynamics Are Regulated by a Minimal Set of Critical Reactions.

Loss of ionic homeostasis during excitotoxic stress depletes ATP levels and activates the AMP-activated protein kinase (AMPK), re-establishing energy production by increased expression of glucose transporters on the plasma membrane. Here, we develop a computational model to test whether this AMPK-mediated glucose import can rapidly restore ATP levels following a transient excitotoxic insult. We demonstrate that a highly compact model, comprising a minimal set of critical reactions, can closely resemble the rapid dynamics and cell-to-cell heterogeneity of ATP levels and AMPK activity, as confirmed by single-cell fluorescence microscopy in rat primary cerebellar neurons exposed to glutamate excitotoxicity. The model further correctly predicted an excitotoxicity-induced elevation of intracellular glucose, and well resembled the delayed recovery and cell-to-cell heterogeneity of experimentally measured glucose dynamics. The model also predicted necrotic bioenergetic collapse and altered calcium dynamics following more severe excitotoxic insults. In conclusion, our data suggest that a minimal set of critical reactions may determine the acute bioenergetic response to transient excitotoxicity and that an AMPK-mediated increase in intracellular glucose may be sufficient to rapidly recover ATP levels following an excitotoxic insult.

Harnessing system models of cell death signalling for cytotoxic chemotherapy: towards personalised medicine approaches?

Most cytotoxic chemotherapeutics are believed to kill cancer cells by inducing apoptosis. Understanding the factors that contribute to impairment of apoptosis in cancer cells is therefore critical for the development of novel therapies that circumvent the widespread chemoresistance. Apoptosis, however, is a complex and tightly controlled process that can be induced by different classes of chemotherapeutics targeting different signalling nodes and pathways. Moreover, apoptosis initiation and apoptosis execution strongly depend on patient-specific, genomic and proteomic signatures. Here, we will review recent translational studies that suggest a critical link between the sensitivity of cancer cells to initiate apoptosis and clinical outcome. Next we will discuss recent advances in the field of system modelling of apoptosis pathways for the prediction of treatment responses. We propose that initiation of mitochondrial apoptosis, defined as the process of mitochondrial outer membrane permeabilisation (MOMP), is a dose-dependent decision process that allows for a prediction of individual therapy responses and therapeutic windows. We provide evidence in contrast that apoptosis execution post-MOMP may be a binary decision that dictates whether apoptosis is executed or not. We will discuss the implications of this concept for the future use of novel adjuvant therapeutics that specifically target apoptosis signalling pathways or which may be used to reduce the impact of cell-to-cell heterogeneity on therapy responses. Finally, we will discuss the technical and regulatory requirements surrounding the use and implications of system-based patient stratification tools for the future of personalised oncology.

Systems analysis of BCL2 protein family interactions establishes a model to predict responses to chemotherapy.

Apoptotic desensitization is a hallmark of cancer cells, but present knowledge of molecular systems controlling apoptosis has yet to provide significant prognostic insights. Here, we report findings from a systems study of the intrinsic pathway of apoptosis by BCL2 family proteins and clinical translation of its findings into a model with applications in colorectal cancer (CRC). By determining absolute protein quantifications in CRC cells and patient tumor samples, we found that BAK and BAX were expressed more highly than their antiapoptotic inhibitors. This counterintuitive finding suggested that sole inhibition of effector BAX and BAK could not be sufficient for systems stability in nonstressed cells. Assuming a model of direct effector activation by BH3-only proteins, we calculated that the amount of stress-induced BH3-only proteins required to activate mitochondrial apoptosis could predict individual death responses of CRC cells to 5-fluorouracil/oxaliplatin. Applying this model predictor to protein profiles in tumor and matched normal tissue samples from 26 patients with CRCs, we found that differences in protein quantities were sufficient to model the increased tumor sensitivity to chemotherapy compared with normal tissue. In addition, these differences were sufficient to differentiate clinical responders from nonresponders with high confidence. Applications of our model, termed DR_MOMP, were used to assess the impact of apoptosis-sensitizing drugs in lowering the necessary dose of state-of-the-art chemotherapy in individual patients. Together, our findings offer a ready clinical tool with the potential to tailor chemotherapy to individual patients.

Systems analysis of cancer cell heterogeneity in caspase-dependent apoptosis subsequent to mitochondrial outer membrane permeabilization.

Deregulation of apoptosis is a hallmark of carcinogenesis. We here combine live cell imaging and systems modeling to investigate caspase-dependent apoptosis execution subsequent to mitochondrial outer membrane permeabilization (MOMP) in several cancer cell lines. We demonstrate that, although most cell lines that underwent MOMP also showed robust and fast activation of executioner caspases and apoptosis, the colorectal cancer cell lines LoVo and HCT-116 Smac(-/-), similar to X-linked inhibitor of apoptosis protein (XIAP)-overexpressing HeLa (HeLa XIAP(Adv)) cells, only showed delayed and often no caspase activation, suggesting apoptosis impairment subsequent to MOMP. Employing APOPTO-CELL, a recently established model of apoptosis subsequent to MOMP, this impairment could be understood by studying the systemic interaction of five proteins that are present in the apoptosis pathway subsequent to MOMP. Using APOPTO-CELL as a tool to study detailed molecular mechanisms during apoptosis execution in individual cell lines, we demonstrate that caspase-9 was the most important regulator in DLD-1, HCT-116, and HeLa cells and identified additional cell line-specific co-regulators. Developing and applying a computational workflow for parameter screening, systems modeling identified that apoptosis execution kinetics are more robust against changes in reaction kinetics in HCT-116 and HeLa than in DLD-1 cells. Our systems modeling study is the first to draw attention to the variability in cell specific protein levels and reaction rates and to the emergent effects of such variability on the efficiency of apoptosis execution and on apoptosis impairment subsequent to MOMP.

Proteasome inhibition can impair caspase-8 activation upon submaximal stimulation of apoptotic tumor necrosis factor-related apoptosis inducing ligand (TRAIL) signaling.

Tumor necrosis factor-related apoptosis inducing ligand (TRAIL) can induce extrinsic apoptosis, resulting in caspase-8 activation, but may also initiate transcription-dependent prosurvival signaling. Proteasome inhibitors were suggested to promote TRAIL signal transduction through the death-inducing signaling complex (DISC) by modulating the relative abundance of core DISC components, thereby enhancing caspase-8 activation and apoptosis. To test this hypothesis, we quantified the changes in DISC protein levels as an early consequence of proteasome inhibition in HeLa cervical cancer cells and, based on these data, mathematically modeled the proapoptotic TRAIL signaling toward caspase-8 activation. Modeling results surprisingly suggested that caspase-8 activation might be delayed in presence of proteasome inhibitors, in particular at submaximal TRAIL doses. Subsequent FRET-based single cell time-lapse imaging at conditions where transcription dependent prosurvival signaling was blocked confirmed this hypothesis: caspase-8 activity was delayed by hours in the presence of proteasome inhibitors epoxomicin or bortezomib. Corresponding delays were detected for effector caspase processing and cell death. Contrary to current models, we therefore provide evidence that synergies between TRAIL and proteasome inhibitors do not result from changes in the levels of core DISC signaling proteins.

Positive feedback in the Akt/mTOR pathway and its implications for growth signal progression in skeletal muscle cells: an analytical study.

The IGF-1 mediated Akt/mTOR pathway has been recently proposed as mediator of skeletal muscle growth and a positive feedback between Akt and mTOR was suggested to induce homogeneous growth signals along the whole spatial extension of such long cells. Here we develop two biologically justified approximations which we study under the presence of four different initial conditions that describe different paradigms of IGF-1 receptor-induced Akt/mTOR activation. In first scenario the activation of the feedback cascade was assumed to be mild or protein turnover considered to be high. In turn, in the second scenario the transcriptional regulation was assumed to maintain defined levels of inactive pro-enzymes. For both scenarios, we were able to obtain closed-form formulas for growth signal progression in time and space and found that a localised initial signal maintains its Gaussian shape, but gets delocalised and exponentially degraded. Importantly, mathematical treatment of the reaction diffusion system revealed that diffusion filtered out high frequencies of spatially periodic initiator signals suggesting that the muscle cell is robust against fluctuations in spatial receptor expression or activation. However, neither scenario was consistent with the presence of stably travelling signal waves. Our study highlights the role of feedback loops in spatiotemporal signal progression and results can be applied to studies in cell proliferation, cell differentiation and cell death in other spatially extended cells.

A structured approach to the study of metabolic control principles in intact and impaired mitochondria.

We devised an approach to extract control principles of cellular bioenergetics for intact and impaired mitochondria from ODE-based models and applied it to a recently established bioenergetic model of cancer cells. The approach used two methods for varying ODE model parameters to determine those model components that, either alone or in combination with other components, most decisively regulated bioenergetic state variables. We found that, while polarisation of the mitochondrial membrane potential (ΔΨ(m)) and, therefore, the protomotive force were critically determined by respiratory complex I activity in healthy mitochondria, complex III activity was dominant for ΔΨ(m) during conditions of cytochrome-c deficiency. As a further important result, cellular bioenergetics in healthy, ATP-producing mitochondria was regulated by three parameter clusters that describe (1) mitochondrial respiration, (2) ATP production and consumption and (3) coupling of ATP-production and respiration. These parameter clusters resembled metabolic blocks and their intermediaries from top-down control analyses. However, parameter clusters changed significantly when cells changed from low to high ATP levels or when mitochondria were considered to be impaired by loss of cytochrome-c. This change suggests that the assumption of static metabolic blocks by conventional top-down control analyses is not valid under these conditions. Our approach is complementary to both ODE and top-down control analysis approaches and allows a better insight into cellular bioenergetics and its pathological alterations.