The Virtual Epileptic Patient: Individualized whole-brain models of epilepsy spread

Controllability of nonlinear epileptic-seizure spreading dynamics in large-scale subject-specific brain networks

Closed-loop electrical stimulation has become an important alternative to resective surgery for control of pharmacologically-resistant focal epileptic seizures. Seizure spread across large-scale brain networks, rather than its focal onset, is what commonly leads to major disruptions in sensorimotor and cognitive processing, as well as loss-of-consciousness, one of the main impairing aspects of the disorder. Electrical stimulation, triggered by early detection of seizure onset in epileptogenic zones (EZs), has been applied to prevent spread and its subsequent effects. Here, we show how linear feedback seizure-spread controllability in subject-specific (white-matter tractography) Epileptor network models is affected by variations in brain excitability, network coupling strength, control latency and gain, and actuation targets. Feedback control can qualitatively change the nonlinear seizure dynamics, and the paths to seizure termination and spread prevention. Notably, control onset latency is a critical parameter leading to a phase transition in spread controllability. Consequently, the efficacy of EZ-only actuation is limited depending on network excitability, coupling strength, and practical latencies for detection and actuation. Additional feedback-stabilization control of theoretically-derived optimal node subsets in the network are necessary for spread prevention. Finally, we contrast our linear-feedback controllability assessment with other measures based on minimum-energy (Gramian) controllability and nonlinear pulse-perturbation approaches.

Precision neurology

Over the past several decades, high-resolution brain imaging, blood and cerebrospinal fluid analyses, and other advanced technologies have changed diagnosis from an exercise depending primarily on the history and physical examination to a computer- and online resource-aided process that relies on larger and larger quantities of data. In addition, randomized controlled trials (RCT) at a population level have led to many new drugs and devices to treat neurological disease, including disease-modifying therapies. We are now at a crossroads. Combinatorially profound increases in data about individuals has led to an alternative to population-based RCTs. Genotyping and comprehensive "deep" phenotyping can sort individuals into smaller groups, enabling precise medical decisions at a personal level. In neurology, precision medicine that includes prediction, prevention and personalization requires that genomic and phenomic information further incorporate imaging and behavioral data. In this article, we review the genomic, phenomic, and computational aspects of precision medicine for neurology. After defining biological markers, we discuss some applications of these "-omic" and neuroimaging measures, and then outline the role of computation and ultimately brain simulation. We conclude the article with a discussion of the relation between precision medicine and value-based care.

Graph Neural Network Learning on the Pediatric Structural Connectome

PURPOSE

Sex classification is a major benchmark of previous work in learning on the structural connectome, a naturally occurring brain graph that has proven useful for studying cognitive function and impairment. While graph neural networks (GNNs), specifically graph convolutional networks (GCNs), have gained popularity lately for their effectiveness in learning on graph data, achieving strong performance in adult sex classification tasks, their application to pediatric populations remains unexplored. We seek to characterize the capacity for GNN models to learn connectomic patterns on pediatric data through an exploration of training techniques and architectural design choices.

METHODS

Two datasets comprising an adult BRIGHT dataset (N = 147 Hodgkin's lymphoma survivors and N = 162 age similar controls) and a pediatric Human Connectome Project in Development (HCP-D) dataset (N = 135 healthy subjects) were utilized. Two GNN models (GCN simple and GCN residual), a deep neural network (multi-layer perceptron), and two standard machine learning models (random forest and support vector machine) were trained. Architecture exploration experiments were conducted to evaluate the impact of network depth, pooling techniques, and skip connections on the ability of GNN models to capture connectomic patterns. Models were assessed across a range of metrics including accuracy, AUC score, and adversarial robustness.

RESULTS

GNNs outperformed other models across both populations. Notably, adult GNN models achieved 85.1% accuracy in sex classification on unseen adult participants, consistent with prior studies. The extension of the adult models to the pediatric dataset and training on the smaller pediatric dataset were sub-optimal in their performance. Using adult data to augment pediatric models, the best GNN achieved comparable accuracy across unseen pediatric (83.0%) and adult (81.3%) participants. Adversarial sensitivity experiments showed that the simple GCN remained the most robust to perturbations, followed by the multi-layer perceptron and the residual GCN.

CONCLUSIONS

These findings underscore the potential of GNNs in advancing our understanding of sex-specific neurological development and disorders and highlight the importance of data augmentation in overcoming challenges associated with small pediatric datasets. Further, they highlight relevant tradeoffs in the design landscape of connectomic GNNs. For example, while the simpler GNN model tested exhibits marginally worse accuracy and AUC scores in comparison to the more complex residual GNN, it demonstrates a higher degree of adversarial robustness.

Precision data-driven modeling of cortical dynamics reveals person-specific mechanisms underpinning brain electrophysiology

Task-free brain activity affords unique insight into the functional structure of brain network dynamics and has been used to identify neural markers of individual differences. In this work, we present an algorithmic optimization framework that directly inverts and parameterizes brain-wide dynamical-systems models involving hundreds of interacting neural populations, from single-subject M/EEG time-series recordings. This technique provides a powerful neurocomputational tool for interrogating mechanisms underlying individual brain dynamics ("precision brain models") and making quantitative predictions. We extensively validate the models' performance in forecasting future brain activity and predicting individual variability in key M/EEG metrics. Last, we demonstrate the power of our technique in resolving individual differences in the generation of alpha and beta-frequency oscillations. We characterize subjects based upon model attractor topology and a dynamical-systems mechanism by which these topologies generate individual variation in the expression of alpha vs. beta rhythms. We trace these phenomena back to global variation in excitatory-inhibitory balance, highlighting the explanatory power of our framework to generate mechanistic insights.

Advances in physiological and clinical relevance of hiPSC-derived brain models for precision medicine pipelines

Precision, or personalized, medicine aims to stratify patients based on variable pathogenic signatures to optimize the effectiveness of disease prevention and treatment. This approach is favorable in the context of brain disorders, which are often heterogeneous in their pathophysiological features, patterns of disease progression and treatment response, resulting in limited therapeutic standard-of-care. Here we highlight the transformative role that human induced pluripotent stem cell (hiPSC)-derived neural models are poised to play in advancing precision medicine for brain disorders, particularly emerging innovations that improve the relevance of hiPSC models to human physiology. hiPSCs derived from accessible patient somatic cells can produce various neural cell types and tissues; current efforts to increase the complexity of these models, incorporating region-specific neural tissues and non-neural cell types of the brain microenvironment, are providing increasingly relevant insights into human-specific neurobiology. Continued advances in tissue engineering combined with innovations in genomics, high-throughput screening and imaging strengthen the physiological relevance of hiPSC models and thus their ability to uncover disease mechanisms, therapeutic vulnerabilities, and tissue and fluid-based biomarkers that will have real impact on neurological disease treatment. True physiological understanding, however, necessitates integration of hiPSC-neural models with patient biophysical data, including quantitative neuroimaging representations. We discuss recent innovations in cellular neuroscience that can provide these direct connections through generative AI modeling. Our focus is to highlight the great potential of synergy between these emerging innovations to pave the way for personalized medicine becoming a viable option for patients suffering from neuropathologies, particularly rare epileptic and neurodegenerative disorders.

A New Perspective in Epileptic Seizure Classification: Applying the Taxonomy of Seizure Dynamotypes to Noninvasive EEG and Examining Dynamical Changes across Sleep Stages

Epilepsy, a neurological disorder characterized by recurrent unprovoked seizures, significantly impacts patient quality of life. Current classification methods focus primarily on clinical observations and electroencephalography (EEG) analysis, often overlooking the underlying dynamics driving seizures. This study uses surface EEG data to identify seizure transitions using a dynamical systems-based framework-the taxonomy of seizure dynamotypes-previously examined only in invasive data. We applied principal component and independent component (IC) analysis to surface EEG recordings from 1,177 seizures in 158 patients with focal epilepsy, decomposing the signals into ICs. The ICs were visually labeled for clear seizure transitions and bifurcation morphologies (BifMs), which were then examined using Bayesian multilevel modeling in the context of clinical factors. Our analysis reveals that certain onset bifurcations (saddle node on invariant circle and supercritical Hopf) are more prevalent during wakefulness compared with their reduced rate during nonrapid eye movement (NREM) sleep, particularly NREM3. We discuss the possible implications of our results in the context of modeling approaches and suggest additional avenues to continue this exploration. Furthermore, we demonstrate the feasibility of automating this classification process using machine learning, achieving high performance in identifying seizure-related ICs and classifying interspike interval changes. Our findings suggest that the noise in surface EEG may obscure certain BifMs, and we suggest technical improvements that could enhance detection accuracy. Expanding the dataset and incorporating long-term biological rhythms, such as circadian and multiday cycles, may provide a more comprehensive understanding of seizure dynamics and improve clinical decision-making.

Homeodynamic feedback inhibition control in whole-brain simulations

Simulations of large-scale brain dynamics are often impacted by overexcitation resulting from heavy-tailed structural network distributions, leading to biologically implausible simulation results. We implement a homeodynamic plasticity mechanism, known from other modeling work, in the widely used Jansen-Rit neural mass model for The Virtual Brain (TVB) simulation framework. We aim at heterogeneously adjusting the inhibitory coupling weights to reach desired dynamic regimes in each brain region. We show that, by using this dynamic approach, we can control the target activity level to obtain biologically plausible brain simulations, including post-synaptic potentials and blood-oxygen-level-dependent functional magnetic resonance imaging (fMRI) activity. We demonstrate that the derived dynamic Feedback Inhibitory Control (dFIC) can be used to enable increased variability of model dynamics. We derive the conditions under which the simulated brain activity converges to a predefined target level analytically and via simulations. We highlight the benefits of dFIC in the context of fitting the TVB model to static and dynamic measures of fMRI empirical data, accounting for global synchronization across the whole brain. The proposed novel method helps computational neuroscientists, especially TVB users, to easily "tune" brain models to desired dynamical regimes depending on the specific requirements of each study. The presented method is a steppingstone towards increased biological realism in brain network models and a valuable tool to better understand their underlying behavior.

Vagus nerve stimulation for epilepsy: A narrative review of factors predictive of response

Vagus nerve stimulation (VNS) is an established therapy for drug-resistant epilepsy. However, there is a lack of reliable predictors of VNS response in clinical use. The identification of factors predictive of VNS response is important for patient selection and stratification as well as tailored stimulation programming. We conducted a narrative review of the existing literature on prognostic markers for VNS response using clinical, demographic, biochemical, and modality-specific information such as from electroencephalography (EEG), magnetoencephalography, and magnetic resonance imaging (MRI). No individual marker demonstrated sufficient predictive power for individual patients, although several have been suggested, with some promising initial findings. Combining markers from underresearched modalities such as T1-weighted MRI morphometrics and EEG may provide better strategies for treatment optimization.

Methods for Identifying Epilepsy Surgery Targets Using Invasive EEG: A Systematic Review

BACKGROUND

The pre-surgical evaluation for drug-resistant epilepsy achieves seizure freedom in only 50-60% of patients. Efforts to identify quantitative intracranial EEG (qEEG) biomarkers of epileptogenicity are needed. This review summarizes and evaluates the design of qEEG studies, discusses barriers to biomarker adoption, and proposes refinements of qEEG study protocols.

METHODS

We included exploratory and prediction prognostic studies from MEDLINE and Scopus published between 2017 and 2023 that investigated qEEG markers for identifying the epileptogenic network as the surgical target. Cohort parameters, ground truth references, and analytical approaches were extracted.

RESULTS

Out of 1789 search results, 128 studies were included. The study designs were highly heterogeneous. Half of the studies included a non-consecutive cohort, with sample sizes ranging from 2 to 166 patients (median of 16). The most common minimum follow-up was one year, and the seizure onset zone was the most common ground truth. Prediction studies were heterogeneous in their analytical approaches, and only 25 studies validated the marker through post-surgical outcome prediction. Outcome prediction performance decreased in larger cohorts. Conversely, longer follow-up periods correlated with higher prediction accuracy, and connectivity-based approaches yielded better predictions. The data and code were available in only 9% of studies.

CONCLUSIONS

To enhance the validation qEEG markers, we propose standardizing study designs to resemble clinical trials. This includes using a consecutive cohort with long-term follow-up, validating against surgical resection as ground truth, and evaluating markers through post-surgical outcome prediction. These considerations would improve the reliability and clinical adoption of qEEG markers.

Secondary thalamic dysfunction underlies abnormal large-scale neural dynamics in chronic stroke

Stroke causes pronounced and widespread slowing of neural activity. Despite decades of work exploring these abnormal neural dynamics and their associated functional impairments, their causes remain largely unclear. To close this gap in understanding, we applied a neurophysiological corticothalamic circuit model to simulate magnetoencephalography (MEG) power spectra recorded from chronic stroke patients. Comparing model-estimated physiological parameters to those of controls, patients demonstrated significantly lower intrathalamic inhibition in the lesioned hemisphere, despite the absence of direct damage to the thalamus itself. We hypothesized that this disinhibition could instead be related to secondary degeneration of the thalamus, for which growing evidence exists in the literature. Further analyses confirmed that spectral slowing correlated significantly with overall secondary degeneration of the ipsilesional thalamus, encompassing decreased thalamic volume, altered tissue microstructure, and decreased blood flow. Crucially, this relationship was mediated by model-estimated thalamic disinhibition, suggesting a causal link between secondary thalamic degeneration and abnormal brain dynamics via thalamic disinhibition. Finally, thalamic degeneration was correlated significantly with poorer cognitive and language outcomes, but not lesion volume, reinforcing that thalamus damage may account for additional individual variability in poststroke disability. Overall, our findings indicate that the frequently observed poststroke slowing reflects a disruption of corticothalamic circuit dynamics due to secondary thalamic dysfunction, and highlights the thalamus as an important target for understanding and potentially treating poststroke brain dysfunction.

Electroencephalography derived connectivity informing epilepsy surgical planning: Towards clinical applications and future perspectives

Epilepsy is one of the most diffused neurological disorders, affecting 50 million people worldwide. Around 30% of patients have drug-resistant epilepsy (DRE), defined as failure of at least two tolerated antiseizure medications (ASMs) to achieve sustained seizure freedom. Brain surgery is an effective therapeutic approach in this group, hinging on the accurate localization of the epileptic focus. The latter task is complex and requires multimodal investigation methods. Epilepsy is also a network disorder and represents one of the best application scenarios of methods leveraging brain functional organization at large scales. Connectivity analysis represents a promising tool for improving surgical assessment, enabling better identification of candidates who could benefit the most from epilepsy surgery. The scalp electroencephalography (EEG) is the most relevant tool to characterize epileptic activity. The EEG has benefited significantly from technological advancement across the last decades. Firstly, electrical source imaging (ESI) allows the reconstruction of electrical activity detected by EEG at the cortex level; secondly, functional connectivity (FC) allows the assessment of functional dependencies across brain areas. The EEG has therefore expanded potential applications in the localization and characterization of the epileptogenic network for surgical planning. As the translation of these methods in clinical practice is little discussed in the literature, we reviewed the investigations using EEG-derived FC. We showed that the FC-informed identification of the epileptic networks improves the localization precision in focal epilepsy. We discussed the heterogeneity in the results and methodology preventing prompt research-to-clinic translation. We finally provided practical suggestions for promoting the applicability of FC-based research in real clinical practice, looking for future research.

Virtual resection evaluation based on sEEG propagation network for drug-resistant epilepsy

Drug-resistant epilepsy with frequent seizures are considered to undergo surgery to become seizure-free, but seizure-free rates have not dramatically improved, partly due to imprecise intervention locations. To address this clinical need, we construct effective connectivity to reveal epilepsy brain dynamics. Based on the propagation path captured by the high order effective connectivity, calculate the control centrality evaluation scheme of the excised area. We used three datasets: simulation dataset, clinical dataset, and public dataset. The epileptogenic propagation network was quantified by calculating high-order effective connection to obtain accurate propagation path, based on this, combined with the outdegree index for virtual resection. By removing electrodes and recalculating control centrality, we quantify each electrode or region's control centrality to evaluate the virtual resection scheme. Three datasets obtained consistent results. We track the accurate propagation path and find the obvious inflection points occurring during the excision process. The minimum intervention targets were obtained by comparing different schemes without recurrence. The clinical data with multiple seizures found that after resection, the brain reaches a stable state and is less likely to continue spreading. By quantitative analysis of control centrality to evaluate the possible excision scheme, finally we obtain the best intervention area for epilepsy, which assist in developing surgical plans.

Brain sodium MRI-derived priors support the estimation of epileptogenic zones using personalized model-based methods in epilepsy

Patients presenting with drug-resistant epilepsy are eligible for surgery aiming to remove the regions involved in the production of seizure activities, the so-called epileptogenic zone network (EZN). Thus the accurate estimation of the EZN is crucial. Data-driven, personalized virtual brain models derived from patient-specific anatomical and functional data are used in Virtual Epileptic Patient (VEP) to estimate the EZN via optimization methods from Bayesian inference. The Bayesian inference approach used in previous VEP integrates priors, based on the features of stereotactic-electroencephalography (SEEG) seizures' recordings. Here, we propose new priors, based on quantitative 23Na-MRI. The 23Na-MRI data were acquired at 7T and provided several features characterizing the sodium signal decay. The hypothesis is that the sodium features are biomarkers of neuronal excitability related to the EZN and will add additional information to VEP estimation. In this paper, we first proposed the mapping from 23Na-MRI features to predict the EZN via a machine learning approach. Then, we exploited these predictions as priors in the VEP pipeline. The statistical results demonstrated that compared with the results from current VEP, the result from VEP based on 23Na-MRI prior has better balanced accuracy, and the similar weighted harmonic mean of the precision and recall.

Spontaneous Brain Activity Emerges from Pairwise Interactions in the Larval Zebrafish Brain

Brain activity is characterized by brainwide spatiotemporal patterns that emerge from synapse-mediated interactions between individual neurons. Calcium imaging provides access to in vivo recordings of whole-brain activity at single-neuron resolution and, therefore, allows the study of how large-scale brain dynamics emerge from local activity. In this study, we use a statistical mechanics approach-the pairwise maximum entropy model-to infer microscopic network features from collective patterns of activity in the larval zebrafish brain and relate these features to the emergence of observed whole-brain dynamics. Our findings indicate that the pairwise interactions between neural populations and their intrinsic activity states are sufficient to explain observed whole-brain dynamics. In fact, the pairwise relationships between neuronal populations estimated with the maximum entropy model strongly correspond to observed structural connectivity patterns. Model simulations also demonstrated how tuning pairwise neuronal interactions drives transitions between observed physiological regimes and pathologically hyperexcitable whole-brain regimes. Finally, we use virtual resection to identify the brain structures that are important for maintaining the brain in a physiological dynamic regime. Together, our results indicate that whole-brain activity emerges from a complex dynamical system that transitions between basins of attraction whose strength and topology depend on the connectivity between brain areas.

Fingerprints of brain disease: connectome identifiability in Alzheimer's disease

Functional connectivity patterns in the human brain, like the friction ridges of a fingerprint, can uniquely identify individuals. Does this "brain fingerprint" remain distinct even during Alzheimer's disease (AD)? Using fMRI data from healthy and pathologically ageing subjects, we find that individual functional connectivity profiles remain unique and highly heterogeneous during mild cognitive impairment and AD. However, the patterns that make individuals identifiable change with disease progression, revealing a reconfiguration of the brain fingerprint. Notably, connectivity shifts towards functional system connections in AD and lower-order cognitive functions in early disease stages. These findings emphasize the importance of focusing on individual variability rather than group differences in AD studies. Individual functional connectomes could be instrumental in creating personalized models of AD progression, predicting disease course, and optimizing treatments, paving the way for personalized medicine in AD management.

The digital twin in neuroscience: from theory to tailored therapy

Digital twins enable simulation, comprehensive analysis and predictions, as virtual representations of physical systems. They are also finding increasing interest and application in the healthcare sector, with a particular focus on digital twins of the brain. We discuss how digital twins in neuroscience enable the modeling of brain functions and pathology as they offer an in-silico approach to studying the brain and illustrating the complex relationships between brain network dynamics and related functions. To showcase the capabilities of digital twinning in neuroscience we demonstrate how the impact of brain tumors on the brain's physical structures and functioning can be modeled in relation to the philosophical concept of plasticity. Against this technically derived backdrop, which assumes that the brain's nonlinear behavior toward improvement and repair can be modeled and predicted based on MRI data, we further explore the philosophical insights of Catherine Malabou. Malabou emphasizes the brain's dual capacity for adaptive and destructive plasticity. We will discuss in how far Malabou's ideas provide a more holistic theoretical framework for understanding how digital twins can model the brain's response to injury and pathology, embracing Malabou's concept of both adaptive and destructive plasticity which provides a framework to address such yet incomputable aspects of neuroscience and the sometimes seemingly unfavorable dynamics of neuroplasticity helping to bridge the gap between theoretical research and clinical practice.

Epileptic network identification: insights from dynamic mode decomposition of sEEG data

Objective.For medically-refractory epilepsy patients, stereoelectroencephalography (sEEG) is a surgical method using intracranial electrode recordings to identify brain networks participating in early seizure organization and propagation (i.e. the epileptogenic zone, EZ). If identified, surgical EZ treatment via resection, ablation or neuromodulation can lead to seizure-freedom. To date, quantification of sEEG data, including its visualization and interpretation, remains a clinical and computational challenge. Given elusiveness of physical laws or governing equations modelling complex brain dynamics, data science offers unique insight into identifying unknown patterns within high-dimensional sEEG data. We apply here an unsupervised data-driven algorithm, dynamic mode decomposition (DMD), to sEEG recordings from five focal epilepsy patients (three with temporal lobe, and two with cingulate epilepsy), who underwent subsequent resective or ablative surgery and became seizure free.Approach.DMD obtains a linear approximation of nonlinear data dynamics, generating coherent structures ('modes') defining important signal features, used to extract frequencies, growth rates and spatial structures. DMD was adapted to produce dynamic modal maps (DMMs) across frequency sub-bands, capturing onset and evolution of epileptiform dynamics in sEEG data. Additionally, we developed a static estimate of EZ-localized electrode contacts, termed the higher-frequency mode-based norm index (MNI). DMM and MNI maps for representative patient seizures were validated against clinical sEEG results and seizure-free outcomes following surgery.Main results.DMD was most informative at higher frequencies, i.e. gamma (including high-gamma) and beta range, successfully identifying EZ contacts. Combined interpretation of DMM/MNI plots best identified spatiotemporal evolution of mode-specific network changes, with strong concordance to sEEG results and outcomes across all five patients. The method identified network attenuation in other contacts not implicated in the EZ.Significance.This is the first application of DMD to sEEG data analysis, supporting integration of neuroengineering, mathematical and machine learning methods into traditional workflows for sEEG review and epilepsy surgical decision-making.

Predicting radiofrequency thermocoagulation surgical outcomes in refractory focal epilepsy patients using functional coupled neural mass model

Objective

The success rate of achieving seizure freedom after radiofrequency thermocoagulation surgery for patients with refractory focal epilepsy is about 20-40%. This study aims to enhance the prediction of surgical outcomes based on preoperative decisions through network model simulation, providing a reference for clinicians to validate and optimize surgical plans.

Methods

Twelve patients with epilepsy who underwent radiofrequency thermocoagulation were retrospectively reviewed in this study. A coupled model based on model subsets of the neural mass model was constructed by calculating partial directed coherence as the coupling matrix from stereoelectroencephalography (SEEG) signals. Multi-channel time-varying model parameters of excitation and inhibitions were identified by fitting the real SEEG signals with the coupled model. Further incorporating these model parameters, the coupled model virtually removed contacts destroyed in radiofrequency thermocoagulation or selected randomly. Subsequently, the coupled model after virtual surgery was simulated.

Results

The identified excitatory and inhibitory parameters showed significant difference before and after seizure onset (p < 0.05), and the trends of parameter changes aligned with the seizure process. Additionally, excitatory parameters of epileptogenic contacts were higher than that of non-epileptogenic contacts, and opposite findings were noticed for inhibitory parameters. The simulated signals of postoperative models to predict surgical outcomes yielded an area under the curve (AUC) of 83.33% and an accuracy of 91.67%.

Conclusion

The multi-channel coupled model proposed in this study with physiological characteristics showed a desirable performance for preoperatively predicting patients' prognoses.

Incorporating slow NMDA-type receptors with nonlinear voltage-dependent magnesium block in a next generation neural mass model: derivation and dynamics

We derive a next generation neural mass model of a population of quadratic-integrate-and-fire neurons, with slow adaptation, and conductance-based AMPAR, GABAR and nonlinear NMDAR synapses. We show that the Lorentzian ansatz assumption can be satisfied by introducing a piece-wise polynomial approximation of the nonlinear voltage-dependent magnesium block of NMDAR current. We study the dynamics of the resulting system for two example cases of excitatory cortical neurons and inhibitory striatal neurons. Bifurcation diagrams are presented comparing the different dynamical regimes as compared to the case of linear NMDAR currents, along with sample comparison simulation time series demonstrating different possible oscillatory solutions. The omission of the nonlinearity of NMDAR currents results in a shift in the range (and possible disappearance) of the constant high firing rate regime, along with a modulation in the amplitude and frequency power spectrum of oscillations. Moreover, nonlinear NMDAR action is seen to be state-dependent and can have opposite effects depending on the type of neurons involved and the level of input firing rate received. The presented model can serve as a computationally efficient building block in whole brain network models for investigating the differential modulation of different types of synapses under neuromodulatory influence or receptor specific malfunction.

Recent Progress in Brain Network Models for Medical Applications: A Review

Importance: Pathological perturbations of the brain often spread via connectome to fundamentally alter functional consequences. By integrating multimodal neuroimaging data with mathematical neural mass modeling, brain network models (BNMs) enable to quantitatively characterize aberrant network dynamics underlying multiple neurological and psychiatric disorders. We delved into the advancements of BNM-based medical applications, discussed the prevalent challenges within this field, and provided possible solutions and future directions. Highlights: This paper reviewed the theoretical foundations and current medical applications of computational BNMs. Composed of neural mass models, the BNM framework allows to investigate large-scale brain dynamics behind brain diseases by linking the simulated functional signals to the empirical neurophysiological data, and has shown promise in exploring neuropathological mechanisms, elucidating therapeutic effects, and predicting disease outcome. Despite that several limitations existed, one promising trend of this research field is to precisely guide clinical neuromodulation treatment based on individual BNM simulation. Conclusion: BNM carries the potential to help understand the mechanism underlying how neuropathology affects brain network dynamics, further contributing to decision-making in clinical diagnosis and treatment. Several constraints must be addressed and surmounted to pave the way for its utilization in the clinic.

Seizure Onset Zone Localization in Focal Epilepsy using stereo-EEG: SCSA and Visibility Graph complementary features

Despite how limiting the symptoms of Epilepsy are in daily life, over a third of epileptic conditions are drug-resistant. These seizures stem from a specific region of the brain: the Seizure Onset Zone (SOZ). The available treatment options include surgical resection or neurostimulation of the SOZ and both options require its accurate localization. In this scope, we aim to analyse stereo-electroencephalograph (s-EEG) data from a feature engineering perspective to propose reliable computer-aided SOZ localization algorithms. This paper investigates Semi-Classical Signal Analysis (SCSA) and Visibility Graphs (VG) features and establishes their complementarity. The study considers brain -scale signals simulated through the Epileptor model using The Virtual Brain (TVB). Results have been undeniably satisfactory and a valid contribution to Seizure Onset Zone localization.Clinical relevance- The paper contributes to the design of frameworks for assisting medical experts in Seizure Onset Zone localization in focal epilepsy.

Parameter estimation in a whole-brain network model of epilepsy: Comparison of parallel global optimization solvers

The Virtual Epileptic Patient (VEP) refers to a computer-based representation of a patient with epilepsy that combines personalized anatomical data with dynamical models of abnormal brain activities. It is capable of generating spatio-temporal seizure patterns that resemble those recorded with invasive methods such as stereoelectro EEG data, allowing for the evaluation of clinical hypotheses before planning surgery. This study highlights the effectiveness of calibrating VEP models using a global optimization approach. The approach utilizes SaCeSS, a cooperative metaheuristic algorithm capable of parallel computation, to yield high-quality solutions without requiring excessive computational time. Through extensive benchmarking on synthetic data, our proposal successfully solved a set of different configurations of VEP models, demonstrating better scalability and superior performance against other parallel solvers. These results were further enhanced using a Bayesian optimization framework for hyperparameter tuning, with significant gains in terms of both accuracy and computational cost. Additionally, we added a scalable uncertainty quantification phase after model calibration, and used it to assess the variability in estimated parameters across different problems. Overall, this study has the potential to improve the estimation of pathological brain areas in drug-resistant epilepsy, thereby to inform the clinical decision-making process.

Individualized epidemic spreading models predict epilepsy surgery outcomes: A pseudo-prospective study

Epilepsy surgery is the treatment of choice for drug-resistant epilepsy patients, but up to 50% of patients continue to have seizures one year after the resection. In order to aid presurgical planning and predict postsurgical outcome on a patient-by-patient basis, we developed a framework of individualized computational models that combines epidemic spreading with patient-specific connectivity and epileptogeneity maps: the Epidemic Spreading Seizure and Epilepsy Surgery framework (ESSES). ESSES parameters were fitted in a retrospective study (N = 15) to reproduce invasive electroencephalography (iEEG)-recorded seizures. ESSES reproduced the iEEG-recorded seizures, and significantly better so for patients with good (seizure-free, SF) than bad (nonseizure-free, NSF) outcome. We illustrate here the clinical applicability of ESSES with a pseudo-prospective study (N = 34) with a blind setting (to the resection strategy and surgical outcome) that emulated presurgical conditions. By setting the model parameters in the retrospective study, ESSES could be applied also to patients without iEEG data. ESSES could predict the chances of good outcome after any resection by finding patient-specific model-based optimal resection strategies, which we found to be smaller for SF than NSF patients, suggesting an intrinsic difference in the network organization or presurgical evaluation results of NSF patients. The actual surgical plan overlapped more with the model-based optimal resection, and had a larger effect in decreasing modeled seizure propagation, for SF patients than for NSF patients. Overall, ESSES could correctly predict 75% of NSF and 80.8% of SF cases pseudo-prospectively. Our results show that individualised computational models may inform surgical planning by suggesting alternative resections and providing information on the likelihood of a good outcome after a proposed resection. This is the first time that such a model is validated with a fully independent cohort and without the need for iEEG recordings.

Treatment effects in epilepsy: a mathematical framework for understanding response over time

Epilepsy is a neurological disorder characterized by recurrent seizures, affecting over 65 million people worldwide. Treatment typically commences with the use of anti-seizure medications, including both mono- and poly-therapy. Should these fail, more invasive therapies such as surgery, electrical stimulation and focal drug delivery are often considered in an attempt to render the person seizure free. Although a significant portion ultimately benefit from these treatment options, treatment responses often fluctuate over time. The physiological mechanisms underlying these temporal variations are poorly understood, making prognosis a significant challenge when treating epilepsy. Here we use a dynamic network model of seizure transition to understand how seizure propensity may vary over time as a consequence of changes in excitability. Through computer simulations, we explore the relationship between the impact of treatment on dynamic network properties and their vulnerability over time that permit a return to states of high seizure propensity. For small networks we show vulnerability can be fully characterised by the size of the first transitive component (FTC). For larger networks, we find measures of network efficiency, incoherence and heterogeneity (degree variance) correlate with robustness of networks to increasing excitability. These results provide a set of potential prognostic markers for therapeutic interventions in epilepsy. Such markers could be used to support the development of personalized treatment strategies, ultimately contributing to understanding of long-term seizure freedom.

Adolescent maturation of cortical excitation-inhibition balance based on individualized biophysical network modeling

The balance of excitation and inhibition is a key functional property of cortical microcircuits which changes through the lifespan. Adolescence is considered a crucial period for the maturation of excitation-inhibition balance. This has been primarily observed in animal studies, yet human in vivo evidence on adolescent maturation of the excitation-inhibition balance at the individual level is limited. Here, we developed an individualized in vivo marker of regional excitation-inhibition balance in human adolescents, estimated using large-scale simulations of biophysical network models fitted to resting-state functional magnetic resonance imaging data from two independent cross-sectional (N = 752) and longitudinal (N = 149) cohorts. We found a widespread relative increase of inhibition in association cortices paralleled by a relative age-related increase of excitation, or lack of change, in sensorimotor areas across both datasets. This developmental pattern co-aligned with multiscale markers of sensorimotor-association differentiation. The spatial pattern of excitation-inhibition development in adolescence was robust to inter-individual variability of structural connectomes and modeling configurations. Notably, we found that alternative simulation-based markers of excitation-inhibition balance show a variable sensitivity to maturational change. Taken together, our study highlights an increase of inhibition during adolescence in association areas using cross sectional and longitudinal data, and provides a robust computational framework to estimate microcircuit maturation in vivo at the individual level.

Virtual brain twins: from basic neuroscience to clinical use

Virtual brain twins are personalized, generative and adaptive brain models based on data from an individual's brain for scientific and clinical use. After a description of the key elements of virtual brain twins, we present the standard model for personalized whole-brain network models. The personalization is accomplished using a subject's brain imaging data by three means: (1) assemble cortical and subcortical areas in the subject-specific brain space; (2) directly map connectivity into the brain models, which can be generalized to other parameters; and (3) estimate relevant parameters through model inversion, typically using probabilistic machine learning. We present the use of personalized whole-brain network models in healthy ageing and five clinical diseases: epilepsy, Alzheimer's disease, multiple sclerosis, Parkinson's disease and psychiatric disorders. Specifically, we introduce spatial masks for relevant parameters and demonstrate their use based on the physiological and pathophysiological hypotheses. Finally, we pinpoint the key challenges and future directions.

Calibrating Deep Learning Classifiers for Patient-Independent Electroencephalogram Seizure Forecasting

The recent scientific literature abounds in proposals of seizure forecasting methods that exploit machine learning to automatically analyze electroencephalogram (EEG) signals. Deep learning algorithms seem to achieve a particularly remarkable performance, suggesting that the implementation of clinical devices for seizure prediction might be within reach. However, most of the research evaluated the robustness of automatic forecasting methods through randomized cross-validation techniques, while clinical applications require much more stringent validation based on patient-independent testing. In this study, we show that automatic seizure forecasting can be performed, to some extent, even on independent patients who have never been seen during the training phase, thanks to the implementation of a simple calibration pipeline that can fine-tune deep learning models, even on a single epileptic event recorded from a new patient. We evaluate our calibration procedure using two datasets containing EEG signals recorded from a large cohort of epileptic subjects, demonstrating that the forecast accuracy of deep learning methods can increase on average by more than 20%, and that performance improves systematically in all independent patients. We further show that our calibration procedure works best for deep learning models, but can also be successfully applied to machine learning algorithms based on engineered signal features. Although our method still requires at least one epileptic event per patient to calibrate the forecasting model, we conclude that focusing on realistic validation methods allows to more reliably compare different machine learning approaches for seizure prediction, enabling the implementation of robust and effective forecasting systems that can be used in daily healthcare practice.

Multi-scale modelling of the epileptic brain: advantages of computational therapy exploration

Objective: Epilepsy is a complex disease spanning across multiple scales, from ion channels in neurons to neuronal circuits across the entire brain. Over the past decades, computational models have been used to describe the pathophysiological activity of the epileptic brain from different aspects. Traditionally, each computational model can aid in optimizing therapeutic interventions, therefore, providing a particular view to design strategies for treating epilepsy. As a result, most studies are concerned with generating specific models of the epileptic brain that can help us understand the certain machinery of the pathological state. Those specific models vary in complexity and biological accuracy, with system-level models often lacking biological details.Approach: Here, we review various types of computational model of epilepsy and discuss their potential for different therapeutic approaches and scenarios, including drug discovery, surgical strategies, brain stimulation, and seizure prediction. We propose that we need to consider an integrated approach with a unified modelling framework across multiple scales to understand the epileptic brain. Our proposal is based on the recent increase in computational power, which has opened up the possibility of unifying those specific epileptic models into simulations with an unprecedented level of detail.Main results: A multi-scale epilepsy model can bridge the gap between biologically detailed models, used to address molecular and cellular questions, and brain-wide models based on abstract models which can account for complex neurological and behavioural observations.Significance: With these efforts, we move toward the next generation of epileptic brain models capable of connecting cellular features, such as ion channel properties, with standard clinical measures such as seizure severity.

Biophysical models applied to dementia patients reveal links between geographical origin, gender, disease duration, and loss of neural inhibition

BACKGROUND

The hypothesis of decreased neural inhibition in dementia has been sparsely studied in functional magnetic resonance imaging (fMRI) data across patients with different dementia subtypes, and the role of social and demographic heterogeneities on this hypothesis remains to be addressed.

METHODS

We inferred regional inhibition by fitting a biophysical whole-brain model (dynamic mean field model with realistic inter-areal connectivity) to fMRI data from 414 participants, including patients with Alzheimer's disease, behavioral variant frontotemporal dementia, and controls. We then investigated the effect of disease condition, and demographic and clinical variables on the local inhibitory feedback, a variable related to the maintenance of balanced neural excitation/inhibition.

RESULTS

Decreased local inhibitory feedback was inferred from the biophysical modeling results in dementia patients, specific to brain areas presenting neurodegeneration. This loss of local inhibition correlated positively with years with disease, and showed differences regarding the gender and geographical origin of the patients. The model correctly reproduced known disease-related changes in functional connectivity.

CONCLUSIONS

Results suggest a critical link between abnormal neural and circuit-level excitability levels, the loss of grey matter observed in dementia, and the reorganization of functional connectivity, while highlighting the sensitivity of the underlying biophysical mechanism to demographic and clinical heterogeneities in the patient population.

Linking fast and slow: The case for generative models

A pervasive challenge in neuroscience is testing whether neuronal connectivity changes over time due to specific causes, such as stimuli, events, or clinical interventions. Recent hardware innovations and falling data storage costs enable longer, more naturalistic neuronal recordings. The implicit opportunity for understanding the self-organised brain calls for new analysis methods that link temporal scales: from the order of milliseconds over which neuronal dynamics evolve, to the order of minutes, days, or even years over which experimental observations unfold. This review article demonstrates how hierarchical generative models and Bayesian inference help to characterise neuronal activity across different time scales. Crucially, these methods go beyond describing statistical associations among observations and enable inference about underlying mechanisms. We offer an overview of fundamental concepts in state-space modeling and suggest a taxonomy for these methods. Additionally, we introduce key mathematical principles that underscore a separation of temporal scales, such as the slaving principle, and review Bayesian methods that are being used to test hypotheses about the brain with multiscale data. We hope that this review will serve as a useful primer for experimental and computational neuroscientists on the state of the art and current directions of travel in the complex systems modelling literature.

A synergetic turn in cognitive neuroscience of brain diseases

Despite significant improvements in our understanding of brain diseases, many barriers remain. Cognitive neuroscience faces four major challenges: complex structure-function associations; disease phenotype heterogeneity; the lack of transdiagnostic models; and oversimplified cognitive approaches restricted to the laboratory. Here, we propose a synergetics framework that can help to perform the necessary dimensionality reduction of complex interactions between the brain, body, and environment. The key solutions include low-dimensional spatiotemporal hierarchies for brain-structure associations, whole-brain modeling to handle phenotype diversity, model integration of shared transdiagnostic pathophysiological pathways, and naturalistic frameworks balancing experimental control and ecological validity. Creating whole-brain models with reduced manifolds combined with ecological measures can improve our understanding of brain disease and help identify novel interventions. Synergetics provides an integrated framework for future progress in clinical and cognitive neuroscience, pushing the boundaries of brain health and disease toward more mature, naturalistic approaches.

Bifurcations and bursting in the Epileptor

The Epileptor is a phenomenological model for seizure activity that is used in a personalized large-scale brain modeling framework, the Virtual Epileptic Patient, with the aim of improving surgery outcomes for drug-resistant epileptic patients. Transitions between interictal and ictal states are modeled as bifurcations, enabling the definition of seizure classes in terms of onset/offset bifurcations. This establishes a taxonomy of seizures grounded in their essential underlying dynamics and the Epileptor replicates the activity of the most common class, as observed in patients with focal epilepsy, which is characterized by square-wave bursting properties. The Epileptor also encodes an additional mechanism to account for interictal spikes and spike and wave discharges. Here we use insights from a more generic model for square-wave bursting, based on the Unfolding Theory approach, to guide the bifurcation analysis of the Epileptor and gain a deeper understanding of the model and the role of its parameters. We show how the Epileptor's parameters can be modified to produce activities for other seizures classes of the taxonomy, as observed in patients, so that the large-scale brain models could be further personalized. Some of these classes have already been described in the literature in the Epileptor, others, predicted by the generic model, are new. Finally, we unveil how the interaction with the additional mechanism for spike and wave discharges alters the bifurcation structure of the main burster.

An exploratory computational analysis in mice brain networks of widespread epileptic seizure onset locations along with potential strategies for effective intervention and propagation control

Mean-field models have been developed to replicate key features of epileptic seizure dynamics. However, the precise mechanisms and the role of the brain area responsible for seizure onset and propagation remain incompletely understood. In this study, we employ computational methods within The Virtual Brain framework and the Epileptor model to explore how the location and connectivity of an Epileptogenic Zone (EZ) in a mouse brain are related to focal seizures (seizures that start in one brain area and may or may not remain localized), with a specific focus on the hippocampal region known for its association with epileptic seizures. We then devise computational strategies to confine seizures (prevent widespread propagation), simulating medical-like treatments such as tissue resection and the application of an anti-seizure drugs or neurostimulation to suppress hyperexcitability. Through selectively removing (blocking) specific connections informed by the structural connectome and graph network measurements or by locally reducing outgoing connection weights of EZ areas, we demonstrate that seizures can be kept constrained around the EZ region. We successfully identified the minimal connections necessary to prevent widespread seizures, with a particular focus on minimizing surgical or medical intervention while simultaneously preserving the original structural connectivity and maximizing brain functionality.

Learn how to interpret and use intracranial EEG findings

Epilepsy surgery is the therapy of choice for many patients with drug-resistant focal epilepsy. Recognizing and describing ictal and interictal patterns with intracranial electroencephalography (EEG) recordings is important in order to most efficiently leverage advantages of this technique to accurately delineate the seizure-onset zone before undergoing surgery. In this seminar in epileptology, we address learning objective "1.4.11 Recognize and describe ictal and interictal patterns with intracranial recordings" of the International League against Epilepsy curriculum for epileptologists. We will review principal considerations of the implantation planning, summarize the literature for the most relevant ictal and interictal EEG patterns within and beyond the Berger frequency spectrum, review invasive stimulation for seizure and functional mapping, discuss caveats in the interpretation of intracranial EEG findings, provide an overview on special considerations in children and in subdural grids/strips, and review available quantitative/signal analysis approaches. To be as practically oriented as possible, we will provide a mini atlas of the most frequent EEG patterns, highlight pearls for its not infrequently challenging interpretation, and conclude with two illustrative case examples. This article shall serve as a useful learning resource for trainees in clinical neurophysiology/epileptology by providing a basic understanding on the concepts of invasive intracranial EEG.

Intelligent digital tools for screening of brain connectivity and dementia risk estimation in people affected by mild cognitive impairment: the AI-Mind clinical study protocol

More than 10 million Europeans show signs of mild cognitive impairment (MCI), a transitional stage between normal brain aging and dementia stage memory disorder. The path MCI takes can be divergent; while some maintain stability or even revert to cognitive norms, alarmingly, up to half of the cases progress to dementia within 5 years. Current diagnostic practice lacks the necessary screening tools to identify those at risk of progression. The European patient experience often involves a long journey from the initial signs of MCI to the eventual diagnosis of dementia. The trajectory is far from ideal. Here, we introduce the AI-Mind project, a pioneering initiative with an innovative approach to early risk assessment through the implementation of advanced artificial intelligence (AI) on multimodal data. The cutting-edge AI-based tools developed in the project aim not only to accelerate the diagnostic process but also to deliver highly accurate predictions regarding an individual's risk of developing dementia when prevention and intervention may still be possible. AI-Mind is a European Research and Innovation Action (RIA H2020-SC1-BHC-06-2020, No. 964220) financed between 2021 and 2026. First, the AI-Mind Connector identifies dysfunctional brain networks based on high-density magneto- and electroencephalography (M/EEG) recordings. Second, the AI-Mind Predictor predicts dementia risk using data from the Connector, enriched with computerized cognitive tests, genetic and protein biomarkers, as well as sociodemographic and clinical variables. AI-Mind is integrated within a network of major European initiatives, including The Virtual Brain, The Virtual Epileptic Patient, and EBRAINS AISBL service for sensitive data, HealthDataCloud, where big patient data are generated for advancing digital and virtual twin technology development. AI-Mind's innovation lies not only in its early prediction of dementia risk, but it also enables a virtual laboratory scenario for hypothesis-driven personalized intervention research. This article introduces the background of the AI-Mind project and its clinical study protocol, setting the stage for future scientific contributions.

Modelling and Controlling System Dynamics of the Brain: An Intersection of Machine Learning and Control Theory

The human brain, as a complex system, has long captivated multidisciplinary researchers aiming to decode its intricate structure and function. This intricate network has driven scientific pursuits to advance our understanding of cognition, behavior, and neurological disorders by delving into the complex mechanisms underlying brain function and dysfunction. Modelling brain dynamics using machine learning techniques deepens our comprehension of brain dynamics from a computational perspective. These computational models allow researchers to simulate and analyze neural interactions, facilitating the identification of dysfunctions in connectivity or activity patterns. Additionally, the trained dynamical system, serving as a surrogate model, optimizes neurostimulation strategies under the guidelines of control theory. In this chapter, we discuss the recent studies on modelling and controlling brain dynamics at the intersection of machine learning and control theory, providing a framework to understand and improve cognitive function, and treat neurological and psychiatric disorders.

High-Density Exploration of Activity States in a Multi-Area Brain Model

To simulate whole brain dynamics with only a few equations, biophysical, mesoscopic models of local neuron populations can be connected using empirical tractography data. The development of mesoscopic mean-field models of neural populations, in particular, the Adaptive Exponential (AdEx mean-field model), has successfully summarized neuron-scale phenomena leading to the emergence of global brain dynamics associated with conscious (asynchronous and rapid dynamics) and unconscious (synchronized slow-waves, with Up-and-Down state dynamics) brain states, based on biophysical mechanisms operating at cellular scales (e.g. neuromodulatory regulation of spike-frequency adaptation during sleep-wake cycles or anesthetics). Using the Virtual Brain (TVB) environment to connect mean-field AdEx models, we have previously simulated the general properties of brain states, playing on spike-frequency adaptation, but have not yet performed detailed analyses of other parameters possibly also regulating transitions in brain-scale dynamics between different brain states. We performed a dense grid parameter exploration of the TVB-AdEx model, making use of High Performance Computing. We report a remarkable robustness of the effect of adaptation to induce synchronized slow-wave activity. Moreover, the occurrence of slow waves is often paralleled with a closer relation between functional and structural connectivity. We find that hyperpolarization can also generate unconscious-like synchronized Up and Down states, which may be a mechanism underlying the action of anesthetics. We conclude that the TVB-AdEx model reveals large-scale properties identified experimentally in sleep and anesthesia.

Enhanced simulations of whole-brain dynamics using hybrid resting-state structural connectomes

The human brain, composed of billions of neurons and synaptic connections, is an intricate network coordinating a sophisticated balance of excitatory and inhibitory activities between brain regions. The dynamical balance between excitation and inhibition is vital for adjusting neural input/output relationships in cortical networks and regulating the dynamic range of their responses to stimuli. To infer this balance using connectomics, we recently introduced a computational framework based on the Ising model, which was first developed to explain phase transitions in ferromagnets, and proposed a novel hybrid resting-state structural connectome (rsSC). Here, we show that a generative model based on the Kuramoto phase oscillator can be used to simulate static and dynamic functional connectomes (FC) with rsSC as the coupling weight coefficients, such that the simulated FC aligns well with the observed FC when compared with that simulated traditional structural connectome.

Altered spread of waves of activities at large scale is influenced by cortical thickness organization in temporal lobe epilepsy: a magnetic resonance imaging-high-density electroencephalography study

Temporal lobe epilepsy is a brain network disorder characterized by alterations at both the structural and the functional levels. It remains unclear how structure and function are related and whether this has any clinical relevance. In the present work, we adopted a novel methodological approach investigating how network structural features influence the large-scale dynamics. The functional network was defined by the spatio-temporal spreading of aperiodic bursts of activations (neuronal avalanches), as observed utilizing high-density electroencephalography in patients with temporal lobe epilepsy. The structural network was modelled as the region-based thickness covariance. Loosely speaking, we quantified the similarity of the cortical thickness of any two brain regions, both across groups and at the individual level, the latter utilizing a novel approach to define the subject-wise structural covariance network. In order to compare the structural and functional networks (at the nodal level), we studied the correlation between the probability that a wave of activity would propagate from a source to a target region and the similarity of the source region thickness as compared with other target brain regions. Building on the recent evidence that large-waves of activities pathologically spread through the epileptogenic network in temporal lobe epilepsy, also during resting state, we hypothesize that the structural cortical organization might influence such altered spatio-temporal dynamics. We observed a stable cluster of structure-function correlation in the bilateral limbic areas across subjects, highlighting group-specific features for left, right and bilateral temporal epilepsy. The involvement of contralateral areas was observed in unilateral temporal lobe epilepsy. We showed that in temporal lobe epilepsy, alterations of structural and functional networks pair in the regions where seizures propagate and are linked to disease severity. In this study, we leveraged on a well-defined model of neurological disease and pushed forward personalization approaches potentially useful in clinical practice. Finally, the methods developed here could be exploited to investigate the relationship between structure-function networks at subject level in other neurological conditions.

Network-motif delay differential analysis of brain activity during seizures

Delay Differential Analysis (DDA) is a nonlinear method for analyzing time series based on principles from nonlinear dynamical systems. DDA is extended here to incorporate network aspects to improve the dynamical characterization of complex systems. To demonstrate its effectiveness, DDA with network capabilities was first applied to the well-known Rössler system under different parameter regimes and noise conditions. Network-motif DDA, based on cortical regions, was then applied to invasive intracranial electroencephalographic data from drug-resistant epilepsy patients undergoing presurgical monitoring. The directional network motifs between brain areas that emerge from this analysis change dramatically before, during, and after seizures. Neural systems provide a rich source of complex data, arising from varying internal states generated by network interactions.

Cell-type specific and multiscale dynamics of human focal seizures in limbic structures

The relationship between clinically accessible epileptic biomarkers and neuronal activity underlying the transition to seizure is complex, potentially leading to imprecise delineation of epileptogenic brain areas. In particular, the pattern of interneuronal firing at seizure onset remains under debate, with some studies demonstrating increased firing and others suggesting reductions. Previous study of neocortical sites suggests that seizure recruitment occurs upon failure of inhibition, with intact feedforward inhibition in non-recruited territories. We investigated whether the same principle applies in limbic structures. We analysed simultaneous electrocorticography (ECoG) and neuronal recordings of 34 seizures in a cohort of 19 patients (10 male, 9 female) undergoing surgical evaluation for pharmacoresistant focal epilepsy. A clustering approach with five quantitative metrics computed from ECoG and multiunit data was used to distinguish three types of site-specific activity patterns during seizures, which at times co-existed within seizures. Overall, 156 single units were isolated, subclassified by cell-type and tracked through the seizure using our previously published methods to account for impacts of increased noise and single-unit waveshape changes caused by seizures. One cluster was closely associated with clinically defined seizure onset or spread. Entrainment of high-gamma activity to low-frequency ictal rhythms was the only metric that reliably identified this cluster at the level of individual seizures (P < 0.001). A second cluster demonstrated multi-unit characteristics resembling those in the first cluster, without concomitant high-gamma entrainment, suggesting feedforward effects from the seizure. The last cluster captured regions apparently unaffected by the ongoing seizure. Across all territories, the majority of both excitatory and inhibitory neurons reduced (69.2%) or ceased firing (21.8%). Transient increases in interneuronal firing rates were rare (13.5%) but showed evidence of intact feedforward inhibition, with maximal firing rate increases and waveshape deformations in territories not fully recruited but showing feedforward activity from the seizure, and a shift to burst-firing in seizure-recruited territories (P = 0.014). This study provides evidence for entrained high-gamma activity as an accurate biomarker of ictal recruitment in limbic structures. However, reduced neuronal firing suggested preserved inhibition in mesial temporal structures despite simultaneous indicators of seizure recruitment, in contrast to the inhibitory collapse scenario documented in neocortex. Further study is needed to determine if this activity is ubiquitous to hippocampal seizures or indicates a 'seizure-responsive' state in which the hippocampus is not the primary driver. If the latter, distinguishing such cases may help to refine the surgical treatment of mesial temporal lobe epilepsy.

Critical-like Brain Dynamics in a Continuum from Second- to First-Order Phase Transition

The classic brain criticality hypothesis postulates that the brain benefits from operating near a continuous second-order phase transition. Slow feedback regulation of neuronal activity could, however, lead to a discontinuous first-order transition and thereby bistable activity. Observations of bistability in awake brain activity have nonetheless remained scarce and its functional significance unclear. Moreover, there is no empirical evidence to support the hypothesis that the human brain could flexibly operate near either a first- or second-order phase transition despite such a continuum being common in models. Here, using computational modeling, we found bistable synchronization dynamics to emerge through elevated positive feedback and occur exclusively in a regimen of critical-like dynamics. We then assessed bistability in vivo with resting-state MEG in healthy adults (7 females, 11 males) and stereo-electroencephalography in epilepsy patients (28 females, 36 males). This analysis revealed that a large fraction of the neocortices exhibited varying degrees of bistability in neuronal oscillations from 3 to 200 Hz. In line with our modeling results, the neuronal bistability was positively correlated with classic assessment of brain criticality across narrow-band frequencies. Excessive bistability was predictive of epileptic pathophysiology in the patients, whereas moderate bistability was positively correlated with task performance in the healthy subjects. These empirical findings thus reveal the human brain as a one-of-a-kind complex system that exhibits critical-like dynamics in a continuum between continuous and discontinuous phase transitions.SIGNIFICANCE STATEMENT In the model, while synchrony per se was controlled by connectivity, increasing positive local feedback led to gradually emerging bistable synchrony with scale-free dynamics, suggesting a continuum between second- and first-order phase transitions in synchrony dynamics inside a critical-like regimen. In resting-state MEG and SEEG, bistability of ongoing neuronal oscillations was pervasive across brain areas and frequency bands and was observed only with concurring critical-like dynamics as the modeling predicted. As evidence for functional relevance, moderate bistability was positively correlated with executive functioning in the healthy subjects, and excessive bistability was associated with epileptic pathophysiology. These findings show that critical-like neuronal dynamics in vivo involves both continuous and discontinuous phase transitions in a frequency-, neuroanatomy-, and state-dependent manner.

Personalized dynamic network models of the human brain as a future tool for planning and optimizing epilepsy therapy

Epilepsy is a common neurological disorder, with one third of patients not responding to currently available antiepileptic drugs. The proportion of pharmacoresistant epilepsies has remained unchanged for many decades. To cure epilepsy and control seizures requires a paradigm shift in the development of new approaches to epilepsy diagnosis and treatment. Contemporary medicine has benefited from the exponential growth of computational modeling, and the application of network dynamics theory to understanding and treating human brain disorders. In epilepsy, the introduction of these approaches has led to personalized epileptic network modeling that can explore the patient's seizure genesis and predict the functional impact of resection on its individual network's propensity to seize. The application of the dynamic systems approach to neurostimulation therapy of epilepsy allows designing stimulation strategies that consider the patient's seizure dynamics and long-term fluctuations in the stability of their epileptic networks. In this article, we review, in a nontechnical fashion suitable for a broad neuroscientific audience, recent progress in personalized dynamic brain network modeling that is shaping the future approach to the diagnosis and treatment of epilepsy.

Computational modeling allows unsupervised classification of epileptic brain states across species

Current advances in epilepsy treatment aim to personalize and responsively adjust treatment parameters to overcome patient heterogeneity in treatment efficiency. For tailoring treatment to the individual and the current brain state, tools are required that help to identify the patient- and time-point-specific parameters of epilepsy. Computational modeling has long proven its utility in gaining mechanistic insight. Recently, the technique has been introduced as a diagnostic tool to predict individual treatment outcomes. In this article, the Wendling model, an established computational model of epilepsy dynamics, is used to automatically classify epileptic brain states in intracranial EEG from patients (n = 4) and local field potential recordings from in vitro rat data (high-potassium model of epilepsy, n = 3). Five-second signal segments are classified to four types of brain state in epilepsy (interictal, preonset, onset, ictal) by comparing a vector of signal features for each data segment to four prototypical feature vectors obtained by Wendling model simulations. The classification result is validated against expert visual assessment. Model-driven brain state classification achieved a classification performance significantly above chance level (mean sensitivity 0.99 on model data, 0.77 on rat data, 0.56 on human data in a four-way classification task). Model-driven prototypes showed similarity with data-driven prototypes, which we obtained from real data for rats and humans. Our results indicate similar electrophysiological patterns of epileptic states in the human brain and the animal model that are well-reproduced by the computational model, and captured by a key set of signal features, enabling fully automated and unsupervised brain state classification in epilepsy.

Virtual brain simulations reveal network-specific parameters in neurodegenerative dementias

Introduction

Neural circuit alterations lay at the core of brain physiopathology, and yet are hard to unveil in living subjects. The Virtual Brain (TVB) modeling, by exploiting structural and functional magnetic resonance imaging (MRI), yields mesoscopic parameters of connectivity and synaptic transmission.

Methods

We used TVB to simulate brain networks, which are key for human brain function, in Alzheimer's disease (AD) and frontotemporal dementia (FTD) patients, whose connectivity and synaptic parameters remain largely unknown; we then compared them to healthy controls, to reveal novel in vivo pathological hallmarks.

Results

The pattern of simulated parameter differed between AD and FTD, shedding light on disease-specific alterations in brain networks. Individual subjects displayed subtle differences in network parameter patterns that significantly correlated with their individual neuropsychological, clinical, and pharmacological profiles.

Discussion

These TVB simulations, by informing about a new personalized set of networks parameters, open new perspectives for understanding dementias mechanisms and design personalized therapeutic approaches.

Altered Relationship between Functional Connectivity and Fiber-Bundle Structure in High-Functioning Male Adults with Autism Spectrum Disorder

Autism spectrum disorder (ASD) is a pervasive neurodevelopmental disorder characterized by abnormalities in structure and function of the brain. However, how ASD affects the relationship between fiber-bundle microstructures and functional connectivity (FC) remains unclear. Here, we analyzed structural and functional images of 26 high-functioning adult males with ASD, alongside 26 age-, gender-, and full-scale IQ-matched typically developing controls (TDCs) from the BNI dataset in the ABIDE database. We utilized fixel-based analysis to extract microstructural information from fiber tracts, which was then used to predict FC using a multilinear model. Our results revealed that the structure-function relationships in both ASD and TDC cohorts were strongly aligned in the primary cortex but decoupled in the high-order cortex, and the ASD patients exhibited reduced structure-function relationships throughout the cortex compared to the TDCs. Furthermore, we observed that the disrupted relationships in ASD were primarily driven by alterations in FC rather than fiber-bundle microstructures. The structure-function relationships in the left superior parietal cortex, right precentral and inferior temporal cortices, and bilateral insula could predict individual differences in clinical symptoms of ASD patients. These findings underscore the significance of altered relationships between fiber-bundle microstructures and FC in the etiology of ASD.

Observers for Phenomenological Models of Epileptic Seizures

Progress towards effective treatment of epileptic seizures has seen much improvement in the past decade. In particular, the emergence of phenomenological models of epileptic seizures specifically designed to capture the electrical seizure dynamics in the Epileptor model is inspiring new approaches to predicting and controlling seizures. These new models present in various forms and contain important but unmeasurable variables that control the occurrence of seizures. These models have been used mostly as nodes in large networks to study the complex brain behaviour of seizures. In order to use this model for the purposes of seizure forecasting or to control seizures through deep brain stimulation, the states of the model will need to be estimated. Although devices such as EEG electrodes can be related to some of the states of the model, most remain unmeasured and would require an observer (as defined in control theory) for their estimation. Additionally, we would like to consider the case for large nodes of systems where the number of electrodes is far smaller than the number of nodes being estimated. In this paper, we provide methods towards obtaining the full states of these phenomenological models using nonlinear observers. In particular, we explore the effectiveness of the Extended Kalman Filter for small networks of nodes of a smoothed sixth order Epileptor model. We show that observer design is possible for this family of systems and identify the difficulties in doing so.Clinical relevance-The methods presented here can be applied with an individual epileptic patient's EEG to reveal previously hidden biomarkers of epilepsy for seizure forecasting.

Subgenual cingulate connectivity as a treatment predictor during low-frequency right dorsolateral prefrontal rTMS: A concurrent TMS-fMRI study

INTRODUCTION

Repetitive transcranial magnetic stimulation (rTMS) to the dorsolateral prefrontal cortex (DLPFC) is effective in alleviating treatment-resistant depression (TRD). It has been proposed that regions within the left DLPFC that are anti-correlated with the right subgenual anterior cingulate cortex (sgACC) may represent optimal individualized target sites for high-frequency left rTMS (HFL).

OBJECTIVE/HYPOTHESIS

This study aimed to explore the effects of low-frequency right rTMS (LFR) on left sgACC connectivity during concurrent TMS-fMRI.

METHODS

34 TRD patients underwent an imaging session that included both a resting-state fMRI run (rs-fMRI0) and a run during which LFR was applied to the right DLPFC (TMS-fMRI). Participants subsequently completed four weeks of LFR treatment. The left sgACC functional connectivity was compared between the rs-fMRI0 run and TMS-fMRI run. Personalized e-fields and a region-of-interest approach were used to calculate overlap of left sgACC functional connectivity at the TMS target and to assess for a relationship with treatment effects.

RESULTS

TMS-fMRI increased left sgACC functional connectivity to parietal regions within the ventral attention network; differences were not significantly associated with clinical improvements. Personalized e-fields were not significant in predicting treatment outcomes (p = 0.18).

CONCLUSION

This was the first study to examine left sgACC anti-correlation with the right DLPFC during an LFR rTMS protocol. In contrast to studies that targeted the left DLPFC, we did not find that higher anti-correlation was associated with clinical outcomes. Our results suggest that the antidepressant mechanism of action of LFR to the right DLPFC may be different than for HFL.

Low-dimensional organization of global brain states of reduced consciousness

Brain states are frequently represented using a unidimensional scale measuring the richness of subjective experience (level of consciousness). This description assumes a mapping between the high-dimensional space of whole-brain configurations and the trajectories of brain states associated with changes in consciousness, yet this mapping and its properties remain unclear. We combine whole-brain modeling, data augmentation, and deep learning for dimensionality reduction to determine a mapping representing states of consciousness in a low-dimensional space, where distances parallel similarities between states. An orderly trajectory from wakefulness to patients with brain injury is revealed in a latent space whose coordinates represent metrics related to functional modularity and structure-function coupling, increasing alongside loss of consciousness. Finally, we investigate the effects of model perturbations, providing geometrical interpretation for the stability and reversibility of states. We conclude that conscious awareness depends on functional patterns encoded as a low-dimensional trajectory within the vast space of brain configurations.