The multistable behavior of neural networks is actively being studied as a landmark of ongoing cerebral activity, reported in both functional Magnetic Resonance Imaging (fMRI) and electro- or magnetoencephalography recordings. This consists of a continuous jumping between different partially synchronized states in the absence of external stimuli. It is thought to be an important mechanism for dealing with sensory novelty and to allow for efficient coding of information in an ever-changing surrounding environment. Many advances have been made to understand how network topology, connection delays, and noise can contribute to building this dynamic. Little or no attention, however, has been paid to the difference between local chaotic and stochastic influences on the switching between different network states. Using a conductance-based neural model that can have chaotic dynamics, we showed that a network can show multistable dynamics in a certain range of global connectivity strength and under deterministic conditions. In the present work, we characterize the multistable dynamics when the networks are, in addition to chaotic, subject to ion channel stochasticity in the form of multiplicative (channel) or additive (current) noise. We calculate the Functional Connectivity Dynamics matrix by comparing the Functional Connectivity (FC) matrices that describe the pair-wise phase synchronization in a moving window fashion and performing clustering of FCs. Moderate noise can enhance the multistable behavior that is evoked by chaos, resulting in more heterogeneous synchronization patterns, while more intense noise abolishes multistability. In networks composed of nonchaotic nodes, some noise can induce multistability in an otherwise synchronized, nonchaotic network. Finally, we found the same results regardless of the multiplicative or additive nature of noise.

Peer Review History: PLOS recognizes the benefits of transparency in the peer review process; therefore, we enable the publication of all of the content of peer review and author responses alongside final, published articles. The editorial history of this article is available here:

No abstract

Chaotic dynamics has been shown in the dynamics of neurons and neural networks, in experimental data and numerical simulations. Theoretical studies have proposed an underlying role of chaos in neural systems. Nevertheless, whether chaotic neural oscillators make a significant contribution to network behaviour and whether the dynamical richness of neural networks is sensitive to the dynamics of isolated neurons, still remain open questions. We investigated synchronization transitions in heterogeneous neural networks of neurons connected by electrical coupling in a small world topology. The nodes in our model are oscillatory neurons that – when isolated – can exhibit either chaotic or non-chaotic behaviour, depending on conductance parameters. We found that the heterogeneity of firing rates and firing patterns make a greater contribution than chaos to the steepness of the synchronization transition curve. We also show that chaotic dynamics of the isolated neurons do not always make a visible difference in the transition to full synchrony. Moreover, macroscopic chaos is observed regardless of the dynamics nature of the neurons. However, performing a Functional Connectivity Dynamics analysis, we show that chaotic nodes can promote what is known as multi-stable behaviour, where the network dynamically switches between a number of different semi-synchronized, metastable states.

The structural connectivity of human brain allows the coexistence of segregated and integrated states of activity. Neuromodulatory systems facilitate the transition between these functional states and recent computational studies have shown how an interplay between the noradrenergic and cholinergic systems define these transitions. However, there is still much to be known about the interaction between the structural connectivity and the effect of neuromodulation, and to what extent the connectome facilitates dynamic transitions. In this work, we use a whole brain model, based on the Jasen and Rit equations plus a human structural connectivity matrix, to find out which structural features of the human connectome network define the optimal neuromodulatory effects. We simulated the effect of the noradrenergic system as changes in filter gain, and studied its effects related to the global-, local-, and meso-scale features of the connectome. At the global-scale, we found that the ability of the network of transiting through a variety of dynamical states is disrupted by randomization of the connection weights. By simulating neuromodulation of partial subsets of nodes, we found that transitions between integrated and segregated states are more easily achieved when targeting nodes with greater connection strengths—local feature—or belonging to the rich club—meso-scale feature. Overall, our findings clarify how the network spatial features, at different levels, interact with neuromodulation to facilitate the switching between segregated and integrated brain states and to sustain a richer brain dynamics.