Ratio control in a cascade model of cell differentiation

We propose a kind of reaction-diffusion equations for cell differentiation, which exhibits the Turing instability. If the diffusivity of some variables is set to be infinity, we get coupled competitive reaction-diffusion equations with a global feedback term. The size ratio of each cell type is controlled by a system parameter in the model. Finally, we extend the model to a cascade model of cell differentiation. A hierarchical spatial structure appears as a result of the cell differentiation. The size ratio of each cell type is also controlled by the system parameter.Comment: 13 pages, 7 figure

Predicting Axonal Response to Molecular Gradients with a Computational Model of Filopodial Dynamics

Axons are often guided to their targets in the developing nervous system by attractive or repulsive molecular concentration gradients. We propose a computational model for gradient sensing and directed movement of the growth cone mediated by filopodia. We show that relatively simple mechanisms are sufficient to generate realistic rajectories for both the short-term response of axons to steep gradients and the long-term response of axons to shallow gradients. The model makes testable predictions for axonal response to attractive and repulsive gradients of different concentrations and steepness, the size of the intracellular amplification of the gradient signal, and the differences in intracellular signaling required for repulsive versus attractive turning

Helical Turing patterns in the Lengyel-Epstein model in thin cylindrical layers

The formation of Turing patterns was investigated in thin cylindrical layers using the Lengyel-Epstein model of the chlorine dioxide-iodine-malonic acid reaction. The influence of the width of the layer W and the diameter D of the inner cylinder on the pattern with intrinsic wavelength l were determined in simulations with initial random noise perturbations to the uniform state for W< l/2 and D l or lower. We show that the geometric constraints of the reaction domain may result in the formation of helical Turing patterns with parameters that give stripes (b ¼ 0.2) or spots (b ¼ 0.37) in two dimensions. For b ¼ 0.2, the helices were composed of lamellae and defects were likely as the diameter of the cylinder increased. With b ¼ 0.37, the helices consisted of semi-cylinders and the orientation of stripes on the outer surface (and hence winding number) increased with increasing diameter until a new stripe appeared

Finding the center reliably: robust patterns of developmental gene expression

We investigate a mechanism for the robust identification of the center of a developing biological system. We assume the existence of two morphogen gradients, an activator emanating from the anterior, and a co-repressor from the posterior. The co-repressor inhibits the action of the activator in switching on target genes. We apply this system to Drosophila embryos, where we predict the existence of a hitherto undetected posterior co-repressor. Using mathematical modelling, we show that a symmetric activator-co-repressor model can quantitatively explain the precise mid-embryo expression boundary of the hunchback gene, and the scaling of this pattern with embryo size.Comment: 4 pages, 3 figure

Morphogen Transport in Epithelia

We present a general theoretical framework to discuss mechanisms of morphogen transport and gradient formation in a cell layer. Trafficking events on the cellular scale lead to transport on larger scales. We discuss in particular the case of transcytosis where morphogens undergo repeated rounds of internalization into cells and recycling. Based on a description on the cellular scale, we derive effective nonlinear transport equations in one and two dimensions which are valid on larger scales. We derive analytic expressions for the concentration dependence of the effective diffusion coefficient and the effective degradation rate. We discuss the effects of a directional bias on morphogen transport and those of the coupling of the morphogen and receptor kinetics. Furthermore, we discuss general properties of cellular transport processes such as the robustness of gradients and relate our results to recent experiments on the morphogen Decapentaplegic (Dpp) that acts in the fruit fly Drosophila

Stability of spikes in the shadow Gierer-Meinhardt system with Robin boundary conditions

We consider the shadow system of the Gierer-Meinhardt system in a smooth bounded domain RN,At=2A−A+,x, t>0, ||t=−||+Ardx, t>0 with the Robin boundary condition +aAA=0, x, where aA>0, the reaction rates (p,q,r,s) satisfy 1<p<()+, q>0, r>0, s0, 1<<+, the diffusion constant is chosen such that 1, and the time relaxation constant is such that 0. We rigorously prove the following results on the stability of one-spike solutions: (i) If r=2 and 1<p<1+4/N or if r=p+1 and 1<p<, then for aA>1 and sufficiently small the interior spike is stable. (ii) For N=1 if r=2 and 1<p3 or if r=p+1 and 1<p<, then for 0<aA<1 the near-boundary spike is stable. (iii) For N=1 if 3<p<5 and r=2, then there exist a0(0,1) and µ0>1 such that for a(a0,1) and µ=2q/(s+1)(p−1)(1,µ0) the near-boundary spike solution is unstable. This instability is not present for the Neumann boundary condition but only arises for the Robin boundary condition. Furthermore, we show that the corresponding eigenvalue is of order O(1) as 0. ©2007 American Institute of Physic

Modelling cell motility and chemotaxis with evolving surface finite elements

We present a mathematical and a computational framework for the modelling of cell motility. The cell membrane is represented by an evolving surface, with the movement of the cell determined by the interaction of various forces that act normal to the surface. We consider external forces such as those that may arise owing to inhomogeneities in the medium and a pressure that constrains the enclosed volume, as well as internal forces that arise from the reaction of the cells' surface to stretching and bending. We also consider a protrusive force associated with a reaction-diffusion system (RDS) posed on the cell membrane, with cell polarization modelled by this surface RDS. The computational method is based on an evolving surface finite-element method. The general method can account for the large deformations that arise in cell motility and allows the simulation of cell migration in three dimensions. We illustrate applications of the proposed modelling framework and numerical method by reporting on numerical simulations of a model for eukaryotic chemotaxis and a model for the persistent movement of keratocytes in two and three space dimensions. Movies of the simulated cells can be obtained from http://homepages.warwick.ac.uk/maskae/CV_Warwick/Chemotaxis.html

Noise-induced inhibitory suppression of malfunction neural oscillators

Motivated by the aim to find new medical strategies to suppress undesirable neural synchronization we study the control of oscillations in a system of inhibitory coupled noisy oscillators. Using dynamical properties of inhibition, we find regimes when the malfunction oscillations can be suppressed but the information signal of a certain frequency can be transmitted through the system. The mechanism of this phenomenon is a resonant interplay of noise and the transmission signal provided by certain value of inhibitory coupling. Analyzing a system of three or four oscillators representing neural clusters, we show that this suppression can be effectively controlled by coupling and noise amplitudes.Comment: 10 pages, 14 figure

Sierpinski signal generates 1/fα1/f^\alpha spectra

We investigate the row sum of the binary pattern generated by the Sierpinski automaton: Interpreted as a time series we calculate the power spectrum of this Sierpinski signal analytically and obtain a unique rugged fine structure with underlying power law decay with an exponent of approximately 1.15. Despite the simplicity of the model, it can serve as a model for $1/f^\alpha$ spectra in a certain class of experimental and natural systems like catalytic reactions and mollusc patterns.Comment: 4 pages (4 figs included). Accepted for publication in Physical Review