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.

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 non-reversibility and thus non-equilibrium in brain dynamics. When applied to large-scale HCP neuroimaging data from close to a thousand participants, the TENET framework suggests that the DMN plays a significant role in orchestrating the hierarchy, i.e. levels of non-reversibility, 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.

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.