Przekaźnictwo dopaminergiczne (DA) odgrywa kluczową rolę w różnych funkcjach poznawczych i behawioralnych, a jego dysfunkcja jest powiązana z wieloma powszechnymi chorobami i zaburzeniami. W związku z tym krytyczne jest opracowanie dokładnych metod ilościowego badania przekaźnictwa DA in vivo, aby ocenić, w jaki sposób jest ono zmienione w kontekście modeli chorób i farmakologii leków. Szybka woltametria cykliczna (FSCV) umożliwia monitorowanie przekaźnictwa DA in vivo z wysoką rozdzielczością przestrzenną i czasową. Choć możliwe jest monitorowanie fizjologicznego przekaźnictwa DA u zwierząt budzących się i swobodnie zachowujących, stymulacja elektryczna wstępujących dróg dopaminergicznych u zwierząt znieczulonych może wywołać silne odpowiedzi DA, które pozwalają na pogłębioną analizę kinetyczną przekaźnictwa DA.
Odpowiedzi DA indukowane stymulacją elektryczną odzwierciedlają dynamiczną interakcję między uwalnianiem a wychwytem DA, a interpretacje tych odpowiedzi opierały się głównie na prostym modelu stymulowanej neurotransmisji DA, zwanym modelem Michaelisa-Menten (M-M).12Model M-M składa się z 3 zmiennych opisujących odpowiedzi DA w zależności od stałej szybkości uwalniania DA oraz stałej wydajności wychwytu zwrotnego (t. j. zależność między szybkością wychwytu zwrotnego DA a pozakomórkowymi stężeniami DA), opisaną przez Równanie 1:

(uwalnianie DA) (wychwyt zwrotny DA)
W równaniu 1 f oznacza częstotliwość stymulacji; [DA]p to szacowany wzrost stężenia DA na jeden impuls stymulacji; Vmax reprezentuje szacowaną maksymalną szybkość wychwytu; a Km jest szacowaną stałą M-M, która teoretycznie odpowiada pozakomórkowemu stężeniu DA nasycającemu 50% transporterów DAT, co prowadzi do połowy maksymalnej szybkości wychwytu. To równanie różniczkowe można scalkulować, aby symulować eksperymentalne odpowiedzi DA poprzez oszacowanie parametrów [DA]p, Vmax oraz Km.
Chociaż model M-M umożliwił znaczące postępy w rozumieniu kinetyki neurotransmisji DA w różnych kontekstach eksperymentalnych, opiera się on na uproszczonych podstawowych założeniach, które ograniczają jego zastosowanie przy modelowaniu odpowiedzi DA wywołanych stymulacją suprafizjologiczną2,13. Na przykład model M-M może jedynie przybliżać kształty odpowiedzi DA, jeśli rosną one w sposób wypukły, ale nie potrafi on wyjaśnić stopniowych (wklęsłych) odpowiedzi wzrostowych stwierdzanych w obszarach grzbietowych prążkowia12. Zatem założenia modelu M-M nie oddają dokładnie dynamicznych procesów uwalniania i wychwytu zwrotnego podczas stymulowanej neurotransmisji DA.
Aby zamodelować stymulowane odpowiedzi DA zgodnie z realistycznymi ramami ilościowymi, opracowano ilościowy model neurobiologiczny (QN) oparty na zasadach kinetyki stymulowanej neurotransmisji wynikających z uzupełniających badań i eksperymentów2. Różne nurty badań nad neurotransmisją wykazują, że (1) stymulowane uwalnianie neuroprzekaźnika jest procesem dynamicznym, którego tempo spada w czasie trwania stymulacji14, (2) uwalnianie trwa w fazie postymulacyjnej z dwufazową kinetyką zaniku15 oraz (3) wydajność wychwytu zwrotnego DA jest stopniowo hamowana w trakcie trwania samej stymulacji2,16. Te trzy koncepcje stanowią fundament modelu QN oraz trzech równań składających się z 12 parametrów opisujących dynamikę uwalniania i wychwytu zwrotnego DA (Tabela 1). Model QN pozwala na precyzyjną symulację heterogenicznych typów eksperymentalnych odpowiedzi DA, a także przewidywanych efektów eksperymentalnych manipulacji parametrami stymulacji i podawania leków2,6. Chociaż dalsze badania są niezbędne do udoskonalenia podejścia do modelowania danych, przyszłe eksperymenty mogą w znacznym stopniu skorzystać z tego ugruntowanego neurobiologicznie podejścia do modelowania, które znacząco rozszerza wnioski wyciągane z paradygmatu stymulowanej neurotransmisji DA.

Tabela 1: Równania i parametry modelowania. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.
Niniejszy samouczek opisuje sposób modelowania danych odpowiedzi DA po stymulacji w celu oszacowania kinetyki uwalniania i ponownego wychwytu DA przy użyciu programu QNsim 1.0. Rzeczywiste gromadzenie i przetwarzanie danych eksperymentalnych nie zostało tutaj opisane; wymagane są jedynie temporalne dane o stężeniu DA. Podstawy teoretyczne i fundamenty struktury QN zostały szeroko opisane wcześniej2, natomiast poniżej przedstawiono praktyczne podejście do zastosowania struktury QN w modelowaniu danych odpowiedzi DA.
Model QN odwzorowuje dynamiczne oddziaływanie pomiędzy: 1) dynamicznym uwalnianiem DA, 2) wychwytem zwrotnym DA oraz 3) wpływem stymulacji suprafizjologicznych na te procesy, aby wyodrębnić istotne informacje kinetyczne z danych dotyczących odpowiedzi DA. Model QN najlepiej nadaje się do modelowania danych FSCV pozyskanych przy użyciu silnych stymulacji suprafizjologicznych o długim czasie trwania (np. stymulacje 60 Hz przez 10 s), które generują silne odpowiedzi DA podatne na analizę kinetyczną. Po dokładnym zamodelowaniu podstawowych procesów uwalniania i wychwytu zwrotnego, parametry modelu mogą zostać wykorzystane do symulacji odpowiedzi DA, która powinna odpowiadać kształtem eksperymentalnej odpowiedzi DA.
Równania modelu QN opisują tempo uwalniania i wychwytu DA w czasie stymulowanych odpowiedzi DA. Model QN opisuje tempo stymulowanego uwalniania DA jako funkcję czasu od rozpoczęcia stymulacji (tstim), podczas której tempo uwalniania DA maleje wykładniczo w trakcie stymulacji. Jest to zgodne z wyczerpywaniem się puli łatwo uwalnialnej, z dodatkowym stałym tempem uwalniania DA (DARss) w celu uwzględnienia uzupełniania pęcherzyków, co jest spójne z innymi doniesieniami (Równanie 2)14,17.

Manipulacje zwiększające tempo uwalniania DA, takie jak zwiększenie ΔDAR, ΔDARτlub DARss, prowadzą do zwiększenia amplitud odpowiedzi na wykresach DA w funkcji czasu. Każdy parametr w różny sposób wpływa na kształt odpowiedzi DA. Zwiększenie DARss i ΔDARτ oba sprawiają, że faza wzrostowa odpowiedzi jest bardziej liniowa (mniej wypukła). Zmniejszenie ΔDARτ sprzyja wypukłości, która jest kontrolowana przez wartość ΔDAR. Na podstawie doświadczeń z modelowania, DARss jest zazwyczaj mniejsze niż 1/5th z ΔDAR; zatem ΔDAR jest parametrem uwalniania, który w głównym stopniu determinuje całkowitą amplitudę odpowiedzi DA.
Szybkość uwalniania DA po stymulacji jest modelowana za pomocą Równania 3 jako kontynuacja stymulowanej szybkości uwalniania DA z momentu zakończenia stymulacji (DARES) w funkcji czasu po stymulacji (tpost). Szybkość uwalniania DA po stymulacji wykazuje dwufazowy wzorzec zaniku, jak opisano wcześniej15, z szybką fazą zaniku wykładniczego i przedłużoną fazą zaniku liniowego, aby wymodelować dwa zależne od wapnia procesy uwalniania neuroprzekaźnika.

(Szybki zanik wykładniczy) (Przedłużony zanik liniowy)
Obecnie nie jest możliwe określenie ilości uwalnianej DA po stymulacji. Ograniczenie to można rozwiązać poprzez systematyczną minimalizację szacunków uwalniania DA po stymulacji oraz walidację parametrów modelu na podstawie zestawu eksperymentalnych odpowiedzi DA, zebranych z tego samego miejsca rejestracji przy zastosowaniu różnych czasów trwania stymulacji. Taka minimalizacja pozwala użytkownikom na dokonywanie konserwatywnych szacunków uwalniania i ponownego wychwytu. Ponieważ stymulacje elektryczne prowadzą do akumulacji wapnia, która sprzyja uwalnianiu neuroprzekaźnika po stymulacji, czas trwania stymulacji wpływa na parametry uwalniania neuroprzekaźnika po stymulacji18,19. Na podstawie doświadczeń w modelowaniu stwierdzono, że wraz ze wzrostem czasu trwania stymulacji wartość τR rośnie, a XR maleje, co jest zgodne z przewidywanymi efektami większej akumulacji wapnia20.
Równanie 4 opisuje szybkość wychwytu zwrotnego DA jako rozszerzenie modelu M-M i wprowadza dynamiczny parametr Km, który zwiększa się podczas stymulacji, aby odwzorować stopniowy spadek wydajności wychwytu zwrotnego spowodowany stymulacjami suprafizjologicznymi2,16. Wartość Km po stymulacji pozostaje stała i odpowiada wartości Km osiągniętej pod koniec stymulacji (KmES).

gdzie

(Podczas stymulacji) (Po stymulacji)
Stymulowane odpowiedzi DA, zwłaszcza w obszarach brzusznych prążkowia, są często niewrażliwe na zmiany początkowej wartości Km (Kmi), co utrudnia zdefiniowanie wartości Kmi. Dlatego, podobnie jak w oryginalnym modelu M-M, Kmi przyjmuje się w przedziale od 0,1 do 0,4 µM dla odpowiedzi DA zarejestrowanych u kontrolnych, nieleczonych zwierząt12. Termin ΔKm określa zakres zmiany wydajności wychwytu podczas stymulacji, który w naszym doświadczeniu wynosi około 20 µM w trakcie 10-sekundowej stymulacji o częstotliwości 60 Hz. Wartości k i Kminf określają, w jaki sposób Km zmienia się w czasie, a zwiększenie któregokolwiek z tych parametrów zwiększa wklęsłość fazy wzrostowej. Vmax to maksymalna szybkość wychwytu, która w części zależy od lokalnej gęstości transporterów DA, wykazującej gradient od strony brzuszzno-przyśrodkowej do grzbietowo-bocznej21. W związku z tym wartości Vmax w grzbietowym prążkowiu (D-Str) są zazwyczaj wyższe niż 30 µM/s, natomiast zazwyczaj niższe niż 30 µM/s w obszarach brzusznych, takich jak jądro półleżące (NAc)6.
Powyższe ogólne wytyczne mogą pomóc w modelowaniu eksperymentalnych danych dotyczących odpowiedzi DA, jednak stworzenie symulacji przybliżającej eksperymentalną odpowiedź DA wymaga iteracyjnej regulacji parametrów modelu. Dokładność parametrów modelu można poprawić poprzez uzyskanie odpowiedzi DA na stymulacje suprafizjologiczne, które zapewniają solidną podstawę do symulacji, a także poprzez uzyskanie i modelowanie wielokrotnych odpowiedzi DA na stymulacje o różnym czasie trwania w tym samym miejscu rejestracji (np. stymulacje 60-Hz o czasie trwania 5 s i 10 s) w celu walidacji dokładności parametrów (patrz: przykładowe dane). W ramach demonstracji do pakietu oprogramowania dołączono zestaw danych zawierający regioselektywne stymulowane odpowiedzi DA zebrane w jądrze półleżącym i grzbiecie prążkowia, przed i po wyzwaniu farmakologicznym, które zostało już zamodelowane przy użyciu struktury QN. W konsekwencji użytkownicy stwierdzą, że tę metodologię można w podobny sposób zastosować do charakterystyki kinetyki neurotransmisji DA w różnych kontekstach chorobowych i manipulacjach farmakologicznych.