Системи реакції-дифузії: як леопард отримує свої плями

Забута стаття Алана Тюрінга

Алана Тюрінга пам'ятають насамперед за його роботу над обчисленнями, криптографією та теоретичними основами інформатики. Але в 1952 році, наприкінці свого життя, він опублікував статтю, яка виявилася настільки ж перетворювальною в зовсім іншій галузі: біології.

«Хімічні основи морфогенезу» пропонували математичний механізм, за яким запліднена яйцеклітина — одна хімічно однорідна клітина — могла б розвинутися в організм зі складною, диференційованою структурою. Тюрінг показав, що якщо дві хімічні речовини з різними швидкостями дифузії реагують одна з одною певним чином, просторові візерунки виникають спонтанно з початково однорідного стану.

Стаття залишалась переважно проігнорованою впродовж десятиліть. Молекулярна біологія тоді перебувала в зародковому стані, а ідея пояснювати біологічні візерунки диференціальними рівняннями здавалась надто абстрактною. До 1970-х і 1980-х років, коли математику розробили детальніше, а експериментальні докази накопичились, механізм Тюрінга визнали однією з найважливіших ідей у біології розвитку.

Принцип активатор-інгібітор

У серці прозріння Тюрінга лежить модель активатор-інгібітор. Уявіть два хімічні види:

• Активатор сприяє власному виробленню (позитивний зворотний зв'язок), а також стимулює вироблення інгібітора.
• Інгібітор пригнічує вироблення активатора (негативний зворотний зв'язок).

Без дифузії це звичайний контур негативного зворотного зв'язку — концентрації осцилюють або встановлюються на однорідній рівновазі. Дивовижний ефект, який відкрив Тюрінг, полягає в тому, що відбувається, коли інгібітор дифундує швидше за активатор.

У невеликій області, де концентрація активатора випадково трохи вища за середню, він стимулює ще більше вироблення активатора — позитивний зворотний зв'язок піднімає локальну концентрацію вгору. Водночас він виробляє інгібітор, але інгібітор швидко розсіюється дифузією, перш ніж встигне пригнітити локальний пік активатора. Концентрація активатора зростає, тоді як концентрація інгібітора зростає повільніше поблизу, але пригнічує вироблення активатора в навколишній області.

Результат: самопідсилюваний пік активатора, оточений пригніченою областю. Помножте це на всю поверхню — і ви отримаєте регулярний візерунок піків і западин — саме такий тип структури, що спостерігається в пігментації тварин, візерунках мушель та ембріональному розвитку.

Рівняння

Загальна система реакції-дифузії активатор-інгібітор описується двома пов'язаними диференціальними рівняннями в частинних похідних:

∂u/∂t = Dᵤ ∇²u + f(u, v)
∂v/∂t = D_v ∇²v + g(u, v)

Де u — концентрація активатора, v — концентрація інгібітора, Dᵤ та D_v — їхні коефіцієнти дифузії (D_v > Dᵤ необхідне для нестійкості Тюрінга), ∇² — оператор Лапласа (що вимірює локальну кривизну концентрації), а f і g — кінетика реакції — конкретні хімічні взаємодії між u та v.

Лапласіанний член реалізує дифузію: хімічні речовини перетікають з областей високої концентрації в області низької концентрації, згладжуючи відмінності. Реакційні члени f та g створюють або знищують хімічні речовини відповідно до локальних концентрацій.

Модель Ґрея-Скотта

Існує багато конкретних моделей реакції-дифузії, але модель Ґрея-Скотта є найбільш візуально багатою й широко дослідженою. Запроваджена П. Ґреєм та С.К. Скоттом у 1984 році, вона моделює автокаталітичну хімічну реакцію:

∂u/∂t = Dᵤ ∇²u − uv² + F(1 − u)
∂v/∂t = D_v ∇²v + uv² − (F + k)v

Тут u — одна хімічна речовина (реагент), v — інша (продукт). Член uv² представляє автокаталітичну реакцію — v каталізує власне вироблення за рахунок u. F — швидкість подачі (швидкість, з якою поповнюється свіжий u), а k — швидкість знищення (швидкість, з якою видаляється v).

Ці два параметри, F та k, є основними регуляторами системи Ґрея-Скотта. У просторі параметрів різні значення (F, k) породжують разюче різні візерунки:

• Плями — ізольовані піки v, що нагадують плями леопарда чи забарвлення гепарда.
• Смуги — видовжені смуги, як у зебри чи деяких видів риб.
• Лабіринтоподібні візерунки — звивисті, лабіринтоподібні канали, поширені в деяких видів коралів.
• Дірки — двоїстість плям: області високого u, оточені фоном v.
• Біжучі хвилі — самопоширювані фронти хімічної активності.
• Солітони — стабільні самовідтворювані плями, що діляться і поширюються.

Багатство простору поведінки моделі Ґрея-Скотта вражає. Дослідник Роберт Мунафо каталогізував десятки різних режимів візерунків, кожен з унікальним якісним характером, і все це — з тих самих двох рівнянь.

Чому ці візерунки з'являються в природі

Механізм Тюрінга тепер підтверджено в численних біологічних системах. Найпереконливіші докази походять із прямого спостереження хімічної динаміки, що породжує візерунки, а не лише кінцевого результату.

Пігментація риб вивчена особливо ретельно. У 2012 році дослідники з університету Осаки безпосередньо спостерігали динаміку активатор-інгібітор у шкірі данио-рерио під час формування смуг. Клітини-активатори (жовті) та клітини-інгібітори (чорні) демонстрували точнісінько ті просторові взаємозв'язки, які передбачали рівняння Тюрінга. Коли співвідношення дифузії змінили генетичною маніпуляцією, смуги перетворилися на плями — точно так, як передбачала математика.

Візерунки на мушлях зростають уздовж краю росту мушлі. У міру росту мушлі пігментоутворювальні клітини вздовж краю формують одновимірну систему реакції-дифузії. Візерунок «друкується» рядок за рядком у міру зростання мушлі, породжуючи математично передбачувані візерунки, що спостерігаються у Conus та інших видів.

Формування пальців в ембріонах — як пальці рук і ніг розвиваються як окремі структури — тепер вважається пов'язаним із динамікою реакції-дифузії в розвитку кінцівкового зачатка. Кількість і відстань між пальцями корелює з довжиною хвилі візерунка Тюрінга, яка залежить від фізичного розміру тканини, що розвивається.

Шкірні гребені, що формують відбитки пальців, також демонструють характеристики, узгоджені з візерункоутворенням реакції-дифузії, хоча конкретні молекулярні компоненти досі дискутуються.

Чисельна реалізація

Модель Ґрея-Скотта чисельно проста в реалізації. Ключові кроки на кожному кроці часу такі:

1. Обчислити лапласіан u та v, використовуючи скінченні різниці або згортку з дискретним лапласіанним ядром.
2. Застосувати реакційні члени: обчислити зміну u та v через автокаталітичну реакцію та члени подачі/знищення.
3. Оновити концентрації u та v, додавши внески дифузії та реакції.
4. Обмежити концентрації допустимим діапазоном (зазвичай від 0 до 1).

for each cell (x, y):
    laplacian_u = discrete_laplacian(u, x, y)
    laplacian_v = discrete_laplacian(v, x, y)

    uvv = u[x,y] * v[x,y] * v[x,y]   // автокаталітичний член

    du = Du * laplacian_u - uvv + F * (1 - u[x,y])
    dv = Dv * laplacian_v + uvv - (F + k) * v[x,y]

    u_new[x,y] = u[x,y] + du * dt
    v_new[x,y] = v[x,y] + dv * dt

Дискретний лапласіан зазвичай обчислюється за допомогою згортки з ядром 3×3 чи 5×5, що усереднює різницю концентрацій між кожною клітиною та її сусідами. Симуляція виконує це оновлення тисячі разів на секунду з використанням обчислень на GPU, що робить можливою взаємодію в реальному часі.

Застосування в мистецтві й дизайні

Окрім біології, системи реакції-дифузії стали потужним інструментом для генеративного мистецтва й процедурного дизайну. Візерунки, які вони породжують, мають якість, що є водночас математичною й органічною — точною, але не механічною.

Дизайнери використовують модель Ґрея-Скотта та споріднені моделі для генерації:

• Текстур для 3D-рендерингу — шкіри тварин, позаземних поверхонь, абстрактних фонів.
• Ювелірних візерунків — плинна, біологічна якість RD-візерунків підходить для декоративного металевиробництва й кераміки.
• Архітектури — параметричних фасадів, натхненних природним візерункоутворенням.
• Візерунків тканини — дизайнів для ткацтва та друку, що відлунюють природну пігментацію.

Ця система також використовується в математичних мистецьких рухах, таких як алгоритмічне мистецтво й генеративне мистецтво, де художник визначає правила й параметри замість малювання окремих елементів.

Відкриті питання

Попри десятиліття досліджень, первісне бачення Тюрінга повної математичної теорії морфогенезу залишається нездійсненим. Багато питань лишаються відкритими:

• Якою є конкретна молекулярна ідентичність активатора та інгібітора в більшості систем візерункоутворення? У багатьох випадках ми можемо спостерігати, що діє RD-динаміка, але точні хімічні речовини невідомі.
• Як візерунки реакції-дифузії координуються з механічними процесами? Клітини не просто змінюють колір — вони рухаються, діляться і змінюють форму. Зв'язок між хімічним візерункоутворенням і механічними силами лише починає вивчатися.
• Як візерунки масштабуються відносно розміру тіла? Ембріон миші та ембріон слона мають різні розміри тканин, але схожі довжини хвиль візерунків відносно тіла. Механізм масштабування залишається незрозумілим.

Стаття Тюрінга 1952 року дала нам математичну рамку. Сімдесят років по тому біологія все ще розкриває глибину того, що він побачив.