A coupled ordinary differential equation lattice model for the simulation of epileptic seizures

A coupled ordinary differential equation lattice model for the CA3 region of the hippocampus (a common location of the epileptic focus) is developed. This model consists of a hexagonal lattice of nodes, each describing a subnetwork consisting of a group of prototypical excitatory pyramidal cells and a group of prototypical inhibitory interneurons connected via on/off excitatory and inhibitory synapses. The nodes communicate using simple rules to simulate the diffusion of extracellular potassium. Both the integration time over which a node's trajectory is integrated before the diffusional event is allowed to occur and the level of inhibition in each node were found to be important parameters. Shorter integration times lead to total synchronization of the lattice (similar to synchronous neural activity occurring during a seizure) whereas longer times cause more random spatiotemporal behavior. Moderately diminished levels of inhibition lead to simple nodal oscillatory behavior. It is postulated that both the lack of inhibition and an alteration in conduction time may be necessary for the development of a behaviorally manifest seizure. (c) 1999 American Institute of Physics.