Rozprzestrzenianie się Turina: wzory reakcji-diffuzyjne i morfogeneza

Turingeva 1952 roku propozycja morfogenezy

Artykuł Alan Turina z 1952 roku „Chemical Basis of Morphogenesis” zabrzmiało na głębokim zagadnieniu: jak może niesymetryczny jądro płodne wytworzyć organizm o precyzyjnej strukturze przestrzennej? Turin postulował morfogene — chemiczne substancje rozpraszalne, które przenoszą informacje pozycyjne poprzez ich gradienty koncentracji. Udowodnił, że system dwóch reagujących i rozpraszających się morfogene może być nieprzystępny dla małych perturbacji przestrzennych, nawet jeśli każdy z indywidualnych procesów (rozpraszanie lub jednolite reakcje) byłby stabilny. Instability Turina jest counterintuitive: rozpraszanie normalnie jest procesem usmiglającym, jednak tutaj tworzy strukturę zamiast jej zniszczać. Artykuł został w większości ignorowany podczas życia Turina i dopiero w latach 70. XX wieku zostało rozpoznane jako fundamentalny.

Mechanizm aktywatora-ihibitora

Ogólny szablon łączy dwie koncentracje, aktywator u i inhibicor v, za pomocą równań reakcji-diffuzji. Krytycznym wymaganiem jest różnorodność diffuzji: inhibicor musi rozprzestrzenić się znacznie szybciej niż aktywator. Gdzie aktywator jest lokalnie podwyższony, on zwiększa się sam na siebie (aktywacja autonomiczna) podczas wysyłania szybszego inhibicora rozpraszającego się wzdłuż do ograniczenia aktywacji w okolicy — aktywacja lokalna, inhibicja laterala — co ustanawia charakterystyczną odległość między elementami wzoru. Model Gierer-Meinhardta (1972), oparty na regeneracji głowy/podeszwy hydry, dostarcza jasnych kinetyki za pomocą niejolinecznej proporcji u²/v, która saturuje wzór przy skończonej amplitudzie.

Model Gray-Scott i jego zoo parametrów

Model Gray-Scott (1984) symuluje reakcję autokatalityczną, w której U konwertuje się na V w obecności V, podczas gdy V wolno zdegradowuje. Dwa parametry — stopa dostawu F (zastępuje U) i stopa zabójcza k (usuwa V) — kontrolują niezwykle bogatą diagram fazowy: (F,k)=(0,035, 0,065) daje punkty, (0,037, 0,060) daje węże, (0,029, 0,057) daje spirale, (0,040, 0,060) daje labiryntowe pasy, a (0,025, 0,055) tworzy „mitozę” — lokalne punkty, które rozszerzają się i podzielają, symulując podział komórkowy.

∂u/∂t = Du·∇²u − u·v² + F·(1 − u)
∂v/∂t = Dv·∇²v + u·v² − (F + k)·v

F=0.035, k=0.065 → spots      F=0.025, k=0.055 → mitosis
F=0.037, k=0.060 → worms      F=0.029, k=0.057 → spirals

Potwierdzone w rzeczywistych zwierzętach

Murray (1988) pokazał, że rozwiązywanie równań reakcji-diffuzyjnych na obrazkach różnych kształtów powtarza wzory skórzanej pasiastej u lwa, pasiastych kikutów konia i przemianę długich pasiastych odcinków ogona na pętle. Najwyraźniejsze potwierdzenie molekularne pochodzi od Sick et al. (2006, Science), którzy zidentyfikowali WNT jako aktywator i DKK1 jako inhibicję przestrzeni kikutów włosów myszy — perturbacja DKK1 zmieniła przestrzeń dokładnie tak jak przewidziano. Nakamasu et al. (2009) później pokazali, że melanofory i xantofory ryby zebra spełniają tę samą regułę aktywacji lokalnej/inhibicji dalekiego zakresu w tworzeniu pasiastych. Tę sama matematyka rządzi również reakcją chemiczną Belousov-Zhabotinskiego, piaszczystymi zboczami i pasami roślinności w ekosystemach suchej stepy.

Często zadawane pytania

Wypróbuj na żywo

Wszystko powyżej działa bezpośrednio w Twojej przeglądarce — otwórz Turing Diffusion i zmieniaj parametry podczas działania. Nic nie jest instalowane ani przesyłane na serwer, cały model działa w jednej karcie.