Підручник МЕДИЧНА ІНФОРМАТИКА

6. Математичне моделювання в імунології.

Імунітет – це складний комплекс реакцій організму на вторгнення антигенів – чужорідних об’єктів або клітин, тканин, білків, що переродилися. Специфічна імунна реакція на молекулярному рівні починається з того, що спеціалізовані плазматичні клітини виробляють у великій кількості білкові молекули – антитіла, які нейтралізують антигени.

Розглянемо модель роботи імунного апарату під час тривалого інфекційного захворювання. Ця модель використовується в клінічній практиці при лікуванні вірусного гепатиту і гострої пневмонії.

Взаємодія антигенів і імунних сил організму в математичній моделі має характер, подібний до взаємодії в системі «хижаки-жертви». «Жертвами» тут є чужорідні об’єкти – антигени з концентрацією Х. «Хижаками» є антитіла з концентрацією N, що створюються плазматичними клітками з концентрацією P. Така модель часто використовується в клінічній практиці при лікуванні вірусного гепатиту, гострої пневмонії та ін. Врахуємо очевидні міркування:

коефіцієнт розмноження антигенів m_р (вірусів, бактерій) має обернено-пропорційну залежність від температури, що пов’язано із стримуючим впливом високої температури на їх розмноження;
природний розпад антитіл і антигенів пропорційний їх концентраціям з коефіцієнтами n_с і m_с , відповідно;
природна загибель плазматичних клітин пропорційна їх концентрації з коефіцієнтом p_с ;
взаємодія антиген-антитіло в реакції аглютинації пропорційно імовірності зустрічі відповідного антитіла з антигеном, тобто добутку Х×N;
надходження антитіл до крові пропорційно концентрації плазматичних кліток P;
швидкість народження плазматичних кліток є функцією концентрації антигенів F(Х) з коефіцієнтом p_p, що має прямо-пропорційну залежить від температури.

З урахування цього можна скласти систему диференціальних рівнянь, що описують тривалий процес інфекційного захворювання:

dХ/dt = m_р×Х – b×Х×N – m_с×Х,

dN/dt = n_р×Х – k×Х×N – n_с×Х,

dP/dt = p_p×F(Х) – p_с×P,

де b – коефіцієнт, що враховує імовірність нейтралізації антигену антитілами при зустрічі з ними; n_р – коефіцієнт розмноження антитіл; k – коефіцієнт, що враховує зменшення числа антитіл за рахунок їх зв’язку з антигенами,

Дослідження характеру рішень математичної моделі дали чотири основні форми перебігу інфекційного захворювання. На рис. 68 показано можливі випадки динаміки імунної реакції (Х – кількість антигенів, t – час).

Рис. 47. Перебіг інфекційного захворювання.

Субклінічна форма (1) – проходить без фізіологічних порушень в організмі і без зовнішніх проявів. Засоби імунного захисту легко знищують антигени, не даючи їм розмножитися до небезпечної межи.

Гостра форма (2) спостерігається коли організм атакується невідомим антигеном у великих кількостях. Спочатку відбувається його інтенсивне розмноження. Коли ж імунна система виробляє проти антигену достатню кількість антитіл, чисельність антигенів різко зменшується.

Хронічна форма (3) – установлюється динамічна рівновага антигенів і антитіл. Виникає стійкий стан хвороби.

Летальна форма (4) – імунна реакція занадто затримується, і велика кількість антигенів приводить до необоротних змін в організмі.

Дослідження математичної моделі полягає у розв’язанні отриманої системи диференціальних рівнянь при відомих коефіцієнтах і з відомими початковими умовами Х(0), N(0), Z(0). Значення коефіцієнтів визначають за результатами спеціальних біохімічних аналізів, для кожної людини вони індивідуальні. Дуже важливо, що така система рівнянь при різних початкових умовах і коефіцієнтах показує абсолютно різну динаміку процесу.

Наприклад, у медичній практиці, лікування деяких інфекційних захворювань проводять методом загострення, тобто переведенням хронічної форми в гостру з наступним одужанням (рис. 69). Для цього потрібно штучно загострити хворобу, тобто ввести в організм у певні моменти часу (t₁, t₂ ) певну кількість Р біостимулятора – антигену, який є конкуруючим, непатогенним і не розмножується. Через деякий час він породжує посилену імунну відповідь, що приводить до швидкого одужання.

Рис. 48. Перебіг захворювання

Дослідження математичної моделі на ЕОМ дає змогу визначити кількість біостимулятора і момент часу його введення в організм хворого, при яких перебіг хвороби набуває потрібної форми.

Переведення хронічної форми в гостру можна здійснити за допомогою температурного ефекту; гіпер- чи гіпотермії, оскільки коефіцієнти m_р i z_р , які відповідають за розмноження антигенів і утворення плазматичних клітин, залежать від температури. Тобто, змінюючи штучно температуру організму чи фізіологічних засобів, які не мають побічної дії на імунну систему, можна досягти потрібного результату. В цьому випадку теж не обійтися без дослідження математичної моделі на ЕОМ. Багаторазове прораховування моделі при різних значеннях температури Т допоможе знайти таку, при якій графік перебігу хвороби набуває потрібної форми.

Содержание