Heterogeneous network flow and Petri nets characterize multilayer complex networks.

Interacting subsystems are commonly described by networks, where multimodal behaviour found in most natural or engineered systems found recent extension in form of multilayer networks. Since multimodal interaction is often not dictated by network topology alone and may manifest in form of cross-layer information exchange, multilayer network flow becomes of relevant further interest. Rationale can be found in most interacting subsystems, where a form of multimodal flow across layers can be observed in e.g., chemical processes, energy networks, logistics, finance, or any other form of conversion process relying on the laws of conservation. To this end, the formal notion of heterogeneous network flow is proposed, as a multilayer flow function aligned with the theory of network flow. Furthermore, dynamic equivalence is established with the framework of Petri nets, as the baseline model of concurrent event systems. Application of the resulting multilayer Laplacian flow and flow centrality is presented, along with graph learning based inference of multilayer relationships over multimodal data. On synthetic data the proposed framework demonstrates benefits of multimodal flow derivation in critical component identification. It also displays applicability in relationship inference (learning based function approximation) on multimodal time series. On real-world data the proposed framework provides, among others, multimodal flow interpretation of U.S. economic activity, uncovering underlying empirical steady state probability distribution, as well as inherent network (economic) robustness.

On convergence for hybrid models of gene regulatory networks under polytopic uncertainties: a Lyapunov approach.

Hybrid models of genetic regulatory networks allow for a simpler analysis with respect to fully detailed quantitative models, still maintaining the main dynamical features of interest. In this paper we consider a piecewise affine model of a genetic regulatory network, in which the parameters describing the production function are affected by polytopic uncertainties. In the first part of the paper, after recalling how the problem of finding a Lyapunov function is solved in the nominal case, we present the considered polytopic uncertain system and then, after describing how to deal with sliding mode solutions, we prove a result of existence of a parameter dependent Lyapunov function subject to the solution of a feasibility linear matrix inequalities problem. In the second part of the paper, based on the previously described Lyapunov function, we are able to determine a set of domains where the system is guaranteed to converge, with the exception of a zero measure set of times, independently from the uncertainty realization. Finally a three nodes network example shows the validity of the results.

Scheduling of energy storage.

The increasing reliance on renewable energy generation means that storage may well play a much greater role in the balancing of future electricity systems. We show how heterogeneous stores, differing in capacity and rate constraints, may be optimally, or nearly optimally, scheduled to assist in such balancing, with the aim of minimizing the total imbalance (unserved energy) over any given period of time. It further turns out that in many cases the optimal policies are such that the optimal decision at each point in time is independent of the future evolution of the supply-demand balance in the system, so that these policies remain optimal in a stochastic environment. This article is part of the theme issue 'The mathematics of energy systems'.

A computational framework for a Lyapunov-enabled analysis of biochemical reaction networks.

Complex molecular biological processes such as transcription and translation, signal transduction, post-translational modification cascades, and metabolic pathways can be described in principle by biochemical reactions that explicitly take into account the sophisticated network of chemical interactions regulating cell life. The ability to deduce the possible qualitative behaviors of such networks from a set of reactions is a central objective and an ongoing challenge in the field of systems biology. Unfortunately, the construction of complete mathematical models is often hindered by a pervasive problem: despite the wealth of qualitative graphical knowledge about network interactions, the form of the governing nonlinearities and/or the values of kinetic constants are hard to uncover experimentally. The kinetics can also change with environmental variations. This work addresses the following question: given a set of reactions and without assuming a particular form for the kinetics, what can we say about the asymptotic behavior of the network? Specifically, it introduces a class of networks that are "structurally (mono) attractive" meaning that they are incapable of exhibiting multiple steady states, oscillation, or chaos by virtue of their reaction graphs. These networks are characterized by the existence of a universal energy-like function called a Robust Lyapunov function (RLF). To find such functions, a finite set of rank-one linear systems is introduced, which form the extremals of a linear convex cone. The problem is then reduced to that of finding a common Lyapunov function for this set of extremals. Based on this characterization, a computational package, Lyapunov-Enabled Analysis of Reaction Networks (LEARN), is provided that constructs such functions or rules out their existence. An extensive study of biochemical networks demonstrates that LEARN offers a new unified framework. Basic motifs, three-body binding, and genetic networks are studied first. The work then focuses on cellular signalling networks including various post-translational modification cascades, phosphotransfer and phosphorelay networks, T-cell kinetic proofreading, and ERK signalling. The Ribosome Flow Model is also studied.

Graph-theoretic characterizations of monotonicity of chemical networks in reaction coordinates.

This paper derives new results for certain classes of chemical reaction networks, linking structural to dynamical properties. In particular, it investigates their monotonicity and convergence under the assumption that the rates of the reactions are monotone functions of the concentrations of their reactants. This is satisfied for, yet not restricted to, the most common choices of the reaction kinetics such as mass action, Michaelis-Menten and Hill kinetics. The key idea is to find an alternative representation under which the resulting system is monotone. As a simple example, the paper shows that a phosphorylation/dephosphorylation process, which is involved in many signaling cascades, has a global stability property. We also provide a global stability result for a more complicated example that describes a regulatory pathway of a prevalent signal transduction module, the MAPK cascade.

Chemical networks with inflows and outflows: a positive linear differential inclusions approach.

Certain mass-action kinetics models of biochemical reaction networks, although described by nonlinear differential equations, may be partially viewed as state-dependent linear time-varying systems, which in turn may be modeled by convex compact valued positive linear differential inclusions. A result is provided on asymptotic stability of such inclusions, and applied to a ubiquitous biochemical reaction network with inflows and outflows, known as the futile cycle. We also provide a characterization of exponential stability of general homogeneous switched systems which is not only of interest in itself, but also plays a role in the analysis of the futile cycle.

Convergence in networks with counterclockwise neural dynamics.

The notion of counterclockwise (ccw) input-output (I-O) dynamics, introduced by Angeli (2006) to deal with questions of multistability in interconnected dynamical systems, is applied and further developed in order to analyze convergence and stability of neural networks. By pursuing a modular approach, we interpret a cellular nonlinear network (CNN) as a positive feedback of a parallel block of single-input-single-output (SISO) dynamical systems, the neurons, and a static multiple-input-multiple-output (MIMO) system that couples them (typically the so-called interconnection matrix). The analysis extends previously known results by enlarging the class of allowed neural dynamics to higher order neurons.

A Petri net approach to the study of persistence in chemical reaction networks.

Persistence is the property, for differential equations in R(n), that solutions starting in the positive orthant do not approach the boundary of the orthant. For chemical reactions and population models, this translates into the non-extinction property: provided that every species is present at the start of the reaction, no species will tend to be eliminated in the course of the reaction. This paper provides checkable conditions for persistence of chemical species in reaction networks, using concepts and tools from Petri net theory, and verifies these conditions on various systems which arise in the modeling of cell signaling pathways.

On predator-prey systems and small-gain theorems.

This paper deals with an almost global convergence result for Lotka-Volterra systems with predator-prey interactions. These systems can be written as (negative) feedback systems. The subsystems of the feedback loop are monotone control systems, possessing particular input-output properties. We use a small-gain theorem, adapted to a context of systems with multiple equilibrium points to obtain the desired almost global convergence result, which provides sufficient conditions to rule out oscillatory or more complicated behavior that is often observed in predator-prey systems.

Detection of multistability, bifurcations, and hysteresis in a large class of biological positive-feedback systems.

It is becoming increasingly clear that bistability (or, more generally, multistability) is an important recurring theme in cell signaling. Bistability may be of particular relevance to biological systems that switch between discrete states, generate oscillatory responses, or "remember" transitory stimuli. Standard mathematical methods allow the detection of bistability in some very simple feedback systems (systems with one or two proteins or genes that either activate each other or inhibit each other), but realistic depictions of signal transduction networks are invariably much more complex. Here, we show that for a class of feedback systems of arbitrary order the stability properties of the system can be deduced mathematically from how the system behaves when feedback is blocked. Provided that this open-loop, feedback-blocked system is monotone and possesses a sigmoidal characteristic, the system is guaranteed to be bistable for some range of feedback strengths. We present a simple graphical method for deducing the stability behavior and bifurcation diagrams for such systems and illustrate the method with two examples taken from recent experimental studies of bistable systems: a two-variable Cdc2/Wee1 system and a more complicated five-variable mitogen-activated protein kinase cascade.