Devlog #62 – Wave 42: Cytoszkielet, Sygnalizacja Komórkowa i Odpowiedź Immunologiczna

Statystyki Wydania

Nowe Symulacje

Cytoszkielet i Ruchliwość Komórki

Pełzanie komórki to jeden z najbardziej złożonych procesów mechanicznych w biologii. Wymaga skoordynowanej aktywności setek filamentów aktynowych, dziesiątek silników molekularnych i setek kompleksów adhezyjnych — wszystko działające w skalach czasowych od milisekund do minut. Symulacja destyluje to do trzech kluczowych procesów: protruzji (polimeryzacja aktyny na krawędzi natarcia), adhezji (tworzenie kompleksu adhezji ogniskowej u podstawy komórki) i retrakcji (skurcz miozyny II włókien naprężeniowych ciągnący tył komórki do przodu).

Treadmilling filamentu aktynowego jest fundamentalnym silnikiem: na końcu kolczastym (plus) podjednostki aktyny-ATP dodają się z prędkością proporcjonalną do lokalnego stężenia monomeru; na końcu szpiczastym (minus) aktyna-ADP odłącza się. Gdy polimeryzacja na końcu kolczastym napotyka opór (błonę komórkową), powstała siła — zapadka polimeryzacyjna — pcha błonę na zewnątrz, generując lamellipodia.

Szczegóły techniczne

• Kontur komórki renderowany jako gładka zamknięta krzywa sklejana z lekkim chwianiem napędzanym tempem protruzji. Lamellipodia wysuwa się na ~30 px w kierunku gradientu sygnału (regulowanego suwakiem).
• 90 segmentów filamentów aktualizowanych na klatkę: każdy ma wierzchołek wzrostu, który polimeryzuje (tip.len += k_on * ATP * dt) i retrakcjonuje (tip.len -= k_off * dt). Filamenty, które osiągną zerową długość, są usuwane; nowe zarodkują na krawędzi natarcia.
• 24 silniki miozyny II interpolują wzdłuż włókien naprężeniowych w cyklu życia 600 klatek, wywierając widoczną siłę ślizgową (pokazaną jako zmiana koloru włókna z bursztynowego na pomarańczowy przy wysokiej aktywności miozyny).
• Prędkość komórki obliczana jako v = v_protruzji − v_retrakcji, podawana w µm/min (jednostkach znormalizowanych). Typowe pełzające komórki poruszają się z prędkością 0,1–1 µm/min in vivo.
• Prawy panel pokazuje wykres szeregu czasowego 175×120 dla prędkości komórki wraz z sześcioma statystykami w czasie rzeczywistym: tempo protruzji, tempo retrakcji, aktywne filamenty, aktywność silników, liczba adhezji i zużycie ATP.

Kaskada Sygnalizacji Komórkowej (MAPK/ERK)

Kaskady transdukcji sygnału osiągają dwa niezwykłe wyczyny: skrajne wzmocnienie (jeden receptor → wiele efektorów) i ostrą odpowiedź progową (przełącznikową). Kaskada MAPK/ERK jest najlepiej zbadanym przykładem i kontroluje proliferację komórek, różnicowanie i przeżycie.

Gdy EGF wiąże się ze swoim receptorem EGFR, receptor dimeryzuje i autofosforyluje się na wielu resztach tyrozynowych. To rekrutuje białka adaptorowe Grb2 i SOS, które katalizują wymianę GDP→GTP na małej GTPazie Ras. Aktywna Ras-GTP rekrutuje Raf do błony, gdzie fosforyluje MEK na dwóch resztach serynowych. Podwójnie ufosforylowany MEK następnie fosforyluje ERK na treoninie i tyrozynie, produkując aktywny ERK. Połączenie dwóch wymaganych etapów fosforylacji przy MEK→ERK tworzy współczynnik Hilla ~2, nadając kaskadzie przełącznikową odpowiedź dawka-efekt.

Szczegóły techniczne

• Stężenia śledzone jako znormalizowane frakcje 0–1. Układ ODE całkowany metodą Eulera z dt = 0,04 na klatkę animacji (efektywnie ~2,4 s/klatkę czasu biologicznego).
• Kinetyka Hilla: [ERK*] = [MEK**]² / (Km² + [MEK**]²) z K_m = 0,4 i wykładnikiem kooperatywności n = 2.
• Białko rusztowania KSR modelowane jako multiplikatywny współczynnik wydajności: rate × (1 + scaffold × 1,2), zgodny z danymi eksperymentalnymi pokazującymi do 2× wzmocnienie tempa.
• Ujemne sprzężenie zwrotne: aktywny ERK fosforyluje i hamuje SOS (Ras-GEF), modelowane jako ligandFrac × (1 − feedback × fbStrength × 0,6).
• Animowane cząstki sygnału podróżują wzdłuż każdej strzałki kaskady z prawdopodobieństwem proporcjonalnym do poziomu aktywności węzła powyżej. Nieaktywne (sprzężenie zwrotne) cząstki są renderowane z 40% przezroczystością.
• Wykres szeregu czasowego 220×240 px pokazuje jednocześnie Ras-GTP, MEK** i ERK**, ułatwiając dostrzeżenie opóźnienia (przesunięcia fazowego) między węzłami wyższego i niższego poziomu.

Odpowiedź Immunologiczna

Symulacja to model agentowy (ABM) działający bezpośrednio w HTML5 Canvas 2D. Agenci należą do pięciu typów — patogen, neutrofil, makrofag, limfocyt T i przeciwciało — każdy z odrębną prędkością, rozmiarem i prawdopodobieństwem zabicia. Patogen podlega wzrostowi logistycznemu ograniczonemu przez zjadliwość; komórki odpornościowe są rekrutowane w trzech fazach czasowych odpowiadających biologii:

• Natychmiastowa (0–80 taktów): Sygnały dopełniacza i komórek tucznych rekrutują neutrofile.
• Wrodzona (80–250 taktów): Przybywają makrofagi, fagocytują patogeny i prezentują antygeny.
• Adaptacyjna (>180 taktów): Limfocyty T namnażają się klonalnie; wydzielane są przeciwciała (pochodzące z limfocytów B) — lub wcześniej, jeśli gospodarz jest zaszczepiony.

Szczegóły techniczne

• Do ~400 żywych agentów w danym momencie. Każda klatka uruchamia przejście ruchu O(n), a następnie przejście kolizji O(n²) (przy dt = 2 klatki dla fizyki, 1 klatka dla renderowania), by wykryć zabicia.
• Prawdopodobieństwo zabicia na kolizję: neutrofil 3–7%, makrofag 2,5–6%, limfocyt T 4–9%, przeciwciało 5–11% (skalowane suwakiem siły odporności).
• Replikacja patogenu: co 40 taktów, jeśli poniżej pojemności nośnej 5 + virulence×95, pojawia się patogen potomny sąsiadujący z agentem rodzicielskim.
• Każdy typ komórki odpornościowej ma czas życia: neutrofile 900 taktów, makrofagi 2000, limfocyty T 2800, przeciwciała 1600 — odzwierciedlając biologiczne okresy półtrwania.
• Szczepienie przyspiesza odporność adaptacyjną: 12 limfocytów T pojawia się przy t = 10 taktów (wobec t = 180 dla niezaszczepionych), a produkcja przeciwciał rozpoczyna się przy t = 60 (wobec 240).
• Prawy panel: wykres historii populacji 185×195 px pokazuje jednocześnie liczbę patogenów, liczbę neutrofili i liczbę limfocytów T na wspólnej, automatycznie skalowanej osi y.

Nowe Kategorie

Wave 42 wprowadza do platformy dwie nowe kategorie symulacji: Biologia Komórki (cytoszkielet, sygnalizacja komórkowa) i Immunologia (odpowiedź immunologiczna). Kategorie te były wcześniej pustymi lukami w spektrum pokrycia platformy. Biologia komórki i immunologia razem stanowią dużą część nowoczesnych badań medycznych, a interaktywne symulacje dają intuicję, której statyczne diagramy nie potrafią dostarczyć — widok posuwającej się do przodu lamellipodii lub obserwowanie pochłaniania patogenów przez neutrofile sprawia, że abstrakcyjne mechanizmy stają się namacalne.

Uwagi Techniczne

Wszystkie trzy symulacje to czysty HTML5/Canvas 2D, zero zewnętrznych zależności. Cytoszkielet i odpowiedź immunologiczna używają pętli requestAnimationFrame przy natywnej częstotliwości odświeżania wyświetlacza. ODE sygnalizacji komórkowej całkuje w mniej niż 0,05 ms na klatkę (Euler, 6 sprzężonych zmiennych). ABM odpowiedzi immunologicznej pomija fizykę co drugą klatkę, by utrzymać 60 fps nawet przy 400 agentach; wykrywanie kolizji O(n²) pozostaje szybkie, ponieważ liczba agentów jest ograniczona, a wewnętrzna pętla wcześnie kończy się dla martwych agentów.

Wszystkie sześć plików (trzy EN + trzy UK) stosuje standardowy układ plików symulacji: ciemny obszar canvas + prawy panel statystyk + suwaki kontrolne + przyciski presetów + karty faktów + pole wzoru + panel informacyjny. Tłumaczenia ukraińskie są kompletne, obejmując wszystkie etykiety interfejsu, nazwy presetów, tekst legendy i panel edukacyjny.

Tagi

Biologia Komórki Cytoszkielet Aktyna Miozyna Treadmilling Lamellipodia MAPK ERK Kaskada Sygnału Kinetyka Hilla Immunologia Neutrofil Makrofag Limfocyt T Model Agentowy Szczepienie Wave 42