Non-linear aspects of dynamic pattern in cellular networks

Effects of noise on the critical points of Turing instability in complex ecosystems

Noise is ubiquitous in natural and artificial systems. In a noisy environment, the interactions among nodes may fluctuate randomly, leading to more complicated interactions. In this paper we focus on the effects of noise and network topology on the Turing pattern of ecological networks with activator-inhibitor structure, which may be interpreted as prey-predator interactions. Based on the stability theory of stochastic differential equations, a sufficient condition for the uniform state is derived. The analytical results indicate that noise is beneficial for the uniform state. When the ratio between the diffusion coefficients of the predator and prey increases, the ecosystems can exhibit a transition from a uniform stable state to a Turing pattern, while when the ratio decreases, the ecosystems transit from a Turing pattern to a uniform stable state. There are two crucial critical points in Turing patterns, forward and backward. We find that both forward and backward critical points increase as the noise intensity increases. This means that noise favors a stable homogeneous state compared to a state with a heterogeneous pattern, which is consistent with the analytical results. In addition, noise can weaken the hysteresis phenomenon and even eliminate it in some cases. Furthermore, we report that network topology plays an important role in modulating the uniform state of ecosystems, such as the size of prey-predator systems, the network connectivity, and the strength of interaction.

Turing patterns in simplicial complexes

The spontaneous emergence of patterns in nature, such as stripes and spots, can be mathematically explained by reaction-diffusion systems. These patterns are often referred as Turing patterns to honor the seminal work of Alan Turing in the early 1950s. With the coming of age of network science, and with its related departure from diffusive nearest-neighbor interactions to long-range links between nodes, additional layers of complexity behind pattern formation have been discovered, including irregular spatiotemporal patterns. Here we investigate the formation of Turing patterns in simplicial complexes, where links no longer connect just pairs of nodes but can connect three or more nodes. Such higher-order interactions are emerging as a new frontier in network science, in particular describing group interaction in various sociological and biological systems, so understanding pattern formation under these conditions is of the utmost importance. We show that a canonical reaction-diffusion system defined over a simplicial complex yields Turing patterns that fundamentally differ from patterns observed in traditional networks. For example, we observe a stable distribution of Turing patterns where the fraction of nodes with reactant concentrations above the equilibrium point is exponentially related to the average degree of 2-simplexes, and we uncover parameter regions where Turing patterns will emerge only under higher-order interactions, but not under pairwise interactions.

Reaction-diffusion models in weighted and directed connectomes

Connectomes represent comprehensive descriptions of neural connections in a nervous system to better understand and model central brain function and peripheral processing of afferent and efferent neural signals. Connectomes can be considered as a distinctive and necessary structural component alongside glial, vascular, neurochemical, and metabolic networks of the nervous systems of higher organisms that are required for the control of body functions and interaction with the environment. They are carriers of functional epiphenomena such as planning behavior and cognition, which are based on the processing of highly dynamic neural signaling patterns. In this study, we examine more detailed connectomes with edge weighting and orientation properties, in which reciprocal neuronal connections are also considered. Diffusion processes are a further necessary condition for generating dynamic bioelectric patterns in connectomes. Based on our high-precision connectome data, we investigate different diffusion-reaction models to study the propagation of dynamic concentration patterns in control and lesioned connectomes. Therefore, differential equations for modeling diffusion were combined with well-known reaction terms to allow the use of connection weights, connectivity orientation and spatial distances.

Three reaction-diffusion systems Gray-Scott, Gierer-Meinhardt and Mimura-Murray were investigated. For this purpose, implicit solvers were implemented in a numerically stable reaction-diffusion system within the framework of neuroVIISAS. The implemented reaction-diffusion systems were applied to a subconnectome which shapes the mechanosensitive pathway that is strongly affected in the multiple sclerosis demyelination disease. It was found that demyelination modeling by connectivity weight modulation changes the oscillations of the target region, i.e. the primary somatosensory cortex, of the mechanosensitive pathway.

In conclusion, a new application of reaction-diffusion systems to weighted and directed connectomes has been realized. Because the implementation were performed in the neuroVIISAS framework many possibilities for the study of dynamic reaction-diffusion processes in empirical connectomes as well as specific randomized network models are available now.

Turing instability in quantum activator–inhibitor systems

Turing instability is a fundamental mechanism of nonequilibrium self-organization. However, despite the universality of its essential mechanism, Turing instability has thus far been investigated mostly in classical systems. In this study, we show that Turing instability can occur in a quantum dissipative system and analyze its quantum features such as entanglement and the effect of measurement. We propose a degenerate parametric oscillator with nonlinear damping in quantum optics as a quantum activator–inhibitor unit and demonstrate that a system of two such units can undergo Turing instability when diffusively coupled with each other. The Turing instability induces nonuniformity and entanglement between the two units and gives rise to a pair of nonuniform states that are mixed due to quantum noise. Further performing continuous measurement on the coupled system reveals the nonuniformity caused by the Turing instability. Our results extend the universality of the Turing mechanism to the quantum realm and may provide a novel perspective on the possibility of quantum nonequilibrium self-organization and its application in quantum technologies.

The qualitative and quantitative relationships between pattern formation and average degree in networked reaction-diffusion systems

The Turing pattern is an important dynamic behavior characteristic of activator-inhibitor systems. Differentiating from traditional assumption of activator-inhibitor interactions in a spatially continuous domain, a Turing pattern in networked reaction-diffusion systems has received much attention during the past few decades. In spite of its great progress, it still fails to evaluate the precise influences of network topology on pattern formation. To this end, we try to promote the research on this important and interesting issue from the point of view of average degree-a critical topological feature of networks. We first qualitatively analyze the influence of average degree on pattern formation. Then, a quantitative relationship between pattern formation and average degree, the exponential decay of pattern formation, is proposed via nonlinear regression. The finding holds true for several activator-inhibitor systems including biology model, ecology model, and chemistry model. The significance of this study lies that the exponential decay not only quantitatively depicts the influence of average degree on pattern formation, but also provides the possibility for predicting and controlling pattern formation.

Turing pattern induced by the directed ER network and delay

Infectious diseases generally spread along with the asymmetry of social network propagation because the asymmetry of urban development and the prevention strategies often affect the direction of the movement. But the spreading mechanism of the epidemic remains to explore in the directed network. In this paper, the main effect of the directed network and delay on the dynamic behaviors of the epidemic is investigated. The algebraic expressions of Turing instability are given to show the role of the directed network in the spread of the epidemic, which overcomes the drawback that undirected networks cannot lead to the outbreaks of infectious diseases. Then, Hopf bifurcation is analyzed to illustrate the dynamic mechanism of the periodic outbreak, which is consistent with the transmission of COVID-19. Also, the discrepancy ratio between the imported and the exported is proposed to explain the importance of quarantine policies and the spread mechanism. Finally, the theoretical results are verified by numerical simulation.

Optimal control of networked reaction-diffusion systems

Patterns in nature are fascinating both aesthetically and scientifically. Alan Turing's celebrated reaction-diffusion model of pattern formation from the 1950s has been extended to an astounding diversity of applications: from cancer medicine, via nanoparticle fabrication, to computer architecture. Recently, several authors have studied pattern formation in underlying networks, but thus far, controlling a reaction-diffusion system in a network to obtain a particular pattern has remained elusive. We present a solution to this problem in the form of an analytical framework and numerical algorithm for optimal control of Turing patterns in networks. We demonstrate our method's effectiveness and discuss factors that affect its performance. We also pave the way for multidisciplinary applications of our framework beyond reaction-diffusion models.

Body size dependent dispersal influences stability in heterogeneous metacommunities

Body size affects key biological processes across the tree of life, with particular importance for food web dynamics and stability. Traits influencing movement capabilities depend strongly on body size, yet the effects of allometrically-structured dispersal on food web stability are less well understood than other demographic processes. Here we study the stability properties of spatially-arranged model food webs in which larger bodied species occupy higher trophic positions, while species’ body sizes also determine the rates at which they traverse spatial networks of heterogeneous habitat patches. Our analysis shows an apparent stabilizing effect of positive dispersal rate scaling with body size compared to negative scaling relationships or uniform dispersal. However, as the global coupling strength among patches increases, the benefits of positive body size-dispersal scaling disappear. A permutational analysis shows that breaking allometric dispersal hierarchies while preserving dispersal rate distributions rarely alters qualitative aspects of metacommunity stability. Taken together, these results suggest that the oft-predicted stabilizing effects of large mobile predators may, for some dimensions of ecological stability, be attributed to increased patch coupling per se, and not necessarily coupling by top trophic levels in particular.

Cooperative metabolic resource allocation in spatially-structured systems

Natural selection has shaped the evolution of cells and multi-cellular organisms such that social cooperation can often be preferred over an individualistic approach to metabolic regulation. This paper extends a framework for dynamic metabolic resource allocation based on the maximum entropy principle to spatiotemporal models of metabolism with cooperation. Much like the maximum entropy principle encapsulates 'bet-hedging' behaviour displayed by organisms dealing with future uncertainty in a fluctuating environment, its cooperative extension describes how individuals adapt their metabolic resource allocation strategy to further accommodate limited knowledge about the welfare of others within a community. The resulting theory explains why local regulation of metabolic cross-feeding can fulfil a community-wide metabolic objective if individuals take into consideration an ensemble measure of total population performance as the only form of global information. The latter is likely supplied by quorum sensing in microbial systems or signalling molecules such as hormones in multi-cellular eukaryotic organisms.

Resilience for stochastic systems interacting via a quasi-degenerate network

A stochastic reaction-diffusion model is studied on a networked support. In each patch of the network, two species are assumed to interact following a non-normal reaction scheme. When the interaction unit is replicated on a directed linear lattice, noise gets amplified via a self-consistent process, which we trace back to the degenerate spectrum of the embedding support. The same phenomenon holds when the system is bound to explore a quasidegenerate network. In this case, the eigenvalues of the Laplacian operator, which governs species diffusion, accumulate over a limited portion of the complex plane. The larger the network, the more pronounced the amplification. Beyond a critical network size, a system deemed deterministically stable, hence resilient, can develop seemingly regular patterns in the concentration amount. Non-normality and quasidegenerate networks may, therefore, amplify the inherent stochasticity and so contribute to altering the perception of resilience, as quantified via conventional deterministic methods.

Pattern Formation over Multigraphs

Two of the most common pattern formation mechanisms are Turing-patterning in reaction-diffusion systems and lateral inhibition of neighboring cells. In this paper, we introduce a broad dynamical model of interconnected modules to study the emergence of patterns, with the above mentioned two mechanisms as special cases. Our results do not restrict the number of modules or their complexity, allow multiple layers of communication channels with possibly different interconnection structure, and do not assume symmetric connections between two connected modules. Leveraging only the static input/output properties of the subsystems and the spectral properties of the interconnection matrices, we characterize the stability of the homogeneous fixed points as well as sufficient conditions for the emergence of spatially non-homogeneous patterns. To obtain these results, we rely on properties of the graphs together with tools from monotone systems theory. As application examples, we consider patterning in neural networks, in reaction-diffusion systems, and contagion processes over random graphs.

Theory of Turing Patterns on Time Varying Networks

The process of pattern formation for a multispecies model anchored on a time varying network is studied. A nonhomogeneous perturbation superposed to an homogeneous stable fixed point can be amplified following the Turing mechanism of instability, solely instigated by the network dynamics. By properly tuning the frequency of the imposed network evolution, one can make the examined system behave as its averaged counterpart, over a finite time window. This is the key observation to derive a closed analytical prediction for the onset of the instability in the time dependent framework. Continuously and piecewise constant periodic time varying networks are analyzed, setting the framework for the proposed approach. The extension to nonperiodic settings is also discussed.

Graph spectral characterization of the XY model on complex networks

There is recent evidence that the XY spin model on complex networks can display three different macroscopic states in response to the topology of the network underpinning the interactions of the spins. In this work we present a way to characterize the macroscopic states of the XY spin model based on the spectral decomposition of time series using topological information about the underlying networks. We use three different classes of networks to generate time series of the spins for the three possible macroscopic states. We then use the temporal Graph Signal Transform technique to decompose the time series of the spins on the eigenbasis of the Laplacian. From this decomposition, we produce spatial power spectra, which summarize the activation of structural modes by the nonlinear dynamics, and thus coherent patterns of activity of the spins. These signatures of the macroscopic states are independent of the underlying network class and can thus be used as robust signatures for the macroscopic states. This work opens avenues to analyze and characterize dynamics on complex networks using temporal Graph Signal Analysis.

Multiple-scale theory of topology-driven patterns on directed networks

Dynamical processes on networks are currently being considered in different domains of cross-disciplinary interest. Reaction-diffusion systems hosted on directed graphs are in particular relevant for their widespread applications, from computer networks to traffic systems. Due to the peculiar spectrum of the discrete Laplacian operator, homogeneous fixed points can turn unstable, on a directed support, because of the topology of the network, a phenomenon which cannot be induced on undirected graphs. A linear analysis can be performed to single out the conditions that underly the instability. The complete characterization of the patterns, which are eventually attained beyond the linear regime of exponential growth, calls instead for a full nonlinear treatment. By performing a multiple time scale perturbative calculation, we here derive an effective equation for the nonlinear evolution of the amplitude of the most unstable mode, close to the threshold of criticality. This is a Stuart-Landau equation the complex coefficients of which appear to depend on the topological features of the embedding directed graph. The theory proves adequate versus simulations, as confirmed by operating with a paradigmatic reaction-diffusion model.

Turing-like instabilities from a limit cycle

The Turing instability is a paradigmatic route to pattern formation in reaction-diffusion systems. Following a diffusion-driven instability, homogeneous fixed points can become unstable when subject to external perturbation. As a consequence, the system evolves towards a stationary, nonhomogeneous attractor. Stable patterns can be also obtained via oscillation quenching of an initially synchronous state of diffusively coupled oscillators. In the literature this is known as the oscillation death phenomenon. Here, we show that oscillation death is nothing but a Turing instability for the first return map of the system in its synchronous periodic state. In particular, we obtain a set of approximated closed conditions for identifying the domain in the parameter space that yields the instability. This is a natural generalization of the original Turing relations, to the case where the homogeneous solution of the examined system is a periodic function of time. The obtained framework applies to systems embedded in continuum space, as well as those defined on a networklike support. The predictive ability of the theory is tested numerically, using different reaction schemes.

Pattern formation in multiplex networks

The advances in understanding complex networks have generated increasing interest in dynamical processes occurring on them. Pattern formation in activator-inhibitor systems has been studied in networks, revealing differences from the classical continuous media. Here we study pattern formation in a new framework, namely multiplex networks. These are systems where activator and inhibitor species occupy separate nodes in different layers. Species react across layers but diffuse only within their own layer of distinct network topology. This multiplicity generates heterogeneous patterns with significant differences from those observed in single-layer networks. Remarkably, diffusion-induced instability can occur even if the two species have the same mobility rates; condition which can never destabilize single-layer networks. The instability condition is revealed using perturbation theory and expressed by a combination of degrees in the different layers. Our theory demonstrates that the existence of such topology-driven instabilities is generic in multiplex networks, providing a new mechanism of pattern formation.

Turing patterns in multiplex networks

The theory of patterns formation for a reaction-diffusion system defined on a multiplex is developed by means of a perturbative approach. The interlayer diffusion constants act as a small parameter in the expansion and the unperturbed state coincides with the limiting setting where the multiplex layers are decoupled. The interaction between adjacent layers can seed the instability of a homogeneous fixed point, yielding self-organized patterns which are instead impeded in the limit of decoupled layers. Patterns on individual layers can also fade away due to cross-talking between layers. Analytical results are compared to direct simulations.

Dispersal-induced destabilization of metapopulations and oscillatory Turing patterns in ecological networks

As shown by Alan Turing in 1952, differential diffusion may destabilize uniform distributions of reacting species and lead to emergence of patterns. While stationary Turing patterns are broadly known, the oscillatory instability, leading to traveling waves in continuous media and sometimes called the wave bifurcation, remains less investigated. Here, we extend the original analysis by Turing to networks and apply it to ecological metapopulations with dispersal connections between habitats. Remarkably, the oscillatory Turing instability does not lead to wave patterns in networks, but to spontaneous development of heterogeneous oscillations and possible extinction of species. We find such oscillatory instabilities for all possible food webs with three predator or prey species, under various assumptions about the mobility of individual species and nonlinear interactions between them. Hence, the oscillatory Turing instability should be generic and must play a fundamental role in metapopulation dynamics, providing a common mechanism for dispersal-induced destabilization of ecosystems.

Chemical computing based on Turing patterns in two coupled cells with equal transport coefficients

Two diffusively coupled reaction cells with a nonlinear reaction are used to perform chemical computing based on targeted perturbations switching between two Turing patterns defining two states of a logical device.

Control of Turing patterns and their usage as sensors, memory arrays, and logic gates

We study a model system of three diffusively coupled reaction cells arranged in a linear array that display Turing patterns with special focus on the case of equal coupling strength for all components. As a suitable model reaction we consider a two-variable core model of glycolysis. Using numerical continuation and bifurcation techniques we analyze the dependence of the system's steady states on varying rate coefficient of the recycling step while the coupling coefficients of the inhibitor and activator are fixed and set at the ratios 100:1, 1:1, and 4:5. We show that stable Turing patterns occur at all three ratios but, as expected, spontaneous transition from the spatially uniform steady state to the spatially nonuniform Turing patterns occurs only in the first case. The other two cases possess multiple Turing patterns, which are stabilized by secondary bifurcations and coexist with stable uniform periodic oscillations. For the 1:1 ratio we examine modular spatiotemporal perturbations, which allow for controllable switching between the uniform oscillations and various Turing patterns. Such modular perturbations are then used to construct chemical computing devices utilizing the multiple Turing patterns. By classifying various responses we propose: (a) a single-input resettable sensor capable of reading certain value of concentration, (b) two-input and three-input memory arrays capable of storing logic information, (c) three-input, three-output logic gates performing combinations of logical functions OR, XOR, AND, and NAND.

Stochastic Turing patterns on a network

The process of stochastic Turing instability on a scale-free network is discussed for a specific case study: the stochastic Brusselator model. The system is shown to spontaneously differentiate into activator-rich and activator-poor nodes outside the region of parameters classically deputed to the deterministic Turing instability. This phenomenon, as revealed by direct stochastic simulations, is explained analytically and eventually traced back to the finite-size corrections stemming from the inherent graininess of the scrutinized medium.

Dissipative structures in a two-cell system: Numerical and experimental approaches

It has been shown that the coupling between the photoreduction of the oxidized form of dichloroindophenol (an artificial electron acceptor) by thylakoids and the incident light intensity can lead to the appearance of multiple steady states when the system is operated under open conditions. In the present work, a numerical study and experimental evidence are presented on the occurrence of dissipative structures in an arrangement of two continuously stirred tank reactors with mutual mass exchange of dichloroindophenol through an inert membrane. The stable spatial structures are generated by the creation of transient internal and external asymmetries. A nontrivial hysteresis effect between symmetric and asymmetric stable steady states has been observed.

A hypothetical mathematical construct explaining the mechanism of biological amplification in an experimental model utilizing picoTesla (PT) electromagnetic fields

We seek to answer the conundrum: What is the fundamental mechanism by which very weak, low frequency Electromagnetic fields influence biosystems? In considering the hydrophobicity of intramembranous protein (IMP) H-bonds which cross the phospholipid bilayer of plasma membranes, and the necessity for photonic recycling in cell surface interactions after dissipation of energetic states, present models lack structure and thermodynamic properties to maintain (DeltaE) sufficient energy sources necessary for amplifications by factors of 10(12). Even though one accepts that the ligand-receptor association alters the conformation of extracellular, extruding portions of IMP's at the cell surface, and that this change can be transmitted to the cytoplasm by the transmembranous helical segments by nonlinear vibrations of proteins with generation of soliton waves, one is still unable to account for repair and balanced function. Indeed, responses of critical molecules to certain magnetic field signals may include enhanced vibrational amplitudes, increased quanta of thermal energies and order inducing interactions. We may accept that microtrabecular reticulum-receptor is associated with actin filaments and ATP molecules which contribute to the activation of the cyclase enzyme system through piezoelectricity. Magnetic fields will pass through the membrane which sharply attenuates the electric field component of an EM field, due to its high impedance. Furthermore, EM oscillations are converted to mechanical vibrations; i.e., photon-phonon transduction, to induce molecular vibrations of frequencies specifically responsible for bioamplifications of weak triggers at the membrane surface, as well as GAP junctions. The hydrogen bonds of considerable importance are those in proteins (10(12)Hz) and DNA (10(11)Hz) and may be viewed as centers of EM radiation emission in the range from the mm microwaves to the far IR. However, classical electrodynamical theory does not yield a model for biomolecular resonant responses which are integrated over time and account for the connection between the phonon field and photons. Jacobson Resonance does supply an initial physical mechanism, as equivalencies in energy to that of Zeeman Resonance (i.e., zero-order magnetic resonance) and cyclotron resonance may be derived from the DeBroglie wave particle equation. For the first time, we view the introduction of Relativity Theory to biology in the expression, mc(2)=BvLq, where m is the mass of a particle in the 'box' or 'string' (molecule in a biosystem), c is the velocity of electromagnetic field in space, independent of its inertial frame of reference, B is the magnetic flux density,v is the velocity of the carrier or 'string' (a one or two dimensional 'box') in which the particle exists, L is its dimension (length) and q represents a unit charge q=1C, by defining electromotive force as energy per unit charge. Equivalencies suggest that qvBL is one of the fundamental expressions of energy of a charged wave-particle in magnetic fields, just as Zeeman and cyclotron resonance energy expressions, gbetaB and qhB/2pim, and is applicable to all charged particles (molecules in biological systems). There may exist spontaneous, independent and incessant interactions of magnetic vector B and particles in biosystems which exert Lorentz forces. Lorentz forces may be transmitted from EM field to gravitational field as a gravity wave which return to the phonon field as microgravitational fluctuations to therein produce quantum vibrational states that increase quanta of thermal energies integrated over time. This may account for the differential of 10(12) between photonic energy of ELF waves and the Boltzman energy kT. Recent data from in vivo controlled studies are included as empirical support for the various hypotheses presented.

The effects of cell density and metabolite flux on cellular dynamics

Density-dependent regulation of cell growth in tissue culture is a well-known phenomenon but the mechanism of regulation remains obscure. Here we explore the effects of cell density and metabolite flux on the collective dynamics of a cell population. The intracellular dynamics are modelled by positive feedback kinetic mechanisms of the kind known to apply to yeast cells. Several experimental observations related to glycolytic oscillations are predicted and it is suggested that the general conclusions may be applicable in a broader context.

Control of somite number during morphogenesis of a vertebrate, Xenopus laevis

During Xenopus embryogenesis and early growth, somite number stays close to a species-typical value for each morphological stage. This remains true even after operations on blastulae which lead to the development of abnormally small but otherwise complete early embryos, involving reduction in number of cells assigned to each somite. Evidence presented suggests that a body-position gradient may be involved, but in rather different ways at different stages, in controlling total somite number.