Морфогенез — процес набуття організмом форми — є однією з центральних нерозв'язаних проблем біології розвитку. Цей симулятор поєднує два фундаментальні механізми: силову механіку клітин, де округлі клітини ростуть, діляться і пакуються в тканину через сили відштовхування та адгезії; та реакційно-дифузійне поле Ґрея–Скотта, що генерує просторово-градієнтне поле морфогенів у тканині. У своїй знаменитій роботі 1952 року Алан Тюрінґ передбачив, що дифундуючі хімічні «морфогени» з різною швидкістю дифузії можуть спонтанно порушити просторову симетрію, породжуючи плями, смуги та лабіринтові патерни на хутрі тварин, шкірі риб та зародках.
Керуйте швидкістю поділу клітин, швидкістю росту, силою адгезії та максимальною кількістю клітин. Увімкніть поле морфогену Ґрея–Скотта (з регульованими швидкістями живлення F і вимирання k) і спостерігайте, як воно «фарбує» кожну клітину за її локальною хімічною долею — перемикаючись між плямами і смугами при зміні F і k через межу нестійкості Тюрінґа.
Що таке модель Ґрея–Скотта і які патерни вона утворює?
Модель Ґрея–Скотта — двокомпонентна реакційно-дифузійна система: ∂U/∂t = D_U∇²U − UV² + F(1−U) та ∂V/∂t = D_V∇²V + UV² − (F+k)V. Залежно від F і k система продукує надзвичайно різноманітні патерни: стабільні плями, плями що діляться «як клітини», смуги, лабіринти та мандрівні хвилі. Ці патерни відповідають реальним біологічним структурам — від текстур коралових рифів до малюнків черепашок конуса.
Як працює поділ клітин у симуляції?
Кожна клітина росте, повільно збільшуючи свій радіус пропорційно до швидкості росту та коефіцієнту скупченості (сусідні клітини сповільнюють ріст). Коли клітина перевищує радіус поділу R_div = 9 пкс, вона з певною ймовірністю ділиться на дві дочірні клітини радіуса R_min = 4, розміщені вздовж випадкової осі. Це реалізує спрощений варіант мітотичного поділу клітин, вловлюючи ключову кінетику: ріст обмежений скупченістю, ймовірність поділу пропорційна слайдеру швидкості поділу.
Що контролює параметр адгезії?
Адгезія задає силу короткодіяльного притяжіння між сусідніми клітинами, що не перетинаються. Висока адгезія — щільні кластери; низька — клітини розходяться після поділу. Це наближення клітинно-клітинної адгезії через E-кадгерин, втрата якого в процесі епітеліально-мезенхімального переходу (ЕМП) дозволяє ракам відокремлюватися від первинної пухлини та метастазувати.
Патерни Тюрінґа виникають, коли активатор (що стимулює і себе, і інгібітор) дифундує повільніше за інгібітор — умова «дифузійно-рушійної нестійкості». Малі випадкові флуктуації концентрації активатора локально підсилюються, але глобально придушуються інгібітором, що швидше поширюється. Це стихійно обирає характерну просторову довжину хвилі. Підтвердженням у реальному світі є патерни меланоцитів у зебрафіша (маніпуляції, що перемикають смуги на плями), відстань між пальцями в кінцівках і розташування волосяних фолікулів.
Кругла тканинна межа в симуляції відповідає біологічно реалістичному обмеженню: у реальних зародках тканини розвиваються в просторових межах, встановлених епітеліальними мембранами та позаклітинним матриксом. Обмеження обмежує доступні довжини хвиль патерну кратними, що «вміщаються» в геометрію тканини, тому кількість і розташування плям частково визначаються геометрично. Ґаструлоїди — синтетичні зародки з стовбурових клітин — стабільно формують вісь перед–зад, демонструючи, що геометрія сама по собі керує морфогенезом.
Ступінь заповнення вимірює, яку частку площі тканини займають тіла клітин (π×r²/площа_тканини). Випадкова упаковка однакових кіл досягає максимуму ~0,906 (гексагональна щільна упаковка), але реальні біологічні тканини зазвичай мають ступінь заповнення 0,7–0,85. В епітеліальних моношарах клітини пакуються у полігони Вороного з характерним розподілом 5-, 6- і 7-кутників — порядок, що виникає з мінімізації механічної енергії.
Морфогени — сигнальні молекули, що утворюють просторові градієнти і вказують клітинам прийняти різні долі залежно від місцевої концентрації. Класичні морфогени: Bicoid (вісь голова–хвіст у дрозофіли), SHH (ідентичність пальців), BMP-4 (дорсо-вентральна вісь). Фактори росту (EGF, FGF, PDGF) переважно стимулюють проліферацію і виживання клітин, а не специфікацію долі. У симуляторі вид V системи Ґрея–Скотта виступає морфогеном, що визначає долю.
В реальних тканинах ріст клітин обмежується контактним гальмуванням — повністю оточені сусідами клітини перестають ділитися. Молекулярна основа — шлях Hippo: при механічній скупченості транскрипційні фактори YAP/TAZ фосфорилюються і інактивуються, вимикаючи гени росту. Ракові клітини часто втрачають контактне гальмування — це характерна ознака злоякісності. У симуляторі скупченість моделюється підрахунком сусідів і пропорційним зниженням швидкості росту.
Характерна довжина хвилі Тюрінґа λ ≈ 2π√(D_activator/degradation_rate) задає відстань між плямами або смугами. Біологічні системи підлаштовують цей масштаб під розміри тканини: при рості зародка однакові параметри морфогену дають пропорційно масштабований патерн. Нещодавні роботи показали, що деякі градієнти морфогенів масштабуються з розміром тканини через петлі зворотного зв'язку, що забезпечує стабільність формоутворення в широкому діапазоні розмірів зародка.
Так. Диференціальне поверхневе натягнення між популяціями клітин рухає сортування клітин (гіпотеза диференціальної адгезії Стайнберґа, 1963); складання тканин при гаструляції та нейруляції рухається апікальним звуженням клітин; стиснення в кишці що розвивається спричиняє її вигин у характерні петлі. Механотрансдукція — перетворення механічних сил у зміни експресії генів — пов'язує хімічний та механічний аспекти морфогенезу, а цитоскелет є ключовим трансдуктором сигналів.
Органоїди — тривимірні міні-органи, вирощені зі стовбурових клітин у гелевому матриксі, — спонтанно самоорганізуються у структури, що нагадують тканину-джерело, завдяки тим самим морфогенним градієнтам типу Тюрінґа, що моделюються тут. Кишкові органоїди мають архітектуру крипта–ворсинки, легеневі — альвеоли. Органоїди революціонізують тестування ліків, персоналізовану медицину (органоїди пухлини пацієнта для прогнозу чутливості до хіміотерапії) та моделювання хвороб.