Uncovering the biological basis of control energy: Structural and metabolic correlates of energy inefficiency in temporal lobe epilepsy

A network control theory pipeline for studying the dynamics of the structural connectome

Network control theory (NCT) is a simple and powerful tool for studying how network topology informs and constrains the dynamics of a system. Compared to other structure-function coupling approaches, the strength of NCT lies in its capacity to predict the patterns of external control signals that may alter the dynamics of a system in a desired way. An interesting development for NCT in the neuroscience field is its application to study behavior and mental health symptoms. To date, NCT has been validated to study different aspects of the human structural connectome. NCT outputs can be monitored throughout developmental stages to study the effects of connectome topology on neural dynamics and, separately, to test the coherence of empirical datasets with brain function and stimulation. Here, we provide a comprehensive pipeline for applying NCT to structural connectomes by following two procedures. The main procedure focuses on computing the control energy associated with the transitions between specific neural activity states. The second procedure focuses on computing average controllability, which indexes nodes' general capacity to control the dynamics of the system. We provide recommendations for comparing NCT outputs against null network models, and we further support this approach with a Python-based software package called 'network control theory for python'. The procedures in this protocol are appropriate for users with a background in network neuroscience and experience in dynamical systems theory.

Atypical connectome topography and signal flow in temporal lobe epilepsy

Temporal lobe epilepsy (TLE) is the most common pharmaco-resistant epilepsy in adults. While primarily associated with mesiotemporal pathology, recent evidence suggests that brain alterations in TLE extend beyond the paralimbic epicenter and impact macroscale function and cognitive functions, particularly memory. Using connectome-wide manifold learning and generative models of effective connectivity, we examined functional topography and directional signal flow patterns between large-scale neural circuits in TLE at rest. Studying a multisite cohort of 95 patients with TLE and 95 healthy controls, we observed atypical functional topographies in the former group, characterized by reduced differentiation between sensory and transmodal association cortices, with most marked effects in bilateral temporo-limbic and ventromedial prefrontal cortices. These findings were consistent across all study sites, present in left and right lateralized patients, and validated in a subgroup of patients with histopathological validation of mesiotemporal sclerosis and post-surgical seizure freedom. Moreover, they were replicated in an independent cohort of 30 TLE patients and 40 healthy controls. Further analyses demonstrated that reduced differentiation related to decreased functional signal flow into and out of temporolimbic cortical systems and other brain networks. Parallel analyses of structural and diffusion-weighted MRI data revealed that topographic alterations were independent of TLE-related cortical thinning but partially mediated by white matter microstructural changes that radiated away from paralimbic circuits. Finally, we found a strong association between the degree of functional alterations and behavioral markers of memory dysfunction. Our work illustrates the complex landscape of macroscale functional imbalances in TLE, which can serve as intermediate markers bridging microstructural changes and cognitive impairment.

A Human Brain Map of Mitochondrial Respiratory Capacity and Diversity

Mitochondrial oxidative phosphorylation (OxPhos) powers brain activity1,2, and mitochondrial defects are linked to neurodegenerative and neuropsychiatric disorders3,4, underscoring the need to define the brain's molecular energetic landscape5-10. To bridge the cognitive neuroscience and cell biology scale gap, we developed a physical voxelization approach to partition a frozen human coronal hemisphere section into 703 voxels comparable to neuroimaging resolution (3×3×3 mm). In each cortical and subcortical brain voxel, we profiled mitochondrial phenotypes including OxPhos enzyme activities, mitochondrial DNA and volume density, and mitochondria-specific respiratory capacity. We show that the human brain contains a diversity of mitochondrial phenotypes driven by both topology and cell types. Compared to white matter, grey matter contains >50% more mitochondria. We show that the more abundant grey matter mitochondria also are biochemically optimized for energy transformation, particularly among recently evolved cortical brain regions. Scaling these data to the whole brain, we created a backward linear regression model integrating several neuroimaging modalities11, thereby generating a brain-wide map of mitochondrial distribution and specialization that predicts mitochondrial characteristics in an independent brain region of the same donor brain. This new approach and the resulting MitoBrainMap of mitochondrial phenotypes provide a foundation for exploring the molecular energetic landscape that enables normal brain functions, relating it to neuroimaging data, and defining the subcellular basis for regionalized brain processes relevant to neuropsychiatric and neurodegenerative disorders.

A Human Brain Map of Mitochondrial Respiratory Capacity and Diversity

Mitochondrial oxidative phosphorylation (OxPhos) powers brain activity1,2, and mitochondrial defects are linked to neurodegenerative and neuropsychiatric disorders3,4, underscoring the need to define the brain's molecular energetic landscape5-10. To bridge the cognitive neuroscience and cell biology scale gap, we developed a physical voxelization approach to partition a frozen human coronal hemisphere section into 703 voxels comparable to neuroimaging resolution (3×3×3 mm). In each cortical and subcortical brain voxel, we profiled mitochondrial phenotypes including OxPhos enzyme activities, mitochondrial DNA and volume density, and mitochondria-specific respiratory capacity. We show that the human brain contains a diversity of mitochondrial phenotypes driven by both topology and cell types. Compared to white matter, grey matter contains >50% more mitochondria. We show that the more abundant grey matter mitochondria also are biochemically optimized for energy transformation, particularly among recently evolved cortical brain regions. Scaling these data to the whole brain, we created a backward linear regression model integrating several neuroimaging modalities11, thereby generating a brain-wide map of mitochondrial distribution and specialization that predicts mitochondrial characteristics in an independent brain region of the same donor brain. This new approach and the resulting MitoBrainMap of mitochondrial phenotypes provide a foundation for exploring the molecular energetic landscape that enables normal brain functions, relating it to neuroimaging data, and defining the subcellular basis for regionalized brain processes relevant to neuropsychiatric and neurodegenerative disorders.

Indirect structural changes and reduced controllability after temporal lobe epilepsy resection

OBJECTIVE

To understand the potential behavioral and cognitive effects of mesial temporal resection for temporal lobe epilepsy (TLE) a method is required to characterize network-wide functional alterations caused by a discrete structural disconnection. The objective of this study was to investigate network-wide alterations in brain dynamics of patients with TLE before and after surgical resection of the seizure focus using average regional controllability (ARC), a measure of the ability of a node to influence network dynamics.

METHODS

Diffusion-weighted imaging (DWI) data were acquired in 27 patients with drug-resistant unilateral mesial TLE who underwent selective amygdalohippocampectomy. Imaging data were acquired before and after surgery and a presurgical and postsurgical structural connectome was generated from whole-brain tractography. Edge-wise strength, node strength, and node ARC were compared before and after surgery. Direct and indirect edge-wise strength changes were identified using patient-specific simulated resections. Direct edges were defined as primary edges disconnected by the resection zone itself. Indirect edges were secondary measured edge strength changes. Changes in node strength and ARC were then related to both direct and indirect edge changes.

RESULTS

We found nodes with significant postsurgical changes in both node strength and ARC surrounding the resection zone (paired t tests, p < .05, Bonferroni corrected). ARC identified additional postsurgical changes in nodes outside of the resection zone within the ipsilateral occipital lobe, which were associated with indirect edge-wise strength changes of the postsurgical network (Fisher's exact test, p < .001). These indirect edge-wise changes were facilitated through the "hub" nodes including the thalamus, putamen, insula, and precuneus.

SIGNIFICANCE

Discrete network disconnection from TLE resection results in widespread structural and functional changes not predicted by disconnection alone. These can be well characterized by dynamic controllability measures such as ARC and may be useful for investigating changes in brain function that may contribute to seizure recurrence and behavioral or cognitive changes after surgery.

Controlling the human connectome with spatially diffuse input signals

The human brain is never at "rest"; its activity is constantly fluctuating over time, transitioning from one brain state-a whole-brain pattern of activity-to another. Network control theory offers a framework for understanding the effort - energy - associated with these transitions. One branch of control theory that is especially useful in this context is "optimal control", in which input signals are used to selectively drive the brain into a target state. Typically, these inputs are introduced independently to the nodes of the network (each input signal is associated with exactly one node). Though convenient, this input strategy ignores the continuity of cerebral cortex - geometrically, each region is connected to its spatial neighbors, allowing control signals, both exogenous and endogenous, to spread from their foci to nearby regions. Additionally, the spatial specificity of brain stimulation techniques is limited, such that the effects of a perturbation are measurable in tissue surrounding the stimulation site. Here, we adapt the network control model so that input signals have a spatial extent that decays exponentially from the input site. We show that this more realistic strategy takes advantage of spatial dependencies in structural connectivity and activity to reduce the energy (effort) associated with brain state transitions. We further leverage these dependencies to explore near-optimal control strategies such that, on a per-transition basis, the number of input signals required for a given control task is reduced, in some cases by two orders of magnitude. This approximation yields network-wide maps of input site density, which we compare to an existing database of functional, metabolic, genetic, and neurochemical maps, finding a close correspondence. Ultimately, not only do we propose a more efficient framework that is also more adherent to well-established brain organizational principles, but we also posit neurobiologically grounded bases for optimal control.

Altered structural connectivity and functional brain dynamics in individuals with heavy alcohol use

Heavy alcohol use and its associated conditions, such as alcohol use disorder (AUD), impact millions of individuals worldwide. While our understanding of the neurobiological correlates of AUD has evolved substantially, we still lack models incorporating whole-brain neuroanatomical, functional, and pharmacological information under one framework. Here, we utilize diffusion and functional magnetic resonance imaging to investigate alterations to brain dynamics in N = 130 individuals with a high amount of current alcohol use. We compared these alcohol using individuals to N = 308 individuals with minimal use of any substances. We find that individuals with heavy alcohol use had less dynamic and complex brain activity, and through leveraging network control theory, had increased control energy to complete transitions between activation states. Further, using separately acquired positron emission tomography (PET) data, we deploy an in silico evaluation demonstrating that decreased D2 receptor levels, as found previously in individuals with AUD, may relate to our observed findings. This work demonstrates that whole-brain, multimodal imaging information can be combined under a network control framework to identify and evaluate neurobiological correlates and mechanisms of AUD.

Network controllability of structural connectomes in the neonatal brain

White matter connectivity supports diverse cognitive demands by efficiently constraining dynamic brain activity. This efficiency can be inferred from network controllability, which represents the ease with which the brain moves between distinct mental states based on white matter connectivity. However, it remains unclear how brain networks support diverse functions at birth, a time of rapid changes in connectivity. Here, we investigate the development of network controllability during the perinatal period and the effect of preterm birth in 521 neonates. We provide evidence that elements of controllability are exhibited in the infant's brain as early as the third trimester and develop rapidly across the perinatal period. Preterm birth disrupts the development of brain networks and altered the energy required to drive state transitions at different levels. In addition, controllability at birth is associated with cognitive ability at 18 months. Our results suggest network controllability develops rapidly during the perinatal period to support cognitive demands but could be altered by environmental impacts like preterm birth.

Using network control theory to study the dynamics of the structural connectome

Network control theory (NCT) is a simple and powerful tool for studying how network topology informs and constrains dynamics. Compared to other structure-function coupling approaches, the strength of NCT lies in its capacity to predict the patterns of external control signals that may alter dynamics in a desired way. We have extensively developed and validated the application of NCT to the human structural connectome. Through these efforts, we have studied (i) how different aspects of connectome topology affect neural dynamics, (ii) whether NCT outputs cohere with empirical data on brain function and stimulation, and (iii) how NCT outputs vary across development and correlate with behavior and mental health symptoms. In this protocol, we introduce a framework for applying NCT to structural connectomes following two main pathways. Our primary pathway focuses on computing the control energy associated with transitioning between specific neural activity states. Our second pathway focuses on computing average controllability, which indexes nodes' general capacity to control dynamics. We also provide recommendations for comparing NCT outputs against null network models. Finally, we support this protocol with a Python-based software package called network control theory for python (nctpy).

Dynamic network properties of the superior temporal gyrus mediate the impact of brain age gap on chronic aphasia severity

Brain structure deteriorates with aging and predisposes an individual to more severe language impairments (aphasia) after a stroke. However, the underlying mechanisms of this relation are not well understood. Here we use an approach to model brain network properties outside the stroke lesion, network controllability, to investigate relations among individualized structural brain connections, brain age, and aphasia severity in 93 participants with chronic post-stroke aphasia. Controlling for the stroke lesion size, we observe that lower average controllability of the posterior superior temporal gyrus (STG) mediates the relation between advanced brain aging and aphasia severity. Lower controllability of the left posterior STG signifies that activity in the left posterior STG is less likely to yield a response in other brain regions due to the topological properties of the structural brain networks. These results indicate that advanced brain aging among individuals with post-stroke aphasia is associated with disruption of dynamic properties of a critical language-related area, the STG, which contributes to worse aphasic symptoms. Because brain aging is variable among individuals with aphasia, our results provide further insight into the mechanisms underlying the variance in clinical trajectories in post-stroke aphasia.

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.

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.

Transitions between cognitive topographies: contributions of network structure, neuromodulation, and disease

Patterns of neural activity underlie human cognition. Transitions between these patterns are orchestrated by the brain's network architecture. What are the mechanisms linking network structure to cognitively relevant activation patterns? Here we implement principles of network control to investigate how the architecture of the human connectome shapes transitions between 123 experimentally defined cognitive activation maps (cognitive topographies) from the NeuroSynth meta-analytic engine. We also systematically incorporate neurotransmitter receptor density maps (18 receptors and transporters) and disease-related cortical abnormality maps (11 neurodegenerative, psychiatric and neurodevelopmental diseases; N = 17 000 patients, N = 22 000 controls). Integrating large-scale multimodal neuroimaging data from functional MRI, diffusion tractography, cortical morphometry, and positron emission tomography, we simulate how anatomically-guided transitions between cognitive states can be reshaped by pharmacological or pathological perturbation. Our results provide a comprehensive look-up table charting how brain network organisation and chemoarchitecture interact to manifest different cognitive topographies. This computational framework establishes a principled foundation for systematically identifying novel ways to promote selective transitions between desired cognitive topographies.

Genetic, individual, and familial risk correlates of brain network controllability in major depressive disorder

Many therapeutic interventions in psychiatry can be viewed as attempts to influence the brain's large-scale, dynamic network state transitions. Building on connectome-based graph analysis and control theory, Network Control Theory is emerging as a powerful tool to quantify network controllability-i.e., the influence of one brain region over others regarding dynamic network state transitions. If and how network controllability is related to mental health remains elusive. Here, from Diffusion Tensor Imaging data, we inferred structural connectivity and inferred calculated network controllability parameters to investigate their association with genetic and familial risk in patients diagnosed with major depressive disorder (MDD, n = 692) and healthy controls (n = 820). First, we establish that controllability measures differ between healthy controls and MDD patients while not varying with current symptom severity or remission status. Second, we show that controllability in MDD patients is associated with polygenic scores for MDD and psychiatric cross-disorder risk. Finally, we provide evidence that controllability varies with familial risk of MDD and bipolar disorder as well as with body mass index. In summary, we show that network controllability is related to genetic, individual, and familial risk in MDD patients. We discuss how these insights into individual variation of network controllability may inform mechanistic models of treatment response prediction and personalized intervention-design in mental health.

Network controllability measures of subnetworks: implications for neurosciences

Objective:Recent progress in network sciences has made it possible to apply key findings from control theory to the study of networks. Referred to as network control theory, this framework describes how the interactions between interconnected system elements and external energy sources, potentially constrained by different optimality criteria, result in complex network behavior. A typical example is the quantification of the functional role certain brain regions or symptoms play in shaping the temporal dynamics of brain activity or the clinical course of a disease, a property that is quantified in terms of the so-called controllability metrics. Critically though, contrary to the engineering context in which control theory was originally developed, a mathematical understanding of the network nodes and connections in neurosciences cannot be assumed. For instance, in the case of psychological systems such as those studied to understand psychiatric disorders, a potentially large set of related variables are unknown. As such, while the measures offered by network control theory would be mathematically correct, in that they can be calculated with high precision, they could have little translational values with respect to their putative role suggested by controllability metrics. It is therefore critical to understand if and how the controllability metrics estimated over subnetworks would deviate, if access to the complete set of variables, as is common in neurosciences, cannot be taken for granted.Approach:In this paper, we use a host of simulations based on synthetic as well as structural magnetic resonance imaging (MRI) data to study the potential deviation of controllability metrics in sub- compared to the full networks. Specifically, we estimate average- and modal-controllability, two of the most widely used controllability measures in neurosciences, in a large number of settings where we systematically vary network type, network size, and edge density.Main results:We find out, across all network types we test, that average and modal controllability are systematically, over- or underestimated depending on the number of nodes in the sub- and full network and the edge density.Significance:Finally, we provide formal theoretical proof that our observations generalize to any network type and discuss the ramifications of this systematic bias and potential solutions to alleviate the problem.

Towards a network control theory of electroconvulsive therapy response

Electroconvulsive Therapy (ECT) is arguably the most effective intervention for treatment-resistant depression. While large interindividual variability exists, a theory capable of explaining individual response to ECT remains elusive. To address this, we posit a quantitative, mechanistic framework of ECT response based on Network Control Theory (NCT). Then, we empirically test our approach and employ it to predict ECT treatment response. To this end, we derive a formal association between Postictal Suppression Index (PSI)-an ECT seizure quality index-and whole-brain modal and average controllability, NCT metrics based on white-matter brain network architecture, respectively. Exploiting the known association of ECT response and PSI, we then hypothesized an association between our controllability metrics and ECT response mediated by PSI. We formally tested this conjecture in N = 50 depressive patients undergoing ECT. We show that whole-brain controllability metrics based on pre-ECT structural connectome data predict ECT response in accordance with our hypotheses. In addition, we show the expected mediation effects via PSI. Importantly, our theoretically motivated metrics are at least on par with extensive machine learning models based on pre-ECT connectome data. In summary, we derived and tested a control-theoretic framework capable of predicting ECT response based on individual brain network architecture. It makes testable, quantitative predictions regarding individual therapeutic response, which are corroborated by strong empirical evidence. Our work might constitute a starting point for a comprehensive, quantitative theory of personalized ECT interventions rooted in control theory.