Traumatic brain injury (TBI) may result in post-traumatic seizures and epilepsy. Approximately 5-7% ofTBI patients suffer from at least one seizure [1]. Thepathophysiological mechanisms are not completelyunderstood, and may also differ between early seizures(< 2 weeks of the injury) and the development of post-traumatic epilepsy. This includes the effects of directphysical trauma, excitotoxicity due to iron released fromthe blood [2] and cytokine TGF-b in blood-brain-barrier-mediated activation of astrocytes [3]. In thiswork we propose that a reduction in network connectiv-ity, as presumed present in several cases of TBI, mayresult in seizures and epilepsy.We use a realistic model of neocortex consisting of sixdifferent, multi-compartmental neurons and Hodgkin-Huxley like ion-channel dynamics [4]. Using physiologi-cally realistic connectivityparameters, we analyze net-works of different sizes; ranging from a microcolumn of656 neurons to a mesocolumn that contains 20k neu-rons. In these networks, small lesions are introduced tosimulate axonal and dendritic damage, thereby limitingaction potential propagation. Furthermore, we analyze alumped model of neocortex that is shown to correspondto the detailed model of the microcolumn [5]. Thismodel consists of a system of two differential equationswith two fixed delays. By using an automated parameterestimation method, parameters are identified for whichthe model’s behavior closely resembles that of the realis-tic model. Subsequently, the dependency and sensitivityon these parameters are studied with bifurcation analy-sis. We generate a mesocolumn by linking severallumped units together. Lesions are then introduced bybreaking or reducing some of the connections betweenthe populations. We also study this case using bothanalytical and numerical bifurcation methods.The ratio between excitatory and inhibitory connec-tions is analytically determined as a function of networksize. It is found that, compared to large networks, smallnetworks tend to have a relatively larger number ofexcitatory connections than inhibitory connections. Thissuggests that a lesion splitting the network into smallersub-networks, could increase the ratio of excitatory andinhibitory?connections in a particular sub-network.Choosing parameters that correspond to a region ofmultistability, as determined by the bifurcation analysis,enables us to create an epileptic focus that spreadsepileptiform activity to neighboring areas.ConclusionsBy using multi-scale modeling, large-scale simulations,bifurcation analysis, and parameter estimation, we studythe effects of small lesions in neocortex. From the large-scale simulations we find that“neuronal peninsulas”,created by these lesions, may evolve into epileptogenicnetworks. By studying bifurcations of a lumped modelwith suitable parameters, regions of multistability areidentified that are hypothesized to correspond withepilepsy.
Read full abstract