scholarly journals High-throughput mathematical analysis identifies Turing networks for patterning with equally diffusing signals

eLife ◽  
2016 ◽  
Vol 5 ◽  
Author(s):  
Luciano Marcon ◽  
Xavier Diego ◽  
James Sharpe ◽  
Patrick Müller

The Turing reaction-diffusion model explains how identical cells can self-organize to form spatial patterns. It has been suggested that extracellular signaling molecules with different diffusion coefficients underlie this model, but the contribution of cell-autonomous signaling components is largely unknown. We developed an automated mathematical analysis to derive a catalog of realistic Turing networks. This analysis reveals that in the presence of cell-autonomous factors, networks can form a pattern with equally diffusing signals and even for any combination of diffusion coefficients. We provide a software (available at http://www.RDNets.com) to explore these networks and to constrain topologies with qualitative and quantitative experimental data. We use the software to examine the self-organizing networks that control embryonic axis specification and digit patterning. Finally, we demonstrate how existing synthetic circuits can be extended with additional feedbacks to form Turing reaction-diffusion systems. Our study offers a new theoretical framework to understand multicellular pattern formation and enables the wide-spread use of mathematical biology to engineer synthetic patterning systems.

2010 ◽  
Vol 7 (1) ◽  
Author(s):  
Paola Lecca ◽  
Adaoha E. C. Ihekwaba ◽  
Lorenzo Dematté ◽  
Corrado Priami

SummaryReaction-diffusion systems are mathematical models that describe how the concentrations of substances distributed in space change under the influence of local chemical reactions, and diffusion which causes the substances to spread out in space. The classical representation of a reaction-diffusion system is given by semi-linear parabolic partial differential equations, whose solution predicts how diffusion causes the concentration field to change with time. This change is proportional to the diffusion coefficient. If the solute moves in a homogeneous system in thermal equilibrium, the diffusion coefficients are constants that do not depend on the local concentration of solvent and solute. However, in nonhomogeneous and structured media the assumption of constant intracellular diffusion coefficient is not necessarily valid, and, consequently, the diffusion coefficient is a function of the local concentration of solvent and solutes. In this paper we propose a stochastic model of reaction-diffusion systems, in which the diffusion coefficients are function of the local concentration, viscosity and frictional forces. We then describe the software tool Redi (REaction-DIffusion simulator) which we have developed in order to implement this model into a Gillespie-like stochastic simulation algorithm. Finally, we show the ability of our model implemented in the Redi tool to reproduce the observed gradient of the bicoid protein in the Drosophila Melanogaster embryo. With Redi, we were able to simulate with an accuracy of 1% the experimental spatio-temporal dynamics of the bicoid protein, as recorded in time-lapse experiments obtained by direct measurements of transgenic bicoidenhanced green fluorescent protein.


2018 ◽  
Vol 11 (04) ◽  
pp. 1850053 ◽  
Author(s):  
Laura Murphy ◽  
Chandrasekhar Venkataraman ◽  
Anotida Madzvamuse

We present a computational framework for isolating spatial patterns arising in the steady states of reaction–diffusion systems. Such systems have been used to model many natural phenomena in areas such as developmental and cancer biology, cell motility and material science. In many of these applications, often one is interested in identifying parameters which will lead to a particular pattern for a given reaction–diffusion model. To attempt to answer this, we compute eigenpairs of the Laplacian on a variety of domains and use linear stability analysis to determine parameter values for the system that will lead to spatially inhomogeneous steady states whose patterns correspond to particular eigenfunctions. This method has previously been used on domains and surfaces where the eigenvalues and eigenfunctions are found analytically in closed form. Our contribution to this methodology is that we numerically compute eigenpairs on arbitrary domains and surfaces. Here we present examples and demonstrate that mode isolation is straightforward especially for low eigenvalues. Additionally, we show that in some cases the inhomogeneous steady state can be a linear combination of eigenfunctions. Finally, we show an example suggesting that pattern formation is robust on similar surfaces in cases that the surface either has or does not have a boundary.


Sign in / Sign up

Export Citation Format

Share Document