Neurones et électricité : La modélisation du potentiel d’action par un dipôle RC

Problématique

En quoi la membrane phospholipidique d’un neurone peut-elle être modélisée par un dipôle électrique RC, et comment l’établissement des équations différentielles de charge et de décharge décrit-il la propagation de l’influx nerveux ?

---

Introduction

Bonjour à toutes et à tous. Pour que nous puissions penser, bouger ou ressentir nos émotions, notre cerveau doit transmettre des milliards d’informations chaque seconde. Cette communication s’effectue le long de cellules hautement spécialisées : les neurones. L’information y circule sous forme de signaux électriques appelés potentiels d’action (ou influx nerveux). Mais comment une cellule vivante, composée d’eau et de molécules organiques, parvient-elle à générer, stocker et propager un courant électrique ?

Pour répondre à cette question, les biologistes et les physiciens ont uni leurs forces. Ils ont découvert que la membrane d’un axone se comporte exactement comme un circuit électrique simple de notre programme de Physique de Terminale : le dipôle RC.

Ayant choisi les spécialités Physique-Chimie et Sciences de la Vie et de la Terre, j’ai trouvé passionnant d’étudier ce croisement entre physique des circuits et neurobiologie. Ma problématique est la suivante : En quoi la membrane phospholipidique d’un neurone peut-elle être modélisée par un dipôle électrique RC, et comment l’établissement des équations différentielles de charge et de décharge décrit-il la propagation de l’influx nerveux ?

Pour y répondre, nous analyserons d’abord l’analogie structurelle entre la membrane cellulaire et les composants d’un dipôle RC. Puis, nous étudierons la modélisation dynamique des variations de tension à l’aide d’équations différentielles.

---

I. L’analogie biologique et électrique : La membrane comme dipôle RC (environ 4 min)

Pour comprendre comment l’électricité circule dans un neurone, examinons la structure de la membrane de son prolongement principal : l’axone.

La membrane cellulaire est constituée d’une bicouche phospholipidique extrêmement fine (d’environ $4$ à $5\text{ nanomètres}$ d’épaisseur). Cette bicouche est hydrophobe et se comporte comme un isolant électrique. Elle sépare deux milieux aqueux riches en ions ($Na^+$, $K^+$, $Cl^-$) qui sont, eux, d’excellents conducteurs électriques : le milieu intracellulaire (cytoplasme) et le milieu extracellulaire.

En physique, un composant constitué de deux plaques conductrices (les armatures) séparées par un isolant (le diélectrique) est un condensateur. La membrane du neurone joue exactement ce rôle. Elle accumule des charges positives sur sa face externe et des charges négatives sur sa face interne, stockant ainsi de l’énergie sous forme de potentiel électrique. La capacité électrique de cette membrane, notée $C$, dépend de sa surface et de son épaisseur.

De plus, cette membrane n’est pas totalement hermétique. Elle est parsemée de protéines spécialisées formant des canaux ioniques (les canaux de fuite et les canaux activés par la tension). Ces canaux permettent aux ions de traverser la membrane sous l’effet du gradient électrochimique. Ce passage d’ions s’oppose à la force de frottement du canal, ce qui équivaut physiquement à une résistance électrique, notée $R$ (ou à une conductance $G = 1/R$).

À l’état de repos, la différence de potentiel transmembranaire (tension membranaire $V_m = V_{\text{int}} - V_{\text{ext}}$) est stable et vaut environ $-70\text{ mV}$. Cette tension est maintenue grâce à l’activité de pompes ioniques actives ($Na^+/K^+$ ATP-ase) qui consomment de l’énergie pour rejeter le sodium et faire entrer le potassium.

---

II. Modélisation dynamique par les équations différentielles (environ 3 min 30s)

Le potentiel d’action est une variation transitoire et stéréotypée de cette tension de repos, divisée en deux phases principales que l’on peut modéliser par les équations différentielles de charge et de décharge d’un condensateur.

La première phase est la dépolarisation. Lorsqu’un stimulus d’intensité suffisante atteint le neurone, il provoque l’ouverture brutale de canaux sodiques dépendants du voltage. Les ions $Na^+$ s’engouffrent massivement dans la cellule. Ce flux de charges positives inverse la polarité de la membrane : la tension $V_m$ remonte rapidement de $-70\text{ mV}$ vers $+40\text{ mV}$.

On modélise cette charge par un circuit série RC alimenté par une force électromotrice $E$ (représentant la tension théorique d’équilibre des ions sodium). En appliquant la loi des mailles :

$$u_R(t) + u_C(t) = E$$

D’après la loi d’Ohm ($u_R = R \cdot i$) et la relation liant l’intensité au flux de charge ($i = C \cdot \frac{du_C}{dt}$), on obtient l’équation différentielle linéaire du premier ordre régissant la tension aux bornes de la membrane :

$$R \cdot C \cdot \frac{du_C(t)}{dt} + u_C(t) = E$$

Posons la constante de temps du neurone $\tau = R \cdot C$. L’équation se récrit :

$$\tau \cdot \frac{du_C(t)}{dt} + u_C(t) = E$$

La solution analytique de cette équation, correspondant à la dépolarisation à partir d’une tension initiale nulle pour simplifier, s’écrit :

$$u_C(t) = E \cdot \left(1 - e^{-t/\tau}\right)$$

La deuxième phase est la repolarisation (et l’hyperpolarisation). Juste après le pic de dépolarisation, les canaux sodiques se ferment et des canaux potassiques dépendants du voltage s’ouvrent, permettant aux ions $K^+$ de sortir de la cellule pour rétablir la polarité initiale.

On modélise ce retour à l’équilibre par la décharge du condensateur dans la résistance, en l’absence de générateur externe ($E = 0$) :

$$\tau \cdot \frac{du_C(t)}{dt} + u_C(t) = 0$$

La solution de cette équation différentielle homogène s’écrit :

$$u_C(t) = U_{\text{max}} \cdot e^{-t/\tau}$$

Où $U_{\text{max}}$ est la tension maximale atteinte au pic du potentiel d’action. La tension redescend de manière exponentielle vers sa valeur de repos.

---

Conclusion

En conclusion, la biophysique montre que le vivant obéit aux mêmes lois fondamentales que nos circuits électroniques. La membrane lipidique d’un neurone agit comme un condensateur plan, tandis que ses canaux ioniques font office de résistances. La résolution des équations différentielles de charge et de décharge du dipôle RC permet de décrire mathématiquement les variations temporelles de la tension lors du potentiel d’action.

Cependant, le modèle RC classique est linéaire car la résistance $R$ y est constante. Dans un neurone réel, la résistance change de façon dynamique et non-linéaire à chaque milliseconde avec l’ouverture et la fermeture des canaux. C’est cette complexité qui a été modélisée en 1952 par Alan Hodgkin et Andrew Huxley, leur valant le prix Nobel pour avoir posé les bases de la neurobiologie quantitative moderne.

Pour mon projet d’orientation, je souhaite intégrer une licence de biophysique ou un double cursus physique-biologie. Je suis convaincue que c’est à l’intersection de ces disciplines que se trouvent les réponses aux grands mystères du cerveau et de la modélisation des réseaux neuronaux artificiels.

Je vous remercie pour votre attention et je me réjouis de répondre à vos questions.

---

Questions Potentielles du Jury

1. Démontrez par analyse dimensionnelle que la constante de temps $\tau = R \cdot C$ est bien homogène à un temps.

Réponse : Pour réaliser l’analyse dimensionnelle, recherchons les dimensions de la résistance $R$ et de la capacité $C$ à partir des lois physiques fondamentales.

1. Dimension de la résistance $R$ : D’après la loi d’Ohm, $U = R \cdot I \implies [R] = [U] \cdot [I]^{-1}$.
2. Dimension de la capacité $C$ : La relation liant la charge d’un condensateur à sa tension est $q = C \cdot U$. De plus, l’intensité du courant est le débit de charge : $I = \frac{dq}{dt} \implies [q] = [I] \cdot [T]$. Par substitution : $$C = \frac{q}{U} \implies [C] = [q] \cdot [U]^{-1} = [I] \cdot [T] \cdot [U]^{-1}$$
3. Dimension du produit $R \cdot C$ : Multiplions les dimensions obtenues : $$[R \cdot C] = [R] \cdot [C] = \left([U] \cdot [I]^{-1}\right) \cdot \left([I] \cdot [T] \cdot [U]^{-1}\right)$$ En simplifiant les termes $[U]$ (tension) et $[I]$ (intensité) : $$[R \cdot C] = [T]$$ Le produit $R \cdot C$ a bien la dimension d’un temps. Son unité dans le système international est la seconde ($\text{s}$).

2. Qu’est-ce que la myélinisation des axones et comment influence-t-elle la vitesse de propagation de l’influx nerveux en modifiant les paramètres électriques de la membrane ?

Réponse : La gaine de myéline est une structure lipidique isolante qui entoure certains axones, interrompue à intervalles réguliers par les nœuds de Ranvier. Sur le plan électrique, la gaine de myéline agit de deux façons :

1. Elle augmente l’épaisseur $e$ de l’isolant : La capacité d’un condensateur plan est inversement proportionnelle à l’épaisseur de l’isolant ($C = \epsilon \cdot \frac{S}{e}$). En augmentant $e$, la myéline divise la capacité membranaire $C$ par un facteur supérieur à $100$.
2. Elle réduit la constante de temps $\tau = RC$ : Puisque la capacité $C$ chute drastiquement, la constante de temps $\tau$ diminue fortement aux zones myélinisées. La membrane se charge et se décharge de façon quasi-instantanée. L’influx nerveux ne se propage pas de façon continue mais « saute » d’un nœud de Ranvier à un autre (où la concentration de canaux ioniques est maximale) : c’est la conduction saltatoire. Cela permet d’augmenter la vitesse de propagation de l’influx de $1\text{ m/s}$ (axone non myélinisé) à plus de $100\text{ m/s}$ (axone myélinisé), tout en réduisant l’énergie consommée par les pompes ioniques.

3. Montrez par le calcul que la fonction $u_C(t) = E \cdot (1 - e^{-t/\tau})$ est l’unique solution de l’équation différentielle $\tau \cdot \frac{du_C(t)}{dt} + u_C(t) = E$ vérifiant $u_C(0) = 0$.

Réponse :

1. Calcul de la dérivée : $$\frac{du_C(t)}{dt} = \frac{d}{dt}\left(E - E \cdot e^{-t/\tau}\right) = -E \cdot \left(-\frac{1}{\tau} \cdot e^{-t/\tau}\right) = \frac{E}{\tau} \cdot e^{-t/\tau}$$
2. Substitution dans le membre de gauche de l’équation différentielle : $$\tau \cdot \frac{du_C(t)}{dt} + u_C(t) = \tau \cdot \left(\frac{E}{\tau} \cdot e^{-t/\tau}\right) + E \cdot \left(1 - e^{-t/\tau}\right)$$ $$\tau \cdot \frac{du_C(t)}{dt} + u_C(t) = E \cdot e^{-t/\tau} + E - E \cdot e^{-t/\tau} = E$$ L’égalité est vérifiée, la fonction est donc bien solution.
3. Vérification de la condition initiale : $$u_C(0) = E \cdot \left(1 - e^{-0/\tau}\right) = E \cdot (1 - 1) = 0$$ La condition initiale $u_C(0) = 0$ est respectée. (L’unicité découle des théorèmes mathématiques de résolution des équations différentielles linéaires du 1er ordre).

4. Qu’est-ce que le modèle de Hodgkin-Huxley et pourquoi est-il supérieur au modèle RC simple pour décrire un neurone réel ?

Réponse : Le modèle de Hodgkin-Huxley (1952) remplace le circuit RC simple par un circuit plus réaliste où :

1. La capacité $C$ représente toujours la membrane lipidique.
2. Les résistances associées aux canaux de sodium ($R_{Na}$) et de potassium ($R_K$) ne sont plus constantes, mais sont modélisées comme des résistances variables contrôlées par le temps et par la tension elle-même ($V_m$).
3. Le modèle intègre des générateurs de tension distincts pour chaque type d’ion ($E_{Na}$, $E_K$, $E_L$), représentant leurs potentiels d’équilibre électrochimique respectifs (calculés par l’équation de Nernst). Ce modèle est supérieur car il s’agit d’un système de 4 équations différentielles non-linéaires couplées qui reproduit le phénomène de seuil d’excitation (la loi du tout ou rien) et la phase d’hyperpolarisation réfractaire, ce qu’un simple circuit RC linéaire à composants fixes est incapable de simuler.

5. Comment les ions sodium et potassium traversent-ils la membrane cellulaire s’ils ne peuvent pas passer à travers la bicouche phospholipidique hydrophobe ?

Réponse : Les ions sodium ($Na^+$) et potassium ($K^+$) sont des particules chargées hydrophiles entourées d’une cage d’hydratation (molécules d’eau). Ils ne peuvent pas traverser spontanément la partie centrale hydrophobe (composée de chaînes d’acides gras) de la bicouche phospholipidique. Pour traverser, ils doivent passer par des protéines transmembranaires spécifiques :

1. Les canaux ioniques passifs : Ce sont des pores hydrophiles sélectifs qui laissent passer les ions par diffusion simple dans le sens de leur gradient électrochimique.
2. Les canaux ioniques dépendants du voltage : Ce sont des protéines qui changent de conformation spatiale (s’ouvrent ou se ferment) sous l’effet de variations de la tension membranaire, jouant le rôle d’interrupteurs électroniques.
3. Les pompes actives (comme la pompe $Na^+/K^+$) : Ce sont des enzymes qui transportent les ions contre leur gradient électrochimique en hydrolysant de l’ATP pour fournir l’énergie nécessaire.