Збудливе середовище ФітцГью-Нагумо · 2D тканина · Спіральний повторний вхід
Клацніть на тканині, щоб збудити область
Мембранний потенціал v ∈ [−2, 2]
2D-симуляція збудливого середовища серцевої тканини за моделлю ФітцХ'ю-Нагумо. Натискайте для стимуляції та спостерігайте хвилі потенціалу дії, спіральний ріентрі та фібриляцію.
Як збудливі клітини (кардіоміоцити) генерують та поширюють електричні імпульси. Модель ФітцХ'ю-Нагумо відтворює деполяризацію, реполяризацію та рефрактерний період лише двома змінними.
Натисніть на тканину для стимуляції. Другий стимул під час рефрактерного хвоста створює спіральний ріентрі — механізм серцевих аритмій. Регулюйте дифузію та часові константи.
Спіральні хвилі ріентрі в серці є провідною причиною шлуночкової тахікардії та фібриляції. Дефібриляція одночасно деполяризує всі клітини, дозволяючи синусовому вузлу відновити контроль.
Це двовимірна модель серцевої тканини Фіцх'ю-Нагумо — ідеалізоване збудливе середовище, яке відтворює, як клітини серця збуджуються та відновлюються. Кожна точка сітки підпорядковується двом зв'язаним звичайним диференціальним рівнянням: швидкій змінній v, схожій на напругу, з кубічною нелінійністю, та повільній змінній відновлення w. Зв'язування сусідніх клітин через дифузійний член перетворює ізольоване збудження на біжучі хвилі, а вчасно поданий другий стимул створює обертальне спіральне повторне збудження (reentry) — абстрактний аналог серцевої аритмії.
Рівняння реакції-дифузії dv/dt = v − v³/3 − w + D·∇²v та dw/dt = ε(v + a − b·w) інтегруються на сітці 120×120 явним прямим методом Ейлера (dt = 0.05) з межами без потоку. Поле напруги v ∈ [−2, 2] кольорово відображається від темно-синього (спокій) до білого (пік деполяризації), а живий графік відображає v у центральній клітині.
Клацніть або торкніться тканини, щоб подати стимул із регульованим радіусом. Скористайтеся кнопками Пауза/Скидання або пресетами «Спіральна хвиля» та «Імпульс». Повзунки задають дифузію D (0.1–2.0, швидкість хвилі), збудливість a (0.05–0.5), швидкість відновлення ε (0.01–0.3, тривалість рефрактерності) та радіус стимулу (2–20 пікс). Статистика показує активні клітини, кількість кроків, напругу в центрі та режим.
Модель Фіцх'ю-Нагумо була виведена у 1961 році як спрощення чотиризмінних рівнянь Ходжкіна-Гакслі до двох змінних: вона відкидає деталі іонних каналів, але зберігає поріг, спайк за принципом «все або нічого» та рефрактерний період, що роблять тканину збудливою.
Це спрощена математична модель збудливої клітини, яка описує поведінку мембрани лише двома змінними замість багатьох іонних струмів у повних рівняннях Ходжкіна-Гакслі. Швидка змінна v передає спайк напруги через кубічний член, а повільна змінна w забезпечує відновлення. Попри свою простоту вона відтворює порогове збудження, потенціал дії за принципом «все або нічого» та рефрактерний період.
Кожна клітина зв'язана зі своїми чотирма сусідами через дифузійний член D·∇²v, який тут обчислюється як дискретний лапласіан. Коли одна ділянка деполяризується, дифузія піднімає напругу сусідів вище порога, тож вони збуджуються по черзі, і збудження поширюється як біжуча хвиля. Збільшення повзунка дифузії D підвищує швидкість хвилі.
Дифузія D масштабує зв'язок між клітинами і відтак швидкість проведення. Збудливість a зсуває поріг збудження динаміки відновлення. Швидкість відновлення ε визначає, наскільки швидко релаксує w, що керує рефрактерним періодом і шириною хвилі. Радіус стимулу задає розмір у пікселях круглої ділянки, на яку впливає кожне клацання.
Спіраль виникає, коли хвиля зустрічає тканину, що з одного боку ще перебуває в рефрактерному стані. Пресет «Спіральна хвиля» налаштовує саме цю умову S1-S2, частково відновлюючи одну половину сітки, тож хвильовий фронт може лише закрутитися навколо рефрактерного краю. Утворена обертальна спіраль є модельним аналогом повторно-збуджувальних аритмій, таких як тахікардія та фібриляція.
Це якісно точно, але не кількісно реалістично. Модель свідомо абстрагується від іонних каналів, реальної геометрії клітин та анізотропії тканини, а її змінні безрозмірні, а не виміряні у мілівольтах чи мілісекундах. Вона чудово підходить для викладання принципів збудливості, поширення та повторного збудження, але клінічна електрофізіологія використовує значно детальніші іонні моделі.
Ця симуляція моделює збудливу тканину серця, використовуючи двозмінну систему реакції-дифузії Фіцх'ю-Нагумо, яка вловлює основну динаміку того, як клітини серцевого м'яза (кардіоміоцити) породжують, підтримують і поширюють електричні імпульси. Швидка змінна v, подібна до напруги, імітує спайк деполяризації мембрани через кубічну нелінійність, тоді як повільна змінна відновлення w забезпечує рефрактерний період; дифузійне зв'язування між клітинами сітки перетворює локальне збудження на біжучі хвилі потенціалу дії через двовимірний шар тканини.
Розуміння серцевих потенціалів дії — біологічна основа електрокардіографії (ЕКГ) та клінічного лікування аритмій — від ритмів водія ритму синоатріального вузла до небезпечних для життя спіральних повторно-збуджувальних патернів, що лежать в основі фібриляції шлуночків.
Серцевий потенціал дії — короткий, стереотипний сигнал напруги, що виникає, коли мембрана клітини серцевого м'яза деполяризується від потенціалу спокою приблизно −90 мВ до піку близько +20 мВ, а потім повільно реполяризується протягом 200-400 мс. Довга фаза плато — унікальна для серцевих клітин порівняно з нейронами — утримує клітину рефрактерною, доки не завершиться механічне скорочення, запобігаючи тетанічному (стійкому) скороченню, яке було б смертельним.
Натисніть кнопку пресету «Спіральна хвиля», щоб завантажити початкову умову S1-S2: нижня половина сітки повністю деполяризована, тоді як права половина має підвищену змінну відновлення w, створюючи частково рефрактерну ділянку. Хвильовий фронт може поширюватися лише там, де тканина збудлива, тож він вигинається навколо рефрактерної межі й починає обертатися. Спіраль також можна створити вручну, клацаючи по сітці, щоб подати два стимули в різний час — ключовим є другий клік під час рефрактерного хвоста першої хвилі.
Колір кодує безрозмірний мембранний потенціал v, який змінюється від −2 (спокій, темно-синій) через проміжні значення (блакитний, зелений, оранжевий) до +2 (пік деполяризації, білий). У біологічних термінах цей діапазон приблизно відповідає потенціалу спокою −90 мВ до піку потенціалу дії +20 мВ. Живий графік у нижньому лівому куті відображає v у центральній клітині сітки за останні 200 кроків симуляції.
Модель визначається двома зв'язаними ЗДР: dv/dt = v - v^3/3 - w + D * Лапласіан(v) та dw/dt = epsilon * (v + a - b*w). Кубічний член v - v^3/3 надає v бістабільний характер з різким порогом — малі збурення згасають, тоді як надпорогові стимули породжують повний спайк. Змінна відновлення w зростає під час деполяризації й повертає v до спокою, створюючи рефрактерний період. Дифузійний член D * Лапласіан(v) зв'язує сусідні клітини, тож збудження поширюється просторово. Це запропонував Річард Фіцх'ю у 1961 році як топологічне спрощення чотиризмінних рівнянь Ходжкіна-Гакслі.
Спіральне повторне збудження виникає, коли хвильовий фронт деполяризації зустрічає ділянку, що все ще рефрактерна, змушуючи його закрутитися в самопідтримувану обертальну спіраль. У реальній серцевій тканині одна стабільна спіраль спричиняє шлуночкову тахікардію (швидку, регулярну, небезпечну); множинні спіралі, що розпадаються на хаотичні патерни, спричиняють фібриляцію шлуночків, при якій скоординоване перекачування руйнується і серце неефективно тремтить. Без дефібриляції або СЛР фібриляція смертельна впродовж хвилин, тому розуміння її механізму є критично важливим для кардіології.
Частково це помилка. Вища дифузія D справді підвищує швидкість проведення, але якщо D у моделі Фіцх'ю-Нагумо підняти занадто високо, лапласіанове згладжування може повністю придушити локалізоване збудження — стимул, занадто малий відносно довжини дифузії, просто розсіється, замість того щоб викликати хвилю. У реальній серцевій тканині аналогічна величина — провідність щілинних контактів; патологічний фіброз, що роз'єднує клітини (знижує ефективне D), уповільнює провідність і сприяє повторному збудженню, тоді як рівномірно висока зв'язаність сприяє стабільним плоским хвилям.
Іонну основу збудливості встановили Алан Ходжкін та Ендрю Гакслі у 1952 році за допомогою експериментів з фіксацією напруги на гігантських аксонах кальмара — робота, за яку вони отримали Нобелівську премію з фізіології або медицини 1963 року. Денис Ноубл адаптував рамки Ходжкіна-Гакслі до серцевих волокон Пуркіньє у 1960 році, створивши першу кількісну модель довгого плато й повільної реполяризації серцевого потенціалу дії. Річард Фіцх'ю звів цю складність до двох змінних у 1961 році; Джін-Ічі Нагумо побудував електронну схему, що реалізує ті самі рівняння, у 1962 році, давши моделі її подвійну назву.
Тісно пов'язані симуляції включають модель нейрона Ходжкіна-Гакслі (повний чотиризмінний прабатько), одноклітинний осцилятор Фіцх'ю-Нагумо та формування патернів реакції-дифузії (патерни Тюрінга використовують подібну математичну структуру). Симуляції течії крові та судинної механіки описують механічний наслідок електричного збудження. Симуляції дифузії ліків і фармакокінетики релевантні, бо багато антиаритмічних препаратів діють, блокуючи певні іонні канали й змінюючи тривалість потенціалу дії чи швидкість проведення.
Обчислювальні серцеві моделі, похідні від традиції Фіцх'ю-Нагумо та Ходжкіна-Гакслі, використовуються в розробці алгоритмів імплантованих кардіовертерів-дефібриляторів (ICD), плануванні катетерної абляції при повторно-збуджувальних аритміях та платформах пацієнт-специфічних цифрових двійників серця, що прогнозують ризик аритмії за даними МРТ. Дифузійно-зв'язана PDE-рамка також лежить в основі електрокардіографічного зображення (ECGI), що неінвазивно реконструює послідовність електричної активації серця з потенціалів на поверхні тіла.
Активні напрями досліджень включають високоточні іонні моделі з понад 50 змінними стану, підігнані під записи людських шлуночкових клітин, GPU-прискорені симуляції цілого серця на анатомічно реалістичних сітках, отриманих з клінічної МРТ, машинно-навчальні сурогати, що замінюють дорогі розв'язання PDE для клінічної підтримки рішень у реальному часі, та оптогенетичне дослідження повторно-збуджувальних контурів, де активовані світлом іонні канали дозволяють дослідникам керувати часом потенціалу дії з мілісекундною точністю в живій тканині.