An official website of the United States government
Boolean networks are a popular modeling framework in computational biology to capture the dynamics of molecular networks, such as gene regulatory networks. It has been observed that many published models of such networks are defined by regulatory rules driving the dynamics that have certain so-called canalizing properties. In this paper, we investigate the dynamics of a random Boolean network with such properties using analytical methods and simulations. From our simulations, we observe that Boolean networks with higher canalizing depth have generally fewer attractors, the attractors are smaller, and the basins are larger, with implications for the stability and robustness of the models. These properties are relevant to many biological applications. Moreover, our results show that, from the standpoint of the attractor structure, high canalizing depth, compared to relatively small positive canalizing depth, has a very modest impact on dynamics. Motivated by these observations, we conduct mathematical study of the attractor structure of a random Boolean network of canalizing depth one (i.e., the smallest positive depth). For every positive integer ℓ , we give an explicit formula for the limit of the expected number of attractors of length ℓ in an n -state random Boolean network as n goes to infinity.
more »
« less
Parity and time reversal elucidate both decision-making in empirical models and attractor scaling in critical Boolean networks
We present new applications of parity inversion and time reversal to the emergence of complex behavior from simple dynamical rules in stochastic discrete models. Our parity-based encoding of causal relationships and time-reversal construction efficiently reveal discrete analogs of stable and unstable manifolds. We demonstrate their predictive power by studying decision-making in systems biology and statistical physics models. These applications underpin a novel attractor identification algorithm implemented for Boolean networks under stochastic dynamics. Its speed enables resolving a long-standing open question of how attractor count in critical random Boolean networks scales with network size and whether the scaling matches biological observations. Via 80-fold improvement in probed network size ( N = 16,384), we find the unexpectedly low scaling exponent of 0.12 ± 0.05, approximately one-tenth the analytical upper bound. We demonstrate a general principle: A system’s relationship to its time reversal and state-space inversion constrains its repertoire of emergent behaviors.
more »
« less
- Award ID(s):
- 1715826
- PAR ID:
- 10333725
- Date Published:
- Journal Name:
- Science Advances
- Volume:
- 7
- Issue:
- 29
- ISSN:
- 2375-2548
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Signaling networks are at the heart of almost all biological processes. Most of these networks contain large number of components, and often either the connections between these components are not known or the rate equations that govern the dynamics of soluble signaling components are not quantified. This uncertainty in network topology and parameters can make it challenging to formulate detailed mathematical models. Boolean networks, in which all components are either on or off, have emerged as viable alternatives to detailed mathematical models that contain rate constants and other parameters. Therefore, open-source platforms of Boolean models for community use are desirable. Here, we present Boolink, a freely available graphical user interface that allows users to easily construct and analyze existing Boolean networks. Boolink can be applied to any Boolean network. We demonstrate its application using a previously published network for abscisic acid (ABA)-driven stomatal closure in Arabidopsis spp. (Arabidopsis thaliana). We also show how Boolink can be used to generate testable predictions by extending the network to include CO2 regulation of stomatal movements. Predictions of the model were experimentally tested, and the model was iteratively modified based on experiments showing that ABA effectively closes Arabidopsis stomata at near-zero CO2 concentrations (1.5-ppm CO2). Thus, Boolink enables public generation and the use of existing Boolean models, including the prior developed ABA signaling model with added CO2 signaling components.more » « less
-
We propose a new algorithm for inference of gene regulatory networks (GRN) from noisy gene expression data based on maximum-likelihood (ML) adaptive filtering and the discrete fish school search algorithm (DFSS). The approach is based on the general partially-observed Boolean dynamical system (POBDS) model, and as such can be used for simultaneous state and parameter estimation for any Boolean dynamical system observed in noise. The proposed DFSS-ML-BKF algorithm combines the ML adaptive Boolean Kalman Filter (ML-BKF) with DFSS, a version of the Fish School Search algorithm tailored for discrete parameter spaces. Results based on synthetic gene expression time-series data using the well-known p53-MDM2 negative-feedback loop GRN demonstrate that DFSS-ML-BKF can infer the network topology accurately and efficiently.more » « less
-
Abstract Summarypystablemotifs is a Python 3 library for analyzing Boolean networks. Its non-heuristic and exhaustive attractor identification algorithm was previously presented in Rozum et al. (2021). Here, we illustrate its performance improvements over similar methods and discuss how it uses outputs of the attractor identification process to drive a system to one of its attractors from any initial state. We implement six attractor control algorithms, five of which are new in this work. By design, these algorithms can return different control strategies, allowing for synergistic use. We also give a brief overview of the other tools implemented in pystablemotifs. Availability and implementationThe source code is on GitHub at https://github.com/jcrozum/pystablemotifs/. Supplementary informationSupplementary data are available at Bioinformatics online.more » « less
-
We consider the problems of asymptotic stability and robustness in large-scale second-order consensus networks and vehicle platoons in the discrete-time domain. First, we develop a graph-theoretic methodology to design the state feedback law for the second-order consensus networks and vehicle platoons in a discrete-time framework. We analyze the stability of such networks based on algebraic properties of the Laplacian matrices of underlying graphs and each vehicle’s update cycle (also known as the time step). We further provide a necessary and sufficient condition of stability of a linear second-order consensus network in the discrete-time domain. Moreover, we evaluate the robustness of the consensus networks by employing the expected value of the steady-state dispersion of the state of the entire network, also known as squared H2-norm, as a performance measure. We show the connection between performance measures with respect to network size, connectivity, and the update cycle. The main contribution of this work is that we provide a formal framework to quantify the relation between scaling performance measures and restrictions of the vehicles’ update cycles. Specifically, we show that denser networks (i.e., networks with more communications/edges) require faster agents (i.e., smaller update cycles) to outperform or achieve the same level of robustness as sparse networks (i.e., networks with fewer communications/edges).more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1126/sciadv.abf8124
