Spirales cardiaques : quand la réentrée électrique entretient l'arythmie

1. Introduction

Chaque battement du cœur résulte d'une onde électrique qui se propage de façon ordonnée à travers le myocarde, déclenchant la contraction musculaire synchronisée qui propulse le sang. Cette onde naît dans le nœud sinusal, traverse les oreillettes, franchit le nœud auriculo-ventriculaire, puis envahit les ventricules en quelques dizaines de millisecondes. Lorsque cette propagation ordonnée est perturbée, une onde peut se refermer sur elle-même et tourner indéfiniment : c'est la réentrée. La réentrée est la cause principale des tachycardies ventriculaires et de la fibrillation, arythmies potentiellement mortelles qui constituent la première cause de mort subite dans les pays industrialisés.

La géométrie de ces ondes de réentrée est celle d'une spirale. Depuis les travaux pionniers de Wiener et Rosenblueth (1946) ^[1] et les expériences de Winfree sur les milieux excitables chimiques (1972) ^[2], puis les simulations numériques de Panfilov et Holden (1990) ^[3], il est établi que le tissu cardiaque se comporte comme un milieu excitable bidimensionnel (ou tridimensionnel) dans lequel des ondes spirales peuvent se former, s'ancrer sur des hétérogénéités anatomiques, dériver, se fragmenter et engendrer une turbulence électrique. Comprendre la physique de ces spirales est devenu un enjeu majeur de la cardiologie computationnelle.

2. Le tissu cardiaque comme milieu excitable

Un milieu excitable est un système qui possède trois propriétés fondamentales : (i) un état de repos stable, (ii) une réponse de tout-ou-rien à une perturbation dépassant un seuil, et (iii) une période réfractaire pendant laquelle le milieu ne peut être ré-excité. Le myocarde réunit ces trois conditions. Chaque cardiomyocyte est couplé électriquement à ses voisins par des jonctions communicantes (connexines 40, 43, 45), formant un syncytium fonctionnel. La propagation de l'onde dépolarisante est assurée par les courants ioniques transmembranaires : courant sodique rapide I_Na (dépolarisation), courant calcique lent I_CaL (plateau), et courant potassique I_K (repolarisation).

La vitesse de conduction longitudinale dans le myocarde ventriculaire est d'environ 0,5–1 m/s, et la vitesse transversale est 3 à 5 fois plus faible en raison de l'anisotropie des fibres musculaires. Cette anisotropie joue un rôle crucial dans la formation et la stabilité des spirales. La longueur d'onde électrique λ = v × DRA (où v est la vitesse de conduction et DRA la durée de la période réfractaire absolue) détermine la taille minimale d'un circuit de réentrée. Pour le ventricule humain, λ ≈ 20–30 cm, ce qui signifie que la réentrée n'est possible que si le circuit est suffisamment grand ou si la longueur d'onde est réduite par une ischémie ou un médicament.

3. Le modèle FitzHugh-Nagumo

Le modèle de FitzHugh-Nagumo (FHN) est la réduction mathématique canonique d'un milieu excitable. Initialement proposé par FitzHugh (1961) ^[4] comme simplification du modèle de Hodgkin-Huxley, il décrit la dynamique d'une variable rapide u (analogue du potentiel de membrane) et d'une variable lente v (variable de récupération, analogue à l'inactivation des canaux) :

D : coefficient de diffusion (couplage spatial) ; α : seuil d'excitation (0 < α < 1) ; ε ≪ 1 : séparation d'échelles de temps ; γ : paramètre de récupération ; I_ext : courant externe.

Dans ce modèle, la nullcline de u (∂u/∂t = 0 sans diffusion) est une cubique en forme de N, et la nullcline de v est une droite. Leur intersection unique définit le point fixe stable (état de repos). Lorsque u dépasse le seuil α, le système effectue une excursion rapide vers la branche haute de la cubique (dépolarisation), puis revient lentement (repolarisation) via la branche basse. La période réfractaire correspond au temps pendant lequel v est encore élevé et empêche une nouvelle excitation.

En dimension deux, avec le terme de diffusion D∇²u, le modèle FHN supporte des ondes planes, des ondes circulaires (cibles), et surtout des ondes spirales. La spirale tourne autour d'un noyau de rotation (le « rotor ») où u et v oscillent périodiquement. La fréquence de rotation de la spirale est supérieure à la fréquence naturelle du tissu, ce qui explique pourquoi la réentrée prend le contrôle du rythme cardiaque.

4. Ondes spirales et réentrée

Une onde spirale dans le tissu cardiaque est une onde de dépolarisation qui tourne autour d'un noyau de phase singulière. Ce noyau est un point (en 2D) ou un filament (en 3D) autour duquel toutes les phases du cycle d'excitation coexistent simultanément. Winfree (1987) ^[5] a montré que ce point singulier est topologiquement inévitable : si une onde est brisée (par une hétérogénéité, une ischémie localisée, ou une stimulation prématurée), ses extrémités libres s'enroulent nécessairement en spirale.

Le mécanisme classique de création d'une spirale est le protocole S1-S2 : une première stimulation S1 crée une onde plane ; une seconde stimulation S2, appliquée pendant la queue réfractaire de S1, crée une onde dont une extrémité est bloquée (tissu encore réfractaire) et l'autre se propage librement. Cette extrémité libre s'enroule en spirale. Ce mécanisme a été démontré expérimentalement par Davidenko et al. (1992) ^[6] sur des tranches de myocarde ventriculaire de chien, en utilisant la cartographie optique par colorants fluorescents sensibles au potentiel.

La forme de la spirale cardiaque est une spirale d'Archimède dans les modèles simples, mais dans les modèles ioniques réalistes elle est plus proche d'une spirale logarithmique dont le pas varie avec la distance au noyau. La vitesse de rotation Ω et le rayon du noyau r₀ sont déterminés par les propriétés locales du tissu (vitesse de conduction, durée du potentiel d'action, pente de la courbe de restitution).

λ : longueur d'onde électrique ; v : vitesse de conduction ; DRA : durée de la période réfractaire absolue ; L_circuit : longueur du circuit anatomique ou fonctionnel.

5. De la tachycardie à la fibrillation

Une spirale stable et solitaire produit une tachycardie ventriculaire monomorphe : le ventricule est activé à la fréquence de rotation de la spirale (typiquement 5–10 Hz, soit 300–600 bpm), bien supérieure à la fréquence sinusale normale (1–2 Hz). La tachycardie est hémodynamiquement mal tolérée car la fréquence élevée réduit le temps de remplissage ventriculaire et donc le débit cardiaque.

La transition vers la fibrillation ventriculaire correspond à la fragmentation de la spirale en de multiples spirales filles, créant une turbulence électrique. Ce phénomène est lié à l'instabilité de la spirale, elle-même gouvernée par la courbe de restitution du potentiel d'action. Karma (1993) ^[7] a montré analytiquement que la spirale devient instable lorsque la pente de la courbe de restitution dépasse 1 :

DPA : durée du potentiel d'action ; DI : intervalle diastolique (temps de récupération entre deux potentiels d'action successifs).

Lorsque la pente de restitution est supérieure à 1, une légère variation de l'intervalle diastolique est amplifiée à chaque cycle, produisant une alternance de potentiels d'action longs et courts (alternance électrique). Cette alternance spatiale crée des hétérogénéités de repolarisation qui fragmentent le front d'onde en de multiples spirales indépendantes. La fibrillation ventriculaire est ainsi un état de turbulence spatio-temporelle dans un milieu excitable.

6. Ancrage, dérive et turbulence

Dans un tissu homogène, le noyau d'une spirale peut dériver lentement (mouvement de translation du rotor) ou rester fixe. En présence d'hétérogénéités anatomiques (cicatrices d'infarctus, zones fibrotiques, orifices valvulaires), la spirale peut s'ancrer sur ces obstacles. L'ancrage stabilise la spirale et rend la tachycardie persistante. Panfilov et Vasiev (1995) ^[8] ont montré que la spirale s'ancre préférentiellement sur des obstacles dont la taille est comparable au rayon du noyau r₀.

La dérive d'une spirale libre dans un milieu légèrement inhomogène peut être décrite par les équations de mouvement du noyau. En présence d'un gradient de paramètre (par exemple un gradient de durée du potentiel d'action dû à une ischémie), la spirale dérive perpendiculairement au gradient (dérive de Doppler) et parallèlement au gradient (dérive de gradient). Ces deux composantes de dérive ont été calculées par Biktashev et Holden (1995) ^[9] en utilisant la théorie des perturbations et les fonctions de réponse de la spirale.

En trois dimensions, l'analogue de la spirale est un scroll wave (onde de défilement), dont le noyau est un filament. Ce filament peut être droit, courbé, ou noué. La dynamique du filament est régie par des équations de courbure-tension : un filament courbé se contracte (tension positive) ou s'étend (tension négative) selon les paramètres du milieu. La tension négative du filament conduit à une instabilité de Winfree-Strogatz qui fragmente le scroll wave en de nombreux filaments enchevêtrés, produisant la turbulence électrique tridimensionnelle observée dans la fibrillation ventriculaire.

7. Modèles ioniques réalistes

Au-delà du modèle FHN, des modèles ioniques réalistes ont été développés pour reproduire fidèlement le potentiel d'action cardiaque humain. Le modèle de Luo-Rudy (1991, 1994) ^{[10] [11]} décrit les courants ioniques ventriculaires de cobaye avec une précision quantitative. Le modèle de Ten Tusscher et Panfilov (2006) ^[12] est le modèle de référence pour le ventricule humain ; il comprend 19 variables d'état et 12 courants ioniques distincts.

Ces modèles sont couplés à des équations de diffusion pour simuler la propagation dans des géométries cardiaques réalistes reconstruites par IRM. La simulation d'un ventricule humain complet à résolution cellulaire (~10⁸ nœuds) requiert des ressources de calcul massives (supercalculateurs). Des approches multi-échelles (modèles de câble, éléments finis, méthodes de domaine fictif) permettent de réduire le coût computationnel tout en préservant la précision des prédictions cliniques.

C_m : capacité membranaire (~1 μF/cm²) ; V_m : potentiel de membrane ; I_ion : somme des courants ioniques transmembranaires ; w : vecteur des variables de porte des canaux ; ρ_i : résistivité intracellulaire ; σ_i : tenseur de conductivité intracellulaire (anisotrope).

8. Défibrillation et contrôle des spirales

La défibrillation électrique consiste à appliquer un choc électrique intense (200–360 J pour un défibrillateur externe) pour interrompre simultanément toutes les spirales en cours. Le mécanisme n'est pas simplement une dépolarisation globale du tissu : Efimov et al. (1998) ^[13] ont montré par cartographie optique que le choc crée des zones de dépolarisation et d'hyperpolarisation virtuelles (électrodes virtuelles) qui génèrent de nouvelles ondes, lesquelles peuvent éliminer les spirales existantes si leur timing est correct.

Des approches alternatives de contrôle des spirales ont été explorées. La stimulation à faible énergie (LEAP — Low-Energy Anti-fibrillation Pacing) utilise une séquence de chocs de faible amplitude pour progressivement synchroniser et éliminer les spirales ^[14]. Le contrôle par rétroaction (feedback control) utilise des mesures en temps réel de l'activité électrique pour appliquer des stimulations ciblées qui déstabilisent le rotor. Ces approches visent à réduire l'énergie du choc de défibrillation de 80 à 90 %, réduisant ainsi les dommages tissulaires et la douleur.

La cartographie des rotors par électrodes endocavitaires (système FIRM — Focal Impulse and Rotor Modulation) permet d'identifier les sites d'ancrage des spirales dans la fibrillation auriculaire et de les ablater par radiofréquence ^[15]. Cette approche illustre comment la compréhension physique des spirales cardiaques se traduit directement en stratégies thérapeutiques.

9. Conclusion

Les ondes spirales cardiaques sont l'une des manifestations les plus dramatiques et les plus étudiées de la physique des milieux excitables. Leur existence résulte d'une combinaison de propriétés universelles — excitabilité, réfractarité, couplage diffusif — et de particularités anatomiques du myocarde. La convergence entre la physique théorique des spirales (modèles FHN, théorie des perturbations, turbulence spatio-temporelle) et la cardiologie clinique (cartographie optique, ablation par radiofréquence, défibrillation à faible énergie) est un exemple remarquable de la façon dont les mathématiques des spirales éclairent des phénomènes vitaux.

La spirale cardiaque n'est pas une curiosité mathématique : c'est une structure dynamique dont la formation, la stabilité et la destruction déterminent si un cœur bat ou s'emballe. Comprendre pourquoi et comment une onde électrique se réentre dans le tissu cardiaque, c'est comprendre l'une des frontières les plus actives entre la physique non linéaire et la médecine.