skip to main content


Title: Cross-scale excitability in networks of quadratic integrate-and-fire neurons
From the action potentials of neurons and cardiac cells to the amplification of calcium signals in oocytes, excitability is a hallmark of many biological signalling processes. In recent years, excitability in single cells has been related to multiple-timescale dynamics through canards , special solutions which determine the effective thresholds of the all-or-none responses. However, the emergence of excitability in large populations remains an open problem. Here, we show that the mechanism of excitability in large networks and mean-field descriptions of coupled quadratic integrate-and-fire (QIF) cells mirrors that of the individual components. We initially exploit the Ott-Antonsen ansatz to derive low-dimensional dynamics for the coupled network and use it to describe the structure of canards via slow periodic forcing. We demonstrate that the thresholds for onset and offset of population firing can be found in the same way as those of the single cell. We combine theoretical analysis and numerical computations to develop a novel and comprehensive framework for excitability in large populations, applicable not only to models amenable to Ott-Antonsen reduction, but also to networks without a closed-form mean-field limit, in particular sparse networks.  more » « less
Award ID(s):
1951099
NSF-PAR ID:
10420009
Author(s) / Creator(s):
; ;
Editor(s):
Blackwell, Kim T.
Date Published:
Journal Name:
PLOS Computational Biology
Volume:
18
Issue:
10
ISSN:
1553-7358
Page Range / eLocation ID:
e1010569
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Synchronization in the gamma band (25–150 Hz) is mediated by PV+ inhibitory interneurons, and evidence is accumulating for the essential role of gamma oscillations in cognition. Oscillations can arise in inhibitory networks via synaptic interactions between individual oscillatory neurons (mean-driven) or via strong recurrent inhibition that destabilizes the stationary background firing rate in the fluctuation-driven balanced state, causing an oscillation in the population firing rate. Previous theoretical work focused on model neurons with Hodgkin's Type 1 excitability (integrators) connected by current-based synapses. Here we show that networks comprised of simple Type 2 oscillators (resonators) exhibit a supercritical Hopf bifurcation between synchrony and asynchrony and a gradual transition via cycle skipping from coupled oscillators to stochastic population oscillator (SPO), as previously shown for Type 1. We extended our analysis to homogeneous networks with conductance rather than current based synapses and found that networks with hyperpolarizing inhibitory synapses were more robust to noise than those with shunting synapses, both in the coupled oscillator and SPO regime. Assuming that reversal potentials are uniformly distributed between shunting and hyperpolarized values, as observed in one experimental study, converting synapses to purely hyperpolarizing favored synchrony in all cases, whereas conversion to purely shunting synapses made synchrony less robust except at very high conductance strengths. In mature neurons the synaptic reversal potential is controlled by chloride cotransporters that control the intracellular concentrations of chloride and bicarbonate ions, suggesting these transporters as a potential therapeutic target to enhance gamma synchrony and cognition.

     
    more » « less
  2. Neural networks with a large number of units ad- mit a mean-field description, which has recently served as a theoretical explanation for the favor- able training properties of “overparameterized” models. In this regime, gradient descent obeys a deterministic partial differential equation (PDE) that converges to a globally optimal solution for networks with a single hidden layer under appro- priate assumptions. In this work, we propose a non-local mass transport dynamics that leads to a modified PDE with the same minimizer. We im- plement this non-local dynamics as a stochastic neuronal birth-death process and we prove that it accelerates the rate of convergence in the mean- field limit. We subsequently realize this PDE with two classes of numerical schemes that converge to the mean-field equation, each of which can easily be implemented for neural networks with finite numbers of units. We illustrate our algorithms with two models to provide intuition for the mech- anism through which convergence is accelerated 
    more » « less
  3. Many cells can generate complementary traveling waves of actin filaments (F-actin) and cytoskeletal regulators. This phenomenon, termed cortical excitability, results from coupled positive and negative feedback loops of cytoskeletal regulators. The nature of these feedback loops, however, remains poorly understood. We assessed the role of the Rho GAP RGA-3/4 in the cortical excitability that accompanies cytokinesis in both frog and starfish. RGA-3/4 localizes to the cytokinetic apparatus, “chases” Rho waves in an F-actin–dependent manner, and when coexpressed with the Rho GEF Ect2, is sufficient to convert the normally quiescent, immature Xenopus oocyte cortex into a dramatically excited state. Experiments and modeling show that changing the ratio of RGA-3/4 to Ect2 produces cortical behaviors ranging from pulses to complex waves of Rho activity. We conclude that RGA-3/4, Ect2, Rho, and F-actin form the core of a versatile circuit that drives a diverse range of cortical behaviors, and we demonstrate that the immature oocyte is a powerful model for characterizing these dynamics.

     
    more » « less
  4. Blackwell, Kim T. (Ed.)
    Transcranial direct current stimulation (tDCS) of the cerebellum has rapidly raised interest but the effects of tDCS on cerebellar neurons remain unclear. Assessing the cellular response to tDCS is challenging because of the uneven, highly stratified cytoarchitecture of the cerebellum, within which cellular morphologies, physiological properties, and function vary largely across several types of neurons. In this study, we combine MRI-based segmentation of the cerebellum and a finite element model of the tDCS-induced electric field (EF) inside the cerebellum to determine the field imposed on the cerebellar neurons throughout the region. We then pair the EF with multicompartment models of the Purkinje cell (PC), deep cerebellar neuron (DCN), and granule cell (GrC) and quantify the acute response of these neurons under various orientations, physiological conditions, and sequences of presynaptic stimuli. We show that cerebellar tDCS significantly modulates the postsynaptic spiking precision of the PC, which is expressed as a change in the spike count and timing in response to presynaptic stimuli. tDCS has modest effects, instead, on the PC tonic firing at rest and on the postsynaptic activity of DCN and GrC. In Purkinje cells, anodal tDCS shortens the repolarization phase following complex spikes (-14.7 ± 6.5% of baseline value, mean ± S.D.; max: -22.7%) and promotes burstiness with longer bursts compared to resting conditions. Cathodal tDCS, instead, promotes irregular spiking by enhancing somatic excitability and significantly prolongs the repolarization after complex spikes compared to baseline (+37.0 ± 28.9%, mean ± S.D.; max: +84.3%). tDCS-induced changes to the repolarization phase and firing pattern exceed 10% of the baseline values in Purkinje cells covering up to 20% of the cerebellar cortex, with the effects being distributed along the EF direction and concentrated in the area under the electrode over the cerebellum. Altogether, the acute effects of tDCS on cerebellum mainly focus on Purkinje cells and modulate the precision of the response to synaptic stimuli, thus having the largest impact when the cerebellar cortex is active. Since the spatiotemporal precision of the PC spiking is critical to learning and coordination, our results suggest cerebellar tDCS as a viable therapeutic option for disorders involving cerebellar hyperactivity such as ataxia. 
    more » « less
  5. Abstract

    Numerical techniques to efficiently model out-of-equilibrium dynamics in interacting quantum many-body systems are key for advancing our capability to harness and understand complex quantum matter. Here we propose a new numerical approach which we refer to as generalized discrete truncated Wigner approximation (GDTWA). It is based on a discrete semi-classical phase space sampling and allows to investigate quantum dynamics in lattice spin systems with arbitraryS ≥ 1/2. We show that the GDTWA can accurately simulate dynamics of large ensembles in arbitrary dimensions. We apply it forS > 1/2 spin-models with dipolar long-range interactions, a scenario arising in recent experiments with magnetic atoms. We show that the method can capture beyond mean-field effects, not only at short times, but it also can correctly reproduce long time quantum-thermalization dynamics. We benchmark the method with exact diagonalization in small systems, with perturbation theory for short times, and with analytical predictions made for models which feature quantum-thermalization at long times. We apply our method to study dynamics in largeS > 1/2 spin-models and compute experimentally accessible observables such as Zeeman level populations, contrast of spin coherence, spin squeezing, and entanglement quantified by single-spin Renyi entropies. We reveal that largeSsystems can feature larger entanglement than correspondingS = 1/2 systems. Our analyses demonstrate that the GDTWA can be a powerful tool for modeling complex spin dynamics in regimes where other state-of-the art numerical methods fail.

     
    more » « less