The effects of cholinergic neuromodulation on neuronal phase-response curves of modeled cortical neurons

Acetylcholine facilitates localized synaptic potentiation and location specific feature binding

Forebrain acetylcholine (ACh) signaling has been shown to drive attention and learning. Recent experimental evidence of spatially and temporally constrained cholinergic signaling has sparked interest to investigate how it facilitates stimulus-induced learning. We use biophysical excitatory-inhibitory (E-I) multi-module neural network models to show that external stimuli and ACh signaling can mediate spatially constrained synaptic potentiation patterns. The effects of ACh on neural excitability are simulated by varying the conductance of a muscarinic receptor-regulated hyperpolarizing slow K+ current (m-current). Each network module consists of an E-I network with local excitatory connectivity and global inhibitory connectivity. The modules are interconnected with plastic excitatory synaptic connections, that change via a spike-timing-dependent plasticity (STDP) rule. Our results indicate that spatially constrained ACh release influences the information flow represented by network dynamics resulting in selective reorganization of inter-module interactions. Moreover the information flow depends on the level of synchrony in the network. For highly synchronous networks, the more excitable module leads firing in the less excitable one resulting in strengthening of the outgoing connections from the former and weakening of its incoming synapses. For networks with more noisy firing patterns, activity in high ACh regions is prone to induce feedback firing of synchronous volleys and thus strengthening of the incoming synapses to the more excitable region and weakening of outgoing synapses. Overall, these results suggest that spatially and directionally specific plasticity patterns, as are presumed necessary for feature binding, can be mediated by spatially constrained ACh release.

Cellular rules underlying psychedelic control of prefrontal pyramidal neurons

Classical psychedelic drugs are thought to increase excitability of pyramidal cells in prefrontal cortex via activation of serotonin 2A receptors (5-HT2ARs). Here, we instead find that multiple classes of psychedelics dose-dependently suppress intrinsic excitability of pyramidal neurons, and that extracellular delivery of psychedelics decreases excitability significantly more than intracellular delivery. A previously unknown mechanism underlies this psychedelic drug action: enhancement of ubiquitously expressed potassium "M-current" channels that is independent of 5-HT2R activation. Using machine-learning-based data assimilation models, we show that M-current activation interacts with previously described mechanisms to dramatically reduce intrinsic excitability and shorten working memory timespan. Thus, psychedelic drugs suppress intrinsic excitability by modulating ion channels that are expressed throughout the brain, potentially triggering homeostatic adjustments that can contribute to widespread therapeutic benefits.

Tuning Neural Synchronization: The Role of Variable Oscillation Frequencies in Neural Circuits

Brain oscillations emerge during sensory and cognitive processes and have been classified into different frequency bands. Yet, even within the same frequency band and between nearby brain locations, the exact frequencies of brain oscillations can differ. These frequency differences (detuning) have been largely ignored and play little role in current functional theories of brain oscillations. This contrasts with the crucial role that detuning plays in synchronization theory, as originally derived in physical systems. Here, we propose that detuning is equally important to understand synchronization in biological systems. Detuning is a critical control parameter in synchronization, which is not only important in shaping phase-locking, but also in establishing preferred phase relations between oscillators. We review recent evidence that frequency differences between brain locations are ubiquitous and essential in shaping temporal neural coordination. With the rise of powerful experimental techniques to probe brain oscillations, the contributions of exact frequency and detuning across neural circuits will become increasingly clear and will play a key part in developing a new understanding of the role of oscillations in brain function.

Modeling cortical synaptic effects of anesthesia and their cholinergic reversal

General anesthetics work through a variety of molecular mechanisms while resulting in the common end point of sedation and loss of consciousness. Generally, the administration of common anesthetics induces reduction in synaptic excitation while promoting synaptic inhibition. Exogenous modulation of the anesthetics' synaptic effects can help determine the neuronal pathways involved in anesthesia. For example, both animal and human studies have shown that exogenously induced increases in acetylcholine in the brain can elicit wakeful-like behavior despite the continued presence of the anesthetic. However, the underlying mechanisms of anesthesia reversal at the cellular level have not been investigated. Here we apply a computational model of a network of excitatory and inhibitory neurons to simulate the network-wide effects of anesthesia, due to changes in synaptic inhibition and excitation, and their reversal by cholinergic activation through muscarinic receptors. We use a differential evolution algorithm to fit model parameters to match measures of spiking activity, neuronal connectivity, and network dynamics recorded in the visual cortex of rodents during anesthesia with desflurane in vivo. We find that facilitating muscarinic receptor effects of acetylcholine on top of anesthetic-induced synaptic changes predicts the reversal of anesthetic suppression of neurons' spiking activity, functional connectivity, as well as pairwise and population interactions. Thus, our model predicts a specific neuronal mechanism for the cholinergic reversal of anesthesia consistent with experimental behavioral observations.

A dynamic clamp protocol to artificially modify cell capacitance

Dynamics of excitable cells and networks depend on the membrane time constant, set by membrane resistance and capacitance. Whereas pharmacological and genetic manipulations of ionic conductances of excitable membranes are routine in electrophysiology, experimental control over capacitance remains a challenge. Here, we present capacitance clamp, an approach that allows electrophysiologists to mimic a modified capacitance in biological neurons via an unconventional application of the dynamic clamp technique. We first demonstrate the feasibility to quantitatively modulate capacitance in a mathematical neuron model and then confirm the functionality of capacitance clamp in in vitro experiments in granule cells of rodent dentate gyrus with up to threefold virtual capacitance changes. Clamping of capacitance thus constitutes a novel technique to probe and decipher mechanisms of neuronal signaling in ways that were so far inaccessible to experimental electrophysiology.

Dynamical Mechanism Underlying Scale-Free Network Reorganization in Low Acetylcholine States Corresponding to Slow Wave Sleep

Sleep is indispensable for most animals' cognitive functions, and is hypothesized to be a major factor in memory consolidation. Although we do not fully understand the mechanisms of network reorganisation driving memory consolidation, available data suggests that sleep-associated neurochemical changes may be important for such processes. In particular, global acetylcholine levels change across the sleep/wake cycle, with high cholinergic tone during wake and REM sleep and low cholinergic tone during slow wave sleep. Furthermore, experimental perturbation of cholinergic tone has been shown to impact memory storage. Through in silico modeling of neuronal networks, we show how spiking dynamics change in highly heterogenous networks under varying levels of cholinergic tone, with neuronal networks under high cholinergic modulation firing asynchronously and at high frequencies, while those under low cholinergic modulation exhibit synchronous patterns of activity. We further examined the network's dynamics and its reorganization mediated via changing levels of acetylcholine within the context of different scale-free topologies, comparing network activity within the hub cells, a small group of neurons having high degree connectivity, and with the rest of the network. We show a dramatic, state-dependent change in information flow throughout the network, with highly active hub cells integrating information in a high-acetylcholine state, and transferring it to rest of the network in a low-acetylcholine state. This result is experimentally corroborated by frequency-dependent frequency changes observed in vivo experiments. Together, these findings provide insight into how new neurons are recruited into memory traces during sleep, a mechanism which may underlie system memory consolidation.

Acetylcholine-gated current translates wake neuronal firing rate information into a spike timing-based code in Non-REM sleep, stabilizing neural network dynamics during memory consolidation

Sleep is critical for memory consolidation, although the exact mechanisms mediating this process are unknown. Combining reduced network models and analysis of in vivo recordings, we tested the hypothesis that neuromodulatory changes in acetylcholine (ACh) levels during non-rapid eye movement (NREM) sleep mediate stabilization of network-wide firing patterns, with temporal order of neurons’ firing dependent on their mean firing rate during wake. In both reduced models and in vivo recordings from mouse hippocampus, we find that the relative order of firing among neurons during NREM sleep reflects their relative firing rates during prior wake. Our modeling results show that this remapping of wake-associated, firing frequency-based representations is based on NREM-associated changes in neuronal excitability mediated by ACh-gated potassium current. We also show that learning-dependent reordering of sequential firing during NREM sleep, together with spike timing-dependent plasticity (STDP), reconfigures neuronal firing rates across the network. This rescaling of firing rates has been reported in multiple brain circuits across periods of sleep. Our model and experimental data both suggest that this effect is amplified in neural circuits following learning. Together our data suggest that sleep may bias neural networks from firing rate-based towards phase-based information encoding to consolidate memories.

We show that neuromodulatory changes during non-rapid eye movement (NREM) sleep generate stable spike timing relationships between neurons, the ordering of which reflects the neurons’ relative firing rates during wake. Learning-dependent ordering of firing in the hippocampus during NREM, acting in tandem with spike timing-dependent plasticity, reconfigures neuronal firing rates across the hippocampal network. This “rescaling” of neuronal firing rates has recently been reported in multiple brain circuits across periods of sleep. Together, our results suggest that the brain is remapping frequency-biased representations of information formed during wake into timing biased-representations during NREM sleep.

Biophysical models of intrinsic homeostasis: Firing rates and beyond

In view of ever-changing conditions both in the external world and in intrinsic brain states, maintaining the robustness of computations poses a challenge, adequate solutions to which we are only beginning to understand. At the level of cell-intrinsic properties, biophysical models of neurons permit one to identify relevant physiological substrates that can serve as regulators of neuronal excitability and to test how feedback loops can stabilize crucial variables such as long-term calcium levels and firing rates. Mathematical theory has also revealed a rich set of complementary computational properties arising from distinct cellular dynamics and even shaping processing at the network level. Here, we provide an overview over recently explored homeostatic mechanisms derived from biophysical models and hypothesize how multiple dynamical characteristics of cells, including their intrinsic neuronal excitability classes, can be stably controlled.

Deciphering how interneuron specific 3 cells control oriens lacunosum-moleculare cells to contribute to circuit function

The wide diversity of inhibitory cells across the brain makes them suitable to contribute to network dynamics in specialized fashions. However, the contributions of a particular inhibitory cell type in a behaving animal are challenging to untangle as one needs to both record cellular activities and identify the cell type being recorded. Thus, using computational modeling and theory to predict and hypothesize cell-specific contributions is desirable. Here, we examine potential contributions of interneuron-specific 3 (I-S3) cells - an inhibitory interneuron found in CA1 hippocampus that only targets other inhibitory interneurons - during simulated theta rhythms. We use previously developed multi-compartment models of oriens lacunosum-moleculare (OLM) cells, the main target of I-S3 cells, and explore how I-S3 cell inputs during in vitro and in vivo scenarios contribute to theta. We find that I-S3 cells suppress OLM cell spiking, rather than engender its spiking via post-inhibitory rebound mechanisms, and contribute to theta frequency spike resonance during simulated in vivo scenarios. To elicit recruitment similar to in vitro experiments, inclusion of disinhibited pyramidal cell inputs is necessary, implying that I-S3 cell firing broadens the window for pyramidal cell disinhibition. Using in vivo virtual networks, we show that I-S3 cells contribute to a sharpening of OLM cell recruitment at theta frequencies. Further, shifting the timing of I-S3 cell spiking due to external modulation shifts the timing of the OLM cell firing and thus disinhibitory windows. We propose a specialized contribution of I-S3 cells to create temporally precise coordination of modulation pathways.

Theta-gamma coupling emerges from spatially heterogeneous cholinergic neuromodulation

Theta and gamma rhythms and their cross-frequency coupling play critical roles in perception, attention, learning, and memory. Available data suggest that forebrain acetylcholine (ACh) signaling promotes theta-gamma coupling, although the mechanism has not been identified. Recent evidence suggests that cholinergic signaling is both temporally and spatially constrained, in contrast to the traditional notion of slow, spatially homogeneous, and diffuse neuromodulation. Here, we find that spatially constrained cholinergic stimulation can generate theta-modulated gamma rhythms. Using biophysically-based excitatory-inhibitory (E-I) neural network models, we simulate the effects of ACh on neural excitability by varying the conductance of a muscarinic receptor-regulated K+ current. In E-I networks with local excitatory connectivity and global inhibitory connectivity, we demonstrate that theta-gamma-coupled firing patterns emerge in ACh modulated network regions. Stable gamma-modulated firing arises within regions with high ACh signaling, while theta or mixed theta-gamma activity occurs at the peripheries of these regions. High gamma activity also alternates between different high-ACh regions, at theta frequency. Our results are the first to indicate a causal role for spatially heterogenous ACh signaling in the emergence of localized theta-gamma rhythmicity. Our findings also provide novel insights into mechanisms by which ACh signaling supports the brain region-specific attentional processing of sensory information.

M-current induced Bogdanov-Takens bifurcation and switching of neuron excitability class

In this work, we consider a general conductance-based neuron model with the inclusion of the acetycholine sensitive, M-current. We study bifurcations in the parameter space consisting of the applied current [Formula: see text], the maximal conductance of the M-current [Formula: see text] and the conductance of the leak current [Formula: see text]. We give precise conditions for the model that ensure the existence of a Bogdanov-Takens (BT) point and show that such a point can occur by varying [Formula: see text] and [Formula: see text]. We discuss the case when the BT point becomes a Bogdanov-Takens-cusp (BTC) point and show that such a point can occur in the three-dimensional parameter space. The results of the bifurcation analysis are applied to different neuronal models and are verified and supplemented by numerical bifurcation diagrams generated using the package MATCONT. We conclude that there is a transition in the neuronal excitability type organised by the BT point and the neuron switches from Class-I to Class-II as conductance of the M-current increases.

Effects of Neuromodulation on Excitatory-Inhibitory Neural Network Dynamics Depend on Network Connectivity Structure

Acetylcholine (ACh), one of the brain's most potent neuromodulators, can affect intrinsic neuron properties through blockade of an M-type potassium current. The effect of ACh on excitatory and inhibitory cells with this potassium channel modulates their membrane excitability, which in turn affects their tendency to synchronize in networks. Here, we study the resulting changes in dynamics in networks with inter-connected excitatory and inhibitory populations (E-I networks), which are ubiquitous in the brain. Utilizing biophysical models of E-I networks, we analyze how the network connectivity structure in terms of synaptic connectivity alters the influence of ACh on the generation of synchronous excitatory bursting. We investigate networks containing all combinations of excitatory and inhibitory cells with high (Type I properties) or low (Type II properties) modulatory tone. To vary network connectivity structure, we focus on the effects of the strengths of inter-connections between excitatory and inhibitory cells (E-I synapses and I-E synapses), and the strengths of intra-connections among excitatory cells (E-E synapses) and among inhibitory cells (I-I synapses). We show that the presence of ACh may or may not affect the generation of network synchrony depending on the network connectivity. Specifically, strong network inter-connectivity induces synchronous excitatory bursting regardless of the cellular propensity for synchronization, which aligns with predictions of the PING model. However, when a network's intra-connectivity dominates its inter-connectivity, the propensity for synchrony of either inhibitory or excitatory cells can determine the generation of network-wide bursting.

Phasic cholinergic signaling promotes emergence of local gamma rhythms in excitatory-inhibitory networks

Recent experimental results have shown that the detection of cues in behavioral attention tasks relies on transient increases of acetylcholine (ACh) release in frontal cortex and cholinergically driven oscillatory activity in the gamma frequency band (Howe et al. Journal of Neuroscience, 2017, 37, 3215). The cue-induced gamma rhythmic activity requires stimulation of M1 muscarinic receptors. Using biophysical computational modeling, we show that a network of excitatory (E) and inhibitory (I) neurons that initially displays asynchronous firing can generate transient gamma oscillatory activity in response to simulated brief pulses of ACh. ACh effects are simulated as transient modulation of the conductance of an M-type K⁺ current which is blocked by activation of muscarinic receptors and has significant effects on neuronal excitability. The ACh-induced effects on the M current conductance, g_Ks , change network dynamics to promote the emergence of network gamma rhythmicity through a Pyramidal-Interneuronal Network Gamma mechanism. Depending on connectivity strengths between and among E and I cells, gamma activity decays with the simulated g_Ks transient modulation or is sustained in the network after the g_Ks transient has completely dissipated. We investigated the sensitivity of the emergent gamma activity to synaptic strengths, external noise and simulated levels of g_Ks modulation. To address recent experimental findings that cholinergic signaling is likely spatially focused and dynamic, we show that localized g_Ks modulation can induce transient changes of cellular excitability in local subnetworks, subsequently causing population-specific gamma oscillations. These results highlight dynamical mechanisms underlying localization of ACh-driven responses and suggest that spatially localized, cholinergically induced gamma may contribute to selectivity in the processing of competing external stimuli, as occurs in attentional tasks.

Emergence of global synchronization in directed excitatory networks of type I neurons

The collective behaviour of neural networks depends on the cellular and synaptic properties of the neurons. The phase-response curve (PRC) is an experimentally obtainable measure of cellular properties that quantifies the shift in the next spike time of a neuron as a function of the phase at which stimulus is delivered to that neuron. The neuronal PRCs can be classified as having either purely positive values (type I) or distinct positive and negative regions (type II). Networks of type 1 PRCs tend not to synchronize via mutual excitatory synaptic connections. We study the synchronization properties of identical type I and type II neurons, assuming unidirectional synapses. Performing the linear stability analysis and the numerical simulation of the extended Kuramoto model, we show that feedforward loop motifs favour synchronization of type I excitatory and inhibitory neurons, while feedback loop motifs destroy their synchronization tendency. Moreover, large directed networks, either without feedback motifs or with many of them, have been constructed from the same undirected backbones, and a high synchronization level is observed for directed acyclic graphs with type I neurons. It has been shown that, the synchronizability of type I neurons depends on both the directionality of the network connectivity and the topology of its undirected backbone. The abundance of feedforward motifs enhances the synchronizability of the directed acyclic graphs.

Network and cellular mechanisms underlying heterogeneous excitatory/inhibitory balanced states

Recent work has explored spatiotemporal relationships between excitatory (E) and inhibitory (I) signaling within neural networks, and the effect of these relationships on network activity patterns. Data from these studies have indicated that excitation and inhibition are maintained at a similar level across long time periods and that excitatory and inhibitory currents may be tightly synchronized. Disruption of this balance-leading to an aberrant E/I ratio-is implicated in various brain pathologies. However, a thorough characterization of the relationship between E and I currents in experimental settings is largely impossible, due to their tight regulation at multiple cellular and network levels. Here, we use biophysical neural network models to investigate the emergence and properties of balanced states by heterogeneous mechanisms. Our results show that a network can homeostatically regulate the E/I ratio through interactions among multiple cellular and network factors, including average firing rates, synaptic weights and average neural depolarization levels in excitatory/inhibitory populations. Complex and competing interactions between firing rates and depolarization levels allow these factors to alternately dominate network dynamics in different synaptic weight regimes. This leads to the emergence of distinct mechanisms responsible for determining a balanced state and its dynamical correlate. Our analysis provides a comprehensive picture of how E/I ratio changes when manipulating specific network properties, and identifies the mechanisms regulating E/I balance. These results provide a framework to explain the diverse, and in some cases, contradictory experimental observations on the E/I state in different brain states and conditions.

Acetylcholine Mediates Dynamic Switching Between Information Coding Schemes in Neuronal Networks

Rate coding and phase coding are the two major coding modes seen in the brain. For these two modes, network dynamics must either have a wide distribution of frequencies for rate coding, or a narrow one to achieve stability in phase dynamics for phase coding. Acetylcholine (ACh) is a potent regulator of neural excitability. Acting through the muscarinic receptor, ACh reduces the magnitude of the potassium M-current, a hyperpolarizing current that builds up as neurons fire. The M-current contributes to several excitability features of neurons, becoming a major player in facilitating the transition between Type 1 (integrator) and Type 2 (resonator) excitability. In this paper we argue that this transition enables a dynamic switch between rate coding and phase coding as levels of ACh release change. When a network is in a high ACh state variations in synaptic inputs will lead to a wider distribution of firing rates across the network and this distribution will reflect the network structure or pattern of external input to the network. When ACh is low, network frequencies become narrowly distributed and the structure of a network or pattern of external inputs will be represented through phase relationships between firing neurons. This work provides insights into how modulation of neuronal features influences network dynamics and information processing across brain states.

Synchronization of Electrically Coupled Resonate-and-Fire Neurons

Electrical coupling between neurons is broadly present across brain areas and is typically assumed to synchronize network activity. However, intrinsic properties of the coupled cells can complicate this simple picture. Many cell types with electrical coupling show a diversity of post-spike subthreshold fluctuations, often linked to subthreshold resonance, which are transmitted through electrical synapses in addition to action potentials. Using the theory of weakly coupled oscillators, we explore the effect of both subthreshold and spike-mediated coupling on synchrony in small networks of electrically coupled resonate-and-fire neurons, a hybrid neuron model with damped subthreshold oscillations and a range of post-spike voltage dynamics. We calculate the phase response curve using an extension of the adjoint method that accounts for the discontinuous post-spike reset rule. We find that both spikes and subthreshold fluctuations can jointly promote synchronization. The subthreshold contribution is strongest when the voltage exhibits a significant post-spike elevation in voltage, or plateau potential. Additionally, we show that the geometry of trajectories approaching the spiking threshold causes a "reset-induced shear" effect that can oppose synchrony in the presence of network asymmetry, despite having no effect on the phase-locking of symmetrically coupled pairs.

How rhythms of the sleeping brain tune memory and synaptic plasticity

Decades of neurobehavioral research has linked sleep-associated rhythms in various brain areas to improvements in cognitive performance. However, it remains unclear what synaptic changes might underlie sleep-dependent declarative memory consolidation and procedural task improvement, and why these same changes appear not to occur across a similar interval of wake. Here we describe recent research on how one specific feature of sleep-network rhythms characteristic of rapid eye movement and non-rapid eye movement-could drive synaptic strengthening or weakening in specific brain circuits. We provide an overview of how these rhythms could affect synaptic plasticity individually and in concert. We also present an overarching hypothesis for how all network rhythms occurring across the sleeping brain could aid in encoding new information in neural circuits.

How single neuron properties shape chaotic dynamics and signal transmission in random neural networks

While most models of randomly connected neural networks assume single-neuron models with simple dynamics, neurons in the brain exhibit complex intrinsic dynamics over multiple timescales. We analyze how the dynamical properties of single neurons and recurrent connections interact to shape the effective dynamics in large randomly connected networks. A novel dynamical mean-field theory for strongly connected networks of multi-dimensional rate neurons shows that the power spectrum of the network activity in the chaotic phase emerges from a nonlinear sharpening of the frequency response function of single neurons. For the case of two-dimensional rate neurons with strong adaptation, we find that the network exhibits a state of "resonant chaos", characterized by robust, narrow-band stochastic oscillations. The coherence of stochastic oscillations is maximal at the onset of chaos and their correlation time scales with the adaptation timescale of single units. Surprisingly, the resonance frequency can be predicted from the properties of isolated neurons, even in the presence of heterogeneity in the adaptation parameters. In the presence of these internally-generated chaotic fluctuations, the transmission of weak, low-frequency signals is strongly enhanced by adaptation, whereas signal transmission is not influenced by adaptation in the non-chaotic regime. Our theoretical framework can be applied to other mechanisms at the level of single neurons, such as synaptic filtering, refractoriness or spike synchronization. These results advance our understanding of the interaction between the dynamics of single units and recurrent connectivity, which is a fundamental step toward the description of biologically realistic neural networks.

Macroscopic phase resetting-curves determine oscillatory coherence and signal transfer in inter-coupled neural circuits

Macroscopic oscillations of different brain regions show multiple phase relationships that are persistent across time and have been implicated in routing information. While multiple cellular mechanisms influence the network oscillatory dynamics and structure the macroscopic firing motifs, one of the key questions is to identify the biophysical neuronal and synaptic properties that permit such motifs to arise. A second important issue is how the different neural activity coherence states determine the communication between the neural circuits. Here we analyse the emergence of phase-locking within bidirectionally delayed-coupled spiking circuits in which global gamma band oscillations arise from synaptic coupling among largely excitable neurons. We consider both the interneuronal (ING) and the pyramidal-interneuronal (PING) population gamma rhythms and the inter coupling targeting the pyramidal or the inhibitory neurons. Using a mean-field approach together with an exact reduction method, we reduce each spiking network to a low dimensional nonlinear system and derive the macroscopic phase resetting-curves (mPRCs) that determine how the phase of the global oscillation responds to incoming perturbations. This is made possible by the use of the quadratic integrate-and-fire model together with a Lorentzian distribution of the bias current. Depending on the type of gamma (PING vs. ING), we show that incoming excitatory inputs can either speed up the macroscopic oscillation (phase advance; type I PRC) or induce both a phase advance and a delay (type II PRC). From there we determine the structure of macroscopic coherence states (phase-locking) of two weakly synaptically-coupled networks. To do so we derive a phase equation for the coupled system which links the synaptic mechanisms to the coherence states of the system. We show that a synaptic transmission delay is a necessary condition for symmetry breaking, i.e. a non-symmetric phase lag between the macroscopic oscillations. This potentially provides an explanation to the experimentally observed variety of gamma phase-locking modes. Our analysis further shows that symmetry-broken coherence states can lead to a preferred direction of signal transfer between the oscillatory networks where this directionality also depends on the timing of the signal. Hence we suggest a causal theory for oscillatory modulation of functional connectivity between cortical circuits.

Large scale brain oscillations emerge from synaptic interactions within neuronal circuits. Over the past years, such macroscopic rhythms have been suggested to play a crucial role in routing the flow of information across cortical regions, resulting in a functional connectome. The underlying mechanism is cortical oscillations that bind together following a well-known motif called phase-locking. While there is significant experimental support for multiple phase-locking modes in the brain, it is still unclear what is the underlying mechanism that permits macroscopic rhythms to phase lock. In the present paper we take up with this issue, and to show that, one can study the emergent macroscopic phase-locking within the mathematical framework of weakly coupled oscillators. We find that under synaptic delays, fully symmetrically coupled networks can display symmetry-broken states of activity, where one network starts to lead in phase the second (also sometimes known as stuttering states). When we analyse how incoming transient signals affect the coupled system, we find that in the symmetry-broken state, the effect depends strongly on which network is targeted (the leader or the follower) as well as the timing of the input. Hence we show how the dynamics of the emergent phase-locked activity imposes a functional directionality on how signals are processed. We thus offer clarification on the synaptic and circuit properties responsible for the emergence of multiple phase-locking patterns and provide support for its functional implication in information transfer.

Population dynamics and entrainment of basal ganglia pacemakers are shaped by their dendritic arbors

The theory of phase oscillators is an essential tool for understanding population dynamics of pacemaking neurons. GABAergic pacemakers in the substantia nigra pars reticulata (SNr), a main basal ganglia (BG) output nucleus, receive inputs from the direct and indirect pathways at distal and proximal regions of their dendritic arbors, respectively. We combine theory, optogenetic stimulation and electrophysiological experiments in acute brain slices to ask how dendritic properties impact the propensity of the various inputs, arriving at different locations along the dendrite, to recruit or entrain SNr pacemakers. By combining cable theory with sinusoidally-modulated optogenetic activation of either proximal somatodendritic regions or the entire somatodendritic arbor of SNr neurons, we construct an analytical model that accurately fits the empirically measured somatic current response to inputs arising from illuminating the soma and various portions of the dendritic field. We show that the extent of the dendritic tree that is illuminated generates measurable and systematic differences in the pacemaker’s phase response curve (PRC), causing a shift in its peak. Finally, we show that the divergent PRCs correctly predict differences in two major features of the collective dynamics of SNr neurons: the fidelity of population responses to sudden step-like changes in inputs; and the phase latency at which SNr neurons are entrained by rhythmic stimulation, which can occur in the BG under both physiological and pathophysiological conditions. Our novel method generates measurable and physiologically meaningful spatial effects, and provides the first empirical demonstration of how the collective responses of SNr pacemakers are determined by the transmission properties of their dendrites. SNr dendrites may serve to delay distal striatal inputs so that they impinge on the spike initiation zone simultaneously with pallidal and subthalamic inputs in order to guarantee a fair competition between the influence of the monosynaptic direct- and polysynaptic indirect pathways.

The substantia nigra pars reticulata (SNr) is a main output nucleus of the basal ganglia (BG), where inputs from the competing direct and indirect pathways converge onto the same neurons. Interestingly, these inputs are differentially distributed with direct and indirect pathway projections arriving at distal and proximal regions of the dendritic arbor, respectively. We employ a novel method combining theory with electrophysiological experiments and optogenetics to study the distinct effects of inputs arriving at different locations along the dendrite. Our approach represents a useful compromise between complexity and reduction in modelling. Our work addresses the question of high fidelity encoding of inputs by networks of neurons in the new context of pacemaking neurons, which are driven to fire by their intrinsic dynamics rather than by a network state. We provide the first empirical demonstration that dendritic delays can introduce latencies in the responses of a population of neurons that are commensurate with synaptic delays, suggesting a new role for SNr dendrites with implications for BG function.

Data-Driven Modeling of Cholinergic Modulation of Neural Microcircuits: Bridging Neurons, Synapses and Network Activity

Neuromodulators, such as acetylcholine (ACh), control information processing in neural microcircuits by regulating neuronal and synaptic physiology. Computational models and simulations enable predictions on the potential role of ACh in reconfiguring network activity. As a prelude into investigating how the cellular and synaptic effects of ACh collectively influence emergent network dynamics, we developed a data-driven framework incorporating phenomenological models of the physiology of cholinergic modulation of neocortical cells and synapses. The first-draft models were integrated into a biologically detailed tissue model of neocortical microcircuitry to investigate the effects of levels of ACh on diverse neuron types and synapses, and consequently on emergent network activity. Preliminary simulations from the framework, which was not tuned to reproduce any specific ACh-induced network effects, not only corroborate the long-standing notion that ACh desynchronizes spontaneous network activity, but also predict that a dose-dependent activation of ACh gives rise to a spectrum of neocortical network activity. We show that low levels of ACh, such as during non-rapid eye movement (nREM) sleep, drive microcircuit activity into slow oscillations and network synchrony, whereas high ACh concentrations, such as during wakefulness and REM sleep, govern fast oscillations and network asynchrony. In addition, spontaneous network activity modulated by ACh levels shape spike-time cross-correlations across distinct neuronal populations in strikingly different ways. These effects are likely due to the regulation of neurons and synapses caused by increasing levels of ACh, which enhances cellular excitability and decreases the efficacy of local synaptic transmission. We conclude by discussing future directions to refine the biological accuracy of the framework, which will extend its utility and foster the development of hypotheses to investigate the role of neuromodulators in neural information processing.

The impact of spike-frequency adaptation on balanced network dynamics

A dynamic balance between strong excitatory and inhibitory neuronal inputs is hypothesized to play a pivotal role in information processing in the brain. While there is evidence of the existence of a balanced operating regime in several cortical areas and idealized neuronal network models, it is important for the theory of balanced networks to be reconciled with more physiological neuronal modeling assumptions. In this work, we examine the impact of spike-frequency adaptation, observed widely across neurons in the brain, on balanced dynamics. We incorporate adaptation into binary and integrate-and-fire neuronal network models, analyzing the theoretical effect of adaptation in the large network limit and performing an extensive numerical investigation of the model adaptation parameter space. Our analysis demonstrates that balance is well preserved for moderate adaptation strength even if the entire network exhibits adaptation. In the common physiological case in which only excitatory neurons undergo adaptation, we show that the balanced operating regime in fact widens relative to the non-adaptive case. We hypothesize that spike-frequency adaptation may have been selected through evolution to robustly facilitate balanced dynamics across diverse cognitive operating states.

Resonance with subthreshold oscillatory drive organizes activity and optimizes learning in neural networks

Networks of neurons need to reliably encode and replay patterns and sequences of activity. In the brain, sequences of spatially coding neurons are replayed in both the forward and reverse direction in time with respect to their order in recent experience. As of yet there is no network-level or biophysical mechanism known that can produce both modes of replay within the same network. Here we propose that resonance, a property of neurons, paired with subthreshold oscillations in neural input facilitate network-level learning of fixed and sequential activity patterns and lead to both forward and reverse replay.

Network oscillations across and within brain areas are critical for learning and performance of memory tasks. While a large amount of work has focused on the generation of neural oscillations, their effect on neuronal populations’ spiking activity and information encoding is less known. Here, we use computational modeling to demonstrate that a shift in resonance responses can interact with oscillating input to ensure that networks of neurons properly encode new information represented in external inputs to the weights of recurrent synaptic connections. Using a neuronal network model, we find that due to an input current-dependent shift in their resonance response, individual neurons in a network will arrange their phases of firing to represent varying strengths of their respective inputs. As networks encode information, neurons fire more synchronously, and this effect limits the extent to which further “learning” (in the form of changes in synaptic strength) can occur. We also demonstrate that sequential patterns of neuronal firing can be accurately stored in the network; these sequences are later reproduced without external input (in the context of subthreshold oscillations) in both the forward and reverse directions (as has been observed following learning in vivo). To test whether a similar mechanism could act in vivo, we show that periodic stimulation of hippocampal neurons coordinates network activity and functional connectivity in a frequency-dependent manner. We conclude that resonance with subthreshold oscillations provides a plausible network-level mechanism to accurately encode and retrieve information without overstrengthening connections between neurons.

Synaptic Impairment and Robustness of Excitatory Neuronal Networks with Different Topologies

Synaptic deficiencies are a known hallmark of neurodegenerative diseases, but the diagnosis of impaired synapses on the cellular level is not an easy task. Nonetheless, changes in the system-level dynamics of neuronal networks with damaged synapses can be detected using techniques that do not require high spatial resolution. This paper investigates how the structure/topology of neuronal networks influences their dynamics when they suffer from synaptic loss. We study different neuronal network structures/topologies by specifying their degree distributions. The modes of the degree distribution can be used to construct networks that consist of rich clubs and resemble small world networks, as well. We define two dynamical metrics to compare the activity of networks with different structures: persistent activity (namely, the self-sustained activity of the network upon removal of the initial stimulus) and quality of activity (namely, percentage of neurons that participate in the persistent activity of the network). Our results show that synaptic loss affects the persistent activity of networks with bimodal degree distributions less than it affects random networks. The robustness of neuronal networks enhances when the distance between the modes of the degree distribution increases, suggesting that the rich clubs of networks with distinct modes keep the whole network active. In addition, a tradeoff is observed between the quality of activity and the persistent activity. For a range of distributions, both of these dynamical metrics are considerably high for networks with bimodal degree distribution compared to random networks. We also propose three different scenarios of synaptic impairment, which may correspond to different pathological or biological conditions. Regardless of the network structure/topology, results demonstrate that synaptic loss has more severe effects on the activity of the network when impairments are correlated with the activity of the neurons.

Qualitative changes in phase-response curve and synchronization at the saddle-node-loop bifurcation

Prominent changes in neuronal dynamics have previously been attributed to a specific switch in onset bifurcation, the Bogdanov-Takens (BT) point. This study unveils another, relevant and so far underestimated transition point: the saddle-node-loop bifurcation, which can be reached by several parameters, including capacitance, leak conductance, and temperature. This bifurcation turns out to induce even more drastic changes in synchronization than the BT transition. This result arises from a direct effect of the saddle-node-loop bifurcation on the limit cycle and hence spike dynamics. In contrast, the BT bifurcation exerts its immediate influence upon the subthreshold dynamics and hence only indirectly relates to spiking. We specifically demonstrate that the saddle-node-loop bifurcation (i) ubiquitously occurs in planar neuron models with a saddle node on invariant cycle onset bifurcation, and (ii) results in a symmetry breaking of the system's phase-response curve. The latter entails an increase in synchronization range in pulse-coupled oscillators, such as neurons. The derived bifurcation structure is of interest in any system for which a relaxation limit is admissible, such as Josephson junctions and chemical oscillators.

Mathematical Frameworks for Oscillatory Network Dynamics in Neuroscience

The tools of weakly coupled phase oscillator theory have had a profound impact on the neuroscience community, providing insight into a variety of network behaviours ranging from central pattern generation to synchronisation, as well as predicting novel network states such as chimeras. However, there are many instances where this theory is expected to break down, say in the presence of strong coupling, or must be carefully interpreted, as in the presence of stochastic forcing. There are also surprises in the dynamical complexity of the attractors that can robustly appear-for example, heteroclinic network attractors. In this review we present a set of mathematical tools that are suitable for addressing the dynamics of oscillatory neural networks, broadening from a standard phase oscillator perspective to provide a practical framework for further successful applications of mathematics to understanding network dynamics in neuroscience.

Interplay between excitability type and distributions of neuronal connectivity determines neuronal network synchronization

While the interplay between neuronal excitability properties and global properties of network topology is known to affect network propensity for synchronization, it is not clear how detailed characteristics of these properties affect spatiotemporal pattern formation. Here we study mixed networks, composed of neurons having type I and/or type II phase response curves, with varying distributions of local and random connections and show that not only average network properties, but also the connectivity distribution statistics, significantly affect network synchrony. Namely, we study networks with fixed networkwide properties, but vary the number of random connections that nodes project. We show that varying node excitability (type I vs type II) influences network synchrony most dramatically for systems with long-tailed distributions of the number of random connections per node. This indicates that a cluster of even a few highly rewired cells with a high propensity for synchronization can alter the degree of synchrony in the network as a whole. We show this effect generally on a network of coupled Kuramoto oscillators and investigate the impact of this effect more thoroughly in pulse-coupled networks of biophysical neurons.

Neurons as oscillators

Regularly spiking neurons can be described as oscillators. In this article we review some of the insights gained from this conceptualization and their relevance for systems neuroscience. First, we explain how a regularly spiking neuron can be viewed as an oscillator and how the phase-response curve (PRC) describes the response of the neuron's spike times to small perturbations. We then discuss the meaning of the PRC for a single neuron's spiking behavior and review the PRCs measured from a variety of neurons in a range of spiking regimes. Next, we show how the PRC can be related to a number of common measures used to quantify neuronal firing, such as the spike-triggered average and the peristimulus histogram. We further show that the response of a neuron to correlated inputs depends on the shape of the PRC. We then explain how the PRC of single neurons can be used to predict neural network behavior. Given the PRC, conduction delays, and the waveform and time course of the synaptic potentials, it is possible to predict neural population behavior such as synchronization. The PRC also allows us to quantify the robustness of the synchronization to heterogeneity and noise. We finally ask how to combine the measured PRCs and the predictions based on PRC to further the understanding of systems neuroscience. As an example, we discuss how the change of the PRC by the neuromodulator acetylcholine could lead to a destabilization of cortical network dynamics. Although all of these studies are grounded in mathematical abstractions that do not strictly hold in biology, they provide good estimates for the emergence of the brain's network activity from the properties of individual neurons. The study of neurons as oscillators can provide testable hypotheses and mechanistic explanations for systems neuroscience.

Jittering waves in rings of pulse oscillators

Rings of oscillators with delayed pulse coupling are studied analytically, numerically, and experimentally. The basic regimes observed in such rings are rotating waves with constant interspike intervals and phase lags between the neighbors. We show that these rotating waves may destabilize leading to the so-called jittering waves. For these regimes, the interspike intervals are no more equal but form a periodic sequence in time. Analytic criterion for the emergence of jittering waves is derived and confirmed by the numerical and experimental data. The obtained results contribute to the hypothesis that the multijitter instability is universal in systems with pulse coupling.

Memory recall and spike-frequency adaptation

The brain can reproduce memories from partial data; this ability is critical for memory recall. The process of memory recall has been studied using autoassociative networks such as the Hopfield model. This kind of model reliably converges to stored patterns that contain the memory. However, it is unclear how the behavior is controlled by the brain so that after convergence to one configuration, it can proceed with recognition of another one. In the Hopfield model, this happens only through unrealistic changes of an effective global temperature that destabilizes all stored configurations. Here we show that spike-frequency adaptation (SFA), a common mechanism affecting neuron activation in the brain, can provide state-dependent control of pattern retrieval. We demonstrate this in a Hopfield network modified to include SFA, and also in a model network of biophysical neurons. In both cases, SFA allows for selective stabilization of attractors with different basins of attraction, and also for temporal dynamics of attractor switching that is not possible in standard autoassociative schemes. The dynamics of our models give a plausible account of different sorts of memory retrieval.

Network burst dynamics under heterogeneous cholinergic modulation of neural firing properties and heterogeneous synaptic connectivity

The characteristics of neural network activity depend on intrinsic neural properties and synaptic connectivity in the network. In brain networks, both of these properties are critically affected by the type and levels of neuromodulators present. The expression of many of the most powerful neuromodulators, including acetylcholine (ACh), varies tonically and phasically with behavioural state, leading to dynamic, heterogeneous changes in intrinsic neural properties and synaptic connectivity properties. Namely, ACh significantly alters neural firing properties as measured by the phase response curve in a manner that has been shown to alter the propensity for network synchronization. The aim of this simulation study was to build an understanding of how heterogeneity in cholinergic modulation of neural firing properties and heterogeneity in synaptic connectivity affect the initiation and maintenance of synchronous network bursting in excitatory networks. We show that cells that display different levels of ACh modulation have differential roles in generating network activity: weakly modulated cells are necessary for burst initiation and provide synchronizing drive to the rest of the network, whereas strongly modulated cells provide the overall activity level necessary to sustain burst firing. By applying several quantitative measures of network activity, we further show that the existence of network bursting and its characteristics, such as burst duration and intraburst synchrony, are dependent on the fraction of cell types providing the synaptic connections in the network. These results suggest mechanisms underlying ACh modulation of brain oscillations and the modulation of seizure activity during sleep states.

Formation and Dynamics of Waves in a Cortical Model of Cholinergic Modulation

Acetylcholine (ACh) is a regulator of neural excitability and one of the neurochemical substrates of sleep. Amongst the cellular effects induced by cholinergic modulation are a reduction in spike-frequency adaptation (SFA) and a shift in the phase response curve (PRC). We demonstrate in a biophysical model how changes in neural excitability and network structure interact to create three distinct functional regimes: localized asynchronous, traveling asynchronous, and traveling synchronous. Our results qualitatively match those observed experimentally. Cortical activity during slow wave sleep (SWS) differs from that during REM sleep or waking states. During SWS there are traveling patterns of activity in the cortex; in other states stationary patterns occur. Our model is a network composed of Hodgkin-Huxley type neurons with a M-current regulated by ACh. Regulation of ACh level can account for dynamical changes between functional regimes. Reduction of the magnitude of this current recreates the reduction in SFA the shift from a type 2 to a type 1 PRC observed in the presence of ACh. When SFA is minimal (in waking or REM sleep state, high ACh) patterns of activity are localized and easily pinned by network inhomogeneities. When SFA is present (decreasing ACh), traveling waves of activity naturally arise. A further decrease in ACh leads to a high degree of synchrony within traveling waves. We also show that the level of ACh determines how sensitive network activity is to synaptic heterogeneity. These regimes may have a profound functional significance as stationary patterns may play a role in the proper encoding of external input as memory and traveling waves could lead to synaptic regularization, giving unique insights into the role and significance of ACh in determining patterns of cortical activity and functional differences arising from the patterns.

On the firing rate dependency of the phase response curve of rat Purkinje neurons in vitro

Synchronous spiking during cerebellar tasks has been observed across Purkinje cells: however, little is known about the intrinsic cellular mechanisms responsible for its initiation, cessation and stability. The Phase Response Curve (PRC), a simple input-output characterization of single cells, can provide insights into individual and collective properties of neurons and networks, by quantifying the impact of an infinitesimal depolarizing current pulse on the time of occurrence of subsequent action potentials, while a neuron is firing tonically. Recently, the PRC theory applied to cerebellar Purkinje cells revealed that these behave as phase-independent integrators at low firing rates, and switch to a phase-dependent mode at high rates. Given the implications for computation and information processing in the cerebellum and the possible role of synchrony in the communication with its post-synaptic targets, we further explored the firing rate dependency of the PRC in Purkinje cells. We isolated key factors for the experimental estimation of the PRC and developed a closed-loop approach to reliably compute the PRC across diverse firing rates in the same cell. Our results show unambiguously that the PRC of individual Purkinje cells is firing rate dependent and that it smoothly transitions from phase independent integrator to a phase dependent mode. Using computational models we show that neither channel noise nor a realistic cell morphology are responsible for the rate dependent shift in the phase response curve.

Changes of Firing Rate Induced by Changes of Phase Response Curve in Bifurcation Transitions

We study dynamical mechanisms responsible for changes of the firing rate during four different bifurcation transitions in the two-dimensional Hindmarsh-Rose (2DHR) neuron model: the saddle node on an invariant circle (SNIC) bifurcation to the supercritical Andronov-Hopf (AH) one, the SNIC bifurcation to the saddle-separatrix loop (SSL) one, the AH bifurcation to the subcritical AH (SAH) one, and the SSL bifurcation to the AH one. For this purpose, we study slopes of the firing rate curve with respect to not only an external input current but also temperature that can be interpreted as a timescale in the 2DHR neuron model. These slopes are mathematically formulated with phase response curves (PRCs), expanding the firing rate with perturbations of the temperature and external input current on the one-dimensional space of the phase [Formula: see text] in the 2DHR oscillator. By analyzing the two different slopes of the firing rate curve with respect to the temperature and external input current, we find that during changes of the firing rate in all of the bifurcation transitions, the calculated slope with respect to the temperature also changes. This is largely dependent on changes in the PRC size that is also related to the slope with respect to the external input current. Furthermore, we find phase transition–like switches of the firing rate with a possible increase of the temperature during the SSL-to-AH bifurcation transition.

Adaptation and shunting inhibition leads to pyramidal/interneuron gamma with sparse firing of pyramidal cells

Gamma oscillations are a prominent phenomenon related to a number of brain functions. Data show that individual pyramidal neurons can fire at rate below gamma with the population showing clear gamma oscillations and synchrony. In one kind of idealized model of such weak gamma, pyramidal neurons fire in clusters. Here we provide a theory for clustered gamma PING rhythms with strong inhibition and weaker excitation. Our simulations of biophysical models show that the adaptation of pyramidal neurons coupled with their low firing rate leads to cluster formation. A partially analytic study of a canonical model shows that the phase response curves with a near zero flat region, caused by the presence of the slow adaptive current, are the key to the formation of clusters. Furthermore we examine shunting inhibition and show that clusters become robust and generic.

Predicting the responses of repetitively firing neurons to current noise

We used phase resetting methods to predict firing patterns of rat subthalamic nucleus (STN) neurons when their rhythmic firing was densely perturbed by noise. We applied sequences of contiguous brief (0.5–2 ms) current pulses with amplitudes drawn from a Gaussian distribution (10–100 pA standard deviation) to autonomously firing STN neurons in slices. Current noise sequences increased the variability of spike times with little or no effect on the average firing rate. We measured the infinitesimal phase resetting curve (PRC) for each neuron using a noise-based method. A phase model consisting of only a firing rate and PRC was very accurate at predicting spike timing, accounting for more than 80% of spike time variance and reliably reproducing the spike-to-spike pattern of irregular firing. An approximation for the evolution of phase was used to predict the effect of firing rate and noise parameters on spike timing variability. It quantitatively predicted changes in variability of interspike intervals with variation in noise amplitude, pulse duration and firing rate over the normal range of STN spontaneous rates. When constant current was used to drive the cells to higher rates, the PRC was altered in size and shape and accurate predictions of the effects of noise relied on incorporating these changes into the prediction. Application of rate-neutral changes in conductance showed that changes in PRC shape arise from conductance changes known to accompany rate increases in STN neurons, rather than the rate increases themselves. Our results show that firing patterns of densely perturbed oscillators cannot readily be distinguished from those of neurons randomly excited to fire from the rest state. The spike timing of repetitively firing neurons may be quantitatively predicted from the input and their PRCs, even when they are so densely perturbed that they no longer fire rhythmically.

Most neurons receive thousands of synaptic inputs per second. Each of these may be individually weak but collectively they shape the temporal pattern of firing by the postsynaptic neuron. If the postsynaptic neuron fires repetitively, its synaptic inputs need not directly trigger action potentials, but may instead control the timing of action potentials that would occur anyway. The phase resetting curve encapsulates the influence of an input on the timing of the next action potential, depending on its time of arrival. We measured the phase resetting curves of neurons in the subthalamic nucleus and used them to accurately predict the timing of action potentials in a phase model subjected to complex input patterns. A simple approximation to the phase model accurately predicted the changes in firing pattern evoked by dense patterns of noise pulses varying in amplitude and pulse duration, and by changes in firing rate. We also showed that the phase resetting curve changes systematically with changes in total neuron conductance, and doing so predicts corresponding changes in firing pattern. Our results indicate that the phase model may accurately represent the temporal integration of complex patterns of input to repetitively firing neurons.

Effect of phase response curve skew on synchronization with and without conduction delays

A central problem in cortical processing including sensory binding and attentional gating is how neurons can synchronize their responses with zero or near-zero time lag. For a spontaneously firing neuron, an input from another neuron can delay or advance the next spike by different amounts depending upon the timing of the input relative to the previous spike. This information constitutes the phase response curve (PRC). We present a simple graphical method for determining the effect of PRC shape on synchronization tendencies and illustrate it using type 1 PRCs, which consist entirely of advances (delays) in response to excitation (inhibition). We obtained the following generic solutions for type 1 PRCs, which include the pulse-coupled leaky integrate and fire model. For pairs with mutual excitation, exact synchrony can be stable for strong coupling because of the stabilizing effect of the causal limit region of the PRC in which an input triggers a spike immediately upon arrival. However, synchrony is unstable for short delays, because delayed inputs arrive during a refractory period and cannot trigger an immediate spike. Right skew destabilizes antiphase and enables modes with time lags that grow as the conduction delay is increased. Therefore, right skew favors near synchrony at short conduction delays and a gradual transition between synchrony and antiphase for pairs coupled by mutual excitation. For pairs with mutual inhibition, zero time lag synchrony is stable for conduction delays ranging from zero to a substantial fraction of the period for pairs. However, for right skew there is a preferred antiphase mode at short delays. In contrast to mutual excitation, left skew destabilizes antiphase for mutual inhibition so that synchrony dominates at short delays as well. These pairwise synchronization tendencies constrain the synchronization properties of neurons embedded in larger networks.

How delays matter in an oscillatory whole-brain spiking-neuron network model for MEG alpha-rhythms at rest

In recent years the study of the intrinsic brain dynamics in a relaxed awake state in the absence of any specific task has gained increasing attention, as spontaneous neural activity has been found to be highly structured at a large scale. This so called resting-state activity has been found to be comprised by nonrandom spatiotemporal patterns and fluctuations, and several Resting-State Networks (RSN) have been found in BOLD-fMRI as well as in MEG signal power envelope correlations. The underlying anatomical connectivity structure between areas of the brain has been identified as being a key to the observed functional network connectivity, but the mechanisms behind this are still underdetermined. Theoretical large-scale brain models for fMRI data have corroborated the importance of the connectome in shaping network dynamics, while the importance of delays and noise differ between studies and depend on the models' specific dynamics. In the current study, we present a spiking neuron network model that is able to produce noisy, distributed alpha-oscillations, matching the power peak in the spectrum of group resting-state MEG recordings. We studied how well the model captured the inter-node correlation structure of the alpha-band power envelopes for different delays between brain areas, and found that the model performs best for propagation delays inside the physiological range (5-10 m/s). Delays also shift the transition from noisy to bursting oscillations to higher global coupling values in the model. Thus, in contrast to the asynchronous fMRI state, delays are important to consider in the presence of oscillation.

Impact of neuronal heterogeneity on correlated colored noise-induced synchronization

Synchronization plays an important role in neural signal processing and transmission. Many hypotheses have been proposed to explain the origin of neural synchronization. In recent years, correlated noise-induced synchronization has received support from many theoretical and experimental studies. However, many of these prior studies have assumed that neurons have identical biophysical properties and that their inputs are well modeled by white noise. In this context, we use colored noise to induce synchronization between oscillators with heterogeneity in both phase-response curves and frequencies. In the low noise limit, we derive novel analytical theory showing that the time constant of colored noise influences correlated noise-induced synchronization and that oscillator heterogeneity can limit synchronization. Surprisingly, however, heterogeneous oscillators may synchronize better than homogeneous oscillators given low input correlations. We also find resonance of oscillator synchronization to colored noise inputs when firing frequencies diverge. Collectively, these results prove robust for both relatively high noise regimes and when applied to biophysically realistic spiking neuron models, and further match experimental recordings from acute brain slices.

A dynamical role for acetylcholine in synaptic renormalization

Although sleep is a fundamental behavior observed in virtually all animal species, its functions remain unclear. One leading proposal, known as the synaptic renormalization hypothesis, suggests that sleep is necessary to counteract a global strengthening of synapses that occurs during wakefulness. Evidence for sleep-dependent synaptic downscaling (or synaptic renormalization) has been observed experimentally, but the physiological mechanisms which generate this phenomenon are unknown. In this study, we propose that changes in neuronal membrane excitability induced by acetylcholine may provide a dynamical mechanism for both wake-dependent synaptic upscaling and sleep-dependent downscaling. We show in silico that cholinergically-induced changes in network firing patterns alter overall network synaptic potentiation when synaptic strengths evolve through spike-timing dependent plasticity mechanisms. Specifically, network synaptic potentiation increases dramatically with high cholinergic concentration and decreases dramatically with low levels of acetylcholine. We demonstrate that this phenomenon is robust across variation of many different network parameters.

Parameterized phase response curves for characterizing neuronal behaviors under transient conditions

Phase response curves (PRCs) are a simple model of how a neuron's spike time is affected by synaptic inputs. PRCs are useful in predicting how networks of neurons behave when connected. One challenge in estimating a neuron's PRCs experimentally is that many neurons do not have stationary firing rates. In this article we introduce a new method to estimate PRCs as a function of firing rate of the neuron. We call the resulting model a parameterized PRC (pPRC). Experimentally, we perturb the neuron applying a current with two parts: 1) a current held constant between spikes but changed at the onset of a spike, used to make the neuron fire at different rates, and 2) a pulse to emulate a synaptic input. A model of the applied constant current and the history is made to predict the interspike interval (ISI). A second model is then made to fit the modulation of the spike time from the expected ISI by the pulsatile stimulus. A polynomial with two independent variables, the stimulus phase and the expected ISI, is used to model the pPRC. The pPRC is validated in a computational model and applied to pyramidal neurons from the CA1 region of the hippocampal slices from rat. The pPRC can be used to model the effect of changing firing rates on network synchrony. It can also be used to characterize the effects of neuromodulators and genetic mutations (among other manipulations) on network synchrony. It can also easily be extended to account for more variables.

Higher-order spike triggered analysis of neural oscillators

For the purpose of elucidating the neural coding process based on the neural excitability mechanism, researchers have recently investigated the relationship between neural dynamics and the spike triggered stimulus ensemble (STE). Ermentrout et al. analytically derived the relational equation between the phase response curve (PRC) and the spike triggered average (STA). The STA is the first cumulant of the STE. However, in order to understand the neural function as the encoder more explicitly, it is necessary to elucidate the relationship between the PRC and higher-order cumulants of the STE. In this paper, we give a general formulation to relate the PRC and the nth moment of the STE. By using this formulation, we derive a relational equation between the PRC and the spike triggered covariance (STC), which is the covariance of the STE. We show the effectiveness of the relational equation through numerical simulations and use the equation to identify the feature space of the rat hippocampal CA1 pyramidal neurons from their PRCs. Our result suggests that the hippocampal CA1 pyramidal neurons oscillating in the theta frequency range are commonly sensitive to inputs composed of theta and gamma frequency components.

When Long-Range Zero-Lag Synchronization is Feasible in Cortical Networks

Many studies have reported long-range synchronization of neuronal activity between brain areas, in particular in the beta and gamma bands with frequencies in the range of 14–30 and 40–80 Hz, respectively. Several studies have reported synchrony with zero phase lag, which is remarkable considering the synaptic and conduction delays inherent in the connections between distant brain areas. This result has led to many speculations about the possible functional role of zero-lag synchrony, such as for neuronal communication, attention, memory, and feature binding. However, recent studies using recordings of single-unit activity and local field potentials report that neuronal synchronization may occur with non-zero phase lags. This raises the questions whether zero-lag synchrony can occur in the brain and, if so, under which conditions. We used analytical methods and computer simulations to investigate which connectivity between neuronal populations allows or prohibits zero-lag synchrony. We did so for a model where two oscillators interact via a relay oscillator. Analytical results and computer simulations were obtained for both type I Mirollo–Strogatz neurons and type II Hodgkin–Huxley neurons. We have investigated the dynamics of the model for various types of synaptic coupling and importantly considered the potential impact of Spike-Timing Dependent Plasticity (STDP) and its learning window. We confirm previous results that zero-lag synchrony can be achieved in this configuration. This is much easier to achieve with Hodgkin–Huxley neurons, which have a biphasic phase response curve, than for type I neurons. STDP facilitates zero-lag synchrony as it adjusts the synaptic strengths such that zero-lag synchrony is feasible for a much larger range of parameters than without STDP.

Intrinsic heterogeneity in oscillatory dynamics limits correlation-induced neural synchronization

Synchronous neural oscillations are found throughout the brain and are thought to contribute to neural coding and the propagation of activity. Several proposed mechanisms of synchronization have gained support through combined theoretical and experimental investigation, including mechanisms based on coupling and correlated input. Here, we ask how correlation-induced synchrony is affected by physiological heterogeneity across neurons. To address this question, we examined cell-to-cell differences in phase-response curves (PRCs), which characterize the response of periodically firing neurons to weak perturbations. Using acute slice electrophysiology, we measured PRCs across a single class of principal neurons capable of sensory-evoked oscillations in vivo: the olfactory bulb mitral cells (MCs). Periodically firing MCs displayed a broad range of PRCs, each of which was well fit by a simple three-parameter model. MCs also displayed differences in firing rate-current relationships and in preferred firing rate ranges. Both the observed PRC heterogeneity and moderate firing rate differences (∼10 Hz) separately reduced the maximum correlation-induced synchrony between MCs by up to 25-30%. Simulations further demonstrated that these components of heterogeneity alone were sufficient to account for the difference in synchronization among heterogeneous vs. homogeneous populations in vitro. Within this simulation framework, independent modulation of specific PRC features additionally revealed which aspects of PRC heterogeneity most strongly impact correlation-induced synchronization. Finally, we demonstrated good agreement of novel mathematical theory with our experimental and simulation results, providing a theoretical basis for the influence of heterogeneity on correlation-induced neural synchronization.

Identifying type I excitability using dynamics of stochastic neural firing patterns

The stochastic firing patterns are simulated near saddle-node bifurcation on an invariant cycle corresponding to type I excitability in stochastic Morris-Lecar model. In absence of external periodic signal, the stochastic firing manifests continuous distribution in ISI histogram (ISIH), whose amplitude at first increases sharply and then decreases exponentially. In presence of the external periodic signal, stochastic firing patterns appear as two cases of integer multiple firing with multiple discrete peaks in ISIH. One manifests perfect exponential decay in all peaks and the other imperfect exponential decay except a lower first peak. These stochastic firing patterns simulated with or without external periodic signal can be demonstrated in the experiments on rat hippocampal CA1 pyramidal neurons. The exponential decay laws in the multiple peaks are also acquired using probability analysis method. The perfect decay law is determined by the independent characteristic within the firing while the imperfect decay law is from the inhibitory effect. In addition, the stochastic firing patterns corresponding to type I excitability are compared to those of type II excitability. The results not only reveal the dynamics of stochastic firing patterns with or without external signal corresponding to type I excitability, but also provide practical indicators to availably identify type I excitability.