Amstein LK, Ackermann J, Hannig J, Đikić I, Fulda S, Koch I. Mathematical modeling of the molecular switch of TNFR1-mediated signaling pathways applying Petri net formalism and in silico knockout analysis.
PLoS Comput Biol 2022;
18:e1010383. [PMID:
35994517 PMCID:
PMC9467317 DOI:
10.1371/journal.pcbi.1010383]
[Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2021] [Revised: 09/12/2022] [Accepted: 07/13/2022] [Indexed: 11/19/2022] Open
Abstract
The paper describes a mathematical model of the molecular switches of cell survival, apoptosis, and necroptosis in cellular signaling pathways initiated by tumor necrosis factor 1. Based on experimental findings in the literature, we constructed a Petri net model based on detailed molecular reactions of the molecular players, protein complexes, post-translational modifications, and cross talk. The model comprises 118 biochemical entities, 130 reactions, and 299 edges. We verified the model by evaluating invariant properties of the system at steady state and by in silico knockout analysis. Applying Petri net analysis techniques, we found 279 pathways, which describe signal flows from receptor activation to cellular response, representing the combinatorial diversity of functional pathways.120 pathways steered the cell to survival, whereas 58 and 35 pathways led to apoptosis and necroptosis, respectively. For 65 pathways, the triggered response was not deterministic and led to multiple possible outcomes. We investigated the in silico knockout behavior and identified important checkpoints of the TNFR1 signaling pathway in terms of ubiquitination within complex I and the gene expression dependent on NF-κB, which controls the caspase activity in complex II and apoptosis induction. Despite not knowing enough kinetic data of sufficient quality, we estimated system’s dynamics using a discrete, semi-quantitative Petri net model.
It is still a challenge to develop mechanistic models for big molecular systems without the knowledge of enough kinetic parameters of sufficient quality. At the same time, more qualitative and semi-quantitative data have been produced in increasing numbers, e.g., by new high-throughput technologies. This has generated demands for new concepts at appropriate abstraction levels. The Petri net formalism enables the integration of qualitative as well as quantitative data and provides algorithms and methods for model verification and model simulation. Moreover, Petri nets exhibit a clear and coherent visualization. Here, we modeled the molecular switches between cell survival, apoptosis, and necroptosis induced by tumor necrosis factor 1. We were interested not only in an exhaustive exploration of all possible signaling pathways, but also in finding the system’s checkpoints. Our Petri net model comprises 118 biochemical entities, 130 reactions, and 299 edges. We found 279 pathways that describe signal flows from receptor activation to cellular response.120 pathways steered the cell to survival, whereas 58 and 35 pathways led to apoptosis and necroptosis, respectively. For 65 pathways, the triggered response was not deterministic, leading to multiple possible outcomes. We applied in silico knockout analyses to the Petri net model and could identify important checkpoints of the tumor necrosis factor 1 signaling pathway.
Collapse