Inferring the dynamical effects of stroke lesions through whole-brain modeling

The localization of coma

Coma and disorders of consciousness (DoC) are common manifestations of acute severe brain injuries. Research into their neuroanatomical basis can be traced from Hippocrates to the present day. Lesions causing DoC have traditionally been conceptualized as decreasing "alertness" from damage to the ascending arousal system, and/or, reducing level of "awareness" due to structural or functional impairment of large-scale brain networks. Within this framework, pharmacological and neuromodulatory interventions to promote recovery from DoC have hitherto met with limited success. This is partly due to inter-individual heterogeneity of brain injury patterns, and an incomplete understanding of brain network properties that characterize consciousness. Advances in multiscale computational modelling of brain dynamics have opened a unique opportunity to explore the causal mechanisms of brain activity at the biophysical level. These models can provide a novel approach for selection and optimization of potential interventions by simulation of brain network dynamics individualized for each patient.

Whole-brain model replicates sleep-like slow-wave dynamics generated by stroke lesions

Focal brain injuries, such as stroke, cause local structural damage as well as alteration of neuronal activity in distant brain regions. Experimental evidence suggests that one of these changes is the appearance of sleep-like slow waves in the otherwise awake individual. This pattern is prominent in areas surrounding the damaged region and can extend to connected brain regions in a way consistent with the individual's specific long-range connectivity patterns. In this paper we present a generative whole-brain model based on (f)MRI data that, in combination with the disconnection mask associated with a given patient, explains the effects of the sleep-like slow waves originated in the vicinity of the lesion area on the distant brain activity. Our model reveals new aspects of their interaction, being able to reproduce functional connectivity patterns of stroke patients and offering a detailed, causal understanding of how stroke-related effects, in particular slow waves, spread throughout the brain. The presented findings demonstrate that the model effectively captures the links between stroke occurrences, sleep-like slow waves, and their subsequent spread across the human brain.

Metric comparison of connectome-based lesion-symptom mapping in post-stroke aphasia

Connectome-based lesion-symptom mapping relates behavioural impairments to disruption of structural brain connectivity. Connectome-based lesion-symptom mapping can be based on different approaches (diffusion MRI versus lesion mask), network scales (whole brain versus regions of interest) and measure types (tract-based, parcel-based, or network-based metrics). We evaluated the similarity of different connectome-based lesion-symptom mapping processing choices and identified factors that influence the results using multiverse analysis-the strategy of conducting and displaying the results of all reasonable processing choices. Metrics derived from lesion masks and diffusion-weighted images were tested for association with Boston Naming Test and Token Test performance in a sample of 50 participants with aphasia following left hemispheric stroke. 'Direct' measures were derived from diffusion-weighted images. 'Indirect' measures were derived by overlaying lesion masks on a white matter atlas. Parcel-based connectomes were constructed for the whole brain and regions of interest (14 language-relevant parcels). Numerous tract-based and network-based metrics were calculated. There was a high discrepancy across processing approaches (diffusion-weighted images versus lesion masks), network scales (whole brain versus regions of interest) and metric types. Results indicate weak correlations and different connectome-based lesion-symptom mapping results across the processing choices. Substantial methodological work is needed to validate the various decision points that arise when conducting connectome-based lesion-symptom mapping analyses. Multiverse analysis is a useful strategy for evaluating the similarity across different processing choices in connectome-based lesion-symptom mapping.

Generative whole-brain dynamics models from healthy subjects predict functional alterations in stroke at the level of individual patients

Computational whole-brain models describe the resting activity of each brain region based on a local model, inter-regional functional interactions, and a structural connectome that specifies the strength of inter-regional connections. Strokes damage the healthy structural connectome that forms the backbone of these models and produce large alterations in inter-regional functional interactions. These interactions are typically measured by correlating the time series of the activity between two brain regions in a process, called resting functional connectivity. We show that adding information about the structural disconnections produced by a patient's lesion to a whole-brain model previously trained on structural and functional data from a large cohort of healthy subjects enables the prediction of the resting functional connectivity of the patient and fits the model directly to the patient's data (Pearson correlation = 0.37; mean square error = 0.005). Furthermore, the model dynamics reproduce functional connectivity-based measures that are typically abnormal in stroke patients and measures that specifically isolate these abnormalities. Therefore, although whole-brain models typically involve a large number of free parameters, the results show that, even after fixing those parameters, the model reproduces results from a population very different than that on which the model was trained. In addition to validating the model, these results show that the model mechanistically captures the relationships between the anatomical structure and the functional activity of the human brain.

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.

The Hopf whole-brain model and its linear approximation

Whole-brain models have proven to be useful to understand the emergence of collective activity among neural populations or brain regions. These models combine connectivity matrices, or connectomes, with local node dynamics, noise, and, eventually, transmission delays. Multiple choices for the local dynamics have been proposed. Among them, nonlinear oscillators corresponding to a supercritical Hopf bifurcation have been used to link brain connectivity and collective phase and amplitude dynamics in different brain states. Here, we studied the linear fluctuations of this model to estimate its stationary statistics, i.e., the instantaneous and lagged covariances and the power spectral densities. This linear approximation-that holds in the case of heterogeneous parameters and time-delays-allows analytical estimation of the statistics and it can be used for fast parameter explorations to study changes in brain state, changes in brain activity due to alterations in structural connectivity, and modulations of parameter due to non-equilibrium dynamics.

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 low dimensional embedding of brain dynamics enhances diagnostic accuracy and behavioral prediction in stroke

Large-scale brain networks reveal structural connections as well as functional synchronization between distinct regions of the brain. The latter, referred to as functional connectivity (FC), can be derived from neuroimaging techniques such as functional magnetic resonance imaging (fMRI). FC studies have shown that brain networks are severely disrupted by stroke. However, since FC data are usually large and high-dimensional, extracting clinically useful information from this vast amount of data is still a great challenge, and our understanding of the functional consequences of stroke remains limited. Here, we propose a dimensionality reduction approach to simplify the analysis of this complex neural data. By using autoencoders, we find a low-dimensional representation encoding the fMRI data which preserves the typical FC anomalies known to be present in stroke patients. By employing the latent representations emerging from the autoencoders, we enhanced patients' diagnostics and severity classification. Furthermore, we showed how low-dimensional representation increased the accuracy of recovery prediction.

Multiscale effects of excitatory-inhibitory homeostasis in lesioned cortical networks: A computational study

Stroke-related disruptions in functional connectivity (FC) often spread beyond lesioned areas and, given the localized nature of lesions, it is unclear how the recovery of FC is orchestrated on a global scale. Since recovery is accompanied by long-term changes in excitability, we propose excitatory-inhibitory (E-I) homeostasis as a driving mechanism. We present a large-scale model of the neocortex, with synaptic scaling of local inhibition, showing how E-I homeostasis can drive the post-lesion restoration of FC and linking it to changes in excitability. We show that functional networks could reorganize to recover disrupted modularity and small-worldness, but not network dynamics, suggesting the need to consider forms of plasticity beyond synaptic scaling of inhibition. On average, we observed widespread increases in excitability, with the emergence of complex lesion-dependent patterns related to biomarkers of relevant side effects of stroke, such as epilepsy, depression and chronic pain. In summary, our results show that the effects of E-I homeostasis extend beyond local E-I balance, driving the restoration of global properties of FC, and relating to post-stroke symptomatology. Therefore, we suggest the framework of E-I homeostasis as a relevant theoretical foundation for the study of stroke recovery and for understanding the emergence of meaningful features of FC from local dynamics.