- 1Department of Psychology, University of Virginia, Charlottesville, VA, United States
- 2Department of Computer Science, University of Virginia, Charlottesville, VA, United States
- 3School of Data Science, University of Virginia, Charlottesville, VA, United States
- 4Biomedical Engineering, University of Virginia, Charlottesville, VA, United States
- 5Department of Physics and Astronomy and the Interdisciplinary Program in Neuroscience, George Mason University, Fairfax, VA, United States
Transient synchronization of bursting activity in neuronal networks, which occurs in patterns of metastable itinerant phase relationships between neurons, is a notable feature of network dynamics observed in vivo. However, the mechanisms that contribute to this dynamical complexity in neuronal circuits are not well understood. Local circuits in cortical regions consist of populations of neurons with diverse intrinsic oscillatory features. In this study, we numerically show that the phenomenon of transient synchronization, also referred to as metastability, can emerge in an inhibitory neuronal population when the neurons’ intrinsic fast-spiking dynamics are appropriately modulated by slower inputs from an excitatory neuronal population. Using a compact model of a mesoscopic-scale network consisting of excitatory pyramidal and inhibitory fast-spiking neurons, our work demonstrates a relationship between the frequency of pyramidal population oscillations and the features of emergent metastability in the inhibitory population. In addition, we introduce a method to characterize collective transitions in metastable networks. Finally, we discuss potential applications of this study in mechanistically understanding cortical network dynamics.
1 Introduction
Computational modeling of neuronal dynamics which incorporates mesoscopic-level connectivity features can potentially offer powerful frameworks for understanding neuronal network mechanisms (Markram et al., 2015; Kopsick et al., 2022). Despite the extensive theoretical work on the emergent dynamics in neuronal networks (Amit, 1989; van Vreeswijk and Sompolinsky, 1996; Rolls and Treves, 1997; Brunel, 2000; Gerstner and Kistler, 2002; Gerstner et al., 2014; Tsuda, 2015), a scalable network representation of key neurodynamical features at the mesoscopic level is currently lacking. A simplified model describing the complexities of between- and within-group neuronal interactions is necessary to advance our multiscale understanding of neuronal systems in typical and atypical conditions (Venkadesh and Van Horn, 2021).
An experimentally observable complexity of in vivo networks is the metastable attractor dynamics at the level of individual neurons. More specifically, an individual neuron can exhibit periodic bursting oscillations that transiently synchronize with the oscillations of other neurons, and these synchronizations can occur at different phases (Skaggs and McNaughton, 1996; Poe et al., 2000; Jones and Wilson, 2005; Lisman, 2005; Ego-Stengel and Wilson, 2007; Vinck et al., 2010). In other words, biological neuronal networks can realize a range of attractor states that are locally observable through different phase-locking modes between individual neurons. Such transient attractor dynamics are also referred to as metastability or itinerancy (Venkadesh et al., 2020), and this phenomenon is hypothesized to be a necessary physical property underlying the coordinated dynamics between spatially distant neuronal populations (Tognoli and Kelso, 2014).
The mathematical conditions for the emergence of transient dynamics have been previously discussed (Tsuda, 2009). More recently, conduction delays between brain areas were associated with metastable oscillations on a macroscopic spatial scale (Cabral et al., 2022). Many possible mechanisms of metastability have been proposed (Kelso, 1996; Shanahan, 2008; Brinkman et al., 2022; Hancock et al., 2023; Mackay et al., 2023). However, the biophysical features of neuronal networks underlying metastability at the level of neurons have not been fully elucidated. While a previous work reported metastable attractor dynamics among intrinsically bursting neurons (Venkadesh et al., 2020), it was unclear if metastability could also be realized in neuronal populations in which bursting dynamics are induced externally. Here, we explore the emergence and characteristics of metastability in a fast-spiking interneuron population that exhibits bursting dynamics due to driving inputs from a population of slower pyramidal neurons. We also examine the quantitative characteristics of collective transitions between attractor states realized in different metastable networks. Finally, we discuss potential applications of our work in understanding network physiology in neurodegenerative disorders.
2 Methods
Networks consisting of 100 pyramidal neurons and 50 fast-spiking interneurons were constructed using Izhikevich neurons (Izhikevich, 2003), where each neuron was governed by Eqs 1–3. The parameters
Network connections were configured as illustrated in Figure 1A. Inter- and intra-population connections were realized through instantiations of our networks using pairwise connection probabilities between individual neurons. We studied the network behavior for a complete range of connection probabilities. Post-synaptic effects were specified using an instantaneous pulse-coupling scheme (Eq. 4), where each presynaptic spike either increased (excitatory) or decreased (inhibitory) the postsynaptic current by the weight parameter (+0.3 for excitatory connections and −0.3 for inhibitory connections).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1302499/fnetp-04-1302499-HTML/image_m/fnetp-04-1302499-g001.jpg)
FIGURE 1. Endogenous itinerancy in fast-spiking interneurons. (A) Left: Network configuration consisting of pyramidal neurons (blue triangles) and fast-spiking interneurons (dark orange ovals). Connections indicated by the arrowhead and circles denote excitatory and inhibitory pulse-coupling, respectively. Right: Spiking frequencies of an isolated excitatory (pyramidal) neuron and an inhibitory (fast-spiking) neuron for various input currents. (B) Raster plot showing spiking activity in pyramidal (blue) and fast-spiking neurons (orange) for one second. (C) (top) Bursting activity (gray) in one representative fast-spiking neuron. The filtered signal, superimposed in black, captures the burst oscillations. (bottom) Burst oscillations of two representative fast-spiking neurons (black and red) showing transitions between in- and anti-phase transiently locked modes.
In Eq. 4,
2.1 Characterizing phase relationships between bursting oscillations
The basic phenomenon of interest is illustrated in Figure 1C bottom, which shows the bursting oscillations of a representative pair of neurons versus time. At first, the neurons are locked in phase, but then they transition to an out-of-phase state. After a short time, they return to the in-phase state, and then switch back to the out-of-phase state. Similar persistent behaviors were described in previous reports (Kasatkin et al., 2019; Venkadesh et al., 2020).
At the population level, the bursting phases tend to cluster and exhibit itinerant switching between clusters. This is illustrated in Figure 2. In the top panel of Figure 2A, the dots represent the instantaneous phases of all the inhibitory neurons at time
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1302499/fnetp-04-1302499-HTML/image_m/fnetp-04-1302499-g002.jpg)
FIGURE 2. (A) Endogenous itinerant states characterized by phase clusters at three time points over a 2 s duration. Neurons (black dots) are ordered by their instantaneous phases in the top panel, and this ordering is preserved in the bottom two panels. Two phase clusters persist over time, and individual neurons spontaneously switch from one cluster to the other. (B) Persistent phase relationships between neurons are illustrated in a polar histogram by plotting the phase differences between all distinct pairs of neurons over the 2 s duration.
In general, the instantaneous phases of the cluster states drift over time. In order to ignore or cancel out this drift, we calculated the phase differences between all possible neuron pairs over all values of discrete time. Plotting these phase differences on a polar histogram as in Figure 2B, we see phase difference clusters at zero and 180°. This is essentially a shift to a reference frame which gives a convenient and consistent way to visualize the phase clusters shown in Figure 2A.
More precisely, we assessed the nature of population-wide phase relationships among interneurons by calculating instantaneous phase differences between randomly selected neuron pairs over the entire duration (10 s). Henceforth, these phase differences will be referred to as
These order parameters characterize clustering patterns in the ensemble as follows. If each ensemble member is plotted as a point on the unit circle in the complex plane at its angle,
2.2 Characterizing exogenously induced transitions
To assess the nature of transitions between clustered states in different metastable networks, the following steps were carried out: A given network was simulated twice with identical initializations. The first simulation was used to establish the baseline instantaneous clustering of neurons, based on their instantaneous phases, at any given time. In the second simulation, a small group of randomly selected interneurons were simultaneously perturbed by single spikes at a certain time. Perturbations were applied via a single spike from an external pyramidal neuron. More specifically, each perturbed interneuron received a single spike that increased its voltage level by a value specified by the perturbation weight. We studied the effect of such perturbations for various numbers of perturbed neurons and for a range of perturbation weights. The perturbations caused some of the neurons to subsequently change their cluster assignments relative to the baseline.
Phase cluster assignments of neurons were obtained using Gaussian mixture models (GMM) (Pedregosa et al., 2011). The GMM group the instantaneous phases of neurons into
To quantify the effect of the perturbation, we calculated the cluster assignment of neurons at a given time in the baseline simulation, and compared this arrangement with how the neurons were distributed among the clusters at the corresponding time in the second simulation (with the perturbation). We measured the similarity of these two clustering arrangements of the neurons using Eqs 7–9.
First, we calculated the well-known “Rand Index” in Eq. 7 (Rand, 1971), where g and ℎ are the numbers of pairs of neurons with the same and different cluster assignments across the two groupings that are being compared, and
We then defined the population’s collective response CR as in Eq. 9 to characterize the extent to which the network (as a whole) generically responds to exogeneous perturbations of a fraction of its elements. The ARIs were measured in 50 simulation trials, each with different realizations of the random connections (
3 Results
The network model incorporates a significant difference in the spiking time scales of the pyramidal and fast-spiking interneurons. The slower excitatory drive from the pyramidal neurons modulates the fast-spiking dynamics of the interneurons, resulting in bursting oscillations in the latter (Figures 1B, C). The inhibitory interactions between these bursting neurons then lead to itinerant metastable behavior. That is, ongoing endogenous transitions occur between different phase-locked (attractor) states. (Figure 1C, also see Section 2.1).
3.1 Interneuron population phase clusters are modulated by the excitatory frequency
The frequency of pyramidal neuron spiking is determined by the value of
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1302499/fnetp-04-1302499-HTML/image_m/fnetp-04-1302499-g003.jpg)
FIGURE 3. The network-wide stability and the number of phase clusters in the fast-spiking interneurons as functions of pyramidal neuron frequency for one representative network. Polar histograms show phase differences
These results show that metastable attractors emerge in the interneuron population when it is appropriately driven by the pyramidal population, and that the nature of these metastable attractors is modulated by the driving frequency. It is interesting to note that the emergence of these clustered states occurs for driving frequencies within the gamma and the high-gamma ranges. This behavior was also present in a larger population of interneurons (see Supplementary Figure S1). Pyramidal neurons have indeed been observed to fire at these frequencies (O’Keefe and Recce, 1993; Dragoi and Buzsáki, 2006; Ray et al., 2008), although such gamma cycles have been observed to be modulated by slower theta ones. Additionally, although it is typical for neurons to generally exhibit aperiodic firing in vivo, selectively recruited pyramidal neurons are known to fire bursts of action potentials in a periodic manner and the phases of these bursts are strongly correlated with an organism’s behavior (O’Keefe and Recce, 1993; Dragoi and Buzsáki, 2006).
The phase-clustering phenomena described above were consistently observed in networks with different connectivity profiles realized using a range of connection probabilities (Figure 4). We analyzed networks instantiated using 50 × 50 values of equally spaced inter- and intra-group connection probabilities in the range [0, 1]. The
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1302499/fnetp-04-1302499-HTML/image_m/fnetp-04-1302499-g004.jpg)
FIGURE 4. The stability of phase clusters versus network connectivity and pyramidal neuron frequency. (A) Stabilities of the maximally stable two-cluster (left) and three-cluster (right) states realized in networks for various excitatory-inhibitory (E-I) and inhibitory-inhibitory (I-I) connection probabilities. (B) Distributions of maximum phase cluster stability are plotted against their respective frequency bins of pyramidal neurons for two-cluster (left) and three-cluster (right) states. Width indicates the density of points for a given maximum stability.
Although the
3.2 Collective response characteristics of metastable networks
The multi-cluster states described above exhibit spontaneous and endogenous itinerant dynamics (see Figures 1C, 2). Motivated by the fact that an organism adapts its internal state to the environment via external stimuli, we sought to measure how a metastable state would be affected by an external perturbation. As described in Methods Section 2.2, the networks were simulated twice with identical initializations. The first simulation established the baseline dynamics as shown in Figure 5A (left column). The panels on the left show the instantaneous phases of individual neurons in the interneuron population at three instants of time, with the neurons ordered along the horizontal axis by increasing phase. Note that the neurons were reordered in each left panel in Figure 5A in order to mask the endogenous itinerant behavior (see Figure 2 for endogenous itinerancy). In the second simulation, an external perturbation was applied. Five randomly selected interneurons each received a single spike at time
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1302499/fnetp-04-1302499-HTML/image_m/fnetp-04-1302499-g005.jpg)
FIGURE 5. Exogenous transitions in the fast-spiking interneuron population. (A) Left and right panels show the temporal evolution of a network in unperturbed and perturbed conditions respectively. Neurons (black dots) in the left panels are sorted by their instantaneous phases at each
We examined the magnitude of collective responses (CR, Eq. 9) in three representative metastable networks that showed high magnitudes of stability in Figure 3: the networks of the two-cluster state that emerged near f = 48 Hz, the three-cluster state that emerged near f = 77Hz, and the two-cluster state that emerged near f = 117 Hz. For each of these three networks, we varied the number of perturbed neurons from 1 to 20 to characterize their collective responses. Additionally, we examined the effect of the perturbation weight (see Section 2.2) by varying it from 0.1 to 1.0, while the number of perturbed neurons was set to a constant value of 5.
Figure 6 shows the CR magnitudes in the three representative networks mentioned above. When only a single neuron was perturbed, the CR remained close to zero, although the perturbed networks began to show small deviations from their respective unperturbed networks at
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1302499/fnetp-04-1302499-HTML/image_m/fnetp-04-1302499-g006.jpg)
FIGURE 6. Magnitude of collective response (Eq. 9) characterized by the exogenous transitions in three representative networks for various numbers of perturbed neurons when the perturbation weight is 0.3 (left panels) and for various perturbation weights in the range [0.1, 1.0], when the number of perturbed neurons is 5 (right panels). In the right panels, solid lines and shaded regions represent means and standard deviations, respectively, across perturbation weights. The pyramidal driving frequencies used are 48 Hz (A), 77 Hz (B), and 117 Hz (C); these correspond to two, three, and two cluster states in the interneuron population, respectively (see Figure 3).
4 Discussion
Understanding the emergent dynamical complexities in a network in terms of its constituent mechanisms remains an important goal. Here, we studied the emergent phenomenon of itinerancy, an instance of metastability, in spiking neuronal network models that represent cortical local circuits. Metastability is characterized by the coexistence of integrated (synchronized states) and segregated (transitions between synchronized states) behaviors in a system of interacting elements (Tognoli and Kelso, 2014). In this study, we characterized the extent to which neurons are synchronized using the Kuramoto-Daido order parameters for various numbers of phase clusters (Eqs 5, 6). We examined the extent to which neurons collectively transition between different synchronized states using a novel metric (Eq. 9). We wish to note that many metrics of metastability in the literature aim to quantify the extent to which both synchronized and desynchronized states coexist in a system of interacting elements. While these metrics provide useful signatures of metastability, they do not capture the synergistic nature of transitions in metastable systems. In this study, once the number of phase clusters in a network was identified using
In a previous study (Venkadesh et al., 2020), we showed that networks of complex periodic oscillators such as intrinsically bursting neurons can realize chaotic itinerancy (Tsuda, 2015). Chaotic itinerancy is an instance of metastability that is realized by the transitions between different basins of attraction through heteroclinic orbits (Hancock et al., 2023). In our current work, we studied how the intrinsic biophysical features of a spiking neuronal network are quantitatively related to metastability using a generalized model of cortical intraregional interactions. We specifically studied how the neuron-level spiking frequencies are associated with the emergent metastability, and how these relationships are governed by the underlying excitatory and inhibitory neuronal connectivity.
We showed that (1) driving input from the pyramidal population to the fast-spiking neurons caused bursting in the latter (see (Pinsky and Rinzel, 1994) for an investigation of the mechanisms of bursting), and that this led (due to inhibitory coupling) to itinerancy in the bursting phase relationships in the fast-spiking population; (2) the complexity of the itinerancy (in the form of the number of clusters) was modulated by the pyramidal driving frequency; (3) this behavior was robust with respect to the network connection probabilities as long as the E-I probability was large enough; (4) the different itinerant states had different responsiveness to perturbations.
It is worth comparing our current model and observations with those of our previous study (Venkadesh et al., 2020). Our previous model of itinerant complexity, which was constructed using only an inhibitory population of intrinsically bursting neurons is simpler than our current model, which includes both an excitatory and an inhibitory population. The construction of our current model was motivated by several considerations: First, intrinsically bursting neurons are biologically rare. For instance, less than 15% of inhibitory neurons are intrinsically “bursting” (multi-periodic) in the rodent hippocampus, and the remaining neurons are intrinsically “spiking” (single-periodic) (Komendantov et al., 2019). Our current model approximation includes slower excitatory and faster inhibitory neurons, both of which are intrinsically single-periodic, and this circuit configuration is ubiquitous in the cortex. Furthermore, as we showed in our previous study, the intrinsically single-periodic neural populations, by themselves, are unable to realize itinerant complexity. On the other hand, the current model describes mechanisms by which itinerant complexity could be realized in such simpler inhibitory neurons that vastly outnumber their intrinsically bursting counterparts in the cortex. Taken together, our studies support the idea that the same dynamics could emerge from multiple network configurations in the cortex (Friston and Price, 2003). Secondly, our previous model did not show that the number of phase clusters within an inhibitory population could change depending on the frequency. In fact, our previous study reported emergent itinerancy at different bursting regimes of the model, where 3-cluster states that only differed in their magnitudes of stability were observed. On the other hand, our current model shows, in addition to the magnitude of stability, the number of phase clusters can also change within an inhibitory population depending on the frequency of excitatory drive. It is worth noting that many intrinsically bursting neurons have a narrow regime for bursting (Komendantov et al., 2019). In other words, the range of excitatory inputs for which an intrinsically bursting neuron actually bursts (as opposed to spike) is narrow. On the other hand, the externally induced bursting is observable for a wide range of excitatory inputs as we have shown in our current study. Finally, our current model suggests that the inhibitory population’s itinerancy could be modulated by electrically stimulating excitatory neurons such as the pyramidal neurons. This is especially important in the context of neuromodulation (Luan et al., 2014; Krishna et al., 2018; Cole et al., 2022). Targeting pyramidal neurons for neuromodulation is currently more feasible than targeting regionally local inhibitory neurons. Because pyramidal neurons project their axons over long distances and their connectivity can be noninvasively mapped using neuroimaging modalities, they are highly feasible targets for neuromodulation. On the other hand, an inhibitory neuronal population’s connectivity is generally restricted to a cortical region, and current neuromodulation methods lack spatial specificity to selectively target such regionally local neuronal populations. Overall, while our previous model was computationally more compact than our current model, this compactness traded off its practical applicability. On the other hand, our current model, while still scalable, has broader applicability and provides a framework for future studies to make testable predictions as discussed in the final paragraph of this section.
The number of connections an inhibitory neuron receives is a crucial parameter that affects the emergence of metastable clusters. Our current model suggests that excitatory connections alone are insufficient to allow the emergence of metastable clusters in an inhibitory population (Supplementary Figure S2). Similarly, inhibitory connections alone are insufficient to allow the emergence of metastable clusters in “non-intrinsically bursting” neurons as we have shown in our previous study. Therefore, it is an appropriate balance between the excitatory and inhibitory connections that allows the emergence of metastable behaviors in our current model. We have shown that such optimal balance could be achieved for various connection probabilities (Figure 4A) and for a larger neuronal population (Supplementary Figure S1). Importantly, the frequency bins at which the 2-cluster and the 3-cluster states achieve maximal stability are consistent across various connection probabilities (Figure 4B) and for a larger neuronal population (Supplementary Figure S1). It is worth mentioning that these frequency bins are characteristic of networks with the specific neuron types considered in our model. These neuron types differ in their intrinsic spiking frequencies (Figure 1). Additionally, the diversity of cortical neuronal populations is not only revealed in their intrinsic spiking frequencies, but also in other temporal features such as delayed spiking, where a neuron elicits spikes following a brief period of quiescence for excitatory inputs. There are many neuronal populations that exhibit such delayed spiking behaviors (Komendantov et al., 2019). This is an important consideration in light of the fact that many studies associated connection delays with emergent metastability (Cabral et al., 2022; Hancock et al., 2023). Future work will investigate if and how the characteristic frequency bins reported in this study are affected by including additional inhibitory neuronal populations with distinct intrinsic spiking profiles and by incorporating a range of connection delays.
The spiking neuronal network models examined in this paper provide compact mesoscopic scale approximations of regional collective dynamics that capture a crucial biological complexity at the level of individual neurons. A limitation of the current work is that it only considered two types of neuronal populations, whereas intraregional neuronal circuits consist of populations with diverse frequency profiles (Komendantov et al., 2019). Nevertheless, the current work demonstrates that the emergent metastability in an interneuron population depends on their connectivity with excitatory neurons and the frequency of the excitatory drive. Their relationships are revealed in the
A previous study reported clustering behaviors in a network of excitatory and inhibitory neurons (Stefanescu and Jirsa, 2008). The authors investigated cluster transitions of neurons under various conditions. It is worth noting that they reported several instances of static clusters in parameter space, where cluster transitions occurred when changing the network parameters such as the connection weights. However, endogenous transitory phenomena have been previously described by Ichiro Tsuda and collaborators using the framework of chaotic itinerancy (Tsuda, 2001; Tsuda et al., 2004). We wish to note that the novelty of our work is in the (a) description of itinerancy that occurs at the level of bursting oscillations, which was motivated by the experimental observations of this phenomenon (see Section 1), (b) association of the emergence of such itinerancy with the underlying biophysical mechanisms, most notably, the intrinsic differences in the membrane excitability levels between the excitatory and inhibitory neural populations, (c) reliable modulation of the itinerant behavior in inhibitory neurons via excitatory neurons, and (d) formulation of a method to characterize the synergy of neural transitions between multiple synchronized states. Furthermore, the Hindmarsh-Rose bursting neurons employed in (Stefanescu and Jirsa, 2008) did not show multiple phase clusters for a single network configuration. The phases of individual bursting neurons were uniformly scattered around a mean phase. The absence of multi-clustered states in their bursting neurons may be because they neglected the connections within the inhibitory subpopulation. In the present model, the bursting neurons not only show multiple distinct phase clusters for a single network configuration, but they also exhibit endogenous cluster transitions without changing the network parameters. It may also be worth pointing out the differences between their network model and the one presented in this study: (1) In (Stefanescu and Jirsa, 2008), the excitatory and the inhibitory neural subpopulations did not differ in their intrinsic dynamics. Their model assumed that in a small cortical volume, the factors affecting the membrane excitability will be similar for all neurons. However, recent efforts to comprehensively characterize the intrinsic physiology of neurons revealed enormous diversity in neural subpopulations even within a small cortical volume (Markram et al., 2015; Komendantov et al., 2019). Motivated by this fact, our model included excitatory and inhibitory subpopulations with sharply different intrinsic excitability levels, which are revealed in their current-frequency curves (Figure 1A). (2) In (Stefanescu and Jirsa, 2008), the connections within the inhibitory subpopulation were neglected. However, connectivity within the inhibitory subpopulations exist in the cortex (Wheeler et al., 2024), and it is an important characterizing feature in data-driven simulations of cortical circuits (Kopsick et al., 2022). As we have shown in this study, connections within the inhibitory subpopulation are crucial for the emergence of itinerant dynamics. (3) Finally, their model utilized a global coupling scheme for the network architecture, where every neuron is connected to every other neuron. In our model, we probabilistically realized connections between neurons, and we explored the network dynamics for the full range of connection probabilities. It is worth noting that cortical neurons are not globally connected, and connection probabilities are important parameters in data-driven simulations of cortical circuits (Markram et al., 2015; Kopsick et al., 2022).
Additionally, the models and the methods of analysis presented in this paper are easily scalable to study interregional interactions in the brain. Empirically estimated anatomical connectivity between brain regions provides useful information for modeling interregional metastability. Note that while the current work suggests a link between two temporal scales (i.e., timescales of neuronal spiking and attractor transitions) of network dynamics, modeling interregional interactions using their anatomical connectivity will further provide a venue to study metastable brain dynamics at a macroscopic spatial scale. Such models are useful mathematical tools that will enable a multiscale understanding of cortical dynamics and illuminate neuronal network mechanisms of pathophysiology in neurological disorders. As an example, we briefly discuss how a brain-wide analysis of metastable dynamics could be useful in elucidating the pathophysiology of Parkinson’s disease (PD). The structural degeneration of cortical and subcortical connectivity is a key factor underlying the motor and cognitive deficits in PD (Cochrane and Ebmeier, 2013; Atkinson-Clement et al., 2017; McGregor and Nelson, 2019). However, the mechanisms by which the structural degeneration alters the brain-wide dynamics in PD are not well understood. A whole-brain model of metastable dynamics constrained by the PD-affected anatomical connectivity can potentially delineate such mechanisms. For instance, amplified synchronization in beta frequency (8 Hz–30 Hz), which has been reported in the cortical and basal ganglia circuits of PD subjects, is hypothesized to contribute to PD symptoms (Hammond et al., 2007; Brittain et al., 2014; McGregor and Nelson, 2019). In particular, longer episodes of beta oscillations in the subthalamic nucleus (STN) (Deffains et al., 2018) and its higher synchronization with the cortical supplementary motor area (SMA) were associated with the freezing of gait in PD (Toledo et al., 2014). It may be hypothesized that the freezing of gait in PD is contributed by a reduction in synergistic transitions between different attractor states realized in these areas. The increased beta synchronization between STN and SMA observed in PD, and the hypothesized reduction in attractor transitions can be validated in whole-brain models by examining their stability (Eq. 6) and collective response characteristics (Eq. 9) respectively. Furthermore, by selectively targeting frequency bands for manipulation, one can investigate causal relationships between specific frequencies of oscillations and the dynamics of transitions between attractor states realized across multiple timescales. Such models could also enable the identification of network mechanisms for restoring optimal brain dynamics in PD and other neurological disorders via neuromodulation.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://github.com/sivaven/Transient-Synchronization.
Author contributions
SV: Conceptualization, Formal Analysis, Investigation, Methodology, Visualization, Writing–original draft. AS: Formal Analysis, Visualization, Writing–review and editing. HS: Investigation, Writing–review and editing. EB: Conceptualization, Funding acquisition, Investigation, Supervision, Visualization, Writing–review and editing. JV: Conceptualization, Funding acquisition, Investigation, Resources, Supervision, Writing–review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported in parts by the Department of Defense grant W81XWH 18 1 0665 PD170037 and an internal funding from the University of Virginia.
Acknowledgments
The authors thank Alev Erisir, Teague Henry, Giselle Petzinger, Michael Jakowec, Dawn Schiehser and Joseph O’Neill for their support, insights, and helpful discussions. The authors appreciate the complementary perspectives of Giorgio Ascoli and Alexander Komendantov and are thankful for their many valuable discussions.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnetp.2024.1302499/full#supplementary-material
References
Amit, D. J. (1989). Modeling brain function: the world of attractor neural networks. Cambridge, United Kingdom: Cambridge University Press.
Atkinson-Clement, C., Pinto, S., Eusebio, A., and Coulon, O. (2017). Diffusion tensor imaging in Parkinson’s disease: review and meta-analysis. Neuroimage Clin. 16, 98–110. doi:10.1016/j.nicl.2017.07.011
Brinkman, B. A. W., Yan, H., Maffei, A., Park, I. M., Fontanini, A., Wang, J., et al. (2022). Metastable dynamics of neural circuits and networks. Appl. Phys. Rev. 9, 011313. doi:10.1063/5.0062603
Brittain, J.-S., Sharott, A., and Brown, P. (2014). The highs and lows of beta activity in cortico-basal ganglia loops. Eur. J. Neurosci. 39, 1951–1959. doi:10.1111/ejn.12574
Brunel, N. (2000). Dynamics of sparsely connected networks of excitatory and inhibitory spiking neurons. J. Comput. Neurosci. 8, 183–208. doi:10.1023/a:1008925309027
Cabral, J., Castaldo, F., Vohryzek, J., Litvak, V., Bick, C., Lambiotte, R., et al. (2022). Metastable oscillatory modes emerge from synchronization in the brain spacetime connectome. Commun. Phys. 5, 184–213. doi:10.1038/s42005-022-00950-y
Cochrane, C. J., and Ebmeier, K. P. (2013). Diffusion tensor imaging in parkinsonian syndromes: a systematic review and meta-analysis. Neurology 80, 857–864. doi:10.1212/WNL.0b013e318284070c
Cole, R. C., Okine, D. N., Yeager, B. E., and Narayanan, N. S. (2022). Neuromodulation of cognition in Parkinson’s disease. Prog. Brain Res. 269, 435–455. doi:10.1016/bs.pbr.2022.01.016
Daido, H. (1996). Onset of cooperative entrainment in limit-cycle oscillators with uniform all-to-all interactions: bifurcation of the order function. Phys. D. Nonlinear Phenom. 91, 24–66. doi:10.1016/0167-2789(95)00260-x
Daoud, M. S., Shehab, M., Al-Mimi, H. M., Abualigah, L., Zitar, R. A., and Shambour, M. K. Y. (2023). Gradient-based optimizer (gbo): a review, theory, variants, and applications. Arch. Comput. Methods Eng. 30, 2431–2449. doi:10.1007/s11831-022-09872-y
Deffains, M., Iskhakova, L., Katabi, S., Israel, Z., and Bergman, H. (2018). Longer β oscillatory episodes reliably identify pathological subthalamic activity in Parkinsonism. Mov. Disord. 33, 1609–1618. doi:10.1002/mds.27418
Dragoi, G., and Buzsáki, G. (2006). Temporal encoding of place sequences by hippocampal cell assemblies. Neuron 50, 145–157. doi:10.1016/j.neuron.2006.02.023
Ego-Stengel, V., and Wilson, M. A. (2007). Spatial selectivity and theta phase precession in CA1 interneurons. Hippocampus 17, 161–174. doi:10.1002/hipo.20253
Friston, K. J., and Price, C. J. (2003). Degeneracy and redundancy in cognitive anatomy. Trends Cogn. Sci. 7, 151–152. doi:10.1016/s1364-6613(03)00054-8
Gerstner, W., and Kistler, W. M. (2002). Spiking neuron models: single neurons, populations, plasticity. Available at: https://books.google.ca/books?hl=en&lr=&id=Rs4oc7HfxIUC&oi=fnd&pg=PR11&ots=2Sh-xSeJT8&sig=x1pk7wNYaCFQ9VbQhrhC_CNeK3E.
Gerstner, W., Kistler, W. M., Naud, R., and Paninski, L. (2014). Neuronal dynamics: from single neurons to networks and models of cognition. Available at: https://books.google.ca/books?hl=en&lr=&id=D4j2AwAAQBAJ&oi=fnd&pg=PR9&ots=-E3ns_xQ-5&sig=tYPm8nWtMgJE-1ZkEp9TQyS4XSs.
Golomb, D., and Hansel, D. (2000). The number of synaptic inputs and the synchrony of large, sparse neuronal networks. Neural comput., 12, 1095–1139. doi:10.1162/089976600300015529
Hammond, C., Bergman, H., and Brown, P. (2007). Pathological synchronization in Parkinson’s disease: networks, models and treatments. Trends Neurosci. 30, 357–364. doi:10.1016/j.tins.2007.05.004
Hancock, F., Rosas, F. E., Zhang, M., Mediano, P. A. M., Luppi, A., Cabral, J., et al. (2023). Metastability demystified — the foundational past, the pragmatic present, and the potential future. Preprints. doi:10.20944/preprints202307.1445.v1
Hubert, L., and Arabie, P. (1985). Comparing partitions. J. Classif. 2, 193–218. doi:10.1007/bf01908075
Izhikevich, E. M. (2003). Simple model of spiking neurons. IEEE Trans. Neural Netw. 14, 1569–1572. doi:10.1109/TNN.2003.820440
Jones, M. W., and Wilson, M. A. (2005). Phase precession of medial prefrontal cortical activity relative to the hippocampal theta rhythm. Hippocampus 15, 867–873. doi:10.1002/hipo.20119
Kasatkin, D. V., Klinshov, V. V., and Nekorkin, V. I. (2019). Itinerant chimeras in an adaptive network of pulse-coupled oscillators. Phys. Rev. E 99, 022203. doi:10.1103/PhysRevE.99.022203
Kelso, J. A. S. (1996). Dynamic patterns: the self-organization of brain and behavior. J. Cogn. Neurosci. 8, 385–386. doi:10.1162/jocn.1996.8.4.385
Komendantov, A. O., Venkadesh, S., Rees, C. L., Wheeler, D. W., Hamilton, D. J., and Ascoli, G. A. (2019). Quantitative firing pattern phenotyping of hippocampal neuron types. Sci. Rep. 9, 17915. doi:10.1038/s41598-019-52611-w
Kopsick, J. D., Tecuatl, C., Moradi, K., Attili, S. M., Kashyap, H. J., Xing, J., et al. (2022). Robust resting-state dynamics in a large-scale spiking neural network model of area CA3 in the mouse Hippocampus. Cogn. Comput. 15, 1190–1210. doi:10.1007/s12559-021-09954-2
Krishna, V., Sammartino, F., and Rezai, A. (2018). A review of the current therapies, challenges, and future directions of transcranial focused ultrasound technology: advances in diagnosis and treatment. JAMA Neurol. 75, 246–254. doi:10.1001/jamaneurol.2017.3129
Lisman, J. (2005). The theta/gamma discrete phase code occurring during the hippocampal phase precession may be a more general brain coding scheme. Hippocampus 15, 913–922. doi:10.1002/hipo.20121
Luan, S., Williams, I., Nikolic, K., and Constandinou, T. G. (2014). Neuromodulation: present and emerging methods. Front. Neuroeng. 7, 27. doi:10.3389/fneng.2014.00027
Mackay, M., Huo, S., and Kaiser, M. (2023). Spatial organisation of the mesoscale connectome: a feature influencing synchrony and metastability of network dynamics. PLoS Comput. Biol. 19, e1011349. doi:10.1371/journal.pcbi.1011349
Markram, H., Muller, E., Ramaswamy, S., Reimann, M. W., Abdellah, M., Sanchez, C. A., et al. (2015). Reconstruction and simulation of neocortical microcircuitry. Cell 163, 456–492. doi:10.1016/j.cell.2015.09.029
McGregor, M. M., and Nelson, A. B. (2019). Circuit mechanisms of Parkinson’s disease. Neuron 101, 1042–1056. doi:10.1016/j.neuron.2019.03.004
O’Keefe, J., and Recce, M. L. (1993). Phase relationship between hippocampal place units and the EEG theta rhythm. Hippocampus 3, 317–330. doi:10.1002/hipo.450030307
Ott, E., Sommerer, J. C., Alexander, J. C., and KanYorke, I. J. A. (1993). Scaling behavior of chaotic systems with riddled basins. Phys. Rev. Lett. 71, 4134–4137. doi:10.1103/PhysRevLett.71.4134
Pedregosa, F., Varoquaux, G., Gramfort, A., Michel, V., Thirion, B., Grisel, O., et al. (2011). Scikit-learn: machine learning in Python. J. Mach. Learn. Res. 12, 2825–2830. doi:10.5555/1953048.2078195
Pinsky, P. F., and Rinzel, J. (1994). Intrinsic and network rhythmogenesis in a reduced Traub model for CA3 neurons. J. Comput. Neurosci. 1, 39–60. doi:10.1007/BF00962717
Poe, G. R., Nitz, D. A., McNaughton, B. L., and Barnes, C. A. (2000). Experience-dependent phase-reversal of hippocampal neuron firing during REM sleep. Brain Res. 855, 176–180. doi:10.1016/s0006-8993(99)02310-0
Rand, W. M. (1971). Objective criteria for the evaluation of clustering methods. J. Am. Stat. Assoc. 66, 846–850. doi:10.1080/01621459.1971.10482356
Ray, S., Crone, N. E., Niebur, E., Franaszczuk, P. J., and Hsiao, S. S. (2008). Neural correlates of high-gamma oscillations (60-200 Hz) in macaque local field potentials and their potential implications in electrocorticography. J. Neurosci. 28, 11526–11536. doi:10.1523/JNEUROSCI.2848-08.2008
Rolls, E., and Treves, A. (1997). Neural networks and brain function. London, England: Oxford University Press.
Selesnick, I. W., and Burrus, C. S. (1998). Generalized digital Butterworth filter design. IEEE Trans. Signal Process. 46, 1688–1694. doi:10.1109/78.678493
Shanahan, M. (2008). Dynamical complexity in small-world networks of spiking neurons. Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 78, 041924. doi:10.1103/PhysRevE.78.041924
Skaggs, W. E., and McNaughton, B. L. (1996). Replay of neuronal firing sequences in rat hippocampus during sleep following spatial experience. Science 271, 1870–1873. doi:10.1126/science.271.5257.1870
Sommerer, J. C., and Ott, E. (1993). A physical system with qualitatively uncertain dynamics. Nature 365, 138–140. doi:10.1038/365138a0
Stefanescu, R. A., and Jirsa, V. K. (2008). A low dimensional description of globally coupled heterogeneous neural networks of excitatory and inhibitory neurons. PLoS Comput. Biol. 4, e1000219. doi:10.1371/journal.pcbi.1000219
Steinley, D. (2004). Properties of the hubert-arabie adjusted Rand index. Psychol. Methods 9, 386–396. doi:10.1037/1082-989X.9.3.386
Stimberg, M., Brette, R., and Goodman, D. F. (2019). Brian 2, an intuitive and efficient neural simulator. Elife 8, e47314. doi:10.7554/eLife.47314
Tognoli, E., and Kelso, J. A. S. (2014). The metastable brain. Neuron 81, 35–48. doi:10.1016/j.neuron.2013.12.022
Toledo, J. B., López-Azcárate, J., Garcia-Garcia, D., Guridi, J., Valencia, M., Artieda, J., et al. (2014). High beta activity in the subthalamic nucleus and freezing of gait in Parkinson’s disease. Neurobiol. Dis. 64, 60–65. doi:10.1016/j.nbd.2013.12.005
Tsuda, I. (2001). Toward an interpretation of dynamic neural activity in terms of chaotic dynamical systems. Behav. Brain Sci. 24, 793–810. ; discussion 810-48. doi:10.1017/s0140525x01000097
Tsuda, I. (2009). Hypotheses on the functional roles of chaotic transitory dynamics. Chaos 19, 015113. doi:10.1063/1.3076393
Tsuda, I. (2015). Chaotic itinerancy and its roles in cognitive neurodynamics. Curr. Opin. Neurobiol. 31, 67–71. doi:10.1016/j.conb.2014.08.011
Tsuda, I., Fujii, H., Tadokoro, S., Yasuoka, T., and Yamaguti, Y. (2004). Chaotic itinerancy as a mechanism of irregular changes between synchronization and desynchronization in a neural network. J. Integr. Neurosci. 3, 159–182. doi:10.1142/s021963520400049x
van Vreeswijk, C., and Sompolinsky, H. (1996). Chaos in neuronal networks with balanced excitatory and inhibitory activity. Science 274, 1724–1726. doi:10.1126/science.274.5293.1724
Venkadesh, S., Barreto, E., and Ascoli, G. A. (2020). Itinerant complexity in networks of intrinsically bursting neurons. Chaos 30, 061106. doi:10.1063/5.0010334
Venkadesh, S., and Van Horn, J. D. (2021). Integrative models of brain structure and dynamics: concepts, challenges, and methods. Front. Neurosci. 15, 752332. doi:10.3389/fnins.2021.752332
Vinck, M., Lima, B., Womelsdorf, T., Oostenveld, R., Singer, W., Neuenschwander, S., et al. (2010). Gamma-phase shifting in awake monkey visual cortex. J. Neurosci. 30, 1250–1257. doi:10.1523/JNEUROSCI.1623-09.2010
Keywords: bursting, synchronization, metastability, emergence, complexity, networks, itinerancy
Citation: Venkadesh S, Shaikh A, Shakeri H, Barreto E and Van Horn JD (2024) Biophysical modulation and robustness of itinerant complexity in neuronal networks. Front. Netw. Physiol. 4:1302499. doi: 10.3389/fnetp.2024.1302499
Received: 26 September 2023; Accepted: 26 February 2024;
Published: 07 March 2024.
Edited by:
Alexander N. Pisarchik, Universidad Politécnica de Madrid, SpainReviewed by:
Simona Olmi, National Research Council (CNR), ItalyLuc Berthouze, University of Sussex, United Kingdom
Copyright © 2024 Venkadesh, Shaikh, Shakeri, Barreto and Van Horn. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Siva Venkadesh, siv30@pitt.edu; Ernest Barreto, ebarreto@gmu.edu
†These authors share senior authorship