-
Classifier construction in Boolean networks using algebraic methods
Authors:
Robert Schwieger,
MatÃas R. Bender,
Heike Siebert,
Christian Haase
Abstract:
We investigate how classifiers for Boolean networks (BNs) can be constructed and modified under constraints. A typical constraint is to observe only states in attractors or even more specifically steady states of BNs. Steady states of BNs are one of the most interesting features for application. Large models can possess many steady states. In the typical scenario motivating this paper we start fro…
▽ More
We investigate how classifiers for Boolean networks (BNs) can be constructed and modified under constraints. A typical constraint is to observe only states in attractors or even more specifically steady states of BNs. Steady states of BNs are one of the most interesting features for application. Large models can possess many steady states. In the typical scenario motivating this paper we start from a Boolean model with a given classification of the state space into phenotypes defined by high-level readout components. In order to link molecular biomarkers with experimental design, we search for alternative components suitable for the given classification task. This is useful for modelers of regulatory networks for suggesting experiments and measurements based on their models. It can also help to explain causal relations between components and phenotypes. To tackle this problem we need to use the structure of the BN and the constraints. This calls for an algebraic approach. Indeed we demonstrate that this problem can be reformulated into the language of algebraic geometry. While already interesting in itself, this allows us to use Groebner bases to construct an algorithm for finding such classifiers. We demonstrate the usefulness of this algorithm as a proof of concept on a model with 25 components.
△ Less
Submitted 19 August, 2021;
originally announced August 2021.
-
Representing Model Ensembles as Boolean Functions
Authors:
Robert Schwieger,
Heike Siebert
Abstract:
Families of ODE models characterized by a common sign structure of their Jacobi matrix are investigated within the formalism of qualitative differential equations. In the context of regulatory networks the sign structure of the Jacobi matrix carries the information about which components of the network inhibit or activate each other. Information about constraints on the behavior of models in this…
▽ More
Families of ODE models characterized by a common sign structure of their Jacobi matrix are investigated within the formalism of qualitative differential equations. In the context of regulatory networks the sign structure of the Jacobi matrix carries the information about which components of the network inhibit or activate each other. Information about constraints on the behavior of models in this family is stored in a so called qualitative state transition graph. We showed previously that a similar approach can be used to analyze a model pool of Boolean functions characterized by a common interaction graph. Here we show that the opposite approach is fruitful as well. We show that the qualitative state transition graph can be reduced to a "skeleton" represented by a Boolean function conserving the reachability properties. This reduction has the advantage that approaches such as model checking and network inference methods can be applied to the "skeleton" within the framework of Boolean networks. Furthermore, our work constitutes an alternative to approach to link Boolean networks and differential equations.
△ Less
Submitted 2 August, 2018;
originally announced August 2018.
-
Basins of Attraction, Commitment Sets and Phenotypes of Boolean Networks
Authors:
Hannes Klarner,
Frederike Heinitz,
Sarah Nee,
Heike Siebert
Abstract:
The attractors of Boolean networks and their basins have been shown to be highly relevant for model validation and predictive modelling, e.g., in systems biology. Yet there are currently very few tools available that are able to compute and visualise not only attractors but also their basins. In the realm of asynchronous, non-deterministic modeling not only is the repertoire of software even more…
▽ More
The attractors of Boolean networks and their basins have been shown to be highly relevant for model validation and predictive modelling, e.g., in systems biology. Yet there are currently very few tools available that are able to compute and visualise 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 ele- ments 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 generalised to sets that are equivalent w.r.t. the long-term behaviours 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 capa- ble of computing the respective sets. All notions are fully integrated as three new modules in our Python package PyBoolNet, including functions for visualising the basins, commitment sets and phenotypes as quotient graphs and pie charts.
△ Less
Submitted 26 July, 2018;
originally announced July 2018.
-
Computing Maximal and Minimal Trap Spaces of Boolean Networks
Authors:
Hannes Klarner,
Alexander Bockmayr,
Heike Siebert
Abstract:
Asymptotic behaviors are often of particular interest when analyzing Boolean networks that represent biological systems such as signal trans- duction or gene regulatory networks. Methods based on a generalization of the steady state notion, the so-called trap spaces, can be exploited to investigate attractor properties as well as for model reduction techniques. In this paper, we propose a novel op…
▽ More
Asymptotic behaviors are often of particular interest when analyzing Boolean networks that represent biological systems such as signal trans- duction or gene regulatory networks. Methods based on a generalization of the steady state notion, the so-called trap spaces, can be exploited to investigate attractor properties as well as for model reduction techniques. In this paper, we propose a novel optimization-based method for com- puting all minimal and maximal trap spaces and motivate their use. In particular, we add a new result yielding a lower bound for the number of cyclic attractors and illustrate the methods with a study of a MAPK pathway model. To test the efficiency and scalability of the method, we compare the performance of the ILP solver gurobi with the ASP solver potassco in a benchmark of random networks.
△ Less
Submitted 25 September, 2015;
originally announced September 2015.
-
Minimization and Equivalence in Multi-valued Logical Models of Regulatory Networks
Authors:
Adam Streck,
Therese Lorenz,
Heike Siebert
Abstract:
Multi-valued logical models can be used to describe biological networks on a high level of abstraction based on the network structure and logical parameters capturing regulatory effects. Interestingly, the dynamics of two distinct models need not necessarily be different, which might hint at either only non-functional characteristics distinguishing the models or at different possible implementatio…
▽ More
Multi-valued logical models can be used to describe biological networks on a high level of abstraction based on the network structure and logical parameters capturing regulatory effects. Interestingly, the dynamics of two distinct models need not necessarily be different, which might hint at either only non-functional characteristics distinguishing the models or at different possible implementations for the same behaviour. Here, we study the conditions allowing for such effects by analysing classes of dynamically equivalent models and both structurally maximal and minimal representatives of such classes. Finally, we present an efficient algorithm that constructs a minimal representative of the respective class of a given multi-valued model.
△ Less
Submitted 27 January, 2016; v1 submitted 23 September, 2015;
originally announced September 2015.