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 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 1000 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.

Multiple sclerosis (MS) is a neurodegenerative disease that affects the central nervous system. Structures affected in MS include the corpus callosum, connecting the hemispheres. Studies have shown that in mammalian brains, structural connectivity is organized according to a conservation principle, an inverse relationship between intra- and interhemispheric connectivity. The aim of this study was to replicate this conservation principle in subjects with MS and to explore how the disease interacts with it. A multicentric dataset has been analyzed including 513 people with MS and 208 healthy controls from 7 different centres. Structural connectivity was quantified through various connectivity measures and graph analysis was used to study the behavior of intra- and interhemispheric connectivity. The association between the intra- and the interhemispheric connectivity showed a similar strength for healthy controls (r = 0.38; p < 0.001) and people with MS (r = 0.35; p < 0.001). Intrahemispheric connectivity was associated with white matter fraction (r = 0.48, p < 0.0001), lesion volume (r = −0.44, p < 0.0001), and the Symbol Digit Modalities Test (SDMT) (r = 0.25, p < 0.0001). Results show that this conservation principle seems to hold for people with MS. These findings support the hypothesis that interhemispheric connectivity decreases at higher cognitive decline and disability levels, while intrahemispheric connectivity increases to maintain the balance. Author Summary In our study, we investigated how multiple sclerosis (MS), a disease that affects the central nervous system, impacts brain connectivity and the conservation principle of connectivity in the brain across hemispheres. By analyzing data from 513 MS patients and 208 healthy individuals, we examined if this conservation principle holds and how it changes due to MS. Our findings revealed that both MS patients and healthy individuals exhibit a similar balance between connections within each hemisphere and between hemispheres. We also observed that as cognitive impairment and disability in MS patients increase, interhemispheric connectivity decreases while intrahemispheric connectivity compensates. This suggests that the brain attempts to maintain balance despite the disease’s progression, highlighting the adaptability of neural connectivity during the course of MS.

Contemplative neuroscience has increasingly explored meditation using neuroimaging. However, the brain mechanisms underlying meditation remain elusive. Here, we implemented a mechanistic framework to explore the spatiotemporal dynamics of expert meditators during meditation and rest, and controls during rest. We first applied a model-free approach by defining a probabilistic metastable substate (PMS) space for each condition, consisting of different probabilities of occurrence from a repertoire of dynamic patterns. Moreover, we implemented a model-based approach by adjusting the PMS of each condition to a whole-brain model, which enabled us to explore in silico perturbations to transition from resting-state to meditation and vice versa. Consequently, we assessed the sensitivity of different brain areas regarding their perturbability and their mechanistic local-global effects. Overall, our work reveals distinct whole-brain dynamics in meditation compared to rest, and how transitions can be induced with localized artificial perturbations. It motivates future work regarding meditation as a practice in health and as a potential therapy for brain disorders. Author Summary Our work explores brain dynamics in a group of expert meditators and controls. First, we characterized meditation and rest with a repertoire of brain patterns, each with its distinct probability of occurrence. Then, we generated whole-brain models of each condition, which enabled us to artificially perturb the systems to induce transitions between rest and meditation. Our results open new avenues in meditation research as a practice in health and disease.

It has been previously shown that traumatic brain injury (TBI) is associated with reductions in metastability in large-scale networks in resting-state fMRI (rsfMRI). However, little is known about how TBI affects the local level of synchronization and how this evolves during the recovery trajectory. Here, we applied a novel turbulent dynamics framework to investigate whole-brain dynamics using an rsfMRI dataset from a cohort of moderate to severe TBI patients and healthy controls (HCs). We first examined how several measures related to turbulent dynamics differ between HCs and TBI patients at 3, 6, and 12 months post-injury. We found a significant reduction in these empirical measures after TBI, with the largest change at 6 months post-injury. Next, we built a Hopf whole-brain model with coupled oscillators and conducted in silico perturbations to investigate the mechanistic principles underlying the reduced turbulent dynamics found in the empirical data. A simulated attack was used to account for the effect of focal lesions. This revealed a shift to lower coupling parameters in the TBI dataset and, critically, decreased susceptibility and information-encoding capability. These findings confirm the potential of the turbulent framework to characterize longitudinal changes in whole-brain dynamics and in the reactivity to external perturbations after TBI. Author Summary Significant advances in nonlinear dynamics and computational modeling have opened up the possibility of studying how whole-brain dynamics may be impacted by traumatic brain injury (TBI). One hypothesis is that whole-brain dynamics may show differential spatiotemporal patterns during the recovery trajectory. We used a novel turbulence framework to examine how several measures related to turbulent dynamics differ between healthy controls and TBI patients at 3, 6, and 12 months post-injury. This framework revealed a significant reduction in these empirical measures after TBI differentially affecting long distances, with the largest change at 6 months post-injury. The alterations observed at network level, however, showed a certain degree of recovery after 1 year. In addition, the Hopf whole-brain model demonstrated decreased susceptibility and information-encoding capability after TBI. The clinical implications of this work are discussed.

Author Summary Here we used a thermodynamics-inspired, deep learning based Temporal Evolution NETwork (TENET) framework to assess the asymmetry in the flow of events, ‘arrow of time’, in human brain signals. This was applied to large-scale HCP neuroimaging data which showed significant changes between the hierarchy of orchestration for the resting state and seven different cognitive tasks. Similarly, the hierarchy of the resting state is significantly different in health compared to neuropsychiatric disorders. This framework provided new insights into the orchestrating of brain dynamics in different brain states. Abstract A promising idea in human cognitive neuroscience is that the default mode network (DMN) is responsible for coordinating the recruitment and scheduling of networks for computing and solving task-specific cognitive problems. This is supported by evidence showing that the physical and functional distance of DMN regions is maximally removed from sensorimotor regions containing environment-driven neural activity directly linked to perception and action, which would allow the DMN to orchestrate complex cognition from the top of the hierarchy. However, discovering the functional hierarchy of brain dynamics requires finding the best way to measure interactions between brain regions. In contrast to previous methods measuring the hierarchical flow of information using, for example, transfer entropy, here we used a thermodynamics-inspired, deep learning based Temporal Evolution NETwork (TENET) framework to assess the asymmetry in the flow of events, ‘arrow of time’, in human brain signals. This provides an alternative way of quantifying hierarchy, given that the arrow of time measures the directionality of information flow that leads to a breaking of the balance of the underlying hierarchy. In turn, the arrow of time is a measure of nonreversibility and thus nonequilibrium in brain dynamics. When applied to large-scale Human Connectome Project (HCP) neuroimaging data from close to a thousand participants, the TENET framework suggests that the DMN plays a significant role in orchestrating the hierarchy, that is, levels of nonreversibility, which changes between the resting state and when performing seven different cognitive tasks. Furthermore, this quantification of the hierarchy of the resting state is significantly different in health compared to neuropsychiatric disorders. Overall, the present thermodynamics-based machine-learning framework provides vital new insights into the fundamental tenets of brain dynamics for orchestrating the interactions between cognition and brain in complex environments.

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. Author Summary Significant progress has been made in understanding the effects of regional heterogeneity on whole-brain dynamics. With imaging technologies, the number of high-resolution reference maps of brain structure and function has been increased, and whole-brain computational models have provided a suitable avenue to investigate the mechanisms supporting the relations between these maps and whole-brain dynamics. Here, we investigate the role of the heterogeneities when synchronous behavior is present in brain dynamics, which we could represent by models capable of oscillating in the presence of a Hopf bifurcation. We found that models with oscillations more faithfully reproduce empirical properties when structural and functional regional heterogeneities are considered, showing that both phenomenological and biophysical models behave similarly at the brink of the Hopf bifurcation.

Beyond the established effects of subthalamic nucleus deep brain stimulation (STN-DBS) in reducing motor symptoms in Parkinson’s disease, recent evidence has highlighted the effect on non-motor symptoms. However, the impact of STN-DBS on disseminated networks remains unclear. This study aimed to perform a quantitative evaluation of network-specific modulation induced by STN-DBS using Leading Eigenvector Dynamics Analysis (LEiDA). We calculated the occupancy of resting-state networks (RSNs) in functional MRI data from 10 patients with Parkinson’s disease implanted with STN-DBS and statistically compared between ON and OFF conditions. STN-DBS was found to specifically modulate the occupancy of networks overlapping with limbic RSNs. STN-DBS significantly increased the occupancy of an orbitofrontal limbic subsystem with respect to both DBS OFF ( p = 0.0057) and 49 age-matched healthy controls ( p = 0.0033). Occupancy of a diffuse limbic RSN was increased with STN-DBS OFF when compared with healthy controls ( p = 0.021), but not when STN-DBS was ON, which indicates rebalancing of this network. These results highlight the modulatory effect of STN-DBS on components of the limbic system, particularly within the orbitofrontal cortex, a structure associated with reward processing. These results reinforce the value of quantitative biomarkers of RSN activity in evaluating the disseminated impact of brain stimulation techniques and the personalization of therapeutic strategies. Author Summary This article addresses a burning question regarding stimulation strategies to rebalance brain network dynamics. Using a rare fMRI dataset of Parkinson’s disease patients implanted with deep brain stimulation, we report evidence of network-specific modulatory effects in the dynamics of resting-state networks. In summary, we found that Leading Eigenvector Dynamics Analysis (LEiDA) successfully identified all seven reference resting-state networks in participants with subthalamic deep brain stimulation (STN-DBS). In particular, STN-DBS increases resting-state orbitofrontal cortex activity. STN-DBS also normalizes wider resting-state limbic network activity, and this correlated with depressive symptoms in pre- and post-operative assessments. The work is limited by a low number of participants, and the retrospective nature of this work, but provides evidence that STN-DBS modulates limbic network occupancy in real time.

Spontaneous brain activity changes across states of consciousness. A particular consciousness-mediated configuration is the anticorrelations between the default mode network and other brain regions. What this antagonistic organization implies about consciousness to date remains inconclusive. In this Perspective Article, we propose that anticorrelations are the physiological expression of the concept of segregation, namely the brain’s capacity to show selectivity in the way areas will be functionally connected. We postulate that this effect is mediated by the process of neural inhibition, by regulating global and local inhibitory activity. While recognizing that this effect can also result from other mechanisms, neural inhibition helps the understanding of how network metastability is affected after disrupting local and global neural balance. In combination with relevant theories of consciousness, we suggest that anticorrelations are a physiological prior that can work as a marker of preserved consciousness. We predict that if the brain is not in a state to host anticorrelations, then most likely the individual does not entertain subjective experience. We believe that this link between anticorrelations and the underlying physiology will help not only to comprehend how consciousness happens, but also conceptualize effective interventions for treating consciousness disorders in which anticorrelations seem particularly affected. Author Summary The fMRI resting paradigm can quantify brain function by surpassing communication and sophisticated setups, hence helping to infer consciousness in individuals who are unable to communicate with their environment. A particular consciousness-mediated rsfMRI configuration is that of functional anticorrelations, that is, the antagonistic relationship between a specific set of brain regions. We suggest that anticorrelations are a key physiological prior, without which consciousness cannot be supported, because the brain cannot segregate how regions get connected. We postulate that segregation is possible thanks to neural inhibition, by regulating global and local inhibitory activity. We believe that the link between anticorrelations and the underlying physiology can help not only to comprehend how consciousness happens, but also conceptualize effective interventions for treating its disorders.

Psychedelic drugs show promise as safe and effective treatments for neuropsychiatric disorders, yet their mechanisms of action are not fully understood. A fundamental hypothesis is that psychedelics work by dose-dependently changing the functional hierarchy of brain dynamics, but it is unclear whether different psychedelics act similarly. Here, we investigated the changes in the brain’s functional hierarchy associated with two different psychedelics (LSD and psilocybin). Using a novel turbulence framework, we were able to determine the vorticity, that is, the local level of synchronization, that allowed us to extend the standard global time-based measure of metastability to become a local-based measure of both space and time. This framework produced detailed signatures of turbulence-based hierarchical change for each psychedelic drug, revealing consistent and discriminate effects on a higher level network, that is, the default mode network. Overall, our findings directly support a prior hypothesis that psychedelics modulate (i.e., “compress”) the functional hierarchy and provide a quantification of these changes for two different psychedelics. Implications for therapeutic applications of psychedelics are discussed. Author Summary Significant progress has been made in understanding the effects of psychedelics on brain function. One of the main hypotheses is that psychedelics work by changing the functional hierarchy of brain dynamics in a dose-dependent manner, modulating the encoding of the precision of priors, beliefs, or assumptions in the brain. We used a novel turbulence framework to investigate the changes in the brain’s functional hierarchy associated with two different psychedelics (LSD and psilocybin). This framework produced detailed signatures of turbulence-based hierarchical change for each psychedelic drug, revealing consistent and discriminate effects on a higher level network, that is, the default mode network.

Previous computational models have related spontaneous resting-state brain activity with local excitatory–inhibitory balance in neuronal populations. However, how underlying neurotransmitter kinetics associated with E–I balance govern resting-state spontaneous brain dynamics remains unknown. Understanding the mechanisms by virtue of which fluctuations in neurotransmitter concentrations, a hallmark of a variety of clinical conditions, relate to functional brain activity is of critical importance. We propose a multiscale dynamic mean field (MDMF) model—a system of coupled differential equations for capturing the synaptic gating dynamics in excitatory and inhibitory neural populations as a function of neurotransmitter kinetics. Individual brain regions are modeled as population of MDMF and are connected by realistic connection topologies estimated from diffusion tensor imaging data. First, MDMF successfully predicts resting-state functional connectivity. Second, our results show that optimal range of glutamate and GABA neurotransmitter concentrations subserve as the dynamic working point of the brain, that is, the state of heightened metastability observed in empirical blood-oxygen-level-dependent signals. Third, for predictive validity the network measures of segregation (modularity and clustering coefficient) and integration (global efficiency and characteristic path length) from existing healthy and pathological brain network studies could be captured by simulated functional connectivity from an MDMF model. Author Summary How changes in neurotransmitter kinetics impact the organization of large-scale neurocognitive networks is an open question in neuroscience. Here, we propose a multiscale dynamic mean field (MDMF) model that incorporates biophysically realistic kinetic parameters of receptor binding in a dynamic mean field model and captures brain dynamics from the “whole brain.” MDMF could reliably reproduce the resting-state brain functional connectivity patterns. Further employing graph theoretic methods, MDMF could qualitatively explain the idiosyncrasies of network integration and segregation measures reported by previous clinical studies.

Neuroimaging techniques are now widely used to study human cognition. The functional associations between brain areas have become a standard proxy to describe how cognitive processes are distributed across the brain network. Among the many analysis tools available, dynamic models of brain activity have been developed to overcome the limitations of original connectivity measures such as functional connectivity. This goes in line with the many efforts devoted to the assessment of directional interactions between brain areas from the observed neuroimaging activity. This opinion article provides an overview of our model-based whole-brain effective connectivity to analyze fMRI data, while discussing the pros and cons of our approach with respect to other established approaches. Our framework relies on the multivariate Ornstein-Uhlenbeck (MOU) process and is thus referred to as MOU-EC. Once tuned, the model provides a directed connectivity estimate that reflects the dynamical state of BOLD activity, which can be used to explore cognition. We illustrate this approach using two applications on task-evoked fMRI data. First, as a connectivity measure, MOU-EC can be used to extract biomarkers for task-specific brain coordination, understood as the patterns of areas exchanging information. The multivariate nature of connectivity measures raises several challenges for whole-brain analysis, for which machine-learning tools present some advantages over statistical testing. Second, we show how to interpret changes in MOU-EC connections in a collective and model-based manner, bridging with network analysis. Our framework provides a comprehensive set of tools that open exciting perspectives to study distributed cognition, as well as neuropathologies. Author Summary Brain connectivity measures have been increasingly used to study cognition and neuropathologies with neuroimaging data. Many methods have been developed with particular objectives, in particular predictability to obtain biomarkers (i.e., which brain regions exhibit changes in their interactions across conditions) and interpretability (relating the connectivity measures back to the biology). In this article we present our framework for whole-brain effective connectivity that aims to find the equilibrium between these two desired properties, relying on a dynamic model that can be fitted to fMRI time series. Meanwhile, we compare it with previous work. Concretely, we show how machine learning can be used to extract biomarkers and how network-oriented analysis can be used to interpret the model estimates in terms of brain dynamics.