The interplay between homeostatic synaptic scaling and homeostatic structural plasticity maintains the robust firing rate of neural networks

eLife Assessment

This valuable study combines experiments and modelling to advance our understanding of the nonlinear nature of homeostatic structural plasticity and its interaction with synaptic scaling. The methodology and findings are solid, although additional work is needed to better link models with experiments and support some of the conclusions drawn. This study will be of interest to theoretical and experimental neuroscientists working in homeostatic plasticity.

https://doi.org/10.7554/eLife.88376.3.sa0

Significance of the findings:

Valuable: Findings that have theoretical or practical implications for a subfield

Landmark
Fundamental
Important
Valuable
Useful

Strength of evidence:

Solid: Methods, data and analyses broadly support the claims with only minor weaknesses

Exceptional
Compelling
Convincing
Solid
Incomplete
Inadequate

During the peer-review process the editor and reviewers write an eLife Assessment that summarises the significance of the findings reported in the article (on a scale ranging from landmark to useful) and the strength of the evidence (on a scale ranging from exceptional to inadequate). Learn more about eLife Assessments

Abstract
Introduction
Results
Discussion
Materials and methods
Data availability
References
Article and author information
Metrics

Abstract

Critical network states and neural plasticity enable adaptive behavior in dynamic environments, supporting efficient information processing and experience-dependent learning. Synaptic-weight-based Hebbian plasticity and homeostatic synaptic scaling are key mechanisms that enable memory while stabilizing network dynamics. However, the role of structural plasticity as a homeostatic mechanism remains less consistently reported, particularly under activity inhibition, leading to an incomplete understanding of its functional impact. In this study, we combined live-cell microscopy of eGFP-labeled neurons in mouse organotypic entorhinal-hippocampal tissue cultures (Thy1-eGFP mice of both sexes) with computational modeling to investigate how synapse-number-based structural plasticity responds to activity perturbations and interacts with homeostatic synaptic scaling. Tracking individual dendritic segments, we found that inhibiting excitatory neurotransmission does not monotonically regulate dendritic spine density. Specifically, inhibition of AMPA receptors with 200 nM 2,3-dioxo-6-nitro-7-sulfamoyl-benzo[f]quinoxaline (NBQX) increased spine density, whereas complete AMPA receptor blockade with 50 μM NBQX reduced it. Motivated by these findings, we developed network simulations incorporating a biphasic structural plasticity rule governing activity-dependent synapse formation. These simulations showed that the biphasic rule maintains neural activity homeostasis under stimulation and permits either synapse formation or synapse loss depending on the degree of activity deprivation. Homeostatic synaptic scaling further modulated recurrent connectivity, network activity, and structural plasticity outcomes. It reduced stimulation-triggered synapse loss by downscaling synaptic weights and rescued silencing-induced synapse loss by upscaling recurrent input, thus reactivating silent neurons. The interaction between these mechanisms provides a mechanistic explanation for divergent findings in the literature. In summary, homeostatic synaptic scaling and homeostatic structural plasticity dynamically compete and compensate for each other, ensuring efficient and robust control of firing rate homeostasis.

Introduction

To survive in dynamic environments, animals must quickly respond to familiar cues, such as those signaling food or predators, while remaining alert to novel stimuli. The former reflects experience-dependent learning, where subtle memory cues can trigger strong responses in the corresponding pathway during recall (Cabeza et al., 1997; Nakazawa et al., 2002; Wade-Bohleber et al., 2019; Bone et al., 2020). The latter requires brain networks to maintain activity in a critical state, allowing efficient information transmission via action potentials and neural avalanches (Beggs and Plenz, 2003; Beggs and Plenz, 2004; Petermann et al., 2009; Shew et al., 2009; Hahn et al., 2010; Shew and Plenz, 2013). These processes represent two seemingly opposing but well-coordinated mechanisms: associative learning and firing rate homeostasis.

The concept of firing rate homeostasis has been hovering in theories since the discovery of long-term synaptic potentiation (LTP; Bliss and Lømo, 1973). LTP is a positive feedback mechanism that adjusts synaptic strength among excitatory neurons (‘neurons that fire together, wire together’, Löwel and Singer, 1992), as postulated by Hebb, 1949. While its associative properties allow synapses to preserve memory traces by increasing synaptic weights, it also risks driving network dynamics toward overexcitation or silence (Bliss and Collingridge, 1993; Figure 1Ai, ii, Bi). In vivo studies in rodents demonstrated firing rate homeostasis in the visual cortex and hippocampus, where neural activity is restored at both the single-neuron (Hengen et al., 2016; Ma et al., 2019; Torrado Pacheco et al., 2019) and population level (Slomowitz et al., 2015; Lu et al., 2019), within days of perturbation. This ability to stabilize network dynamics without erasing memories underscores the robustness of the brain’s complex networks (Alon et al., 1999; von Dassow et al., 2000; Carlson and Doyle, 2002; Kitano, 2002; Kitano, 2004; Aoki et al., 2019).

Figure 1

Download asset Open asset

Overview of Hebbian and homeostatic plasticity.

(Ai) In a static network, neural activity is driven by external input. (**Aii**) Hebbian plasticity amplifies network response via a positive feedback mechanism, increasing sensitivity to external input. (**Aiii**) Homeostatic plasticity restores activity toward a setpoint through negative feedback, supporting firing rate homeostasis. (Bi) Hebbian plasticity is synapse-specific, strengthening recurrent connectivity by potentiating the weights of activated synapses. In contrast, synaptic depression can occur either in the same synapse (homosynaptic depression) or in neighboring synapses (heterosynaptic depression) depending on the stimulation protocol. (**Bii**) Homeostatic synaptic scaling is a cell-autonomous process involving proportional up- or down-scaling of all input synaptic weights in response to prolonged changes in postsynaptic activity. (**Biii**) Homeostatic structural plasticity is likewise cell-autonomous and involves compensatory spine loss during chronic excitation. In contrast, chronic inhibition can lead to divergent and sometimes non-homeostatic changes in spine density.

Robustness is a fundamental property of complex biological and engineering systems, enabling adaptation to external perturbations by restoring an original attractor state or shifting to a new one while preserving function (Kitano, 2002; Kitano, 2004). Key features of robustness include negative feedback control, redundancy, and heterogeneity (Alon et al., 1999; von Dassow et al., 2000; Carlson and Doyle, 2002; Kitano, 2002; Kitano, 2004; Aoki et al., 2019). For instance, homeostatic synaptic scaling exemplifies a negative feedback mechanism (Turrigiano et al., 1998; Figure 1Aiii, Bii), where neurons proportionally adjust their synaptic weights to compensate for long-term changes in activity (Turrigiano et al., 1998), membrane potential (Leslie et al., 2001), or intracellular calcium concentration (Slutsky et al., 2004). Redundancy and heterogeneity describe the presence of multiple mechanisms capable of achieving the same outcome, ensuring system robustness against failure. Although redundancy and heterogeneity have been recognized in biological systems in the form of genetic buffering and convergent molecular circuits (Csete and Doyle, 2002; Kitano, 2002; Kitano, 2004), their representation in activity-dependent plasticity has been less explicitly discussed. Mechanisms such as homosynaptic long-term depression (LTD; Lynch et al., 1977) under the Bienenstock–Cooper–Munro (BCM) rule (Bienenstock et al., 1982) or spike-timing-dependent plasticity (STDP; Markram et al., 1997; Bi and Poo, 1998; Gütig et al., 2003; Izhikevich and Desai, 2003), and heterosynaptic LTD via synaptic tagging and capture (Redondo and Morris, 2011) help stabilize firing rates. Inhibitory plasticity also restores the excitation-inhibition balance (Froemke, 2015; Hennequin et al., 2017; Hiratani and Fukai, 2017). These mechanisms, whether modulating synaptic transmission locally or globally, demonstrate functional redundancy and implementational heterogeneity, aligning with homeostatic synaptic scaling to maintain network stability.

However, it remains unclear whether structural plasticity exhibits similar redundancy. Structural plasticity encompasses changes in dendritic spine sizes and numbers, axonal boutons, synapse count, and network connectivity–—all critical for synaptic transmission and memory capacity (Emina and Kropff, 2022). The homeostatic structural plasticity (HSP) model, which posits that synapse formation and elimination are homeostatically regulated (Butz et al., 2009a; Butz et al., 2009b; Butz and van Ooyen, 2013; Butz et al., 2014; Diaz Pier et al., 2016), has expanded theoretical models of network reorganization following diverse perturbations (Gallinaro and Rotter, 2018; Lu et al., 2019; Gallinaro et al., 2022; Lu et al., 2022; Anil et al., 2023). Nonetheless, experimental findings present a more complex picture (Figure 1Biii): synaptic downscaling is often associated with spine loss, while chronic activity deprivation leads to inconsistent and mostly non-homeostatic spine loss (Moulin et al., 2022). This complexity highlights the need for a deeper understanding of activity-dependent structural plasticity, particularly the rules governing spine number changes. A potentially overlooked factor is the degree of inhibition. Given the variability in experimental conditions—ranging from monocular deprivation to lesions and pharmacological treatments—we hypothesized that: (i) The relationship between neural activity and spine numbers is non-monotonic, leading to divergent responses under different levels of inhibition. (ii) Structural responses may be influenced by interactions with other plasticity mechanisms, especially homeostatic synaptic scaling, which is crucial during chronic activity deprivation, as evidenced by both computational models and experimental data.

In this study, we applied two concentrations of the competitive AMPA-receptor antagonist NBQX to gradually suppress neural activity. Using whole-cell patch-clamp recordings and time-lapse imaging, we found that activity deprivation affects dendritic spine numbers and sizes in a non-monotonic manner: a low concentration of NBQX (200 nM) increased spine numbers, while a higher concentration (50 μM) reduced them. These findings guided the development of a spiking neural network model with a biphasic, synapse-number-based HSP rule, which successfully reproduced both homeostatic and non-homeostatic spine changes under different levels of activity perturbations. Additionally, we implemented a monotonic, synaptic-weight-based homeostatic synaptic scaling (HSS) rule, as described in the literature (van Rossum et al., 2000). While largely redundant to the biphasic HSP rule, the HSS rule facilitates the transition from non-homeostatic synapse loss to homeostatic synapse regeneration by adjusting recurrent connectivity in silenced networks. Our results highlight the redundancy and complementarity between HSP and HSS rules, underscoring their combined role in ensuring robust and efficient adaptation of network activity (Yi et al., 2000; Briat et al., 2016; Aoki et al., 2019).

Results

Experimental evidence supports a non-monotonic structural plasticity rule under graded activity blockade

Variability in experimental results regarding structural plasticity under chronic activity deprivation (Moulin et al., 2022) may arise from (i) non-monotonic, activity-dependent regulation of spine numbers, (ii) concurrent changes in synaptic weights and firing rates due to synaptic scaling, or (iii) a combination of both. To systematically address this, we designed computer simulations that incorporate a synapse-number-based structural plasticity rule alone and in combination with synaptic-weight-based synaptic scaling to examine their interactions. In our implementation, these activity-dependent models act as integral feedback mechanisms, both governed by intracellular calcium concentration $C (t)$ and a setpoint value $ϵ$ (Figure 2A), which monitor neural activity and influence synaptic weights or synapse numbers. This approach eliminates the need for manual time-scale alignment. However, structural plasticity models have different variants—some proposing that neural activity influences neurite growth and retraction in a monotonic way (Gallinaro and Rotter, 2018; Lu et al., 2019; Gallinaro et al., 2022) or some in a non-monotonic manner (Butz et al., 2009a; Butz and van Ooyen, 2013; Butz et al., 2014; Diaz Pier et al., 2016). This highlights the need for additional experimental data to refine model selection.

Figure 2

Download asset Open asset

Non-monotonic spine density changes under NBQX inhibition support a biphasic structural plasticity rule.

(A) Integral feedback control framework for homeostatic synaptic scaling and structural plasticity. Both mechanisms use intracellular calcium concentration ( $[C a^{2 +}]_{i}$ , $C (t)$ ) to track neural activity (AP, action potential, $S (t)$ ). Calcium increases with each postsynaptic spike via calcium influx ( $β_{C a}$ ) and decays with a time constant $τ_{C a}$ . Synaptic scaling adjusts synaptic weights $w (t)$ multiplicatively using a scaling factor $ρ$ , based on the deviation of calcium from a setpoint $ϵ$ . Structural plasticity similarly uses $C (t)$ to regulate the growth and retraction of axonal boutons and dendritic spines. Two variants of the structural plasticity rule are illustrated: one with a linear and the other with a non-linear (biphasic) calcium dependency. (B) Example CA1 pyramidal neuron recorded and *post hoc* identified in an entorhinal-hippocampal tissue culture. Scale bar: 500 μm. (**Ci-Cii**) Group data of AMPA receptor-mediated sEPSC from three experimental groups (control: N=18; 200 nM NBQX: N=18; 50 μM NBQX: N=12). (D) Representative Thy1-eGFP culture and time-lapse imaging of dendritic segments in the stratum radiatum (rad.) before and after three-day treatment. Scale bars: 200 μm, 5 μm. (E) Spine density at baseline and after treatment. Values are normalized to the respective baseline per segment. Light-shaded lines show individual trajectories (solid = increase; dashed = decrease). Dark-shaded lines with error bars represent group means and SEM for each group (control: N=19; 200 nM NBQX: N=24; 50 μM NBQX: N=33).

To test whether spine density responds monotonically to deprivation in neural activity, we used two concentrations of the AMPA receptor antagonist NBQX, 200 nM and 50 μM, and studied CA1 pyramidal neurons in mouse organotypic entorhinal-hippocampal tissue cultures (Figure 2B). Whole-cell patch-clamp recording of AMPA receptor-mediated spontaneous excitatory postsynaptic currents (sEPSCs) revealed a significant reduction in mean sEPSC amplitude with 200 nM NBQX and a near-to-complete blockade of excitatory synaptic transmission with 50 μM NBQX (Figure 2Ci–Cii). Both NBQX concentrations were applied to Thy1-eGFP cultures for three days to induce chronic activity inhibition. Individual dendritic segments (side branches of apical dendrites in the stratum radiatum) were imaged before and after treatment using time-lapse microscopy (Figure 2D). The analysis showed that 200 nM NBQX increased spine density ( $p = 0.003$ , Wilcoxon test), while 50 μM NBQX reduced spine density ( $p = 0.008$ , Wilcoxon test), compared to baseline (Figure 2E). No significant changes in spine density were observed in control cultures ( $p = 0.06$ , Wilcoxon test). These results demonstrated a non-monotonic relationship between neural activity and spine numbers under activity deprivation: partial inhibition increased spine numbers, whereas complete inhibition decreased them.

Stabilizing and characterizing a biphasic synapse-number-based structural plasticity rule

Since changes in spine numbers or density are functionally linked to synapse formation, loss, and rewiring, we aimed to establish a synapse-number-based structural plasticity rule in point neurons and study synaptic rewiring in a large, topology-free spiking neuron network to examine the first hypothesis, thereby bypassing the complexities of neuronal morphology (Figure 3A). For comparison, we first simulated the classical linear growth rule (Figure 3Bi). To capture the observed non-monotonic dependency in a simplified manner, we adopted a Gaussian-shaped growth rule for synaptic element numbers with two setpoints. By setting the first setpoint at either $η = 0$ or $η > 0$ , we generated two variants of the Gaussian rule (Figure 3Bii–Biii). We then grew three neural networks in which synapses among excitatory neurons were governed by these three rules, respectively (Figure 3C). As reported in previous studies (Diaz Pier et al., 2016; Gallinaro and Rotter, 2018; Lu et al., 2019; Lu et al., 2022), the linear-rule-guided network developed smoothly into a homogeneous, sparsely connected network (10% of connection probability, equivalent to approximately 1000 excitatory synapses per neuron) with an average firing rate near the target value of $ϵ = 7.9 H z$ (Figure 3Di). The Gaussian rule with a zero setpoint ( $η = 0$ ) also guided the network to a similar equilibrium state (Figure 3Dii). However, the Gaussian rule with two non-zero setpoints ( $η = 0.7$ and $ϵ = 7.9$ ) initially failed to stabilize network development, resulting in network fragmentation and activity silencing (Figure 3Diii).

Figure 3 with 1 supplement see all

Download asset Open asset

Establishing a spiking neural network with three distinct structural plasticity rules.

(A) A network of point neurons was used to study structural plasticity. Dendritic spines are visualized as pink sticks on the soma, and axonal boutons are represented by half-circles–––solid or empty depending on their growth state. An empty half-circle with a dashed line denotes a retracting axon. In this model, intracellular calcium concentration $C (t)$ is linearly correlated with neural firing rate, so neural activity is used throughout this study as a proxy for ‘firing rate homeostasis’, given its physical unit compared to the arbitrary unit (a.u.) of calcium concentration. (**Bi-Biii**) Three structural plasticity growth rules regulate changes in synaptic element numbers. (Bi) Linear rule with a single setpoint ( $ϵ = 7.9$ , orange line). (**Bii**) Gaussian rule with two setpoints, one at zero ( $η = 0$ , grey line; $ϵ = 7.9$ , orange line). (**Biii**) Gaussian rule with two non-zero setpoints ( $η = 0.7$ , grey line; $ϵ = 7.9$ , orange line). Three shades represent 100%, 50%, or 10% of the original growth rate ( $ν$ ), with positive and negative values indicating the rate of synaptic element growth and loss, respectively. (C) Network architecture based on the Brunel network model, consisting of 10000 excitatory (E, blue) and 2500 inhibitory neurons (I, red), receiving external Poisson input (P). Synapses among I-I, I-E, and E-I populations are hard-wired with a fixed 10% probability. E-E synapses are subject to structural plasticity rules. (**Di-Diii**) Temporal evolution of neural activity and network connectivity ( $Γ$ ) during network growth, guided by the three plasticity rules. Unless otherwise stated, the curves and shaded areas represent the mean and standard deviation of firing rates in E and I populations. An equilibrium state ( $Γ = 10 %$ ) was reached under the linear and zero-Gaussian rules (**Di, Dii**), but not with the biphasic Gaussian rule (**Diii**). Firing rate distributions and network connectivity matrices at the marked time points (solid triangles) are shown in Figure 3—figure supplement 1 for **Di-Dii**, and in Figure 4 for **Diii**.

Closer examination of the neural firing rates and network connectivity at the late growth stage of this network (Figure 3Diii) revealed that half of the excitatory neurons were silent (Figure 4A), and the network exhibited inhomogeneous connectivity, with nearly half of the neurons isolated (Figure 4Bi–Bii). Correlating the firing rates with synapse numbers of individual excitatory neurons confirmed that the silent neurons were isolated, while neurons with firing rates close to the target possessed approximately 1000 synapses, similar to the linear rule (Figure 4C). These results align with the properties of the Gaussian growth rules: one setpoint is stable ( $ϵ = 7.9$ ), while the second setpoint can be zero if the synapse numbers remain unchanged after silencing ( $η = 0$ ) or non-zero if disconnection occurs after silencing ( $η > 0$ ). We refer to the Gaussian rule with one zero-valued setpoint as the zero Gaussian rule, and the variant with two non-zero setpoints is termed the biphasic Gaussian rule. After applying appropriate damping current, the biphasic Gaussian rule produced network statistics comparable to those observed with the linear and zero Gaussian rules (Figure 4D, E and F), thereby allowing a systematic comparison of their effects following network development.

Figure 4

Download asset Open asset

Silent neurons remain isolated under the biphasic Gaussian rule.

(A) Histogram of firing rates for excitatory and inhibitory neurons sampled at the time point indicated in Figure 3Diii. Nearly half of the excitatory neurons remained silent. The blue vertical line marks the mean firing rate of non-silent excitatory neurons, while the orange dashed line indicates the target firing rate ( $ϵ = 7.9 H z$ ). (**Bi-Bii**) Network connectivity matrix and distribution of synapse numbers across individual excitatory neurons. (C) Correlation heatmap showing the relationship between neural activity and synapse number per excitatory neuron. Silent neurons failed to form synapses, whereas active neurons fired around the target rate and formed approximately 1000 synapses with other active excitatory neurons. (D) Network architecture after introducing a damping facilitating current ( $I_{f a c i l i t a t i n g}$ ) to promote network development. (E) Temporal dynamics of firing rate (mean and standard deviation) and network connectivity after current injection. The inset displays firing rate distributions for excitatory and inhibitory neurons at the time point indicated by the solid triangle. The facilitating current started at 750 pA and decayed linearly to zero over 4000 s. (F) Different initial amplitudes of the facilitating current led to distinct levels of network connectivities. An initial value of 750 pA was used consistently throughout the study.

The biphasic Gaussian rule reconciles homeostatic and non-homeostatic structural changes following activity perturbation

To examine the first hypothesis—that synapse-number-based structural plasticity responds non-monotonically to activity perturbations—we systematically modulated the input strength to a subgroup of excitatory neurons (10% of the entire excitatory population) after network development. The input strength varied from 0% to 200% fold of the original intensity (FOI, Figure 5A), mimicking both input deprivation and stimulation. Figure 5Bi–Biii and C illustrates three example protocols demonstrating neural activity and network connectivity under the biphasic Gaussian rule. As expected, stimulation increased the activity of the targeted subpopulation (S) (left panel in Bii) and triggered a reduction in connectivity both among the stimulated neurons (S-S) and between the stimulated and non-stimulated excitatory neurons (S-E, left panel in Biii). The connectivity matrices in panel C display the final network connectivity. Ultimately, homeostatic disconnection restored the activity levels of the stimulated neurons to the setpoint value (orange line in Bii). In contrast, the outcome of activity deprivation depended on its severity. Both mild deprivation and silencing reduced the neural activity of the affected subpopulation (middle and right panels in Bii). However, only mild deprivation allowed the network to restore activity to the setpoint through a homeostatic increase in network connectivity (middle panels in Biii and C). In the case of complete silencing, the affected neurons disconnected from the network (right panels in Biii and C).

Figure 5

Download asset Open asset

Divergent regulation of network connectivity under stimulation and deprivation via three structural plasticity rules.

(A) A subpopulation (S) comprising 10% excitatory neurons (E) was subject to activity perturbation. All E-E, S-S, and E-S synapses were governed by the biphasic Gaussian rule. (Bi) Activity perturbation protocol. Three levels of input intensity, expressed as a percentage of the original Poisson input (fold of intensity, FOI), were used: *stimulation* (110% FOI), *mild deprivation* (95% FOI), and *silencing* (0% FOI). (**Bii**) Temporal dynamics of mean firing rates and standard deviation for the S, E, and inhibitory population (I) under each protocol. (**Biii**) Temporal evolution of overall network connectivity ( $Γ_{o v e r a l l}$ ) and subpopulation-specific connectivity ( $Γ_{s u b s}$ ) under each condition. Due to the symmetric application of plasticity rules to spines and boutons, the connection probabilities from E to S and S to E are equivalent. Therefore, only S-E traces are shown unless otherwise noted. (C) Final network connectivity matrices following stimulation, mild deprivation, or silencing. (D) Average number of incoming synapses to neurons in S subpopulation across a range of FOIs. Empty green circles indicate values from networks exposed to extreme stimulation or inhibition, where both activity and network connectivity were unstable. (E) Examples of two neurons receiving identical external input but with different calcium decay time constants ( $τ_{C a}$ ). The upper panel shows membrane potential traces; the middle panel displays spike trains, and the lower panel illustrates integrated calcium concentration over time. (F) Connectivity trajectories of silenced subnetwork under three different conditions. The dashed circle highlights a region shown at higher magnification.

A systematic analysis of incoming synapse numbers confirmed the biphasic dependency of the biphasic Gaussian rule (dark yellow curve in Figure 5D), in contrast to the monotonic behavior of the linear rule (green curve). The zero Gaussian rule, which features a zero-valued setpoint ( $η = 0$ ), exhibited an intermediate profile: external stimulation and mild deprivation triggered homeostatic reductions or increases in synapse numbers, respectively, but strong deprivation left the network silent yet structurally intact (light yellow curve). Among the three rules, the biphasic Gaussian rule best captured the homeostatic characteristics of the linear rule while also allowing for non-homeostatic, silencing-induced spine loss, consistent with our experimental findings and previous studies (Keck et al., 2008; Vuksic et al., 2011; Vlachos et al., 2012; Vlachos et al., 2013; Bissen et al., 2021).

Lesion- or denervation-induced plasticity is of particular interest in the context of activity-dependent structural plasticity and holds clinical relevance (Desmurget et al., 2007; Di Pino et al., 2014; Sampaio-Baptista et al., 2018). The biphasic Gaussian rule provides a useful model for studying denervation-induced degeneration and subsequent regeneration. However, we observed a small, biologically unrealistic increase in connectivity immediately following silencing, possibly caused by a high growth rate or lingering effects of prior activity on intracellular calcium concentrations. In Figure 5E, we show that two neurons receiving identical external input displayed identical membrane potential dynamics and spike trains, but a difference in the calcium time constant ( $τ_{C a}$ ) led to a divergence in accumulated calcium concentration. Since synapse-number dynamics in our model depend on intracellular calcium concentration—used here as a proxy for firing rate—a shorter time constant ( $τ_{C a} = 1 s$ , orange curve in Figure 5F) carried less activity history, resulting in a smoother reduction in network connectivity. A similar effect was achieved by reducing the structural plasticity growth rate to 10% of the original value (light blue dashed curve, Figure 5F). Our computer simulations demonstrated that the time scales of calcium dynamics and structural plasticity interact. The decay time constant of intracellular calcium concentration is difficult to estimate, given the diversity of calcium signals and their temporal and spatial compartmentalization. Moreover, structural plasticity typically operates over longer timescales—minutes to days. Therefore, we used a slow growth rate and the original calcium time constant in all subsequent simulations. The faster growth rate employed in Figure 6D–F was used solely to accelerate spine turnover and reduce simulation time to simulate long-term effects.

Figure 6

Download asset Open asset

Activity perturbation and recurrent connectivity jointly shape the evolution of network connectivity.

(A) Structural plasticity rules applied to both axonal boutons (A, light brown curve) and dendritic spines (D, pink curve), using identical growth functions with $η_{A} = η_{D}$ (upper inset). Variations with asymmetric parameters (different $η$ values could be used for axonal and dendritic elements, such as $η_{A} > η_{D}$ or $η_{A} < η_{D}$ ) are conceptually possible. A silencing protocol was implemented using a fold of original input intensity (FOI). (B) Temporal dynamics of neural activity (mean and standard deviation) for the silenced subpopulation (S), excitatory neurons (E), and inhibitory neurons (I) under two different conditions. (C) Final network connectivity matrices at the end of the silencing protocol under the same two conditions shown in B. (D) Experimental protocols used to assess the combined effects of recurrent connectivity and external stimulation. External stimulation was introduced to the deprived subnetwork under conditions of varying recurrent connectivity. In Protocol 2 (Prtcl. 2), the external stimulation intensity was doubled relative to Protocol 1 (Prtcl. 1). The same structural plasticity rule was applied to both axonal and dendritic elements ( $η_{A} = η_{D}$ ). (E) Temporal trajectories of neural activity (mean and standard deviation) following silencing and external stimulation. The network in the left panel was simulated with a structural growth rate 10 times faster than that in the middle and right panel, resulting in lower recurrent connectivity. (F) Corresponding evolution of overall connectivity ( $Γ_{o v e r a l l}$ , top) and subnetwork-specific connectivity ( $Γ_{s u b s}$ , bottom) for excitatory-excitatory (E–E), subpopulation (S–S), and cross-population (S–E) connections.

Recurrent connectivity and activity perturbation shape synapse-number dynamics during structural recovery

While silencing-induced disconnection captures the characteristic spine loss observed during input deprivation, complete neural isolation is rarely observed in vivo. Instead, spine recovery often follows an initial loss, consistent with experimental observations and previous modeling studies (Keck et al., 2008; Vlachos et al., 2012; Bissen et al., 2021). Prior simulations using Gaussian rules have also reported such ‘physiological’ recovery, assuming distinct growth rules for axonal boutons and dendritic spines (Butz and van Ooyen, 2013; Butz et al., 2014). However, when we implemented different $η$ values for axonal bouton and dendritic spine growth in a topology-free network (without distance-dependent connectivity), synapse numbers did not recover after silencing (Figure 6A and B). Instead, we observed only asymmetric connectivity between input and output synapses of the deprived subpopulation (Figure 6C). Closer examination of the network revealed the critical role of recurrent input. In networks with distance-dependent connectivity, neurons at the edge of the deprived region remained partially active due to input from non-deprived neighbors, whereas neurons at the center became fully isolated—even though external input was uniformly removed.

This observation prompted further investigation into the role of recurrent connectivity in synaptic recovery. In our topology-free network, we systematically modulated both external input and internal connectivity. After applying the silencing protocol, we introduced external stimulation to the deprived subnetwork under conditions of varying recurrent connectivity. In Figure 6E and F (Protocol 1), the network shown in the left panel was simulated with a ten-fold faster timescale than the middle panel, resulting in lower recurrent connectivity at the time stimulation was applied. To isolate the effect of recurrent connectivity, we avoided varying stimulation parameters—such as intensity, duration, interval, or session number—which have been addressed in prior work (Lu et al., 2019; Manos et al., 2021; Anil et al., 2023). Instead, we kept stimulation constant and compared two intensities as a proof of concept. The results demonstrated that network recovery depends on the product of external input and recurrent connectivity. In the left panel of Figure 6E, the deprived subnetwork remained silent due to low recurrent connectivity, despite receiving stimulation. In contrast, higher recurrent connectivity in the middle panel allowed the same stimulation to reactivate the subnetwork. Connectivity traces in Figure 6F confirmed that residual recurrent input amplified the stimulation, enabling synapse regeneration and rewiring. However, when the external stimulation was too strong (protocol 2), the amplification effect of recurrent connectivity pushed the system into a different regime: from silencing-induced degeneration to over-excitation-induced homeostatic degeneration (Figure 6E and F, right panels). These findings suggested that recurrent connectivity serves as a critical amplifier of external input, shaping structural recovery or degeneration following activity perturbations.

Synaptic scaling modulates structural plasticity and reshapes connectivity during activity perturbation

So far, we have confirmed that the biphasic Gaussian rule reconciled both homeostatic and non-homeostatic regulation of spine numbers under activity deprivation, as proposed in the first hypothesis, and that recurrent connectivity plays a critical role in shaping the expression of this rule. In line with this framework, homeostatic synaptic scaling may dynamically modulate neural activity and reshape recurrent connectivity by adjusting functional synaptic transmission—forming the basis of the second hypothesis. Specifically, downscaling of excitatory synapses may reduce heightened neural activity, potentially masking structural changes if it acts rapidly and effectively. Conversely, homeostatic strengthening of excitatory synapses under silencing may restore activity, thereby shifting structural plasticity from non-homeostatic synapse loss to homeostatic rewiring—similar to the effects observed with external stimulation or distance-dependent connectivity.

To test this hypothesis, we implemented a monotonic, activity-dependent synaptic scaling rule, guided by the same intracellular calcium concentration as the biphasic Gaussian structural plasticity rule. To ensure consistent baseline connectivity and avoid developmental confounds, we first allowed the network to grow using the biphasic Gaussian structural rule until it reached equilibrium. Synaptic scaling was then activated at the onset of the activity perturbation protocol (Figure 7A, Figure 7—figure supplement 1A). To disentangle the contributions of structural plasticity and synaptic scaling, we analyzed two connectivity measures: (1) Structural Connectivity ( $Γ_{s t r u c .}$ ), determined solely by synapse number, and (2) Effective Connectivity ( $Γ_{e f f e c .}$ ), defined as the product of synapse numbers and synaptic weights. Without synaptic scaling, synaptic weights remained uniform, so structural and effective connectivity were identical (Figure 7C, left panels; Figure 7—figure supplement 1Bi). Intuitively, when neural activity increased through stimulation, both structural plasticity and synaptic scaling acted redundantly in a homeostatic manner. In this case, synaptic downscaling attenuated homeostatic disconnection, making spine loss less pronounced (Figure 7—figure supplement 1Bii).

Figure 7 with 1 supplement see all

Download asset Open asset

Homeostatic synaptic scaling (HSS) shapes effective connectivity and interacts with homeostatic structural plasticity (HSP).

(A) Protocol of silencing (0% FOI) was applied to a subpopulation of excitatory neurons (S) and synaptic scaling was activated with three different strengths: $ρ = 0$ (no scaling), $ρ = 0.01$ (weak scaling), and $ρ = 0.02$ (strong scaling). (**B-C**) Time courses of network activity (mean and standard deviation) and network connectivity. $Γ_{s t r u c .}$ denotes structural connectivity and $Γ_{e f f e c .}$ refers to effective connectivity. (**Di-Dii**) Structural and effective connectivity matrices of the entire network at time point $t 2$ . (**Ei-Eii**) Structural and effective connectivity matrices at time point $t 3$ . (F) Raster plots of spiking activity for 100 representative neurons at the time points $t 1$ , $t 2$ , and $t 3$ , including inhibitory (I, red), excitatory (E, blue), and silenced excitatory neurons (S, grey).

In contrast, activity deprivation revealed a more complex interaction. With weak synaptic scaling (Figure 7B and C, $ρ = 0.01$ ), effective connectivity to the deprived neurons (E-S) increased gradually but had minimal impact on firing rates and failed to restore synapse numbers and network connectivity. When scaling strength was increased ( $ρ = 0.02$ ; Figure 7C, right panel), effective E-S connectivity increased more substantially, reactivating previously silent neurons and initiating synapse regeneration and rewiring. We further examined connectivity matrices at two time points: $t 2$ (Figure 7Di, ii) and $t 3$ (Figure 7Ei, ii). Notably, at $t 2$ , input synapse numbers to the deprived subgroup decreased, but the remaining synapses receiving input from non-deprived neurons increased in synaptic weight. This led to the reactivation of silent neurons (Figure 7B, right panels), and by $t 3$ , input synapse numbers began to recover. These results confirm that monotonic synaptic scaling can oppose the biphasic structural plasticity rule under strong deprivation, facilitating a transition from synapse loss to activity-dependent rewiring. Intriguingly, with sufficient scaling strength, both synapse numbers and average neural activity were restored after silencing, and the network’s temporal dynamics were altered. The previously deprived subpopulation showed high synchronization after recovery (Figure 7F), suggesting mechanisms relevant to pathological states such as post-traumatic epilepsy. Overall, our simulations suggested that biphasic HSP and monotonic homeostatic synaptic scaling play complementary and partially redundant roles in maintaining firing rate homeostasis. Their interaction provides a biologically plausible mechanism for reconciling stable activity control with flexible network reorganization.

Spine size dynamics reveal an interaction between structural plasticity and synaptic scaling

The altered dynamics in the deprived neural population after recovery may result from dominant input originating from non-deprived excitatory neurons. To test this, we analyzed the firing rates and input synaptic weights of representative neurons from both the deprived (grey) and non-deprived (blue) excitatory populations (Figure 8A). This analysis revealed a marked increase in synaptic weights from non-deprived to deprived neurons (Figure 8A, lower panel), further supported by a global shift in synaptic weight distribution (Figure 8B). Additionally, the distribution analysis indicated changes in input synapse numbers after silencing (see values above Figure 8B), prompting us to assess both synaptic weights and input synapse numbers for a representative neuron (solid line) and its peers (dashed line in Figure 8C) to gain further insight into structural plasticity. Consistent with experimental findings, input synapse numbers of deprived neurons decreased after silencing but gradually recovered. Interestingly, recovery patterns diverged: synapses from reactivated neurons exceeded their initial numbers (upper panel), while those from non-deprived neurons (middle panel) and total synapse counts (lower panel) returned to slightly lower levels than baseline. These results demonstrated that an appropriate combination of homeostatic synaptic scaling and structural plasticity can replicate the process of spine loss and recovery following activity deprivation. Importantly, the model captures a commonly observed outcome of chronic activity silencing: fewer but stronger synapses, corresponding to fewer but larger spines.

Figure 8 with 1 supplement see all

Download asset Open asset

Spine size analysis suggests interactions between functional and structural plasticity.

(A) Activity trajectories and average incoming excitatory synaptic weights of two representative excitatory neurons: one from the deprived subpopulation (N1, grey) and one from the non-deprived group (N2, blue). The gray dashed line in the upper panel indicates the target activity setpoint. Shaded areas represent the standard deviation of incoming synaptic weights. Synapses were categorized based on their source: from the deprived subpopulation (S, grey) or from non-deprived excitatory neurons (E, blue). Triangles in the middle panel indicate the time points used for synaptic weight distribution analyzed in B. (B) Distributions of synaptic weights for neuron N1 *before* and *after* activity deprivation. Numbers above the plots indicate the total synapse count per source type at each time point. (C) Normalized synapse numbers over time of neuron N1 (solid line) and the averaged across all deprived neurons (dashed line). (Di-Diii) Baseline spine size distributions for three groups: control (N=489), 200 nM-NBQX treated (N=736), and 50 μM-NBQX treated (N=675). (E) Cumulative distribution functions (cdf) of spine sizes before and after the three-day treatment. The inset shows average normalized spine sizes across conditions. (F) Normalized spine size changes grouped by initial spine size. The $x$ -axis denotes the upper limit of each bin. Data points above the dashed line represent spine enlargement; points below indicate shrinkage. Error bars in the inset of E and in F represent SEM for each group.

To test whether this principle also holds in our NBQX experiments—and to assess potential dose dependency—we analyzed spine sizes before and after chronic NBQX treatment in mouse entorhinal-hippocampal slice cultures. As expected, initial spine size distributions in all three groups showed a long-tail pattern, with a predominance of small spines (Figure 8Di–Diii), consistent with previous reports (Arellano et al., 2007; Buzsáki and Mizuseki, 2014). After 3 days of treatment, the distribution of spine sizes changes ( $Δ s p i n e s i z e$ ) approximated a normal distribution (Figure 8—figure supplement 1Bi, iii). At the population level (Figure 8E), significant increases in raw spine sizes were observed following 200 nM NBQX treatment ( $p < 0.001$ , Wilcoxon test; 99.9% CI = [9.09, 25.8], LMM), while 50 μM NBQX treatment led to significant reductions ( $p < 0.001$ , Wilcoxon test; 99.9% CI = [9.09, 25.8], LMM), indicating a non-monotonic, activity-dependent response. When size changes were normalized by their initial spine sizes, a general increase was observed across all three groups (Figure 8E inset), likely driven by the growth of many small spines. Sorting spine size changes by initial size (Figure 8F) revealed that in control segments (black curve), small spines generally enlarged while large spines tended to shrink. Under partial inhibition with 200 nM NBQX (orange curve), this relationship was shifted upward, promoting spine enlargement regardless of initial size. However, complete inhibition with 50 μM NBQX (light blue curve), not only reversed the typical growth of small spines but also altered the fate of a subset of large spines—favoring further enlargement instead of the typical shrinkage observed in medium-sized spines.

These findings suggested a homeostatic adjustment of spine sizes in response to partial and complete inhibition. Notably, under complete inhibition, the spine size changes were strongly dependent on initial spine sizes, paralleling both experimental and simulation results in which complete activity blockade resulted in fewer but enlarged spines. However, this pattern deviates from the expected synaptic weight changes predicted by the synaptic scaling rule. This discrepancy will be further examined in the discussion, where we integrate electrophysiological and imaging data for a more comprehensive interpretation.

Hybrid combinations of synaptic scaling and structural plasticity reconcile spine density inconsistencies

Having obtained both spine density and size data for individual segments, we explored whether a unifying principle might link changes in spine size and density. We quantified the average spine size change ( $\bar{Δ} s p i n e s i z e$ ) for each segment and correlated these values with baseline spine densities (Figure 9A) and with changes in spine density ( $Δ s p i n e d e n s i t y$ , Figure 9B). A significant positive correlation between spine size changes and baseline spine densities was observed in the 200nM NBQX treatment group (Figure 9A, $r (24) = 0.65$ , $p = 0.00064$ , Pearson’s correlation; Figure 9—figure supplement 1A), but not in the control or 50 μM NBQX groups. When we correlated spine size changes with spine density changes across all segments, weak inhibition (200 nM NBQX) typically resulted in coordinated increases in both size and density (see orange dots clustered in the upper right quadrant in Figure 9B). This suggested that spine-dense segments may be more responsive to weak inhibition, reflecting a synergistic interaction between structural and functional plasticity. In contrast, under strong inhibition (with 50 μM NBQX), the size-density relationships were highly variable (light blue dots in Figure 9B), consistent with the diverse spine remodeling outcomes observed. Although our imaging data did not reveal a universal rule, they suggested that responses to activity deprivation are highly segment-specific. This led us to hypothesize that different combinations of homeostatic synaptic scaling and structural plasticity operate at the dendritic segment level, or may even vary between cell types (Figure 9C). Such diversity could help explain inconsistent findings across brain regions or between apical and basal dendrites of the same neuron (Lenz et al., 2023b).

Figure 9 with 3 supplements see all

Download asset Open asset

Systematic study of the interaction between homeostatic synaptic scaling and structural plasticity in response to activity silencing.

(A) Correlation between each dendritic segment’s average change in spine size and its initial spine density. Marker size reflects the net change in spine density over the three-day period. (B) Correlation between each dendritic segment’s average changes in spine size and its change in spine density. (C) Conceptual hypothesis: Different combinations of homeostatic synaptic scaling (HSS) and homeostatic structural plasticity (HSP) may form a spectrum, allowing individual neuron types or even distinct dendritic segments within the same neuron to be governed by unique rule combinations. This framework could account for the variability observed in experimental studies of structural plasticity. Sample neuron reconstructions are based on CA1 pyramidal cells (left two) and dentate gyrus granule cells (right two) recorded in our lab. (D) Simulation protocol used to explore the interaction between homeostatic synaptic scaling and structural plasticity. The growth rate of the HSP rule ( $ν$ ) and the strength of the HSS rule ( $ρ$ ) were systematically varied. (**Ei-Eii**) Example time course of structural connection probability and neural activity in the deprived subpopulation (S) under selected combinations of ( $ν$ ) and ( $ρ$ ). (**Fi-Fii**) Quantification of deviations in connectivity and firing rate under different parameter combinations. Deviations were computed as the area between the actual trajectory and either the equilibrium connection probability (10%) or the target firing rate ( $ϵ$ ) over the silencing period. Results represent the mean from 11 independent trials.

To test this hypothesis, we simulated the same silencing protocol while varying the structural plasticity growth rate ( $ν$ ) and synaptic scaling strength ( $ρ$ ) (Figure 9D). Example trajectories of neural activity and structural connectivity in the deprived subpopulation under different conditions are shown in Figure 9Ei, ii. To quantify rewiring efficacy, we calculated the cumulative deviation of structural connectivity from the equilibrium value (10%) in the deprived subnetwork over the silencing period (Figure 9Fi). Cold colors indicate failed connectivity restoration; warm colors indicate successful or excessive rewiring. Observationally, slower growth rates preserved more residual synapses, allowing even weak synaptic scaling to restore connectivity. In contrast, faster growth resulted in more extensive spine loss, requiring stronger synaptic scaling for recovery.

We also tracked deviations of firing rates from target (7.9 Hz) over the same period (Figure 9Fii). Parameter combinations that failed to restore network connectivity (Figure 7Fi) showed persistent silencing (dark blue in Figure 7F). However, successful synapse regeneration (warm colors in Figure 7Fi) did not always ensure optimal firing rate recovery. Best results were achieved with moderate synapse regeneration (light pinkish colors in panel Fii), whereas strong structural growth combined with strong synaptic scaling could induce excessive activity (bright colors) or delayed rebound activity (dark blue colors). In conclusion, our systematic simulations demonstrated that calcium-based homeostatic synaptic scaling and structural plasticity can interact in a dynamic, context-dependent manner. Their interplay—shaped by timing and strength—offers a plausible explanation for the diverse spine remodeling outcomes observed experimentally (Moulin et al., 2022).

Discussion

The mammalian brain is a complex system composed of billions of neurons and non-neuron cells, organized into distinct but interconnected layers, regions, and circuits (Herculano-Houzel, 2012). It is remarkably robust—information overload does not result in catastrophic forgetting (French, 1999), nor does sensory deprivation disrupt overall network dynamics (Hengen et al., 2016; Ma et al., 2019; Torrado Pacheco et al., 2019). Hebbian plasticity and homeostatic synaptic scaling (HSS), two well-characterized forms of functional plasticity, are central to preserving this stability. In contrast, structural plasticity has received comparatively less attention, owing to inconclusive experimental findings (Moulin et al., 2022) and a lack of comprehensive computational models (Butz et al., 2009b). Under conditions of activity deprivation, both homeostatic and non-homeostatic changes in spine numbers have been reported (Moulin et al., 2022), challenging the explanatory power of the HSP framework. In this study, we confirmed a biphasic HSP rule with experiments and examined its interaction with a monotonic HSS mechanism during activity deprivation using computer simulations. Our time-lapse imaging experiments confirmed the non-monotonic, activity-dependent regulation of spine numbers, supporting the proposed biphasic HSP model. Complementary simulations further demonstrated that this rule reconciles opposing structural outcomes—spine gain and loss—and interacts with HSS to jointly regulate synapse number and synaptic strength. We argue that the redundancy and heterogeneity between HSP and HSS mechanisms are critical for maintaining firing rate homeostasis and ensuring long-term network stability.

Our research significantly advances the understanding of structural plasticity by providing direct evidence of its activity-dependent and differential regulation. While structural plasticity is frequently linked to cortical reorganization following injury, amputation, or disease-related cognitive decline, its theoretical modeling has lagged behind that of functional plasticity—largely due to challenges of capturing high-dimensional morphological dynamics. Among existing mathematical frameworks, some have focused on dendritic tree development (Vormberg et al., 2017; Ferreira Castro et al., 2020; Cuntz et al., 2021), while others emphasized synapse formation and rewiring as key distinctions from functional plasticity models. However, these early models generally adopted a synapse-specific, Hebbian-like approach, relying on use-dependent growth and pruning of individual synapses (Bourjaily and Miller, 2011; Fauth et al., 2015; Navlakha et al., 2015; Spiess et al., 2016; Wang et al., 2021), which often proved insufficient to stabilize network dynamics. In contrast, Butz and van Ooyen introduced a cell-autonomous HSP rule, grounded in the earlier work of van Ooyen and van Pelt’s work (van Ooyen and van Pelt, 1994) and conceptually aligned with ideas proposed by Dammasch, 1989. This model uses a target firing rate—or intracellular calcium concentration—as a setpoint to regulate spine and bouton formation or elimination (van Ooyen, 2011; Butz and van Ooyen, 2013). As a synapse-number-based mechanism, it serves as a structural analog to synaptic scaling for maintaining firing rate homeostasis. At the same time, it retains the associative nature of correlation-based models (Gallinaro and Rotter, 2018; Lu et al., 2019; Gallinaro et al., 2022) and promotes self-organized criticality in developing neural networks (van Ooyen and Butz-Ostendorf, 2019). Whether implemented in linear or Gaussian form, this rule has long been proposed as a means of integrating Hebbian and homeostatic plasticity, capturing a range of phenomena observed during development and pathological conditions. Our findings support the biological plausibility of a biphasic variant of this rule, substantiated by both experimental and computational evidence.

We used gradual activity blockade combined with time-lapse imaging to test the hypothesis that spine numbers follow a biphasic relationship with neural activity, thereby addressing longstanding inconsistencies in spine density findings across the literature (Moulin et al., 2022). Many experimental paradigms—such as lesion-based models or pharmacological inhibition—introduce confounding factors, including neuroinflammatory cytokine production (Vlachos et al., 2012; Keck et al., 2013; Vlachos et al., 2013; Barnes et al., 2017), which can independently influence plasticity mechanisms (Turrigiano, 2006; Becker et al., 2013; Kleidonas and Vlachos, 2021). Non-traumatic deprivation models like dark rearing (Wallace and Bear, 2004) or whisker trimming (Trachtenberg et al., 2002) result in non-uniform input deprivation across the dendritic arbor, complicating interpretations of spine density changes across compartments. Even precise manipulation, such as pre-synaptic blockade (Quinn et al., 2019) or glutamate uncaging near individual spines (Tong et al., 2021), may be limited by their homogeneity of effects or lack of circuit specificity. Importantly, even studies using drugs that uniformly target postsynaptic neurons have reported inconsistent effects on spine density (Bacci et al., 2001; Thiagarajan et al., 2005; Wierenga et al., 2005; Fishbein and Segal, 2007; Mitra et al., 2012; Hobbiss et al., 2018), suggesting that the intensity of deprivation may be a key factor. Variables such as drug type (non-competitive vs. competitive), experimental context (in vivo vs. in vitro), and circuit integrity strongly modulate the degree of inhibition experienced by target neurons. Our previous work in CA1 pyramidal neurons revealed pathway-specific differences in structural, cellular, and molecular responses to synaptic deprivation, highlighting the importance of input origin in shaping plasticity (Lenz et al., 2023a). Based on these observations, we hypothesized a biphasic activity-dependent rule: partial activity suppression promotes spine growth, whereas complete inhibition leads to spine loss. Indeed, our time-lapse imaging confirmed this prediction. Spine density increased following partial AMPA receptor blockade (200 nM NBQX) and decreased under complete inhibition (50 μM NBQX). However, we also observed complex changes in spine size under both conditions, suggesting an interaction between functional and structural plasticity, which we discuss in subsequent sections.

The biphasic structural plasticity model, implemented as a Gaussian-shaped growth curve, successfully reconciled the diverse, activity-dependent changes in spine density observed in our experimental data and simulations. With two non-zero setpoints ( $η > 0$ ), the model predicted a range of synapse number outcomes—including increases, no changes, or reductions—depending on the degree of activity perturbation. Notably, silencing-induced synapse loss, which cannot be captured by linear (monotonic) models, falls within the predictive scope of the biphasic rule. A key emergent feature of this model is its ability to induce phase transition through modulation of external input or intrinsic connectivity. For example, appropriately timed and scaled stimulation can reverse silencing-induced synapse loss, triggering homeostatic rewiring. This mechanism may have clinical relevance, potentially informing non-invasive brain stimulation strategies aimed at mitigating stroke-induced network disconnection (Bai et al., 2022). In parallel, HSS contributes to synaptic recovery by adjusting synaptic weights, consistent with empirical observations of both synapse numbers and synaptic strength. Similar forms of rewiring have been described in previous work by Butz and van Ooyen using the same Gaussian HSP framework (Butz et al., 2009a; Butz and van Ooyen, 2013; Butz et al., 2014). Their studies incorporated distance-dependent connectivity and emphasized the need for distinct growth rules for axonal boutons and dendritic spines. Our ability to replicate such rewiring in a topology-free, homogeneous network supports and extends their findings. It underscores the importance of intrinsic connectivity in governing synapse loss and recovery, even in the absence of predefined spatial constraints. Importantly, our co-simulation of structural and functional plasticity revealed a synergistic interaction between these mechanisms—an aspect often overlooked in computational studies, yet critical for interpreting biological data. Additionally, inhibitory plasticity may further contribute to the physiological regulation of HSP, as suggested by numerical studies (Sinha et al., 2015); however, this remains beyond the scope of the present work and warrants future investigation.

Nevertheless, our model suggests a potential mechanism by which injury-related seizures could emerge following otherwise adaptive recovery. Stroke is a leading cause of seizures and the development of epilepsy in the elderly population (Beghi et al., 2011; Guekht and Bornstein, 2012; Galovic et al., 2021) with higher seizures incidence reported after hemorrhagic compared to ischemic stroke (Bladin et al., 2000). According to our model, cell death and resulting input deprivation initiate a cascade of plastic reorganization aimed at restoring network activity. While this reorganization supports functional recovery, it may also create a substrate for epileptiform activity. In turn, seizure activity can exacerbate neural damage, creating a vicious cycle of degeneration and maladaptive plasticity. A closer examination of the epidemiological data reveals a more nuanced relationship between plasticity and post-stroke seizures. Seizures occurring after stroke are typically classified by onset time: acute symptomatic seizures (or early seizures) occur within seven days of infarct and are linked to metabolic or toxic effects of stroke, whereas remote symptomatic seizures (or late seizures) arise after one week and carry a higher risk of developing into chronic epilepsy (Galovic et al., 2021). Although structural and functional plasticity occur within minutes to days, the timing of late seizures aligns more closely with prolonged, plasticity-driven network remodeling. Additional risk factors—such as younger age (<65 years) (Phan et al., 2022) and cortical involvement (Galovic et al., 2018)—have been associated with both increased seizure susceptibility and heightened neural plasticity. These insights suggest that plasticity itself may contribute to post-stroke epileptogenesis. Accordingly, interventions that modulate neural plasticity such as antiepilepsy drugs (Chang et al., 2022; Ryu et al., 2024) or non-invasive brain stimulation (Adeyemo et al., 2012; Aderinto et al., 2024) hold promise not only for post-stroke rehabilitation but also for seizure prevention in at-risk populations.

In this context, we examined how the combined action of HSP and HSS shapes the interpretation of experimental results. In our NBQX-treated slice cultures, we observed a simultaneous increase in spine size and density after three days of partial inhibition (200 nM NBQX). This enlargement of spines likely enhanced the stability and detectability of small spines in imaging, suggesting that the apparent increase in spine density reflects a combination of spinogenesis, increased spine stability, and potentiation of previously small (and silent) spines. Similar interpretational challenges arise in studies of synaptic scaling results. For example, Turrigiano et al., 1998 reported changes in miniature EPSC amplitudes (but not frequencies) following chronic inhibition and excitation, inferring synaptic weight scaling. However, other studies have documented concurrent changes in both amplitude and frequency. One explanation is the unsilencing of synapses, whereby formerly silent synapses become active (Nakayama et al., 2005), increasing event frequency but lower mean EPSC amplitude due to their relatively small weights. This process may reflect synapse-number-based structural plasticity, which may not be easily detectable with standard imaging methods, depending on variables such as laser intensity, z-resolution, or exposure time. Additionally, under conditions of complete inhibition (50 μM NBQX), we observed distinct and opposing size changes: small spines tended to shrink or disappear, while large spines further enlarged, resulting in an overall increase in sEPSC amplitudes. Such differential response could explain discrepancies between electrophysiological and imaging data. Together, these findings suggest that spine numbers and synaptic weights represent complementary facets of a broader, integrated concept of ‘effective synaptic transmission’. Relying on only one of these measures risks misrepresenting the true functional state of the network.

Our systematic study successfully replicated the empirical observation that chronic activity inhibition leads to spine loss, with the remaining spines often exhibiting enlarged heads (Thiagarajan et al., 2005; Fishbein and Segal, 2007). This phenomenon can be explained by a shared calcium signaling mechanism underlying both HSP and HSS rules. While HSP consistently uses intracellular calcium concentration as a setpoint, most HSS models instead rely on firing rate setpoint (Tetzlaff et al., 2011; Tetzlaff et al., 2012; Tetzlaff et al., 2013; Herpich and Tetzlaff, 2019), synaptic weight normalization (Sullivan and de Sa, 2006), or presynaptic mechanisms (Liu and Buonomano, 2009; Liu, 2011). Notably, one earlier model introduced an integral feedback control algorithm (van Rossum et al., 2000), using a variable with an exponential kernel to track neural spiking activity and adjust synaptic weights. At that time, intracellular calcium had not yet been identified as the principal variable tracking neural activity (Grewe et al., 2010). Building on this idea, we reinterpreted the model’s control variable using a calcium-based variant, allowing us to examine the interaction between the HSP and HSS more mechanistically. In this way, our study also served as a Gedankenexperiment—a conceptual experiment—to explore how structural and functional plasticity interact. If both rules are redundantly monotonic, the one with the faster timescale will dominate. However, when a biphasic HSP rule is combined with a monotonic HSS rule, our simulations showed divergent behavior under strong inhibition: neurons exhibit reduced spine density and enlarged spine heads, mimicking experimental observations. From an energy efficiency standpoint, synapse formation and spine growth are constrained by protein synthesis and ATP availability. It is thus energetically more favorable to potentiate existing synapses than to form new ones. This trade-off supports the complementary roles of HSP and HSS in maintaining network activity. This interplay is also supported by recent theoretical work on dendritic trafficking (Aljaberi et al., 2021) and studies of morphological changes in Alzheimer’s disease (Bachmann et al., 2020). Our findings highlight the interdependency between HSP and HSS, which may underlie divergent results across brain regions, cell types, dendritic compartments, and experimental settings. In particular, the hybrid combination of structural and functional plasticity across dendritic segments supports the emerging view of dendritic compartmentalization (Sutton et al., 2006; Branco et al., 2008; Makara et al., 2009), a principle that is critical to understanding plasticity but remains underrepresented in computational models. More computational work addressing these multiscale interactions is both warranted and encouraged.

While our study addressed both HSS and HSP, it is important to note that the use of NBQX is not commonly used in experiments focused specifically on synaptic scaling. NBQX, a competitive AMPA receptor antagonist, enables precise modulation of network activity via dose titration (Wrathall et al., 1994), offering greater flexibility than tetrodotoxin (TTX), which globally blocks action potentials by inhibiting sodium channels. This property makes NBQX a valuable tool for activity modulation, as demonstrated in our electrophysiological recordings and supported by prior in vivo and in vitro studies (Wrathall et al., 1994; Follett et al., 2000). Notably, a previous study (McKinney et al., 1999) reported a reduction in spine density after applying 20 μM NBQX to hippocampal slices for seven days, with surviving spines predominantly exhibiting a stubby-like morphology. This aligns with our observation of decreased spine density and reduced spine size following 50 μM NBQX treatment and raises the possibility of investigating the functional role of surviving spines under chronic activity blockade (Obashi et al., 2021). However, the specific binding affinity of NBQX for GluA2-containing AMPA receptors merits further attention. These receptors are implicated in both TTX-induced synaptic scaling (Gainey et al., 2009) and glutamate-induced spine protrusion in the presence of TTX (Richards et al., 2005). This selectivity introduces complexity when interpreting NBQX effects, particularly in the absence of TTX. Although there is no definitive evidence that NBQX alone disrupts AMPA receptors synthesis or membrane insertion, its broader impact on synaptic scaling and structural plasticity remains to be fully elucidated. Further studies are warranted to clarify these mechanisms and refine the interpretation of NBQX-based manipulations.

In summary, the use of intracellular calcium concentration as a shared control variable in both HSP and HSS rules closely resembles an integral feedback controller (Figure 9—figure supplement 2), a widely used engineering approach to achieve robust performance despite external perturbations. In both rules, calcium concentration serves as the control signal to drive connectivity kernel modification via negative feedback, which integrates neural firing rates over time and returns a filtered signal responding to its activity. However, the HSP and HSS rules differ in their detailed mechanism. The HSP rule more closely aligns with a Proportional-Integral-Differential (PID) controller, while the HSS rule operates more like a Proportional-Integral (PI) controller. Besides sharing the integral (I) component, the HSP rule’s proportional (P) component is linked to the growth rate of synaptic elements (Diaz Pier et al., 2016), which modifies the number of synaptic elements created or deleted per time unit based on a magnitude relative to the difference between the expected and actual calcium concentration. The scaling factor ( $ρ$ ) in the HSS rule regulates the proportional effects of synaptic scaling. Uniquely, the HSP rule includes a derivative (D) component, determined by the steepness of the growth curve, which is determined by the steepness of the growth curve and dictates how quickly or slowly the system approaches the desired neural activity. This derivative term plays a crucial role in ensuring a smooth and robust response to perturbations, preventing biologically implausible overshooting and abrupt oscillations. The effect is especially pronounced in the Gaussian implementation, where the distance to the target firing rate accelerates differently depending on the position on the calcium concentration axis. Although the exact shape of the growth rule in biological HSP remains to be determined, we propose that a biphasic curve with variable slope offers a redundant and heterogeneous structural backup to synaptic scaling, supporting robust firing rate homeostasis across a range of conditions. This framework may also help anticipate interactions with other plasticity mechanisms, particularly those involving sub- or supra-threshold calcium dynamics (Graupner and Brunel, 2012; Luboeinski and Tetzlaff, 2023).

Materials and methods

Key resources table

Reagent type (species) or resource	Designation	Source or reference	Identifiers	Additional information
Strain, strain background (mice)	C57BL/6 J mouse	Jackson Laboratory	RRID:IMSR_JAX:000664
Strain, strain background (mice)	Thy1-eGFP mouse	Jackson Laboratory	RRID:IMSR_JAX:007788
Peptide, recombinant protein	Streptavidin*	Invitrogen	Cat#: S32354 RRID:AB_2315383	Post hoc labeling (1:1000)
Chemical compound, drug	Paraformaldehyde (PFA)	Carl Roth	Cat#: 0335.3	Final concentration: 4% (w/v) in PB or PBS
Chemical compound, drug	NBQX	Tocris Bioscience	Cat#: 1044	Final concentrations: 200 nM and 50 µM
Chemical compound, drug	DAPI (1 mg/ml in water)	Thermo Fisher Scientific	Cat#: 62248	Post hoc labeling (1:2000)
Software, algorithm	Python	Python Software Foundation	RRID:SCR_008394
Software, algorithm	NEST simulator	NEST Initiative	RRID:SCR_002963
Software, algorithm	Clampfit^†	Molecular Devices	RRID:SCR_011323
Software, algorithm	Image J	NIH ^‡	RRID:SCR_003070
Software, algorithm	Prism (Graphpad)	GraphPad Software	RRID:SCR_002798
Other	OCT ^§	Method described in Del Turco and Deller, 2007	N/A

*

Alexa Fluor 488-Conjugate.
†

pClamp software package.
‡

National Institutes of Health.
§

Organotypic entorhinal-hippocampal tissue cultures.

	plasticity rules	$η$	growth rate	$I_{f a c i l i t a t i n g}$	FOI	$τ_{C a}$	$ρ$
Figure 3Bi, Di	linear	-	100%, 50%, 10%			10 s	-
Figure 3Bii, Dii	Gaussian	0	100%, 50%, 10%			10 s	-
Figure 3Biii, Diii	Gaussian	0.0007	100%, 50%			10 s	-
Figure 4A-C	Gaussian	0.0007	100%	750 pA		10 s
Figure 4E	Gaussian	0.0007	100%	750 pA		10 s	-
Figure 4F	Gaussian	0.0007	100%	0, 200, 500, 750 pA		10 s	-
Figure 5B and C	Gaussian	0.0007	100%	750 pA	110%,95%, 0%	10 s	-
Figure 5D	Gaussian	0.0007	100%	750 pA	[0%, 200%]	10 s	-
Figure 5E	-	-	-	-	-	1 s, 10 s	-
Figure 5F, orange	Gaussian	0.0007	100%	750 pA	0%	1 s	-
Figure 5F, dark blue	Gaussian	0.0007	100%	750 pA	0%	10 s	-
Figure 5F, light blue	Gaussian	0.0007	10%	750 pA	0%	10 s	-
Figure 6B and C	Gaussian	0.0007 for D,0.0004 or0.001 for A	10%	750 pA	0%	10 s	-
Figure 6E and F	Gaussian	0.0007	10%, 100%	750 pA	0% to100% or0% to200%	10 s	-
Figure 7B–F	Gaussian +scaling	0.0007	10%	750 pA	0%	10 s	0, 0.01, 0.02
Figure 8C–F	Gaussian +scaling	0.0007	[0%, 100%]	750 pA	0%	10 s	[0, 0.1]

Share this article

Cite this article

Overview of Hebbian and homeostatic plasticity.

Non-monotonic spine density changes under NBQX inhibition support a biphasic structural plasticity rule.

Establishing a spiking neural network with three distinct structural plasticity rules.

Silent neurons remain isolated under the biphasic Gaussian rule.

Divergent regulation of network connectivity under stimulation and deprivation via three structural plasticity rules.

Activity perturbation and recurrent connectivity jointly shape the evolution of network connectivity.

Homeostatic synaptic scaling (HSS) shapes effective connectivity and interacts with homeostatic structural plasticity (HSP).

Spine size analysis suggests interactions between functional and structural plasticity.

Systematic study of the interaction between homeostatic synaptic scaling and structural plasticity in response to activity silencing.

Parameters for the point neuron model.

Parameters for the network model.

Parameters for the linear structural plasticity model.

Parameters for the zero Gaussian structural plasticity model.

Parameters for the biphasic Gaussian structural plasticity model.

Parameters for the synaptic scaling model.

Protocols for numeric stimulation.

Author details

Han Lu

Present address

Contribution

For correspondence

Competing interests

Sandra Diaz-Pier

Contribution

Competing interests

Maximilian Lenz

Contribution

Competing interests

Andreas Vlachos

Contribution

For correspondence

Competing interests

Citations by DOI

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Categories and tags

Research organism