Dynamics underlying epileptic seizures span multiple scales in space and time, therefore, understanding seizure mechanisms requires identifying the relations between seizure components within and across these scales, together with the analysis of their dynamical repertoire. In this view, mathematical models have been developed, ranging from single neuron to neural population. In this study, we consider a neural mass model able to exactly reproduce the dynamics of heterogeneous spiking neural networks. We combine mathematical modeling with structural information from non invasive brain imaging, thus building large-scale brain network models to explore emergent dynamics and test the clinical hypothesis. We provide a comprehensive study on the effect of external drives on neuronal networks exhibiting multistability, in order to investigate the role played by the neuroanatomical connectivity matrices in shaping the emergent dynamics. In particular, we systematically investigate the conditions under which the network displays a transition from a low activity regime to a high activity state, which we identify with a seizure-like event. This approach allows us to study the biophysical parameters and variables leading to multiple recruitment events at the network level. We further exploit topological network measures in order to explain the differences and the analogies among the subjects and their brain regions, in showing recruitment events at different parameter values. We demonstrate, along with the example of diffusion-weighted magnetic resonance imaging (dMRI) connectomes of 20 healthy subjects and 15 epileptic patients, that individual variations in structural connectivity, when linked with mathematical dynamic models, have the capacity to explain changes in spatiotemporal organization of brain dynamics, as observed in network-based brain disorders. In particular, for epileptic patients, by means of the integration of the clinical hypotheses on the epileptogenic zone (EZ), i.e., the local network where highly synchronous seizures originate, we have identified the sequence of recruitment events and discussed their links with the topological properties of the specific connectomes. The predictions made on the basis of the implemented set of exact mean-field equations turn out to be in line with the clinical pre-surgical evaluation on recruited secondary networks.

The structural connectivity of human brain allows the coexistence of segregated and integrated states of activity. Neuromodulatory systems facilitate the transition between these functional states and recent computational studies have shown how an interplay between the noradrenergic and cholinergic systems define these transitions. However, there is still much to be known about the interaction between the structural connectivity and the effect of neuromodulation, and to what extent the connectome facilitates dynamic transitions. In this work, we use a whole brain model, based on the Jasen and Rit equations plus a human structural connectivity matrix, to find out which structural features of the human connectome network define the optimal neuromodulatory effects. We simulated the effect of the noradrenergic system as changes in filter gain, and studied its effects related to the global-, local-, and meso-scale features of the connectome. At the global-scale, we found that the ability of the network of transiting through a variety of dynamical states is disrupted by randomization of the connection weights. By simulating neuromodulation of partial subsets of nodes, we found that transitions between integrated and segregated states are more easily achieved when targeting nodes with greater connection strengths—local feature—or belonging to the rich club—meso-scale feature. Overall, our findings clarify how the network spatial features, at different levels, interact with neuromodulation to facilitate the switching between segregated and integrated brain states and to sustain a richer brain dynamics.

A great deal of research has been devoted on the investigation of neural dynamics in various network topologies. However, only a few studies have focused on the influence of autapses, synapses from a neuron onto itself via closed loops, on neural synchronization. Here, we build a random network with adaptive exponential integrate-and-fire neurons coupled with chemical synapses, equipped with autapses, to study the effect of the latter on synchronous behavior. We consider time delay in the conductance of the pre-synaptic neuron for excitatory and inhibitory connections. Interestingly, in neural networks consisting of both excitatory and inhibitory neurons, we uncover that synchronous behavior depends on their synapse type. Our results provide evidence on the synchronous and desynchronous activities that emerge in random neural networks with chemical, inhibitory and excitatory synapses where neurons are equipped with autapses.