Potencjał czynnościowy
Kalmar Atlantycki Loligo pealei Potencjał czynnościowy polega na krótkotrwałej depolaryzacji
błony komórkowej. Wczesne doświadczenia (K.C. Cole i H. J.
Curtis, 1939) pokazały, że błona komórkowa staje się spolaryzowana dodatnio (ok. +50 mV) podczas maksimum potencjału czynnościowego. Gdyby powodował go jedynie chwilowy wzrost przepuszczalności dla wszystkich jonów, błona osiągnęła by 0 mV, lecz nie więcej. Obiektem do badań potencjału czynnościowego był akson Kalmara Atlantyckiego.
Potencjał czynnościowy – impuls sodowy
Zależność potencjału
czynnościowego od stężenia sodu. A i B: Maksimum potencjału
czynnościowego maleje wraz maleniem stężenia Na+ w płynie zewnątrzkomórkowym. Silna zależność wartości maksimum od stężenia Na+ wskazuje na duża przepuszczalność błony dla tych jonów w trakcie impulsu.
Alan Hodgkin i Bernard Katz odkryli, że amplituda potencjału czynnościowego zależy od koncentracji Na
+na zewnątrz komórki. Postawili hipotezę, że chwilowa zmiana przepuszczalności i wpływ jonów Na
+do wnętrza komórki powoduje
potencjał czynnościowy. Potwierdzeniem tej hipotezy była obserwacja, że
maksimum potencjału czynnościowego wynosi +55mV, co jest bliskie wartości
potencjału równowagi dla sodu. Ich eksperymenty wskazały również, że zanik
potencjału czynnościowego może być związany ze wzrostem przepuszczalności
dla jonów K
+i ich wypływem z komórki.
Eksperyment HH
A. Elektroda mierząca potencjał wewnątrzkomórkowy w aksonie (średnica 500 um). Obraz z przodu i z boku jest możliwy dzięki systemowi luster.
B. Pierwszy wewnątrzkomórkowy zapis potencjału czynnościowego. Sinusoida o częstości 500 Hz była używana jako znacznik czasu.
Hodgkin AL & Huxley AF (1939). Action potentials recorded from inside a nerve fibre. Nature 144, 710–711.
wzrost g
NaPotencjał czynnościowy – wszystko albo nic
depolaryza -cja błony
napływ Na
+‘Wybuchowa’ natura impulsu jest związana z kanałami sodowymi o przepuszczalności zależnej od
napięcia na błonie i sprzężeniem zwrotnym dodatnim z depolaryzacją błony. Repolaryzacja jest
związana z otwarciem kanałów potasowych i wypływem potasu z komórki.
Voltage clamp
Technika voltage clamp była opracowana przez Kenneth’a Cole’a w 1949 r. Alan Hodgkin i Andrew Huxley wykorzystał ją w serii eksperymentów (1952) nad mechanizmem generacji potencjału czynnościowego. Voltage clamp pozwala mierzyć wpływ zmian potencjału błonowego na przewodnictwa jonowe.
Voltage clamp działa na zasadzie ujemnego sprzężenia zwrotnego. Potencjał błonowy jest mierzony przez wzmacniacz (membrane potential amplifier) podłączony do elektrody zewnątrz i wewnątrzkomórkowej. Jest on przekazywany do wzmacniacza (feedback amplifier). Drugie wejście do wzmacniacza stanowi potencjał z generatora ustalany przez eksperymentatora (command potential). Wzmacniacz oblicza różnicę napięć i przekazuje sygnał na elektrodę biegnącą wewnątrz komórki. Prąd potrzebny do utrzymania napięcia na zadanym poziomie jest miarą prądu błonowego płynącego przez kanały jonowe.
Eksperyment voltage clamp – wstępne wyniki
Mała depolaryzacja wywołuje prąd kondensatora Ic = C dV/dt oraz leak Il.
Większa depolaryzacja wywołuje większy prąd kondensatora Ic oraz Il oraz dodatkowo prąd dokomórkowy, a następnie
odkomórkowy.
Depolaryzacja w obecności tetrodoxyny (TTX) blokującej kanały Na+, a następnie w obecności tetraethyloammonium (TEA) blokującego kanał K+, pozwala zobaczyć
‘czysty’ prąd IK i INa, po odjęciu Ic oraz Il.
I m = C m dV
dt + I K + I Na + I L I m = C m dV
dt + g K (V, t)(V - E K ) + g Na ( V, t)(V - E Na ) + g L ( V - E L )
Okładka programu wręczenia Nagrody Nobla w 1963 r.
g K ( V,t) = I K ( V,t) (V - E K ) g Na ( V, t) = I Na ( V, t)
( V - E Na )
Eksperyment voltage clamp - wyniki
Prawo Ohma (prąd jest proporcjonalny do
napięcia):
Znając I
K, I
Na, E
K, E
Na, oraz V można obliczyć g
Ki g
Na. I
K, I
Namożna otrzymać z pomiarów voltage clamp, E
K, E
Na-
stałe, V – ustala eksperymentator.
Pomiary voltage clamp dla różnych wartości V pozwoliły HH postawić hipotezę, że kanał Na
+posiada bramkę aktywacyjną i bramkę inaktywacyjną. Obie muszą być otwarte by kanał mógł przewodzić jony.
Bramka aktywacyjna jest zamknięta gdy błona znajduje się poniżej potencjału spoczynkowego i otwiera się szybko przy depolaryzacji. Bramka inaktywacyjna jest otwarta przy potencjale spoczynkowym i wolno zamyka się w wyniku depolaryzacji. Kanał K
+posiada tylko bramkę aktywacyjną otwierającą się wolno w wyniku depolaryzacji.
HH model - bramki
Zachowanie pojedynczych kanałów może być rejestrowane za pomocą patch clamp.
W zapisach widać szybkie otwieranie i zamykanie
pojedynczych kanałów. Ich suma daje gładki przebieg wartości prądu
Terminologia:
Włączenie prądu: aktywacja, wyłączenie prądu: dezaktywacja: oba procesy związane z otwarciem/zamknięciem bramki aktywacyjnej.
Gdy prąd się włącza, a następnie wyłącza pomimo stałego poziomu napięcia na błonie, oznacza to, że posiada on bramkę inaktywacją.
Procesy związane z zamknięciem/otwarciem bramki inaktywacyjnej to inaktywacja prądu/deinaktywacja prądu.
Aby mógł płynąć prąd, obie bramki musza być otwarte.
HH model - bramki
1 - y y
Otwarty Zamknięty
y - prawdopodobieństwo, że bramka jest w stanie otwartym, 1-y – że w stanie zamkniętym, , – stałe szybkości.
Uwaga: , zależą do V.
Zakładamy kinetykę reakcji pierwszego rzędu:
y dt y
dy = ( 1 - ) -
W stanie ustalonym:
-
-
=
= y y
dt
dy 0 ( 1 )
Stąd:
= + y
Podstawiając do równania:
) )(
(
) (
) (
) (
) 1
(
y y
y y
y y
dt y dy
- +
=
+ -
+
= +
-
= -
-
=
Model bramki (gate model – Hodgkin i Huxley (1952))
dy(-1) ( y
- y)
ò = (-1) ( ò + ) dt
ln( y
- y) = -( + )t+ C
Model bramki (gate model – Hodgkin i Huxley (1952)) Całkując dostajemy:
Czyli:
y= y
- Ae
-(
+
)tWidać, że y dochodzi do stanu ustalonego eksponencjalnie, ze stałą czasową , gdzie:
y
(V) = (V)
(V)+ (V)
(V) = 1
(V)+ (V)
stan ustalony
stała czasowa
Model bramki (gate model – Hodgkin i Huxley (1952))
y
0= y
- A A = y
- y
0Definiując stan poczatkowy:
y 0 = y(t= 0)
Dostajemy:
Ostatecznie:
y= y
-[(y
- y
0)e
-(
+
)t]
Lub:
y= y
0+[(y
- y
0)(1 - e
-(
+
)t)]
K K
K
K V t Y V t g n g
g ( , ) = ( , ) = 4
Na Na
Na
Na V t Y V t g m h g g ( , ) = ( , ) = 3
HH zauważyli, że g
Ki g
Nanie są funkcjami exp(-t/) lecz raczej potęgami funkcji eksponencjalnych.
Zaproponowali:
Gdzie n, m i h opisane są modelem bramki:
n n
n n
n n n
n
n n n
dt dn
= +
= + -
-
=
1
,
, )
1 (
m m
m m
m m m
m
m m m
dt dm
= +
= + -
-
=
1
,
, )
1 (
h h
h h
h h h
h
h h h
dt dh
= +
= + -
-
=
1
,
, )
1 (
HH model
Wzrost i spadek przewodnictwa K+
HH model
Przewodnictwo potasu Przewodnictwo sodu
Punkt doświadczalny Krzywa teoretyczna
g
Na= g
Nam
3[1-exp(-t/
m)]
3h
0exp(-t/
h) g
K= g
K[n
-(n
-n
0)exp(-t/
h)]
4Obserwując zmiany przewodnictwa dla różnych wartości depolaryzacji :
HH zaproponowali:
HH model
n = n
n , n = 1- n
n
m = m
m , m = 1- m
m
h = h
h , h = 1 - h
h
A następnie z funkcji
n,
m,
h, n
, m
, h
n,
n,
m,
m,
h,
hna podstawie zależności:
Z dopasowanych przebiegów g
Ki g
NaHH wyznaczyli funkcje:
,
, , ,
,
m hn
m
h
n
wyznaczyli
Następnie, do
n,
n,
m,
m,
h,
h
n(V) = 0.01(-V +10) / (e
-V+10
10
-1)
n(V) = 0.125e
-V 80
m( V) = 0.1(-V + 25) / (e
-V+25
10
-1)
m( V) = 4e
-V 18
h(V) = 0.07e
-V 20
h(V) =1/ (e
-V+30
10
-1)
HH model
dopasowali empirycznie znalezione funkcje:
Stałe szybkości aktywacji i dezaktywacji przewodnictwa K+
Stałe szybkości aktywacji i dezaktywacji przewodnictwa Na+
Stałe szybkości inaktywacji i deinaktywacji przewodnictwa Na+
Podsumowując, model Hodgkina i Huxleya jest opisany następującym ukladem rownan:
HH model
gdzie:
oraz:
I
m= C
mdV
dt + g
Kn
4( V - E
K) + g
Nam
3h(V - E
Na) + g
L(V - E
L)
HH model
HH rozwiazali numerycznie swoje równania na mechanicznym kalkulatorze i otrzymali bardzo dobrą zgodność z doświadczeniem. Model HH jest wciąż
uznawany za największy sukces w ilościowym modelowaniu mózgu a nawet i w całych naukach biologicznych.
model
eksperyment
A. Obliczony (góra) i rzeczywisty (dół) potencjał czynnościowy w aksonie kalmara.
B. Kalkulator mechaniczny (Brunsviga Maschinenwerke, Grimme, Natalis & Co) w laboratorium Alana Hodgkina w
Cambridge. Na tej maszynie Huxley rozwiązał równania HH numerycznie. Obceny w roku 1951 komupter EDSAC 1 był wyłączony z uzytkowania na 6 miesiecy. Obliczenia generacji potencjału czynnosciowego na mechanicznej maszynie liczącej trwały 3 tygodnie.
Generacja potencjału czynnościowego - podsumowanie
-100 -80 -60 -40 -20 0
0.2 0.4 0.6 0.8 1
activation gate
mV
-100 -80 -60 -40 -20 0
0.2 0.4 0.6 0.8 1
inactivation gate
mV
-100 -80 -60 -40 -20 0
5 10
15 tau activation
mV
ms
-100 -80 -60 -40 -20 0
100 200 300
tau recovery and inactivation
mV
ms