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

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

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

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

Time-resolved coupling between connectome harmonics and subjective experience under the psychedelic DMT.

Exploring the intricate relationship between brain's structure and function, and how this affects subjective experience is a fundamental pursuit in neuroscience. Psychedelic substances offer a unique insight into the influences of specific neurotransmitter systems on perception, cognition and consciousness. Specifically, their impact on brain function propagates across the structural connectome - a network of white matter pathways linking different regions. To comprehensively grasp the effects of psychedelic compounds on brain function, we used a theoretically rigorous framework known as connectome harmonic decomposition. This framework provides a robust method to characterize how brain function intricately depends on the organized network structure of the human connectome. We show that the connectome harmonic repertoire under DMT is reshaped in line with other reported psychedelic compounds - psilocybin, LSD and ketamine. Furthermore, we show that the repertoire entropy of connectome harmonics increases under DMT, as with those other psychedelics. Importantly, we demonstrate for the first time that measures of energy spectrum difference and repertoire entropy of connectome harmonics indexes the intensity of subjective experience of the participants in a time-resolved manner reflecting close coupling between connectome harmonics and subjective experience.

The flattening of spacetime hierarchy of the N,N-dimethyltryptamine brain state is characterized by harmonic decomposition of spacetime (HADES) framework.

The human brain is a complex system, whose activity exhibits flexible and continuous reorganization across space and time. The decomposition of whole-brain recordings into harmonic modes has revealed a repertoire of gradient-like activity patterns associated with distinct brain functions. However, the way these activity patterns are expressed over time with their changes in various brain states remains unclear. Here, we investigate healthy participants taking the serotonergic psychedelic N,N-dimethyltryptamine (DMT) with the Harmonic Decomposition of Spacetime (HADES) framework that can characterize how different harmonic modes defined in space are expressed over time. HADES demonstrates significant decreases in contributions across most low-frequency harmonic modes in the DMT-induced brain state. When normalizing the contributions by condition (DMT and non-DMT), we detect a decrease specifically in the second functional harmonic, which represents the uni- to transmodal functional hierarchy of the brain, supporting the leading hypothesis that functional hierarchy is changed in psychedelics. Moreover, HADES' dynamic spacetime measures of fractional occupancy, life time and latent space provide a precise description of the significant changes of the spacetime hierarchical organization of brain activity in the psychedelic state.

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.

Functional hierarchies in brain dynamics characterized by signal reversibility in ferret cortex.

Brain signal irreversibility has been shown to be a promising approach to study neural dynamics. Nevertheless, the relation with cortical hierarchy and the influence of different electrophysiological features is not completely understood. In this study, we recorded local field potentials (LFPs) during spontaneous behavior, including awake and sleep periods, using custom micro-electrocorticographic (µECoG) arrays implanted in ferrets. In contrast to humans, ferrets remain less time in each state across the sleep-wake cycle. We deployed a diverse set of metrics in order to measure the levels of complexity of the different behavioral states. In particular, brain irreversibility, which is a signature of non-equilibrium dynamics, captured by the arrow of time of the signal, revealed the hierarchical organization of the ferret's cortex. We found different signatures of irreversibility and functional hierarchy of large-scale dynamics in three different brain states (active awake, quiet awake, and deep sleep), showing a lower level of irreversibility in the deep sleep stage, compared to the other. Irreversibility also allowed us to disentangle the influence of different cortical areas and frequency bands in this process, showing a predominance of the parietal cortex and the theta band. Furthermore, when inspecting the embedded dynamic through a Hidden Markov Model, the deep sleep stage was revealed to have a lower switching rate and lower entropy production. These results suggest functional hierarchies in organization that can be revealed through thermodynamic features and information theory metrics.

A low dimensional embedding of brain dynamics enhances diagnostic accuracy and behavioral prediction in stroke.

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

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.

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

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

Multi-modal and multi-model interrogation of large-scale functional brain networks.

Existing whole-brain models are generally tailored to the modelling of a particular data modality (e.g., fMRI or MEG/EEG). We propose that despite the differing aspects of neural activity each modality captures, they originate from shared network dynamics. Building on the universal principles of self-organising delay-coupled nonlinear systems, we aim to link distinct features of brain activity - captured across modalities - to the dynamics unfolding on a macroscopic structural connectome. To jointly predict connectivity, spatiotemporal and transient features of distinct signal modalities, we consider two large-scale models - the Stuart Landau and Wilson and Cowan models - which generate short-lived 40 Hz oscillations with varying levels of realism. To this end, we measure features of functional connectivity and metastable oscillatory modes (MOMs) in fMRI and MEG signals - and compare them against simulated data. We show that both models can represent MEG functional connectivity (FC), functional connectivity dynamics (FCD) and generate MOMs to a comparable degree. This is achieved by adjusting the global coupling and mean conduction time delay and, in the WC model, through the inclusion of balance between excitation and inhibition. For both models, the omission of delays dramatically decreased the performance. For fMRI, the SL model performed worse for FCD and MOMs, highlighting the importance of balanced dynamics for the emergence of spatiotemporal and transient patterns of ultra-slow dynamics. Notably, optimal working points varied across modalities and no model was able to achieve a correlation with empirical FC higher than 0.4 across modalities for the same set of parameters. Nonetheless, both displayed the emergence of FC patterns that extended beyond the constraints of the anatomical structure. Finally, we show that both models can generate MOMs with empirical-like properties such as size (number of brain regions engaging in a mode) and duration (continuous time interval during which a mode appears). Our results demonstrate the emergence of static and dynamic properties of neural activity at different timescales from networks of delay-coupled oscillators at 40 Hz. Given the higher dependence of simulated FC on the underlying structural connectivity, we suggest that mesoscale heterogeneities in neural circuitry may be critical for the emergence of parallel cross-modal functional networks and should be accounted for in future modelling endeavours.

The lack of temporal brain dynamics asymmetry as a signature of impaired consciousness states.

Life is a constant battle against equilibrium. From the cellular level to the macroscopic scale, living organisms as dissipative systems require the violation of their detailed balance, i.e. metabolic enzymatic reactions, in order to survive. We present a framework based on temporal asymmetry as a measure of non-equilibrium. By means of statistical physics, it was discovered that temporal asymmetries establish an arrow of time useful for assessing the reversibility in human brain time series. Previous studies in human and non-human primates have shown that decreased consciousness states such as sleep and anaesthesia result in brain dynamics closer to the equilibrium. Furthermore, there is growing interest in the analysis of brain symmetry based on neuroimaging recordings and since it is a non-invasive technique, it can be extended to different brain imaging modalities and applied at different temporo-spatial scales. In the present study, we provide a detailed description of our methodological approach, paying special attention to the theories that motivated this work. We test, for the first time, the reversibility analysis in human functional magnetic resonance imaging data in patients suffering from disorder of consciousness. We verify that the tendency of a decrease in the asymmetry of the brain signal together with the decrease in non-stationarity are key characteristics of impaired consciousness states. We expect that this work will open the way for assessing biomarkers for patients' improvement and classification, as well as motivating further research on the mechanistic understanding underlying states of impaired consciousness.

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

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

Distributed harmonic patterns of structure-function dependence orchestrate human consciousness.

A central question in neuroscience is how consciousness arises from the dynamic interplay of brain structure and function. Here we decompose functional MRI signals from pathological and pharmacologically-induced perturbations of consciousness into distributed patterns of structure-function dependence across scales: the harmonic modes of the human structural connectome. We show that structure-function coupling is a generalisable indicator of consciousness that is under bi-directional neuromodulatory control. We find increased structure-function coupling across scales during loss of consciousness, whether due to anaesthesia or brain injury, capable of discriminating between behaviourally indistinguishable sub-categories of brain-injured patients, tracking the presence of covert consciousness. The opposite harmonic signature characterises the altered state induced by LSD or ketamine, reflecting psychedelic-induced decoupling of brain function from structure and correlating with physiological and subjective scores. Overall, connectome harmonic decomposition reveals how neuromodulation and the network architecture of the human connectome jointly shape consciousness and distributed functional activation across scales.

Dynamic sensitivity analysis: Defining personalised strategies to drive brain state transitions via whole brain modelling.

Traditionally, in neuroimaging, model-free analyses are used to find significant differences between brain states via signal detection theory. Depending on the a priori assumptions about the underlying data, different spatio-temporal features can be analysed. Alternatively, model-based techniques infer features from the data and compare significance from model parameters. However, to assess transitions from one brain state to another remains a challenge in current paradigms. Here, we introduce a "Dynamic Sensitivity Analysis" framework that quantifies transitions between brain states in terms of stimulation ability to rebalance spatio-temporal brain activity towards a target state such as healthy brain dynamics. In practice, it means building a whole-brain model fitted to the spatio-temporal description of brain dynamics, and applying systematic stimulations in-silico to assess the optimal strategy to drive brain dynamics towards a target state. Further, we show how Dynamic Sensitivity Analysis extends to various brain stimulation paradigms, ultimately contributing to improving the efficacy of personalised clinical interventions.

Understanding brain states across spacetime informed by whole-brain modelling.

In order to survive in a complex environment, the human brain relies on the ability to flexibly adapt ongoing behaviour according to intrinsic and extrinsic signals. This capability has been linked to specific whole-brain activity patterns whose relative stability (order) allows for consistent functioning, supported by sufficient intrinsic instability needed for optimal adaptability. The emergent, spontaneous balance between order and disorder in brain activity over spacetime underpins distinct brain states. For example, depression is characterized by excessively rigid, highly ordered states, while psychedelics can bring about more disordered, sometimes overly flexible states. Recent developments in systems, computational and theoretical neuroscience have started to make inroads into the characterization of such complex dynamics over space and time. Here, we review recent insights drawn from neuroimaging and whole-brain modelling motivating using mechanistic principles from dynamical system theory to study and characterize brain states. We show how different healthy and altered brain states are associated to characteristic spacetime dynamics which in turn may offer insights that in time can inspire new treatments for rebalancing brain states in disease. This article is part of the theme issue 'Emergent phenomena in complex physical and socio-technical systems: from cells to societies'.

Effects of classic psychedelic drugs on turbulent signatures in brain dynamics.

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.

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.

Metastable oscillatory modes emerge from synchronization in the brain spacetime connectome.

A rich repertoire of oscillatory signals is detected from human brains with electro- and magnetoencephalography (EEG/MEG). However, the principles underwriting coherent oscillations and their link with neural activity remain under debate. Here, we revisit the mechanistic hypothesis that transient brain rhythms are a signature of metastable synchronization, occurring at reduced collective frequencies due to delays between brain areas. We consider a system of damped oscillators in the presence of background noise - approximating the short-lived gamma-frequency oscillations generated within neuronal circuits - coupled according to the diffusion weighted tractography between brain areas. Varying the global coupling strength and conduction speed, we identify a critical regime where spatially and spectrally resolved metastable oscillatory modes (MOMs) emerge at sub-gamma frequencies, approximating the MEG power spectra from 89 healthy individuals at rest. Further, we demonstrate that the frequency, duration, and scale of MOMs - as well as the frequency-specific envelope functional connectivity - can be controlled by global parameters, while the connectome structure remains unchanged. Grounded in the physics of delay-coupled oscillators, these numerical analyses demonstrate how interactions between locally generated fast oscillations in the connectome spacetime structure can lead to the emergence of collective brain rhythms organized in space and time.

Stable biomarker identification for predicting schizophrenia in the human connectome.

Schizophrenia, as a psychiatric disorder, has recognized brain alterations both at the structural and at the functional magnetic resonance imaging level. The developing field of connectomics has attracted much attention as it allows researchers to take advantage of powerful tools of network analysis in order to study structural and functional connectivity abnormalities in schizophrenia. Many methods have been proposed to identify biomarkers in schizophrenia, focusing mainly on improving the classification performance or performing statistical comparisons between groups. However, the stability of biomarkers selection has been for long overlooked in the connectomics field. In this study, we follow a machine learning approach where the identification of biomarkers is addressed as a feature selection problem for a classification task. We perform a recursive feature elimination and support vector machines (RFE-SVM) approach to identify the most meaningful biomarkers from the structural, functional, and multi-modal connectomes of healthy controls and patients. Furthermore, the stability of the retrieved biomarkers is assessed across different subsamplings of the dataset, allowing us to identify the affected core of the pathology. Considering our technique altogether, it demonstrates a principled way to achieve both accurate and stable biomarkers while highlighting the importance of multi-modal approaches to brain pathology as they tend to reveal complementary information.

A brief exposure to rightward prismatic adaptation changes resting-state network characteristics of the ventral attentional system.

A brief session of rightward prismatic adaptation (R-PA) has been shown to alleviate neglect symptoms in patients with right hemispheric damage, very likely by switching hemispheric dominance of the ventral attentional network (VAN) from the right to the left and by changing task-related activity within the dorsal attentional network (DAN). We have investigated this very rapid change in functional organisation with a network approach by comparing resting-state connectivity before and after a brief exposure i) to R-PA (14 normal subjects; experimental condition) or ii) to plain glasses (12 normal subjects; control condition). A whole brain analysis (comprising 129 regions of interest) highlighted R-PA-induced changes within a bilateral, fronto-temporal network, which consisted of 13 nodes and 11 edges; all edges involved one of 4 frontal nodes, which were part of VAN. The analysis of network characteristics within VAN and DAN revealed a R-PA-induced decrease in connectivity strength between nodes and a decrease in local efficiency within VAN but not within DAN. These results indicate that the resting-state connectivity configuration of VAN is modulated by R-PA, possibly by decreasing its modularity.