Computational modeling of high frequency oscillations recorded with clinical intracranial macroelectrodes

High Frequency Oscillations (HFOs) are a potential biomarker of epileptogenic regions. They have been extensively investigated in terms of automatic detection, classification and feature extraction. However, the mechanisms governing the generation of HFOs as well as the observability conditions on clinical intracranial macroelectrodes remain elusive. In this paper, we propose a novel physiologically-relevant macroscopic model for accurate simulation of HFOs as invasively recorded in epileptic patients. This model accounts for both the temporal and spatial properties of the cortical patch at the origin of epileptiform activity. Indeed, neuronal populations are combined with a 3D geometrical representation to simulate an extended epileptic source. Then, by solving the forward problem, the contributions of neuronal population signals are projected onto intracerebral electrode contacts. The obtained signals are qualitatively and quantitatively compared to real HFOs, and a relationship is drawn between macroscopic model parameters such as synchronization and spatial extent on the one hand, and HFO features such as the wave and fast ripple (200-600 Hz) components, on the other hand.