Whole-Brain Multimodal Neuroimaging Model Using Serotonin Receptor Maps Explains Non-linear Functional Effects of LSD

A computational approach to evaluate how molecular mechanisms impact large-scale brain activity

Assessing the impact of pharmaceutical compounds on brain activity is a critical issue in contemporary neuroscience. Currently, no systematic approach exists for evaluating these effects in whole-brain models, which typically focus on macroscopic phenomena, while pharmaceutical interventions operate at the molecular scale. Here we address this issue by presenting a computational approach for brain simulations using biophysically grounded mean-field models that integrate membrane conductances and synaptic receptors, showcased in the example of anesthesia. We show that anesthetics targeting GABAA and NMDA receptors can switch brain activity to generalized slow-wave patterns, as observed experimentally in deep anesthesia. To validate our models, we demonstrate that these slow-wave states exhibit reduced responsiveness to external stimuli and functional connectivity constrained by anatomical connectivity, mirroring experimental findings in anesthetized states across species. Our approach, founded on mean-field models that incorporate molecular realism, provides a robust framework for understanding how molecular-level drug actions impact whole-brain dynamics.

Network control energy reductions under DMT relate to serotonin receptors, signal diversity, and subjective experience

Psychedelics offer a profound window into the human brain through their robust effects on perception, subjective experience, and brain activity patterns. The serotonergic psychedelic N,N-dimethyltryptamine (DMT) induces a profoundly immersive altered state of consciousness lasting under 20 min, allowing the entire experience to be captured during a single functional magnetic resonance imaging (fMRI) scan. Using network control theory, we map energy trajectories of 14 individuals undergoing fMRI during DMT and placebo. We find that global control energy is reduced after DMT injection compared to placebo. Longitudinal trajectories of global control energy correlate with longitudinal trajectories of electroencephalography (EEG) signal diversity (a measure of entropy) and subjective drug intensity ratings. At the regional level, spatial patterns of DMT's effects on these metrics correlate with serotonin 2a receptor density from positron emission tomography (PET) data. Using receptor distribution and pharmacokinetic information, we recapitulate DMT's effects on global control energy trajectories, demonstrating control models can predict pharmacological effects on brain dynamics.

Transient destabilization of whole brain dynamics induced by N,N-Dimethyltryptamine (DMT)

The transition towards the brain state induced by psychedelic drugs is frequently neglected in favor of a static description of their acute effects. We use a time-dependent whole-brain model to reproduce large-scale brain dynamics measured with fMRI from 15 volunteers under 20 mg intravenous N,N-Dimethyltryptamine (DMT), a short-acting psychedelic. To capture its transient effects, we parametrize the proximity to a global bifurcation using a pharmacokinetic equation. Simulated perturbations reveal a transient of heightened reactivity concentrated in fronto-parietal regions and visual cortices, correlated with serotonin 5HT2a receptor density, the primary target of psychedelics. These advances suggest a mechanism to explain key features of the psychedelic state and also predicts that the temporal evolution of these features aligns with pharmacokinetics. Our results contribute to understanding how psychedelics introduce a transient where minimal perturbations can achieve a maximal effect, shedding light on how short psychedelic episodes may extend an overarching influence over time.

Modulatory Neurotransmitter Genotypes Shape Dynamic Functional Connectome Reconfigurations

Dynamic reconfigurations of the functional connectome across different connectivity states are highly heritable, predictive of cognitive abilities, and linked to mental health. Despite their established heritability, the specific polymorphisms that shape connectome dynamics are largely unknown. Given the widespread regulatory impact of modulatory neurotransmitters on functional connectivity, we comprehensively investigated a large set of single nucleotide polymorphisms (SNPs) of their receptors, metabolic enzymes, and transporters in 674 healthy adult subjects (347 females) from the Human Connectome Project. Preregistered modulatory neurotransmitter SNPs and dynamic connectome features entered a Stability Selection procedure with resampling. We found that specific subsets of these SNPs explain individual differences in temporal phenotypes of fMRI-derived connectome dynamics for which we previously established heritability. Specifically, noradrenergic polymorphisms explained Fractional Occupancy, i.e., the proportion of time spent in each connectome state, and cholinergic polymorphisms explained Transition Probability, i.e., the probability to transition between state pairs, respectively. This work identifies specific genetic effects on connectome dynamics via the regulatory impact of modulatory neurotransmitter systems. Our observations highlight the potential of dynamic connectome features as endophenotypes for neurotransmitter-focused precision psychiatry.

Neural mass modeling for the masses: Democratizing access to whole-brain biophysical modeling with FastDMF

Different whole-brain computational models have been recently developed to investigate hypotheses related to brain mechanisms. Among these, the Dynamic Mean Field (DMF) model is particularly attractive, combining a biophysically realistic model that is scaled up via a mean-field approach and multimodal imaging data. However, an important barrier to the widespread usage of the DMF model is that current implementations are computationally expensive, supporting only simulations on brain parcellations that consider less than 100 brain regions. Here, we introduce an efficient and accessible implementation of the DMF model: the FastDMF. By leveraging analytical and numerical advances-including a novel estimation of the feedback inhibition control parameter and a Bayesian optimization algorithm-the FastDMF circumvents various computational bottlenecks of previous implementations, improving interpretability, performance, and memory use. Furthermore, these advances allow the FastDMF to increase the number of simulated regions by one order of magnitude, as confirmed by the good fit to fMRI data parcellated at 90 and 1,000 regions. These advances open the way to the widespread use of biophysically grounded whole-brain models for investigating the interplay between anatomy, function, and brain dynamics and to identify mechanistic explanations of recent results obtained from fine-grained neuroimaging recordings.

Synergistic, multi-level understanding of psychedelics: three systematic reviews and meta-analyses of their pharmacology, neuroimaging and phenomenology

Serotonergic psychedelics induce altered states of consciousness and have shown potential for treating a variety of neuropsychiatric disorders, including depression and addiction. Yet their modes of action are not fully understood. Here, we provide a novel, synergistic understanding of psychedelics arising from systematic reviews and meta-analyses of three hierarchical levels of analysis: (1) subjective experience (phenomenology), (2) neuroimaging and (3) molecular pharmacology. Phenomenologically, medium and high doses of LSD yield significantly higher ratings of visionary restructuralisation than psilocybin on the 5-dimensional Altered States of Consciousness Scale. Our neuroimaging results reveal that, in general, psychedelics significantly strengthen between-network functional connectivity (FC) while significantly diminishing within-network FC. Pharmacologically, LSD induces significantly more inositol phosphate formation at the 5-HT2A receptor than DMT and psilocin, yet there are no significant between-drug differences in the selectivity of psychedelics for the 5-HT2A, 5-HT2C, or D2 receptors, relative to the 5-HT1A receptor. Our meta-analyses link DMT, LSD, and psilocybin to specific neural fingerprints at each level of analysis. The results show a highly non-linear relationship between these fingerprints. Overall, our analysis highlighted the high heterogeneity and risk of bias in the literature. This suggests an urgent need for standardising experimental procedures and analysis techniques, as well as for more research on the emergence between different levels of psychedelic effects.

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.

Psychedelic 5-HT2A receptor agonism: neuronal signatures and altered neurovascular coupling

Psychedelics hold therapeutic promise for mood disorders due to rapid, sustained results. Human neuroimaging studies have reported dramatic serotonin-2A receptor-(5-HT2AR)-dependent changes in functional brain reorganization that presumably reflect neuromodulation. However, the potent vasoactive effects of serotonin have been overlooked. We found psilocybin-mediated alterations to fMRI-HRFs in humans, suggesting potentially altered NVC. To assess the neuronal, hemodynamic, and neurovascular coupling (NVC) effects of the psychedelic 5-HT2AR agonist, 2,5-Dimethoxy-4-iodoamphetamine (DOI), wide-field optical imaging (WFOI) was used in awake Thy1-jRGECO1a mice during stimulus-evoked and resting-state conditions. While DOI partially altered tasked-based NVC, more pronounced NVC alterations occurred under resting-state conditions and were strongest in association regions. Further, calcium and hemodynamic activity reported different accounts of RSFC changes under DOI. Co-administration of DOI and the 5-HT2AR antagonist, MDL100907, reversed many of these effects. Dissociation between neuronal and hemodynamic signals emphasizes a need to consider neurovascular effects of psychedelics when interpreting blood-oxygenation-dependent neuroimaging measures.

The Algorithmic Agent Perspective and Computational Neuropsychiatry: From Etiology to Advanced Therapy in Major Depressive Disorder

Major Depressive Disorder (MDD) is a complex, heterogeneous condition affecting millions worldwide. Computational neuropsychiatry offers potential breakthroughs through the mechanistic modeling of this disorder. Using the Kolmogorov theory (KT) of consciousness, we developed a foundational model where algorithmic agents interact with the world to maximize an Objective Function evaluating affective valence. Depression, defined in this context by a state of persistently low valence, may arise from various factors-including inaccurate world models (cognitive biases), a dysfunctional Objective Function (anhedonia, anxiety), deficient planning (executive deficits), or unfavorable environments. Integrating algorithmic, dynamical systems, and neurobiological concepts, we map the agent model to brain circuits and functional networks, framing potential etiological routes and linking with depression biotypes. Finally, we explore how brain stimulation, psychotherapy, and plasticity-enhancing compounds such as psychedelics can synergistically repair neural circuits and optimize therapies using personalized computational models.

Revisiting the role of computational neuroimaging in the era of integrative neuroscience

Computational models have become integral to human neuroimaging research, providing both mechanistic insights and predictive tools for human cognition and behavior. However, concerns persist regarding the ecological validity of lab-based neuroimaging studies and whether their spatiotemporal resolution is not sufficient for capturing neural dynamics. This review aims to re-examine the utility of computational neuroimaging, particularly in light of the growing prominence of alternative neuroscientific methods and the growing emphasis on more naturalistic behaviors and paradigms. Specifically, we will explore how computational modeling can both enhance the analysis of high-dimensional imaging datasets and, conversely, how neuroimaging, in conjunction with other data modalities, can inform computational models through the lens of neurobiological plausibility. Collectively, this evidence suggests that neuroimaging remains critical for human neuroscience research, and when enhanced by computational models, imaging can serve an important role in bridging levels of analysis and understanding. We conclude by proposing key directions for future research, emphasizing the development of standardized paradigms and the integrative use of computational modeling across neuroimaging techniques.

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.

Neuroimaging Biomarkers for Drug Discovery and Development in Schizophrenia

Schizophrenia is a chronic mental illness that affects up to 1% of the population. While efficacious therapies are available for positive symptoms, effective treatment of cognitive and negative symptoms remains an unmet need after decades of research. New developments in the field of neuroimaging are accelerating our knowledge gain regarding the underlying pathophysiology of symptoms in schizophrenia and psychosis spectrum disorders, inspiring new targets for drug development. However, no validated and qualified biomarkers are currently available to support the development of new therapeutics. This review summarizes the current use of neuroimaging technology in clinical drug development for psychotic disorders. As exemplified by drug development programs that target NMDA receptor hypofunction, neuroimaging results play a critical role in target discovery and establishing target engagement and dose selection. Furthermore, pharmacological neuroimaging may provide response biomarkers that allow for early decision making in proof-of-concept studies that leverage pharmacological challenge models in healthy volunteers. That said, while response and predictive biomarkers are starting to be evaluated in patient populations, they continue to play a limited role. Novel approaches to neuroimaging data acquisition and analysis may aid the establishment of biomarkers that are predictive at the individual level in the future. Nevertheless, various gaps in knowledge need to be addressed and biomarkers need to be validated to establish them as "fit for purpose" in drug development.

EEG microstate transition cost correlates with task demands

The ability to solve complex tasks relies on the adaptive changes occurring in the spatio-temporal organization of brain activity under different conditions. Altered flexibility in these dynamics can lead to impaired cognitive performance, manifesting for instance as difficulties in attention regulation, distraction inhibition, and behavioral adaptation. Such impairments result in decreased efficiency and increased effort in accomplishing goal-directed tasks. Therefore, developing quantitative measures that can directly assess the effort involved in these transitions using neural data is of paramount importance. In this study, we propose a framework to associate cognitive effort during the performance of tasks with electroencephalography (EEG) activation patterns. The methodology relies on the identification of discrete dynamical states (EEG microstates) and optimal transport theory. To validate the effectiveness of this framework, we apply it to a dataset collected during a spatial version of the Stroop task, a cognitive test in which participants respond to one aspect of a stimulus while ignoring another, often conflicting, aspect. The Stroop task is a cognitive test where participants must respond to one aspect of a stimulus while ignoring another, often conflicting, aspect. Our findings reveal an increased cost linked to cognitive effort, thus confirming the framework's effectiveness in capturing and quantifying cognitive transitions. By utilizing a fully data-driven method, this research opens up fresh perspectives for physiologically describing cognitive effort within the brain.

Beyond the usual suspects: multi-factorial computational models in the search for neurodegenerative disease mechanisms

From Alzheimer's disease to amyotrophic lateral sclerosis, the molecular cascades underlying neurodegenerative disorders remain poorly understood. The clinical view of neurodegeneration is confounded by symptomatic heterogeneity and mixed pathology in almost every patient. While the underlying physiological alterations originate, proliferate, and propagate potentially decades before symptomatic onset, the complexity and inaccessibility of the living brain limit direct observation over a patient's lifespan. Consequently, there is a critical need for robust computational methods to support the search for causal mechanisms of neurodegeneration by distinguishing pathogenic processes from consequential alterations, and inter-individual variability from intra-individual progression. Recently, promising advances have been made by data-driven spatiotemporal modeling of the brain, based on in vivo neuroimaging and biospecimen markers. These methods include disease progression models comparing the temporal evolution of various biomarkers, causal models linking interacting biological processes, network propagation models reproducing the spatial spreading of pathology, and biophysical models spanning cellular- to network-scale phenomena. In this review, we discuss various computational approaches for integrating cross-sectional, longitudinal, and multi-modal data, primarily from large observational neuroimaging studies, to understand (i) the temporal ordering of physiological alterations, i(i) their spatial relationships to the brain's molecular and cellular architecture, (iii) mechanistic interactions between biological processes, and (iv) the macroscopic effects of microscopic factors. We consider the extents to which computational models can evaluate mechanistic hypotheses, explore applications such as improving treatment selection, and discuss how model-informed insights can lay the groundwork for a pathobiological redefinition of neurodegenerative disorders.

Whole brain modelling for simulating pharmacological interventions on patients with disorders of consciousness

Disorders of consciousness (DoC) represent a challenging and complex group of neurological conditions characterised by profound disturbances in consciousness. The current range of treatments for DoC is limited. This has sparked growing interest in developing new treatments, including the use of psychedelic drugs. Nevertheless, clinical investigations and the mechanisms behind them are methodologically and ethically constrained. To tackle these limitations, we combined biologically plausible whole-brain models with deep learning techniques to characterise the low-dimensional space of DoC patients. We investigated the effects of model pharmacological interventions by including the whole-brain dynamical consequences of the enhanced neuromodulatory level of different neurotransmitters, and providing geometrical interpretation in the low-dimensional space. Our findings show that serotonergic and opioid receptors effectively shifted the DoC models towards a dynamical behaviour associated with a healthier state, and that these improvements correlated with the mean density of the activated receptors throughout the brain. These findings mark an important step towards the development of treatments not only for DoC but also for a broader spectrum of brain diseases. Our method offers a promising avenue for exploring the therapeutic potential of pharmacological interventions within the ethical and methodological confines of clinical research.

Contributions of network structure, chemoarchitecture and diagnostic categories to transitions between cognitive topographies

The mechanisms linking the brain's network structure to cognitively relevant activation patterns remain largely unknown. Here, by leveraging principles of network control, we show how the architecture of the human connectome shapes transitions between 123 experimentally defined cognitive activation maps (cognitive topographies) from the NeuroSynth meta-analytic database. Specifically, we systematically integrated large-scale multimodal neuroimaging data from functional magnetic resonance imaging, diffusion tractography, cortical morphometry and positron emission tomography to simulate how anatomically guided transitions between cognitive states can be reshaped by neurotransmitter engagement or by changes in cortical thickness. Our model incorporates neurotransmitter-receptor density maps (18 receptors and transporters) and maps of cortical thickness pertaining to a wide range of mental health, neurodegenerative, psychiatric and neurodevelopmental diagnostic categories (17,000 patients and 22,000 controls). The results provide a comprehensive look-up table charting how brain network organization and chemoarchitecture interact to manifest different cognitive topographies, and establish a principled foundation for the systematic identification of ways to promote selective transitions between cognitive topographies.

Breathwork-induced psychedelic experiences modulate neural dynamics

Breathwork is an understudied school of practices involving intentional respiratory modulation to induce an altered state of consciousness (ASC). We simultaneously investigate the phenomenological and neural dynamics of breathwork by combining Temporal Experience Tracing, a quantitative methodology that preserves the temporal dynamics of subjective experience, with low-density portable EEG devices. Fourteen novice participants completed a course of up to 28 breathwork sessions-of 20, 40, or 60 min-in 28 days, yielding a neurophenomenological dataset of 301 breathwork sessions. Using hypothesis-driven and data-driven approaches, we found that "psychedelic-like" subjective experiences were associated with increased neural Lempel-Ziv complexity during breathwork. Exploratory analyses showed that the aperiodic exponent of the power spectral density-but not oscillatory alpha power-yielded similar neurophenomenological associations. Non-linear neural features, like complexity and the aperiodic exponent, neurally map both a multidimensional data-driven composite of positive experiences, and hypothesis-driven aspects of psychedelic-like experience states such as high bliss.

Unravelling consciousness and brain function through the lens of time, space, and information

Disentangling how cognitive functions emerge from the interplay of brain dynamics and network architecture is among the major challenges that neuroscientists face. Pharmacological and pathological perturbations of consciousness provide a lens to investigate these complex challenges. Here, we review how recent advances about consciousness and the brain's functional organisation have been driven by a common denominator: decomposing brain function into fundamental constituents of time, space, and information. Whereas unconsciousness increases structure-function coupling across scales, psychedelics may decouple brain function from structure. Convergent effects also emerge: anaesthetics, psychedelics, and disorders of consciousness can exhibit similar reconfigurations of the brain's unimodal-transmodal functional axis. Decomposition approaches reveal the potential to translate discoveries across species, with computational modelling providing a path towards mechanistic integration.

In-vivo neuronal dysfunction by Aβ and tau overlaps with brain-wide inflammatory mechanisms in Alzheimer's disease

The molecular mechanisms underlying neuronal dysfunction in Alzheimer's disease (AD) remain uncharacterized. Here, we identify genes, molecular pathways and cellular components associated with whole-brain dysregulation caused by amyloid-beta (Aβ) and tau deposits in the living human brain. We obtained in-vivo resting-state functional MRI (rs-fMRI), Aβ- and tau-PET for 47 cognitively unimpaired and 16 AD participants from the Translational Biomarkers in Aging and Dementia cohort. Adverse neuronal activity impacts by Aβ and tau were quantified with personalized dynamical models by fitting pathology-mediated computational signals to the participant's real rs-fMRIs. Then, we detected robust brain-wide associations between the spatial profiles of Aβ-tau impacts and gene expression in the neurotypical transcriptome (Allen Human Brain Atlas). Within the obtained distinctive signature of in-vivo neuronal dysfunction, several genes have prominent roles in microglial activation and in interactions with Aβ and tau. Moreover, cellular vulnerability estimations revealed strong association of microglial expression patterns with Aβ and tau's synergistic impact on neuronal activity (q < 0.001). These results further support the central role of the immune system and neuroinflammatory pathways in AD pathogenesis. Neuronal dysregulation by AD pathologies also associated with neurotypical synaptic and developmental processes. In addition, we identified drug candidates from the vast LINCS library to halt or reduce the observed Aβ-tau effects on neuronal activity. Top-ranked pharmacological interventions target inflammatory, cancer and cardiovascular pathways, including specific medications undergoing clinical evaluation in AD. Our findings, based on the examination of molecular-pathological-functional interactions in humans, may accelerate the process of bringing effective therapies into clinical practice.

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.

Psychedelics and disorders of consciousness: the current landscape and the path forward

Modern medicine has been shaken by the surge of psychedelic science that proposes a new approach to mitigate mental disorders, such as depression and post-traumatic stress disorder. Clinical trials to investigate whether psychedelic substances can treat psychiatric conditions are now underway, yet less discussion gravitates around their use in neurological disorders due to brain injury. One suggested implementation of brain-complexity enhancing psychedelics is to treat people with post-comatose disorders of consciousness (DoC). In this article, we discuss the rationale of this endeavour, examining possible outcomes of such experiments by postulating the existence of an optimal level of complexity. We consider the possible counterintuitive effects of both psychedelics and DoC on the functional connectivity of the default mode network and its possible impact on selfhood. We also elaborate on the role of computational modelling in providing complementary information to experimental studies, both contributing to our understanding of the treatment mechanisms and providing a path towards personalized medicine. Finally, we update the discourse surrounding the ethical considerations, encompassing clinical and scientific values.

The Thermodynamics of Mind

To not only survive, but also thrive, the brain must efficiently orchestrate distributed computation across space and time. This requires hierarchical organisation facilitating fast information transfer and processing at the lowest possible metabolic cost. Quantifying brain hierarchy is difficult but can be estimated from the asymmetry of information flow. Thermodynamics has successfully characterised hierarchy in many other complex systems. Here, we propose the 'Thermodynamics of Mind' framework as a natural way to quantify hierarchical brain orchestration and its underlying mechanisms. This has already provided novel insights into the orchestration of hierarchy in brain states including movie watching, where the hierarchy of the brain is flatter than during rest. Overall, this framework holds great promise for revealing the orchestration of cognition.

Gaming expertise induces meso‑scale brain plasticity and efficiency mechanisms as revealed by whole-brain modeling

Video games are a valuable tool for studying the effects of training and neural plasticity on the brain. However, the underlying mechanisms related to plasticity-associated brain structural changes and their impact on brain dynamics are unknown. Here, we used a semi-empirical whole-brain model to study structural neural plasticity mechanisms linked to video game expertise. We hypothesized that video game expertise is associated with neural plasticity-mediated changes in structural connectivity that manifest at the meso‑scale level, resulting in a more segregated functional network topology. To test this hypothesis, we combined structural connectivity data of StarCraft II video game players (VGPs, n = 31) and non-players (NVGPs, n = 31), with generic fMRI data from the Human Connectome Project and computational models, to generate simulated fMRI recordings. Graph theory analysis on simulated data was performed during both resting-state conditions and external stimulation. VGPs' simulated functional connectivity was characterized by a meso‑scale integration, with increased local connectivity in frontal, parietal, and occipital brain regions. The same analyses at the level of structural connectivity showed no differences between VGPs and NVGPs. Regions that increased their connectivity strength in VGPs are known to be involved in cognitive processes crucial for task performance such as attention, reasoning, and inference. In-silico stimulation suggested that differences in FC between VGPs and NVGPs emerge in noisy contexts, specifically when the noisy level of stimulation is increased. This indicates that the connectomes of VGPs may facilitate the filtering of noise from stimuli. These structural alterations drive the meso‑scale functional changes observed in individuals with gaming expertise. Overall, our work sheds light on the mechanisms underlying structural neural plasticity triggered by video game experiences.

Personalized whole-brain neural mass models reveal combined Aβ and tau hyperexcitable influences in Alzheimer's disease

Neuronal dysfunction and cognitive deterioration in Alzheimer's disease (AD) are likely caused by multiple pathophysiological factors. However, mechanistic evidence in humans remains scarce, requiring improved non-invasive techniques and integrative models. We introduce personalized AD computational models built on whole-brain Wilson-Cowan oscillators and incorporating resting-state functional MRI, amyloid-β (Aβ) and tau-PET from 132 individuals in the AD spectrum to evaluate the direct impact of toxic protein deposition on neuronal activity. This subject-specific approach uncovers key patho-mechanistic interactions, including synergistic Aβ and tau effects on cognitive impairment and neuronal excitability increases with disease progression. The data-derived neuronal excitability values strongly predict clinically relevant AD plasma biomarker concentrations (p-tau217, p-tau231, p-tau181, GFAP) and grey matter atrophy obtained through voxel-based morphometry. Furthermore, reconstructed EEG proxy quantities show the hallmark AD electrophysiological alterations (theta band activity enhancement and alpha reductions) which occur with Aβ-positivity and after limbic tau involvement. Microglial activation influences on neuronal activity are less definitive, potentially due to neuroimaging limitations in mapping neuroprotective vs detrimental activation phenotypes. Mechanistic brain activity models can further clarify intricate neurodegenerative processes and accelerate preventive/treatment interventions.

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.

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.

The unique neural signature of your trip: Functional connectome fingerprints of subjective psilocybin experience

The emerging neuroscientific frontier of brain fingerprinting has recently established that human functional connectomes (FCs) exhibit fingerprint-like idiosyncratic features, which map onto heterogeneously distributed behavioral traits. Here, we harness brain-fingerprinting tools to extract FC features that predict subjective drug experience induced by the psychedelic psilocybin. Specifically, in neuroimaging data of healthy volunteers under the acute influence of psilocybin or a placebo, we show that, post psilocybin administration, FCs become more idiosyncratic owing to greater intersubject dissimilarity. Moreover, whereas in placebo subjects idiosyncratic features are primarily found in the frontoparietal network, in psilocybin subjects they concentrate in the default mode network (DMN). Crucially, isolating the latter revealed an FC pattern that predicts subjective psilocybin experience and is characterized by reduced within-DMN and DMN-limbic connectivity, as well as increased connectivity between the DMN and attentional systems. Overall, these results contribute to bridging the gap between psilocybin-mediated effects on brain and behavior, while demonstrating the value of a brain-fingerprinting approach to pharmacological neuroimaging.

Brain dynamics predictive of response to psilocybin for treatment-resistant depression

Psilocybin therapy for depression has started to show promise, yet the underlying causal mechanisms are not currently known. Here, we leveraged the differential outcome in responders and non-responders to psilocybin (10 and 25 mg, 7 days apart) therapy for depression-to gain new insights into regions and networks implicated in the restoration of healthy brain dynamics. We used large-scale brain modelling to fit the spatiotemporal brain dynamics at rest in both responders and non-responders before treatment. Dynamic sensitivity analysis of systematic perturbation of these models enabled us to identify specific brain regions implicated in a transition from a depressive brain state to a healthy one. Binarizing the sample into treatment responders (>50% reduction in depressive symptoms) versus non-responders enabled us to identify a subset of regions implicated in this change. Interestingly, these regions correlate with in vivo density maps of serotonin receptors 5-hydroxytryptamine 2a and 5-hydroxytryptamine 1a, which psilocin, the active metabolite of psilocybin, has an appreciable affinity for, and where it acts as a full-to-partial agonist. Serotonergic transmission has long been associated with depression, and our findings provide causal mechanistic evidence for the role of brain regions in the recovery from depression via psilocybin.

Neural Geometrodynamics, Complexity, and Plasticity: A Psychedelics Perspective

We explore the intersection of neural dynamics and the effects of psychedelics in light of distinct timescales in a framework integrating concepts from dynamics, complexity, and plasticity. We call this framework neural geometrodynamics for its parallels with general relativity's description of the interplay of spacetime and matter. The geometry of trajectories within the dynamical landscape of "fast time" dynamics are shaped by the structure of a differential equation and its connectivity parameters, which themselves evolve over "slow time" driven by state-dependent and state-independent plasticity mechanisms. Finally, the adjustment of plasticity processes (metaplasticity) takes place in an "ultraslow" time scale. Psychedelics flatten the neural landscape, leading to heightened entropy and complexity of neural dynamics, as observed in neuroimaging and modeling studies linking increases in complexity with a disruption of functional integration. We highlight the relationship between criticality, the complexity of fast neural dynamics, and synaptic plasticity. Pathological, rigid, or "canalized" neural dynamics result in an ultrastable confined repertoire, allowing slower plastic changes to consolidate them further. However, under the influence of psychedelics, the destabilizing emergence of complex dynamics leads to a more fluid and adaptable neural state in a process that is amplified by the plasticity-enhancing effects of psychedelics. This shift manifests as an acute systemic increase of disorder and a possibly longer-lasting increase in complexity affecting both short-term dynamics and long-term plastic processes. Our framework offers a holistic perspective on the acute effects of these substances and their potential long-term impacts on neural structure and function.

A role for the serotonin 2A receptor in the expansion and functioning of human transmodal cortex

Integrating independent but converging lines of research on brain function and neurodevelopment across scales, this article proposes that serotonin 2A receptor (5-HT2AR) signalling is an evolutionary and developmental driver and potent modulator of the macroscale functional organization of the human cerebral cortex. A wealth of evidence indicates that the anatomical and functional organization of the cortex follows a unimodal-to-transmodal gradient. Situated at the apex of this processing hierarchy-where it plays a central role in the integrative processes underpinning complex, human-defining cognition-the transmodal cortex has disproportionately expanded across human development and evolution. Notably, the adult human transmodal cortex is especially rich in 5-HT2AR expression and recent evidence suggests that, during early brain development, 5-HT2AR signalling on neural progenitor cells stimulates their proliferation-a critical process for evolutionarily-relevant cortical expansion. Drawing on multimodal neuroimaging and cross-species investigations, we argue that, by contributing to the expansion of the human cortex and being prevalent at the apex of its hierarchy in the adult brain, 5-HT2AR signalling plays a major role in both human cortical expansion and functioning. Owing to its unique excitatory and downstream cellular effects, neuronal 5-HT2AR agonism promotes neuroplasticity, learning and cognitive and psychological flexibility in a context-(hyper)sensitive manner with therapeutic potential. Overall, we delineate a dual role of 5-HT2ARs in enabling both the expansion and modulation of the human transmodal cortex.

Psychedelia: The interplay of music and psychedelics

Music and psychedelics have been intertwined throughout the existence of Homo sapiens, from the early shamanic rituals of the Americas and Africa to the modern use of psychedelic-assisted therapy for a variety of mental health conditions. Across such settings, music has been highly prized for its ability to guide the psychedelic experience. Here, we examine the interplay between music and psychedelics, starting by describing their association with the brain's functional hierarchy that is relied upon for music perception and its psychedelic-induced manipulation, as well as an exploration of the limited research on their mechanistic neural overlap. We explore music's role in Western psychedelic therapy and the use of music in indigenous psychedelic rituals, with a specific focus on ayahuasca and the Santo Daime Church. Furthermore, we explore work relating to the evolution and onset of music and psychedelic use. Finally, we consider music's potential to lead to altered states of consciousness in the absence of psychedelics as well as the development of psychedelic music. Here, we provide an overview of several perspectives on the interaction between psychedelic use and music-a topic with growing interest given increasing excitement relating to the therapeutic efficacy of psychedelic interventions.

Towards a biologically annotated brain connectome

The brain is a network of interleaved neural circuits. In modern connectomics, brain connectivity is typically encoded as a network of nodes and edges, abstracting away the rich biological detail of local neuronal populations. Yet biological annotations for network nodes - such as gene expression, cytoarchitecture, neurotransmitter receptors or intrinsic dynamics - can be readily measured and overlaid on network models. Here we review how connectomes can be represented and analysed as annotated networks. Annotated connectomes allow us to reconceptualize architectural features of networks and to relate the connection patterns of brain regions to their underlying biology. Emerging work demonstrates that annotated connectomes help to make more veridical models of brain network formation, neural dynamics and disease propagation. Finally, annotations can be used to infer entirely new inter-regional relationships and to construct new types of network that complement existing connectome representations. In summary, biologically annotated connectomes offer a compelling way to study neural wiring in concert with local biological features.

Gaming expertise induces meso-scale brain plasticity and efficiency mechanisms as revealed by whole-brain modeling Gaming expertise, neuroplasticity and functional dynamics

Video games are a valuable tool for studying the effects of training and neural plasticity on the brain. However, the underlaying mechanisms related to plasticity-induced brain structural changes and their impact in brain dynamics are unknown. Here, we used a semi-empirical whole-brain model to study structural neural plasticity mechanisms linked to video game expertise. We hypothesized that video game expertise is associated with neural plasticity-mediated changes in structural connectivity that manifest at the meso-scale level, resulting in a more segregated functional network topology. To test this hypothesis, we combined structural connectivity data of StarCraft II video game players (VGPs, n = 31) and non-players (NVGPs, n = 31), with generic fMRI data from the Human Connectome Project and computational models, with the aim of generating simulated fMRI recordings. Graph theory analysis on simulated data was performed during both resting-state conditions and external stimulation. VGPs' simulated functional connectivity was characterized by a meso-scale integration, with increased local connectivity in frontal, parietal and occipital brain regions. The same analyses at the level of structural connectivity showed no differences between VGPs and NVGPs. Regions that increased their connectivity strength in VGPs are known to be involved in cognitive processes crucial for task performance such as attention, reasoning, and inference. In-silico stimulation suggested that differences in FC between VGPs and NVGPs emerge in noisy contexts, specifically when the noisy level of stimulation is increased. This indicates that the connectomes of VGPs may facilitate the filtering of noise from stimuli. These structural alterations drive the meso-scale functional changes observed in individuals with gaming expertise. Overall, our work sheds light into the mechanisms underlying structural neural plasticity triggered by video game experiences.

Structural basis of envelope and phase intrinsic coupling modes in the cerebral cortex

Intrinsic coupling modes (ICMs) can be observed in ongoing brain activity at multiple spatial and temporal scales. Two families of ICMs can be distinguished: phase and envelope ICMs. The principles that shape these ICMs remain partly elusive, in particular their relation to the underlying brain structure. Here we explored structure-function relationships in the ferret brain between ICMs quantified from ongoing brain activity recorded with chronically implanted micro-ECoG arrays and structural connectivity (SC) obtained from high-resolution diffusion MRI tractography. Large-scale computational models were used to explore the ability to predict both types of ICMs. Importantly, all investigations were conducted with ICM measures that are sensitive or insensitive to volume conduction effects. The results show that both types of ICMs are significantly related to SC, except for phase ICMs when using measures removing zero-lag coupling. The correlation between SC and ICMs increases with increasing frequency which is accompanied by reduced delays. Computational models produced results that were highly dependent on the specific parameter settings. The most consistent predictions were derived from measures solely based on SC. Overall, the results demonstrate that patterns of cortical functional coupling as reflected in both phase and envelope ICMs are both related, albeit to different degrees, to the underlying structural connectivity in the cerebral cortex.

Modeling the role of the thalamus in resting-state functional connectivity: Nature or structure

The thalamus is a central brain structure that serves as a relay station for sensory inputs from the periphery to the cortex and regulates cortical arousal. Traditionally, it has been regarded as a passive relay that transmits information between brain regions. However, recent studies have suggested that the thalamus may also play a role in shaping functional connectivity (FC) in a task-based context. Based on this idea, we hypothesized that due to its centrality in the network and its involvement in cortical activation, the thalamus may also contribute to resting-state FC, a key neurological biomarker widely used to characterize brain function in health and disease. To investigate this hypothesis, we constructed ten in-silico brain network models based on neuroimaging data (MEG, MRI, and dwMRI), and simulated them including and excluding the thalamus, and raising the noise into thalamus to represent the afferences related to the reticular activating system (RAS) and the relay of peripheral sensory inputs. We simulated brain activity and compared the resulting FC to their empirical MEG counterparts to evaluate model's performance. Results showed that a parceled version of the thalamus with higher noise, able to drive damped cortical oscillators, enhanced the match to empirical FC. However, with an already active self-oscillatory cortex, no impact on the dynamics was observed when introducing the thalamus. We also demonstrated that the enhanced performance was not related to the structural connectivity of the thalamus, but to its higher noisy inputs. Additionally, we highlighted the relevance of a balanced signal-to-noise ratio in thalamus to allow it to propagate its own dynamics. In conclusion, our study sheds light on the role of the thalamus in shaping brain dynamics and FC in resting-state and allowed us to discuss the general role of criticality in the brain at the mesoscale level.

Revealing the combined roles of Aβ and tau in Alzheimer's disease via a pathophysiological activity decoder

Neuronal dysfunction and cognitive deterioration in Alzheimer's disease (AD) are likely caused by multiple pathophysiological factors. However, evidence in humans remains scarce, necessitating improved non-invasive techniques and integrative mechanistic models. Here, we introduce personalized brain activity models incorporating functional MRI, amyloid-β (Aβ) and tau-PET from AD-related participants ( N = 132 ) . Within the model assumptions, electrophysiological activity is mediated by toxic protein deposition. Our integrative subject-specific approach uncovers key patho-mechanistic interactions, including synergistic Aβ and tau effects on cognitive impairment and neuronal excitability increases with disease progression. The data-derived neuronal excitability values strongly predict clinically relevant AD plasma biomarker concentrations (p-tau217, p-tau231, p-tau181, GFAP). Furthermore, our results reproduce hallmark AD electrophysiological alterations (theta band activity enhancement and alpha reductions) which occur with Aβ-positivity and after limbic tau involvement. Microglial activation influences on neuronal activity are less definitive, potentially due to neuroimaging limitations in mapping neuroprotective vs detrimental phenotypes. Mechanistic brain activity models can further clarify intricate neurodegenerative processes and accelerate preventive/treatment interventions.

Molecular-enriched functional connectivity in the human brain using multiband multi-echo simultaneous ASL/BOLD fMRI

Receptor-enriched analysis of functional connectivity by targets (REACT) is a strategy to enrich functional MRI (fMRI) data with molecular information on the neurotransmitter distribution density in the human brain, providing a biological basis to the functional connectivity (FC) analysis. Although this approach has been used in BOLD fMRI studies only so far, extending its use to ASL imaging would provide many advantages, including the more direct link of ASL with neuronal activity compared to BOLD and its suitability for pharmacological MRI studies assessing drug effects on baseline brain function. Here, we applied REACT to simultaneous ASL/BOLD resting-state fMRI data of 29 healthy subjects and estimated the ASL and BOLD FC maps related to six molecular systems. We then compared the ASL and BOLD FC maps in terms of spatial similarity, and evaluated and compared the test-retest reproducibility of each modality. We found robust spatial patterns of molecular-enriched FC for both modalities, moderate similarity between BOLD and ASL FC maps and comparable reproducibility for all but one molecular-enriched functional networks. Our findings showed that ASL is as informative as BOLD in detecting functional circuits associated with specific molecular pathways, and that the two modalities may provide complementary information related to these circuits.

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.

Gradients of neurotransmitter receptor expression in the macaque cortex

Dynamics and functions of neural circuits depend on interactions mediated by receptors. Therefore, a comprehensive map of receptor organization across cortical regions is needed. In this study, we used in vitro receptor autoradiography to measure the density of 14 neurotransmitter receptor types in 109 areas of macaque cortex. We integrated the receptor data with anatomical, genetic and functional connectivity data into a common cortical space. We uncovered a principal gradient of receptor expression per neuron. This aligns with the cortical hierarchy from sensory cortex to higher cognitive areas. A second gradient, driven by serotonin 5-HT1A receptors, peaks in the anterior cingulate, default mode and salience networks. We found a similar pattern of 5-HT1A expression in the human brain. Thus, the macaque may be a promising translational model of serotonergic processing and disorders. The receptor gradients may enable rapid, reliable information processing in sensory cortical areas and slow, flexible integration in higher cognitive areas.

From neurotransmitters to networks: Transcending organisational hierarchies with molecular-informed functional imaging

The human brain exhibits complex interactions across micro, meso-, and macro-scale organisational principles. Recent synergistic multi-modal approaches have begun to link micro-scale information to systems level dynamics, transcending organisational hierarchies and offering novel perspectives into the brain's function and dysfunction. Specifically, the distribution of micro-scale properties (such as receptor density or gene expression) can be mapped onto macro-scale measures from functional MRI to provide novel neurobiological insights. Methodological approaches to enrich functional imaging analyses with molecular information are rapidly evolving, with several streams of research having developed relatively independently, each offering unique potential to explore the trans-hierarchical functioning of the brain. Here, we address the three principal streams of research - spatial correlation, molecular-enriched network, and in-silico whole brain modelling analyses - to provide a critical overview of the different sources of molecular information, how this information can be utilised within analyses of fMRI data, the merits and pitfalls of each methodology, and, through the use of key examples, highlight their promise to shed new light on key domains of neuroscientific inquiry.

Spatial Correspondence of LSD-Induced Variations on Brain Functioning at Rest With Serotonin Receptor Expression

BACKGROUND

Lysergic acid diethylamide (LSD) is an atypical psychedelic compound that exerts its effects through pleiotropic actions, mainly involving 1A/2A serotoninergic (5-HT) receptor subtypes. However, the mechanisms by which LSD promotes a reorganization of the brain's functional activity and connectivity are still partially unknown.

METHODS

Our study analyzed resting-state functional magnetic resonance imaging data acquired from 15 healthy volunteers undergoing LSD single-dose intake. A voxelwise analysis investigated the alterations of the brain's intrinsic functional connectivity and local signal amplitude induced by LSD or by a placebo. Quantitative comparisons assessed the spatial overlap between these 2 indices of functional reorganization and the topography of receptor expression obtained from a publicly available collection of in vivo, whole-brain atlases. Finally, linear regression models explored the relationships between changes in resting-state functional magnetic resonance imaging and behavioral aspects of the psychedelic experience.

RESULTS

LSD elicited modifications of the cortical functional architecture that spatially overlapped with the distribution of serotoninergic receptors. Local signal amplitude and functional connectivity increased in regions belonging to the default mode and attention networks associated with high expression of 5-HT2A receptors. These functional changes correlate with the occurrence of simple and complex visual hallucinations. At the same time, a decrease in local signal amplitude and intrinsic connectivity was observed in limbic areas, which are dense with 5-HT1A receptors.

CONCLUSIONS

This study provides new insights into the neural processes underlying the brain network reconfiguration induced by LSD. It also identifies a topographical relationship between opposite effects on brain functioning and the spatial distribution of different 5-HT receptors.

The impact of regional heterogeneity in whole-brain dynamics in the presence of oscillations

Large variability exists across brain regions in health and disease, considering their cellular and molecular composition, connectivity, and function. Large-scale whole-brain models comprising coupled brain regions provide insights into the underlying dynamics that shape complex patterns of spontaneous brain activity. In particular, biophysically grounded mean-field whole-brain models in the asynchronous regime were used to demonstrate the dynamical consequences of including regional variability. Nevertheless, the role of heterogeneities when brain dynamics are supported by synchronous oscillating state, which is a ubiquitous phenomenon in brain, remains poorly understood. Here, we implemented two models capable of presenting oscillatory behavior with different levels of abstraction: a phenomenological Stuart-Landau model and an exact mean-field model. The fit of these models informed by structural- to functional-weighted MRI signal (T1w/T2w) allowed us to explore the implication of the inclusion of heterogeneities for modeling resting-state fMRI recordings from healthy participants. We found that disease-specific regional functional heterogeneity imposed dynamical consequences within the oscillatory regime in fMRI recordings from neurodegeneration with specific impacts on brain atrophy/structure (Alzheimer's patients). Overall, we found that models with oscillations perform better when structural and functional regional heterogeneities are considered, showing that phenomenological and biophysical models behave similarly at the brink of the Hopf bifurcation.

In vivo mapping of pharmacologically induced functional reorganization onto the human brain's neurotransmitter landscape

To understand how pharmacological interventions can exert their powerful effects on brain function, we need to understand how they engage the brain's rich neurotransmitter landscape. Here, we bridge microscale molecular chemoarchitecture and pharmacologically induced macroscale functional reorganization, by relating the regional distribution of 19 neurotransmitter receptors and transporters obtained from positron emission tomography, and the regional changes in functional magnetic resonance imaging connectivity induced by 10 different mind-altering drugs: propofol, sevoflurane, ketamine, lysergic acid diethylamide (LSD), psilocybin, N,N-Dimethyltryptamine (DMT), ayahuasca, 3,4-methylenedioxymethamphetamine (MDMA), modafinil, and methylphenidate. Our results reveal a many-to-many mapping between psychoactive drugs' effects on brain function and multiple neurotransmitter systems. The effects of both anesthetics and psychedelics on brain function are organized along hierarchical gradients of brain structure and function. Last, we show that regional co-susceptibility to pharmacological interventions recapitulates co-susceptibility to disorder-induced structural alterations. Collectively, these results highlight rich statistical patterns relating molecular chemoarchitecture and drug-induced reorganization of the brain's functional architecture.

Neuromodulatory control of complex adaptive dynamics in the brain

How is the massive dimensionality and complexity of the microscopic constituents of the nervous system brought under sufficiently tight control so as to coordinate adaptive behaviour? A powerful means for striking this balance is to poise neurons close to the critical point of a phase transition, at which a small change in neuronal excitability can manifest a nonlinear augmentation in neuronal activity. How the brain could mediate this critical transition is a key open question in neuroscience. Here, I propose that the different arms of the ascending arousal system provide the brain with a diverse set of heterogeneous control parameters that can be used to modulate the excitability and receptivity of target neurons-in other words, to act as control parameters for mediating critical neuronal order. Through a series of worked examples, I demonstrate how the neuromodulatory arousal system can interact with the inherent topological complexity of neuronal subsystems in the brain to mediate complex adaptive behaviour.

The current and future contribution of neuroimaging to the understanding of disorders of consciousness

Patients with disorders of consciousness (DoC) represent a group of severely brain-injured patients with varying capacities for consciousness in terms of both wakefulness and awareness. The current state-of-the-art for assessing these patients is through standardised behavioural examinations, but inaccuracies are commonplace. Neuroimaging and electrophysiological techniques have revealed vast insights into the relationships between neural alterations, andcognitive and behavioural features of consciousness in patients with DoC. This has led to the establishment of neuroimaging paradigms for the clinical assessment of DoC patients. Here, we review selected neuroimaging findings on the DoC population, outlining key findings of the dysfunction underlying DoC and presenting the current clinical utility of neuroimaging tools. We discuss that whilst individual brain areas play instrumental roles in generating and supporting consciousness, activation of these areas alone is not sufficient for conscious experience. Instead, for consciousness to arise, we need preserved thalamo-cortical circuits, in addition to sufficient connectivity between distinctly differentiated brain networks, underlined by connectivity both within, and between such brain networks. Finally, we present recent advances and future perspectives in computational methodologies applied to DoC, supporting the notion that progress in the science of DoC will be driven by a symbiosis of these data-driven analyses, and theory-driven research. Both perspectives will work in tandem to provide mechanistic insights contextualised within theoretical frameworks which ultimately inform the practice of clinical neurology.

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.

Time-resolved network control analysis links reduced control energy under DMT with the serotonin 2a receptor, signal diversity, and subjective experience

Psychedelics offer a profound window into the functioning of the human brain and mind through their robust acute effects on perception, subjective experience, and brain activity patterns. In recent work using a receptor-informed network control theory framework, we demonstrated that the serotonergic psychedelics lysergic acid diethylamide (LSD) and psilocybin flatten the brain's control energy landscape in a manner that covaries with more dynamic and entropic brain activity. Contrary to LSD and psilocybin, whose effects last for hours, the serotonergic psychedelic N,N-dimethyltryptamine (DMT) rapidly induces a profoundly immersive altered state of consciousness lasting less than 20 minutes, allowing for the entirety of the drug experience to be captured during a single resting-state fMRI scan. Using network control theory, which quantifies the amount of input necessary to drive transitions between functional brain states, we integrate brain structure and function to map the energy trajectories of 14 individuals undergoing fMRI during DMT and placebo. Consistent with previous work, we find that global control energy is reduced following injection with DMT compared to placebo. We additionally show longitudinal trajectories of global control energy correlate with longitudinal trajectories of EEG signal diversity (a measure of entropy) and subjective ratings of drug intensity. We interrogate these same relationships on a regional level and find that the spatial patterns of DMT's effects on these metrics are correlated with serotonin 2a receptor density (obtained from separately acquired PET data). Using receptor distribution and pharmacokinetic information, we were able to successfully recapitulate the effects of DMT on global control energy trajectories, demonstrating a proof-of-concept for the use of control models in predicting pharmacological intervention effects on brain dynamics.

A whole-brain model of the neural entropy increase elicited by psychedelic drugs

Psychedelic drugs, including lysergic acid diethylamide (LSD) and other agonists of the serotonin 2A receptor (5HT2A-R), induce drastic changes in subjective experience, and provide a unique opportunity to study the neurobiological basis of consciousness. One of the most notable neurophysiological signatures of psychedelics, increased entropy in spontaneous neural activity, is thought to be of relevance to the psychedelic experience, mediating both acute alterations in consciousness and long-term effects. However, no clear mechanistic explanation for this entropy increase has been put forward so far. We sought to do this here by building upon a recent whole-brain model of serotonergic neuromodulation, to study the entropic effects of 5HT2A-R activation. Our results reproduce the overall entropy increase observed in previous experiments in vivo, providing the first model-based explanation for this phenomenon. We also found that entropy changes were not uniform across the brain: entropy increased in all regions, but the larger effect were localised in visuo-occipital regions. Interestingly, at the whole-brain level, this reconfiguration was not well explained by 5HT2A-R density, but related closely to the topological properties of the brain's anatomical connectivity. These results help us understand the mechanisms underlying the psychedelic state and, more generally, the pharmacological modulation of whole-brain activity.

Reduced emergent character of neural dynamics in patients with a disrupted connectome

High-level brain functions are widely believed to emerge from the orchestrated activity of multiple neural systems. However, lacking a formal definition and practical quantification of emergence for experimental data, neuroscientists have been unable to empirically test this long-standing conjecture. Here we investigate this fundamental question by leveraging a recently proposed framework known as "Integrated Information Decomposition," which establishes a principled information-theoretic approach to operationalise and quantify emergence in dynamical systems - including the human brain. By analysing functional MRI data, our results show that the emergent and hierarchical character of neural dynamics is significantly diminished in chronically unresponsive patients suffering from severe brain injury. At a functional level, we demonstrate that emergence capacity is positively correlated with the extent of hierarchical organisation in brain activity. Furthermore, by combining computational approaches from network control theory and whole-brain biophysical modelling, we show that the reduced capacity for emergent and hierarchical dynamics in severely brain-injured patients can be mechanistically explained by disruptions in the patients' structural connectome. Overall, our results suggest that chronic unresponsiveness resulting from severe brain injury may be related to structural impairment of the fundamental neural infrastructures required for brain dynamics to support emergence.