Homeostasis in a feed forward loop gene regulatory motif

Dynamics and bifurcations in genetic circuits with fibration symmetries

Circuit building blocks of gene regulatory networks (GRN) have been identified through the fibration symmetries of the underlying biological graph. Here, we analyse analytically six of these circuits that occur as functional and synchronous building blocks in these networks. Of these, the lock-on, toggle switch, Smolen oscillator, feed-forward fibre and Fibonacci fibre circuits occur in living organisms, notably Escherichia coli; the sixth, the repressilator, is a synthetic GRN. We consider synchronous steady states determined by a fibration symmetry (or balanced colouring) and determine analytic conditions for local bifurcation from such states, which can in principle be either steady-state or Hopf bifurcations. We identify conditions that characterize the first bifurcation, the only one that can be stable near the bifurcation point. We model the state of each gene in terms of two variables: mRNA and protein concentration. We consider all possible 'admissible' models-those compatible with the network structure-and then specialize these general results to simple models based on Hill functions and linear degradation. The results systematically classify using graph symmetries the complexity and dynamics of these circuits, which are relevant to understand the functionality of natural and synthetic cells.

Classification of 2-node excitatory-inhibitory networks

We classify connected 2-node excitatory-inhibitory networks under various conditions. We assume that, as well as for connections, there are two distinct node-types, excitatory and inhibitory. In our classification we consider four different types of excitatory-inhibitory networks: restricted, partially restricted, unrestricted and completely unrestricted. For each type we give two different classifications. Using results on ODE-equivalence and minimality, we classify the ODE-classes and present a minimal representative for each ODE-class. We also classify all the networks with valence ≤2. These classifications are up to renumbering of nodes and the interchange of 'excitatory' and 'inhibitory' on nodes and arrows. These classifications constitute a first step towards analysing dynamics and bifurcations of excitatory-inhibitory networks. The results have potential applications to biological network models, especially neuronal networks, gene regulatory networks, and synthetic gene networks.

Homeostasis in networks with multiple inputs

Homeostasis, also known as adaptation, refers to the ability of a system to counteract persistent external disturbances and tightly control the output of a key observable. Existing studies on homeostasis in network dynamics have mainly focused on 'perfect adaptation' in deterministic single-input single-output networks where the disturbances are scalar and affect the network dynamics via a pre-specified input node. In this paper we provide a full classification of all possible network topologies capable of generating infinitesimal homeostasis in arbitrarily large and complex multiple inputs networks. Working in the framework of 'infinitesimal homeostasis' allows us to make no assumption about how the components are interconnected and the functional form of the associated differential equations, apart from being compatible with the network architecture. Remarkably, we show that there are just three distinct 'mechanisms' that generate infinitesimal homeostasis. Each of these three mechanisms generates a rich class of well-defined network topologies-called homeostasis subnetworks. More importantly, we show that these classes of homeostasis subnetworks provides a topological basis for the classification of 'homeostasis types': the full set of all possible multiple inputs networks can be uniquely decomposed into these special homeostasis subnetworks. We illustrate our results with some simple abstract examples and a biologically realistic model for the co-regulation of calcium ( Ca ) and phosphate ( PO 4 ) in the rat. Furthermore, we identify a new phenomenon that occurs in the multiple input setting, that we call homeostasis mode interaction, in analogy with the well-known characteristic of multiparameter bifurcation theory.

Classification of infinitesimal homeostasis in four-node input-output networks

An input-output network has an input node [Formula: see text], an output node o, and regulatory nodes [Formula: see text]. Such a network is a core network if each [Formula: see text] is downstream from [Formula: see text] and upstream from o. Wang et al. (J Math Biol 82:62, 2021. https://doi.org/10.1007/s00285-021-01614-1 ) show that infinitesimal homeostasis can be classified in biochemical networks through infinitesimal homeostasis in core subnetworks. Golubitsky and Wang (J Math Biol 10:1-23, 2020) show that there are three types of 3-node core networks and three types of infinitesimal homeostasis in 3-node core networks. This paper uses the theory developed in Wang et al. (2021) to show that there are twenty types of 4-node core networks (Theorem 1.3) and seventeen types of infinitesimal homeostasis in 4-node core networks (Theorem 1.7). Biological contexts illustrate the classification theorems and show that the theory can be an aid when calculating homeostasis in specific biochemical networks.

A homeostasis criterion for limit cycle systems based on infinitesimal shape response curves

Homeostasis occurs in a control system when a quantity remains approximately constant as a parameter, representing an external perturbation, varies over some range. Golubitsky and Stewart (J Math Biol 74(1-2):387-407, 2017) developed a notion of infinitesimal homeostasis for equilibrium systems using singularity theory. Rhythmic physiological systems (breathing, locomotion, feeding) maintain homeostasis through control of large-amplitude limit cycles rather than equilibrium points. Here we take an initial step to study (infinitesimal) homeostasis for limit-cycle systems in terms of the average of a quantity taken around the limit cycle. We apply the "infinitesimal shape response curve" (iSRC) introduced by Wang et al. (SIAM J Appl Dyn Syst 82(7):1-43, 2021) to study infinitesimal homeostasis for limit-cycle systems in terms of the mean value of a quantity of interest, averaged around the limit cycle. Using the iSRC, which captures the linearized shape displacement of an oscillator upon a static perturbation, we provide a formula for the derivative of the averaged quantity with respect to the control parameter. Our expression allows one to identify homeostasis points for limit cycle systems in the averaging sense. We demonstrate in the Hodgkin-Huxley model and in a metabolic regulatory network model that the iSRC-based method provides an accurate representation of the sensitivity of averaged quantities.

The structure of infinitesimal homeostasis in input-output networks

Homeostasis refers to a phenomenon whereby the output [Formula: see text] of a system is approximately constant on variation of an input [Formula: see text]. Homeostasis occurs frequently in biochemical networks and in other networks of interacting elements where mathematical models are based on differential equations associated to the network. These networks can be abstracted as digraphs [Formula: see text] with a distinguished input node [Formula: see text], a different distinguished output node o, and a number of regulatory nodes [Formula: see text]. In these models the input-output map [Formula: see text] is defined by a stable equilibrium [Formula: see text] at [Formula: see text]. Stability implies that there is a stable equilibrium [Formula: see text] for each [Formula: see text] near [Formula: see text] and infinitesimal homeostasis occurs at [Formula: see text] when [Formula: see text]. We show that there is an [Formula: see text] homeostasis matrix [Formula: see text] for which [Formula: see text] if and only if [Formula: see text]. We note that the entries in H are linearized couplings and [Formula: see text] is a homogeneous polynomial of degree [Formula: see text] in these entries. We use combinatorial matrix theory to factor the polynomial [Formula: see text] and thereby determine a menu of different types of possible homeostasis associated with each digraph [Formula: see text]. Specifically, we prove that each factor corresponds to a subnetwork of [Formula: see text]. The factors divide into two combinatorially defined classes: structural and appendage. Structural factors correspond to feedforward motifs and appendage factors correspond to feedback motifs. Finally, we discover an algorithm for determining the homeostasis subnetwork motif corresponding to each factor of [Formula: see text] without performing numerical simulations on model equations. The algorithm allows us to classify low degree factors of [Formula: see text]. There are two types of degree 1 homeostasis (negative feedback loops and kinetic or Haldane motifs) and there are two types of degree 2 homeostasis (feedforward loops and a degree two appendage motif).

Infinitesimal homeostasis in three-node input-output networks

Homeostasis occurs in a system where an output variable is approximately constant on an interval on variation of an input variable [Formula: see text]. Homeostasis plays an important role in the regulation of biological systems, cf.  Ferrell (Cell Syst 2:62-67, 2016), Tang and McMillen (J Theor Biol 408:274-289, 2016), Nijhout et al. (BMC Biol 13:79, 2015), and Nijhout et al. (Wiley Interdiscip Rev Syst Biol Med 11:e1440, 2018). A method for finding homeostasis in mathematical models is given in the control theory literature as points where the derivative of the output variable with respect to [Formula: see text] is identically zero. Such points are called perfect homeostasis or perfect adaptation. Alternatively, Golubitsky and Stewart (J Math Biol 74:387-407, 2017) use an infinitesimal notion of homeostasis (namely, the derivative of the input-output function is zero at an isolated point) to introduce singularity theory into the study of homeostasis. Reed et al. (Bull Math Biol 79(9):1-24, 2017) give two examples of infinitesimal homeostasis in three-node chemical reaction systems: feedforward excitation and substrate inhibition. In this paper we show that there are 13 different three-node networks leading to 78 three-node input-output network configurations, under the assumption that there is one input node, one output node, and they are distinct. The different configurations are based on which node is the input node and which node is the output node. We show nonetheless that there are only three basic mechanisms for three-node input-output networks that lead to infinitesimal homeostasis and we call them structural homeostasis, Haldane homeostasis, and null-degradation homeostasis. Substantial parts of this classification are given in Ma et al. (Cell 138:760-773, 2009) and Ferrell (2016) among others. Our contributions include giving a complete classification using general admissible systems (Golubitsky and Stewart in Bull Am Math Soc 43:305-364, 2006) rather than specific biochemical models, relating the types of infinitesimal homeostasis to the graph theoretic existence of simple paths, and providing the basis to use singularity theory to study higher codimension homeostasis singularities such as the chair singularities introduced in Nijhout and Reed (Integr Comp Biol 54(2):264-275, 2014. https://doi.org/10.1093/icb/icu010) and Nijhout et al. (Math Biosci 257:104-110, 2014). See Golubitsky and Stewart (2017). The first two of these mechanisms are illustrated by feedforward excitation and substrate inhibition. Structural homeostasis occurs only when the network has a feedforward loop as a subnetwork; that is, when there are two distinct simple paths connecting the input node to the output node. Moreover, when the network is just the feedforward loop motif itself, one of the paths must be excitatory and one inhibitory to support infinitesimal homeostasis. Haldane homeostasis occurs when there is a single simple path from the input node to the output node and then only when one of the couplings along this path has strength 0. Null-degradation homeostasis is illustrated by a biochemical example from Ma et al. (2009); this kind of homeostasis can occur only when the degradation constant of the third node is 0. The paper ends with an analysis of Haldane homeostasis infinitesimal chair singularities.

Large-Scale Functional Analysis of CRP-Mediated Feed-Forward Loops

Feed-forward loops (FFLs) represent an important and basic network motif to understand specific biological functions. Cyclic-AMP (cAMP) receptor protein (CRP), a transcription factor (TF), mediates catabolite repression and regulates more than 400 genes in response to changes in intracellular concentrations of cAMP in Escherichia coli. CRP participates in some FFLs, such as araBAD and araFGH operons and adapts to fluctuating environmental nutrients, thereby enhancing the survivability of E. coli. Although computational simulations have been conducted to explore the potential functionality of FFLs, a comprehensive study on the functions of all structural types on the basis of in vivo data is lacking. Moreover, the regulatory role of CRP-mediated FFLs (CRP-FFLs) remains obscure. We identified 393 CRP-FFLs in E. coli using EcoCyc and RegulonDB. Dose⁻response genomic microarray of E. coli revealed dynamic gene expression of each target gene of CRP-FFLs in response to a range of cAMP dosages. All eight types of FFLs were present in CRP regulon with various expression patterns of each CRP-FFL, which were further divided into five functional groups. The microarray and reported regulatory relationships identified 202 CRP-FFLs that were directly regulated by CRP in these eight types of FFLs. Interestingly, 34% (147/432) of genes were directly regulated by CRP and CRP-regulated TFs, which indicates that these CRP-regulated genes were also regulated by other CRP-regulated TFs responding to environmental signals through CRP-FFLs. Furthermore, we applied gene ontology annotation to reveal the biological functions of CRP-FFLs.