AbstractUsing The Virtual Brain (TVB, thevirtualbrian.org) simulation platform, we explored for 50 individual adult human brains (ages 18-80), how personalized connectome based brain network modelling captures various empirical observations as measured by functional magnetic resonance imaging (fMRI) and electroencephalography (EEG). We compare simulated activity based on individual structural connectomes (SC) inferred from diffusion weighted imaging with fMRI and EEG in the resting state. We systematically explore the role of the following model parameters: conduction velocity, global coupling and graph theoretical features of individual SC. First, a subspace of the parameter space is identified for each subject that results in realistic brain activity, i.e. reproducing the following prominent features of empirical EEG-fMRI activity: topology of resting-state fMRI functional connectivity (FC), functional connectivity dynamics (FCD), electrophysiological oscillations in the delta (3-4 Hz) and alpha (8-12 Hz) frequency range and their bimodality, i.e. low and high energy modes. Interestingly, FCD fit, bimodality and static FC fit are highly correlated. They all show their optimum in the same range of global coupling. In other words, only when our local model is in a bistable regime we are able to generate switching of modes in our global network. Second, our simulations reveal the explicit network mechanisms that lead to electrophysiological oscillations, their bimodal behaviour and inter-regional differences. Third, we discuss biological interpretability of the Stefanescu-Jirsa-Hindmarsh-Rose-3D model when embedded inside the large-scale brain network and mechanisms underlying the emergence of bimodality of the neural signal.With the present study, we set the cornerstone for a systematic catalogue of spatiotemporal brain activity regimes generated with the connectome-based brain simulation platform The Virtual Brain.Author SummaryIn order to understand brain dynamics we use numerical simulations of brain network models. Combining the structural backbone of the brain, that is the white matter fibres connecting distinct regions in the grey matter, with dynamical systems describing the activity of neural populations we are able to simulate brain function on a large scale. In order to make accurate prediction with this network, it is crucial to determine optimal model parameters. We here use an explorative approach to adjust model parameters to individual brain activity, showing that subjects have their own optimal point in the parameter space, depending on their brain structure and function. At the same time, we investigate the relation between bistable phenomena on the scale of neural populations and the changed in functional connectivity on the brain network scale. Our results are important for future modelling approaches trying to make accurate predictions of brain function.
Probabilistic models for neural populations that naturally capture global coupling and criticality