1
|
Vegué M, Allard A, Desrosiers P. Firing rate distributions in plastic networks of spiking neurons. Netw Neurosci 2025; 9:447-474. [PMID: 40161997 PMCID: PMC11949577 DOI: 10.1162/netn_a_00442] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/19/2024] [Accepted: 01/05/2025] [Indexed: 04/02/2025] Open
Abstract
In recurrent networks of leaky integrate-and-fire neurons, the mean-field theory has been instrumental in capturing the statistical properties of neuronal activity, like firing rate distributions. This theory has been applied to networks with either homogeneous synaptic weights and heterogeneous connections per neuron or vice versa. Our work expands mean-field models to include networks with both types of structural heterogeneity simultaneously, particularly focusing on those with synapses that undergo plastic changes. The model introduces a spike trace for each neuron, a variable that rises with neuron spikes and decays without activity, influenced by a degradation rate r p and the neuron's firing rate ν. When the ratio α = ν/r p is significantly high, this trace effectively estimates the neuron's firing rate, allowing synaptic weights at equilibrium to be determined by the firing rates of connected neurons. This relationship is incorporated into our mean-field formalism, providing exact solutions for firing rate and synaptic weight distributions at equilibrium in the high α regime. However, the model remains accurate within a practical range of degradation rates, as demonstrated through simulations with networks of excitatory and inhibitory neurons. This approach sheds light on how plasticity modulates both activity and structure within neuronal networks, offering insights into their complex behavior.
Collapse
Affiliation(s)
- Marina Vegué
- Departament de Matemàtiques, Universitat Politècnica de Catalunya, Barcelona, Spain
| | - Antoine Allard
- Département de Physique, de Génie Physique et d’Optique, Université Laval, Québec, Canada
- Centre Interdisciplinaire en Modélisation Mathématique, Université Laval, Québec, Canada
| | - Patrick Desrosiers
- Département de Physique, de Génie Physique et d’Optique, Université Laval, Québec, Canada
- Centre Interdisciplinaire en Modélisation Mathématique, Université Laval, Québec, Canada
- CERVO Brain Research Center, Québec, Canada
| |
Collapse
|
2
|
Ocker GK, Litwin-Kumar A, Doiron B. Self-Organization of Microcircuits in Networks of Spiking Neurons with Plastic Synapses. PLoS Comput Biol 2015; 11:e1004458. [PMID: 26291697 PMCID: PMC4546203 DOI: 10.1371/journal.pcbi.1004458] [Citation(s) in RCA: 46] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2014] [Accepted: 07/19/2015] [Indexed: 11/18/2022] Open
Abstract
The synaptic connectivity of cortical networks features an overrepresentation of certain wiring motifs compared to simple random-network models. This structure is shaped, in part, by synaptic plasticity that promotes or suppresses connections between neurons depending on their joint spiking activity. Frequently, theoretical studies focus on how feedforward inputs drive plasticity to create this network structure. We study the complementary scenario of self-organized structure in a recurrent network, with spike timing-dependent plasticity driven by spontaneous dynamics. We develop a self-consistent theory for the evolution of network structure by combining fast spiking covariance with a slow evolution of synaptic weights. Through a finite-size expansion of network dynamics we obtain a low-dimensional set of nonlinear differential equations for the evolution of two-synapse connectivity motifs. With this theory in hand, we explore how the form of the plasticity rule drives the evolution of microcircuits in cortical networks. When potentiation and depression are in approximate balance, synaptic dynamics depend on weighted divergent, convergent, and chain motifs. For additive, Hebbian STDP these motif interactions create instabilities in synaptic dynamics that either promote or suppress the initial network structure. Our work provides a consistent theoretical framework for studying how spiking activity in recurrent networks interacts with synaptic plasticity to determine network structure.
Collapse
Affiliation(s)
- Gabriel Koch Ocker
- Department of Neuroscience, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America
- Center for the Neural Basis of Cognition, University of Pittsburgh and Carnegie Melon University, Pittsburgh, Pennsylvania, United States of America
| | - Ashok Litwin-Kumar
- Center for the Neural Basis of Cognition, University of Pittsburgh and Carnegie Melon University, Pittsburgh, Pennsylvania, United States of America
- Department of Mathematics, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America
- Center for Theoretical Neuroscience, Columbia University, New York, New York, United States of America
| | - Brent Doiron
- Center for the Neural Basis of Cognition, University of Pittsburgh and Carnegie Melon University, Pittsburgh, Pennsylvania, United States of America
- Department of Mathematics, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America
| |
Collapse
|
3
|
Mikkelsen K, Imparato A, Torcini A. Sisyphus effect in pulse-coupled excitatory neural networks with spike-timing-dependent plasticity. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2014; 89:062701. [PMID: 25019808 DOI: 10.1103/physreve.89.062701] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/24/2014] [Indexed: 06/03/2023]
Abstract
The collective dynamics of excitatory pulse-coupled neural networks with spike-timing-dependent plasticity (STDP) is studied. Depending on the model parameters stationary states characterized by high or low synchronization can be observed. In particular, at the transition between these two regimes, persistent irregular low frequency oscillations between strongly and weakly synchronized states are observable, which can be identified as infraslow oscillations with frequencies ≃0.02-0.03 Hz. Their emergence can be explained in terms of the Sisyphus effect, a mechanism caused by a continuous feedback between the evolution of the coherent population activity and of the average synaptic weight. Due to this effect, the synaptic weights have oscillating equilibrium values, which prevents the neuronal population from relaxing into a stationary macroscopic state.
Collapse
Affiliation(s)
- Kaare Mikkelsen
- Department of Physics and Astronomy, University of Aarhus, Ny Munkegade, Building 1520, DK-8000 Aarhus C, Denmark
| | - Alberto Imparato
- Department of Physics and Astronomy, University of Aarhus, Ny Munkegade, Building 1520, DK-8000 Aarhus C, Denmark
| | - Alessandro Torcini
- Department of Physics and Astronomy, University of Aarhus, Ny Munkegade, Building 1520, DK-8000 Aarhus C, Denmark and CNR-Consiglio Nazionale delle Ricerche-Istituto dei Sistemi Complessi, via Madonna del Piano 10, I-50019 Sesto Fiorentino, Italy and INFN Sez. Firenze, via Sansone 1, I-50019 Sesto Fiorentino, Italy
| |
Collapse
|
4
|
Mikkelsen K, Imparato A, Torcini A. Emergence of slow collective oscillations in neural networks with spike-timing dependent plasticity. PHYSICAL REVIEW LETTERS 2013; 110:208101. [PMID: 25167453 DOI: 10.1103/physrevlett.110.208101] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/01/2013] [Revised: 03/19/2013] [Indexed: 06/03/2023]
Abstract
The collective dynamics of excitatory pulse coupled neurons with spike-timing dependent plasticity is studied. The introduction of spike-timing dependent plasticity induces persistent irregular oscillations between strongly and weakly synchronized states, reminiscent of brain activity during slow-wave sleep. We explain the oscillations by a mechanism, the Sisyphus Effect, caused by a continuous feedback between the synaptic adjustments and the coherence in the neural firing. Due to this effect, the synaptic weights have oscillating equilibrium values, and this prevents the system from relaxing into a stationary macroscopic state.
Collapse
Affiliation(s)
- Kaare Mikkelsen
- Department of Physics and Astronomy, University of Aarhus, Ny Munkegade, Building 1520, DK-8000 Aarhus C, Denmark
| | - Alberto Imparato
- Department of Physics and Astronomy, University of Aarhus, Ny Munkegade, Building 1520, DK-8000 Aarhus C, Denmark
| | - Alessandro Torcini
- Department of Physics and Astronomy, University of Aarhus, Ny Munkegade, Building 1520, DK-8000 Aarhus C, Denmark and CNR-Consiglio Nazionale delle Ricerche-Istituto dei Sistemi Complessi, via Madonna del Piano 10, I-50019 Sesto Fiorentino, Italy and INFN Sezione di Firenze, via Sansone, 1-I-50019 Sesto Fiorentino, Italy
| |
Collapse
|
5
|
Hashizume Y, Araki O. Analytical condition for synchrony in a neural network with two periodic inputs. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2013; 87:012713. [PMID: 23410365 DOI: 10.1103/physreve.87.012713] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/18/2012] [Revised: 08/22/2012] [Indexed: 06/01/2023]
Abstract
In this study, we apply a mean-field theory to the neural network model with two periodic inputs in order to clarify the conditions of synchronies. This mean-field theory yields a self-consistent condition for the synchrony and enables us to study the effects of synaptic connections for the behavior of neural networks. Then, we obtain a condition of synaptic connections for the synchrony with the cycle time T. The neurons in neural networks receive sensory inputs and top-down inputs from outside of the network. When the network neurons receive two or more inputs, their synchronization depends on the conditions of inputs. We also analyze this case using the mean-field theory. As a result, we clarify the following points: (i) The stronger synaptic connections enhance the shorter synchrony cycle of neurons. (ii) The cycle of the synchrony becomes longer as the cycle of external inputs becomes longer. (iii) The relationships among synaptic weights, the properties of input trains, and the cycle of synchrony are expressed by one equation, and there are two areas for asynchrony. With regard to the third point, the yielded equation is so simple for calculation that it can easily provide us with feasible and infeasible conditions for synchrony.
Collapse
Affiliation(s)
- Yoichiro Hashizume
- Department of Applied Physics, Tokyo University of Science, Kagurazaka 1-3, Tokyo.
| | | |
Collapse
|
6
|
Chen CC, Jasnow D. Event-driven simulations of a plastic, spiking neural network. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2011; 84:031908. [PMID: 22060404 DOI: 10.1103/physreve.84.031908] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/29/2011] [Indexed: 05/31/2023]
Abstract
We consider a fully connected network of leaky integrate-and-fire neurons with spike-timing-dependent plasticity. The plasticity is controlled by a parameter representing the expected weight of a synapse between neurons that are firing randomly with the same mean frequency. For low values of the plasticity parameter, the activities of the system are dominated by noise, while large values of the plasticity parameter lead to self-sustaining activity in the network. We perform event-driven simulations on finite-size networks with up to 128 neurons to find the stationary synaptic weight conformations for different values of the plasticity parameter. In both the low- and high-activity regimes, the synaptic weights are narrowly distributed around the plasticity parameter value consistent with the predictions of mean-field theory. However, the distribution broadens in the transition region between the two regimes, representing emergent network structures. Using a pseudophysical approach for visualization, we show that the emergent structures are of "path" or "hub" type, observed at different values of the plasticity parameter in the transition region.
Collapse
Affiliation(s)
- Chun-Chung Chen
- Physics Division, National Center for Theoretical Sciences, Hsinchu, Taiwan 300, Republic of China
| | | |
Collapse
|