Prediction and control of focal seizure spread: Random walk with restart on heterogeneous brain networks

Modeling the interplay between regional heterogeneity and critical dynamics underlying brain functional networks

The human brain exhibits heterogeneity across regions and network connectivity patterns; However, how these heterogeneities contribute to whole-brain network functions and cognitive capacities remains unclear. In this study, we focus on the regional heterogeneity reflected in local dynamics and study how it contributes to the emergence of functional connectivity patterns, network ignition dynamics of the empirical brains. We find that the level of synchrony among voxelwise neural activities measured from the fMRI data is significantly correlated with the transcriptional variations in excitatory and inhibitory receptor gene expression. Consequently, we construct heterogeneous whole-brain network models with nodal excitability calibrated by the synchronization measure of regional dynamics. We demonstrate that as the extent of heterogeneity increases, the models operating around the critical point between order and disorder generate simulated functional connectivity networks increasingly similar to empirical resting-state or working memory task-evoked function connectivity networks. Furthermore, the heterogeneous models can predict individual differences in resting-state and task-evoked reconfiguration of the functional connectivity, as well as the comparative causal effect of empirical brain networks-that is, how the dynamics of one brain region affect whole-brain synchronization. Overall, this work demonstrates the viability of using regional heterogeneous functional signals to improve the performance of the whole-brain models, and illustrates how regional heterogeneity in human brains interplays with structural connections and critical dynamics to contribute to the emergence of functional connectivity networks.

Exploration of interictal to ictal transition in epileptic seizures using a neural mass model

An epileptic seizure can usually be divided into three stages: interictal, preictal, and ictal. However, the seizure underlying the transition from interictal to ictal activities in the brain involves complex interactions between inhibition and excitation in groups of neurons. To explore this mechanism at the level of a single population, this paper employed a neural mass model, named the complete physiology-based model (cPBM), to reconstruct electroencephalographic (EEG) signals and to infer the changes in excitatory/inhibitory connections related to excitation-inhibition (E-I) balance based on an open dataset recorded for ten epileptic patients. Since epileptic signals display spectral characteristics, spectral dynamic causal modelling (DCM) was applied to quantify these frequency characteristics by maximizing the free energy in the framework of power spectral density (PSD) and estimating the cPBM parameters. In addition, to address the local maximum problem that DCM may suffer from, a hybrid deterministic DCM (H-DCM) approach was proposed, with a deterministic annealing-based scheme applied in two directions. The H-DCM approach adjusts the temperature introduced in the objective function by gradually decreasing the temperature to obtain relatively good initialization and then gradually increasing the temperature to search for a better estimation after each maximization. The results showed that (i) reconstructed EEG signals belonging to the three stages together with their PSDs can be reproduced from the estimated parameters of the cPBM; (ii) compared to DCM, traditional D-DCM and anti D-DCM, the proposed H-DCM shows higher free energies and lower root mean square error (RMSE), and it provides the best performance for all stages (e.g., the RMSEs between the reconstructed PSD computed from the reconstructed EEG signal and the sample PSD obtained from the real EEG signal are 0.33 ± 0.08, 0.67 ± 0.37 and 0.78 ± 0.57 in the interictal, preictal and ictal stages, respectively); and (iii) the transition from interictal to ictal activity can be explained by an increase in the connections between pyramidal cells and excitatory interneurons and between pyramidal cells and fast inhibitory interneurons, as well as a decrease in the self-loop connection of the fast inhibitory interneurons in the cPBM. Moreover, the E-I balance, defined as the ratio between the excitatory connection from pyramidal cells to fast inhibitory interneurons and the inhibitory connection with the self-loop of fast inhibitory interneurons, is also significantly increased during the epileptic seizure transition.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11571-023-09976-6.

The distribution and heterogeneity of excitability in focal epileptic network potentially contribute to the seizure propagation

Introduction

Existing dynamical models can explain the transmigration mechanisms involved in seizures but are limited to a single modality. Combining models with networks can reproduce scaled epileptic dynamics. And the structure and coupling interactions of the network, as well as the heterogeneity of both the node and network activities, may influence the final state of the network model.

Methods

We built a fully connected network with focal nodes prominently interacting and established a timescale separated epileptic network model. The factors affecting epileptic network seizure were explored by varying the connectivity patterns of focal network nodes and modulating the distribution of network excitability.

Results

The whole brain network topology as the brain activity foundation affects the consistent delayed clustering seizure propagation. In addition, the network size and distribution heterogeneity of the focal excitatory nodes can influence seizure frequency. With the increasing of the network size and averaged excitability level of focal network, the seizure period decreases. In contrast, the larger heterogeneity of excitability for focal network nodes can lower the functional activity level (average degree) of focal network. There are also subtle effects of focal network topologies (connection patterns of excitatory nodes) that cannot be ignored along with non-focal nodes.

Discussion

Unraveling the role of excitatory factors in seizure onset and propagation can be used to understand the dynamic mechanisms and neuromodulation of epilepsy, with profound implications for the treatment of epilepsy and even for the understanding of the brain.