Циркадний ритм: Мережі генів годинника – Моделювання клітинного осцилятора

Основний механізм транскрипції-перекладу

Мелатонічний годинник у ссавців базується на механізмі транскрипції-перекладу (ТПЗ). Фактори транскрипції CLOCK та BMAL1 утворюють гетеродимер, який стимулює експресію генів Per (Period) та Cry (Cryptochrome). Білки Per та Cry накопичуються в цитоплазмі, димеризуються, повертаються в ядро і безпосередньо пригнічують комплекс CLOCK:BMAL1 – негативний зворотний зв’язок, побудований повністю на рівні генного вираження, а не швидкої ферментативної кінетики.

CLOCK:BMAL1 → активує транскрипцію Per, Cry Per, Cry (білок) накопичуються → димеризуються → повертаються в ядро Per:Cry → пригнічує активність CLOCK:BMAL1 → зменшується транскрипція Per, Cry деградація (за допомогою фосфорилювання CK1δ/ε-медіатором → протеасом) придушення знято → цикл відновлюється – повний період циклу ≈ 24 години у ссавців Другий, переплетений цикл (ROR/REV-ERB, що діють на Bmal1 безпосередньо) додає стійкості, але цикл PER/CRY-CLOCK/BMAL1 є мінімальним основним механізмом, який математична модель повинна відтворити 24-годинний ритм.

CLOCK:BMAL1 → activates Per, Cry transcription
PER, CRY (protein) accumulate → dimerise → re-enter nucleus
PER:CRY → represses CLOCK:BMAL1 activity → transcription falls
PER, CRY degrade (via CK1δ/ε-mediated phosphorylation → proteasome)
repression lifts → cycle restarts     — full cycle period ≈ 24 h in mammals

Гівінський осцилятор: мінімальна зворотна модель

Модель Бріана Гівінса 1965 року спрощує будь-який такий цикл до трьох змінних: мРНК (X), проміжний білок (Y) і фінальний білковий репресор (Z), при цьому Z інактивує власну транскрипцію мРНК (X):

dX/dt = v₁ / (1 + (Z/K)^n) − d₁X // мРНК (напр., Per/Cry) dY/dt = v₂X − d₂Y // цитоплазматичний білок dZ/dt = v₃Y − d₃Z // ядерний репресор (комплекс PER:CRY) Термін інактивування є Hill-функцією з коефіцієнтом кооперативності n. Гівінс показав, що підтримання осциляції в цій системі з трьох змінних ОДУ вимагає n > 8 — надмірно високого рівня молекулярної кооперативності, оскільки реалізація димеризації та інактивування PER/CRY комплексу відбувається набагато більш м'яким нелінійністю.

dX/dt = v₁ / (1 + (Z/K)^n) − d₁X    // mRNA (e.g. Per/Cry)
dY/dt = v₂X − d₂Y                    // cytoplasmic protein
dZ/dt = v₃Y − d₃Z                    // nuclear repressor (PER:CRY complex)

Чому затримка має значення: DDE vs. ODE

Реальний молекулярний годинник досягає надійної 24-годинної осциляції з меншою кооперативністю, ніж вимагає модель Goodwin's ODE, оскільки біологія вводить явне часове запізнення: транскрипція, ядерна експортуваність, переклад, фолдинг білка, димерування та імпорт PER:CRY займають кілька годин — не миттєво, як це припускається в ODE. Заміна миттєвого зворотного зв’язку на запізнення відмінне рівняння (DDE) дозволяє це врахувати безпосередньо:

dX/dt = v₁ / (1 + (Z(t − τ)/K)^n) − d₁X // пригнічення дії X за допомогою значення Z's τ годин у минулому Завдяки затримці τ декількох годин, стабільні осциляції з’являються навіть при n = 1 (простий, не кооперативний Michaelis-Menten механізм пригнічення) — що відповідає помірній кооперативності, яка спостерігається при пригніченні CLOCK:BMAL1 PER:CRY. Затримка, а не надчутливість, є ключовою рисою, яка робить біологічний годинник стійким: явне транскрипційне та перекладачне запізнення замінює кооперативність у генеруванні обмеженого циклу.

dX/dt = v₁ / (1 + (Z(t − τ)/K)^n) − d₁X
// repression acts on X using Z's value τ hours in the past

Обмежені цикли та біфуркація Хопфа

Зі збільшенням затримки τ (або, еквівалентно, коефіцієнта Хілла n) за критичного значення, стабільний нерухомий пункт у системі втрачає стабільність внаслідок біфуркації Хопфа, і народжується стабільний обмежений цикл — замклена, самопідтримувана орбіта в фазовому просторі (X, Y, Z), в яку система орієнтується незалежно від невеликих збурень від свого початкового стану. Під порогом нерухомий пункт є стабільним, а рівні mRNA/білків встановлюються на константи без ритму; під час біфуркації нерухомий пункт втрачає стабільність; за його межами з'являється самопідтримувана ~24-годинна осциляція, амплітуда якої визначається нелінійністю, а не початковими умовами; і щоденний світловий вплив фазово зміщує цей обмежений цикл, щоб він збігався з 24-годинним сонячним днем (ентрейнмент).

Ця незалежність від початкових умов є математичною ознакою міцного біологічного годинника: клітина, що починається з будь-якого стартового рівня PER/CRY, сходиться до однієї ритмічної орбіти, тому ізольовані фібробласти всі встановлюються в порівняно з нею ~24-годинними циклами, незважаючи на те, що вони починаються з різних транскрипційних станів.

Зв’язні клітини: від одного осцилятора до ритму тканини

Одиночний годинник клітини сам по собі є шумним осцилятором — його період коливається приблизно на ±1–2 години від клітини до клітини через генетичну експресію. У супрахиазматичному ядрі, близько 20 000 таких шумних одиночних осциляторів зв’язані між собою за допомогою нейропептидного сигналізу (переважно VIP, що діє на рецептори VPAC2), синхронізуючи їх фази в єдиний, точний тканини рівень ритму — такий же принцип зв’язного осцилятора використовується для зв’язаних підвісних систем або моделі Курамото.

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

Часті запитання

Спробуйте наживо

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