Whole-brain Computational Modeling Research Articles

Turbulent dynamics and whole-brain modeling: toward new clinical applications for traumatic brain injury.

Traumatic Brain Injury (TBI) is a prevalent disorder mostly characterized by persistent impairments in cognitive function that poses a substantial burden on caregivers and the healthcare system worldwide. Crucially, severity classification is primarily based on clinical evaluations, which are non-specific and poorly predictive of long-term disability. In this Mini Review, we first provide a description of our model-free and model-based approaches within the turbulent dynamics framework as well as our vision on how they can potentially contribute to provide new neuroimaging biomarkers for TBI. In addition, we report the main findings of our recent study examining longitudinal changes in moderate-severe TBI (msTBI) patients during a one year spontaneous recovery by applying the turbulent dynamics framework (model-free approach) and the Hopf whole-brain computational model (model-based approach) combined with in silico perturbations. Given the neuroinflammatory response and heightened risk for neurodegeneration after TBI, we also offer future directions to explore the association with genomic information. Moreover, we discuss how whole-brain computational modeling may advance our understanding of the impact of structural disconnection on whole-brain dynamics after msTBI in light of our recent findings. Lastly, we suggest future avenues whereby whole-brain computational modeling may assist the identification of optimal brain targets for deep brain stimulation to promote TBI recovery.

Task-relevant brain dynamics among cognitive subsystems induced by regional stimulation in a whole-brain computational model.

Cognition involves the global integration of distributed brain regions that are known to work cohesively as cognitive subsystems during brain functioning. Empirical evidence has suggested that spatiotemporal phase relationships between brain regions, measured as synchronization and metastability, may encode important task-relevant information. However, it remains largely unknown how phase relationships aggregate at the level of cognitive subsystems under different cognitive processing. Here, we probe this question by simulating task-relevant brain dynamics through regional stimulation of a whole-brain dynamical network model operating in the resting-state dynamical regime. The model is constructed with structurally embedded Stuart-Laudon oscillators and then fitted with human resting-state functional magnetic resonance imaging data. Based on this framework, we first demonstrate the plausibility of introducing the cognitive system partition into the modeling analysis framework by showing that the clustering of regions across functional networks is better circumscribed by the predefined partition. At the cognitive subsystem level, we focus on how task-relevant phase dynamics are organized in terms of synchronization and metastability. We found that patterns of cognitive synchronization are more task specific, whereas patterns of cognitive metastability are more consistent across different states, suggesting it may encode a more task-general property during cognitive processing, an inherent property conferred by brain organization. This consistent network architecture in cognitive metastability may be related to the distinct functional responses of realistic cognitive systems. We also provide empirical evidence to partially support our computational results. Our paper may provide insights for the mechanisms underlying task-relevant brain dynamics, and establish a model-based link between brain structure, dynamics, and cognition, a fundamental step for computationally aided brain interventions.

Neural fingerprints of gambling disorder using diffusion tensor imaging

Gambling disorder (GD) is a behavioral addiction associated with personal, social and occupational consequences. Thus, examining GD's clinical relationship with its neural substrates is critical. We compared neural fingerprints using diffusion tensor imaging (DTI) in GD subjects undergoing treatment relative to healthy volunteers (HV). Fifty-three (25 GD, 28 age-matched HV) males were scanned with structural magnetic resonance imaging (MRI) and DTI. We applied probabilistic tractography based on DTI scanning data, preprocessed and analyzed using permutation testing of individual connectivity weights between regions for group comparison. Permutation-based comparisons between group-averaged connectomes highlighted significant structural differences. The GD group demonstrated increased connectivity, and striatal network reorganisation, contrasted by reduced connectivity within and to frontal lobe nodes. Modularity analysis revealed that the GD group had fewer hubs integrating information across the brain. We highlight GD neural changes involved in controlling risk-seeking behaviors. The observed striatal restructuring converges with previous research, and the increased connectivity affects subnetworks highly active in gambling situations, although these findings are not significant when correcting for multiple comparisons. Modularity analysis underlines that, despite connectivity increases, the GD connectome loses hubs, impeding its neuronal network coherence. Together, these results demonstrate the feasibility of using whole-brain computational modeling in assessing GD.

Strength-dependent perturbation of whole-brain model working in different regimes reveals the role of fluctuations in brain dynamics.

Despite decades of research, there is still a lack of understanding of the role and generating mechanisms of the ubiquitous fluctuations and oscillations found in recordings of brain dynamics. Here, we used whole-brain computational models capable of presenting different dynamical regimes to reproduce empirical data's turbulence level. We showed that the model's fluctuations regime fitted to turbulence more faithfully reproduces the empirical functional connectivity compared to oscillatory and noise regimes. By applying global and local strength-dependent perturbations and subsequently measuring the responsiveness of the model, we revealed each regime's computational capacity demonstrating that brain dynamics is shifted towards fluctuations to provide much-needed flexibility. Importantly, fluctuation regime stimulation in a brain region within a given resting state network modulates that network, aligned with previous empirical and computational studies. Furthermore, this framework generates specific, testable empirical predictions for human stimulation studies using strength-dependent rather than constant perturbation. Overall, the whole-brain models fitted to the level of empirical turbulence together with functional connectivity unveil that the fluctuation regime best captures empirical data, and the strength-dependent perturbative framework demonstrates how this regime provides maximal flexibility to the human brain.

High-order functional redundancy in ageing explained via alterations in the connectome in a whole-brain model.

The human brain generates a rich repertoire of spatio-temporal activity patterns, which support a wide variety of motor and cognitive functions. These patterns of activity change with age in a multi-factorial manner. One of these factors is the variations in the brain's connectomics that occurs along the lifespan. However, the precise relationship between high-order functional interactions and connnectomics, as well as their variations with age are largely unknown, in part due to the absence of mechanistic models that can efficiently map brain connnectomics to functional connectivity in aging. To investigate this issue, we have built a neurobiologically-realistic whole-brain computational model using both anatomical and functional MRI data from 161 participants ranging from 10 to 80 years old. We show that the differences in high-order functional interactions between age groups can be largely explained by variations in the connectome. Based on this finding, we propose a simple neurodegeneration model that is representative of normal physiological aging. As such, when applied to connectomes of young participant it reproduces the age-variations that occur in the high-order structure of the functional data. Overall, these results begin to disentangle the mechanisms by which structural changes in the connectome lead to functional differences in the ageing brain. Our model can also serve as a starting point for modeling more complex forms of pathological ageing or cognitive deficits.

A multiscale brain network model links Alzheimer’s disease-mediated neuronal hyperactivity to large-scale oscillatory slowing

BackgroundNeuronal hyperexcitability and inhibitory interneuron dysfunction are frequently observed in preclinical animal models of Alzheimer’s disease (AD). This study investigates whether these microscale abnormalities explain characteristic large-scale magnetoencephalography (MEG) activity in human early-stage AD patients.MethodsTo simulate spontaneous electrophysiological activity, we used a whole-brain computational network model comprised of 78 neural masses coupled according to human structural brain topology. We modified relevant model parameters to simulate six literature-based cellular scenarios of AD and compare them to one healthy and six contrast (non-AD-like) scenarios. The parameters include excitability, postsynaptic potentials, and coupling strength of excitatory and inhibitory neuronal populations. Whole-brain spike density and spectral power analyses of the simulated data reveal mechanisms of neuronal hyperactivity that lead to oscillatory changes similar to those observed in MEG data of 18 human prodromal AD patients compared to 18 age-matched subjects with subjective cognitive decline.ResultsAll but one of the AD-like scenarios showed higher spike density levels, and all but one of these scenarios had a lower peak frequency, higher spectral power in slower (theta, 4–8Hz) frequencies, and greater total power. Non-AD-like scenarios showed opposite patterns mainly, including reduced spike density and faster oscillatory activity. Human AD patients showed oscillatory slowing (i.e., higher relative power in the theta band mainly), a trend for lower peak frequency and higher total power compared to controls. Combining model and human data, the findings indicate that neuronal hyperactivity can lead to oscillatory slowing, likely due to hyperexcitation (by hyperexcitability of pyramidal neurons or greater long-range excitatory coupling) and/or disinhibition (by reduced excitability of inhibitory interneurons or weaker local inhibitory coupling strength) in early AD.ConclusionsUsing a computational brain network model, we link findings from different scales and models and support the hypothesis of early-stage neuronal hyperactivity underlying E/I imbalance and whole-brain network dysfunction in prodromal AD.

Effective connectivity extracts clinically relevant prognostic information from resting state activity in stroke.

Recent resting-state functional MRI studies in stroke patients have identified two robust biomarkers of acute brain dysfunction: a reduction of inter-hemispheric functional connectivity between homotopic regions of the same network, and an abnormal increase of ipsi-lesional functional connectivity between task-negative and task-positive resting-state networks. Whole-brain computational modelling studies, at the individual subject level, using undirected effective connectivity derived from empirically measured functional connectivity, have shown a reduction of measures of integration and segregation in stroke as compared to healthy brains. Here we employ a novel method, first, to infer whole-brain directional effective connectivity from zero-lagged and lagged covariance matrices, then, to compare it to empirically measured functional connectivity for predicting stroke versus healthy status, and patient performance (zero, one, multiple deficits) across neuropsychological tests. We also investigated the accuracy of functional connectivity versus model effective connectivity in predicting the long-term outcome from acute measures. Both functional and effective connectivity predicted healthy from stroke individuals significantly better than the chance-level; however, accuracy for the effective connectivity was significantly higher than for functional connectivity at 1- to 2-week, 3-month and 1-year post-stroke. Predictive functional connections mainly included those reported in previous studies (within-network inter-hemispheric and between task-positive and -negative networks intra-hemispherically). Predictive effective connections included additional between-network links. Effective connectivity was a better predictor than functional connectivity of the number of behavioural domains in which patients suffered deficits, both at 2-week and 1-year post-onset of stroke. Interestingly, patient deficits at 1-year time-point were better predicted by effective connectivity values at 2 weeks rather than at 1-year time-point. Our results thus demonstrate that the second-order statistics of functional MRI resting-state activity at an early stage of stroke, derived from a whole-brain effective connectivity, estimated in a model fitted to reproduce the propagation of neuronal activity, has pertinent information for clinical prognosis.

Perturbations in dynamical models of whole-brain activity dissociate between the level and stability of consciousness.

Consciousness transiently fades away during deep sleep, more stably under anesthesia, and sometimes permanently due to brain injury. The development of an index to quantify the level of consciousness across these different states is regarded as a key problem both in basic and clinical neuroscience. We argue that this problem is ill-defined since such an index would not exhaust all the relevant information about a given state of consciousness. While the level of consciousness can be taken to describe the actual brain state, a complete characterization should also include its potential behavior against external perturbations. We developed and analyzed whole-brain computational models to show that the stability of conscious states provides information complementary to their similarity to conscious wakefulness. Our work leads to a novel methodological framework to sort out different brain states by their stability and reversibility, and illustrates its usefulness to dissociate between physiological (sleep), pathological (brain-injured patients), and pharmacologically-induced (anesthesia) loss of consciousness.

Cholinergic neuromodulation of inhibitory interneurons facilitates functional integration in whole-brain models.

Segregation and integration are two fundamental principles of brain structural and functional organization. Neuroimaging studies have shown that the brain transits between different functionally segregated and integrated states, and neuromodulatory systems have been proposed as key to facilitate these transitions. Although whole-brain computational models have reproduced this neuromodulatory effect, the role of local inhibitory circuits and their cholinergic modulation has not been studied. In this article, we consider a Jansen & Rit whole-brain model in a network interconnected using a human connectome, and study the influence of the cholinergic and noradrenergic neuromodulatory systems on the segregation/integration balance. In our model, we introduce a local inhibitory feedback as a plausible biophysical mechanism that enables the integration of whole-brain activity, and that interacts with the other neuromodulatory influences to facilitate the transition between different functional segregation/integration regimes in the brain.

Data-driven method to infer the seizure propagation patterns in an epileptic brain from intracranial electroencephalography.

Surgical interventions in epileptic patients aimed at the removal of the epileptogenic zone have success rates at only 60-70%. This failure can be partly attributed to the insufficient spatial sampling by the implanted intracranial electrodes during the clinical evaluation, leading to an incomplete picture of spatio-temporal seizure organization in the regions that are not directly observed. Utilizing the partial observations of the seizure spreading through the brain network, complemented by the assumption that the epileptic seizures spread along the structural connections, we infer if and when are the unobserved regions recruited in the seizure. To this end we introduce a data-driven model of seizure recruitment and propagation across a weighted network, which we invert using the Bayesian inference framework. Using a leave-one-out cross-validation scheme on a cohort of 45 patients we demonstrate that the method can improve the predictions of the states of the unobserved regions compared to an empirical estimate that does not use the structural information, yet it is on the same level as the estimate that takes the structure into account. Furthermore, a comparison with the performed surgical resection and the surgery outcome indicates a link between the inferred excitable regions and the actual epileptogenic zone. The results emphasize the importance of the structural connectome in the large-scale spatio-temporal organization of epileptic seizures and introduce a novel way to integrate the patient-specific connectome and intracranial seizure recordings in a whole-brain computational model of seizure spread.

Increased sensitivity to strong perturbations in a whole-brain model of LSD

Lysergic acid diethylamide (LSD) is a potent psychedelic drug, which has seen a revival in clinical and pharmacological research within recent years. Human neuroimaging studies have shown fundamental changes in brain-wide functional connectivity and an expansion of dynamical brain states, thus raising the question about a mechanistic explanation of the dynamics underlying these alterations. Here, we applied a novel perturbational approach based on a whole-brain computational model, which opens up the possibility to externally perturb different brain regions in silico and investigate differences in dynamical stability of different brain states, i.e. the dynamical response of a certain brain region to an external perturbation. After adjusting the whole-brain model parameters to reflect the dynamics of functional magnetic resonance imaging (fMRI) BOLD signals recorded under the influence of LSD or placebo, perturbations of different brain areas were simulated by either promoting or disrupting synchronization in the regarding brain region. After perturbation offset, we quantified the recovery characteristics of the brain area to its basal dynamical state with the Perturbational Integration Latency Index (PILI) and used this measure to distinguish between the two brain states. We found significant changes in dynamical complexity with consistently higher PILI values after LSD intake on a global level, which indicates a shift of the brain's global working point further away from a stable equilibrium as compared to normal conditions. On a local level, we found that the largest differences were measured within the limbic network, the visual network and the default mode network. Additionally, we found a higher variability of PILI values across different brain regions after LSD intake, indicating higher response diversity under LSD after an external perturbation. Our results provide important new insights into the brain-wide dynamical changes underlying the psychedelic state - here provoked by LSD intake - and underline possible future clinical applications of psychedelic drugs in particular psychiatric disorders.

Signature of consciousness in brain-wide synchronization patterns of monkey and human fMRI signals

During the sleep-wake cycle, the brain undergoes profound dynamical changes, which manifest subjectively as transitions between conscious experience and unconsciousness. Yet, neurophysiological signatures that can objectively distinguish different consciousness states based are scarce. Here, we show that differences in the level of brain-wide signals can reliably distinguish different stages of sleep and anesthesia from the awake state in human and monkey fMRI resting state data. Moreover, a whole-brain computational model can faithfully reproduce changes in global synchronization and other metrics such as functional connectivity, structure-function relationship, integration and segregation across vigilance states. We demonstrate that the awake brain is close to a Hopf bifurcation, which naturally coincides with the emergence of globally correlated fMRI signals. Furthermore, simulating lesions of individual brain areas highlights the importance of connectivity hubs in the posterior brain and subcortical nuclei for maintaining the model in the awake state, as predicted by graph-theoretical analyses of structural data.

Whole-brain computational modeling reveals disruption of microscale brain dynamics in HIV infected individuals.

MRI‐based neuroimaging techniques have been used to investigate brain injury associated with HIV‐infection. Whole‐brain cortical mean‐field dynamic modeling provides a way to integrate structural and functional imaging outcomes, allowing investigation of microscale brain dynamics. In this study, we adopted the relaxed mean‐field dynamic modeling to investigate structural and functional connectivity in 42 HIV‐infected subjects before and after 12‐week of combination antiretroviral therapy (cART) and compared them with 46 age‐matched healthy subjects. Microscale brain dynamics were modeled by a set of parameters including two region‐specific microscale brain properties, recurrent connection strengths, and subcortical inputs. We also analyzed the relationship between the model parameters (i.e., the recurrent connection and subcortical inputs) and functional network topological characterizations, including smallworldness, clustering coefficient, and network efficiency. The results show that untreated HIV‐infected individuals have disrupted local brain dynamics that in part correlate with network topological measurements. Notably, after 12 weeks of cART, both the microscale brain dynamics and the network topological measurements improved and were closer to those in the healthy brain. This was also associated with improved cognitive performance, suggesting that improvement in local brain dynamics translates into clinical improvement.

Data augmentation based on dynamical systems for the classification of brain states

Abstract The application of machine learning algorithms to neuroimaging data shows great promise for the classification of physiological and pathological brain states. However, classifiers trained on high dimensional data are prone to overfitting, especially for a low number of training samples. We describe the use of whole-brain computational models for data augmentation in brain state classification. Our low dimensional model is based on nonlinear oscillators coupled by the empirical structural connectivity of the brain. We use this model to enhance a dataset consisting of functional magnetic resonance imaging recordings acquired during all stages of the human wake-sleep cycle. After fitting the model to the average functional connectivity of each state, we show that the synthetic data generated by the model yields classification accuracies comparable to those obtained from the empirical data. We also show that models fitted to individual subjects generate surrogates with enough information to train classifiers that present significant transfer learning accuracy to the whole sample. Whole-brain computational modeling represents a useful tool to produce large synthetic datasets for data augmentation in the classification of certain brain states, with potential applications to computer-assisted diagnosis and prognosis of neuropsychiatric disorders.

Adaptive frequency-based modeling of whole-brain oscillations: Predicting regional vulnerability and hazardousness rates.

Whole-brain computational modeling based on structural connectivity has shown great promise in successfully simulating fMRI BOLD signals with temporal coactivation patterns that are highly similar to empirical functional connectivity patterns during resting state. Importantly, previous studies have shown that spontaneous fluctuations in coactivation patterns of distributed brain regions have an inherent dynamic nature with regard to the frequency spectrum of intrinsic brain oscillations. In this modeling study, we introduced frequency dynamics into a system of coupled oscillators, where each oscillator represents the local mean-field model of a brain region. We first showed that the collective behavior of interacting oscillators reproduces previously shown features of brain dynamics. Second, we examined the effect of simulated lesions in gray matter by applying an in silico perturbation protocol to the brain model. We present a new approach to map the effects of vulnerability in brain networks and introduce a measure of regional hazardousness based on mapping of the degree of divergence in a feature space.

Perturbation of whole-brain dynamics in silico reveals mechanistic differences between brain states

Human neuroimaging research has revealed that wakefulness and sleep involve very different activity patterns. Yet, it is not clear why brain states differ in their dynamical complexity, e.g. in the level of integration and segregation across brain networks over time. Here, we investigate the mechanisms underlying the dynamical stability of brain states using a novel off-line in silico perturbation protocol. We first adjust a whole-brain computational model to the basal dynamics of wakefulness and deep sleep recorded with fMRI in two independent human fMRI datasets. Then, the models of sleep and awake brain states are perturbed using two distinct multifocal protocols either promoting or disrupting synchronization in randomly selected brain areas. Once perturbation is halted, we use a novel measure, the Perturbative Integration Latency Index (PILI), to evaluate the recovery back to baseline. We find a clear distinction between models, consistently showing larger PILI in wakefulness than in deep sleep, corroborating previous experimental findings. In the models, larger recoveries are associated to a critical slowing down induced by a shift in the model's operation point, indicating that the awake brain operates further from a stable equilibrium than deep sleep. This novel approach opens up for a new level of artificial perturbative studies unconstrained by ethical limitations allowing for a deeper investigation of the dynamical properties of different brain states.

Effect of Field Spread on Resting-State Magneto Encephalography Functional Network Analysis: A Computational Modeling Study.

A popular way to analyze resting-state electroencephalography (EEG) and magneto encephalography (MEG) data is to treat them as a functional network in which sensors are identified with nodes and the interaction between channel time series and the network connections. Although conceptually appealing, the network-theoretical approach to sensor-level EEG and MEG data is challenged by the fact that EEG and MEG time series are mixtures of source activity. It is, therefore, of interest to assess the relationship between functional networks of source activity and the ensuing sensor-level networks. Since these topological features are of high interest in experimental studies, we address the question of to what extent the network topology can be reconstructed from sensor-level functional connectivity (FC) measures in case of MEG data. Simple simulations that consider only a small number of regions do not allow to assess network properties; therefore, we use a diffusion magnetic resonance imaging-constrained whole-brain computational model of resting-state activity. Our motivation lies behind the fact that still many contributions found in the literature perform network analysis at sensor level, and we aim at showing the discrepancies between source- and sensor-level network topologies by using realistic simulations of resting-state cortical activity. Our main findings are that the effect of field spread on network topology depends on the type of interaction (instantaneous or lagged) and leads to an underestimation of lagged FC at sensor level due to instantaneous mixing of cortical signals, instantaneous interaction is more sensitive to field spread than lagged interaction, and discrepancies are reduced when using planar gradiometers rather than axial gradiometers. We, therefore, recommend using lagged interaction measures on planar gradiometer data when investigating network properties of resting-state sensor-level MEG data.

Novel Intrinsic Ignition Method Measuring Local-Global Integration Characterizes Wakefulness and Deep Sleep.

A precise definition of a brain state has proven elusive. Here, we introduce the novel local-global concept of intrinsic ignition characterizing the dynamical complexity of different brain states. Naturally occurring intrinsic ignition events reflect the capability of a given brain area to propagate neuronal activity to other regions, giving rise to different levels of integration. The ignitory capability of brain regions is computed by the elicited level of integration for each intrinsic ignition event in each brain region, averaged over all events. This intrinsic ignition method is shown to clearly distinguish human neuroimaging data of two fundamental brain states (wakefulness and deep sleep). Importantly, whole-brain computational modelling of this data shows that at the optimal working point is found where there is maximal variability of the intrinsic ignition across brain regions. Thus, combining whole brain models with intrinsic ignition can provide novel insights into underlying mechanisms of brain states.

Uncovering the underlying mechanisms and whole-brain dynamics of deep brain stimulation for Parkinson\u2019s disease

Deep brain stimulation (DBS) for Parkinson’s disease is a highly effective treatment in controlling otherwise debilitating symptoms. Yet the underlying brain mechanisms are currently not well understood. Whole-brain computational modeling was used to disclose the effects of DBS during resting-state functional Magnetic Resonance Imaging in ten patients with Parkinson’s disease. Specifically, we explored the local and global impact that DBS has in creating asynchronous, stable or critical oscillatory conditions using a supercritical bifurcation model. We found that DBS shifts global brain dynamics of patients towards a Healthy regime. This effect was more pronounced in very specific brain areas such as the thalamus, globus pallidus and orbitofrontal regions of the right hemisphere (with the left hemisphere not analyzed given artifacts arising from the electrode lead). Global aspects of integration and synchronization were also rebalanced. Empirically, we found higher communicability and coherence brain measures during DBS-ON compared to DBS-OFF. Finally, using our model as a framework, artificial in silico DBS was applied to find potential alternative target areas for stimulation and whole-brain rebalancing. These results offer important insights into the underlying large-scale effects of DBS as well as in finding novel stimulation targets, which may offer a route to more efficacious treatments.

Increased Stability and Breakdown of Brain Effective Connectivity During Slow-Wave Sleep: Mechanistic Insights from Whole-Brain Computational Modelling

Recent research has found that the human sleep cycle is characterised by changes in spatiotemporal patterns of brain activity. Yet, we are still missing a mechanistic explanation of the local neuronal dynamics underlying these changes. We used whole-brain computational modelling to study the differences in global brain functional connectivity and synchrony of fMRI activity in healthy humans during wakefulness and slow-wave sleep. We applied a whole-brain model based on the normal form of a supercritical Hopf bifurcation and studied the dynamical changes when adapting the bifurcation parameter for all brain nodes to best match wakefulness and slow-wave sleep. Furthermore, we analysed differences in effective connectivity between the two states. In addition to significant changes in functional connectivity, synchrony and metastability, this analysis revealed a significant shift of the global dynamic working point of brain dynamics, from the edge of the transition between damped to sustained oscillations during wakefulness, to a stable focus during slow-wave sleep. Moreover, we identified a significant global decrease in effective interactions during slow-wave sleep. These results suggest a mechanism for the empirical functional changes observed during slow-wave sleep, namely a global shift of the brain’s dynamic working point leading to increased stability and decreased effective connectivity.