Reliability and subject specificity of personalized whole-brain dynamical models

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.

Viscous dynamics associated with hypoexcitation and structural disintegration in neurodegeneration via generative whole-brain modeling

INTRODUCTION

Alzheimer's disease (AD) and behavioral variant frontotemporal dementia (bvFTD) lack mechanistic biophysical modeling in diverse, underrepresented populations. Electroencephalography (EEG) is a high temporal resolution, cost-effective technique for studying dementia globally, but lacks mechanistic models and produces non-replicable results.

METHODS

We developed a generative whole-brain model that combines EEG source-level metaconnectivity, anatomical priors, and a perturbational approach. This model was applied to Global South participants (AD, bvFTD, and healthy controls).

RESULTS

Metaconnectivity outperformed pairwise connectivity and revealed more viscous dynamics in patients, with altered metaconnectivity patterns associated with multimodal disease presentation. The biophysical model showed that connectome disintegration and hypoexcitability triggered altered metaconnectivity dynamics and identified critical regions for brain stimulation. We replicated the main results in a second subset of participants for validation with unharmonized, heterogeneous recording settings.

DISCUSSION

The results provide a novel agenda for developing mechanistic model-inspired characterization and therapies in clinical, translational, and computational neuroscience settings.

The Commonality and Individuality of Human Brains When Performing Tasks

It is imperative to study individual brain functioning toward understanding the neural bases responsible for individual behavioral and clinical traits. The complex and dynamic brain activity varies from area to area and from time to time across the entire brain, and BOLD-fMRI measures this spatiotemporal activity at large-scale systems level. We present a novel method to investigate task-evoked whole brain activity that varies not only from person to person but also from task trial to trial within each task type, offering a means of characterizing the individuality of human brains when performing tasks. For each task trial, the temporal correlation of task-evoked ideal time signal with the time signal of every point in the brain yields a full spatial map that characterizes the whole brain's functional co-activity (FC) relative to the task-evoked ideal response. For any two task trials, regardless of whether they are the same task or not, the spatial correlation of their corresponding two FC maps over the entire brain quantifies the similarity between these two maps, offering a means of investigating the variation in the whole brain activity trial to trial. The results demonstrated a substantially varied whole brain activity from trial to trial for each task category. The degree of this variation was task type-dependent and varied from subject to subject, showing a remarkable individuality of human brains when performing tasks. It demonstrates the potential of using the presented method to investigate the relationship of the whole brain activity with individual behavioral and clinical traits.

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.

Moving from phenomenological to predictive modelling: Progress and pitfalls of modelling brain stimulation in-silico

Brain stimulation is an increasingly popular neuromodulatory tool used in both clinical and research settings; however, the effects of brain stimulation, particularly those of non-invasive stimulation, are variable. This variability can be partially explained by an incomplete mechanistic understanding, coupled with a combinatorial explosion of possible stimulation parameters. Computational models constitute a useful tool to explore the vast sea of stimulation parameters and characterise their effects on brain activity. Yet the utility of modelling stimulation in-silico relies on its biophysical relevance, which needs to account for the dynamics of large and diverse neural populations and how underlying networks shape those collective dynamics. The large number of parameters to consider when constructing a model is no less than those needed to consider when planning empirical studies. This piece is centred on the application of phenomenological and biophysical models in non-invasive brain stimulation. We first introduce common forms of brain stimulation and computational models, and provide typical construction choices made when building phenomenological and biophysical models. Through the lens of four case studies, we provide an account of the questions these models can address, commonalities, and limitations across studies. We conclude by proposing future directions to fully realise the potential of computational models of brain stimulation for the design of personalized, efficient, and effective stimulation strategies.

Whole-brain modeling explains the context-dependent effects of cholinergic neuromodulation

Integration and segregation are two fundamental principles of brain organization. The brain manages the transitions and balance between different functional segregated or integrated states through neuromodulatory systems. Recently, computational and experimental studies suggest a pro-segregation effect of cholinergic neuromodulation. Here, we studied the effects of the cholinergic system on brain functional connectivity using both empirical fMRI data and computational modeling. First, we analyzed the effects of nicotine on functional connectivity and network topology in healthy subjects during resting-state conditions and during an attentional task. Then, we employed a whole-brain neural mass model interconnected using a human connectome to simulate the effects of nicotine and investigate causal mechanisms for these changes. The drug effect was modeled decreasing both the global coupling and local feedback inhibition parameters, consistent with the known cellular effects of acetylcholine. We found that nicotine incremented functional segregation in both empirical and simulated data, and the effects are context-dependent: observed during the task, but not in the resting state. In-task performance correlates with functional segregation, establishing a link between functional network topology and behavior. Furthermore, we found in the empirical data that the regional density of the nicotinic acetylcholine α₄β₂ correlates with the decrease in functional nodal strength by nicotine during the task. Our results confirm that cholinergic neuromodulation promotes functional segregation in a context-dependent fashion, and suggest that this segregation is suited for simple visual-attentional tasks.

A Holistic Analysis of Individual Brain Activity Revealed the Relationship of Brain Areal Activity with the Entire Brain's Activity

The relationship between brain areal activity and the entire brain's activity is unknown, and understanding this relationship is imperative for understanding the neural mechanisms of human brain function at systems level. The complex activity of human brains varies from area to area and from time to time across the whole brain. BOLD-fMRI measures this spatiotemporal activity at a large-scale systems level. The BOLD time signal of an area reflects a collective neuronal activity of over one million neurons under that area, and the temporal correlation of this time signal with that of every point in the brain yields a full spatial map that characterizes the entire brain's functional co-activity (FC) relative to that area's activity. Here we show a quantitative relationship between brain areal activity and the activity of the entire brain. The temporal correlation coefficient r of the signal time courses of two areas quantifies the degree of co-activity between the two areas, and the spatial correlation coefficient R of their corresponding two FC maps quantifies the co-activity between these two maps. We found that a modified sigmoid function quantified this R with r, i.e., Rr=1+ra-1-ra1+ra+1-ra, revealing a relationship between the activity of brain areas and that of the entire brain. The parameter a in this equation was found to be associated with the mean degree of the temporal co-activity among all brain areas, and its value was brain functional state dependent too. Our study demonstrated a novel approach for analyzing fMRI data to holistically characterize the entire brain's activity quantitatively for any brain functional state in individual humans.

Whole-brain dynamical modelling for classification of Parkinson's disease

Simulated whole-brain connectomes demonstrate enhanced inter-individual variability depending on the data processing and modelling approach. By considering the human brain connectome as an individualized attribute, we investigate how empirical and simulated whole-brain connectome-derived features can be utilized to classify patients with Parkinson's disease against healthy controls in light of varying data processing and model validation. To this end, we applied simulated blood oxygenation level-dependent signals derived by a whole-brain dynamical model simulating electrical signals of neuronal populations to reveal differences between patients and controls. In addition to the widely used model validation via fitting the dynamical model to empirical neuroimaging data, we invented a model validation against behavioural data, such as subject classes, which we refer to as behavioural model fitting and show that it can be beneficial for Parkinsonian patient classification. Furthermore, the results of machine learning reported in this study also demonstrated that the performance of the patient classification can be improved when the empirical data are complemented by the simulation results. We also showed that the temporal filtering of blood oxygenation level-dependent signals influences the prediction results, where filtering in the low-frequency band is advisable for Parkinsonian patient classification. In addition, composing the feature space of empirical and simulated data from multiple brain parcellation schemes provided complementary features that improved prediction performance. Based on our findings, we suggest that combining the simulation results with empirical data is effective for inter-individual research and its clinical application.