Le modèle de Lotka-Volterra

1. Le modèle de Lotka-Volterra

2. Version Volterra • Croissance exponentielle proies x en absence de prédateurs • Décroissance exponentielle des prédateurs y en absence de proies • L’ « interaction » proies-prédateurs est proportionnelle au nombre de rencontres entre les deux espèces, donc à leurs densités respectives(les coeff. B et D sont différents car la disparition d’une proie n’entraîne pas de facto la naissance d’un prédateur)

3. Modèle de Lotka-Volterra k1 Z + X 2X k2 Y + X 2Y Version Lotka k3 Y W k2 = B k3 = C Si [Z] >> [X] alors k1Z=cste=A

5. Etats stationnaires 1ère solution : Premier état stationnaire : u0=0, v0=0 2ème solution : Deuxième état stationnaire : u0=1, v0=1

6. Analyse des modes normaux Matrice jacobienne Eléments de la jacobienne Jacobienne à l’état stationnaire 1 (u0=0, v0=0) Jacobienne à l’état stationnaire 2 (u0=1, v0=1)

7. Jacobienne à l’état stationnaire 1 (u0=0, v0=0) Analyse des modes normaux Trace de la jacobienne : Déterminant de la jacobienne : Equation caractéristique : Discriminant de l’équation caractéristique : La constante a étant toujours positive, D est négatif, ce qui signifie que les deux racines R1 et R2 sont de signes contraires. L'état stationnaire est donc instable et de type point de selle. Les trajectoires s'écartent de cet état. 5 6

8. Jacobienne à l’état stationnaire 2 (u0=1, v0=1) Analyse des modes normaux Trace de la jacobienne : Déterminant de la jacobienne : Equation caractéristique : Discriminant de l’équation caractéristique : Le discriminant étant toujours négatif, l'équation caractéristique admet deux racines toujours imaginaires. La constante a étant toujours positive, D est positif. La trace est toujours nulle. L'état stationnaire est un centre (cas très particulier). Les trajectoires sont (dans le système linéarisé) des cercles centrés sur l'état stationnaire. 3 2

10. Sensibilité du modèle de Lotka-Volterra aux conditions initiales

12. Oscillations glycolytiques

13. Pye, 1971 Extrait S. carlbergensis Milieu de culture : tréhalose

15. Influence de la vitesse d’injection du substrat

17. + S P vPFK (a) (g) a

18. + S P vin vPFK vout

20. Etat(s) stationnaire(s) du système soit avec dont l’unique racine positive est :

21. Analyse des modes normaux

22. Analyse des modes normaux : déterminant de la jacobienne L’état stationnaire (a*, g*) ne peut donc être en aucun cas du type point de selle.

23. Analyse des modes normaux : trace de la jacobienne Condition d’instabilité : Tr* > 0 soit

24. Analyse des modes normaux : discriminant La valeur stationnaire du discriminant peut être positive ou négative. L’état stationnaire peut donc être, a priori, stable ou instable, de type « nœud » ou « foyer », selon la valeur des paramètres.

25. 1. 0 < µ < µ1(d) = 0.032144 tr* < 0, d* > 0 nœud stable 2. µ1(d) = 0.032144 < µ < µ1(tr) = 0.03598 tr* < 0, d* < 0 foyer stable 3. µ1(tr) = 0.03598 < µ < µ2(d) = 0.04098 tr* > 0, d* < 0 foyer instable 4. µ2(d) = 0.04098 < µ < µ3(d) = 0.3574 tr* > 0, d* > 0 nœud instable 5. µ3(d) = 0.3574 < µ < µ2(tr) = 0.4572 tr* > 0, d* < 0 foyer instable 6. µ2(tr) = 0.4572 < µ < µ4(d) = 0.5543 tr* < 0, d* < 0 foyer stable 7. µ4(d) = 0.5543 < µ tr* < 0, d* > 0 nœud stable µ = 0.36 3

26. 1. 0 < µ < µ1(d) = 0.032144 tr* < 0, d* > 0 nœud stable 2. µ1(d) = 0.032144 < µ < µ1(tr) = 0.03598 tr* < 0, d* < 0 foyer stable 3. µ1(tr) = 0.03598 < µ < µ2(d) = 0.04098 tr* > 0, d* < 0 foyer instable 4. µ2(d) = 0.04098 < µ < µ3(d) = 0.3574 tr* > 0, d* > 0 nœud instable 5. µ3(d) = 0.3574 < µ < µ2(tr) = 0.4572 tr* > 0, d* < 0 foyer instable 6. µ2(tr) = 0.4572 < µ < µ4(d) = 0.5543 tr* < 0, d* < 0 foyer stable 7. µ4(d) = 0.5543 < µ tr* < 0, d* > 0 nœud stable µ = 0.45 3

27. 1. 0 < µ < µ1(d) = 0.032144 tr* < 0, d* > 0 nœud stable 2. µ1(d) = 0.032144 < µ < µ1(tr) = 0.03598 tr* < 0, d* < 0 foyer stable 3. µ1(tr) = 0.03598 < µ < µ2(d) = 0.04098 tr* > 0, d* < 0 foyer instable 4. µ2(d) = 0.04098 < µ < µ3(d) = 0.3574 tr* > 0, d* > 0 nœud instable 5. µ3(d) = 0.3574 < µ < µ2(tr) = 0.4572 tr* > 0, d* < 0 foyer instable 6. µ2(tr) = 0.4572 < µ < µ4(d) = 0.5543 tr* < 0, d* < 0 foyer stable 7. µ4(d) = 0.5543 < µ tr* < 0, d* > 0 nœud stable µ = 0.55 2

28. Méthode des isoclines nulles

35. Diffusion Première loi de Fick : le flux matériel J (cellules, concentration, nombre d’individus appartenant à une population, etc.) est proportionnel au gradient de « concentration » de ce matériel, ce qui, dans un espace à une dimension, s’exprime mathématiquement par : D ( coefficient de diffusion) a pour dimension [longueur]2[temps]-1 (cm2 sec-1 par exemple). Deuxième loi de Fick : Equation générale d’un processus de réaction-diffusion : (diffusion à une dimension)

36. Diffusion Réaction – diffusion à une dimension : Formalisme plus général : avec, dans le cas d’une diffusion monodirectionnelle : Cas général d’une diffusion à 3 dimensions :

40. Scilab cd ..\g

41. Qu‘est-ce que le phénomène d’excitabilité ? • Dans quelles conditions paramétriques peut-il être observé ? • Pourquoi n’observe-t-on jamais qu’une seule vague (et non, par exemple, un train d’oscillations) ? • Pourquoi lorsqu’ils se rencontrent, les deux fronts de vague s’annihilent-ils (ce qui a pour conséquence de faire que chaque point du cortex est soumis une et une seule fois au passage de la vague) ? • Quel est le mécanisme expliquant la transmission d’une unité à l’autre de la fluctuation initiale (phénomène de relais) ? La transmission est-elle soumise à des contraintes paramétriques particulières ?

42. L’excitabilité est une propriété d’un état stationnaire stable d’un système qui posséde au moins une isocline nulle ayant une branche de pente négative (forme en S). Dans la plupart des zones de pente positive, l’état stationnaire est très stable : le système ne peut amplifier un pulse de concentration. Cependant, à mesure que l’état stationnaire s’approche de l’extremum de la branche stable passant par l’origine (0, 0), le système devient excitable, c’est-à-dire capable d’amplifier un pulse de concentration (ici g). La zone d’excitabilité est en général très étroite.

43. Condition à l’excitation effective d’un système excitable : Pour que la fluctuation soit amplifiée, il faut que son amplitude dépasse un seuil. Le seuil diminue à mesure que l’état stationnaire soumis à la fluctuation approche le domaine instable. L’excitabilité se définit ainsi comme la capacité d’un système à amplifier, au-delà d’un seuil, un pulse d’une de ses variables, à partir d’un état stationnaire stable.

44. Trajectoire dans le plan de phase, suite à une excitation supraliminaire

45. A t= 0, Dg=0.4

46. Perturbations et période réfractaire Dg = 2.0 Dg = 0.4 Dg = 0.2 g seuil

47. Période réfractaire annihilation de la vague

48. Le phénomène de relais • Transmission de l’excitation d’une unité à l’autre, via la diffusion du calcium cytosolique Programme VagueCaRelais (nb_cycles = 100) Existence d’une limite supérieure et d’une limite inférieure à la diffusion entre lesquelles le signal d’excitation amplifié est transmis entre unités corticales. Détermination des valeurs limites de d pour lesquelles le relais se produit

49. De l'excitabilité aux oscillations µ = 0.2 µ = 0.0333

50. Relais d'oscillations Tous les compartiments basculent dans un état oscillant : synchronisation des oscillations, en phase dans tous les compartiments Un seul compartiment bascule dans un état oscillant : ondes concentriques en provenance d’un émetteur Une rangée de compartiments basculent dans un état oscillant : vagues successives Dans les 2 derniers cas, comment sont relayées les oscillations par les unités restées dans l’état excitable ? Programme VagueCaRelais