Навіщо симулювати молекулярну біологію?
Молекулярна біологія працює в масштабах, невидимих для будь-якого мікроскопа, який можна навести на шкільний зразок. Одна вилка реплікації ДНК рухається зі швидкістю ~1 000 пар основ за секунду; рибосома транслює 10–20 амінокислот за секунду; CRISPR-Cas9 знаходить свою 20-нуклеотидну мішень у геномі з 3 мільярдів пар основ. Ці числа мають значення, бо розкривають механізм життя — але вибудувати навколо них інтуїцію лише з діаграм складно.
Анімація перекриває цей розрив. Коли можна спостерігати, як хеліказа розкручує подвійну спіраль, бачити, як фрагменти Оказакі з'являються позаду полімерази відстаючого ланцюга, а потім сповільнити анімацію, щоб розглянути кожен фермент по черзі, підручникова діаграма перетворюється на процес, який ви розумієте, а не на картинку, яку ви завчаєте.
Цей огляд проведе вас крізь шість симуляцій молекулярної біології в колекції, згрупованих за масштабом: спочатку сама молекула ДНК, потім клітина в цілому, потім спадковість через покоління, потім ферментативна хімія і, нарешті, редагування генів.
Шар 1: молекула ДНК
Реплікація ДНК
Кожного разу, коли клітина ділиться, вона спочатку повинна зробити точну копію всього свого геному. Симуляція реплікації ДНК анімує вилку реплікації: хеліказа розкручує й розділяє подвійну спіраль; примаза закладає короткі РНК-праймери; ДНК-полімераза III безперервно рухається вздовж лідируючого ланцюга й синтезує відстаючий ланцюг окремими фрагментами Оказакі; а ДНК-лігаза зшиває розриви, щойно кожен фрагмент Оказакі завершено.
Ключові рівняння: напівконсервативна реплікація
Комплементарне спарювання основ: A ↔ T (2 водневі зв'язки) G ↔ C (3 водневі зв'язки) Напрямок полімеризації: тільки 5′ → 3′ Лідируючий ланцюг: безперервний синтез у напрямку вилки Відстаючий ланцюг: перервний — фрагменти Оказакі ~200 нт (еукаріоти) Частота помилок Pol III (до коректури): ~10⁻⁵ Частота помилок після коректури 3′→5′ екзонуклеазою: ~10⁻⁷ Кінцева геномна частота помилок (+ репарація невідповідностей): ~10⁻¹⁰
Симуляція дозволяє налаштувати швидкість відтворення й наблизитись до окремих пар основ. Увімкніть "режим пар основ", щоб побачити кожну пару A–T і G–C у міру проходження хеліказы. Лічильник прогресу підраховує реплікові пари основ у реальному часі.
Транскрипція ДНК
Щойно послідовність ДНК існує, клітина зчитує обрані гени й перетворює їх на матричну РНК — процес, який називається транскрипцією. Симуляція транскрипції ДНК показує рух РНК-полімерази вздовж матричного ланцюга, збираючи комплементарні нуклеотиди мРНК (A→U, T→A, G→C, C→G). Коли полімераза рухається, подвійна спіраль ДНК знову формується позаду неї, а зростаючий ланцюг мРНК відшаровується.
Правила та швидкості транскрипції
Спарювання основ ДНК → мРНК (матричний ланцюг зчитується 3′→5′):
A → U T → A G → C C → G
Швидкість елонгації Pol II еукаріотів: ~20–60 нт/с
Швидкість елонгації Pol прокаріотів: ~40–80 нт/с
Центральна догма (спрощено):
ДНК ──транскрипція──▶ мРНК ──трансляція──▶ Білок
Кодон: триплет основ мРНК; 4³ = 64 можливі кодони
(кодують 20 амінокислот + 3 стоп-кодони)
Реплікація ДНК
Хеліказа, Pol III, примаза, лігаза — анімована вилка реплікації з фрагментами Оказакі та режимом пар основ.
Транскрипція ДНК
Відстеження РНК-полімерази вздовж матричного ланцюга; збірка мРНК із кольоровим позначенням пуринів/піримідинів.
Шар 2: клітина
Мітоз і мейоз
Читання геному — лише одна з вимог. Клітини також повинні точно ділитися, щоб кожна дочірня клітина отримала повний геном. Симуляція мітозу та мейозу анімує обидва типи поділу через усі фази: інтерфаза → профаза → метафаза → анафаза → телофаза → цитокінез для мітозу; а також додаткові специфічні для мейозу стадії, що зменшують кількість хромосом удвічі для статевого розмноження.
Мітоз проти мейозу — ключові відмінності
Мітоз (7 фаз): Інтерфаза → Профаза → Метафаза → Анафаза → Телофаза → Цитокінез Дочірні клітини: 2 × диплоїдні (2n) Мета: ріст, репарація Мейоз (10 фаз, 2 поділи): Мейоз I (редукційний): гомологи розходяться → 2 × гаплоїдні (n) Мейоз II (еквіваційний): хроматиди розходяться → 4 × гаплоїдні (n) Ключова подія — кросинговер: обмін несестринськими хроматидами під час профази I → генетична рекомбінація Приклад кількості хромосом (людина): 2n = 46 → Мейоз → n = 23 (гамети)
Симуляція включає панель фазових перемикачів, тож можна перейти до будь-якої стадії, крокувати вручну або запустити анімацію з налаштовуваною швидкістю. Веретено поділу, ядерна оболонка, хромосоми та борозна поділу — усе відтворено на полотні з реалістичними пропорціями.
Міжкатегорійний зв'язок: мейоз — це також місце, де починається менделівська генетика — випадкове розходження гомологічних хромосом у різні гамети є фізичним механізмом, що лежить в основі Першого закону Менделя про сегрегацію.
Шар 3: спадковість через покоління
Менделівська генетика та рівновага Харді-Вайнберга
Експерименти Менделя з горохом у 1860-х встановили правила спадковості ще до того, як хтось знав, що таке ДНК. Симуляція менделівської генетики реалізує моно- та дигібридну решітку Пеннета, рівняння рівноваги Харді-Вайнберга та генератор Монте-Карло, що витягує 200 нащадків із заданих користувачем частот генотипів.
Ключові закони та рівняння
Перший закон Менделя (сегрегація): Кожен організм несе два алелі; вони сегрегують рівномірно в гамети → p(A) + p(a) = 1 Другий закон Менделя (незалежне комбінування): Алелі різних генів комбінуються незалежно (коли розташовані на негомологічних хромосомах) Рівновага Харді-Вайнберга (без відбору, мутацій, дрейфу): p² + 2pq + q² = 1 де p = частота(A), q = частота(a) → частоти генотипів: AA = p², Aa = 2pq, aa = q² Умова відхилення від Х-В: ΔГет = 2pq(поточне) − 2pq(Х-В) ≠ 0 (коефіцієнт інбридингу F вимірює надлишок гомозиготності)
П'ять пресетів ознак (колір пелюсток, форма насіння, група крові ABO, кодомінантність зросту та серпоподібноклітинна анемія) охоплюють як прості менделівські, так і складніші сценарії кодомінантності та неповного домінування. Пресет групи крові коректно реалізує трьохалельну систему IA/IB/i.
Шар 4: ферментативна хімія
Ферментативна кінетика — Міхаеліс-Ментен
Клітини не просто копіюють інформацію — вони постійно проводять хімічні реакції, всі каталізовані ферментами. Симуляція ферментативної кінетики візуалізує модель Міхаеліса-Ментен взаємодії фермент-субстрат, будує подвійно-обернений графік Лайнуївера-Берка й демонструє три класи інгібування: конкурентне, неконкурентне та безконкурентне.
Кінетика Міхаеліса-Ментен
Схема реакції: E + S ⇌ ES → E + P
k₁ k₂
k₋₁
Рівняння Міхаеліса-Ментен:
v = Vmax · [S] / (Km + [S])
Vmax = k₂ · [E]загальна (максимальна швидкість при насиченні субстратом)
Km = (k₋₁ + k₂) / k₁ (концентрація субстрату при ½ Vmax)
Лайнуївер-Берк (подвійна обернена):
1/v = (Km/Vmax) · (1/[S]) + 1/Vmax
→ перетин з y: 1/Vmax; перетин з x: −1/Km
Типи інгібування (приклад конкурентного):
v = Vmax · [S] / (Km·α + [S]) де α = 1 + [I]/Ki
→ уявне Km зростає; Vmax не змінюється
Вигляд кривої виснаження субстрату показує, як продукт накопичується з часом у міру споживання субстрату — корисне доповнення до звичного стаціонарного вигляду v проти [S], що зазвичай подається в підручниках.
Шар 5: редагування генів
Редагування генів CRISPR-Cas9
CRISPR-Cas9 — найзначніший розвиток біотехнології з часів ПЛР. Симуляція CRISPR анімує молекулярний механізм: направляюча РНК (гРНК) шукає в геномі комплементарну 20-нуклеотидну протоспейсерну послідовність; Cas9 фіксується, розкручує локальну ДНК, перевіряє ділянку PAM (NGG для SpCas9) і робить дволанцюговий розрив. Потім клітина відновлює розрив або негомологічним з'єднанням кінців (NHEJ, схильне до помилок), або гомологічно-спрямованою репарацією (HDR, точна).
Механізм CRISPR-Cas9
Компоненти:
гРНК = 20-нт спейсер (відповідає мішені) + каркасна РНК
Cas9 = ендонуклеаза з доменами нуклеази RuvC і HNH
PAM = мотив, суміжний з протоспейсером (5′-NGG-3′ для SpCas9)
Розпізнавання мішені:
Спейсер гРНК спарюється з неспецифічним ланцюгом
PAM має бути безпосередньо 3′ від мішені на неспецифічному ланцюгу
Cas9 розкручує ~20 п.о. і перевіряє комплементарність пар основ
Розщеплення:
Домен HNH розрізає неспецифічний ланцюг (+3 нт від PAM)
Домен RuvC розрізає матричний ланцюг
→ тупокінцевий дволанцюговий розрив (DSB)
Шляхи репарації:
NHEJ: схильне до помилок з'єднання → інделі → нокаут гена
HDR: надається гомологічний шаблон → точне редагування
вимагає перебування клітини у фазі S/G2
Симуляція включає накладення оцінки позацільової дії: невідповідності в насіннєвій ділянці (позиції 1–10, рахуючи від PAM) зважуються вище, ніж проксимальні невідповідності, відображаючи емпіричні дані про точність Cas9. Змініть послідовність гРНК і дивіться, як прогнозована оцінка позацільової дії змінюється в реальному часі.
Повна колекція молекулярної біології
Реплікація ДНК
Анімація вилки реплікації: хеліказа, Pol III, фрагменти Оказакі, лігаза. Режим пар основ.
Транскрипція ДНК
РНК-полімераза рухається вздовж матричного ланцюга; збірка мРНК з A–T–G–C.
Мітоз і мейоз
Усі 7 фаз мітозу + усі 10 фаз мейозу. Фазові перемикачі, кросинговер, анімація веретена.
Менделівська генетика
Решітки Пеннета, рівновага Харді-Вайнберга, нащадки Монте-Карло, 5 пресетів ознак.
Ферментативна кінетика
Крива v проти [S] Міхаеліса-Ментен, Лайнуївер-Берк, конкурентне/неконкурентне інгібування.
CRISPR-Cas9
Пошук мішені гРНК, перевірка PAM, дволанцюговий розрив, шляхи репарації NHEJ проти HDR.
Зв'язки з іншими колекціями
Молекулярна біологія не існує ізольовано. Кілька симуляцій в інших категоріях безпосередньо пов'язані з темами цього огляду:
- Синаптична передача — білки, що формують іонні канали й машинерію злиття везикул, транслюються з генів тим самим механізмом, що показує симуляція транскрипції ДНК.
- Фармакокінетика — метаболізм ліків переважно керується ферментами цитохрому P450; їхня кінетика слідує кривим Міхаеліса-Ментен, ідентичним тим, що в симуляції ферментативної кінетики.
- Броунівський рух — теплові флуктуації на масштабі окремих молекул спричиняють випадкові блукання, що зводять ферменти й субстрати в контакт; той самий коефіцієнт дифузії з'являється в обох симуляціях.
- Динаміка хижак-жертва — рівняння Лотки-Вольтерри, що описують екологічні цикли популяцій, математично еквівалентні рівнянням ферментативної кінетики (обидва є законами швидкості за принципом діючих мас).