Turing-Mechanismus

Der Turing-Mechanismus ist ein von dem britischen Mathematiker Alan Turing, einem der einflussreichsten Theoretiker der frühen Informatik, beschriebener Mechanismus, wie Reaktions-Diffusions-Systeme spontan Strukturen bilden können. Dieser Prozess steht auch noch heute im Mittelpunkt vieler chemisch-biologischer Strukturbildungstheorien, er erklärt zum Beispiel die Morphogenese von farbigen Mustern auf dem Fell von Tieren wie Zebra, Giraffe oder Kudu.

Turing hatte sich von 1952 bis zu seinem frühen Tod im Jahre 1954 mit Problemen der Theoretischen Biologie beschäftigt. In seiner 1952 veröffentlichten Arbeit zum Thema The Chemical Basis of Morphogenesis wurde dieser heute als Turing-Mechanismus bekannte Prozess erstmals beschrieben^[1]. Spätere Arbeiten Turings, unter anderem über die Bedeutung der Fibonacci-Zahlen für die morphologische Struktur von Pflanzen, blieben unveröffentlicht. Wegen Turings früherer Mitarbeit an nachrichtendienstlichen Projekten, wie der Dechiffrierung des Enigma-Codes, wurden seine gesammelten Werke erst 1992 zur Veröffentlichung freigegeben.

Turing-Modell für zwei Chemikalien

Ein Turing-Modell für zwei Chemikalien in dimensionsloser Form ist zum Beispiel gegeben durch das System von partiellen Differenzialgleichungen

$\left. \begin{matrix} u_{t} &=& \Delta u+\gamma f(u,v) \\ v_{t} &=& d\Delta v+\gamma g(u,v) \end{matrix} \right\} \mbox{ auf } B \subseteq \mathbb{R}^n$

mit Neumann-Randbedingungen $n\cdot\nabla u=n\cdot\nabla v=0 \mbox{ auf } \partial B$ und den Anfangsdaten u(.,t = 0) = u₀(.), v(.,t = 0) = v₀(.).

Der Vektor n ist der äußere Einheitsnormalvektor an $\partial B$, die Konstante d beschreibt das Verhältnis der Diffusionskoeffizienten der zwei Substanzen, u ist die Konzentration der aktivierenden Substanz (Aktivator) und v ist die Konzentration der deaktivierenden Substanz (Deaktivator oder Inhibitor). Die Konstante γ kann als Größe des Gebiets B (eine gewisse Potenz davon, abhängig von der Dimension n) oder als relative Stärke der Reaktionsterme im Vergleich zu den Diffusionseffekten interpretiert werden.

Die zentrale Idee der Turing-Instabilität besteht darin, ein räumlich homogenes, linear stabiles System zu betrachten, das im inhomogenen Fall, wenn sich Diffusion beschrieben durch geeignet gewählte Diffusionskoeffizienten auswirkt (d. h. $d \neq 1$), instabil ist („diffusionsgetriebene Instabilität“). Dieses Konzept war neuartig, da Diffusion auf dem Gebiet der partiellen Differenzialgleichungen im allgemeinen als stabilisierender Faktor aufgefasst wird.

Es basiert darauf, dass der Inhibitor schneller diffundiert als der Aktivator. Zuerst einmal entsteht dort, wo sich viel Aktivatorsubstanz befindet, auch viel Inhibitor. Dies führt allerdings nicht zum Verschwinden der Aktivatorsubstanz, da sich der Inhibitor aufgrund der schnellen Diffusion rasch verflüchtigt. Vor allem in Modellkonfigurationen auf beschränktem Gebiet mit verschwindender Neumann-Randbedingung sammeln sich vielmehr vergleichsweise hohe Inhibitorkonzentrationen an solchen abgelegenen Punkten des Gebiets an, wo sich wenig Aktivator befindet. Dort können sie nämlich die Zunahme der Aktivatorsubstanz erfolgreich unterbinden. An diesem Punkt zeigt sich, dass Modelle auf unbeschränktem Gebiet ein qualitativ anderes Verhalten zeigen, da die Inhibitorsubstanz tendenziell nach Unendlich abdiffundiert.

Für die spezielle Wahl B = [0,L] mit L > 0 bestimmen wir den Turing-Raum, d. h. die Menge an Parameterwerten, für die wir Turing-Instabilität beobachten können.

Lineare Stabilität des homogenen Systems

Sei (u₀,v₀) ein stationärer Zustand, d. h. f(u₀,v₀) = g(u₀,v₀) = 0, dann ist das obige System in Abwesenheit von Diffusionseffekten linear stabil an diesem stationären Zustand, wenn gilt

f_u + g_v < 0 und f_ug_v − f_vg_u > 0,

wobei wir hier und im folgenden die partiellen Ableitungen von f and g am stationären Zustand auswerten.

Instabilität räumlicher Störungen

Es sei W_k(r) eine Lösung des Eigenwertproblems ΔW_k + k²W = 0, d. h. $W_{k}(r)\propto\cos\left(\frac{n\pi}{L}r\right)$, wobei $n\in\mathbb{Z}$ und $k:=\frac{n\pi}{L}$ ist die sog. Wellenzahl. Mit dem Ansatz

$w(r,t)=\sum_{k}c_{k}\exp(\lambda(k)t)W_{k}(r)\;$

zeigt sich, dass das System linear instabil ist, falls $\Re\lambda(k)>0$ für ein $k\neq0$ gilt. Der Ausdruck λ(k) wird Dispersionsrelation genannt und nimmt nur dann positive Werte an, wenn gilt

df_u + g_v > 0 und $(df_{u}+g_{v})^{2}-4d(f_{u}g_{v}-f_{v}g_{u})>0\;.$

Aus der ersten Ungleichung ergibt sich insbesondere d > 1. Abhängig von der konkreten Wahl der Parameter hat das folgende Intervall von Wellenzahlen positiven Realteil,

$\gamma\frac{(df_{u}+g_{v})-\sqrt{(df_{u}+g_{v})^{2}-4d\det(A)}}{2d}<k^{2}<\gamma \frac{(df_{u}+g_{v})+\sqrt{(df_{u}+g_{v})^{2}-4d\det(A)}}{2d}\;,$

wobei $A:=\begin{pmatrix} f_u & f_v \\ g_u & g_v \end{pmatrix}$. Man beachte, dass $\gamma\propto L^{2}$. Das oben genannte Intervall heißt instabiles Intervall. Die diesen Wellenzahlen entsprechenden Wellenlängen nehmen im Laufe der Zeitentwicklung an Amplitude zu, während andere Wellenlängen gedämpft werden. Diese instabilen Mode beschreiben die verstärkten Muster. Da k nur diskrete Werte annimmt, gibt es nur eine endliche Zahl an verstärkten Wellenlängen.

Beispiele für geeignete Reaktionsterme sind

$f(u,v)=a-bu+\frac{u^{2}}{v}$ und $g(u,v)=u^{2}-v\;$ (Gierer, Meinhardt; 1972),

sowie

f(u,v) = a − u − h(u,v) und g(u,v) = α(b − v) − h(u,v), wobei $h(u,v):=\frac{\rho uv}{1+u+Ku^{2}}$, (Thomas, 1975).

Die Konstanten a, b, α, ρ und K sind positive Parameter, die geeignet gewählt werden müssen, damit das System die obigen Voraussetzungen erfüllt.

Musterbildung

Unbeschränkte Gebiete B entsprechen Modellen, die für Situationen relevant sind, in denen der Embryo weit größer ist als die Größenordnung der auszubildenden Muster und daher der Rand des Gebiets nichts zur Bevorzugung bestimmter Wellenlängen beitragen kann. Die Analysis ist in diesem Fall etwas einfacher. Im allgemeinen gibt es keine endliche Zahl verstärkter Wellenlängen, sondern eine gewisse Wellenzahl, die den größten Eigenwert aufweist und deren Muster letztlich ausgebildet wird.

Falls das Gebiet B größer wird im Laufe der Zeit, zum Beispiel dann, wenn der Embryo wächst, dann nimmt der Wert γ zu, und an gewissen Bifurkationspunkten werden aus verstärkten Moden gedämpfte, d. h., sie fallen aus dem instabilen Intervall, oder höhere Wellenzahlen, die bisher noch stabil waren, werden instabil. Dieser Prozess wird Moden-Selektion genannt und erklärt die komplexe Entwicklung von Mustern während der Morphogenese.

Literatur

• J. D. Murray: Mathematical Biology. Band 2: Spatial Models and Biomedical Applications. 3rd edition. Springer, New York NY u. a. 2003, ISBN 0-387-95228-4 (Interdisciplinary applied mathematics 18).

Einzelnachweise

1. ↑ Alan Turing: The chemical basis of morphogenesis (PDF). Phil. Trans. R. Soc. London B 237 pp 37-72 (1952). Originalartikel.