Чому працює математична біологія
Живі системи виглядають ззовні неймовірно складними. Помірний ліс містить тисячі взаємодіючих видів на десятках трофічних рівнів, вбудованих у мінливі абіотичні умови. І все ж напрочуд прості математичні моделі — системи пов'язаних звичайних диференціальних рівнянь — вловлюють основну якісну поведінку: коливальні цикли хижак-жертва, каскад, що йде за зникненням ключового виду, збіжність до сталого стану, коли ріст і смертність урівноважуються.
Причина в тому, що екологія, як і фізика, керується законами збереження та процесами швидкостей. Енергія надходить як первинна продукція, підіймається трофічними рівнями з ефективністю приблизно 10% на крок і зрештою розсіюється у вигляді тепла. Математика, що описує ці потоки, та сама, що застосовується в хімічній кінетиці й гідродинаміці. Ця структурна подібність робить екологічне моделювання потужним — і саме тому симуляції, побудовані на тих самих рівняннях, здатні розкривати справжнє екологічне розуміння.
Цикли хижак-жертва: рівняння Лотки-Вольтерри
Основоположну модель популяційної екології незалежно опублікували Альфред Лотка у 1925 році та Віто Вольтерра у 1926-му. Вольтерру спонукала загадка з даних адріатичного рибальства: під час Першої світової війни, коли вилов риби припинився, кількість акул зросла, а кормової риби — зменшилася, всупереч тому, що очікувала більшість людей. Його рівняння пояснили чому.
Модель хижак-жертва Лотки-Вольтерри
dN/dt = αN − βNP (здобич: народжуваність − хижацтво)
dP/dt = δNP − γP (хижак: конверсія − смертність)
N = популяція здобичі, P = популяція хижака
α = народжуваність здобичі, β = інтенсивність хижацтва
δ = ефективність конверсії хижака, γ = смертність хижака
Нерухома точка: N* = γ/δ, P* = α/β (нейтрально стійкий центр)
Період коливань ≈ 2π / √(αγ)
Модель Лотки-Вольтерри передбачає нейтрально стійкі коливання: популяції циклічно змінюються навколо нерухомої точки, не сходячись і не розходячись. Реальні системи демонструють затухаючі коливання через ємність середовища та щільнозалежні ефекти, але базовий цикл — чисельність здобичі стимулює зростання хижака, чисельність хижака спричиняє спад здобичі, спад здобичі веде до спаду хижака, що дає змогу здобичі відновитися — стійко спостерігається в природі. Класичні дані отримано з архівів хутряної торгівлі компанії Гудзонової затоки щодо популяцій рисі та зайця-біляка з 1845 по 1935 рік.
Хижак-жертва (Лотка-Вольтерра)
Агентна модель Лотки-Вольтерри на 2D-сітці: кролики пасуться й розмножуються, лисиці полюють і гинуть. Графіки популяцій оновлюються в реальному часі поряд із просторовою сіткою. Спостерігайте класичний цикл у фазовому просторі між чисельністю здобичі й хижака.
Харчова мережа
Екосистема з шести видів (Трава / Кущі / Кролик / Олень / Лисиця / Вовк) на основі розширеної Лотки-Вольтерри та RK4. Мережевий граф, де розмір вузла = популяція, панель часового ряду, перемикання видів клацанням. Пресети: стабільність / колапс хижаків / бум рослин / трофічний каскад / конкурентне витіснення.
Цикл рись-заєць триває приблизно 10 років. Популяції зайців сягають піку приблизно раз на десятиліття; популяції рисі відстають на 1–2 роки. На піку зайців може бути до 2 000 на км² — у мінімумі менше 10. Механізм частково трофічний (хижацтво рисі) і частково висхідний (перевипас кормових рослин під час піку чисельності зайців спричиняє харчовий колапс, що вбиває зайців ще до спаду популяції рисі). Симуляція відтворює обидва ефекти, коли увімкнено ємність рослинного середовища.
Трофічні каскади: як вовки змінюють русла річок
Трофічний каскад — це непрямий ефект, що поширюється харчовою мережею зверху вниз (від хижаків) або знизу вгору (від первинної продуктивності). Найвідоміший задокументований приклад — повернення сірих вовків до Єллоустоунського національного парку в 1995 році. Вовки скоротили чисельність оленів-вапіті і, що важливіше, змінили їхню поведінку — вапіті перестали переїдати прибережну рослинність, побоюючись засідки. Прибережна рослинність відновилася, береги річок стабілізувалися завдяки кореневим системам, а річкові русла звузилися й поглибшали. Повернулися бобри. Повернулися співочі птахи. Вовки, змінивши лише одну взаємодію "головний хижак — травоїдний", змінили фізичну географію парку.
Математика трофічних каскадів включає трирівневі системи Лотки-Вольтерри: рослина–травоїдний–м'ясоїдний. Рівновага рослинного рівня може контролюватися головним хижаком через непряме пригнічення травоїдних — явище, яке екологи називають "вивільненням мезохижака", коли головного хижака усувають.
Трофічний каскад
Трирівнева Лотка-Вольтерра: Рослини → Травоїдні → М'ясоїдні. Увімкніть видалення хижака, щоб побачити вибух чисельності травоїдних і колапс рослинності. Сценарій "вивільнення мезохижака" відтворює задокументовану Єллоустоунську динаміку з живими кривими часового ряду.
Зграя Boids
Емерджентна поведінка зграї з трьох локальних правил: розділення, вирівнювання, згуртованість. Алгоритм Рейнольдса (1986), 150 агентів, налаштовувані ваги та радіуси сусідства. Спостерігайте, як скоординований груповий рух виникає з суто локальних взаємодій — мікрокосм того, як виникає колективна поведінка в групах тварин.
Трирівнева Лотка-Вольтерра (трофічний каскад)
dV/dt = r·V(1 − V/K) − a₁·V·H (рослини)
dH/dt = e₁·a₁·V·H − d₁·H − a₂·H·C (травоїдні)
dC/dt = e₂·a₂·H·C − d₂·C (м'ясоїдні)
V = рослинність, H = травоїдні, C = м'ясоїдні
r = темп росту рослин, K = ємність середовища
a₁, a₂ = темпи атаки; e₁, e₂ = ефективності конверсії
d₁, d₂ = природна смертність
Циркадні ритми: біологічні годинники
Майже всі живі організми — від ціанобактерій до людей — підтримують внутрішній годинник, що відстежує ~24-годинний цикл доби. Ці циркадні ритми регулюють експресію генів, секрецію гормонів, температуру тіла, харчову поведінку та цикли сон-неспання. Порушення циркадних ритмів (позмінна робота, джетлаг, штучне освітлення вночі) корелює з підвищеним ризиком метаболічного синдрому, раку та серцево-судинних захворювань.
Молекулярний механізм розкрили Голл, Росбаш і Янг (Нобелівська премія з фізіології або медицини, 2017). Основа — петля зворотного зв'язку транскрипції-трансляції: білковий комплекс CLOCK/BMAL1 активує транскрипцію генів PER і CRY; білки PER і CRY накопичуються, пригнічують CLOCK/BMAL1, знижують власну транскрипцію, деградують, і цикл повторюється. Математичний опис — осцилятор Гудвіна: петля негативного зворотного зв'язку із затримкою в часі, що підтримує коливання.
Осцилятор циркадного ритму
Осцилятор Гудвіна, що моделює мелатонін, кортизол і базову температуру тіла протягом 24-годинного циклу. Чотири пресети: норма / позмінна робота / джетлаг на схід / джетлаг на захід. Спостерігайте, як світлове захоплення зсуває фазу годинника і як розбалансування спричиняє фазовий дрейф.
Броунівський рух
Дифузія частинок під дією теплового шуму: D = k_B T / (6πηr). Анімовані частинки, живе середньоквадратичне зміщення проти теоретичного прогнозу, повзунки температури та в'язкості. Співвідношення Ейнштейна, що пов'язує дифузію з дисипацією, — фундаментальне для молекулярної біології.
Осцилятор Гудвіна (циркадний базовий контур)
dX/dt = v₁ / (K₁ + Zⁿ) − v₂·X / (K₂ + X)
dY/dt = v₃·X − v₄·Y / (K₄ + Y)
dZ/dt = v₅·Y − v₆·Z / (K₆ + Z)
X = мРНК, Y = цитоплазматичний білок, Z = ядерний білок
n = коефіцієнт Хілла (≥ 9 для сталих коливань)
Період ≈ 24 год, задається темпами деградації v₂, v₄, v₆
Захоплення: зовнішнє світло змінює v₁ залежно від фази
Осцилятор Гудвіна — прототиповий осцилятор із негативним зворотним зв'язком: генний продукт пригнічує власне вироблення через ланцюг проміжних сполук. Стійкі коливання вимагають коефіцієнта Хілла n ≥ 9 — вкрай кооперативного нелінійного пригнічення. Симуляція дозволяє дослідити, як зміна темпів деградації зсуває період і як криві фазового відгуку на світловий імпульс (PRC) описують захоплення ритму.
Броунівський рух і дифузія в живих системах
У 1827 році ботанік Роберт Броун спостерігав хаотичний рух зерен пилку у воді. У 1905 році Ейнштейн вивів коефіцієнт дифузії D = k_B T / (6πηr) з перших принципів, використовуючи кінетичну теорію тепла, — і показав, що середньоквадратичне зміщення зростає лінійно з часом: ⟨r²⟩ = 6Dt для 3D-дифузії, 4Dt для 2D, 2Dt для 1D. Це дало перше кількісне свідчення молекулярної реальності.
У живих клітинах броунівський рух — не просто курйоз, а основний механізм транспорту на нанорівні. Молекули менші за ~5 нм дифундують достатньо швидко, щоб дослідити весь об'єм клітини за секунди. Більші комплекси й органели потребують активного транспорту молекулярними моторами (кінезин, динеїн) уздовж цитоскелетних доріжок. Межа переходу між пасивною дифузією й активним транспортом визначається балансом теплової енергії k_B T та енергії, доступної з однієї молекули АТФ (~80 пН·нм).
Клітинна сигналізація спирається на каскади дифузії. Коли гормон зв'язується з рецептором на поверхні клітини, сигнал поширюється всередину клітини через каскад фосфорилювання — кожна кіназа активує наступну. Швидкість цього каскаду зрештою обмежена дифузією активованої кінази до її субстрату. Саме тому сигнальні молекули зазвичай малі (швидка дифузія), а клітина компактна (короткі відстані дифузії). Затримка у 3 хвилини між сплеском адреналіну та піком мобілізації глікогену в м'язових клітинах здебільшого визначається дифузійним транспортом.
Екологічні мережі: від простого до складного
Реальні екосистеми — це не прості пари хижак-жертва. Харчова мережа помірного озера може містити сотні видів, з'єднаних тисячами трофічних зв'язків. Екологи характеризують такі мережі за зв'язністю C (частка можливих зв'язків, що існують), середньою довжиною ланцюга та часткою верхніх, проміжних і базових видів.
Парадокс Роберта Мея 1972 року математично показав, що випадкові мережі стають менш стійкими (у сенсі лінійної стійкості) зі збільшенням розміру і зв'язності — що суперечило тогочасному пануючому уявленню про те, що різноманіття сприяє стійкості. Розв'язання полягає в тому, що реальні харчові мережі не є випадковими: вони мають структуру (переважають слабкі зв'язки, сильні взаємодії рідкісні), яка дозволяє великим складним мережам бути динамічно стійкими.
Епідемічна модель SIR
Компартментна модель Сприйнятливі–Інфіковані–Одужалі: dS/dt = −βSI/N, dI/dt = βSI/N − γI. Живі криві та сітка агентів. Чітко видно поріг колективного імунітету R₀ = β/γ — необхідна частка вакцинації p ≥ 1 − 1/R₀.
Фуражування мурашиної колонії
Оптимізація мурашиної колонії (ACO) на 2D-арені фуражування: мурахи залишають феромонні сліди, випаровування запобігає заклинюванню, позитивний зворотний зв'язок концентрує рух на найкоротшому шляху. Емерджентний колективний інтелект із простих локальних правил.
Потік енергії та трофічна ефективність
"Правило 10%" в екології стверджує, що приблизно 10% енергії на одному трофічному рівні доступне на наступному. Первинний продуцент вловлює сонячну радіацію; травоїдна тварина, що поїдає цю рослину, засвоює близько 10% у вигляді біомаси; хижак, що поїдає травоїдну тварину, засвоює близько 10% від того, тобто лише ~1% початкової сонячної енергії. Саме тому харчові ланцюги рідко довші за 5–6 ланок: далі просто бракує енергії для підтримки життєздатної популяції хижаків.
Причина неефективності термодинамічна: більша частина енергії їжі розсіюється у вигляді тепла під час дихання, втрачається з екскрементами або витрачається на неростові метаболічні процеси (витрати на підтримку). Екологічна ефективність = ефективність засвоєння × ефективність чистої продукції × ефективність споживання. Для теплокровних хребетних витрати на підтримку дуже високі (ендотермія коштовна), що знижує екологічну ефективність до 1–5%. Для ектотермів (комахи, риби) ефективність може сягати 15–20%.
Трофічна ефективність і потік енергії
Екологічна ефективність: ε = P_n / P_{n-1} ≈ 0.10
Продукція на рівні n: P_n = ε^n · P_0
P_0 = чиста первинна продукція (NPP, г C м⁻² рік⁻¹)
NPP помірного лісу ≈ 600 г C м⁻² рік⁻¹
NPP океану ≈ 130 г C м⁻² рік⁻¹
Доступно головному хижаку (5 рівнів): ~0.001% від NPP
Ефективність засвоєння (травоїдні): ~60–80%
Ефективність чистої продукції (ендотерми): ~1–3%