Measuring resetting of brain dynamics at epileptic seizures: application of global optimization and spatial synchronization techniques

Epileptic seizures are manifestations of intermittent spatiotemporal transitions of the human brain from chaos to order. Measures of chaos, namely maximum Lyapunov exponents (STL _max ), from dynamical analysis of the electroencephalograms (EEGs) at critical sites of the epileptic brain, progressively converge (diverge) before (after) epileptic seizures, a phenomenon that has been called dynamical synchronization (desynchronization). This dynamical synchronization/desynchronization has already constituted the basis for the design and development of systems for long-term (tens of minutes), on-line, prospective prediction of epileptic seizures. Also, the criterion for the changes in the time constants of the observed synchronization/desynchronization at seizure points has been used to show resetting of the epileptic brain in patients with temporal lobe epilepsy (TLE), a phenomenon that implicates a possible homeostatic role for the seizures themselves to restore normal brain activity. In this paper, we introduce a new criterion to measure this resetting that utilizes changes in the level of observed synchronization/desynchronization. We compare this criterion’s sensitivity of resetting with the old one based on the time constants of the observed synchronization/desynchronization. Next, we test the robustness of the resetting phenomena in terms of the utilized measures of EEG dynamics by a comparative study involving STL _max , a measure of phase (φ _max ) and a measure of energy (E) using both criteria (i.e. the level and time constants of the observed synchronization/desynchronization). The measures are estimated from intracranial electroencephalographic (iEEG) recordings with subdural and depth electrodes from two patients with focal temporal lobe epilepsy and a total of 43 seizures. Techniques from optimization theory, in particular quadratic bivalent programming, are applied to optimize the performance of the three measures in detecting preictal entrainment. It is shown that using either of the two resetting criteria, and for all three dynamical measures, dynamical resetting at seizures occurs with a significantly higher probability (α=0.05) than resetting at randomly selected non-seizure points in days of EEG recordings per patient. It is also shown that dynamical resetting at seizures using time constants of STL _max synchronization/desynchronization occurs with a higher probability than using the other synchronization measures, whereas dynamical resetting at seizures using the level of synchronization/desynchronization criterion is detected with similar probability using any of the three measures of synchronization. These findings show the robustness of seizure resetting with respect to measures of EEG dynamics and criteria of resetting utilized, and the critical role it might play in further elucidation of ictogenesis, as well as in the development of novel treatments for epilepsy.