Modelo de potencial de accin de Hodking & Huxley

Modelo Matemtico HodgkinHuxley

Reporte Final
Bogot, Colombia

Abstract: Con este modelo se describir cmo

se inicia y se transmite el potencial de accin a
travs de las neuronas. Para esto, es necesario
profundizar en el comportamiento de las
corrientes de los iones de potasio, sodio y
cloro (ion para la corriente de fuga) a partir de
la variacin del voltaje que se genera alrededor
de la membrana del axn durante el
experimento.
Palabras claves: Bifurcacin de Hopf
I. INTRODUCCIN
El modelo de Hodkin-Huxley (HH) es la
representacin matemtica del proceso, mediante el
cual se estudia el comportamiento del transporte del
potencial de accin largo de una neurona. Durante
este proceso se realizan intercambios progresivos
de iones en forma de corrientes entre el interior y el
exterior de la membrana a travs de los canales que
se abren para permitir su paso.
Durante el experimento de Hodkin y Huxley se
identificaron tres corrientes de iones diferentes, la de
potasio, la de sodio, y una corriente de fuga
compuesta en su mayora por iones de cloro. La
conduccin de stas corrientes no se produce
instantneamente ya que la membrana opone una
resistencia que disminuye con el aumento del
voltaje, es decir, que la conductancia de sta
aumenta como el comportamiento de la carga de
un condensador.
Las caractersticas descritas anteriormente llevaron
a que Hodkin y Huxley desarrollaran un modelo de
cuatro dimensiones con variables de estado como el
voltaje a travs de la membrana, la probabilidad de
que se abran los canales de sodio (m,h) y de
potasio (n).
II. MODELO DE HODKIN-HUXLEY (HH)
El punto de inicio del modelo de HH est
representado por un circuito elctrico equivalente

En el cual se encuentran diferentes tipos de

corrientes inicas como sodio (INa), potasio Ik y una
corriente que Hodking y Huxley denominaron
corriente de fuga, representada por IL.
Cm e Im estn relacionadas con el comportamiento
como tal de la membrana.
A su vez, el modelo est representado por 4
ecuaciones no lineales que Hodking- Huxley
emplearon para representar su modelo.

En donde, y son los coeficientes que

representan las partculas de compuertas que
determinan la velocidad en la que se genera el
cambio de las partculas de los canales de iones
del estado abierto al estado cerrado, las cuales
estn definidas, para cada uno de los iones, por
las siguientes ecuaciones.
Partculas variables de activacin y desactivacin
de las compuertas de sodio:

h=
e

30v
10

h=0,07 e

v
20

+1

m=4 e 18 m=

25v
v

10e 80 1
Partculas variables de activacin y desactivacin
de las compuertas de potasio:

Modelo de potencial de accin de Hodking & Huxley

v

n=0,125 e 80 n=

10v
100e

10v
10

Adems, para evaluar el modelo se debe de tener

en cuenta las siguientes constantes:

potasio que representan a las partculas de

compuerta y para ste modelo, estn dadas
por las siguientes constantes.

Y a su vez, maneja los mismos valores para las

constantes VNa, Vk, VL, Na, K y L anteriormente
mencionados en el modelo no simplificado.
Para este modelo realizaron un procedimiento
matemtico, donde, mediante la comprobacin de
los teoremas, se demostr que existe una
bifurcacin de Hopf. A raz de sta, se pudo hallar
que existen tres puntos crticos para cada uno de
los parmetros (I,VE).

Para las cuales Cm es la capacitancia que se

produce a travs de la membrana de la neurona,
VNa, Vk y VL son los voltajes de equilibrio (No existe
intercambio de iones entre el interior y el exterior
de la membrana) del sodio, el potasio y otros
iones (conformados mayormente por iones de
cloro), respectivamente, adems tenemos las
conductancias mximas a travs de los canales
de sodio y potasio y la conductancia mxima de
fuga estn representadas por Na, K y L,
respectivamente.

Si tomamos el parmetro VE=0 para evaluar la

respuesta del sistema solo dependiente de I,
obtenemos las siguientes graficas:

III. MODELO SIMPLIFICADO DE HODKINGHUXLEY

Como podemos inferir, el modelo HH posee un
grado de complejidad elevado as que se ha
buscado realizar una simplificacin a dos
dimensiones, el cual est descrito por las siguientes
dos ecuaciones

Donde V es el voltaje a travs de la membrana,

I es la corriente de entrada a la neurona, V E el
voltaje del campo elctrico que se est generando
entre el interior y el exterior de la membrana, n la
probabilidad de que se abra un canal de potasio,
1, 2, 1 y 2 que son los coeficientes que
determinan la velocidad en la que se genera el
cambio de las partculas de los canales de iones
del estado abierto al estado cerrado.
En el cual, a diferencia del modelo no simplificado,
las contantes halladas por medio de una
linealizacin a las tasas de cambio de estado
abierto a cerrado de las partculas del canal de

Figura. 1
Como podemos ver en la Fig.1 (a) y (b), que
representa la respuesta para 0<I<2.4082A/cm2, el
sistema no gener un potencial de accin, pues de
acuerdo a la grfica el voltaje slo genera una
oscilacin antes de estabilizarse en un valor fijo, al
igual
que
en
(c)
y
(d),

Modelo de potencial de accin de Hodking & Huxley

2.4083A/cm2<I<5.2647A/cm2, se generan un
nmero limitado de picos antes de tender a una
constante. Si se sigue aumentando el valor de I se
obtienen las siguientes respuestas:

*Imagen ampliada al final del documento [Anexo 1].

Se realiz las grficas variando el valor de la
corriente para poder observar la similitud del
modelo simplificado con el modelo sin simplificar.
El cdigo que se muestra a continuacin fue
empleado en cada una de la elaboracin de las
grficas variando el valor de la corriente que se
indica en cada caso.

Grfica Sim1,
I=2.4283uA

para

una

corriente

%Graficas I=2.4283uA
Figura. 2
Para la Fig.2 se incrementaron los valores de I,
como podemos ver para el intervalo de
5.2648A/cm2<I<454A/cm2
obtenemos
la
respuesta deseada de la generacin peridica
estable de picos de potencial de accin como se
aprecia en (a) y (b), pero cuando I aumenta mucho
(454A/cm2<I) el sistema se vuelve inestable.
No se enfatizar el procedimiento mediante el cual
se estudi la respuesta del sistema a la variable
VE con I=0, ya que para el modelo sin simplificar,
el cual se trabajar, maneja como nica entrada
del sistema, la corriente y no es tenido en cuenta
el voltaje del campo elctrico que se est
generando entre el interior y el exterior de la
membrana, el cual est representado por esta
variable (VE).

plot(V0)
title('Respuesta del sistema')
xlabel('Tiempo(ms)')
ylabel('V(t)')
grid on
figure
plot(V0.data,n.data)
title('n(t) Vs V(t)')
xlabel('V(t)')
ylabel('n(t)')
grid on

IV. SIMULACIN DE ECUACIONES

Mediante el uso de la herramienta de MatlabSimulink, se desarroll el diagrama de bloques que
comprende al modelo de 4 dimensiones
relacionadas con el modelo de HH.

Figuras Sim1

Grfica Sim2, para una corriente I=5.25

uA

Modelo de potencial de accin de Hodking & Huxley

trminos de v, ya se puede hallar los puntos de
equilibrio para cada una de stas, los cuales
corresponden a los siguientes datos:

Figuras Sim2

Grfica Sim3,
I=5.2648 uA

para

una

corriente

m=0.0532
h=0.5945
n=0.3183
V=0.462
stos sern los puntos de equilibrio los cuales se
reemplazaran en los trminos derivados para
finalmente hallar la matriz linealizada

Figuras Sim3

Grfica Sim4, para una corriente I=180

uA

Matriz Linealizada.
La ecuacin de salida con respecto a la variable
de estado V.

Figuras Sim4

Grfica Sim5,
I=600uA

para

una

corriente

Figuras Sim5
Como se observa en las simulaciones realizadas
con el modelo no simplificado, las cuales guardan
gran similitud con el los resultados obtenidos con
el modelo simplificado, se puede inferir que la
validez de esta simplificacin es de gran
importancia ya que reduce la complejidad del
modelo, sin embargo, ste empieza a generar
gran discrepancia con el modelo original de HH
cuando la corriente empieza a aumentar.
V. MODELO LINEALIZADO
Una vez teniendo las 4 ecuaciones que describen
el modelo, habiendo derivado cada una de las
componentes que componen la matriz A y hallando
las ecuaciones correspondientes para m, n y h en

Published with MATLAB R2013b

VI. COMPARACIN Y ANLISIS DE
RESULTADOS

Hodking y huxley realizaron diferentes anlisis

para lograr un completo entendimiento del
potencial de accin que se generaba alrededor de
la membrana del axn del calamar. Como ya se ha
mencionado, para el anlisis del modelo sin
simplificar no se tuvo en cuenta el comportamiento
que el campo elctrico, V E, generaba en sta, se
analizar tomando VE=0 y variando los intervalos
de la corriente.

Modelo de potencial de accin de Hodking & Huxley

Para el modelo simplificado, el primer intervalo a
estudiar fue cuando la corriente I=2,4082 A y para
el caso de la simulacin realizada en simulink el
primer intervalo fue de I=2,4283 A. Observando
las grficas correspondientes a estos valores, las
cuales son las imgenes a) y b) de la figura 1 y la
figura sim1, podemos observar su similitud, en
las primeras grficas, en que no alcanza a haber
un pulso para luego la seal mantenerse
constante. Para el caso de las fases de stas
mismas, la fase de la imagen b) es coherente con
la grfica a) ya que en sta se puede ver que no se
alcanza a completar este pulso, a diferencia de la
fase para Sim1, la cual es una grficas cerrada,
entendindose como pulso completo.
Las siguientes grficas a analizar son las
correspondientes a las imgenes a) y b) de la
figura 2 las cuales corresponden a una corriente
I=5,2648 A y la imagen Sim3 desarrollada a una
corriente I=5,2648 A las cuales concuerdan
completamente ya que en este punto las respuesta
del sistema es oscilatoria estable ya que se
generan pulsos peridicamente que en ningn
momento tienden a un valor infinito y el diagrama
de fase tiende despus oscilar a un punto fijo, lo
cual es muy coherente ya que los valores de
corriente tomados para analizar cada grfica, son
exactamente iguales.
Para las terceras grficas a comparar las cuales
corresponden a las imgenes e) y f) dela figura 1
comparndolas con la imagen nombrada como
sim2, las cuales manejan una corriente I=5,2647
A, para el primer caso e I=5,25 A, para el
segundo caso, en las cuales se generan una
cantidad finita de oscilaciones antes de volver a un
estado de equilibrio.
Y para finalizar, se comparan, con corrientes
iguales a I=600A para el primer caso relacionado
a la imagen sim5, e I=455A para las imgenes
e) y f) de la figura 2 en las cuales, debido al gran
aumento en la magnitud de la corriente el sistema
se vuelve inestable y no genera respuesta
apreciable del sistema.
Cuando se tiene un campo elctrico excitado a
corrientes muy pequeas, las fibras nerviosas no
generan potencial de accin, por lo tanto, el voltaje
en la membrana llega a equilibrio despus de
pequeas oscilaciones. En cambio, cuando se

aumenta un poco la corriente suceden tres

cambios notables: el nmero de parmetro
aumenta, la amplitud del voltaje de la membrana
disminuye y aumenta la frecuencia de oscilacin
de la membrana.

VII. CONCLUSIONES

Cuando se quiere realizar el muestreo y

visualizar bien los pulsos que se generan
al observar el comportamiento de la
transferencia de potencial de accin a
travs de las neuronas, se deben tomar
valores de corrientes dentro de un rango
apropiado (donde la respuesta del
sistema oscila y es estable) o si no, se
visualizar una seal constante lo cual no
permitir tomar valores y conclusiones
muy precisas.

El modelo simplificado de Hodking y

Huxley que emplea 2 dimensiones,
hallado a partir de un proceso matemtico
y una linealizacin, guarda gran relacin
con el modelo de cuatro dimensiones,
debido a que las respuestas obtenidas
con cada uno de estos modelos, a un
mismo estmulo, son muy parecidas.

Este modelo presentado de cuatro

dimensiones de Hodking y Huxley, ofrece
un cuadro de trabajo para las personas
interesadas en la excitabilidad de las
neuronas cuando se le inyecta una
corriente al axn, en este caso, del
calamar.

La complejidad de las cuatro ecuaciones

no lineales hace que existan diferentes
tipos de resultados, por lo cual en
ocasiones es mejor tener en cuenta el
modelo simplificado de dos dimensiones
de HH.

Este tipo de modelos experimentales es

muy til y es una tcnica muy popular
para modelar los diferentes tipos de
canales existentes.

Estos modelos son aplicables al estudio

del sistema nervioso ya que realizan una

Modelo de potencial de accin de Hodking & Huxley

aproximacin
muy
realista
al
comportamiento de la neurona segn
diferentes estmulos.

VIII. Referencias
Sterrat, D., Graham, B., Gillies, A., &
Willshaw,
D. (2011). Principles of Computational
Modelling in Neuroscience. Cambridge

Wang, H., Yu, Y., Wang, M., & Zhao, R. (2010). The
Dynamical Analysis of a Two-dimensional
Simplified Hodgkin-Huxley Model. Sixth
International Conference on Natural
Computation, (pgs. 814-818).