Effects of physiological parameter evolution on the dynamics of tonic-clonic seizures

• Published in: PloS one Vol. 15; no. 4; p. e0230510
• Main Authors: Deeba, F, Sanz-Leon, P, Robinson, P A
• Format: Journal Article
• Language: English
• Published: United States Public Library of Science 02.04.2020; Public Library of Science (PLoS)
• Subjects: Analysis; Biology and Life Sciences; Brain; Cerebral cortex; Computer and Information Sciences; Cortex (temporal); Hopf bifurcation; Medicine and Health Sciences; Oscillations; Physical Sciences; Physics; Physiological aspects; Physiological effects; Physiology; Power (Philosophy); Research and Analysis Methods; Seizures; Seizures (Medicine); Spectra; Spectral analysis; Spectrum analysis; Strength; Thalamus; Time; Australia
• ISSN: 1932-6203; 1932-6203
• DOI: 10.1371/journal.pone.0230510

The temporal and spectral characteristics of tonic-clonic seizures are investigated using a neural field model of the corticothalamic system in the presence of a temporally varying connection strength between the cerebral cortex and thalamus. Increasing connection strength drives the system into ∼ 10 Hz seizure oscillations once a threshold is passed and a subcritical Hopf bifurcation occurs. In this study, the spectral and temporal characteristics of tonic-clonic seizures are explored as functions of the relevant properties of physiological connection strengths, such as maximum strength, time above threshold, and the ramp rate at which the strength increases or decreases. Analysis shows that the seizure onset time decreases with the maximum connection strength and time above threshold, but increases with the ramp rate. Seizure duration and offset time increase with maximum connection strength, time above threshold, and rate of change. Spectral analysis reveals that the power of nonlinear harmonics and the duration of the oscillations increase as the maximum connection strength and the time above threshold increase. A secondary limit cycle at ∼ 18 Hz, termed a saddle-cycle, is also seen during seizure onset and becomes more prominent and robust with increasing ramp rate. If the time above the threshold is too small, the system does not reach the 10 Hz limit cycle, and only exhibits 18 Hz saddle-cycle oscillations. It is also seen that the time to reach the saturated large amplitude limit-cycle seizure oscillation from both the instability threshold and from the end of the saddle-cycle oscillations is inversely proportional to the square root of the ramp rate.