The attractors of Boolean networks and their basins have been shown to be highly relevant for model validation and predictive modeling, e.g., in systems biology. Yet, there are currently very few tools available that are able to compute and visualize not only attractors but also their basins. In the realm of asynchronous, non-deterministic modeling not only is the repertoire of software even more limited, but also the formal notions for basins of attraction are often lacking. In this setting, the difficulty both for theory and computation arises from the fact that states may be elements of several distinct basins. In this paper, we address this topic by partitioning the state space into sets that are committed to the same attractors. These commitment sets can easily be generalized to sets that are equivalent w.r.t. the long-term behaviors of pre-selected nodes which leads us to the notions of markers and phenotypes which we illustrate in a case study on bladder tumorigenesis. For every concept, we propose equivalent CTL model checking queries and an extension of the state of the art model checking software NuSMV is made available that is capable of computing the respective sets. All notions are fully integrated as three new modules in our Python package PyBoolNet, including functions for visualizing the basins, commitment sets, and phenotypes as quotient graphs and pie charts.

Motivation: The goal of this project is to provide a simple interface to working with Boolean networks. Emphasis is put on easy access to a large number of common tasks including the generation and manipulation of networks, attractor and basin computation, model checking and trap space computation, execution of established graph algorithms as well as graph drawing and layouts. Results: P y B ool N et is a Python package for working with Boolean networks that supports simple access to model checking via N u SMV, standard graph algorithms via N etwork X and visualization via dot . In addition, state of the art attractor computation exploiting P otassco ASP is implemented. The package is function-based and uses only native Python and N etwork X data types. Availability and Implementation: https://github.com/hklarner/PyBoolNet. Contact: [email protected].

This paper is concerned with the analysis of labeled Thomas networks using discrete time series. It focuses on refining the given edge labels and on assessing the data quality. The results are aimed at being exploitable for experimental design and include the prediction of new activatory or inhibitory effects of given interactions and yet unobserved oscillations of specific components in between specific sampling intervals. On the formal side, we generalize the concept of edge labels and introduce a discrete time series interpretation. This interpretation features two original concepts: 1) Incomplete measurements are admissible, and 2) it allows qualitative assumptions about the changes in gene expression by means of monotonicity. On the computational side, we provide a Python script, erda.py, that automates the suggested workflow by model checking and constraint satisfaction. We illustrate the workflow by investigating the yeast network IRMA.