A model of the onset of the senescence associated secretory phenotype after DNA damage induced senescence.

Cells and tissues are exposed to stress from numerous sources. Senescence is a protective mechanism that prevents malignant tissue changes and constitutes a fundamental mechanism of aging. It can be accompanied by a senescence associated secretory phenotype (SASP) that causes chronic inflammation. We present a Boolean network model-based gene regulatory network of the SASP, incorporating published gene interaction data. The simulation results describe current biological knowledge. The model predicts different in-silico knockouts that prevent key SASP-mediators, IL-6 and IL-8, from getting activated upon DNA damage. The NF-κB Essential Modulator (NEMO) was the most promising in-silico knockout candidate and we were able to show its importance in the inhibition of IL-6 and IL-8 following DNA-damage in murine dermal fibroblasts in-vitro. We strengthen the speculated regulator function of the NF-κB signaling pathway in the onset and maintenance of the SASP using in-silico and in-vitro approaches. We were able to mechanistically show, that DNA damage mediated SASP triggering of IL-6 and IL-8 is mainly relayed through NF-κB, giving access to possible therapy targets for SASP-accompanied diseases.

ViSiBooL-visualization and simulation of Boolean networks with temporal constraints.

Summary: Mathematical models and their simulation are increasingly used to gain insights into cellular pathways and regulatory networks. Dynamics of regulatory factors can be modeled using Boolean networks (BNs), among others. Text-based representations of models are precise descriptions, but hard to understand and interpret. ViSiBooL aims at providing a graphical way of modeling and simulating networks. By providing visualizations of static and dynamic network properties simultaneously, it is possible to directly observe the effects of changes in the network structure on the behavior. In order to address the challenges of clear design and a user-friendly graphical user interface (GUI), ViSiBooL implements visual representations of BNs. Additionally temporal extensions of the BNs for the modeling of regulatory time delays are incorporated. The GUI of ViSiBooL allows to model, organize, simulate and visualize BNs as well as corresponding simulation results such as attractors. Attractor searches are performed in parallel to the modeling process. Hence, changes in the network behavior are visualized at the same time. Availability and Implementation: ViSiBooL (Java 8) is freely available at http://sysbio.uni-ulm.de/?Software:ViSiBooL . Contact: hans.kestler@uni-ulm.de.

Boolean modeling identifies Greatwall/MASTL as an important regulator in the AURKA network of neuroblastoma.

Aurora Kinase A (AURKA) is often overexpressed in neuroblastoma (NB) with poor outcome. The causes of AURKA overexpression in NB are unknown. Here, we describe a gene regulatory network consisting of core regulators of AURKA protein expression and activation during mitosis to identify potential causes. This network was transformed to a dynamic Boolean model. Simulated activation of the serine/threonine protein kinase Greatwall (GWL, encoded by MASTL) that attenuates the pivotal AURKA inhibitor PP2A, predicted stabilization of AURKA. Consistent with this notion, gene set enrichment analysis showed enrichment of mitotic spindle assembly genes and MYCN target genes in NB with high GWL/MASTL expression. In line with the prediction of GWL/MASTL enhancing AURKA, elevated expression of GWL/MASTL was associated with NB risk factors and poor survival of patients. These results establish Boolean network modeling of oncogenic pathways in NB as a useful means for guided discovery in this enigmatic cancer.

Predicting Variabilities in Cardiac Gene Expression with a Boolean Network Incorporating Uncertainty.

Gene interactions in cells can be represented by gene regulatory networks. A Boolean network models gene interactions according to rules where gene expression is represented by binary values (on / off or {1, 0}). In reality, however, the gene's state can have multiple values due to biological properties. Furthermore, the noisy nature of the experimental design results in uncertainty about a state of the gene. Here we present a new Boolean network paradigm to allow intermediate values on the interval [0, 1]. As in the Boolean network, fixed points or attractors of such a model correspond to biological phenotypes or states. We use our new extension of the Boolean network paradigm to model gene expression in first and second heart field lineages which are cardiac progenitor cell populations involved in early vertebrate heart development. By this we are able to predict additional biological phenotypes that the Boolean model alone is not able to identify without utilizing additional biological knowledge. The additional phenotypes predicted by the model were confirmed by published biological experiments. Furthermore, the new method predicts gene expression propensities for modelled but yet to be analyzed genes.