Whole-Brain Network Models: From Physics to Bedside

A biologically inspired computational model of human ventral temporal cortex

Our minds represent miscellaneous objects in the physical world metaphorically in an abstract and complex high-dimensional object space, which is implemented in a two-dimensional surface of the ventral temporal cortex (VTC) with topologically organized object selectivity. Here we investigated principles guiding the topographical organization of object selectivities in the VTC by constructing a hybrid Self-Organizing Map (SOM) model that harnesses a biologically inspired algorithm of wiring cost minimization and adheres to the constraints of the lateral wiring span of human VTC neurons. In a series of in silico experiments with functional brain neuroimaging and neurophysiological single-unit data from humans and non-human primates, the VTC-SOM predicted the topographical structure of fine-scale category-selective regions (face-, tool-, body-, and place-selective regions) and the boundary in large-scale abstract functional maps (animate vs. inanimate, real-word small-size vs. big-size, central vs. peripheral), with no significant loss in functionality (e.g., categorical selectivity and view-invariant representations). In addition, when the same principle was applied to V1 orientation preferences, a pinwheel-like topology emerged, suggesting the model's broad applicability. In summary, our study illustrates that the simple principle of wiring cost minimization, coupled with the appropriate biological constraint of lateral wiring span, is able to implement the high-dimensional object space in a two-dimensional cortical surface.

The virtual multiple sclerosis patient

Multiple sclerosis (MS) diagnosis typically involves assessing clinical symptoms, MRI findings, and ruling out alternative explanations. While myelin damage broadly affects conduction speeds, traditional tests focus on specific white-matter tracts, which may not reflect overall impairment accurately. In this study, we integrate diffusion tensor immaging (DTI) and magnetoencephalography (MEG) data into individualized virtual brain models to estimate conduction velocities for MS patients and controls. Using Bayesian inference, we demonstrated a causal link between empirical spectral changes and inferred slower conduction velocities in patients. Remarkably, these velocities proved superior predictors of clinical disability compared to structural damage. Our findings underscore a nuanced relationship between conduction delays and large-scale brain dynamics, suggesting that individualized velocity alterations at the whole-brain level contribute causatively to clinical outcomes in MS.

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.

Review of the Brain's Behaviour after Injury and Disease for Its Application in an Agent-Based Model (ABM)

The brain is the most complex organ in the human body and, as such, its study entails great challenges (methodological, theoretical, etc.). Nonetheless, there is a remarkable amount of studies about the consequences of pathological conditions on its development and functioning. This bibliographic review aims to cover mostly findings related to changes in the physical distribution of neurons and their connections-the connectome-both structural and functional, as well as their modelling approaches. It does not intend to offer an extensive description of all conditions affecting the brain; rather, it presents the most common ones. Thus, here, we highlight the need for accurate brain modelling that can subsequently be used to understand brain function and be applied to diagnose, track, and simulate treatments for the most prevalent pathologies affecting the brain.

Energy landscape analysis of brain network dynamics in Alzheimer's disease

Background

Alzheimer's disease (AD) is a common neurodegenerative dementia, characterized by abnormal dynamic functional connectivity (DFC). Traditional DFC analysis, assuming linear brain dynamics, may neglect the complexity of the brain's nonlinear interactions. Energy landscape analysis offers a holistic, nonlinear perspective to investigate brain network attractor dynamics, which was applied to resting-state fMRI data for AD in this study.

Methods

This study utilized resting-state fMRI data from 60 individuals, comparing 30 Alzheimer's patients with 30 controls, from the Alzheimer's Disease Neuroimaging Initiative. Energy landscape analysis was applied to the data to characterize the aberrant brain network dynamics of AD patients.

Results

The AD group stayed in the co-activation state for less time than the healthy control (HC) group, and a positive correlation was identified between the transition frequency of the co-activation state and behavior performance. Furthermore, the AD group showed a higher occurrence frequency and transition frequency of the cognitive control state and sensory integration state than the HC group. The transition between the two states was positively correlated with behavior performance.

Conclusion

The results suggest that the co-activation state could be important to cognitive processing and that the AD group possibly raised cognitive ability by increasing the occurrence and transition between the impaired cognitive control and sensory integration states.

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.

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.

Modeling brain network flexibility in networks of coupled oscillators: a feasibility study

Modeling the functionality of the human brain is a major goal in neuroscience for which many powerful methodologies have been developed over the last decade. The impact of working memory and the associated brain regions on the brain dynamics is of particular interest due to their connection with many functions and malfunctions in the brain. In this context, the concept of brain flexibility has been developed for the characterization of brain functionality. We discuss emergence of brain flexibility that is commonly measured by the identification of changes in the cluster structure of co-active brain regions. We provide evidence that brain flexibility can be modeled by a system of coupled FitzHugh-Nagumo oscillators where the network structure is obtained from human brain Diffusion Tensor Imaging (DTI). Additionally, we propose a straightforward and computationally efficient alternative macroscopic measure, which is derived from the Pearson distance of functional brain matrices. This metric exhibits similarities to the established patterns of brain template flexibility that have been observed in prior investigations. Furthermore, we explore the significance of the brain's network structure and the strength of connections between network nodes or brain regions associated with working memory in the observation of patterns in networks flexibility. This work enriches our understanding of the interplay between the structure and function of dynamic brain networks and proposes a modeling strategy to study brain flexibility.

The time-evolving epileptic brain network: concepts, definitions, accomplishments, perspectives

Epilepsy is now considered a network disease that affects the brain across multiple levels of spatial and temporal scales. The paradigm shift from an epileptic focus-a discrete cortical area from which seizures originate-to a widespread epileptic network-spanning lobes and hemispheres-considerably advanced our understanding of epilepsy and continues to influence both research and clinical treatment of this multi-faceted high-impact neurological disorder. The epileptic network, however, is not static but evolves in time which requires novel approaches for an in-depth characterization. In this review, we discuss conceptual basics of network theory and critically examine state-of-the-art recording techniques and analysis tools used to assess and characterize a time-evolving human epileptic brain network. We give an account on current shortcomings and highlight potential developments towards an improved clinical management of epilepsy.

Patient-specific modeling for guided rehabilitation of stroke patients: the BrainX3 use-case

BrainX3 is an interactive neuroinformatics platform that has been thoughtfully designed to support neuroscientists and clinicians with the visualization, analysis, and simulation of human neuroimaging, electrophysiological data, and brain models. The platform is intended to facilitate research and clinical use cases, with a focus on personalized medicine diagnostics, prognostics, and intervention decisions. BrainX3 is designed to provide an intuitive user experience and is equipped to handle different data types and 3D visualizations. To enhance patient-based analysis, and in keeping with the principles of personalized medicine, we propose a framework that can assist clinicians in identifying lesions and making patient-specific intervention decisions. To this end, we are developing an AI-based model for lesion identification, along with a mapping of tract information. By leveraging the patient's lesion information, we can gain valuable insights into the structural damage caused by the lesion. Furthermore, constraining whole-brain models with patient-specific disconnection masks can allow for the detection of mesoscale excitatory-inhibitory imbalances that cause disruptions in macroscale network properties. Finally, such information has the potential to guide neuromodulation approaches, assisting in the choice of candidate targets for stimulation techniques such as Transcranial Ultrasound Stimulation (TUS), which modulate E-I balance, potentiating cortical reorganization and the restoration of the dynamics and functionality disrupted due to the lesion.

A phenomenological model of whole brain dynamics using a network of neural oscillators with power-coupling

We present a general, trainable oscillatory neural network as a large-scale model of brain dynamics. The model has a cascade of two stages - an oscillatory stage and a complex-valued feedforward stage - for modelling the relationship between structural connectivity and functional connectivity from neuroimaging data under resting brain conditions. Earlier works of large-scale brain dynamics that used Hopf oscillators used linear coupling of oscillators. A distinctive feature of the proposed model employs a novel form of coupling known as power coupling. Oscillatory networks based on power coupling can accurately model arbitrary multi-dimensional signals. Training the lateral connections in the oscillator layer is done by a modified form of Hebbian learning, whereas a variation of the complex backpropagation algorithm does training in the second stage. The proposed model can not only model the empirical functional connectivity with remarkable accuracy (correlation coefficient between simulated and empirical functional connectivity- 0.99) but also identify default mode network regions. In addition, we also inspected how structural loss in the brain can cause significant aberration in simulated functional connectivity and functional connectivity dynamics; and how it can be restored with optimized model parameters by an in silico perturbational study.

Identifying good practices for detecting inter-regional linear functional connectivity from EEG

Aggregating voxel-level statistical dependencies between multivariate time series is an important intermediate step when characterising functional connectivity (FC) between larger brain regions. However, there are numerous ways in which voxel-level data can be aggregated into inter-regional FC, and the advantages of each of these approaches are currently unclear. In this study we generate ground-truth data and compare the performances of various pipelines that estimate directed and undirected linear phase-to-phase FC between regions. We test the ability of several existing and novel FC analysis pipelines to identify the true regions within which connectivity was simulated. We test various inverse modelling algorithms, strategies to aggregate time series within regions, and connectivity metrics. Furthermore, we investigate the influence of the number of interactions, the signal-to-noise ratio, the noise mix, the interaction time delay, and the number of active sources per region on the ability of detecting phase-to-phase FC. Throughout all simulated scenarios, lowest performance is obtained with pipelines involving the absolute value of coherency. Further, the combination of dynamic imaging of coherent sources (DICS) beamforming with directed FC metrics that aggregate information across multiple frequencies leads to unsatisfactory results. Pipelines that show promising results with our simulated pseudo-EEG data involve the following steps: (1) Source projection using the linearly-constrained minimum variance (LCMV) beamformer. (2) Principal component analysis (PCA) using the same fixed number of components within every region. (3) Calculation of the multivariate interaction measure (MIM) for every region pair to assess undirected phase-to-phase FC, or calculation of time-reversed Granger Causality (TRGC) to assess directed phase-to-phase FC. We formulate recommendations based on these results that may increase the validity of future experimental connectivity studies. We further introduce the free ROIconnect plugin for the EEGLAB toolbox that includes the recommended methods and pipelines that are presented here. We show an exemplary application of the best performing pipeline to the analysis of EEG data recorded during motor imagery.

Alzheimer's Disease: Insights from Large-Scale Brain Dynamics Models

Alzheimer's disease (AD) is a degenerative brain disease, and the condition is difficult to assess. In the past, numerous brain dynamics models have made remarkable contributions to neuroscience and the brain from the microcosmic to the macroscopic scale. Recently, large-scale brain dynamics models have been developed based on dual-driven multimodal neuroimaging data and neurodynamics theory. These models bridge the gap between anatomical structure and functional dynamics and have played an important role in assisting the understanding of the brain mechanism. Large-scale brain dynamics have been widely used to explain how macroscale neuroimaging biomarkers emerge from potential neuronal population level disturbances associated with AD. In this review, we describe this emerging approach to studying AD that utilizes a biophysically large-scale brain dynamics model. In particular, we focus on the application of the model to AD and discuss important directions for the future development and analysis of AD models. This will facilitate the development of virtual brain models in the field of AD diagnosis and treatment and add new opportunities for advancing clinical neuroscience.

Computational modelling in disorders of consciousness: Closing the gap towards personalised models for restoring consciousness

Disorders of consciousness are complex conditions characterised by persistent loss of responsiveness due to brain injury. They present diagnostic challenges and limited options for treatment, and highlight the urgent need for a more thorough understanding of how human consciousness arises from coordinated neural activity. The increasing availability of multimodal neuroimaging data has given rise to a wide range of clinically- and scientifically-motivated modelling efforts, seeking to improve data-driven stratification of patients, to identify causal mechanisms for patient pathophysiology and loss of consciousness more broadly, and to develop simulations as a means of testing in silico potential treatment avenues to restore consciousness. As a dedicated Working Group of clinicians and neuroscientists of the international Curing Coma Campaign, here we provide our framework and vision to understand the diverse statistical and generative computational modelling approaches that are being employed in this fast-growing field. We identify the gaps that exist between the current state-of-the-art in statistical and biophysical computational modelling in human neuroscience, and the aspirational goal of a mature field of modelling disorders of consciousness; which might drive improved treatments and outcomes in the clinic. Finally, we make several recommendations for how the field as a whole can work together to address these challenges.

A structural magnetic resonance imaging review of clinical motor outcomes from deep brain stimulation in movement disorders

Patients with movement disorders treated by deep brain stimulation do not always achieve successful therapeutic alleviation of motor symptoms, even in cases where surgery is without complications. Magnetic resonance imaging (MRI) offers methods to investigate structural brain-related factors that may be predictive of clinical motor outcomes. This review aimed to identify features which have been associated with variability in clinical post-operative motor outcomes in patients with Parkinson's disease, dystonia, and essential tremor from structural MRI modalities. We performed a literature search for articles published between 1 January 2000 and 1 April 2022 and identified 5197 articles. Following screening through our inclusion criteria, we identified 60 total studies (39 = Parkinson's disease, 11 = dystonia syndromes and 10 = essential tremor). The review captured a range of structural MRI methods and analysis techniques used to identify factors related to clinical post-operative motor outcomes from deep brain stimulation. Morphometric markers, including volume and cortical thickness were commonly identified in studies focused on patients with Parkinson's disease and dystonia syndromes. Reduced metrics in basal ganglia, sensorimotor and frontal regions showed frequent associations with reduced motor outcomes. Increased structural connectivity to subcortical nuclei, sensorimotor and frontal regions was also associated with greater motor outcomes. In patients with tremor, increased structural connectivity to the cerebellum and cortical motor regions showed high prevalence across studies for greater clinical motor outcomes. In addition, we highlight conceptual issues for studies assessing clinical response with structural MRI and discuss future approaches towards optimizing individualized therapeutic benefits. Although quantitative MRI markers are in their infancy for clinical purposes in movement disorder treatments, structural features obtained from MRI offer the powerful potential to identify candidates who are more likely to benefit from deep brain stimulation and provide insight into the complexity of disorder pathophysiology.

White-matter degradation and dynamical compensation support age-related functional alterations in human brain

Structural connectivity of the brain at different ages is analyzed using diffusion-weighted magnetic resonance imaging (MRI) data. The largest decrease of streamlines is found in frontal regions and for long inter-hemispheric links. The average length of the tracts also decreases, but the clustering is unaffected. From functional MRI we identify age-related changes of dynamic functional connectivity (dFC) and spatial covariation features of functional connectivity (FC) links captured by metaconnectivity. They indicate more stable dFC, but wider range and variance of MC, whereas static features of FC did not show any significant differences with age. We implement individual connectivity in whole-brain models and test several hypotheses for the mechanisms of operation among underlying neural system. We demonstrate that age-related functional fingerprints are only supported if the model accounts for: (i) compensation of the individual brains for the overall loss of structural connectivity and (ii) decrease of propagation velocity due to the loss of myelination. We also show that with these 2 conditions, it is sufficient to decompose the time-delays as bimodal distribution that only distinguishes between intra- and inter-hemispheric delays, and that the same working point also captures the static FC the best, and produces the largest variability at slow time-scales.