Whole brain functional connectivity: Insights from next generation neural mass modelling incorporating electrical synapses.

The ready availability of brain connectome data has both inspired and facilitated the modelling of whole brain activity using networks of phenomenological neural mass models that can incorporate both interaction strength and tract length between brain regions. Recently, a new class of neural mass model has been developed from an exact mean field reduction of a network of spiking cortical cell models with a biophysically realistic model of the chemical synapse. Moreover, this new population dynamics model can naturally incorporate electrical synapses. Here we demonstrate the ability of this new modelling framework, when combined with data from the Human Connectome Project, to generate patterns of functional connectivity (FC) of the type observed in both magnetoencephalography and functional magnetic resonance neuroimaging. Some limited explanatory power is obtained via an eigenmode description of frequency-specific FC patterns, obtained via a linear stability analysis of the network steady state in the neigbourhood of a Hopf bifurcation. However, direct numerical simulations show that empirical data is more faithfully recapitulated in the nonlinear regime, and exposes a key role of gap junction coupling strength in generating empirically-observed neural activity, and associated FC patterns and their evolution. Thereby, we emphasise the importance of maintaining known links with biological reality when developing multi-scale models of brain dynamics. As a tool for the study of dynamic whole brain models of the type presented here we further provide a suite of C++ codes for the efficient, and user friendly, simulation of neural mass networks with multiple delayed interactions.