A framework for quantifying the coupling between brain connectivity and heartbeat dynamics: Insights into the disrupted network physiology in Parkinson's disease.

Parkinson's disease (PD) often shows disrupted brain connectivity and autonomic dysfunctions, progressing alongside with motor and cognitive decline. Recently, PD has been linked to a reduced sensitivity to cardiac inputs, that is, cardiac interoception. Altogether, those signs suggest that PD causes an altered brain-heart connection whose mechanisms remain unclear. Our study aimed to explore the large-scale network disruptions and the neurophysiology of disrupted interoceptive mechanisms in PD. We focused on examining the alterations in brain-heart coupling in PD and their potential connection to motor symptoms. We developed a proof-of-concept method to quantify relationships between the co-fluctuations of brain connectivity and cardiac sympathetic and parasympathetic activities. We quantified the brain-heart couplings from electroencephalogram and electrocardiogram recordings from PD patients on and off dopaminergic medication, as well as in healthy individuals at rest. Our results show that the couplings of fluctuating alpha and gamma connectivity with cardiac sympathetic dynamics are reduced in PD patients, as compared to healthy individuals. Furthermore, we show that PD patients under dopamine medication recover part of the brain-heart coupling, in proportion with the reduced motor symptoms. Our proposal offers a promising approach to unveil the physiopathology of PD and promoting the development of new evaluation methods for the early stages of the disease.

Causality of cortical and cardiovascular activity during cyclic alternating pattern in non-rapid eye movement sleep.

The dynamic interplay between central and autonomic nervous system activities plays a pivotal role in orchestrating sleep. Macrostructural changes such as sleep-stage transitions or phasic, brief cortical events elicit fluctuations in neural outflow to the cardiovascular system, but the causal relationships between cortical and cardiovascular activities underpinning the microstructure of sleep are largely unknown. Here, we investigate cortical-cardiovascular interactions during the cyclic alternating pattern (CAP) of non-rapid eye movement sleep in a diverse set of overnight polysomnograms. We determine the Granger causality in both 507 CAP and 507 matched non-CAP sequences to assess the causal relationships between electroencephalography (EEG) frequency bands and respiratory and cardiovascular variables (heart period, respiratory period, pulse arrival time and pulse wave amplitude) during CAP. We observe a significantly stronger influence of delta activity on vascular variables during CAP sequences where slow, low-amplitude EEG activation phases (A1) dominate than during non-CAP sequences. We also show that rapid, high-amplitude EEG activation phases (A3) provoke a more pronounced change in autonomic activity than A1 and A2 phases. Our analysis provides the first evidence on the causal interplay between cortical and cardiovascular activities during CAP. Granger causality analysis may also be useful for probing the level of decoupling in sleep disorders. This article is part of the theme issue 'Advanced computation in cardiovascular physiology: new challenges and opportunities'.

Assessment of Outliers and Detection of Artifactual Network Segments Using Univariate and Multivariate Dispersion Entropy on Physiological Signals.

Network physiology has emerged as a promising paradigm for the extraction of clinically relevant information from physiological signals by moving from univariate to multivariate analysis, allowing for the inspection of interdependencies between organ systems. However, for its successful implementation, the disruptive effects of artifactual outliers, which are a common occurrence in physiological recordings, have to be studied, quantified, and addressed. Within the scope of this study, we utilize Dispersion Entropy (DisEn) to initially quantify the capacity of outlier samples to disrupt the values of univariate and multivariate features extracted with DisEn from physiological network segments consisting of synchronised, electroencephalogram, nasal respiratory, blood pressure, and electrocardiogram signals. The DisEn algorithm is selected due to its efficient computation and good performance in the detection of changes in signals for both univariate and multivariate time-series. The extracted features are then utilised for the training and testing of a logistic regression classifier in univariate and multivariate configurations in an effort to partially automate the detection of artifactual network segments. Our results indicate that outlier samples cause significant disruption in the values of extracted features with multivariate features displaying a certain level of robustness based on the number of signals formulating the network segments from which they are extracted. Furthermore, the deployed classifiers achieve noteworthy performance, where the percentage of correct network segment classification surpasses 95% in a number of experimental setups, with the effectiveness of each configuration being affected by the signal in which outliers are located. Finally, due to the increase in the number of features extracted within the framework of network physiology and the observed impact of artifactual samples in the accuracy of their values, the implementation of algorithmic steps capable of effective feature selection is highlighted as an important area for future research.

Information Transfer in Linear Multivariate Processes Assessed through Penalized Regression Techniques: Validation and Application to Physiological Networks.

The framework of information dynamics allows the dissection of the information processed in a network of multiple interacting dynamical systems into meaningful elements of computation that quantify the information generated in a target system, stored in it, transferred to it from one or more source systems, and modified in a synergistic or redundant way. The concepts of information transfer and modification have been recently formulated in the context of linear parametric modeling of vector stochastic processes, linking them to the notion of Granger causality and providing efficient tools for their computation based on the state-space (SS) representation of vector autoregressive (VAR) models. Despite their high computational reliability these tools still suffer from estimation problems which emerge, in the case of low ratio between data points available and the number of time series, when VAR identification is performed via the standard ordinary least squares (OLS). In this work we propose to replace the OLS with penalized regression performed through the Least Absolute Shrinkage and Selection Operator (LASSO), prior to computation of the measures of information transfer and information modification. First, simulating networks of several coupled Gaussian systems with complex interactions, we show that the LASSO regression allows, also in conditions of data paucity, to accurately reconstruct both the underlying network topology and the expected patterns of information transfer. Then we apply the proposed VAR-SS-LASSO approach to a challenging application context, i.e., the study of the physiological network of brain and peripheral interactions probed in humans under different conditions of rest and mental stress. Our results, which document the possibility to extract physiologically plausible patterns of interaction between the cardiovascular, respiratory and brain wave amplitudes, open the way to the use of our new analysis tools to explore the emerging field of Network Physiology in several practical applications.

Information Dynamics of the Brain, Cardiovascular and Respiratory Network during Different Levels of Mental Stress.

In this study, an analysis of brain, cardiovascular and respiratory dynamics was conducted combining information-theoretic measures with the Network Physiology paradigm during different levels of mental stress. Starting from low invasive recordings of electroencephalographic, electrocardiographic, respiratory, and blood volume pulse signals, the dynamical activity of seven physiological systems was probed with one-second time resolution measuring the time series of the δ , θ , α and ß brain wave amplitudes, the cardiac period (RR interval), the respiratory amplitude, and the duration of blood pressure wave propagation (pulse arrival time, PAT). Synchronous 5-min windows of these time series, obtained from 18 subjects during resting wakefulness (REST), mental stress induced by mental arithmetic (MA) and sustained attention induced by serious game (SG), were taken to describe the dynamics of the nodes composing the observed physiological network. Network activity and connectivity were then assessed in the framework of information dynamics computing the new information generated by each node, the information dynamically stored in it, and the information transferred to it from the other network nodes. Moreover, the network topology was investigated using directed measures of conditional information transfer and assessing their statistical significance. We found that all network nodes dynamically produce and store significant amounts of information, with the new information being prevalent in the brain systems and the information storage being prevalent in the peripheral systems. The transition from REST to MA was associated with an increase of the new information produced by the respiratory signal time series (RESP), and that from MA to SG with a decrease of the new information produced by PAT. Each network node received a significant amount of information from the other nodes, with the highest amount transferred to RR and the lowest transferred to δ , θ , α and ß . The topology of the physiological network underlying such information transfer was node- and state-dependent, with the peripheral subnetwork showing interactions from RR to PAT and between RESP and RR, PAT consistently across states, the brain subnetwork resulting more connected during MA, and the subnetwork of brain-peripheral interactions involving different brain rhythms in the three states and resulting primarily activated during MA. These results have both physiological relevance as regards the interpretation of central and autonomic effects on cardiovascular and respiratory variability, and practical relevance as regards the identification of features useful for the automatic distinction of different mental states.

The emergence of identity, agency and consciousness from the temporal dynamics of neural elaboration.

Identity-differentiating self from external reality-and agency-being the author of one's acts-are generally considered intrinsic properties of awareness and looked at as mental constructs generated by consciousness. Here a different view is proposed. All physiological systems display complex time-dependent regulations to adapt or anticipate external changes. To interact with rapid changes, an animal needs a nervous system capable of modelling and predicting (not simply representing) it. Different algorithms must be employed to predict the momentary location of an object based on sensory information (received with a delay), or to design in advance and direct the trajectory of movement. Thus, the temporal dynamics of external events and action must be handled in differential ways, thereby generating the distinction between self and non-self ("identity") as an intrinsic computational construct in neuronal elaboration. Handling time is not what neurons are designed for. Neuronal circuits are based on parallel processing: each bit of information diverges on many neurons, each of which combines it with many other data. Spike firing reports the likelihood that the specific pattern the neuron is designed to respond to is present in the incoming data. This organization seems designed to process synchronous datasets. However, since neural networks can introduce delays in processing, time sequences can be transformed into simultaneous patterns and analysed as such. This way predictive algorithms can be implemented, and continually improved through neuronal plasticity. To successfully interact with the external reality, the nervous system must model and predict, but also differentially handle perceptual functions or motor activity, by putting in register information that becomes available at different time moments. Also, to learn through positive/negative reinforcement, modelling must establish a causal relation between motor control and its consequences: the contrast between phase lag in perception and phase lead (and control) in motor programming produces the emergence of identity (discerning self from surrounding) and agency (control on actions) as necessary computational constructs to model reality. This does not require any form of awareness. In a brain, capable of producing awareness, these constructs may evolve from mere computational requirements into mental (conscious) constructs.

Feedback control of collective dynamics in an oscillator population with time-dependent connectivity.

We present a numerical study of pulsatile feedback-based control of synchrony level in a highly-interconnected oscillatory network. We focus on a nontrivial case when the system is close to the synchronization transition point and exhibits collective rhythm with strong amplitude modulation. We pay special attention to technical but essential steps like causal real-time extraction of the signal of interest from a noisy measurement and estimation of instantaneous phase and amplitude. The feedback loop's parameters are tuned automatically to suppress synchrony. Though the study is motivated by neuroscience, the results are relevant to controlling oscillatory activity in ensembles of various natures and, thus, to the rapidly developing field of network physiology.

Muscle Recruitment Strategies in a Redundant Task: Age Differences Through Network Analyses.

There are numerous studies comparing young and old adults in terms of muscle coordination in standard tasks (e.g., walking, reaching) and small variations of them. These tasks might hide differences: individuals would converge to similar behavior as they practice these throughout life. Also, we are unaware of studies that considered the muscle recruitment nested dynamics. For this reason, our study evaluated how young and old women coordinate and control the movement system while performing an unusual redundant motor control task through the network physiology approach. We acquired electromyographic signals from nine leg muscles of the dominant and non-dominant limbs during maximum voluntary isometric contractions (knee extension and flexion) and co-contraction bouts. Our results showed that young participants presented higher peak torque output, with similar EMG variability, compared to older participants. Considering firing rate frequencies, old and young women demonstrated different traits for network clustering and efficiency for the task. Age seems to affect muscle coordination at higher frequencies, even with a similar number of muscle synergies, indicating that younger women might have more integrated synergies than older women. The findings also point to differential muscle coordination adaptability.

Physiological network approach to prognosis in cirrhosis: A shifting paradigm.

Decompensated liver disease is complicated by multi-organ failure and poor prognosis. The prognosis of patients with liver failure often dictates clinical management. Current prognostic models have focused on biomarkers considered as individual isolated units. Network physiology assesses the interactions among multiple physiological systems in health and disease irrespective of anatomical connectivity and defines the influence or dependence of one organ system on another. Indeed, recent applications of network mapping methods to patient data have shown improved prediction of response to therapy or prognosis in cirrhosis. Initially, different physical markers have been used to assess physiological coupling in cirrhosis including heart rate variability, heart rate turbulence, and skin temperature variability measures. Further, the parenclitic network analysis was recently applied showing that organ systems connectivity is impaired in patients with decompensated cirrhosis and can predict mortality in cirrhosis independent of current prognostic models while also providing valuable insights into the associated pathological pathways. Moreover, network mapping also predicts response to intravenous albumin in patients hospitalized with decompensated cirrhosis. Thus, this review highlights the importance of evaluating decompensated cirrhosis through the network physiologic prism. It emphasizes the limitations of current prognostic models and the values of network physiologic techniques in cirrhosis.

Interoception, network physiology and the emergence of bodily self-awareness.

The interplay between the brain and interoceptive signals is key in maintaining internal balance and orchestrating neural dynamics, encompassing influences on perceptual and self-awareness. Central to this interplay is the differentiation between the external world, others and the self, a cornerstone in the construction of bodily self-awareness. This review synthesizes physiological and behavioral evidence illustrating how interoceptive signals can mediate or influence bodily self-awareness, by encompassing interactions with various sensory modalities. To deepen our understanding of the basis of bodily self-awareness, we propose a network physiology perspective. This approach explores complex neural computations across multiple nodes, shifting the focus from localized areas to large-scale neural networks. It examines how these networks operate in parallel with and adapt to changes in visceral activities. Within this framework, we propose to investigate physiological factors that disrupt bodily self-awareness, emphasizing the impact of interoceptive pathway disruptions, offering insights across several clinical contexts. This integrative perspective not only can enhance the accuracy of mental health assessments but also paves the way for targeted interventions.

Application of early warning signs to physiological contexts: a comparison of multivariate indices in patients on long-term hemodialysis.

Early warnings signs (EWSs) can anticipate abrupt changes in system state, known as "critical transitions," by detecting dynamic variations, including increases in variance, autocorrelation (AC), and cross-correlation. Numerous EWSs have been proposed; yet no consensus on which perform best exists. Here, we compared 15 multivariate EWSs in time series of 763 hemodialyzed patients, previously shown to present relevant critical transition dynamics. We calculated five EWSs based on AC, six on variance, one on cross-correlation, and three on AC and variance. We assessed their pairwise correlations, trends before death, and mortality predictive power, alone and in combination. Variance-based EWSs showed stronger correlations (r = 0.663 ± 0.222 vs. 0.170 ± 0.205 for AC-based indices) and a steeper increase before death. Two variance-based EWSs yielded HR95 > 9 (HR95 standing for a scale-invariant metric of hazard ratio), but combining them did not improve the area under the receiver-operating curve (AUC) much compared to using them alone (AUC = 0.798 vs. 0.796 and 0.791). Nevertheless, the AUC reached 0.825 when combining 13 indices. While some indicators did not perform overly well alone, their addition to the best performing EWSs increased the predictive power, suggesting that indices combination captures a broader range of dynamic changes occurring within the system. It is unclear whether this added benefit reflects measurement error of a unified phenomenon or heterogeneity in the nature of signals preceding critical transitions. Finally, the modest predictive performance and weak correlations among some indices call into question their validity, at least in this context.

Information theory reveals physiological manifestations of COVID-19 that correlate with symptom density of illness.

Algorithms for the detection of COVID-19 illness from wearable sensor devices tend to implicitly treat the disease as causing a stereotyped (and therefore recognizable) deviation from healthy physiology. In contrast, a substantial diversity of bodily responses to SARS-CoV-2 infection have been reported in the clinical milieu. This raises the question of how to characterize the diversity of illness manifestations, and whether such characterization could reveal meaningful relationships across different illness manifestations. Here, we present a framework motivated by information theory to generate quantified maps of illness presentation, which we term "manifestations," as resolved by continuous physiological data from a wearable device (Oura Ring). We test this framework on five physiological data streams (heart rate, heart rate variability, respiratory rate, metabolic activity, and sleep temperature) assessed at the time of reported illness onset in a previously reported COVID-19-positive cohort (N = 73). We find that the number of distinct manifestations are few in this cohort, compared to the space of all possible manifestations. In addition, manifestation frequency correlates with the rough number of symptoms reported by a given individual, over a several-day period prior to their imputed onset of illness. These findings suggest that information-theoretic approaches can be used to sort COVID-19 illness manifestations into types with real-world value. This proof of concept supports the use of information-theoretic approaches to map illness manifestations from continuous physiological data. Such approaches could likely inform algorithm design and real-time treatment decisions if developed on large, diverse samples.

Network approach reveals preferential T-cell and macrophage association with α-linked ß-cells in early stage of insulitis in NOD mice.

One of the challenges in studying islet inflammation-insulitis-is that it is a transient phenomenon. Traditional reporting of the insulitis progression is based on cumulative, donor-averaged values of leucocyte density in the vicinity of pancreatic islets, that hinder intra- and inter-islet heterogeneity of disease progression. Here, we aimed to understand why insulitis is non-uniform, often with peri-insulitis lesions formed on one side of an islet. To achieve this, we demonstrated the applicability of network theory in detangling intra-islet multi-cellular interactions during insulitis. Specifically, we asked the question "What is unique about regions of the islet that interact with immune cells first". This study utilized the non-obese diabetic mouse model of type one diabetes and examined the interplay among α-, ß-, T-cells, myeloid cells, and macrophages in pancreatic islets during the progression of insulitis. Disease evolution was tracked based on the T/ß cell ratio in individual islets. In the early stage, we found that immune cells are preferentially interacting with α-cell-rich regions of an islet. At the islet periphery α-linked ß-cells were found to be targeted significantly more compared to those without α-cell neighbors. Additionally, network analysis revealed increased T-myeloid, and T-macrophage interactions with all ß-cells.

Age-related breakdown in networks of inter-muscular coordination.

Assessing inter-muscular coordination in older adults is crucial, as it directly impacts an individual's ability for independent functioning, injury prevention, and active engagement in daily activities. However, the precise mechanisms by which distinct muscle fiber types synchronize their activity across muscles to generate coordinated movements in older adults remain unknown. Our objective is to investigate how distinct muscle groups dynamically synchronize with each other in young and older adults during exercise. Thirty-five young adults and nine older adults performed one bodyweight squat set until exhaustion. Simultaneous surface electromyography (sEMG) recordings were taken from the left and right vastus lateralis, and left and right erector spinae. To quantify inter-muscular coordination, we first obtained ten time series of sEMG band power for each muscle, representing the dynamics of different muscle fiber types. Next, we calculated the bivariate equal-time Pearson's cross-correlation for each pair of sEMG band power time series across all leg and back muscles. The main results show (i) an overall reduction in the degree of inter-muscular coordination, and (ii) increased stratification of the inter-muscular network in older adults compared to young adults. These findings suggest that as individuals age, the global inter-muscular network becomes less flexible and adaptable, hindering its ability to reorganize effectively in response to fatigue or other stimuli. This network approach opens new avenues for developing novel network-based markers to characterize multilevel inter-muscular interactions, which can help target functional deficits and potentially reduce the risk of falls and neuro-muscular injuries in older adults.

Emergence of metastability in frustrated oscillatory networks: the key role of hierarchical modularity.

Oscillatory complex networks in the metastable regime have been used to study the emergence of integrated and segregated activity in the brain, which are hypothesised to be fundamental for cognition. Yet, the parameters and the underlying mechanisms necessary to achieve the metastable regime are hard to identify, often relying on maximising the correlation with empirical functional connectivity dynamics. Here, we propose and show that the brain's hierarchically modular mesoscale structure alone can give rise to robust metastable dynamics and (metastable) chimera states in the presence of phase frustration. We construct unweighted 3-layer hierarchical networks of identical Kuramoto-Sakaguchi oscillators, parameterized by the average degree of the network and a structural parameter determining the ratio of connections between and within blocks in the upper two layers. Together, these parameters affect the characteristic timescales of the system. Away from the critical synchronization point, we detect the emergence of metastable states in the lowest hierarchical layer coexisting with chimera and metastable states in the upper layers. Using the Laplacian renormalization group flow approach, we uncover two distinct pathways towards achieving the metastable regimes detected in these distinct layers. In the upper layers, we show how the symmetry-breaking states depend on the slow eigenmodes of the system. In the lowest layer instead, metastable dynamics can be achieved as the separation of timescales between layers reaches a critical threshold. Our results show an explicit relationship between metastability, chimera states, and the eigenmodes of the system, bridging the gap between harmonic based studies of empirical data and oscillatory models.

Measures of the coupling between fluctuating brain network organization and heartbeat dynamics.

In recent years, there has been an increasing interest in studying brain-heart interactions. Methodological advancements have been proposed to investigate how the brain and the heart communicate, leading to new insights into some neural functions. However, most frameworks look at the interaction of only one brain region with heartbeat dynamics, overlooking that the brain has functional networks that change dynamically in response to internal and external demands. We propose a new framework for assessing the functional interplay between cortical networks and cardiac dynamics from noninvasive electrophysiological recordings. We focused on fluctuating network metrics obtained from connectivity matrices of EEG data. Specifically, we quantified the coupling between cardiac sympathetic-vagal activity and brain network metrics of clustering, efficiency, assortativity, and modularity. We validate our proposal using open-source datasets: one that involves emotion elicitation in healthy individuals, and another with resting-state data from patients with Parkinson's disease. Our results suggest that the connection between cortical network segregation and cardiac dynamics may offer valuable insights into the affective state of healthy participants, and alterations in the network physiology of Parkinson's disease. By considering multiple network properties, this framework may offer a more comprehensive understanding of brain-heart interactions. Our findings hold promise in the development of biomarkers for diagnostic and cognitive/motor function evaluation.

Resilience of the slow component in timescale-separated synchronized oscillators.

Physiological networks are usually made of a large number of biological oscillators evolving on a multitude of different timescales. Phase oscillators are particularly useful in the modelling of the synchronization dynamics of such systems. If the coupling is strong enough compared to the heterogeneity of the internal parameters, synchronized states might emerge where phase oscillators start to behave coherently. Here, we focus on the case where synchronized oscillators are divided into a fast and a slow component so that the two subsets evolve on separated timescales. We assess the resilience of the slow component by, first, reducing the dynamics of the fast one using Mori-Zwanzig formalism. Second, we evaluate the variance of the phase deviations when the oscillators in the two components are subject to noise with possibly distinct correlation times. From the general expression for the variance, we consider specific network structures and show how the noise transmission between the fast and slow components is affected. Interestingly, we find that oscillators that are among the most robust when there is only a single timescale, might become the most vulnerable when the system undergoes a timescale separation. We also find that layered networks seem to be insensitive to such timescale separations.

Chaos control in cardiac dynamics: terminating chaotic states with local minima pacing.

Current treatments of cardiac arrhythmias like ventricular fibrillation involve the application of a high-energy electric shock, that induces significant electrical currents in the myocardium and therefore involves severe side effects like possible tissue damage and post-traumatic stress. Using numerical simulations on four different models of 2D excitable media, this study demonstrates that low energy pulses applied shortly after local minima in the mean value of the transmembrane potential provide high success rates. We evaluate the performance of this approach for ten initial conditions of each model, ten spatially different stimuli, and different shock amplitudes. The investigated models of 2D excitable media cover a broad range of dominant frequencies and number of phase singularities, which demonstrates, that our findings are not limited to a specific kind of model or parameterization of it. Thus, we propose a method that incorporates the dynamics of the underlying system, even during pacing, and solely relies on a scalar observable, which is easily measurable in numerical simulations.

Application of physiological network mapping in the prediction of survival in critically ill patients with acute liver failure.

Reduced functional connectivity of physiological systems is associated with poor prognosis in critically ill patients. However, physiological network analysis is not commonly used in clinical practice and awaits quantitative evidence. Acute liver failure (ALF) is associated with multiorgan failure and mortality. Prognostication in ALF is highly important for clinical management but is currently dependent on models that do not consider the interaction between organ systems. This study aims to examine whether physiological network analysis can predict survival in patients with ALF. Data from 640 adult patients admitted to the ICU for paracetamol-induced ALF were extracted from the MIMIC-III database. Parenclitic network analysis was performed on the routine biomarkers using 28-day survivors as reference population and network clusters were identified for survivors and non-survivors using k-clique percolation method. Network analysis showed that liver function biomarkers were more clustered in survivors than in non-survivors. Arterial pH was also found to cluster with serum creatinine and bicarbonate in survivors compared with non-survivors, where it clustered with respiratory nodes indicating physiologically distinctive compensatory mechanism. Deviation along the pH-bicarbonate and pH-creatinine axes significantly predicts mortality independent of current prognostic indicators. These results demonstrate that network analysis can provide pathophysiologic insight and predict survival in critically ill patients with ALF.

Data-driven and equation-free methods for neurological disorders: analysis and control of the striatum network.

The striatum as part of the basal ganglia is central to both motor, and cognitive functions. Here, we propose a large-scale biophysical network for this part of the brain, using modified Hodgkin-Huxley dynamics to model neurons, and a connectivity informed by a detailed human atlas. The model shows different spatio-temporal activity patterns corresponding to lower (presumably normal) and increased cortico-striatal activation (as found in, e.g., obsessive-compulsive disorder), depending on the intensity of the cortical inputs. By applying equation-free methods, we are able to perform a macroscopic network analysis directly from microscale simulations. We identify the mean synaptic activity as the macroscopic variable of the system, which shows similarity with local field potentials. The equation-free approach results in a numerical bifurcation and stability analysis of the macroscopic dynamics of the striatal network. The different macroscopic states can be assigned to normal/healthy and pathological conditions, as known from neurological disorders. Finally, guided by the equation-free bifurcation analysis, we propose a therapeutic close loop control scheme for the striatal network.