Abstract
Cell-to-cell variability during TNFα stimulated Tumor Necrosis Factor Receptor 1 (TNFR1) signaling can lead to single-cell level pro-survival and apoptotic responses. This variability stems from the heterogeneity in signal flow through intracellular signaling entities that regulate the balance between these two phenotypes. Using systematic Boolean dynamic modeling of a TNFR1 signaling network, we demonstrate that the signal flow path variability can be modulated to enable cells favour apoptosis. We developed a computationally efficient approach "Boolean Modeling based Prediction of Steady-state probability of Phenotype Reachability (BM-ProSPR)" to accurately predict the network's ability to settle into different phenotypes. Model analysis juxtaposed with the experimental observations revealed that NFκB and PI3K transient responses guide the XIAP behaviour to coordinate the crucial dynamic cross-talk between the pro-survival and apoptotic arms at the single-cell level. Model predicted the experimental observations that ~31% apoptosis increase can be achieved by arresting Comp1 - IKK* activity which regulates the NFκB and PI3K dynamics. Arresting Comp1 - IKK* activity causes signal flow path re-wiring towards apoptosis without significantly compromising NFκB levels, which govern adequate cell survival. Priming an ensemble of cancerous cells with inhibitors targeting the specific interaction involving Comp1 and IKK* prior to TNFα exposure could enable driving them towards apoptosis.