Matière et Révolution

Periodic Stimulation of Biological Oscillators

Robert Paris — 2019-11-09T09:24:17Z

Synchronization and Rhythmic Processes in Physiology

Multiple Oscilators Provide Metastblity in Rhythm

Phase-locking

Zones of phase locking

Respiratory System

Cardiac Phase-locked Rhythms

Cardiac Phase Resetting

Phase-locking in Parasystole Rhythms

Van der Pol Oscillator

Phase-locking of the van der Pol Oscillator

Dynamical Features od Neurons and the Brain

Periodic Stimulation of Biological Oscillators

Periodic stimulation of spontaneously oscillating physiological rhythms has powerful effects on the intrinsic rhythm. As the frequency and amplitude of the periodic stimulus are varied, a variety of different coupling patterns are set up between the stimulus and the spontaneous oscillator. In some situations the spontaneous rhythm is entrained or phase locked to the forcing stimulus so that for each N cycles of the stimulus there are M cycles of the spontaneous rhythm, and the spontaneous oscillation occurs at fixed phase (or phases) of the periodic stimulus (N:M phase locking). In addition to phase-locked rhythms, it is also possible to observe irregular or aperiodic rhythms in which fixed phase relationships and regular repeating cyclic patterns are not observed. We will discuss the main experimental observations in phase locking experiments, with particular reference to experiments on the cardiac and respiratory systems. We will develop the mathematical concepts needed to analyse phase locking in integrate and fire models and in limit cycle models. We will discuss several phenomena entailing phase locking in humans.

One experimental paradigm for studying physiological oscillators is to subject the oscillator to periodic stimuli while maintaining physiological conditions as constant as possible. We do not attempr to give a complete summary of this large body of work but discuss two illustrative systems: the mechanical ventilation of animals, and the periodic stimulation of cardiac oscillations using an intracellular microelectrode. We draw generalizations from these systems which are broadly applicable in a wide variety of experimental systems.

It has been known since the time of Hering and Breuer that periodic lung inflation by a mechanical ventilator can lead to phase locking between the ventilator and the intrinsic respiratory rhythm in mammals. The entrainment is believed to be mediated by the Hering-Breuer reflexes in which expansion of the lungs inhibits inspiration ans prolongs expiration. Afferent activity from stretch receptors in the lungs is carried by afferent fibers in the vagus nerve. Respiratory entrainment can be studied by mechanically ventilating an animal at different volume and frequencies while attempting to maintain constant levels of anesthesia, body temperature, and blood gases.

Experiments were performed on paralyzed, pentobarbital-anesthetized adult cats, and central respiratory activity was monitored by recording from a branch of the phrenic nerve. In a normal unparalyzed animal, the phrenic nerve nerve innervates the diaphragm, and phrenic nerve activity causes the diaphragm to contract, thus leading to inspiration. In the paralyzed animal, neuromuscular transmission between the phrenic nerve and the diaphragm is blocked, and lung inflation is due solely to the mechanical ventilator. However, afferent activity from stretch receptors in the lung is still carried by the vagus nerve…

As the ventilator volume and frequency are varied, a number of different rhythms are established between the ventilator and phrenic activity. These different rhythms are organized in an orderly fashion in the ventilator volume-ventilator frequency plane. Insets show representative traces of ventilator volumes and phrenic nerve activity corresponding to different stable phase-locked rhythms, as well as non-phase-locked rhythms (which occur in the shaded regions). In these experiments, low ventilation must be maintained. Likewise, very high volumes and frequencies could not be studied because of the mechanical limitations of the ventilator and the limited lung capacity of the cat.

The respiratory rhythm is generated in a complex network of neurons in the brain stem. In contrast, the cardiac rhythm is generated in a specialized region of electrically coupled cells in the SA node that act as a pacemaker for the heart. We now describe experiments in which the effects of pulsatile electrical stimuli delivered to spontaneously beating cells derived from the ventricles of embryonic chick heart were determined.

As the frequency and current intensity of the electrical stimuli are varied, a variety of different rhythms between the stimulator and the heart cells are established… The insets represent the different observed phase-locking patterns, and the solid lines represent the results of theoretical computations based on phase-resetting experiments using single pulses. In addition to the stable phase-locked rhythms, there are a number of irregular rhythms…

Although thses results have been obtained from two very different physiological systems using different types of periodic stimulation, there are certain striking features common to both. The following generalizations are applicable to a large number of experiments of periodic forcing of biological oscillators.

1. The stable zones of phase locking most commonly observed correspond to low-order ratios between the number of cycles of the forcing stimulus and the intrinsic rhythm (i.e., 2:1 , 3:2 , 1:1 , 2:3 , 1:2). Although other N:M ratios with larger values of N and M can also be observed, these occupy smaller areas in the frequency- amplitude parameter space, and they are consequently easily overlooked or obscured by noise.

2. The stable rhythms are organized in the frequency-amplitude plane in an orderly fashion. It is common to associate a rotation number p = M/N with an N:M rhythm. Then, as the stimulation frequency increases at fixed stimulus amplitude, p decreases.

3. At very low stimulation amplitudes, it is difficult to maintain stable phase locking.

4. If the regions of frequency-amplitude parameter space between stable phase-locking zones are studied, then it is generally possible to find stimulation parameters that give rise to irregular dynamics.

Despite their similarities, there are differences between the cardiac and respiratory systems. For example, the 2:2 region observed in the periodically stimulated cardiac cells was not observed in the mechanically ventilated cats. The goal of theoretical studies of phase locking is to provide a theoretical basis for understanding the similarities and differences between the different preparations. Ideally, one would like to be able to make predictions about the phase locking as stimulation parameters vary, based on the mechanisms of rhythmogenesis and the coupling of the stimulator to the intrinsic rhythm. In practice, the mechanical analysis of periodically forced nonlinear oscillators is an extremely difficult problem, and detailed quantitative understanding of dynamics has been obtained in only a few special situations.

The theoretical studies make it clear that although there are many similarities between the dynamics in different systems, there will also generally be differences if the dynamics are studied in sufficient detail. In order to establish this assertion, we describe the dynamics resulting from periodic stimulation in a number of different model systems. We also discuss entrainment of biological oscillators in a number of clinical contexts. (…)

One motivation for the analysis of periodic forcing of integrate and fire models comes from the experimental studies on the mechanical ventilation of cats. An examination of the insets, in particular the 1:2 and 3:2 rhythms, shows large lung inflations coincident with the premature termination of the inspiratory activity. This gives a striking visual suggestion for the utility of an integrate and fire model with modulated threshold for this system. As we already discussed, the Hering-Breuer reflexes can be modeled by integrate and fire models. (…)

In addition to studies of integrate and fire models in which the thresholds are sinusoidal, there have been several studies of the properties of integrate and fire models in which the thresholds are piecewise linear. Such studies allow a much more complete analysis of the dynamics than is possible using sinusoidal thresholds. For example, Lasota and Mackey proved that such models display chaotic dynamics in some parameter ranges. The significance of this finding lies in the observation of chaos even in extremely simplified models for the periods forcing of oscillations. In view of this finding, it seems likely that the appearance of chaotic dynamics at least over some range of stimulation parameters in periodically forced physiological oscillators may be very common.

Attemps to develop more realistic integrate and fire models for the entrainment of the respiratory rhythm have been made by Petrillo and Glass. They assume that two thresholds are modulated by the volume of the mechanical ventilator. The timing of inspiration and expiration is represented by activities that oscillate between the two thresholds. A brief delay occurs between the time the inspiratory onset threshold was reached and the beginning of inspiration. This model has five parameters, all of which could be determined from the data. Numerical simulations of the model show good agreement with the experimental results. However, a detailed study of the bifurcations in this model was not carried through.

Another attempt to develop integrate and fire models for the entrainment of biological rhythms was made by Daan and coworkers for the circadian rhythm. They assumed that there were two sinusoidally modulated thresholds. With this model they were able to provide a partial explanation for data on sleep duration as a function of the time of onset entrainment (corresponding to the normal circadian rhythm), or of other possible phase-locking zones as a function of the parameters.

In conclusion, the assumption that there are periodic inputs to integrate and fire models provides a conceptually simple means to model phase locking in biological systems. Even the simplest possible models lead to extremely complicated dynamics, which are only partially understood.

Many biological rhythms are best represented mathematically as limit cycle oscillations in differential equations. Because these rhythms interact with each other and because there is a periodic stimulation from the external environment, it is important to understand the effects of periodic forcing on a limit cycle oscillation. One prototypical model for periodically forced limit cycles is the sinusoidally forced van der Pol equation. J.H. Jensen and coworkers have shown that sinusoidal forcing also gives rise to chaotic dynamics in mathematical models of excitable neural and cardiac tissue, and Aihara and coworkers have demonstrated strange attractors from sinusoidally forced squid axons.

We now consider the effects of a periodic train of short pulsatile stimuli on limit cycle oscillations. In the event that the limit cycle is rapidly restablished following a single stimulus, it is straightforward to compute the effect of periodic stimuli once the effects of a single stimulus are understood.

The main idea can be developed from a consideration of the effects of a periodic stimulation of the Poincaré oscillator...”

Zones of phase locking from Glass and Petrillo (1984)

Theoretically computed phase-locking zones from Glass, Guevara and Shrier (1987)

Regions of phase locking of the van der Pol oscillator from Hayashi (1964)

Theoretical model of phase-locking from Petrillo and Glass (1984)

Phase locking of respiration from Bramble (1983)

Phase-locking zones in a mathematical model of parasystole from Moe et al. (1977)

Phase locked rhythms set up between a mechanical ventilator and the respiratory rhythms in an anesthetized human being from Graves et al. (1986)

Glass-Mackey, « The Rythms of Life »

Rythmes non-linéaires auto-organisés

Robert Paris — 2008-09-30T19:29:53Z

Réactions biochimiques, créateurs de formes et de rythmes, le film

L'oscillateur Van der Pol

Phase-locking ou interaction entre rythmes auto-organisés

Simon Diner

Extrait de « Les voies du chaos dans l'école russe », tiré de l'ouvrage collectif « Chaos et déterminisme » :

NON LINEARITE DES RYTHMES AUTO-ORGANISES (DISSIPATIFS ET ENTRETENUS)

« Dans les oscillations non-linéaires, l'ordre et le désordre se côtoient, se relaient, se confortent, voilà la surprise. (…) C'est l'instauration d'une véritable conception dialectique de l'ordre et du désordre qui n'a pas fini de nous étonner. »
Ce que prolonge Dahan Dalmedico dans « Retour sur l'histoire de la philosophie » du même ouvrage :
« L'étude des systèmes dynamiques chaotiques exige une véritable dialectique entre l'instabilité d'un système dynamique chaotique et sa stabilité structurelle. »

(…)

« L'école de Mandelstham-Andronov et le paradigme des auto-oscillations

Parmi tous les mouvements mécaniques et physiques, les oscillations occupent une place à part Ce sont les mouvements ou les changements d'état qui présentent un certain degré de répétitivité ou de périodicité. Dans le cas le plus simple, celui des petits mouvements d'un pendule balançoire, ou d'un ressort, la force responsable du mouvement reste simplement proportionnelle aux déplacements du système. Ce sont les oscillations linéaires dont l'oscillateur harmonique est le modèle universel. A vrai dire, toutes les oscillations qui existent réellement dans la nature sont plus ou moins non linéaires. Les oscillations linéaires ne sont qu'un modèle mathématique approché dont l'importance est liée au rôle mathématique joué par les fonctions périodiques (analyse de Fourier). Au 19ème siècle, le modèle de l'oscillateur linéaire s'est imposé à travers le développement de l'étude des ondes en optique, en acoustique et en électromagnétisme. Ce caractère universel du modèle d'oscillateur linéaire triomphe dans le célèbre livre de Lord Rayleigh : « The theory of sound » (1877).
L'oscillateur linéaire va aussi sous-tendre et structurer toute la physique quantique. De par l'utilisation qu'elle fait de la théorie des opérateurs linéaires dans l'espace vectoriel des états (espace de Hilbert), la mécanique quantique est comme l'apothéose d'un paradigme en développement depuis deux siècles. (…) En fait, le modèle d'oscillateur linéaire convient parfaitement aux phénomènes stationnaires, ceux dont l'évolution présente l'harmonie et la régularité d'un équilibre mobile. Mais dès qu'il s'agit d'évolutions temporelles dramatiques, en particulier de phénomènes de transition d'un état d'équilibre à un autre, de phénomène de création ou de disparition de mouvements, la non-linéarité devient une propriété organique essentielle. C'est précisément le cas lorsqu'on engendre et entretient des oscillations à partir de phénomènes non oscillatoires : chute d'un poids dans une horloge, frottement de l'archet dans un violon, souffle de l'instrumentaliste dans une flûte, émission de la voix humaine.
L'horlogerie est le domaine privilégié des oscillations non linéaires. (…) Lord Rayleigh est le premier à distinguer les traits caractéristiques des systèmes susceptibles d'engendrer des oscillations non amorties, en particulier la non-linéarité des équations du mouvement.
En fait, l'intérêt pour les études théoriques des oscillations non linéaires ne va pas venir de l'horlogerie mais de deux autres domaines techniques. Dans la seconde moitié du 19ème siècle, la construction de régulateurs pour les machines devient un problème technologique essentiel. L'enjeu est d'empêcher l'apparitions d'oscillations. Worms et Romilly (1872) et L.A. Vichnegradski (1876) reconnaissent la nécessité du frottement pour la stabilisation des régulateurs. H. Léauté (1885) montre le rôle essentiel joué par la non-linéarité dans certains types de régulateurs que l'on ne peut étudier par linéarisation comme le faisaient I.A. Vichnegradski et A. Stodola (1893).
Au début du 20ème siècle, commence la réalisation de dispositifs radiotechniques, pour l'émission et la réception des ondes électromagnétiques. Il s'agit là d'engendrer des oscillations. Dans les dispositifs radiotechniques, tenir compte de la non-linéarité s'avère essentiel, mais les approches restèrent longtemps ad hoc. C'est dans cet esprit que s'effectuèrent en particulier les travaux fondamentaux du radiophysicien hollandais B. Van der Pol. (…)
Au moment même où naît la mécanique quantique, A.A Andronov (élève de l'école russe de L.I. Mandelstham, qui a choisi comme domaine de recherche l'étude des vibrations non linéaires contrairement à la physique quantique) participe à l'émergence d'un nouveau paradigme dont l'acte fondateur, son travail de diplôme, paraît en français dans les « Comptes rendus de l'Académie des sciences » du 14 octobre 1929 : « Les cycles-limites de Poincaré et la théorie des oscillations auto-entretenues ». Andronov y reconnaît pour la première fois que dans un oscillateur de la radiophysique comme celui de Van der Pol, système non-conservatif (dissipatif), dont les oscillations sont entretenues en puisant de l'énergie à des sources non vibratoires, le mouvement dans l'espace des phases est du type « cycle limite », notion introduite par Poincaré en 1880, dans un contexte purement mathématique. Il reconnaît d'emblée la caractère très général de ces « auto-oscillations » comme il les nomme, les voyant intervenir en acoustique, en radiophysique, en chimie (réactions périodiques) et en biologie.
Les auto-oscillations ont des caractéristiques spécifiques :
– amplitude et fréquence indépendantes des conditions initiales
– apparition en l'absence d'excitation périodique extérieure
– contrôle par rétroaction, de la source d'énergie pour compenser la dissipation, sans influer sur l'amplitude et la fréquence (…)

L'auto-oscilateur reçoit une définition très générale : système engendrant des oscillations non amorties, entretenue par une source d'énergie extérieure, dans un dispositif non-linéaire dissipatif, et dont l'aspect et les propriétés sont déterminés par le système lui-même sans dépendre des conditions initiales. Dans ces conditions, les auto-oscillations peuvent être non seulement périodiques mais quasi périodiques et même stochastiques. Andronov a eu en effet le mérite de montrer pour la prmeière fois l'existence physique d'un attracteur qui ne soit pas un point d'équilibre. Le cycle limite est en effet un attracteur périodique. Par la suite, la notion d'attracteur sera élargie, jusqu'à l'apparition du concept d' »attracteur étrange », forme mathématique des auto-oscillations stochastiques.
On attribue souvent à E. Lorentz la découverte en 1963 du premier « mouvement chaotique sur un attracteur étrange ». Découverte qui passa inaperçue et ne commença à être reconnue que dans la seconde moitié des années 70. On ne sait pas que dans les années 50, les travaux de l'école de Gorki, sous la direction d'un élève d'Andronov, Yu. I. Neimark, ont mis en évidence l'existence d'auto-oscillations schochastiques , par l'application de la méthode des transformations ponctuelles. (…) Avec les systèmes auto-oscillants comme cible favorite, Andronov et ses élèves vont baliser tout le champ des vibrations non linéaires et créer les outils et les concepts fondamentaux de la physique non linéaire. Dès 1933, dans son rapport à la première conférence soviétique sur les vibrations, Andronov développe le thème de la théorie des bifurcations, c'est-à-dire du changement de caractère qualitatif du portrait de phase d'un système dynamique, lors de la variation des paramètres du système. (…) En 1937 paraît la bible des vibrations non-linéaires : « La théorie des vibrations » de A. A. Andronov, A. A. Vitt et S. E. Khaïkir. La signature de A. A. Vitt en a disparu dans la seconde édition, car il a été assassiné lors des grandes purges staliniennes.
(…) Dès le début des années quarante, Kolmogorov s'intéresse en probabiliste à la turbulence. (…) Dans les années cinquante, il passe à l'étude des systèmes dynamiques. (…) Dans son exposé au congrès international de mathématiques d'Amsterdam de 1954, il présente une splendide synthèse des résultats obtenus depuis H. Poincaré. (…) Et Kolmogorov formule la première version du résultat fondamental qui va devenir, quelques années plus tard, le théorème de Kolmogorov, Arnold et Moser (théorème KAM) sur la préservation des mouvements quasi périodiques dans les systèmes hamiltoniens. »