Phase-dependent stimulation effects on bursting activity in a neural network cortical simulation

Purpose: A neural network simulation with realistic cortical architecture has been used to study synchronized bursting as a seizure representation. This model has the property that bursting epochs arise and cease spontaneously, and bursting epochs can be induced by external stimulation. We have used this simulation to study the time-frequency properties of the evolving bursting activity, as well as effects due to network stimulation. Methods: The model represents a cortical region of 1.6 mm × 1.6 mm, and includes seven neuron classes organized by cortical layer, inhibitory or excitatory properties, and electrophysiological characteristics. There are a total of 65,536 modeled single compartment neurons that operate according to a version of Hodgkin-Huxley dynamics. The intercellular wiring is based on histological studies and our previous modeling efforts. Results: The bursting phase is characterized by a flat frequency spectrum. Stimulation pulses are applied to this modeled network, with an electric field provided by a 1 mm radius circular electrode represented mathematically in the simulation. A phase dependence to the post-stimulation quiescence is demonstrated, with local relative maxima in efficacy occurring before or during the network depolarization phase in the underlying activity. Brief periods of network insensitivity to stimulation are also demonstrated. The phase dependence was irregular and did not reach statistical significance when averaged over the full 2.5 s of simulated bursting investigated. This result provides comparison with previous in vivo studies which have also demonstrated increased efficacy of stimulation when pulses are applied at the peak of the local field potential during cortical after discharges. The network bursting is synchronous when comparing the different neuron classes represented up to an uncertainty of 10 ms. Studies performed with an excitatory chandelier cell component demonstrated increased synchronous bursting in the model, as predicted from experimental work. Conclusions: This large-scale multi-neuron neural network simulation reproduces many aspects of evolving cortical bursting behavior as well as the timing-dependent effects of electrical stimulation on that bursting.