🧬 Біологія · Обчислювальна біологія
📅 Березень 2026⏱ 12 хв🟡 Середній
·
Останнє оновлення: 28 травня 2026 р.
Системна біологія: мережі, схеми та емерджентне життя
Клітина містить ~25 000 генів, ~100 000 різних білків і мільйони метаболітів, що взаємодіють через десятки тисяч хімічних реакцій одночасно. Системна біологія застосовує інженерну інтуїцію та математичне моделювання, щоб зрозуміти не окремі молекули, а колективну поведінку цих мереж — прагнучи пояснити, як життя постає з хімії.
1. Типи біологічних мереж
Різні типи молекулярних взаємодій утворюють різні класи мереж, кожен із характерною топологією та функцією:
- Генні регуляторні мережі (GRN): Орієнтований граф, у якому фактори транскрипції (TF) зв'язуються з промоторами та активують або пригнічують експресію цільових генів. Вузли = гени, ребра = регуляторні взаємодії (+/−). GRN людини: ~1500 TF, що регулюють ~25 000 мішеней. Ключова властивість: безмасштабна топологія (більшість генів мають мало зв'язків, кілька «вузлових» TF — тисячі).
- Мережі білок-білкових взаємодій (PPI): Неорієнтований граф фізичного зв'язування білків. PPI-мережа людини: ~300 000 взаємодій (оцінка, дуже неповна). Вузли часто є життєво важливими генами — вимкнення їх → загибель клітини.
- Метаболічні мережі: Орієнтований дводольний граф (метаболіти ↔ реакції). Стехіометрична матриця S: S·v = 0 у стаціонарному стані (S = стехіометрія, v = вектор потоків). Аналіз балансу потоків (FBA) оптимізує цільову функцію (швидкість росту) за обмежень на потоки.
- Сигнальні мережі: Білки передають сигнали від поверхні клітини до ядра через посттрансляційні модифікації (фосфорилювання, убіквітинування). Динамічні, швидкі (від секунд до хвилин).
2. Мережеві мотиви: будівельні блоки схем
Алон (2007) — мережеві мотиви: патерни, що трапляються набагато частіше
в біологічних мережах, ніж у випадкових мережах із тим самим розподілом
ступенів. Ключові мотиви: 1. Негативна саморегуляція
(авторепресія): TF пригнічує власну транскрипцію. ~50% TF E. coli
авторепресовані. Функція: пришвидшує час відгуку, зменшує шум
експресії (дисперсію). Лабораторія Алона: концентрація білка-відповіді зростає та
менше «перестрибує» за авторепресії → швидша й точніша адаптація.
2. Петля прямого зв'язку (FFL) — 8 підтипів, 2 когерентні / 2 некогерентні основні
типи: X активує Y, X активує Z, і Y активує Z: когерентна FFL
типу 1 (більшість у генних мережах E. coli): Z вмикається лише тоді, коли
присутні ОБИДВА — X і Y. Функція: знакочутлива затримка — Z реагує
на стійкий сигнал X, але ігнорує короткі імпульси X (бо Y має
повільну динаміку). ~ фільтр низьких частот / фільтр шуму. Некогерентна
FFL типу 1: X активує Z, але Y пригнічує Z. X також активує Y. Функція:
генератор імпульсів — Z вмикається, коли активується X, потім вимикається, коли Y
накопичується. Перехідний відгук. 3. Позитивний зворотний зв'язок / бістабільний
перемикач: X активує Y, Y активує X. Два стабільні стани: ОБИДВА вимкнені або
ОБИДВА увімкнені. Бістабільність → ухвалення рішень, фіксація долі клітини. Приклад:
Mos/MAP-кіназа при дозріванні ооцита Xenopus.
3. ОДУ-моделі експресії генів
Проста модель експресії гена: dm/dt = α_m − β_m · m (продукування мРНК −
деградація) dp/dt = α_p · m − β_p · p (синтез білка −
деградація) де m = концентрація мРНК, p = концентрація білка
α_m = швидкість транскрипції (з активатором / без нього) β_m = швидкість
деградації мРНК (t½ мРНК ~3-10 min у E. coli; ~30-60 min у ссавців)
α_p = швидкість трансляції β_p = швидкість деградації білка (t½ білка
~години-дні; у E. coli ~20 годин у середньому) Стаціонарний стан: m_ss =
α_m/β_m, p_ss = α_p · α_m / (β_m · β_p) Час відгуку (на стрибкоподібну зміну
транскрипції): t_response ≈ ln(2) · max(1/β_m, 1/β_p) Визначається
найповільнішою швидкістю деградації (зазвичай білка). Функція Гілла
(активація транскрипції): f(X) = X^n / (K_d^n + X^n) (активатор)
f(X) = 1 / (1 + (X/K_d)^n) (репресор) K_d = константа дисоціації
(концентрація при напівмаксимальній активації) n = коефіцієнт Гілла
(кооперативність; сигмоїдальний для n > 1) n = 1: Міхаеліс-Ментен (без
кооперативності), градуйований відгук n = 2-4: перемикальний відгук; n = ∞:
ідеальний бінарний перемикач Бістабільність потребує: n ≥ 2 у петлі
позитивного зворотного зв'язку (приблизно)
4. Сигнальні каскади: шлях MAPK
Каскад MAPK (мітоген-активована протеїнкіназа) — це консервований сигнальний модуль, що контролює проліферацію, диференціювання клітин і стресові відповіді. Його архітектура у клітинах ссавців:
- Рівень 1 — Ras (ГТФаза): Активована рецепторна тирозин- кіназа (напр., EGFR) залучає GEF (Sos) → Ras-GTP. Інактивується білками GAP.
- Рівень 2 — Raf (MAPKKK): Ras-GTP залучає та активує кіназу Raf. Raf фосфорилює MEK.
- Рівень 3 — MEK (MAPKK): Кіназа подвійної специфічності — фосфорилює ERK по залишках і Thr, і Tyr (незвична подвійна вимога).
- Рівень 4 — ERK (MAPK): Двічі фосфорильована ERK (ppERK) є активною. Мішені: фактори транскрипції (Elk-1, c-Fos), рибосомні ферменти (RSK), регулятори клітинного циклу (циклін D1).
Ультрачутливість і бістабільність: Трирівнева
архітектура каскаду в поєднанні з вимогою подвійного фосфорилювання
для активації ERK створює дуже ультрачутливий (крутий
вхід-вихід) відгук. Голдбетер і Кошланд (1981) показали, що
протилежно спрямовані цикли кінази/фосфатази можуть породжувати позірно кооперативну
поведінку (видимий коефіцієнт Гілла >1) без фактичної
кооперативності в окремих реакціях. Ця «ультрачутливість
нульового порядку» може бути поширеною в сигналінгу й пояснює, як клітини
ухвалюють перемикальні рішення.
5. Осцилятор p53-Mdm2
Осцилятор p53 («охоронець геному») після пошкодження ДНК: Структура
зворотного зв'язку: p53 → активує → Mdm2 (транскрипційно) Mdm2 → пригнічує
→ p53 (убіквітин-опосередкована деградація) Пошкодження ДНК → порушує
опосередковану Mdm2 деградацію p53. ОДУ-модель (спрощена, Geva-Zatorsky
та ін., 2006): dp53/dt = α_p53 − β_p53 · p53 · Mdm2 ... (продукування −
опосередкована Mdm2 деградація) dMdm2/dt = α_Mdm2 · p53 − β_Mdm2 · Mdm2
(транскрипція через p53 − деградація) [Плюс член затримки для ядерного
транспорту Mdm2 ~30-40 min] Результат: імпульси p53 та Mdm2 кожні ~5-6 годин
після пошкодження ДНК. Спостерігалося у: одноклітинній прижиттєвій візуалізації флуоресцентного
p53 у клітинах MCF7. Кожен імпульс відповідає одній «спробі» оцінити
пошкодження та запустити апоптоз. Цифровий відгук p53: Кількість імпульсів
кодує тяжкість пошкодження. Незначне пошкодження: 1-2 імпульси → репарація ДНК,
виживання. Тяжке пошкодження: багато імпульсів → рішення про апоптоз. Значення для
раку: p53 мутований у ~50% усіх ракових пухлин людини (найпоширеніша
мутація в раку). Втрата p53 → клітини оминають апоптоз після пошкодження
ДНК → неконтрольований ріст. MDM2 надекспресований у ~10% ракових пухлин →
надмірна деградація p53.
6. Булеві мережі та стани-атрактори
Коли кінетичні параметри невідомі, булеві мережі пропонують грубішу, але обчислювально доступну модель. Кожен ген перебуває або в стані ON (1), або OFF (0), із логічними правилами оновлення:
Булева мережа: Стан: вектор (x₁, x₂, ..., x_N) ∈ {0,1}^N Оновлення:
x_i(t+1) = f_i(x₁(t), ..., x_N(t)) — логічна функція 2^N можливих
станів. За послідовних або синхронних оновлень траєкторії у просторі
станів детерміновані → сходяться до атракторів: Атрактор-нерухома
точка: один стан (стабільний тип клітини / онтогенетична доля) Граничний
цикл: коливальна траєкторія (клітинний цикл, циркадний осцилятор)
NK-модель Кауфмана (1969): N генів, кожен із K входами (призначеними
випадково). K=2: упорядкований режим — небагато атракторів малої довжини
(стійкі до збурень) K>2: хаотичний режим — експоненційна кількість
довгих атракторів (крихкі) Реальні GRN: K ~ 2 (розріджена зв'язність) →
стійка регуляція генів. Приклад — мережа сегментної полярності Drosophila
(Albert & Othmer, 2003): 15-генна булева мережа відтворює правильні
патерни експресії генів для всіх сегментів ембріона мухи. 7 атракторів
відповідають 7 спостережуваним типам клітин. Стійка до >90% одногенних
збурень.
7. Застосування синтетичної біології
Синтетична біологія застосовує інженерні принципи проєктування — стандартизовані компоненти, модульний дизайн, абстракцію — для побудови нових біологічних схем:
- Тригерний перемикач (Gardner та ін., 2000): Два репресори, кожен пригнічує промотор іншого. Два стабільні стани: R₁-on/R₂-off ТА R₁-off/R₂-on. Перемикання шляхом тимчасової індукції. Перша сконструйована бістабільна схема. Реалізована в E. coli.
- Репресилятор (Elowitz & Leibler, 2000): Три гени-репресори (lacI, tetR, cI) у циклі — кожен пригнічує наступний. Результат: коливання рівнів білків із періодом ~150 хвилин у E. coli. Перший сконструйований генетичний осцилятор. Продемонстрував, що періодичну динаміку можна спроєктувати синтетично.
- Біосенсори: Клітини, сконструйовані для виявлення забруднювачів, метаболітів чи патогенів. Фактор транскрипції, активований цільовою речовиною → запускає репортерний ген. Застосування: виявлення забруднення води, діагностика.
- Метаболічна інженерія: Артемізинін (протималярійний препарат, спершу з рослини A. annua) виробляється у сконструйованих S. cerevisiae у промисловому масштабі (лабораторія Кіслінга / Amyris). Фарнезил- пірофосфат → артемізинінова кислота через 10 гетерологічних генів. Знижує вартість виробництва вдесятеро.
- Інженерія CAR-T-клітин: Т-клітини, сконструйовані з химерними антигенними рецепторами — синтетичними рецепторними білками, що зв'язуються з пухлинними антигенами → активують знищення Т-клітинами. Схвалені FDA для певних видів раку крові. CAR із логічними вентилями (вентилі AND, NOT із двома типами рецепторів) розробляються для солідних пухлин.