Modèle de Hodgkin–Huxley

Le modèle Hodgkin-Huxley est un modèle mathématique qui décrit comment les potentiels d'action des neurones sont initiés et propagés. Il s'agit d'un ensemble d'équations différentielles non linéaires qui se rapprochent des caractéristiques électriques de cellules excitables telles que les neurones ou les myocytes cardiaques. C'est un système dynamique en temps continu.

Composants de base des modèles de type Hodgkin-Huxley. Les modèles de type Hodgkin-Huxley représentent la caractéristique biophysique des membranes cellulaires. La bicouche lipidique est représentée par une capacité (C _m). Les canaux ioniques voltage-dépendants et de fuite sont représentés par des conductances non linéaires (g _n) et linéaires (g _L), respectivement. Les gradients électrochimiques entraînant le flux d'ions sont représentés par des batteries (E), et les pompes et échangeurs d'ions sont représentés par des sources de courant (I _p).

Originellement le modèle fut décrit en 1952 par Alan Hodgkin et Andrew Huxley pour expliquer les mécanismes ioniques sous-tendant l'initiation et la propagation des potentiels d'action observé dans l'axone géant du calmar. Le prix Nobel de physiologie ou médecine de 1963 leur fut attribué pour ce travail.

Alors que les expériences originales ne traitaient que des canaux sodium et potassium, le modèle Hodgkin-Huxley peut également être étendu pour tenir compte d'autres espèces de canaux ioniques.

Composants de base

Le modèle de Hodgkin-Huxley modélise l'influence des canaux ionique de potassium et sodium sur le potentiel membranaire.

Le modèle de Hodgkin-Huxley traite chaque composant d'une cellule excitable comme un élément électrique (comme le montre la figure). La bicouche lipidique est représentée par une capacité $C_{m}$. Les canaux ioniques voltage-dépendants sont représentés par les conductances $g_{n}$, soumises aux gradients électrochimiques entraînant le flux d'ions $E_{n}$, (où n est le canal ionique spécifique) fonction de la tension et du temps. Les canaux de fuite sont représentés par des conductances linéaires $g_{L}$. Les valeurs de $E_{n}$dont les valeurs sont déterminées par le rapport entre les concentrations intra- et extracellulaires des espèces ioniques d'intérêt. Enfin, les pompes à ions sont représentées par des sources de courant ($I_{p}$). Le potentiel membranaire est noté $I_{c}$.

Ainsi, le courant traversant la bicouche lipidique est donné par:

$I_{c}=C_{m}{\frac {{\mathrm {d} }V_{m}}{{\mathrm {d} }t}}$

et le courant traversant un canal ionique est donné par :

$I_{i}={g_{i}}(V_{m}-V_{i})\;$

où $V_{i}$ est le potentiel d'inversion du i-ième canal ionique. Ainsi, pour une cellule possédant des canaux sodium et potassium, le courant total traversant la membrane est donné par :

$I=C_{m}{\frac {{\mathrm {d} }V_{m}}{{\mathrm {d} }t}}+{\bar {g}}_{\text{K}}n^{4}(V_{m}-V_{K})+{\bar {g}}_{\text{Na}}m^{3}h(V_{m}-V_{Na})+{\bar {g}}_{l}(V_{m}-V_{l})$,

où I est le courant membranaire total par unité de surface, $C_{m}$ est la capacité membranaire par unité de surface, $g_{K}$ et $g_{Na}$ sont les conductances de potassium et de sodium par unité de surface, respectivement, $V_{K}$ et $V_{Na}$ sont les potentiels d'inversion du potassium et du sodium, respectivement, et $g_{l}$et $V_{l}$ sont respectivement la conductance de fuite par unité de surface et le potentiel d'inversion de fuite. Les éléments dépendant du temps de cette équation sont $V_{m}$, $g_{Na}$ et $g_{K}$, où les deux dernières conductances dépendent également explicitement de la tension.

Caractérisation du courant ionique

Dans les canaux ioniques voltage-dépendants, la conductance du canal $g_{l}$ est fonction à la fois du temps et de la tension ($g_{n}(t,V)$ dans la figure), tandis que dans les canaux de fuite $g_{l}$ est une constante ($g_{L}$ dans la figure). Le courant généré par les pompes ioniques dépend des espèces ioniques spécifiques à cette pompe. Les sections suivantes décriront ces formulations plus en détail.

Canaux ioniques voltage-dépendants

À l'aide d'une série d'expériences de voltage imposé et en faisant varier les concentrations extracellulaires de sodium et de potassium, Hodgkin et Huxley ont développé un modèle dans lequel les propriétés d'une cellule excitable sont décrites par un ensemble de quatre équations différentielles ordinaires. Ce modèle est développé originellement pour l'axones de calmar qui ne contiennent que des canaux de sodium et potassium mais peut être étendu à d'autres axones possédant plus de canaux. Ainsi en plus l'équation du courant total donné ci-dessus, on a aussi :

$I_{\mathrm {Na} }(t)={\bar {g}}_{\mathrm {Na} }m(V_{m})^{3}h(V_{m})(V_{m}-E_{\mathrm {Na} })$,

${\frac {dn}{dt}}=\alpha _{n}(V_{m})(1-n)-\beta _{n}(V_{m})n$

${\frac {dm}{dt}}=\alpha _{m}(V_{m})(1-m)-\beta _{m}(V_{m})m$

${\frac {dh}{dt}}=\alpha _{h}(V_{m})(1-h)-\beta _{h}(V_{m})h$

où I est le courant par unité de surface, et $\alpha _{i}$ et $\beta _{i}$ sont des constantes de vitesse pour le i- ième canal ionique, qui dépendent de la tension mais pas du temps. ${\bar {g}}_{n}$ est la valeur maximale de la conductance. n, m et h sont des quantités sans dimension entre 0 et 1 qui sont respectivement associées à l'activation des canaux potassiques, à l'activation des canaux sodiques et à l'inactivation des canaux sodiques. Pour $p=(n,m,h)$, $\alpha _{p}$ et $\beta _{p}$ prendre la forme :

$\alpha _{p}(V_{m})=p_{\infty }(V_{m})/\tau _{p}$,

$\beta _{p}(V_{m})=(1-p_{\infty }(V_{m}))/\tau _{p}$.

$p_{\infty }$ et $(1-p_{\infty })$ sont les valeurs d'état stable pour l'activation et l'inactivation, respectivement, et sont généralement représentées par les équations de Boltzmann en tant que fonctions de $V_{m}$. Dans l'article original de Hodgkin et Huxley, les fonctions $\alpha$ et $\beta$ sont donnés par :

${\begin{array}{lll}\alpha _{n}(V_{m})={\frac {0.01(10+V)}{\exp {\big (}{\frac {10+V}{10}}{\big )}-1}}&\alpha _{m}(V_{m})={\frac {0.1(25+V)}{\exp {\big (}{\frac {25+V}{10}}{\big )}-1}}&\alpha _{h}(V_{m})=0.07\exp {\bigg (}{\frac {V}{20}}{\bigg )}\\\beta _{n}(V_{m})=0.125\exp {\bigg (}{\frac {V}{80}}{\bigg )}&\beta _{m}(V_{m})=4\exp {\bigg (}{\frac {V}{18}}{\bigg )}&\beta _{h}(V_{m})={\frac {1}{\exp {\big (}{\frac {30+V}{10}}{\big )}+1}}\end{array}}$

où $V=V_{rest}-V_{m}$ désigne la dépolarisation négative en mV.

Alors que dans de nombreux logiciels actuels, les modèles de type Hodgkin-Huxley généralisent $\alpha$ et $\beta$ à :

${\frac {A_{p}(V_{m}-B_{p})}{\exp {\big (}{\frac {V_{m}-B_{p}}{C_{p}}}{\big )}-D_{p}}}$

Afin de caractériser les canaux voltage-dépendants, les équations sont ajustées aux données obtenu sous une tension d'entrée constante. Pour une dérivation des équations de Hodgkin-Huxley sous tension d'entrée constante, voir. En bref, pour un potentiel membranaire maintenu à une valeur constante les équations de déclenchement non linéaires se réduisent aux équations de la forme :

$m(t)=m_{0}-[(m_{0}-m_{\infty })(1-e^{-t/\tau _{m}})]$,

$h(t)=h_{0}-[(h_{0}-h_{\infty })(1-e^{-t/\tau _{h}})]$,

$n(t)=n_{0}-[(n_{0}-n_{\infty })(1-e^{-t/\tau _{n}})]$.

Ainsi, pour chaque valeur de potentiel de membrane $V_{m}$ les courants de sodium et de potassium peuvent être décrits par :

$I_{\mathrm {K} }(t)={\bar {g}}_{\mathrm {K} }n(V_{m})^{4}(V_{m}-E_{\mathrm {K} })$.

Afin d'arriver à la solution complète pour un potentiel d'action propagé, on écrit le courant I du côté gauche de la première équation différentielle en fonction de V, de sorte que l'équation soit pour la tension seule. La relation entre I et V peut être dérivée de la théorie des câbles et est donnée par :

$I={\frac {a}{2R}}{\frac {\partial ^{2}V}{\partial x^{2}}}$,

où a est le rayon de l'axone, R est la résistance spécifique de l'axoplasme et x est la position le long de la fibre nerveuse.

L'algorithme de Levenberg-Marquardt est souvent utilisé pour ajuster ces équations aux données de voltage-clamp.

Canaux de fuite

Les canaux de fuite tiennent compte de la perméabilité naturelle de la membrane aux ions et prennent la forme de l'équation des canaux voltage-dépendants, où la conductance $g_{leak}$ est une constante. Ainsi, le courant de fuite dû aux canaux ioniques de fuite passifs dans le formalisme de Hodgkin-Huxley est $I_{l}=g_{leak}(V-V_{leak})$.

Pompes et échangeurs

Le potentiel membranaire dépend du maintien des gradients de concentration ionique à travers elle. Le maintien de ces gradients de concentration nécessite un transport actif d'espèces ioniques. Les échangeurs sodium-potassium et sodium-calcium sont les plus connus d'entre eux. Certaines des propriétés de base de l'échangeur Na/Ca sont déjà bien établies : la stoechiométrie d'échange est de 3 Na ⁺ : 1 Ca ²⁺ et l'échangeur est électrogène et sensible à la tension. L'échangeur Na/K a également été décrit en détail, avec une stoechiométrie 3 Na ⁺ : 2 K ⁺ .

Propriétés mathématiques

Le modèle Hodgkin-Huxley peut être considéré comme un système d'équations différentielles avec quatre variables d'état, $V_{m}(t),n(t),m(t)$, et $h(t)$, qui changent par rapport au temps $t$. Le système est difficile à étudier car il s'agit d'un système non linéaire et ne peut pas être résolu analytiquement. Cependant, il existe de nombreuses méthodes numériques disponibles pour analyser le système. Certaines propriétés et comportements généraux, tels que les cycles limites, peuvent être prouvés.

Une simulation du modèle Hodgkin-Huxley dans l'espace des phases, en termes de tension v(t) et de variable de déclenchement de potassium n(t). La courbe fermée est connue sous le nom de cycle limite.

Collecteur central

Comme il existe quatre variables d'état, il peut être difficile de visualiser le chemin dans l'espace des phases. Habituellement, deux variables sont choisies, la tension $V_{m}(t)$ et la variable de contrôle du potassium $n(t)$, permettant de visualiser le cycle limite. Cependant, il faut être prudent car il s'agit d'une méthode ad hoc de visualisation du système à 4 dimensions. Cela ne prouve pas l'existence du cycle limite.

Une meilleure projection peut être construite à partir d'une analyse minutieuse du Jacobien du système, évalué au point d'équilibre. Plus précisément, les valeurs propres du Jacobien sont indicatives de l'existence de la variété centrale. De même, les vecteurs propres du Jacobien révèlent l'orientation de la variété centrale. Le modèle Hodgkin-Huxley a deux valeurs propres négatives et deux valeurs propres complexes avec des parties réelles légèrement positives. Les vecteurs propres associés aux deux valeurs propres négatives se réduiront à zéro à mesure que le temps t augmente. Les deux vecteurs propres complexes restants définissent la variété centrale. En d'autres termes, le système à 4 dimensions s'effondre sur un plan à 2 dimensions. Toute solution partant du collecteur central se désintégrera vers le collecteur central. De plus, le cycle limite est contenu sur le collecteur central.

La tension v (t) (en millivolts) du modèle Hodgkin-Huxley, représentée graphiquement sur 50 millisecondes. Le courant injecté varie de -5 nanoampères à 12 nanoampères. Le graphique passe par trois étapes : une étape d'équilibre, une étape à pic unique et une étape de cycle limite.

Bifurcations

Si le courant injecté $I$ ont été utilisés comme paramètre de bifurcation, alors le modèle Hodgkin-Huxley subit une bifurcation de Hopf. Comme avec la plupart des modèles neuronaux, l'augmentation du courant injecté augmentera le taux de décharge du neurone. Une conséquence de la bifurcation de Hopf est qu'il y a une cadence de tir minimale. Cela signifie que soit le neurone ne se déclenche pas du tout (correspondant à une fréquence nulle), soit se déclenche à la vitesse de déclenchement minimale. En raison du principe tout ou rien, il n'y a pas d'augmentation régulière de l'amplitude du potentiel d'action, mais plutôt un "saut" soudain d'amplitude. La transition qui en résulte est connue sous le nom de canard.

Améliorations et modèles alternatifs

Le modèle Hodgkin-Huxley est considéré comme l'une des grandes réalisations de la biophysique du XXe siècle. Néanmoins, les modèles modernes de type Hodgkin-Huxley ont été étendus de plusieurs manières importantes :

• Des populations de canaux ioniques supplémentaires ont été incorporées sur la base de données expérimentales.
• Le modèle Hodgkin-Huxley a été modifié pour incorporer la théorie des états de transition et produire des modèles thermodynamiques Hodgkin-Huxley.
• Les modèles incorporent souvent des géométries très complexes de dendrites et d'axones, souvent basées sur des données de microscopie.
• Modèles stochastiques du comportement des canaux ioniques, conduisant à des systèmes hybrides stochastiques.
• Le modèle Poisson-Nernst-Planck (PNP) est basé sur une approximation du champ moyen des interactions ioniques et des descriptions de continuum de concentration et de potentiel électrostatique.

Plusieurs modèles neuronaux simplifiés ont également été développés (comme le modèle FitzHugh-Nagumo), facilitant une simulation efficace à grande échelle de groupes de neurones, ainsi qu'un aperçu mathématique de la dynamique de génération de potentiel d'action.

Notes et références

(en) Cet article est partiellement ou en totalité issu de l’article de Wikipédia en anglais intitulé « Hodgkin–Huxley model » (voir la liste des auteurs).

Voir aussi

Articles connexes

• Potentiel d'action
• Modèles du neurone biologique

Bibliographie

• A. L. Hodgkin et A. F. Huxley, « Currents carried by sodium and potassium ions through the membrane of the giant axon of Loligo », The Journal of Physiology, vol. 116, n^o 4,‎ 28 avril 1952, p. 449–472 (ISSN 0022-3751, PMID 14946713, PMCID 1392213, lire en ligne, consulté le 4 octobre 2021).
• A. L. Hodgkin et A. F. Huxley, « The components of membrane conductance in the giant axon of Loligo », The Journal of Physiology, vol. 116, n^o 4,‎ 28 avril 1952, p. 473–496 (ISSN 0022-3751, PMID 14946714, PMCID 1392209, lire en ligne, consulté le 4 octobre 2021).
• A. L. Hodgkin et A. F. Huxley, « The components of membrane conductance in the giant axon of Loligo », The Journal of Physiology, vol. 116, n^o 4,‎ 28 avril 1952, p. 473–496 (ISSN 0022-3751, PMID 14946714, PMCID 1392209, lire en ligne, consulté le 4 octobre 2021).
• A. L. Hodgkin et A. F. Huxley, « A quantitative description of membrane current and its application to conduction and excitation in nerve », The Journal of Physiology, vol. 117, n^o 4,‎ 28 août 1952, p. 500–544 (ISSN 0022-3751, PMID 12991237, PMCID 1392413, lire en ligne, consulté le 4 octobre 2021).
• A. L. Hodgkin, A. F. Huxley et B. Katz, « Measurement of current-voltage relations in the membrane of the giant axon of Loligo », The Journal of Physiology, vol. 116, n^o 4,‎ 28 avril 1952, p. 424–448 (ISSN 0022-3751, PMID 14946712, PMCID 1392219, lire en ligne, consulté le 4 octobre 2021).

Liens externes

• Simulation Javascript interactive du modèle HH — Fonctionne dans n'importe quel navigateur compatible HTML5. Permet de modifier les paramètres du modèle et l'injection de courant.
• Applet Java interactif du modèle HH — Les paramètres du modèle peuvent être modifiés ainsi que les paramètres d'excitation et les tracés de l'espace des phases de toutes les variables sont possibles.
• Lien direct vers le modèle Hodgkin-Huxley et une description dans la base de données BioModels.
• Impulsions neuronales : le potentiel d'action en action par Garrett Neske, The Wolfram Demonstrations Project.
• Modèle interactif Hodgkin-Huxley par Shimon Marom, The Wolfram Demonstrations Project.
• ModelDB — Une base de données de code source en neurosciences computationnelles contenant 4 versions (dans différents simulateurs) du modèle Hodgkin-Huxley original et des centaines de modèles qui appliquent le modèle Hodgkin-Huxley à d'autres canaux dans de nombreux types de cellules électriquement excitables.
• Plusieurs articles sur la version stochastique du modèle et son lien avec l'original.