skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: A compositional account of motifs, mechanisms, and dynamics in biochemical regulatory networks
Regulatory networks depict promoting or inhibiting interactions between molecules in a biochemical system. We introduce a category-theoretic formalism for regulatory networks, using signed graphs to model the networks and signed functors to describe occurrences of one network in another, especially occurrences of network motifs. With this foundation, we establish functorial mappings between regulatory networks and other mathematical models in biochemistry. We construct a functor from reaction networks, modeled as Petri nets with signed links, to regulatory networks, enabling us to precisely define when a reaction network could be a physical mechanism underlying a regulatory network. Turning to quantitative models, we associate a regulatory network with a Lotka-Volterra system of differential equations, defining a functor from the category of signed graphs to a category of parameterized dynamical systems. We extend this result from closed to open systems, demonstrating that Lotka-Volterra dynamics respects not only inclusions and collapsings of regulatory networks, but also the process of building up complex regulatory networks by gluing together simpler pieces. Formally, we use the theory of structured cospans to produce a lax double functor from the double category of open signed graphs to that of open parameterized dynamical systems. Throughout the paper, we ground the categorical formalism in examples inspired by systems biology.  more » « less
Award ID(s):
1839968
PAR ID:
10513253
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
Arxiv
Date Published:
Journal Name:
Compositionality
Volume:
6
ISSN:
2631-4444
Page Range / eLocation ID:
2
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; (, Physical Review E)
    We study the stochastic spatial Lotka-Volterra model for predator-prey interaction subject to a periodically varying carrying capacity. The Lotka-Volterra model with on-site lattice occupation restrictions (i.e., finite local carrying capacity) that represent finite food resources for the prey population exhibits a continuous active-to-absorbing phase transition. The active phase is sustained by the existence of spatio-temporal patterns in the form of pursuit and evasion waves. Monte Carlo simulations on a two-dimensional lattice are utilized to investigate the effect of seasonal variations of the environment on species coexistence. The results of our simulations are also compared to a mean-field analysis in order to specifically delineate the impact of stochastic fluctuations and spatial correlations. We find that the parameter region of predator and prey coexistence is enlarged relative to the stationary situation when the carrying capacity varies periodically. The (quasi-)stationary regime of our periodically varying Lotka-Volterra predator-prey system shows qualitative agreement between the stochastic model and the mean-field approximation. However, under periodic carrying capacity-switching environments, the mean-field rate equations predict period-doubling scenarios that are washed out by internal reaction noise in the stochastic lattice model. Utilizing visual representations of the lattice simulations and dynamical correlation functions, we study how the pursuit and evasion waves are affected by ensuing resonance effects. Correlation function measurements indicate a time delay in the response of the system to sudden changes in the environment. Resonance features are observed in our simulations that cause prolonged persistent spatial correlations. Different effective static environments are explored in the extreme limits of fast and slow periodic switching. The analysis of the mean-field equations in the fast-switching regime enables a semi-quantitative description of the (quasi-)stationary state. 
    more » « less
  2. ; ; ; ; ; ; ; ; ; (, PMLR)
    Complex systems are characterized by intricate interactions between entities that evolve dynamically over time. Accurate inference of these dynamic relationships is crucial for understanding and predicting system behavior. In this paper, we propose Regulatory Temporal Interaction Network Inference (RiTINI) for inferring time-varying interaction graphs in complex systems using a novel combination of space-and-time graph attentions and graph neural ordinary differential equations (ODEs). RiTINI leverages time-lapse signals on a graph prior, as well as perturbations of signals at various nodes in order to effectively capture the dynamics of the underlying system. This approach is distinct from traditional causal inference networks, which are limited to inferring acyclic and static graphs. In contrast, RiTINI can infer cyclic, directed, and time-varying graphs, providing a more comprehensive and accurate representation of complex systems. The graph attention mechanism in RiTINI allows the model to adaptively focus on the most relevant interactions in time and space, while the graph neural ODEs enable continuous-time modeling of the system’s dynamics. We evaluate RiTINI’s performance on simulations of dynamical systems, neuronal networks, and gene regulatory networks, demonstrating its state-of-the-art capability in inferring interaction graphs compared to previous methods. 
    more » « less
  3. ; ; (, Inverse Problems)
    Abstract Our recent work on linear and affine dynamical systems has laid out a general framework for inferring the parameters of a differential equation model from a discrete set of data points collected from a system being modeled. It introduced a new class of inverse problems where qualitative information about the parameters and the associated dynamics of the system is determined for regions of the data space, rather than just for isolated experiments. Rigorous mathematical results have justified this approach and have identified common features that arise for certain classes of integrable models. In this work we present a thorough numerical investigation that shows that several of these core features extend to a paradigmatic linear-in-parameters model, the Lotka–Volterra (LV) system, which we consider in the conservative case as well as under the addition of terms that perturb the system away from this regime. A central construct for this analysis is a concise representation of parameter and dynamical features in the data space that we call thePn-diagram, which is particularly useful for visualization of the qualitative dependence of the system dynamics on data for low-dimensional (smalln) systems. Our work also exposes some new properties related to non-uniqueness that arise for these LV systems, with non-uniqueness manifesting as a multi-layered structure in the associatedP2-diagrams. 
    more » « less
  4. ; ; (, Proceedings of the ... American Control Conference)
    A hybrid control architecture for nonlinear dynamical systems is described which combines the advantages of model based control with those of real-time learning. The idea is to generate input-output data from an error system involving the plant and a proposed model. A discretized Chen-Fliess functional series is then identified from this data and used in conjunction with the model for predictive control. This method builds on the authors’ previous work on model-free control of a single-input, single-output Lotka-Volterra system.The problem is revisited here, but now with the introduction of a model for the dynamics. The single-input, multiple-output version of the problem is also investigated as a way to enhance closed-loop performance 
    more » « less
  5. ; ; (, Journal of Animal Ecology)
    Abstract Resilience is broadly understood as the ability of an ecological system to resist and recover from perturbations acting on species abundances and on the system's structure. However, one of the main problems in assessing resilience is to understand the extent to which measures of recovery and resistance provide complementary information about a system. While recovery from abundance perturbations has a strong tradition under the analysis of dynamical stability, it is unclear whether this same formalism can be used to measure resistance to structural perturbations (e.g. perturbations to model parameters).Here, we provide a framework grounded on dynamical and structural stability in Lotka–Volterra systems to link recovery from small perturbations on species abundances (i.e. dynamical indicators) with resistance to parameter perturbations of any magnitude (i.e. structural indicators). We use theoretical and experimental multispecies systems to show that the faster the recovery from abundance perturbations, the higher the resistance to parameter perturbations.We first use theoretical systems to show that the return rate along the slowest direction after a small random abundance perturbation (what we call full recovery) is negatively correlated with the largest random parameter perturbation that a system can withstand before losing any species (what we call full resistance). We also show that the return rate along the second fastest direction after a small random abundance perturbation (what we call partial recovery) is negatively correlated with the largest random parameter perturbation that a system can withstand before at most one species survives (what we call partial resistance). Then, we use a dataset of experimental microbial systems to confirm our theoretical expectations and to demonstrate that full and partial components of resilience are complementary.Our findings reveal that we can obtain the same level of information about resilience by measuring either a dynamical (i.e. recovery) or a structural (i.e. resistance) indicator. Irrespective of the chosen indicator (dynamical or structural), our results show that we can obtain additional information by separating the indicator into its full and partial components. We believe these results can motivate new theoretical approaches and empirical analyses to increase our understanding about risk in ecological systems. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email