Abstract

This paper introduces the notion of approximating asynchronous attractors of Boolean networks by minimal trap spaces. We define three criteria for determining the quality of an approximation: “faithfulness” which requires that the oscillating variables of all attractors in a trap space correspond to their dimensions, “univocality” which requires that there is a unique attractor in each trap space, and “completeness” which requires that there are no attractors outside of a given set of trap spaces. Each is a reachability property for which we give equivalent model checking queries. Whereas faithfulness and univocality can be decided by model checking the corresponding subnetworks, the naive query for completeness must be evaluated on the full state space. Our main result is an alternative approach which is based on the iterative refinement of an initially poor approximation. The algorithm detects so-called autonomous sets in the interaction graph, variables that contain all their regulators, and considers their intersection and extension in order to perform model checking on the smallest possible state spaces. A benchmark, in which we apply the algorithm to 18 published Boolean networks, is given. In each case, the minimal trap spaces are faithful, univocal, and complete, which suggests that they are in general good approximations for the asymptotics of Boolean networks.

Highlights

  • Boolean and multi-valued networks are frequently used to model the dynamics of biological processes that involve gene regulation and signal transduction

  • We developed the notion of an approximation of attractors of a Boolean network

  • Minimal trap spaces are approximations that can be computed for networks with hundreds of variables using Answer Set Programing (ASP) solvers

Read more

Summary

Introduction

Boolean and multi-valued networks are frequently used to model the dynamics of biological processes that involve gene regulation and signal transduction. The dynamics of such models is captured by the state transition graph, a directed graph that relates states to potential successor states. The long-term behaviors correspond to the minimal trap sets of the state transition graph which are called its attractors. We have suggested to compute the minimal trap spaces of a network to obtain an approximation for its cyclic attractors (Klarner et al, 2014) and proposed an efficient, Answer Set Programing (ASP)-based method for their computation. There is a Supplementary Material that contains proofs for the formal statements in the main text

Background
The Attractor Detection Problem
Approximating Attractors by Subspaces
Completeness
Deciding Completeness by Iterative Refinement
Results
Conclusion and Outlook
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call