Девлог #70 – Хвиля 50: каскад комплементу, прецесія спіна & гальванічний елемент

Хвиля 50 одним поглядом

Каскад комплементу — три шляхи до імунітету

Система комплементу складається з понад 30 білків плазми крові, що формують самопідсилювальний протеолітичний каскад. Усі три шляхи активації — класичний, альтернативний та лектиновий — сходяться на розщепленні C3:

• Класичний шлях — запускається комплексами антиген-антитіло, що зв'язують C1q; вимагає розпізнавання адаптивною імунною системою.
• Альтернативний шлях — спонтанний низькорівневий гідроліз C3 (tick-over), що підсилюється на поверхнях мікроорганізмів, позбавлених регуляторних білків.
• Лектиновий шлях — манозозв'язувальний лектин (MBL) або фіколіни зв'язують вуглеводні патерни на патогенах і активують серинові протеази MASP-1 та MASP-2.

Кожен шлях утворює C3-конвертазу, яка розщеплює C3 на C3a (анафілатоксин) і C3b. C3b депонується на поверхні патогена (опсонізація) і за допомогою додаткових етапів розщеплення формує C5-конвертазу. Розщеплення C5 вивільняє C5a (потужний хемоатрактант) і C5b, який ініціює утворення мембраноатакувального комплексу (MAC: C5b-6-7-8-9). MAC вбудовується в бактеріальну мембрану як трансмембранна пора діаметром ~10 нм, спричиняючи осмотичний лізис.

Архітектура анімації

Полотно розділене на три візуальні зони: плаваючий пул білків C3 та частинок-ініціаторів у рідкій фазі, тіло бактерії в центрі-праворуч, що представляє цільову поверхню, та шар дрейфу, куди анафілатоксини відлітають угору після розщеплення. Переходи станів частинок підпорядковуються ймовірнісному скінченному автомату:

• вільна → зв'язування, коли частинка C3/ініціатора досягає поверхні бактерії і випадковий кидок кубика падає нижче ймовірності активації (модулюється шляхом, температурою та рівнем інгібітора DAF/CD55).
• зв'язування → c3b після короткої затримки збірки; частинка перемальовується як фрагмент опсоніну, прикріплений до бактерії.
• Після накопичення достатньої кількості депозитів C3b на контурі бактерії з'являється індикатор пори MAC.
• Частинки C5a спавняться і відлітають, підраховуючись на панелі статистики.

Температура впливає через мультиплікативний коефіцієнт k_T = exp(0.035 × (T − 37)), застосований до всіх швидкостей реакцій, що відображає приблизне значення Q⊂10; ≈ 1.3 для більшості реакцій комплементу.

Прецесія спіна — рівняння Блоха та фізика МРТ

Ядерний магнітний резонанс (ЯМР) та МРТ спираються на єдиний квантово-механічний факт: ядра з ненульовим спіном мають магнітний момент μ, що прецесує навколо прикладеного поля B₀ з ларморівською частотою:

ω₀ = γ B₀

Для ¹H (протона), γ = 42.577 МГц/Тл, тож клінічний томограф 3 Тл резонує на 127.7 МГц. Інші ядра охоплюють широкий діапазон: ⁰4¹3C з 10.708 МГц/Тл використовується для спектроскопії вуглецю-13; ¹9F (40.078) та ³¹P (17.235) важливі відповідно для досліджень метаболізму ліків та фосфатної енергетики.

Рівняння Блоха

Макроскопічний вектор намагніченості M підпорядковується феноменологічним рівнянням Блоха, які додають експоненційну релаксацію до рівняння крутного моменту:

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

dM_xy/dt = −M_xy / T₂

T₁ (поздовжня, спін-ґраткова) керує відновленням M_z до рівноважного значення M₀; T₂ (поперечна, спін-спінова) керує розфазуванням обертової компоненти M_xy. У тканинах T₁ та T₂ варіюються достатньо між типами тканин, щоб створювати контраст МРТ: сіра речовина має T₁≈1.3 с, біла речовина T₁≈0.95 с; спинномозкова рідина має T₂≈2 с, тоді як м'яз — T₂≈50 мс.

Ізометрична сфера Блоха та FID

Складність візуалізації сфери Блоха на площинному полотні полягає в розбірливому відображенні тривимірних осей без контексту WebGL. Симуляція використовує просту ізометричну проєкцію: тривимірна точка (x, y, z) відображається в координати полотна через cx = px + (x − z) × cos(30°), cy = py − (x + z) × sin(30°) + y × 0.85. Слід вектора M зберігається як кільцевий буфер з 80 точок і малюється зі спадною прозорістю, простежуючи спіраль прецесії, поки M_xy спадає. Після РЧ-імпульсу панель вільної індукційної релаксації (FID) показує затухаючий синусоїдальний сигнал, який томографи МРТ оцифровують і перетворюють за Фур'є для отримання спектра в частотній області.

Гальванічний елемент — рівняння Нернста та динаміка електродів

Гальванічний (вольтаїчний) елемент перетворює спонтанну хімічну енергію окисно-відновної реакції на електричну роботу. У класичному елементі Даніеля цинк окислюється на аноді, а мідь осаджується на катоді:

Zn(тв) → Zn²⁺(вод) + 2e⁻   E° = −0.76 В

Cu²⁺(вод) + 2e⁻ → Cu(тв)   E° = +0.34 В

Стандартний потенціал елемента E°_ел = E°_катод − E°_анод = 1.10 В. За нестандартних концентрацій рівняння Нернста коригує значення з урахуванням реакційного коефіцієнта Q:

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

де R = 8.314 Дж/(моль·К), T — абсолютна температура, n — кількість переданих електронів, а F = 96485 Кл/моль — стала Фарадея. Для елемента Даніеля з Q = [Zn²⁺]/[Cu²⁺]:

• Розведення катода (низька [Cu²⁺]) зменшує E (елемент розряджається)
• Збільшення концентрації на аноді (висока [Zn²⁺]) також зменшує E (принцип Ле Шательє)
• Вища температура збільшує RT/nF, посилюючи поправку Нернста

Вільна енергія Гіббса ΔG = −nFE; коли ΔG < 0, реакція є спонтанною. Симулятор показує як ΔG° (стандартне значення), так і поточне значення E за Нернстом та струм I = E/(R_наван.).

Стратегія анімації частинок

Два напівелементи намальовані поруч із U-подібним сольовим містком між ними. Аніони (NO⊂3;&sup0; із сольового містка) дрейфують ліворуч, а катіони (K¹) дрейфують праворуч, підтримуючи електронейтральність у міру перебігу реакції. Іони металу, що розчиняються з анода, дрейфують праворуч через розчин; іони, що осідають на катоді, дрейфують ліворуч і зникають при контакті з поверхнею електрода. Електрони представлені як маленькі заряджені частинки, що пульсують уздовж зовнішнього дроту від анода до катода. Швидкість анімації масштабується з обчисленим струмом I, тож навантаження з низьким опором дає візуально швидший потік електронів.

Нотатки з реалізації

Комплемент: керування станом каскаду

Симуляція підтримує глобальний об'єкт state, що відстежує c3bCount, macCount, c5aCount та активні збірки конвертаз. На кожному кадрі анімації всі вільні частинки перевіряються відповідно до поточної ймовірності активації: P_bind = базова_швидкість × k_T × (1 − частка_інгібітора). Це дозволяє каскаду залишатися чутливим до змін повзунків без перезапуску симуляції.

Прецесія спіна: чисельне інтегрування

Рівняння Блоха інтегруються простим методом Ейлера з кроком dt, обмеженим 5 мс, щоб запобігти нестабільності при високих значеннях ω₀. Кут прецесії зростає на Δφ = ω₀ × dt на кожному кадрі, тоді як M_z та M_xy аналітично слідують своїм відповідним експоненційним відгукам у межах кожного кроку.

Гальванічний елемент: обмеження Нернста

За дуже низьких концентрацій (Q → 0) ln Q розходиться. Симулятор обмежує Q мінімальним значенням 10³⁽&sup0;⁾, щоб уникнути значень NaN у розрахунку Нернста, що відображає практичну, обмежену дифузією концентрацію поблизу поверхонь електродів.

Наступні кроки

Хвиля 51 продовжить розширювати недопредставлені категорії. Кандидати включають додаткові симуляції з імунології (активація Т-клітин, кінетика антитіл), додаткові теми з електрохімії (електроліз, концентраційні елементи) та квантовий магнетизм (обмінна взаємодія, модель Гайзенберга). Як завжди, пропозиції через сторінку контактів вітаються.