Fluctuation-driven Turing patterns

Computation of random time-shift distributions for stochastic population models

Even in large systems, the effect of noise arising from when populations are initially small can persist to be measurable on the macroscale. A deterministic approximation to a stochastic model will fail to capture this effect, but it can be accurately approximated by including an additional random time-shift to the initial conditions. We present a efficient numerical method to compute this time-shift distribution for a large class of stochastic models. The method relies on differentiation of certain functional equations, which we show can be effectively automated by deriving rules for different types of model rates that arise commonly when mass-action mixing is assumed. Explicit computation of the time-shift distribution can be used to build a practical tool for the efficient generation of macroscopic trajectories of stochastic population models, without the need for costly stochastic simulations. Full code is provided to implement the calculations and we demonstrate the method on an epidemic model and a model of within-host viral dynamics.

Computing macroscopic reaction rates in reaction-diffusion systems using Monte Carlo simulations

Stochastic reaction-diffusion models are employed to represent many complex physical, biological, societal, and ecological systems. The macroscopic reaction rates describing the large-scale, long-time kinetics in such systems are effective, scale-dependent renormalized parameters that need to be either measured experimentally or computed by means of a microscopic model. In a Monte Carlo simulation of stochastic reaction-diffusion systems, microscopic probabilities for specific events to happen serve as the input control parameters. To match the results of any computer simulation to observations or experiments carried out on the macroscale, a mapping is required between the microscopic probabilities that define the Monte Carlo algorithm and the macroscopic reaction rates that are experimentally measured. Finding the functional dependence of emergent macroscopic rates on the microscopic probabilities (subject to specific rules of interaction) is a very difficult problem, and there is currently no systematic, accurate analytical way to achieve this goal. Therefore, we introduce a straightforward numerical method of using lattice Monte Carlo simulations to evaluate the macroscopic reaction rates by directly obtaining the count statistics of how many events occur per simulation time step. Our technique is first tested on well-understood fundamental examples, namely, restricted birth processes, diffusion-limited two-particle coagulation, and two-species pair annihilation kinetics. Next we utilize the thus gained experience to investigate how the microscopic algorithmic probabilities become coarse-grained into effective macroscopic rates in more complex model systems such as the Lotka-Volterra model for predator-prey competition and coexistence, as well as the rock-paper-scissors or cyclic Lotka-Volterra model and its May-Leonard variant that capture population dynamics with cyclic dominance motifs. Thereby we achieve a more thorough and deeper understanding of coarse graining in spatially extended stochastic reaction-diffusion systems and the nontrivial relationships between the associated microscopic and macroscopic model parameters, with a focus on ecological systems. The proposed technique should generally provide a useful means to better fit Monte Carlo simulation results to experimental or observational data.

Patterns governed by chemotaxis and time delay

In this paper, sufficient conditions of Turing instability are established for general delayed reaction-diffusion-chemotaxis models with no-flux boundary conditions. In particular, we address the difficulty brought about by the time delay in investigating the Turing instability. These models involve the time delay parameter τ and the general density (concentration) function in chemotaxis terms with the chemotaxis parameter χ. Theoretical results reveal that the time delay parameter τ could determine the stability of positive equilibrium for ordinary differential equations, while the chemotaxis parameter χ could describe the stability of positive equilibrium for partial differential equations. In this fashion, the general conditions of Turing instability are presented. To confirm their validity, the delayed chemotaxis-type predator-prey model and the phytoplankton-zooplankton model are considered. It is found that these two models admit Turing instability, and the numerical results are in good agreement with the theoretical analysis. The obtained results are helpful in the application of the Turing pattern in models with time delay and chemotaxis.

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.

Symmetry-simplicity, broken symmetry-complexity

Complex phenomena are made possible when: (i) fundamental physical symmetries are broken and (ii) from the set of broken symmetries historically selected ground states are applied to performing mechanical work and storing adaptive information. Over the course of several decades Philip Anderson enumerated several key principles that can follow from broken symmetry in complex systems. These include emergence, frustrated random functions, autonomy and generalized rigidity. I describe these as the four Anderson Principles all of which are preconditions for the emergence of evolved function. I summarize these ideas and discuss briefly recent extensions that engage with the related concept of functional symmetry breaking, inclusive of information, computation and causality.

Numerical and renormalization group analysis of the phase diagram of a stochastic cubic autocatalytic reaction-diffusion system

The renormalization group is a set of tools that can be used to incorporate the effect of fluctuations in a dynamical system as a rescaling of the system's parameters. Here, we apply the renormalization group to a pattern-forming stochastic cubic autocatalytic reaction-diffusion model and compare its predictions with numerical simulations. Our results demonstrate a good agreement within the range of validity of the theory and show that external noise can be used as a control parameter in such systems.

Constrained Langevin approximation for the Togashi-Kaneko model of autocatalytic reactions

The Togashi Kaneko model (TK model) is a simple stochastic reaction network that displays discreteness-induced transitions between meta-stable patterns. Here we study a constrained Langevin approximation (CLA) of this model. This CLA, derived under the classical scaling, is an obliquely reflected diffusion process on the positive orthant and hence respects the constraint that chemical concentrations are never negative. We show that the CLA is a Feller process, is positive Harris recurrent and converges exponentially fast to the unique stationary distribution. We also characterize the stationary distribution and show that it has finite moments. In addition, we simulate both the TK model and its CLA in various dimensions. For example, we describe how the TK model switches between meta-stable patterns in dimension six. Our simulations suggest that, when the volume of the vessel in which all of the reactions that take place is large, the CLA is a good approximation of the TK model in terms of both the stationary distribution and the transition times between patterns.

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.

Stochastic Model for Quasi-One-Dimensional Transitional Turbulence with Streamwise Shear Interactions

The transition to turbulence in wall-bounded shear flows is typically subcritical, with a poorly understood interplay between spatial fluctuations, pattern formation, and stochasticity near the critical Reynolds number. Here, we present a spatially extended stochastic minimal model for the energy budget in transitional pipe flow, which successfully recapitulates the way localized patches of turbulence (puffs) decay, split, and grow, respectively, as the Reynolds number increases through the laminar-turbulent transition. Our approach takes into account the flow geometry, as we demonstrate by extending the model to quasi-one-dimensional Taylor-Couette flow, reproducing the observed directed percolation pattern of turbulent patches in space and time.

Synthetic spatial patterning in bacteria: advances based on novel diffusible signals

Engineering multicellular patterning may help in the understanding of some fundamental laws of pattern formation and thus may contribute to the field of developmental biology. Furthermore, advanced spatial control over gene expression may revolutionize fields such as medicine, through organoid or tissue engineering. To date, foundational advances in spatial synthetic biology have often been made in prokaryotes, using artificial gene circuits. In this review, engineered patterns are classified into four levels of increasing complexity, ranging from spatial systems with no diffusible signals to systems with complex multi-diffusor interactions. This classification highlights how the field was held back by a lack of diffusible components. Consequently, we provide a summary of both previously characterized and some new potential candidate small-molecule signals that can regulate gene expression in Escherichia coli. These diffusive signals will help synthetic biologists to successfully engineer increasingly intricate, robust and tuneable spatial structures.

Turing pattern design principles and their robustness

Turing patterns have morphed from mathematical curiosities into highly desirable targets for synthetic biology. For a long time, their biological significance was sometimes disputed but there is now ample evidence for their involvement in processes ranging from skin pigmentation to digit and limb formation. While their role in developmental biology is now firmly established, their synthetic design has so far proved challenging. Here, we review recent large-scale mathematical analyses that have attempted to narrow down potential design principles. We consider different aspects of robustness of these models and outline why this perspective will be helpful in the search for synthetic Turing-patterning systems. We conclude by considering robustness in the context of developmental modelling more generally. This article is part of the theme issue 'Recent progress and open frontiers in Turing's theory of morphogenesis'.

The design principles of discrete turing patterning systems

The formation of spatial structures lies at the heart of developmental processes. However, many of the underlying gene regulatory and biochemical processes remain poorly understood. Turing patterns constitute a main candidate to explain such processes, but they appear sensitive to fluctuations and variations in kinetic parameters, raising the question of how they may be adopted and realised in naturally evolved systems. The vast majority of mathematical studies of Turing patterns have used continuous models specified in terms of partial differential equations. Here, we complement this work by studying Turing patterns using discrete cellular automata models. We perform a large-scale study on all possible two-species networks and find the same Turing pattern producing networks as in the continuous framework. In contrast to continuous models, however, we find these Turing pattern topologies to be substantially more robust to changes in the parameters of the model. We also find that diffusion-driven instabilities are substantially weaker predictors for Turing patterns in our discrete modelling framework in comparison to the continuous case, in the sense that the presence of an instability does not guarantee a pattern emerging in simulations. We show that a more refined criterion constitutes a stronger predictor. The similarity of the results for the two modelling frameworks suggests a deeper underlying principle of Turing mechanisms in nature. Together with the larger robustness in the discrete case this suggests that Turing patterns may be more robust than previously thought.

Stochastic self-assembly of reaction-diffusion patterns in synaptic membranes

Synaptic receptor and scaffold molecules self-assemble into membrane protein domains, which play an important role in signal transmission across chemical synapses. Experiment and theory have shown that the formation of receptor-scaffold domains of the characteristic size observed in nerve cells can be understood from the receptor and scaffold reaction and diffusion processes suggested by experiments. We employ here kinetic Monte Carlo (KMC) simulations to explore the self-assembly of synaptic receptor-scaffold domains in a stochastic lattice model of receptor and scaffold reaction-diffusion dynamics. For reaction and diffusion rates within the ranges of values suggested by experiments we find, in agreement with previous mean-field calculations, self-assembly of receptor-scaffold domains of a size similar to that observed in experiments. Comparisons between the results of our KMC simulations and mean-field solutions suggest that the intrinsic noise associated with receptor and scaffold reaction and diffusion processes accelerates the self-assembly of receptor-scaffold domains, and confers increased robustness to domain formation. In agreement with experimental observations, our KMC simulations yield a prevalence of scaffolds over receptors in receptor-scaffold domains. Our KMC simulations show that receptor and scaffold reaction-diffusion dynamics can inherently give rise to plasticity in the overall properties of receptor-scaffold domains, which may be utilized by nerve cells to regulate the receptor number at chemical synapses.

Turing's Diffusive Threshold in Random Reaction-Diffusion Systems

Turing instabilities of reaction-diffusion systems can only arise if the diffusivities of the chemical species are sufficiently different. This threshold is unphysical in most systems with N=2 diffusing species, forcing experimental realizations of the instability to rely on fluctuations or additional nondiffusing species. Here, we ask whether this diffusive threshold lowers for N>2 to allow "true" Turing instabilities. Inspired by May's analysis of the stability of random ecological communities, we analyze the probability distribution of the diffusive threshold in reaction-diffusion systems defined by random matrices describing linearized dynamics near a homogeneous fixed point. In the numerically tractable cases N⩽6, we find that the diffusive threshold becomes more likely to be smaller and physical as N increases, and that most of these many-species instabilities cannot be described by reduced models with fewer diffusing species.

Plastic systemic inhibition controls amplitude while allowing phase pattern in a stochastic neural field model

We investigate oscillatory phase pattern formation and amplitude control for a linearized stochastic neuron field model by simulating Mexican-hat-coupled stochastic processes. We find, for several choices of parameters, that spatial pattern formation in the temporal phases of the coupled processes occurs if and only if their amplitudes are allowed to grow unrealistically large. Stimulated by recent work on homeostatic inhibitory plasticity, we introduce static and plastic (adaptive) systemic inhibitory mechanisms to keep the amplitudes stochastically bounded. We find that systems with static inhibition exhibited bounded amplitudes but no sustained phase patterns. With plastic systemic inhibition, on the other hand, the resulting systems exhibit both bounded amplitudes and sustained phase patterns. These results demonstrate that plastic inhibitory mechanisms in neural field models can dynamically control amplitudes while allowing patterns of phase synchronization to develop. Similar mechanisms of plastic systemic inhibition could play a role in regulating oscillatory functioning in the brain.

Stochastic Turing Pattern Formation in a Model with Active and Passive Transport

We investigate Turing pattern formation in a stochastic and spatially discretized version of a reaction-diffusion-advection (RDA) equation, which was previously introduced to model synaptogenesis in C. elegans. The model describes the interactions between a passively diffusing molecular species and an advecting species that switches between anterograde and retrograde motor-driven transport (bidirectional transport). Within the context of synaptogenesis, the diffusing molecules can be identified with the protein kinase CaMKII and the advecting molecules as glutamate receptors. The stochastic dynamics evolves according to an RDA master equation, in which advection and diffusion are both modeled as hopping reactions along a one-dimensional array of chemical compartments. Carrying out a linear noise approximation of the RDA master equation leads to an effective Langevin equation, whose power spectrum provides a means of extending the definition of a Turing instability to stochastic systems, namely in terms of the existence of a peak in the power spectrum at a nonzero spatial frequency. We thus show how noise can significantly extend the range over which spontaneous patterns occur, which is consistent with previous studies of RD systems.

Precision and dissipation of a stochastic Turing pattern

Spontaneous pattern formation is a fundamental scientific problem that has received much attention since the seminal theoretical work of Turing on reaction-diffusion systems. In molecular biophysics, this phenomenon often takes place under the influence of large fluctuations. It is then natural to inquire about the precision of such pattern. In particular, spontaneous pattern formation is a nonequilibrium phenomenon, and the relation between the precision of a pattern and the thermodynamic cost associated with it remains largely unexplored. Here, we analyze this relation with a paradigmatic stochastic reaction-diffusion model, i.e., the Brusselator in one spatial dimension. We find that the precision of the pattern is maximized for an intermediate thermodynamic cost, i.e., increasing the thermodynamic cost beyond this value makes the pattern less precise. Even though fluctuations get less pronounced with an increase in thermodynamic cost, we argue that larger fluctuations can also have a positive effect on the precision of the pattern.

Pattern formation mechanisms of self-organizing reaction-diffusion systems

Embryonic development is a largely self-organizing process, in which the adult body plan arises from a ball of cells with initially nearly equal potency. The reaction-diffusion theory first proposed by Alan Turing states that the initial symmetry in embryos can be broken by the interplay between two diffusible molecules, whose interactions lead to the formation of patterns. The reaction-diffusion theory provides a valuable framework for self-organized pattern formation, but it has been difficult to relate simple two-component models to real biological systems with multiple interacting molecular species. Recent studies have addressed this shortcoming and extended the reaction-diffusion theory to realistic multi-component networks. These efforts have challenged the generality of previous central tenets derived from the analysis of simplified systems and guide the way to a new understanding of self-organizing processes. Here, we discuss the challenges in modeling multi-component reaction-diffusion systems and how these have recently been addressed. We present a synthesis of new pattern formation mechanisms derived from these analyses, and we highlight the significance of reaction-diffusion principles for developmental and synthetic pattern formation.

Coloured Noise from Stochastic Inflows in Reaction-Diffusion Systems

In this paper, we present a framework for investigating coloured noise in reaction-diffusion systems. We start by considering a deterministic reaction-diffusion equation and show how external forcing can cause temporally correlated or coloured noise. Here, the main source of external noise is considered to be fluctuations in the parameter values representing the inflow of particles to the system. First, we determine which reaction systems, driven by extrinsic noise, can admit only one steady state, so that effects, such as stochastic switching, are precluded from our analysis. To analyse the steady-state behaviour of reaction systems, even if the parameter values are changing, necessitates a parameter-free approach, which has been central to algebraic analysis in chemical reaction network theory. To identify suitable models, we use tools from real algebraic geometry that link the network structure to its dynamical properties. We then make a connection to internal noise models and show how power spectral methods can be used to predict stochastically driven patterns in systems with coloured noise. In simple cases, we show that the power spectrum of the coloured noise process and the power spectrum of the reaction-diffusion system modelled with white noise multiply to give the power spectrum of the coloured noise reaction-diffusion system.

Interface Fingering Instability Triggered by a Density-Coupled Oscillatory Chemical Reaction via Precipitation

A density fingering hydrodynamic instability is triggered by a chemical reaction at the interface between two fluids. The density instability is controlled by the density gradient between both solutions, while the excitability of the bubble-free Belousov-Zhabotinsky-1,4-cyclohexanedione (BZ-CHD) oscillatory chemical reaction controls the importance of the chemistry in the system. Both parameters are thoroughly analyzed, and the mechanism underlying the instability is unveiled. The experimental observations lead us to modify the existing and accepted models for the BZ-CHD reaction within this context. The important role played by precipitation is considered in this context and included into the model. The modified kinetic model once coupled with fluid dynamics along with the precipitation mechanism was able to reproduce the experimental observations.

Noise sharing and Mexican-hat coupling in a stochastic neural field

A diffusion-type coupling operator that is biologically significant in neuroscience is a difference of Gaussian functions (Mexican-hat operator) used as a spatial-convolution kernel. We are interested in pattern formation by stochastic neural field equations, a class of space-time stochastic differential-integral equations using the Mexican-hat kernel. We explore quantitatively how the parameters that control the shape of the coupling kernel, the coupling strength, and aspects of spatially smoothed space-time noise influence the pattern in the resulting evolving random field. We confirm that a spatial pattern that is damped in time in a deterministic system may be sustained and amplified by stochasticity. We find that spatially smoothed noise alone causes pattern formation even without direct spatial coupling. Our analysis of the interaction between coupling and noise sharing allows us to determine parameter combinations that are optimal for the formation of spatial pattern.

Stochastic fluctuations and quasipattern formation in reaction-diffusion systems with anomalous transport

Many approaches to modeling reaction-diffusion systems with anomalous transport rely on deterministic equations which ignore fluctuations arising due to finite particle numbers. Starting from an individual-based model we use a generating-functional approach to derive a Gaussian approximation for this intrinsic noise in subdiffusive systems. This results in corrections to the deterministic fractional reaction-diffusion equations. Using this analytical approach, we study the onset of noise-driven quasipatterns in reaction-subdiffusion systems. We find that subdiffusion can be conducive to the formation of both deterministic and stochastic patterns. Our analysis shows that the combination of subdiffusion and intrinsic stochasticity can reduce the threshold ratio of the effective diffusion coefficients required for pattern formation to a greater degree than either effect on its own.

Spatial Interactions and Oscillatory Tragedies of the Commons

A tragedy of the commons (TOC) occurs when individuals acting in their own self-interest deplete commonly held resources, leading to a worse outcome than had they cooperated. Over time, the depletion of resources can change incentives for subsequent actions. Here, we investigate long-term feedback between game and environment across a continuum of incentives in an individual-based framework. We identify payoff-dependent transition rules that lead to oscillatory TOCs in stochastic simulations and the mean field limit. Further extending the stochastic model, we find that spatially explicit interactions can lead to emergent, localized dynamics, including the propagation of cooperative wave fronts and cluster formation of both social context and resources. These dynamics suggest new mechanisms underlying how TOCs arise and how they might be averted.

Analysing diffusion and flow-driven instability using semidefinite programming

Diffusion and flow-driven instability, or transport-driven instability, is one of the central mechanisms to generate inhomogeneous gradient of concentrations in spatially distributed chemical systems. However, verifying the transport-driven instability of reaction-diffusion-advection systems requires checking the Jacobian eigenvalues of infinitely many Fourier modes, which is computationally intractable. To overcome this limitation, this paper proposes mathematical optimization algorithms that determine the stability/instability of reaction-diffusion-advection systems by finite steps of algebraic calculations. Specifically, the stability/instability analysis of Fourier modes is formulated as a sum-of-squares optimization program, which is a class of convex optimization whose solvers are widely available as software packages. The optimization program is further extended for facile computation of the destabilizing spatial modes. This extension allows for predicting and designing the shape of the concentration gradient without simulating the governing equations. The streamlined analysis process of self-organized pattern formation is demonstrated with a simple illustrative reaction model with diffusion and advection.

Noise⁻Seeded Developmental Pattern Formation in Filamentous Cyanobacteria

Under nitrogen-poor conditions, multicellular cyanobacteria such as Anabaena sp. PCC 7120 undergo a process of differentiation, forming nearly regular, developmental patterns of individual nitrogen-fixing cells, called heterocysts, interspersed between intervals of vegetative cells that carry out photosynthesis. Developmental pattern formation is mediated by morphogen species that can act as activators and inhibitors, some of which can diffuse along filaments. We survey the limitations of the classical, deterministic Turing mechanism that has been often invoked to explain pattern formation in these systems, and then, focusing on a simpler system governed by birth-death processes, we illustrate pedagogically a recently proposed paradigm that provides a much more robust description of pattern formation: stochastic Turing patterns. We emphasize the essential role that cell-to-cell differences in molecular numbers—caused by inevitable fluctuations in gene expression—play, the so called demographic noise, in seeding the formation of stochastic Turing patterns over a much larger region of parameter space, compared to their deterministic counterparts.

Stochastic Turing patterns in a synthetic bacterial population

The origin of biological morphology and form is one of the deepest problems in science, underlying our understanding of development and the functioning of living systems. In 1952, Alan Turing showed that chemical morphogenesis could arise from a linear instability of a spatially uniform state, giving rise to periodic pattern formation in reaction-diffusion systems but only those with a rapidly diffusing inhibitor and a slowly diffusing activator. These conditions are disappointingly hard to achieve in nature, and the role of Turing instabilities in biological pattern formation has been called into question. Recently, the theory was extended to include noisy activator-inhibitor birth and death processes. Surprisingly, this stochastic Turing theory predicts the existence of patterns over a wide range of parameters, in particular with no severe requirement on the ratio of activator-inhibitor diffusion coefficients. To explore whether this mechanism is viable in practice, we have genetically engineered a synthetic bacterial population in which the signaling molecules form a stochastic activator-inhibitor system. The synthetic pattern-forming gene circuit destabilizes an initially homogenous lawn of genetically engineered bacteria, producing disordered patterns with tunable features on a spatial scale much larger than that of a single cell. Spatial correlations of the experimental patterns agree quantitatively with the signature predicted by theory. These results show that Turing-type pattern-forming mechanisms, if driven by stochasticity, can potentially underlie a broad range of biological patterns. These findings provide the groundwork for a unified picture of biological morphogenesis, arising from a combination of stochastic gene expression and dynamical instabilities.

Robust stochastic Turing patterns in the development of a one-dimensional cyanobacterial organism

Under nitrogen deprivation, the one-dimensional cyanobacterial organism Anabaena sp. PCC 7120 develops patterns of single, nitrogen-fixing cells separated by nearly regular intervals of photosynthetic vegetative cells. We study a minimal, stochastic model of developmental patterns in Anabaena that includes a nondiffusing activator, two diffusing inhibitor morphogens, demographic fluctuations in the number of morphogen molecules, and filament growth. By tracking developing filaments, we provide experimental evidence for different spatiotemporal roles of the two inhibitors during pattern maintenance and for small molecular copy numbers, justifying a stochastic approach. In the deterministic limit, the model yields Turing patterns within a region of parameter space that shrinks markedly as the inhibitor diffusivities become equal. Transient, noise-driven, stochastic Turing patterns are produced outside this region, which can then be fixed by downstream genetic commitment pathways, dramatically enhancing the robustness of pattern formation, also in the biologically relevant situation in which the inhibitors' diffusivities may be comparable.

Coevolution Maintains Diversity in the Stochastic "Kill the Winner" Model

The "kill the winner" hypothesis is an attempt to address the problem of diversity in biology. It argues that host-specific predators control the population of each prey, preventing a winner from emerging and thus maintaining the coexistence of all species in the system. We develop a stochastic model for the kill the winner paradigm and show that the stable coexistence state of the deterministic kill the winner model is destroyed by demographic stochasticity, through a cascade of extinction events. We formulate an individual-level stochastic model in which predator-prey coevolution promotes the high diversity of the ecosystem by generating a persistent population flux of species.

Distribution of randomly diffusing particles in inhomogeneous media

Diffusion can be conceptualized, at microscopic scales, as the random hopping of particles between neighboring lattice sites. In the case of diffusion in inhomogeneous media, distinct spatial domains in the system may yield distinct particle hopping rates. Starting from the master equations (MEs) governing diffusion in inhomogeneous media we derive here, for arbitrary spatial dimensions, the deterministic lattice equations (DLEs) specifying the average particle number at each lattice site for randomly diffusing particles in inhomogeneous media. We consider the case of free (Fickian) diffusion with no steric constraints on the maximum particle number per lattice site as well as the case of diffusion under steric constraints imposing a maximum particle concentration. We find, for both transient and asymptotic regimes, excellent agreement between the DLEs and kinetic Monte Carlo simulations of the MEs. The DLEs provide a computationally efficient method for predicting the (average) distribution of randomly diffusing particles in inhomogeneous media, with the number of DLEs associated with a given system being independent of the number of particles in the system. From the DLEs we obtain general analytic expressions for the steady-state particle distributions for free diffusion and, in special cases, diffusion under steric constraints in inhomogeneous media. We find that, in the steady state of the system, the average fraction of particles in a given domain is independent of most system properties, such as the arrangement and shape of domains, and only depends on the number of lattice sites in each domain, the particle hopping rates, the number of distinct particle species in the system, and the total number of particles of each particle species in the system. Our results provide general insights into the role of spatially inhomogeneous particle hopping rates in setting the particle distributions in inhomogeneous media.

Motif analysis for small-number effects in chemical reaction dynamics

The number of molecules involved in a cell or subcellular structure is sometimes rather small. In this situation, ordinary macroscopic-level fluctuations can be overwhelmed by non-negligible large fluctuations, which results in drastic changes in chemical-reaction dynamics and statistics compared to those observed under a macroscopic system (i.e., with a large number of molecules). In order to understand how salient changes emerge from fluctuations in molecular number, we here quantitatively define small-number effect by focusing on a "mesoscopic" level, in which the concentration distribution is distinguishable both from micro- and macroscopic ones and propose a criterion for determining whether or not such an effect can emerge in a given chemical reaction network. Using the proposed criterion, we systematically derive a list of motifs of chemical reaction networks that can show small-number effects, which includes motifs showing emergence of the power law and the bimodal distribution observable in a mesoscopic regime with respect to molecule number. The list of motifs provided herein is helpful in the search for candidates of biochemical reactions with a small-number effect for possible biological functions, as well as for designing a reaction system whose behavior can change drastically depending on molecule number, rather than concentration.

Stochastic lattice model of synaptic membrane protein domains

Neurotransmitter receptor molecules, concentrated in synaptic membrane domains along with scaffolds and other kinds of proteins, are crucial for signal transmission across chemical synapses. In common with other membrane protein domains, synaptic domains are characterized by low protein copy numbers and protein crowding, with rapid stochastic turnover of individual molecules. We study here in detail a stochastic lattice model of the receptor-scaffold reaction-diffusion dynamics at synaptic domains that was found previously to capture, at the mean-field level, the self-assembly, stability, and characteristic size of synaptic domains observed in experiments. We show that our stochastic lattice model yields quantitative agreement with mean-field models of nonlinear diffusion in crowded membranes. Through a combination of analytic and numerical solutions of the master equation governing the reaction dynamics at synaptic domains, together with kinetic Monte Carlo simulations, we find substantial discrepancies between mean-field and stochastic models for the reaction dynamics at synaptic domains. Based on the reaction and diffusion properties of synaptic receptors and scaffolds suggested by previous experiments and mean-field calculations, we show that the stochastic reaction-diffusion dynamics of synaptic receptors and scaffolds provide a simple physical mechanism for collective fluctuations in synaptic domains, the molecular turnover observed at synaptic domains, key features of the observed single-molecule trajectories, and spatial heterogeneity in the effective rates at which receptors and scaffolds are recycled at the cell membrane. Our work sheds light on the physical mechanisms and principles linking the collective properties of membrane protein domains to the stochastic dynamics that rule their molecular components.

Giant Amplification of Noise in Fluctuation-Induced Pattern Formation

The amplitude of fluctuation-induced patterns might be expected to be proportional to the strength of the driving noise, suggesting that such patterns would be difficult to observe in nature. Here, we show that a large class of spatially extended dynamical systems driven by intrinsic noise can exhibit giant amplification, yielding patterns whose amplitude is comparable to that of deterministic Turing instabilities. The giant amplification results from the interplay between noise and nonorthogonal eigenvectors of the linear stability matrix, yielding transients that grow with time, and which, when driven by the ever-present intrinsic noise, lead to persistent large amplitude patterns. This mechanism shows that fluctuation-induced Turing patterns are observable, and are not strongly limited by the amplitude of demographic stochasticity nor by the value of the diffusion coefficients.

Intrinsic noise in systems with switching environments

We study individual-based dynamics in finite populations, subject to randomly switching environmental conditions. These are inspired by models in which genes transition between on and off states, regulating underlying protein dynamics. Similarly, switches between environmental states are relevant in bacterial populations and in models of epidemic spread. Existing piecewise-deterministic Markov process approaches focus on the deterministic limit of the population dynamics while retaining the randomness of the switching. Here we go beyond this approximation and explicitly include effects of intrinsic stochasticity at the level of the linear-noise approximation. Specifically, we derive the stationary distributions of a number of model systems, in good agreement with simulations. This improves existing approaches which are limited to the regimes of fast and slow switching.

Turing Patterning Using Gene Circuits with Gas-Induced Degradation of Quorum Sensing Molecules

The Turing instability was proposed more than six decades ago as a mechanism leading to spatial patterning, but it has yet to be exploited in a synthetic biology setting. Here we characterize the Turing instability in a specific gene circuit that can be implemented in vitro or in populations of clonal cells producing short-range activator N-Acyl homoserine lactone (AHL) and long-range inhibitor hydrogen peroxide (H2O2) gas. Slowing the production rate of the AHL-degrading enzyme, AiiA, generates stable fixed states, limit cycle oscillations and Turing patterns. Further tuning of signaling parameters determines local robustness and controls the range of unstable wavenumbers in the patterning regime. These findings provide a roadmap for optimizing spatial patterns of gene expression based on familiar quorum and gas sensitive E. coli promoters. The circuit design and predictions may be useful for (re)programming spatial dynamics in synthetic and natural gene expression systems.

Positive feedback can lead to dynamic nanometer-scale clustering on cell membranes

Clustering of molecules on biological membranes is a widely observed phenomenon. A key example is the clustering of the oncoprotein Ras, which is known to be important for signal transduction in mammalian cells. Yet, the mechanism by which Ras clusters form and are maintained remains unclear. Recently, it has been discovered that activated Ras promotes further Ras activation. Here we show using particle-based simulation that this positive feedback is sufficient to produce persistent clusters of active Ras molecules at the nanometer scale via a dynamic nucleation mechanism. Furthermore, we find that our cluster statistics are consistent with experimental observations of the Ras system. Interestingly, we show that our model does not support a Turing regime of macroscopic reaction-diffusion patterning, and therefore that the clustering we observe is a purely stochastic effect, arising from the coupling of positive feedback with the discrete nature of individual molecules. These results underscore the importance of stochastic and dynamic properties of reaction diffusion systems for biological behavior.

Small-scale properties of a stochastic cubic-autocatalytic reaction-diffusion model

We investigate the small-scale properties of a stochastic cubic-autocatalytic reaction-diffusion (CARD) model using renormalization techniques. We renormalize noise-induced ultraviolet divergences and obtain β functions for the decay rate and coupling at one loop. Assuming colored (power-law) noise, our results show that the behavior of both decay rate and coupling with scale depends crucially on the noise exponent. Interpreting the CARD model as a proxy for a (very simple) living system, our results suggest that power-law correlations in environmental fluctuations can both decrease or increase the growth of structures at smaller scales.

Self-assembly and plasticity of synaptic domains through a reaction-diffusion mechanism

Signal transmission across chemical synapses relies crucially on neurotransmitter receptor molecules, concentrated in postsynaptic membrane domains along with scaffold and other postsynaptic molecules. The strength of the transmitted signal depends on the number of receptor molecules in postsynaptic domains, and activity-induced variation in the receptor number is one of the mechanisms of postsynaptic plasticity. Recent experiments have demonstrated that the reaction and diffusion properties of receptors and scaffolds at the membrane, alone, yield spontaneous formation of receptor-scaffold domains of the stable characteristic size observed in neurons. On the basis of these experiments we develop a model describing synaptic receptor domains in terms of the underlying reaction-diffusion processes. Our model predicts that the spontaneous formation of receptor-scaffold domains of the stable characteristic size observed in experiments depends on a few key reactions between receptors and scaffolds. Furthermore, our model suggests novel mechanisms for the alignment of pre- and postsynaptic domains and for short-term postsynaptic plasticity in receptor number. We predict that synaptic receptor domains localize in membrane regions with an increased receptor diffusion coefficient or a decreased scaffold diffusion coefficient. Similarly, we find that activity-dependent increases or decreases in receptor or scaffold diffusion yield a transient increase in the number of receptor molecules concentrated in postsynaptic domains. Thus, the proposed reaction-diffusion model puts forth a coherent set of biophysical mechanisms for the formation, stability, and plasticity of molecular domains on the postsynaptic membrane.

Path-integral methods for analyzing the effects of fluctuations in stochastic hybrid neural networks

We consider applications of path-integral methods to the analysis of a stochastic hybrid model representing a network of synaptically coupled spiking neuronal populations. The state of each local population is described in terms of two stochastic variables, a continuous synaptic variable and a discrete activity variable. The synaptic variables evolve according to piecewise-deterministic dynamics describing, at the population level, synapses driven by spiking activity. The dynamical equations for the synaptic currents are only valid between jumps in spiking activity, and the latter are described by a jump Markov process whose transition rates depend on the synaptic variables. We assume a separation of time scales between fast spiking dynamics with time constant [Formula: see text] and slower synaptic dynamics with time constant τ. This naturally introduces a small positive parameter [Formula: see text], which can be used to develop various asymptotic expansions of the corresponding path-integral representation of the stochastic dynamics. First, we derive a variational principle for maximum-likelihood paths of escape from a metastable state (large deviations in the small noise limit [Formula: see text]). We then show how the path integral provides an efficient method for obtaining a diffusion approximation of the hybrid system for small ϵ. The resulting Langevin equation can be used to analyze the effects of fluctuations within the basin of attraction of a metastable state, that is, ignoring the effects of large deviations. We illustrate this by using the Langevin approximation to analyze the effects of intrinsic noise on pattern formation in a spatially structured hybrid network. In particular, we show how noise enlarges the parameter regime over which patterns occur, in an analogous fashion to PDEs. Finally, we carry out a [Formula: see text]-loop expansion of the path integral, and use this to derive corrections to voltage-based mean-field equations, analogous to the modified activity-based equations generated from a neural master equation.

Theoretical analysis of discreteness-induced transition in autocatalytic reaction dynamics

Transitions in the qualitative behavior of chemical reaction dynamics with a decrease in molecule number have attracted much attention. Here, a method based on a Markov process with a tridiagonal transition matrix is applied to the analysis of this transition in reaction dynamics. The transition to bistability due to the small-number effect and the mean switching time between the bistable states are analytically calculated in agreement with numerical simulations. In addition, a novel transition involving the reversal of the chemical reaction flow is found in the model under an external flow, and also in a three-component model. The generality of this transition and its correspondence to biological phenomena are also discussed.

Stochastic Turing patterns: analysis of compartment-based approaches

Turing patterns can be observed in reaction-diffusion systems where chemical species have different diffusion constants. In recent years, several studies investigated the effects of noise on Turing patterns and showed that the parameter regimes, for which stochastic Turing patterns are observed, can be larger than the parameter regimes predicted by deterministic models, which are written in terms of partial differential equations (PDEs) for species concentrations. A common stochastic reaction-diffusion approach is written in terms of compartment-based (lattice-based) models, where the domain of interest is divided into artificial compartments and the number of molecules in each compartment is simulated. In this paper, the dependence of stochastic Turing patterns on the compartment size is investigated. It has previously been shown (for relatively simpler systems) that a modeler should not choose compartment sizes which are too small or too large, and that the optimal compartment size depends on the diffusion constant. Taking these results into account, we propose and study a compartment-based model of Turing patterns where each chemical species is described using a different set of compartments. It is shown that the parameter regions where spatial patterns form are different from the regions obtained by classical deterministic PDE-based models, but they are also different from the results obtained for the stochastic reaction-diffusion models which use a single set of compartments for all chemical species. In particular, it is argued that some previously reported results on the effect of noise on Turing patterns in biological systems need to be reinterpreted.

Path-integral calculation for the emergence of rapid evolution from demographic stochasticity

Genetic variation in a population can sometimes arise so fast as to modify ecosystem dynamics. Such phenomena have been observed in natural predator-prey systems and characterized in the laboratory as showing unusual phase relationships in population dynamics, including a π phase shift between predator and prey (evolutionary cycles) and even undetectable prey oscillations compared to those of the predator (cryptic cycles). Here we present a generic individual-level stochastic model of interacting populations that includes a subpopulation of low nutritional value to the predator. Using a master equation formalism and by mapping to a coherent state path integral solved by a system-size expansion, we show that evolutionary and cryptic quasicycles can emerge generically from the combination of intrinsic demographic fluctuations and clonal mutations alone, without additional biological mechanisms.

Characterization of spiraling patterns in spatial rock-paper-scissors games

The spatiotemporal arrangement of interacting populations often influences the maintenance of species diversity and is a subject of intense research. Here, we study the spatiotemporal patterns arising from the cyclic competition between three species in two dimensions. Inspired by recent experiments, we consider a generic metapopulation model comprising "rock-paper-scissors" interactions via dominance removal and replacement, reproduction, mutations, pair exchange, and hopping of individuals. By combining analytical and numerical methods, we obtain the model's phase diagram near its Hopf bifurcation and quantitatively characterize the properties of the spiraling patterns arising in each phase. The phases characterizing the cyclic competition away from the Hopf bifurcation (at low mutation rate) are also investigated. Our analytical approach relies on the careful analysis of the properties of the complex Ginzburg-Landau equation derived through a controlled (perturbative) multiscale expansion around the model's Hopf bifurcation. Our results allow us to clarify when spatial "rock-paper-scissors" competition leads to stable spiral waves and under which circumstances they are influenced by nonlinear mobility.

Feedback, receptor clustering, and receptor restriction to single cells yield large Turing spaces for ligand-receptor-based Turing models

Turing mechanisms can yield a large variety of patterns from noisy, homogenous initial conditions and have been proposed as patterning mechanism for many developmental processes. However, the molecular components that give rise to Turing patterns have remained elusive, and the small size of the parameter space that permits Turing patterns to emerge makes it difficult to explain how Turing patterns could evolve. We have recently shown that Turing patterns can be obtained with a single ligand if the ligand-receptor interaction is taken into account. Here we show that the general properties of ligand-receptor systems result in very large Turing spaces. Thus, the restriction of receptors to single cells, negative feedbacks, regulatory interactions among different ligand-receptor systems, and the clustering of receptors on the cell surface all greatly enlarge the Turing space. We further show that the feedbacks that occur in the FGF10-SHH network that controls lung branching morphogenesis are sufficient to result in large Turing spaces. We conclude that the cellular restriction of receptors provides a mechanism to sufficiently increase the size of the Turing space to make the evolution of Turing patterns likely. Additional feedbacks may then have further enlarged the Turing space. Given their robustness and flexibility, we propose that receptor-ligand-based Turing mechanisms present a general mechanism for patterning in biology.

Coexistence of productive and non-productive populations by fluctuation-driven spatio-temporal patterns

Cooperative interactions, their stability and evolution, provide an interesting context in which to study the interface between cellular and population levels of organization. Here we study a public goods model relevant to microorganism populations actively extracting a growth resource from their environment. Cells can display one of two phenotypes - a productive phenotype that extracts the resources at a cost, and a non-productive phenotype that only consumes the same resource. Both proliferate and are free to move by diffusion; growth rate and diffusion coefficient depend only weakly phenotype. We analyze the continuous differential equation model as well as simulate stochastically the full dynamics. We find that the two sub-populations, which cannot coexist in a well-mixed environment, develop spatio-temporal patterns that enable long-term coexistence in the shared environment. These patterns are purely fluctuation-driven, as the corresponding continuous spatial system does not display Turing instability. The average stability of coexistence patterns derives from a dynamic mechanism in which the producing sub-population equilibrates with the environmental resource and holds it close to an extinction transition of the other sub-population, causing it to constantly hover around this transition. Thus the ecological interactions support a mechanism reminiscent of self-organized criticality; power-law distributions and long-range correlations are found. The results are discussed in the context of general pattern formation and critical behavior in ecology as well as in an experimental context.

Linear noise approximation for stochastic oscillations of intracellular calcium

A stochastic model of intracellular calcium oscillations is analytically studied. The governing master equation is expanded under the linear noise approximation and a closed prediction for the power spectrum of fluctuations analytically derived. A peak in the obtained power spectrum profile signals the presence of stochastic, noise induced oscillations which extend also outside the region where a deterministic limit cycle is predicted to occur.

A GRAPH PARTITIONING APPROACH TO PREDICTING PATTERNS IN LATERAL INHIBITION SYSTEMS

We analyze spatial patterns on networks of cells where adjacent cells inhibit each other through contact signaling. We represent the network as a graph where each vertex represents the dynamics of identical individual cells and where graph edges represent cell-to-cell signaling. To predict steady-state patterns we find equitable partitions of the graph vertices and assign them into disjoint classes. We then use results from monotone systems theory to prove the existence of patterns that are structured in such a way that all the cells in the same class have the same final fate. To study the stability properties of these patterns, we rely on the graph partition to perform a block decomposition of the system. Then, to guarantee stability, we provide a small-gain type criterion that depends on the input-output properties of each cell in the reduced system. Finally, we discuss pattern formation in stochastic models. With the help of a modal decomposition we show that noise can enhance the parameter region where patterning occurs.