Browsing by Author "Mischaikow, Konstantin"

Now showing 1 - 8 of 8

Combinatorial Representation of Parameter Space for Switching Networks
(2016-11) Cummins, Bree; Harker, Shaun; Mischaikow, Konstantin; Mok, Kafung; Gedeon, Tomas
We describe the theoretical and computational framework for the Dynamic Signatures Generated by Regulatory Networks (DSGRN) database. The motivation stems from an urgent need to understand the global dynamics of biologically relevant signal transduction/gene regulatory networks that have at least 5 to 10 nodes, involve multiple interactions, and have decades of parameters. The input to the database computations is a regulatory network, i.e., a directed graph with edges indicating up or down regulation. A computational model based on switching networks is generated from the regulatory network. The phase space dimension of this model equals the number of nodes and the associated parameter space consists of one parameter for each node (a decay rate) and three parameters for each edge (low level of expression, high level of expression, and threshold at which expression levels change). Since the nonlinearities of switching systems are piecewise constant, there is a natural decomposition of phase space into cells from which the dynamics can be described combinatorially in terms of a state transition graph. This in turn leads to a compact representation of the global dynamics called an annotated Morse graph that identifies recurrent and nonrecurrent dynamics. The focus of this paper is on the construction of a natural computable finite decomposition of parameter space into domains where the annotated Morse graph description of dynamics is constant. We use this decomposition to construct an SQL database that can be effectively searched for dynamical signatures such as bistability, stable or unstable oscillations, and stable equilibria. We include two simple 3-node networks to provide small explicit examples of the type of information stored in the DSGRN database. To demonstrate the computational capabilities of this system we consider a simple network associated with p53 that involves 5 nodes and a 29-dimensional parameter space.
DSGRN: Examining the Dynamics of Families of Logical Models
(2018-06) Cummins, Bree; Gedeon, Tomas; Harker, Shaun; Mischaikow, Konstantin
We present a computational tool DSGRN for exploring the dynamics of a network by computing summaries of the dynamics of switching models compatible with the network across all parameters. The network can arise directly from a biological problem, or indirectly as the interaction graph of a Boolean model. This tool computes a finite decomposition of parameter space such that for each region, the state transition graph that describes the coarse dynamical behavior of a network is the same. Each of these parameter regions corresponds to a different logical description of the network dynamics. The comparison of dynamics across parameters with experimental data allows the rejection of parameter regimes or entire networks as viable models for representing the underlying regulatory mechanisms. This in turn allows a search through the space of perturbations of a given network for networks that robustly fit the data. These are the first steps toward discovering a network that optimally matches the observed dynamics by searching through the space of networks.
Experimental guidance for discovering genetic networks through hypothesis reduction on time series
(Public Library of Science, 2022-10) Cummins, Breschine; Motta, Francis C.; Moseley, Robert C.; Deckard, Anastasia; Campione, Sophia; Gedeon, Tomáš; Mischaikow, Konstantin; Haase, Steven B.
Large programs of dynamic gene expression, like cell cyles and circadian rhythms, are controlled by a relatively small “core” network of transcription factors and post-translational modifiers, working in concerted mutual regulation. Recent work suggests that system-independent, quantitative features of the dynamics of gene expression can be used to identify core regulators. We introduce an approach of iterative network hypothesis reduction from time-series data in which increasingly complex features of the dynamic expression of individual, pairs, and entire collections of genes are used to infer functional network models that can produce the observed transcriptional program. The culmination of our work is a computational pipeline, Iterative Network Hypothesis Reduction from Temporal Dynamics (Inherent dynamics pipeline), that provides a priority listing of targets for genetic perturbation to experimentally infer network structure. We demonstrate the capability of this integrated computational pipeline on synthetic and yeast cell-cycle data.
Extending Combinatorial Regulatory Network Modeling to Include Activity Control and Decay Modulation
(Society for Industrial & Applied Mathematics, 2022-09) Cummins, Bree; Gameiro, Marcio; Gedeon, Tomas; Kepley, Shane; Mischaikow, Konstantin; Zhang, Lun
Understanding how the structure of within-system interactions affects the dynamics of the system is important in many areas of science. We extend a network dynamics modeling platform DSGRN, which combinatorializes both dynamics and parameter space to construct finite but accurate summaries of network dynamics, to new types of interactions. While the standard DSGRN assumes that each network edge controls the rate of abundance of the target node, the new edges may control either activity level or a decay rate of its target. While motivated by processes of post-transcriptional modification and ubiquitination in systems biology, our extension is applicable to the dynamics of any signed directed network.
Global dynamics for steep nonlinearities in two dimensions
(2017-01) Gedeon, Tomas; Harker, Shaun; Kokubu, Hiroshi; Mischaikow, Konstantin; Oka, Hiroe
This paper discusses a novel approach to obtaining mathematically rigorous results on the global dynamics of ordinary differential equations. We study switching models of regulatory networks. To each switching network we associate a Morse graph, a computable object that describes a Morse decomposition of the dynamics. In this paper we show that all smooth perturbations of the switching system share the same Morse graph and we compute explicit bounds on the size of the allowable perturbation. This shows that computationally tractable switching systems can be used to characterize dynamics of smooth systems with steep nonlinearities.
Global dynamics for switching systems and their extensions by linear differential equations
(2018-11) Huttinga, Zane; Cummins, Bree; Gedeon, Tomas; Mischaikow, Konstantin
Switching systems use piecewise constant nonlinearities to model gene regulatory networks. This choice provides advantages in the analysis of behavior and allows the global description of dynamics in terms of Morse graphs associated to nodes of a parameter graph. The parameter graph captures spatial characteristics of a decomposition of parameter space into domains with identical Morse graphs. However, there are many cellular processes that do not exhibit threshold-like behavior and thus are not well described by a switching system. We consider a class of extensions of switching systems formed by a mixture of switching interactions and chains of variables governed by linear differential equations. We show that the parameter graphs associated to the switching system and any of its extensions are identical. For each parameter graph node, there is an order-preserving map from the Morse graph of the switching system to the Morse graph of any of its extensions. We provide counterexamples that show why possible stronger relationships between the Morse graphs are not valid.
Rational design of complex phenotype via network models
(Public Library of Science, 2021-07) Gameiro, Marcio; Gedeon, Tomáš; Kepley, Shane; Mischaikow, Konstantin
We demonstrate a modeling and computational framework that allows for rapid screening of thousands of potential network designs for particular dynamic behavior. To illustrate this capability we consider the problem of hysteresis, a prerequisite for construction of robust bistable switches and hence a cornerstone for construction of more complex synthetic circuits. We evaluate and rank most three node networks according to their ability to robustly exhibit hysteresis where robustness is measured with respect to parameters over multiple dynamic phenotypes. Focusing on the highest ranked networks, we demonstrate how additional robustness and design constraints can be applied. We compare our results to more traditional methods based on specific parameterization of ordinary differential equation models and demonstrate a strong qualitative match at a small fraction of the computational cost.
Rational design of complex phenotype via network models
(Public Library of Science, 2021-07) Gameiro, Marcio; Gedeon, Tomáš; Kepley, Shane; Mischaikow, Konstantin
We demonstrate a modeling and computational framework that allows for rapid screening of thousands of potential network designs for particular dynamic behavior. To illustrate this capability we consider the problem of hysteresis, a prerequisite for construction of robust bistable switches and hence a cornerstone for construction of more complex synthetic circuits. We evaluate and rank most three node networks according to their ability to robustly exhibit hysteresis where robustness is measured with respect to parameters over multiple dynamic phenotypes. Focusing on the highest ranked networks, we demonstrate how additional robustness and design constraints can be applied. We compare our results to more traditional methods based on specific parameterization of ordinary differential equation models and demonstrate a strong qualitative match at a small fraction of the computational cost.