Нейроваскулярне з’єднання та сигнали БОЛД

Активні нейрони потребують доставки крові

Нейрони мають майже безперервний запас палива — вони покладаються на постійне постачання глюкози та кисню з крові. Коли ділянка кортикси стає активною, її енергетичні потреби зростають протягом мілісекунд, і відбувається координаційний сигнал, який залучає більше крові до цієї конкретної ділянки протягом декількох секунд: поруч розташовані артеріоли розширюються, локальний кровотік збільшується та об’єм крові в цій області зростає. Ця ланцюгова реакція — від збудженого нейрона до розширеного кровопостачання — є нейровіскулярним зв’язком, і вона опосередкована «нейровіскулярною одиницею» нейронів, астроцитів та стінки судної посудини, з оксидом азоту як одним із ключових вазодилататорних сигналів, що виділяються в рамках відповіді.

Чому сигнал перевищує межі

Ось неінтуїтивне пояснення, яке робить можливим fMRI. Активні нейрони споживають більше кисню, що само по собі збільшує локальну дезоксигемоглобін — а дезоксигемоглобін є парамагнітним і спотворює локальне магнітне поле, темніючи сигнал МРТ. Однак підвищена кровопостачання, викликана нейровіскулярною взаємозв’язністю, доставляє значно більшу кількість киснем насиченої крові, ніж це використовує тканина. В результаті концентрація дезоксигемоглобіну в активній області фактично падає, і сигнал МРТ у цій області стає яскравішим. Це підвищення яскравості, залежне від рівня кисню (BOLD), яке вимірюють майже всі дослідження на основі завдань.

Модель повітряної кулі

Стандартна кількісна модель цього процесу, представлена Бауксоном та колегами у 1998 році та розширена Фрістоном та колегами до широко використовуваної моделі «мішень-Віндкесел», розглядає місцеву венозну камеру як еластичну кулю, наповнену припливом і злиту відтоком, що залежить нелінійно від того, наскільки повністю заповнена куля:

dv/dt = ( f_in(t) - f_out(v) ) / τ // об’єм крові кулі f_out(v) = v^(1/α) // відтік зростає з повнотою dq/dt = ( f_in(t) * E(f_in)/E0 - f_out(v) * q/v ) / τ // деоксигемоглобін BOLD(t) ≈ V0 * ( k1*(1 - q) + k2*(1 - q/v) + k3*(1 - v) ) Об’єм v збільшується, коли приплив f_in перевищує відтік, і відтік зростає швидше, ніж лінійно, коли куля вже досить повна — властивість, запозичена з моделі Віндкеселя артеріальної відповідності. Деоксигемоглобін q підвищується з вилученням кисню, але вимивається відтоком, і останій сигнал BOLD є зваженою комбінацією того, наскільки змінився об’єм крові та вміст деоксигемоглобіну відносно базового рівня, із ваговими коефіцієнтами (k1, k2, k3), які залежать від потужності сканера та послідовності імпульсів.

dv/dt = ( f_in(t) - f_out(v) ) / tau         // blood volume balloon

f_out(v) = v^(1/alpha)                        // outflow grows with fullness

dq/dt = ( f_in(t) * E(f_in)/E0 - f_out(v) * q/v ) / tau   // deoxyhemoglobin

BOLD(t) ≈ V0 * ( k1*(1 - q) + k2*(1 - q/v) + k3*(1 - v) )

Чому сигнал відстає та перевищує

Оскільки модель містить дві закорельовані змінні стану (об’єм і деоксигемоглобін) з різними часовими константами, BOLD-відповідь на короткий імпульс нейронної активності є повільним, плавним підвищенням — рулонкова гемодинамічна функція — яке досягає піку приблизно через 4–6 секунд після нейронного події і не повертається до початкового рівня протягом 15–20 секунд, що значно перевищує мілісекундний часовий масштаб основних спалахів. Також добре задокументований пост-стимуляційний відкат — сигнал тимчасово падає нижче базового рівня після основного піку — випливає безпосередньо з моделі: об’єм крові повертається до базового рівня повільніше, ніж кровотік, тому на кілька секунд деоксигемоглобін підвищений порівняно з все ще роздутим кулька, тимчасово затемнюючи сигнал перед тим, як все стабілізується.

This lag is a fundamental consequence of the model’s simplified representation of vascular dynamics and neuronal activity. The slow relaxation of blood volume contributes to the overall delay in the BOLD response.

Що це означає для інтерпретації карти фМРТ

Оскільки BOLD є гемодинамічним проксі, а не прямим електричним вимірюванням, його часовий роздільний здатність фундаментально обмежений багатосекундним часом вазогенного відповіді, незалежно від того, наскільки швидко сканер може семплювати. Це також означає, що амплітуда BOLD відображає сукупність нейронної активності та вазогенної реактивності, які можуть відрізнятися між мозковими областями, людьми та станами здоров’я — це застереження, яке необхідно враховувати кожній фМРТ-дослідженню при порівнянні інтенсивності сигналу в різних місцях або популяціях.

Frequently asked questions

Спробуйте наживо

Усе, що вище, працює прямо у вашому браузері — відкрийте Neurovascular Coupling і змінюйте параметри під час роботи. Нічого не встановлюється, нічого не завантажується на сервер, уся модель живе в одній вкладці.