Молекулярна клітинна біологія — це біофізика в найдетальнішому прояві: окремі молекули, що прикладають сили в кілька піконьютонів, дифундуючи й реагуючи в компартменті заледве мікрона завширшки, підкоряючись термодинамічним законам, що не терплять винятків. Одномолекулярні методи — оптичний пінцет, FRET, крио-ЕМ, атомно-силова мікроскопія — дозволили спостерігати роботу окремих молекулярних машин у реальному часі, розкриваючи механізми, що водночас елегантні й вражаючі.
1. Реплікація ДНК і реплісома
Подвійна спіраль ДНК зберігає генетичну інформацію в послідовності пар основ: аденін з'єднується з тиміном (2 водневі зв'язки), гуанін із цитозином (3 водневі зв'язки). Перш ніж клітина поділиться, весь геном має бути скопійований із надзвичайною точністю. У людини ~6 мільярдів пар основ реплікуються за ~8 годин скоординованим молекулярним комплексом, що зветься реплісомою, досягаючи рівня похибки приблизно 1 на 10^{10} після коректорського читання й репарації неспарованих основ.
Кінетика та точність реплікації ДНК
Стабільність пар основ Уотсона-Кріка (ΔG за 37°C, у воді):
Пара AT: ΔG ≈ −1 ккал/моль (2 водневі зв'язки)
Пара GC: ΔG ≈ −3 ккал/моль (3 водневі зв'язки)
Стекінг сусідніх пар додає ~−1 до −3 ккал/моль
Температура плавлення дуплексу ДНК T_m:
T_m = 81,5 + 16,6 log[Na+] + 0,41(%GC) − 675/n (правило Уоллеса, n = довжина в п.о.)
T_m олігонуклеотиду: зазвичай 55–65 °C для 20-мірних праймерів у ПЛР
Реплікаційний апарат (реплісома E. coli):
Геліказа (DnaB): розкручує дволанцюгову ДНК зі швидкістю ~1000 п.о./с, використовує гідроліз АТФ
Праймаза (DnaG): синтезує короткі РНК-праймери (~10 нт) для кожного фрагмента Оказакі
ДНК-полімераза III: процесивний синтез, ~1000 нт/с, коректорська екзонуклеаза 3'→5'
Лідируючий ланцюг (безперервний): ~1000 нт/с; відстаючий ланцюг (переривчастий): фрагменти Оказакі 1–2 т.п.о.
Лігаза (LigA): зшиває розриви після дозрівання фрагментів Оказакі
Точність реплікації:
Частота помилок вставки «голої» полімерази: ~10^{-5} на основу
Коректорська екзонуклеаза 3'→5': знижує до ~10^{-7}
Репарація неспарованих основ (MutS/MutL/MutH): знижує до ~10^{-10}
Геном людини 3×10^9 п.о.: ~0,3 помилки на поділ (корисно для еволюції!)
Реплікація в людини:
~30 000 точок початку реплікації на гаплоїдний геном
Середня відстань між точками початку: ~100 т.п.о.
Швидкість вилки: ~1 т.п.о./хв (проти 60 т.п.о./хв у бактерій)
Тривалість S-фази: ~8 год
2. Згортання білків і енергетичний ландшафт
Поліпептидний ланцюг із N амінокислот має ~3^N можливих конформацій основного ланцюга — для білка зі 100 залишків це більше конформацій, ніж атомів у видимому Всесвіті. Проте більшість білків спонтанно згортаються до унікального нативного стану за мікросекунди-секунди. Парадокс Левінталя (1969) вказував, що випадковий пошук неможливий; має існувати спрямований механізм згортання. Сучасна відповідь — енергетична воронка: конформаційний простір має форму нерівної воронки, а частково згорнуті структури мають нижчу вільну енергію, ніж повністю розгорнуті.
Термодинаміка згортання білків та AlphaFold
Вільна енергія згортання: ΔG_fold = ΔH_fold − TΔS_fold ΔH_fold: водневі зв'язки, ван-дер-ваальсові сили, дисульфідні зв'язки (− кілька ккал/моль кожен) ΔS_fold: втрата ентропії ланцюга (несприятливо) + гідрофобний ефект (сприятливо) Чиста ΔG_fold ≈ −5 до −15 ккал/моль для типового стабільного білка Теплова енергія k_B T ≈ 0,6 ккал/моль при 300 К Гідрофобний ефект (домінівна рушійна сила): Неполярні залишки заховані в ядрі → впорядкована водна оболонка вивільняється → великий приріст ентропії Перенесення неполярної поверхні: ΔG ≈ 25 кал/(моль Ų) Ядро складається ~на 50% з гідрофобних залишків (Лей, Вал, Іле, Фен, Мет, Ала) Швидкості згортання (кореляція Тірумалая 1995, Плаксо 1998): Двостанові фолдери: ln k_f корелює з σ (відносний порядок контактів) Швидкі фолдери: мкс–мс (напр., головний домен вілліну 35 а.к. згортається за ~5 мкс) Повільні фолдери: хвилини–години (пріони, великі β-структурні білки) Молекулярні шаперони: Hsp70 (DnaK): зв'язує відкриті гідрофобні ділянки, АТФ-залежний GroEL/GroES (Hsp60): клітка ізолює неправильно згорнутий білок, 2 цикли гідролізу АТФ ~30% новосинтезованих білків E. coli взаємодіють з шаперонами AlphaFold2 (Jumper та ін., Nature 2021): Архітектура трансформера, навчена на 170 000 структур PDB Передбачено ~200 млн структур у базі AlphaFold DB (охоплює більшість відомих родин білків) Точність: медіанний TM-score ~0,92 на мішенях CASP14 Ключові інновації: Evoformer (коеволюційна увага), інваріантна точкова увага Залишкові виклики: внутрішньо невпорядковані ділянки, комплекси білок-РНК, мембранна динаміка
Неправильне згортання білків лежить в основі амілоїдних хвороб: Альцгеймера (Aβ/tau), Паркінсона (α-синуклеїн), діабету 2 типу (IAPP) та пріонних хвороб (PrP). У кожному випадку нормально розчинний білок набуває багатого на β-структуру фібрилярного стану, який термодинамічно стабільніший за нативну структуру, але біохімічно токсичний.
3. АТФ-синтаза: роторний двигун, керований протонами
АТФ-синтаза (F_0F_1-АТФаза) — найважливіший фермент біології: вона синтезує ~90% АТФ, що використовується практично кожною еукаріотичною та більшістю бактеріальних клітин. Вона працює за рахунок протон-рушійної сили — поєднання градієнта pH та мембранного потенціалу на внутрішній мітохондріальній мембрані — уперше описаної Пітером Мітчеллом у його хемиосмотичній гіпотезі (Нобелівська премія 1978). Фермент — буквально роторний двигун: ротор c-кільце обертається зі швидкістю 100–400 об/хв, а Пол Бойєр (Нобелівська премія 1997) показав, що всі три каталітичні сайти в голівці F_1 працюють послідовно через механізм зміни зв'язування.
Термодинаміка та обертання АТФ-синтази
Протон-рушійна сила (ПРС): ΔG_pmf = ΔG_p + ΔG_pH ΔG_p = F Δψ (електрична складова; Δψ ≈ −150 до −180 мВ у мітохондріях) ΔG_pH = −2,3 RT ΔpH (ΔpH ≈ 0,5–1,0 одиниці ≈ 30–60 мВ еквівалент) Загальна ПРС ≈ −200 мВ ≈ −4,6 ккал/моль на протон Термодинаміка синтезу АТФ: АДФ + Фн → АТФ + H_2O ΔG°' = +7,3 ккал/моль (стандартна), у клітині ΔG ≈ +12–14 ккал/моль Протонів на одну АТФ: n = ΔG_ATP / ΔG_pmf ≈ 3 (у ссавців: c-кільце має 8 субодиниць, 8/3 ≈ 2,7 H+/АТФ) Стехіометрія ротора F_0 (кількість субодиниць c-кільця, залежить від організму): Мітохондрії бика: 8 c-субодиниць (2,67 H+/АТФ) Хлоропласт CF_1CF_0: 14 c-субодиниць (4,67 H+/АТФ) E. coli: 10 c-субодиниць (3,33 H+/АТФ) Покрокове обертання (одномолекулярні експерименти, Ноджі 1997): F_1 обертається кроками по 120° (3-кратна симетрія кільця α_3β_3) Кожен крок 120° = одна синтезована АТФ (або гідролізована у зворотному напрямку) Підкрок: 80° (зв'язування АТФ) + 40° (вивільнення Фн) Гальмівний крутний момент: ~40 пН·нм; ефективність ~100% (близько до термодинамічної оборотності) Виробництво АТФ на молекулу глюкози (повне окиснення): Гліколіз: 2 АТФ (чисто) Декарбоксилювання пірувату + цикл Кребса: 2 АТФ + 8 НАДН + 2 ФАДН_2 Окисне фосфорилювання: 32 АТФ (теоретично; ~30 спостерігається) Загалом: ~36–38 АТФ на молекулу глюкози
4. Динаміка цитоскелета: актин і мікротрубочки
Цитоскелет — структурний і механічний каркас клітини. Він складається з трьох типів філаментів: актинові мікрофіламенти (~7 нм), проміжні філаменти (~10 нм, переважно структурні) та мікротрубочки (~25 нм). Актин і мікротрубочки перебувають у стані постійної динамічної рівноваги — полімеризуючись на одному кінці й деполімеризуючись на іншому — явище, зване тредмілінгом, а для мікротрубочок — динамічною нестабільністю.
Кінетика полімеризації та динамічна нестабільність
Кінетика полімеризації актину: k_on (плюс-кінець) = 11,6 мкМ^-1с^-1; k_off = 1,4 с^-1 Критична концентрація (плюс-кінець): c_crit = k_off/k_on ≈ 0,12 мкМ (АТФ-актин) Тредмілінг: чисте додавання на плюс-кінці, чиста втрата на мінус-кінці Гідроліз АТФ після полімеризації: t_½ ≈ 2 с → АДФ-актин слабший Профілін прискорює додавання на плюс-кінці; кофілін розрізає АДФ-актин Динамічна нестабільність мікротрубочок (Мітчісон і Кіршнер, 1984): ГТФ-тубулін полімеризується на плюс-кінці, утворює ГТФ-шапку Гідроліз ГТФ до ГДФ: затримка ~100 с (стохастична) Втрата ГТФ-шапки → катастрофа (швидка деполімеризація зі швидкістю ~30 нм/с) Рятування: перехід назад до росту (→ ~1 нм/с) Частота катастроф: ~0,05 с^-1; частота рятувань: ~0,15 с^-1 Генерація сили полімеризацією (модель броунівського храповика): Сила актинового філамента: F = k_B T / δ · ln(k_on [G-актин] / k_off) δ = розмір мономера (~2,7 нм для актину) Гальмівна сила: ~1–2 пН (передній край актину, ламелоподії) Мікротрубочка: до ~5 пН (проти кінетохори під час сегрегації хромосом) Розмір кроку та сила кінезину: Кінезин-1 ходить, перебираючи «руками»: кроки по 8 нм уздовж мікротрубочки (один димер тубуліну) Швидкість: ~800 нм/с, витрачаючи 1 АТФ на крок Гальмівна сила: ~7 пН (виміряно оптичним пінцетом, Свобода та ін., 1993) Процесивність: ~1 км на «шосе» (сотні послідовних кроків до від'єднання)
5. CRISPR-Cas9: точне редагування геному
CRISPR-Cas9 (короткі паліндромні повтори, регулярно згруповані кластерами) було відкрито як бактеріальну адаптивну імунну систему, що зберігає фрагменти вірусної ДНК і використовує їх для розпізнавання та розрізання вірусних загарбників, що повертаються. Дженніфер Дудна та Еммануель Шарпантьє адаптували її в програмований редактор геному (Нобелівська премія 2020), для якого потрібна лише направляюча РНК (гРНК), щоб спрямувати ендонуклеазу Cas9 до будь-якої послідовності, що відповідає 20-нуклеотидному протоспейсеру, суміжному з мотивом PAM.
Механізм CRISPR-Cas9 та кінетика пошуку
Структура направляючої РНК: crRNA (CRISPR РНК): 20-нт спейсер + повторюваний каркас tracrRNA: активуючий каркас, злитий із crRNA у конструкції sgRNA sgRNA (єдина направляюча РНК): злиті crRNA+tracrRNA, ~100 нт загалом Розпізнавання PAM (NGG для SpCas9): 3D-дифузія: Cas9-sgRNA обшукує ~10^7 сайтів PAM у геномі людини (~60 мільйонів NGG) 1D-ковзання: Cas9 перевіряє 2–3 п.о., суміжні з кожним PAM Час перебування на нецільовому сайті: ~0,1 с; на цільовому: ~години (утворюється стабільна R-петля) Загальний час пошуку: ~6 год (модель Штернберга та ін., 2014) Утворення R-петлі (вторгнення ланцюга): sgRNA спарюється з комплементарним нематричним ланцюгом (3'→5') Витіснений ланцюг утворює R-петлю; домен HNH розрізає матричний ланцюг Домен RuvC розрізає нематричний ланцюг Розріз утворює тупокінцевий подвійний розрив (DSB) на 3 п.о. вище PAM Наслідки репарації ДНК: NHEJ (негомологічне з'єднання кінців): швидке, але схильне до помилок; вставки/делеції (indels) → вимкнення гена HDR (репарація, спрямована гомологією): потребує донорської матриці; точні редагування, поодинокі нуклеотиди Ефективність HDR: 0,1–5% у клітинах людини (переважає в фазах S/G2) Інструменти CRISPR наступного покоління: Базові редактори (Комор, 2016): дезаміназа, злита з dCas9; C→T або A→G без DSB Прайм-редагування (Анзалоне, 2019): зворотна транскриптаза + pegRNA; вставки, делеції, усі 12 трансверсій/транзицій CRISPRi / CRISPRa: dCas9, злитий із репресором/активатором; модуляція експресії без редагування Схвалені терапії: exagamglogene autotemcel (exa-cel) для серпоподібноклітинної анемії та β-таласемії (FDA, 2023)
Перша терапія на основі CRISPR, схвалена FDA (грудень 2023, Casgevy), лікує серпоподібноклітинну анемію та β-таласемію, редагуючи гемопоетичні стовбурові клітини для реактивації фетального гемоглобіну. Дані лікування 2023 року показали, що 28 з 29 пацієнтів із серпоподібноклітинною анемією були вільні від вазооклюзійних кризів через 12 місяців після лікування — віха в історії медицини.
Пов'язані симуляції
Редагування генів CRISPR
Покрокова візуалізація направляючої РНК, зв'язування Cas9, розпізнавання PAM та розрізання ланцюга.
Потенціал дії серця
Модель іонних каналів Ходжкіна-Хакслі: стробування Na+, K+, Ca2+ та динаміка аритмії.
Молекулярна динаміка
Частинки Леннарда-Джонса: рівноважний стан, фазові переходи, коефіцієнти дифузії.
Нейронні коливання
Модель зв'язаних осциляторів мозкових ритмів: синхронізація тета-, альфа- та гама-діапазонів.