Genetyka i biologia molekularna — DNA, dziedziczenie mendlowskie, podział komórki i CRISPR

Dlaczego symulować biologię molekularną?

Biologia molekularna działa w skalach niewidocznych dla żadnego mikroskopu, jaki można skierować na szkiełko w klasie. Pojedyncze widełki replikacyjne DNA poruszają się z prędkością ~1000 par zasad na sekundę; rybosom tłumaczy 10–20 aminokwasów na sekundę; CRISPR-Cas9 znajduje swój 20-nukleotydowy cel w genomie liczącym 3 miliardy par zasad. Te liczby mają znaczenie, ponieważ ujawniają maszynerię życia — ale trudno zbudować wokół nich intuicję, patrząc tylko na diagramy.

Animacja wypełnia tę lukę. Gdy możesz obserwować, jak helikaza rozwija podwójną helisę, widzieć fragmenty Okazaki pojawiające się za polimerazą nici opóźnionej, a następnie spowolnić animację, by zbadać każdy enzym po kolei, diagram z podręcznika staje się procesem, który rozumiesz, a nie obrazkiem, który zapamiętujesz.

Ten spotlight przechodzi przez sześć symulacji biologii molekularnej w kolekcji, pogrupowanych według skali: najpierw sama cząsteczka DNA, potem cała komórka, następnie dziedziczenie na przestrzeni pokoleń, potem chemia enzymatyczna, a na końcu edycja genów.

Warstwa 1: Cząsteczka DNA

Replikacja DNA

Za każdym razem, gdy komórka się dzieli, musi najpierw wykonać dokładną kopię całego swojego genomu. Symulacja replikacji DNA animuje widełki replikacyjne: helikaza rozwija i rozdziela podwójną helisę; prymaza układa krótkie startery RNA; polimeraza DNA III podąża wzdłuż nici prowadzącej w sposób ciągły i syntetyzuje nić opóźnioną w postaci dyskretnych fragmentów Okazaki; a ligaza DNA zszywa nacięcia po ukończeniu każdego fragmentu Okazaki.

Komplementarne parowanie zasad:
  A ↔ T  (2 wiązania wodorowe)
  G ↔ C  (3 wiązania wodorowe)

Kierunek polimeryzacji: tylko 5′ → 3′
Nić prowadząca: synteza ciągła w kierunku widełek
Nić opóźniona: nieciągła — fragmenty Okazaki ~200 nt (eukarionty)

Wskaźnik błędów Pol III (przed korektą): ~10⁻⁵
Wskaźnik błędów po korekcie 3′→5′ egzonukleazą: ~10⁻⁷
Ostateczny genomowy wskaźnik błędów (+ naprawa niedopasowań): ~10⁻¹⁰

Symulacja pozwala dostosować prędkość odtwarzania i przybliżyć poszczególne pary zasad. Włącz „tryb par zasad", by zobaczyć każdą parę A–T i G–C podczas przechodzenia helikazy. Tracker postępu liczy zreplikowane pary zasad w czasie rzeczywistym.

Transkrypcja DNA

Gdy sekwencja DNA już istnieje, komórka odczytuje wybrane geny i przekształca je w informacyjny RNA — proces zwany transkrypcją. Symulacja transkrypcji DNA pokazuje postęp polimerazy RNA wzdłuż nici matrycowej, składającej komplementarne nukleotydy mRNA (A→U, T→A, G→C, C→G). Podczas przesuwania się polimerazy podwójna helisa DNA odtwarza się za nią, a rosnąca nić mRNA odrywa się.

Parowanie zasad DNA → mRNA (nić matrycowa czytana 3′→5′):
  A → U      T → A      G → C      C → G

Tempo elongacji eukariotycznej Pol II: ~20–60 nt/s
Tempo elongacji prokariotycznej Pol:   ~40–80 nt/s

Centralny dogmat (uproszczony):
  DNA ──transkrypcja──▶ mRNA ──translacja──▶ Białko

Kodon: trójka zasad mRNA; 4³ = 64 możliwe kodony
       (kodujące 20 aminokwasów + 3 kodony stop)

Warstwa 2: Komórka

Mitoza i mejoza

Odczytywanie genomu to tylko jeden z wymogów. Komórki muszą też dokładnie się dzielić, tak by każda komórka potomna otrzymała kompletny genom. Symulacja mitozy i mejozy animuje oba typy podziału przez wszystkie fazy: interfaza → profaza → metafaza → anafaza → telofaza → cytokineza dla mitozy; oraz dodatkowe etapy specyficzne dla mejozy, które zmniejszają liczbę chromosomów o połowę na potrzeby rozmnażania płciowego.

Mitoza (7 faz):
  Interfaza → Profaza → Metafaza → Anafaza
  → Telofaza → Cytokineza
  Komórki potomne: 2 × diploidalne (2n)  Cel: wzrost, naprawa

Mejoza (10 faz, 2 podziały):
  Mejoza I (redukcyjna): homologi się rozdzielają → 2 × haploidalne (n)
  Mejoza II (ekwacyjna): chromatydy się rozdzielają → 4 × haploidalne (n)
  Kluczowe zdarzenie — crossing-over: wymiana chromatyd niesiostrzanych
  podczas profazy I → rekombinacja genetyczna

Przykład liczby chromosomów (człowiek):
  2n = 46  →  Mejoza  →  n = 23 (gamety)

Symulacja zawiera panel z chipami faz, dzięki czemu można przejść do dowolnego etapu, przechodzić przez fazy ręcznie lub uruchomić animację z regulowaną prędkością. Wrzeciono, otoczka jądrowa, chromosomy i bruzda podziałowa są renderowane na canvasie z realistycznymi proporcjami.

Warstwa 3: Dziedziczenie na przestrzeni pokoleń

Genetyka mendlowska i równowaga Hardy'ego-Weinberga

Eksperymenty Mendla na roślinach grochu w latach 60. XIX wieku ustaliły zasady dziedziczenia, zanim ktokolwiek wiedział, czym jest DNA. Symulacja genetyki mendlowskiej implementuje zarówno kwadrat Punnetta dla krzyżówki monohybrydowej i dihybrydowej, równania równowagi Hardy'ego-Weinberga, jak i generator Monte Carlo losujący 200 potomków z zdefiniowanych przez użytkownika częstości genotypów.

Pierwsze prawo Mendla (segregacji):
  Każdy organizm niesie dwa allele; segregują się
  równo do gamet  →  p(A) + p(a) = 1

Drugie prawo Mendla (niezależnej dystrybucji):
  Allele różnych genów rozdzielają się niezależnie
  (gdy znajdują się na niehomologicznych chromosomach)

Równowaga Hardy'ego-Weinberga (bez selekcji, mutacji, dryfu):
  p² + 2pq + q² = 1      gdzie p = częst.(A), q = częst.(a)
  → częstości genotypów: AA = p², Aa = 2pq, aa = q²

Warunek odstępstwa od H-W:
  ΔHet = 2pq(obecne) − 2pq(H-W) ≠ 0
  (współczynnik chowu wsobnego F mierzy nadmiar homozygotyczności)

Pięć presetów cech (kolor płatków, kształt nasion, grupa krwi ABO, kodominacja wysokości i anemia sierpowata) obejmuje zarówno proste scenariusze mendlowskie, jak i bardziej złożone przypadki kodominacji i niepełnej dominacji. Preset grupy krwi poprawnie implementuje trójallelowy system IA/IB/i.

Warstwa 4: Chemia enzymatyczna

Kinetyka enzymatyczna — Michaelis-Menten

Komórki nie tylko kopiują informację; nieustannie prowadzą reakcje chemiczne, wszystkie katalizowane przez enzymy. Symulacja kinetyki enzymatycznej wizualizuje model Michaelisa-Menten interakcji enzym-substrat, rysuje wykres podwójnej odwrotności Lineweavera-Burka i demonstruje trzy klasy inhibicji: kompetycyjną, niekompetycyjną i mieszaną (typu non-competitive).

Schemat reakcji:   E + S ⇌ ES → E + P
                        k₁     k₂
                        k₋₁

Równanie Michaelisa-Menten:
  v = Vmax · [S] / (Km + [S])

  Vmax = k₂ · [E]całkowite    (maks. prędkość przy nasyceniu substratem)
  Km = (k₋₁ + k₂) / k₁      (stężenie substratu przy ½ Vmax)

Lineweaver-Burk (podwójna odwrotność):
  1/v = (Km/Vmax) · (1/[S]) + 1/Vmax
  → punkt przecięcia z osią y: 1/Vmax; z osią x: −1/Km

Typy inhibicji (przykład kompetycyjnej):
  v = Vmax · [S] / (Km·α + [S])   gdzie α = 1 + [I]/Ki
  → pozorne Km rośnie; Vmax bez zmian

Widok krzywej zubożenia substratu pokazuje, jak produkt gromadzi się w czasie, w miarę zużywania substratu — przydatne uzupełnienie typowego dla podręczników widoku v w funkcji [S] w stanie stacjonarnym.

Warstwa 5: Edycja genów

Edycja genów CRISPR-Cas9

CRISPR-Cas9 to najważniejsze osiągnięcie biotechnologiczne od czasu PCR. Symulacja CRISPR animuje mechanizm molekularny: RNA prowadzący (gRNA) przeszukuje genom w poszukiwaniu komplementarnej 20-nukleotydowej sekwencji protospacer; Cas9 zablokowuje się, rozwija lokalny fragment DNA, weryfikuje miejsce PAM (NGG dla SpCas9) i wykonuje dwuniciowe pęknięcie. Komórka następnie naprawia pęknięcie albo przez niehomologiczne łączenie końców (NHEJ, podatne na błędy), albo przez naprawę kierowaną homologią (HDR, precyzyjną).

Składniki:
  gRNA  = 20-nt spacer (dopasowany do celu) + RNA rusztowania
  Cas9  = endonukleaza z domenami nukleazowymi RuvC i HNH
  PAM   = motyw sąsiadujący z protospacerem (5′-NGG-3′ dla SpCas9)

Rozpoznanie celu:
  Spacer gRNA paruje się z nicią niematrycową
  PAM musi znajdować się bezpośrednio 3′ od celu na nici niematrycowej
  Cas9 rozwija ~20 pz i sprawdza komplementarność par zasad

Cięcie:
  Domena HNH tnie nić niematrycową (+3 nt od PAM)
  Domena RuvC tnie nić matrycową
  → tępo zakończone dwuniciowe pęknięcie (DSB)

Ścieżki naprawy:
  NHEJ:  ligacja podatna na błędy → indele → knockout genu
  HDR:   dostarczona matryca homologiczna → precyzyjna edycja
         wymaga, by komórka była w fazie S/G2

Symulacja zawiera nakładkę punktacji celów pozagenomowych (off-target): niedopasowania w regionie zalążkowym (pozycje 1–10, licząc od PAM) są ważone wyżej niż niedopasowania bliższe. odzwierciedlając empiryczne dane wierności Cas9. Dostosuj sekwencję gRNA i obserwuj, jak przewidywany wynik off-target zmienia się w czasie rzeczywistym.

Pełna kolekcja biologii molekularnej

Powiązania między kolekcjami

Biologia molekularna nie istnieje w izolacji. Kilka symulacji w innych kategoriach łączy się bezpośrednio z tym, co omawia ten spotlight:

• Transmisja synaptyczna — białka tworzące kanały jonowe i maszynerię fuzji pęcherzyków są tłumaczone z genów przez tę samą maszynerię, którą pokazuje symulacja transkrypcji DNA.
• Farmakokinetyka — metabolizm leków jest napędzany głównie przez enzymy cytochromu P450; ich kinetyka podąża krzywymi Michaelisa-Menten identycznymi jak w symulacji kinetyki enzymatycznej.
• Ruch Browna — fluktuacje termiczne w skali pojedynczych cząsteczek napędzają błądzenia losowe, które zbliżają enzymy i substraty; ten sam współczynnik dyfuzji pojawia się w obu symulacjach.
• Dynamika drapieżnik-ofiara — równania Lotki-Volterry opisujące ekologiczne cykle populacyjne są matematycznie równoważne równaniom kinetyki enzymatycznej (oba to prawa szybkości działania mas).

Algorytmy i metody w tej kolekcji