Devlog #70 – Fala 50: kaskada dopełniacza, precesja spinu & ogniwo galwaniczne

Fala 50 w skrócie

Kaskada dopełniacza — trzy drogi do odporności

Układ dopełniacza tworzy ponad 30 białek osocza, które składają się na samowzmacniającą się kaskadę proteolityczną. Wszystkie trzy drogi aktywacji — klasyczna, alternatywna i lektynowa — zbiegają się w rozszczepieniu C3:

• Droga klasyczna — wyzwalana przez kompleksy antygen-przeciwciało wiążące C1q; wymaga rozpoznania przez odporność adaptacyjną.
• Droga alternatywna — spontaniczna niskopoziomowa hydroliza C3 (tzw. tick-over), wzmacniana na powierzchniach drobnoustrojów pozbawionych białek regulatorowych.
• Droga lektynowa — lektyna wiążąca mannozę (MBL) lub fikoliny wiążą wzorce węglowodanowe na patogenach i aktywują proteazy serynowe MASP-1 i MASP-2.

Każda droga tworzy konwertazę C3, która rozszczepia C3 na C3a (anafilatoksyna) i C3b. C3b osadza się na powierzchni patogenu (opsonizacja) i w kolejnych etapach rozszczepienia prowadzi do złożenia konwertazy C5. Rozszczepienie C5 uwalnia C5a (silny czynnik chemotaktyczny) oraz C5b, który zapoczątkowuje powstanie kompleksu atakującego błonę (MAC: C5b-6-7-8-9). MAC wbudowuje się w błonę bakteryjną jako transbłonowy por o średnicy ~10 nm, powodując lizę osmotyczną.

Architektura animacji

Płótno podzielone jest na trzy strefy wizualne: unoszącą się pulę białek C3 i cząstek inicjujących w fazie płynnej, ciało bakterii w środkowo-prawej części reprezentujące powierzchnię docelową oraz warstwę dryfu, w której anafilatoksyny uciekają ku górze po rozszczepieniu. Przejścia stanów cząstek podlegają probabilistycznemu automatowi skończonemu:

• free (wolna) → binding (wiążąca się), gdy cząstka C3/inicjatora dociera do powierzchni bakterii, a losowy rzut kostką wypada poniżej prawdopodobieństwa aktywacji (modulowanego drogą aktywacji, temperaturą i poziomem inhibitora DAF/CD55).
• binding → c3b po krótkim opóźnieniu związanym ze złożeniem kompleksu; cząstka jest przerysowywana jako fragment opsoniny przytwierdzony do bakterii.
• Po nagromadzeniu wystarczającej liczby depozytów C3b na konturze bakterii pojawia się wskaźnik porów MAC.
• Cząstki C5a są generowane i odpływają, a ich liczba jest zliczana na panelu statystyk.

Temperatura wchodzi do modelu poprzez multiplikatywny czynnik k_T = exp(0.035 × (T − 37)), stosowany do wszystkich szybkości reakcji, co odzwierciedla przybliżoną wartość Q⊂10; ≈ 1.3 dla większości reakcji dopełniacza.

Precesja spinu — równania Blocha i fizyka MRI

Jądrowy rezonans magnetyczny (NMR) i MRI opierają się na jednym kwantowo-mechanicznym fakcie: jądra o niezerowym spinie mają moment magnetyczny μ, który precesuje wokół przyłożonego pola B₀ z częstością Larmora:

ω₀ = γ B₀

Dla ¹H (protonu) γ = 42.577 MHz/T, więc kliniczny skaner 3 T rezonuje na 127.7 MHz. Inne jądra obejmują szeroki zakres: ⁰4¹3C przy 10.708 MHz/T jest wykorzystywany w spektroskopii węgla-13; ¹9F (40.078) i ³¹P (17.235) mają istotne zastosowanie odpowiednio w badaniach metabolizmu leków i energetyki fosforanowej.

Równania Blocha

Makroskopowy wektor namagnesowania M podlega fenomenologicznym równaniom Blocha, które dodają eksponencjalną relaksację do równania momentu obrotowego:

dM_z/dt = −(M_z − M₀) / T₁

dM_xy/dt = −M_xy / T₂

T₁ (podłużna, spin-sieć) opisuje powrót M_z do wartości równowagowej M₀; T₂ (poprzeczna, spin-spin) opisuje rozfazowanie obracającej się składowej M_xy. W tkankach T₁ i T₂ różnią się na tyle między rodzajami tkanek, że generują kontrast MRI: istota szara ma T₁≈1.3 s, istota biała T₁≈0.95 s; płyn mózgowo-rdzeniowy ma T₂≈2 s, podczas gdy mięsień T₂≈50 ms.

Izometryczna sfera Blocha i FID

Wyzwaniem przy wizualizacji sfery Blocha na płaskim płótnie jest czytelne przedstawienie trójwymiarowych osi bez kontekstu WebGL. Symulacja wykorzystuje prostą projekcję izometryczną: punkt 3D (x, y, z) jest mapowany na współrzędne płótna poprzez cx = px + (x − z) × cos(30°), cy = py − (x + z) × sin(30°) + y × 0.85. Ślad wektora M jest przechowywany jako bufor cykliczny 80 punktów i rysowany z malejącą przezroczystością, śledząc spiralę precesji w miarę zanikania M_xy. Po impulsie RF panel swobodnego zaniku indukcji (FID) pokazuje tłumiony sygnał sinusoidalny, który skanery MRI digitalizują i przekształcają transformatą Fouriera, aby uzyskać widmo w dziedzinie częstotliwości.

Ogniwo galwaniczne — równanie Nernsta i dynamika elektrod

Ogniwo galwaniczne (woltaiczne) przekształca spontaniczną energię chemiczną reakcji redoks w pracę elektryczną. W klasycznym ogniwie Daniella cynk utlenia się na anodzie, a miedź osadza się na katodzie:

Zn(s) → Zn²⁺(aq) + 2e⁻   E° = −0.76 V

Cu²⁺(aq) + 2e⁻ → Cu(s)   E° = +0.34 V

Standardowy potencjał ogniwa E°_ogniwa = E°_katoda − E°_anoda = 1.10 V. Przy niestandardowych stężeniach równanie Nernsta koryguje wartość o iloraz reakcji Q:

E = E° − (RT / nF) · ln Q

gdzie R = 8.314 J/(mol·K), T to temperatura bezwzględna, n to liczba przenoszonych elektronów, a F = 96485 C/mol to stała Faradaya. Dla ogniwa Daniella z Q = [Zn²⁺]/[Cu²⁺]:

• Rozcieńczenie katody (niskie [Cu²⁺]) zmniejsza E (ogniwo się rozładowuje)
• Zwiększenie stężenia na anodzie (wysokie [Zn²⁺]) również zmniejsza E (reguła przekory Le Chateliera)
• Wyższa temperatura zwiększa RT/nF, wzmacniając poprawkę Nernsta

Energia swobodna Gibbsa ΔG = −nFE; gdy ΔG < 0, reakcja jest spontaniczna. Symulator wyświetla zarówno ΔG° (wartość standardową), jak i bieżące E z równania Nernsta oraz prąd I = E/(R_obc.).

Strategia animacji cząstek

Dwa półogniwa są narysowane obok siebie z mostkiem solnym w kształcie litery U pomiędzy nimi. Aniony (NO⊂3;&sup0; z mostka solnego) dryfują w lewo, a kationy (K¹) dryfują w prawo, utrzymując elektroobojętność w miarę postępu reakcji. Jony metalu rozpuszczające się z anody dryfują w prawo przez roztwór; jony osadzające się na katodzie dryfują w lewo i znikają po kontakcie z powierzchnią elektrody. Elektrony są przedstawione jako małe naładowane cząstki pulsujące wzdłuż zewnętrznego przewodu od anody do katody. Prędkość animacji skaluje się z obliczonym prądem I, tak że obciążenie o niskiej oporności daje wizualnie szybszy przepływ elektronów.

Uwagi implementacyjne

Dopełniacz: zarządzanie stanem kaskady

Symulacja utrzymuje globalny obiekt state, który śledzi c3bCount, macCount, c5aCount oraz aktywne zespoły konwertaz. W każdej klatce animacji wszystkie wolne cząstki są sprawdzane względem bieżącego prawdopodobieństwa aktywacji: P_bind = szybkość_bazowa × k_T × (1 − frakcja_inhibitora). Dzięki temu kaskada pozostaje responsywna na zmiany suwaków bez konieczności restartowania symulacji.

Precesja spinu: całkowanie numeryczne

Równania Blocha są całkowane prostą metodą Eulera z krokiem dt ograniczonym do 5 ms, aby zapobiec niestabilności przy wysokich wartościach ω₀. Kąt precesji zwiększa się o Δφ = ω₀ × dt w każdej klatce, natomiast M_z i M_xy podążają analitycznie za swoimi odpowiednimi odpowiedziami eksponencjalnymi w obrębie każdego kroku.

Ogniwo galwaniczne: ograniczenie Nernsta

Przy bardzo niskich stężeniach (Q → 0) ln Q dąży do nieskończoności. Symulator ogranicza Q do minimalnej wartości 10³⁽&sup0;⁾, aby uniknąć wartości NaN w obliczeniach Nernsta, co odzwierciedla praktyczne, ograniczone dyfuzją stężenie w pobliżu powierzchni elektrod.

Następne kroki

Fala 51 będzie kontynuować rozbudowę słabo reprezentowanych kategorii. Kandydaci obejmują dodatkowe symulacje z immunologii (aktywacja limfocytów T, kinetyka przeciwciał), dodatkowe tematy z elektrochemii (elektroliza, ogniwa stężeniowe) oraz magnetyzm kwantowy (oddziaływanie wymienne, model Heisenberga). Jak zawsze, sugestie za pośrednictwem strony kontaktowej są mile widziane.