IGEM:Cambridge/2008/Turing Pattern Formation/Modelling
Modelling ReactionDiffusion SystemsMathematically, reactiondiffusion systems are coupled nonlinear differential equations that can be solved numerically. Our first step will be to implement/model a simple twocomponent system originally proposed by Turing. After validating our numerical method and computating our first patterns (img to the right), we will be thinking about a more realistic system describing the behaviour of the activator and inhibitor system that we intend to engineer with B.subtilis. This will include an analysis of enzyme kinetics and we hope to deduce the parameter ranges, in which our B.subtilis construct will be able to form patterns, thus feeding back into our design decisions (promoter strength, rbs choice etc.). We will also take an investigative approach and ask whether Turinglike patterns can originate from a system that is simpler than the one envisioned, consisting of fewer components, e.g. depending on a single signalling molecule only. IntroductionHow does pattern formation occur? Turing considered diffusion to be the crucial component. He showed that pattern formation occurs of the system has a stable steady state in the absence of diffusion, but allows unstable states to develop when the diffusion term is added. This results in divergence (or, in a biological context, gene expression/cell differentiation) and subsequent pattern formation. In the following sections, we shall deduce conditions for pattern formation for this simple system. We will work with a twocomponent system and we nondimensionalise the variables. The resulting general form is:
where [math]\gamma[/math] (determining the scale) is a constant and d stands for the diffusion ratio. We note that nondimensionalisation has the added advantage that we can now map a specific pattern onto a wide range of biological parameters, the easiest example being that two different pairs of diffusion rates will result in the same pattern if the respective diffusion ratios remain unchanged. Turing System(also known as Schnakenberg reaction) Case: Without diffusion, need stable steady stateConsider the following system without diffusion terms:
We determine the steady state solution and add a small pertubation [math]\tilde{A}, \tilde{B} [/math] to linearise the system about the steady state in order to determine its stability. In Matrix notation, the linearised system can be written as:
where the Jacobian is evaluated at the steady states of A and B. This is a set of coupled firstorder ODEs and solutions are proportional to [math]exp(\lambda t)[/math]. In order to have a stable steady state, we need the real part of [math]\lambda[/math] to be negative. This requires [math]f_A + g_B \lt 0[/math] (need ve trace) and [math]f_A g_Bf_B g_A \gt 0[/math] (need +ve determinant), which we note as our first two conditions for pattern formation. Case: With diffusion, need unstable steady statesNow we add diffusion to our system and require the resulting steady state to be unstable. After linearisation, the equation can be expressed as follows:
In the following steps, we consider solutions of this particular form [math]\mathbf{w}(\mathbf{x},t) = \sum_{k} c_{k} e^{\lambda t} \mathbf{v}_k(\mathbf{x}) [/math] with [math] \mathbf{v}_k [/math] satisfying [math] \nabla^2 \mathbf{v}_k + k^2\mathbf{v}_k = 0[/math].
Given a (twocomponent) system, we can now infer which parameters to choose in order for the system to develop patterns. This saves us a lot of time as reaction diffusion systems are extremely sensitive to parameter changes and simply trying out sets of parameters will usually result in no patterning at all. Bacillus and its signalling systemsWe will first look at the specific components of our construct and models describing the behaviour of the agr and lux quorumsensing systems will be introduced. We will then model our envisioned activatorinhibitor system and deduce under which conditions pattern formation will occur. Peptidesignalling systemReceiverWe have created a simple model describing the agrreceiver device and it exhibits a twostate behaviour, akin to a switch. This is typical for quorumsensing systems: once concentrations of AIP (our signalling molecule) pass a certain threshold value, the agrreceiver will jump from its previously "off"state (negligible levels of PoPS from the P2 promoter) into an "on"state (significant levels of PoPS). In nature, once the bacterial density (and thus AIP concentration) surpasses a critical value, the bacteria will e.g. sporulate, become virulent or fluoresce  behaviour that is biologically expensive (however, critical) and only advantageous at high cell densities. The equations of this model read as follows:
A rigorous justification (i.e. all the steps leading to this set of equations) should probably be on here as well. SenderInteresting note: There is experimental evidence that AIP production rates are not depended on the basal expression level of agrB and agrD. AHLsignalling systemAccording to our experiments, B.subtilis does not degrade AHL. Model behaviour, look at Bangalore 2007 & Canton et al. (2008), assay luxsystem in B.subtilis. A Repressible promoterActivatorinhibitor systemLocal activation, lateral inhibition. Beyond TuringPattern formation with a single signalling molecule only?

<html>
<script language="JAVASCRIPT">
Today = new Date();
Jamboree = new Date("November 8, 2008");
msInADay = 1000 * 60 * 60 * 24;
display = Math.floor((Jamboree.getTime()  Today.getTime())/msInADay);
document.write("There are " + display +" days left until the Jamboree!");
</script>
</html>