Modelos cuantitativos del potencial de acción.

En neurofisiología se han desarrollado varios modelos matemáticos del potencial de acción , que se dividen en dos tipos básicos. El primer tipo busca modelar cuantitativamente los datos experimentales, es decir, reproducir exactamente las mediciones de corriente y voltaje. El renombrado modelo de Hodgkin-Huxley del axón del calamar Loligo ejemplifica tales modelos. ^[1] Aunque cualitativamente correcto, el modelo HH no describe con precisión cada tipo de membrana excitable, ya que considera solo dos iones (sodio y potasio), cada uno con un solo tipo de canal sensible al voltaje. Sin embargo, otros iones como el calcio pueden ser importantes y existe una gran diversidad de canales para todos los iones. ^[2] Como ejemplo, el potencial de acción cardíaco ilustra cómo se pueden generar potenciales de acción de diferentes formas en membranas con canales de calcio sensibles al voltaje y diferentes tipos de canales de sodio/potasio. El segundo tipo de modelo matemático es una simplificación del primer tipo; El objetivo no es reproducir los datos experimentales, sino comprender cualitativamente el papel de los potenciales de acción en los circuitos neuronales. Para tal fin, los modelos fisiológicos detallados pueden ser innecesariamente complicados y oscurecer el "bosque de los árboles". El modelo de FitzHugh-Nagumo es típico de esta clase, que a menudo se estudia por su comportamiento de arrastre . ^[3] El arrastre se observa comúnmente en la naturaleza, por ejemplo en la iluminación sincronizada de las luciérnagas , que está coordinada por una explosión de potenciales de acción; ^[4] el arrastre también se puede observar en neuronas individuales. ^[5] Ambos tipos de modelos pueden usarse para comprender el comportamiento de pequeñas redes neuronales biológicas , como los generadores de patrones centrales responsables de algunas acciones reflejas automáticas. ^[6] Tales redes pueden generar un patrón temporal complejo de potenciales de acción que se utiliza para coordinar las contracciones musculares, como las involucradas en la respiración o la natación rápida para escapar de un depredador. ^[7]

Modelo de Hodgkin-Huxley

Circuito eléctrico equivalente para el modelo de potencial de acción de Hodgkin-Huxley. I _m y V _m representan la corriente y el voltaje a través de un pequeño parche de membrana, respectivamente. La C _m representa la capacitancia del parche de membrana, mientras que las cuatro g representan las conductancias de cuatro tipos de iones. Las dos conductancias de la izquierda, para el potasio (K) y el sodio (Na), se muestran con flechas para indicar que pueden variar con el voltaje aplicado, correspondientes a los canales iónicos sensibles al voltaje .

En 1952, Alan Lloyd Hodgkin y Andrew Huxley desarrollaron un conjunto de ecuaciones para ajustar sus datos experimentales de fijación de voltaje en la membrana axonal. ^{[1] [8]} El modelo supone que la capacitancia de la membrana C es constante; por lo tanto, el voltaje transmembrana V cambia con la corriente transmembrana total I _tot según la ecuación

${\displaystyle C{\frac {dV}{dt}}=I_{\mathrm {tot} }=I_{\mathrm {ext} }+I_{\mathrm {Na} }+I_{\mathrm {K} }+I_{\mathrm {L} }}$

donde I _Na , I _K e IL son corrientes transportadas a través de los canales locales de sodio, los canales de potasio y los canales _de "fuga" (un comodín), respectivamente. El término inicial I _ext representa la corriente que llega de fuentes externas, como los potenciales postsinápticos excitadores de las dendritas o el electrodo de un científico.

El modelo supone además que un canal iónico determinado está completamente abierto o cerrado; si está cerrado, su conductancia es cero, mientras que si está abierto, su conductancia es un valor constante g . Por lo tanto, la corriente neta a través de un canal iónico depende de dos variables: la probabilidad p de que el canal esté _abierto y la diferencia de voltaje con respecto al voltaje de equilibrio de ese ión, V − V _eq . Por ejemplo, la corriente a través del canal de potasio se puede escribir como

${\displaystyle I_{\mathrm {K} }=g_{\mathrm {K} }\left(V-E_{\mathrm {K} }\right)p_{\mathrm {open,K} }}$

que es equivalente a la ley de Ohm . Por definición, no fluye corriente neta ( I _K = 0) cuando el voltaje transmembrana es igual al voltaje de equilibrio de ese ion (cuando V = E _K ).

Para ajustar sus datos con precisión, Hodgkin y Huxley supusieron que cada tipo de canal iónico tenía múltiples "puertas", de modo que el canal estaba abierto sólo si todas las puertas estaban abiertas y, en caso contrario, se cerraban. También asumieron que la probabilidad de que una puerta estuviera abierta era independiente de que las otras puertas estuvieran abiertas; Esta suposición fue posteriormente validada para la puerta de inactivación. ^[9] Hodgkin y Huxley modelaron el canal de potasio sensible al voltaje con cuatro puertas; Si p _n denota la probabilidad de que una sola puerta esté abierta, la probabilidad de que todo el canal esté abierto es el producto de cuatro de dichas probabilidades, es decir, p _{abierta, K} = n ⁴ . De manera similar, se modeló la probabilidad del canal de sodio sensible al voltaje para que tuviera tres puertas similares de probabilidad my una cuarta puerta, asociada con la inactivación, de probabilidad h ; por tanto, p _{abierto, Na} = m ³ h . Se supone que las probabilidades de cada puerta obedecen a una cinética de primer orden.

${\displaystyle {\frac {dm}{dt}}=-{\frac {m-m_{\mathrm {eq} }}{\tau _{m}}}}$

donde tanto el valor de equilibrio m _eq como la constante de tiempo de relajación τ _m dependen del voltaje instantáneo V a través de la membrana. Si V cambia en una escala de tiempo más lentamente que τ _m , la probabilidad m siempre será aproximadamente igual a su valor de equilibrio m _eq ; sin embargo, si V cambia más rápidamente, entonces m quedará rezagado respecto de m _eq . Al ajustar sus datos de fijación de voltaje, Hodgkin y Huxley pudieron modelar cómo estos valores de equilibrio y constantes de tiempo variaban con la temperatura y el voltaje transmembrana. ^[1] Las fórmulas son complejas y dependen exponencialmente del voltaje y la temperatura. Por ejemplo, la constante de tiempo para la probabilidad de activación del canal de sodio h varía como 3 ^{(θ−6.3)/10} con la temperatura Celsius θ y con el voltaje V como

${\displaystyle {\frac {1}{\tau _{h}}}=0.07e^{-V/20}+{\frac {1}{1+e^{3-V/10}}}.}$

En resumen, las ecuaciones de Hodgkin-Huxley son ecuaciones diferenciales ordinarias complejas y no lineales en cuatro variables independientes : el voltaje transmembrana V y las probabilidades m , h y n . ^[10] No se ha descubierto ninguna solución general de estas ecuaciones. Un método menos ambicioso pero generalmente aplicable para estudiar tales sistemas dinámicos no lineales es considerar su comportamiento en las proximidades de un punto fijo . ^[11] Este análisis muestra que el sistema Hodgkin-Huxley sufre una transición desde la inactividad estable hasta las oscilaciones explosivas a medida que la corriente estimulante I _ext aumenta gradualmente; Sorprendentemente, el axón vuelve a estar establemente inactivo a medida que la corriente estimulante aumenta aún más. ^[12] También se ha llevado a cabo un estudio más general de los tipos de comportamiento cualitativo de los axones predichos por las ecuaciones de Hodgkin-Huxley. ^[10]

Modelo FitzHugh-Nagumo

Figura FHN: Para imitar el potencial de acción, el modelo de FitzHugh-Nagumo y sus parientes utilizan una función g ( V ) con resistencia diferencial negativa (una pendiente negativa en la gráfica I vs. V ). A modo de comparación, una resistencia normal tendría una pendiente positiva, según la ley de Ohm I = GV , donde la conductancia G es la inversa de la resistencia G =1/ R .

Debido a la complejidad de las ecuaciones de Hodgkin-Huxley, se han desarrollado varias simplificaciones que exhiben un comportamiento cualitativamente similar. ^{[3] [13]} El modelo de FitzHugh-Nagumo es un ejemplo típico de un sistema tan simplificado. ^{[14] [15]} Basado en el diodo túnel , el modelo FHN tiene solo dos variables independientes, pero exhibe un comportamiento de estabilidad similar al de las ecuaciones completas de Hodgkin-Huxley. ^[16] Las ecuaciones son

${\displaystyle C{\frac {dV}{dt}}=I-g(V),}$

${\displaystyle L{\frac {dI}{dt}}=E-V-RI}$

donde g(V) es función del voltaje V que tiene una región de pendiente negativa en el medio, flanqueada por un máximo y un mínimo (Figura FHN). Un caso simple muy estudiado del modelo de FitzHugh-Nagumo es el modelo nervioso de Bonhoeffer-van der Pol, que se describe mediante las ecuaciones ^[17]

${\displaystyle C{\frac {dV}{dt}}=I-\epsilon \left({\frac {V^{3}}{3}}-V\right),}$

${\displaystyle L{\frac {dI}{dt}}=-V}$

donde se supone que el coeficiente ε es pequeño. Estas ecuaciones se pueden combinar en una ecuación diferencial de segundo orden.

${\displaystyle C{\frac {d^{2}V}{dt^{2}}}+\epsilon \left(V^{2}-1\right){\frac {dV}{dt}}+{\frac {V}{L}}=0.}$

Esta ecuación de van der Pol ha estimulado muchas investigaciones en las matemáticas de sistemas dinámicos no lineales . Keener ha desarrollado circuitos de amplificador operacional que realizan los modelos FHN y van der Pol del potencial de acción. ^[18]

Morris y Lecar desarrollaron un híbrido de los modelos Hodgkin-Huxley y FitzHugh-Nagumo en 1981 y lo aplicaron a la fibra muscular de los percebes . ^[19] Fiel a la fisiología del percebe, el modelo de Morris-Lecar reemplaza la corriente de sodio dependiente de voltaje del modelo de Hodgkin-Huxley con una corriente de calcio dependiente de voltaje. No hay inactivación (no hay variable h ) y la corriente de calcio se equilibra instantáneamente, de modo que nuevamente, solo hay dos variables dependientes del tiempo: el voltaje transmembrana V y la probabilidad de puerta de potasio n . Se han estudiado en detalle la explosión, el arrastre y otras propiedades matemáticas de este modelo. ^[20]

Los modelos más simples del potencial de acción son los modelos de "vaciar y llenar" (también llamados modelos de "integrar y disparar"), en los que la señal de entrada se suma (la fase de "llenado") hasta que alcanza un umbral, disparando una pulso y restablecer la suma a cero (la fase de "limpieza"). ^{[3] [21] [22]} Todos estos modelos son capaces de exhibir arrastre , que se observa comúnmente en los sistemas nerviosos. ^[3]

Potenciales y corrientes extracelulares.

Mientras que los modelos anteriores simulan el voltaje y la corriente transmembrana en un solo parche de membrana, otros modelos matemáticos se refieren a los voltajes y corrientes en la solución iónica que rodea la neurona. ^[23] Estos modelos son útiles para interpretar datos de electrodos extracelulares, que eran comunes antes de la invención del electrodo de pipeta de vidrio que permitía el registro intracelular. ^[24] El medio extracelular puede modelarse como una solución iónica isotrópica normal ; En tales soluciones, la corriente sigue las líneas del campo eléctrico , según la forma continua de la ley de Ohm.

${\displaystyle \mathbf {j} =\sigma \mathbf {E} }$

donde j y E son vectores que representan la densidad de corriente y el campo eléctrico , respectivamente, y donde σ es la conductividad . Por lo tanto, j se puede encontrar a partir de E , que a su vez se puede encontrar usando las ecuaciones de Maxwell . Las ecuaciones de Maxwell pueden reducirse a un problema relativamente simple de electrostática , ya que las concentraciones iónicas cambian demasiado lentamente (en comparación con la velocidad de la luz ) para que los efectos magnéticos sean importantes. El potencial eléctrico φ( x ) en cualquier punto extracelular x se puede resolver utilizando las identidades de Green ^[23]

${\displaystyle \phi (\mathbf {x} )={\frac {1}{4\pi \sigma _{\mathrm {outside} }}}\oint _{\mathrm {membrane} }{\frac {\partial }{\partial n}}{\frac {1}{\left|\mathbf {x} -{\boldsymbol {\xi }}\right|}}\left[\sigma _{\mathrm {outside} }\phi _{\mathrm {outside} }({\boldsymbol {\xi }})-\sigma _{\mathrm {inside} }\phi _{\mathrm {inside} }({\boldsymbol {\xi }})\right]dS}$

donde la integración es sobre toda la superficie de la membrana; es una posición en la membrana, σ _dentro y φ _dentro son la conductividad y el potencial justo dentro de la membrana, y σ _fuera y φ _fuera los valores correspondientes justo fuera de la membrana. Por lo tanto, dados estos valores de σ y φ en la membrana, el potencial extracelular φ ( x ) se puede calcular para cualquier posición x ; a su vez, el campo eléctrico E y la densidad de corriente j se pueden calcular a partir de este campo potencial. ^[25] ${\displaystyle {\boldsymbol {\xi }}}$

Ver también

• Modelos de neuronas biológicas.
• Ecuación actual de GHK
• Modelos de computación neuronal.
• Conducción saltatoria
• Bioelectrónica
• Teoría del cable

Referencias

1. Hodgkin AL , Huxley AF , Katz B (1952). "Corrientes transportadas por iones de sodio y potasio a través de la membrana del axón gigante de Loligo". Revista de fisiología . 116 (4): 424–448. doi : 10.1113/jphysiol.1952.sp004717. PMC  1392213 . PMID  14946713.{{cite journal}}: Mantenimiento CS1: varios nombres: lista de autores ( enlace )
   Hodgkin AL , Huxley AF (1952). "Corrientes transportadas por iones de sodio y potasio a través de la membrana del axón gigante de Loligo". Revista de fisiología . 116 (4): 449–472. doi : 10.1113/jphysiol.1952.sp004717. PMC  1392213 . PMID  14946713.
   Hodgkin AL , Huxley AF (1952). "Los componentes de la conductancia de membrana en el axón gigante de Loligo". J Physiol . 116 (4): 473–496. doi : 10.1113/jphysiol.1952.sp004718. PMC  1392209 . PMID  14946714.
   Hodgkin AL , Huxley AF (1952). "El doble efecto del potencial de membrana sobre la conductancia del sodio en el axón gigante de Loligo". J Physiol . 116 (4): 497–506. doi : 10.1113/jphysiol.1952.sp004719. PMC  1392212 . PMID  14946715.
   Hodgkin AL , Huxley AF (1952). "Una descripción cuantitativa de la corriente de membrana y su aplicación a la conducción y excitación en los nervios". J Physiol . 117 (4): 500–544. doi : 10.1113/jphysiol.1952.sp004764. PMC  1392413 . PMID  12991237.
2. Yamada WM, Koch C , Adams PR (1989). "Múltiples canales y dinámica del calcio". En C. Koch , I Segev (ed.). Métodos de modelado neuronal: de las sinapsis a las redes . Cambridge, Massachusetts: Libro de Bradford, The MIT Press. págs. 97-133. ISBN 978-0-262-11133-1.{{cite book}}: Mantenimiento CS1: varios nombres: lista de autores ( enlace )
3. Hoppensteadt FC (1986). Una introducción a las matemáticas de las neuronas . Cambridge: Prensa de la Universidad de Cambridge. ISBN 978-0-521-31574-6.
4. Hanson, FE; Caso, JF; Dólar, E.; Dólar, J. (1971). "Sincronía y arrastre de destellos en una luciérnaga de Nueva Guinea". Ciencia . 174 (4005): 161–164. Código Bib : 1971 Ciencia... 174.. 161H. doi :10.1126/ciencia.174.4005.161. PMID  17742039. S2CID  41760422.
5. Guttman R, Feldman L, Jacobsson E (1980). "Arrastre de frecuencia del axón de calamar". J. Miembro. Biol . 56 (1): 9–18. doi :10.1007/BF01869347. PMID  7441721. S2CID  10775478.
6. Obtener PA (1989). "Reconstrucción de pequeñas redes neuronales". En C Koch y I Segev (ed.). Métodos de modelado neuronal: de las sinapsis a las redes . Cambridge, Massachusetts: Libro de Bradford, The MIT Press. págs. 171-194. ISBN 978-0-262-11133-1.
7. Hooper, Scott L. "Generadores de patrones centrales". Embryonic ELS (1999) http://www.els.net/elsonline/figpage/I0000206.html ^{[ enlace muerto permanente ]} (2 de 2) [6/2/2001 11:42:28 a. m.] En línea: consultado el 27 de noviembre 2007 [1].
8. Nelson ME, Rinzel J (1994). "El modelo Hodgkin-Huxley" (PDF) . En Bower J, Beeman D (eds.). El libro de GÉNESIS: Exploración de modelos neuronales realistas con el sistema de simulación neuronal general . Nueva York: Springer Verlag. págs. 29–49.
9. Armstrong CM, Bezanilla F, Rojas E (1973). "Destrucción de la inactivación de la conductancia del sodio en axones de calamar perfundidos con pronasa". J. Gen. Physiol . 62 (4): 375–391. doi :10.1085/jgp.62.4.375. PMC 2226121 . PMID  4755846. 
   Rojas E, Rudy B (1976). "Destrucción de la inactivación de la conductancia del sodio por una proteasa específica en fibras nerviosas perfundidas de Loligo". J. Physiol . 262 (2): 501–531. doi :10.1113/jphysiol.1976.sp011608. PMC  1307656 . PMID  994046.
10. Sato S, Fukai H, Nomura T, Doi S (2005). "Análisis de bifurcación de las ecuaciones de Hodgkin-Huxley". En Reeke GN, Poznanski RR, Lindsay KA, Rosenberg JR, Sporns O (eds.). Modelado en neurociencias: de los sistemas biológicos a la robótica neuromimética (2ª ed.). Boca Ratón: CRC Press. págs. 459–478. ISBN 978-0-415-32868-5.
11. Guckenheimer J , Holmes P (1986). Oscilaciones no lineales, sistemas dinámicos y bifurcaciones de campos vectoriales (segunda impresión, edición revisada y corregida). Nueva York: Springer Verlag. págs. 12-16. ISBN 978-0-387-90819-9.
12. Sabah NH, Spangler RA (1970). "Respuesta repetitiva del modelo de Hodgkin-Huxley para el axón gigante del calamar". Revista de Biología Teórica . 29 (2): 155-171. Código Bib : 1970JThBi..29..155S. doi :10.1016/0022-5193(70)90017-2. PMID  5500466.
    Evans JW (1972). "Ecuaciones de axones nerviosos. I. Aproximaciones lineales". Universidad de Indiana. Matemáticas. J.​ 21 (9): 877–885. doi : 10.1512/iumj.1972.21.21071 .
    Evans JW, Feroe J (1977). "Teoría de la estabilidad local del impulso nervioso". Matemáticas. Biosci . 37 (1–2): 23–50. doi :10.1016/0025-5564(77)90076-1.
13. FitzHugh R (1960). "Umbrales y mesetas en las ecuaciones nerviosas de Hodgkin-Huxley". J. Gen. Physiol . 43 (5): 867–896. doi :10.1085/jgp.43.5.867. PMC 2195039 . PMID  13823315. 
    Kepler TB, Abbott LF, Marder E (1992). "Reducción de modelos neuronales basados ​​en conductancia". Cibernética biológica . 66 (5): 381–387. doi :10.1007/BF00197717. PMID  1562643. S2CID  6789007.
14. FitzHugh R (1961). "Impulsos y estados fisiológicos en modelos teóricos de membrana nerviosa". Revista Biofísica . 1 (6): 445–466. Código bibliográfico : 1961BpJ......1..445F. doi :10.1016/S0006-3495(61)86902-6. PMC 1366333 . PMID  19431309. 
15. Nagumo J, Arimoto S, Yoshizawa S (1962). "Una línea de transmisión de pulso activa que simula el axón nervioso". Actas del IRE . 50 (10): 2061-2070. doi :10.1109/JRPROC.1962.288235. S2CID  51648050.
16. FitzHugh R (1969). "Modelos matemáticos de excitación y propagación en nervio". En HP Schwann (ed.). Ingeniería biológica . Nueva York: McGraw-Hill. págs. 1–85.
17. Bonhoeffer KF (1948). "ACTIVACIÓN DEL HIERRO PASIVO COMO MODELO DE EXCITACIÓN DEL NERVIO". J. Gen. Physiol . 32 (1): 69–91. doi :10.1085/jgp.32.1.69. PMC 2213747 . PMID  18885679. 
    Bonhoeffer KF (1953). "Modelle der Nervenerregung". Naturwissenschaften . 40 (11): 301–311. Código Bib : 1953NW..... 40.. 301B. doi :10.1007/BF00632438. S2CID  19149460.
    van der Pol B (1926). "Sobre las oscilaciones de relajación". Revista Filosófica . 2 : 978–992.
    van der Pol B , van der Mark J (1928). "El latido del corazón considerado como una oscilación de relajación y un modelo eléctrico del corazón". Revista Filosófica . 6 : 763–775.
    van der Pol B , van der Mark J (1929). "El latido del corazón considerado como una oscilación de relajación y un modelo eléctrico del corazón". Arco. Neerl. Fisiol . 14 : 418–443.
18. Keener JP (1983). "Circuitos analógicos para las ecuaciones de van der Pol y FitzHugh-Nagumo". Transacciones IEEE sobre sistemas, hombre y cibernética . 13 (5): 1010–1014. doi :10.1109/TSMC.1983.6313098. S2CID  20077648.
19. Morris C, Lecar H (1981). "Oscilaciones de voltaje en la fibra muscular gigante de percebe". Revista Biofísica . 35 (1): 193–213. Código Bib : 1981BpJ....35..193M. doi :10.1016/S0006-3495(81)84782-0. PMC 1327511 . PMID  7260316. 
20. Rinzel J, Ermentrout GB (1989). "Análisis de oscilaciones y excitabilidad neuronal". En C. Koch , I Segev (ed.). Métodos de modelado neuronal: de las sinapsis a las redes . Cambridge, Massachusetts: Libro de Bradford, The MIT Press. págs. 135-169. ISBN 978-0-262-11133-1.
21. Keener JP, Hoppensteadt FC, Rinzel J (1981). "Modelos de integración y activación de respuesta de la membrana nerviosa a la entrada oscilatoria". Revista SIAM de Matemática Aplicada . 41 (3): 503–517. doi :10.1137/0141042.
22. Cejnar, Pavel; Vyšata, Oldřich; Kukal, Jaromír; Beránek, Martín; Vališ, Martín; Procházka, Aleš (2020). "El modelo simple de interruptor de condensador de neurona excitadora e inhibidora con todas las partes explicadas biológicamente permite oscilaciones caóticas dependientes del patrón de disparo de entrada". Informes científicos . 10 (1): 7353. Código bibliográfico : 2020NatSR..10.7353C. doi : 10.1038/s41598-020-63834-7 . PMC 7192907 . PMID  32355185. 
23. Stevens, Charles F. (1966). Neurofisiología: introducción . Nueva York: John Wiley and Sons. págs. 161-173. LCCN  66015872.
24. Ling G, Gerard RW (1949). "El potencial de membrana normal de las fibras de rana sartorio". J. Celda. comp. Fisiol . 34 (3): 383–396. doi :10.1002/jcp.1030340304. PMID  15410483.
25. Lorente de N°R (1947). "Un estudio de fisiología nerviosa". Semental. Instituto Rockefeller. Medicina. Investigación . 132 : Cap. dieciséis.
    Mauro A (1960). "Propiedades de generadores delgados pertenecientes a potenciales electrofisiológicos en conductores de volumen". J. Neurofisiol . 23 (2): 132-143. doi :10.1152/jn.1960.23.2.132.
    Woodbury JW (1965). "Capítulo 3: Potenciales en un conductor de volumen". En TC Ruch; HD Patton (eds.). Fisiología y Biofísica . Filadelfia: WB Saunders Co.

Otras lecturas

• Vidrio L, Mackey MC (1988). De los relojes al caos: los ritmos de la vida . Princeton, Nueva Jersey: Universidad de Princeton. ISBN 978-0-691-08496-1.