Figures and data

Learning Mechanisms in direct pathway Striatal Projection Neurons (dSPNs) for the Nonlinear Feature Binding Problem (NFBP)
(A) : Inputs and assumed supralinearity that could solve the NFBP: The NFBP is represented with an example from visual feature binding. In the simplest form of the NFBP, a stimulus has two features, here shape and form, each with two possible values, strawberry and banana, and red and yellow, respectively. In the NFBP the neuron should learn to respond by spiking to two of the feature combinations, representing the relevant stimuli (red strawberry and yellow banana), while remaining silent for the other two feature combinations which represent the irrelevant stimuli (yellow strawberry and red banana). Assuming that each feature is represented with locally clustered synapses, a solution of the NFBP can be achieved when the co-active clusters on a single dendrite, corresponding to a relevant stimulus, evoke a plateau potential, thus superlinearly exciting the soma. Conversely, co-activation of synaptic clusters for the irrelevant combinations should not evoke plateau potentials.
(B) : Synaptic clustering in dendrites: Illustration of how synaptic plasticity in SPNs may contribute to solving the NFBP for a pre-existing arrangement of synaptic clusters on two dendrites. A plasticity rule which strengthens only synaptic clusters representing relevant feature combinations, so that they produce robust supralinear responses, while weakening synapses activated by irrelevant feature combinations, could solve the NFBP.
(C) : Dopamine (Da) feedback: dopaminergic inputs from the midbrain to the striatum (Str) guides the learning process, differentiating between positive feedback for relevant stimuli and negative feedback for irrelevant stimuli. Positive feedback represented by dopamine peaks is necessary for LTP, and negative feedback represented by a dopamine pause is necessary for LTD.
(D) : Signaling pathways underlying synaptic plasticity in dSPNs: Illustrations of signaling components at the corticostriatal synapse that modify synaptic strength (redrawn from Shen et al., 2008). NMDA calcium influx, followed by stimulation of D1 dopamine receptors (D1Rs), triggers LTP (while inhibiting the LTD cascade). L-type calcium influx and activation of metabotropic glutamate receptors (mGluRs) when D1Rs are free of dopamine triggers LTD (while counteracting the LTP cascade).

Characterization of dendritic plateau potentials in the model.
(A) : Somatic voltage, spine voltage, NMDA calcium, and L-type calcium evoked by a cluster varying in size from 1 to 20 synapses. A plateau potential is evoked when glutamate spillover is activated, here triggered when 8 synapses with a weight of 0.25 each are coactivated (corresponding to the “baseline” weights in C).
(B) : Schematic of the neuron morphology with an arrow indicating the stimulated dendritic branch shown in panel A.
(C) : The mean maximal amplitude, with standard deviation shown in bars, of the measures in A averaged over 10 different dendrites, and 10 trials per dendrite. The curves represent clusters with different synaptic weights: baseline (0.25, 0.625 nS), strengthened (0.4, 1 nS) and weakened weights (0.2, 0.5 nS). Synaptic background noise is used in all simulations to elevate the membrane potential to ranges seen in vivo (Reig and Silberberg, 2014).

Example of setup and learning-induced synaptic plasticity
(A) : Illustration of input configuration. The upper panel illustrates the arrangement of features across two dendrites. Each dendrite has synaptic clusters for three features, allowing the representation of each of the four stimuli if seen from the whole neuron’s perspective (but only one relevant stimulus per dendrite). The bottom panel illustrates the stimulation protocol used in the simulation, where stimulus presentation is followed by a dopaminergic feedback signal only if the neuron spikes. Two stimuli are spaced 800 ms apart, to allow for the calcium dynamics to reach baseline levels.
(B) : Example voltage in the soma and the middle of dendrite 1 (d1) and dendrite 2 (d2) before and after learning. Each dendrite here stops responding to the respective irrelevant stimuli during learning.
(C), (D): Evolution of synaptic conductances during learning (top row) and examples of peak Ca²⁺ levels (dots) alongside kernel dynamics in single synapses (three lower rows). Panels in (C) depict clustered synapses in dendrite 1 (d1), where ‘yellow’ and ‘banana’ generally undergo LTP, while ‘red’ undergoes LTD. Panels in (D) show distributed, feature-unspecific synapses. Among these, the purple and blue traces represent synapses that are weakened, while the green trace exemplifies a synapse near one of the clusters that, by chance, achieves a sufficiently high local NMDA calcium level for LTP to dominate. The initial synaptic conductances are drawn from a normal distribution with a mean of 0.625 and a standard deviation of 0.375 nS. MP-K and LTP-K refer to metaplasticity and LTP kernels, respectively, while the solid line represents the midpoint of the kernels, where LTP is strongest. ‘Max’ indicates the peak NMDA calcium during a single stimulus.

Impact of input configurations and synaptic cluster locations on NFBP learning performance
(A): Illustration of creating different input configurations where two, three or four features are placed in two dendritic locations.
(B) : Performance trajectories for all 31 unique input configurations over training epochs. Each line represents the learning performance of a single configuration, with blue traces indicating setups where no dendrite receives more than three features and deep red traces representing setups where at least one dendrite contains all four features. The inset shows the number of configurations in each group that reached the NFBP learning criterion: 18/18 for the blue group and 5/13 for the deep red group.
(C) : Performance on the last 160 stimuli for the two groups of configurations: (C1) setups where each dendrite has at most three features and (C2) setups where at least one dendrite contains all four features.
(D) and (E): Performance in a three-feature configuration as a function of cluster location. (D) illustrates the distribution of synaptic clusters across dendritic locations, and (E) shows performance over the last 160 stimulus presentations as a function of somatic distance for these configurations. A total of 60 unique configurations were tested, with synapse clusters randomly assigned to different dendritic locations. The solid line represents a quadratic fit to the data.

Effects of inhibitory inputs on performance
(A) Dendritic input configuration with inhibitory synapses added. The setup of the excitatory and feature-unspecific synapses is the same as in Figure 3A (the feature-unspecific synapses are not shown). Inhibitory synaptic connections for each of the four features are added in both dendrites, with one synapse per feature.
(B) : Displays average performance for 31 unique input configurations of 2, 3, or 4 features on two dendrites as a comparison between the setup with (orange) and without inhibitory plasticity (blue). The bar plot inset shows the number of configurations that successfully solved the NFBP with and without inhibition.
(C) : Performance over the last 160 stimulus presentations as a function of somatic distance of the synaptic clusters for the input configuration in (A). A total of 60 unique configurations were tested, with synapse clusters randomly assigned to different dendritic locations. The solid line represents a quadratic fit to the data. The dashed line is the corresponding quadratic fit from Figure 4E.
(D) : Peak voltage-gated calcium (left panel, dots) and plasticity thresholds dynamics (lines) during learning (number of training examples). The upper threshold (CTH) moves towards peak calcium levels while the lower threshold (CTL) moves towards the next highest level. (Right) Inhibitory synaptic conductances for the synapses in dendrite 1. Strengthened inhibitory synapses prevent excitation by the corresponding features (here ‘red’ and ‘strawberry’), while weakened inhibitory synapses allow excitation by their corresponding features (here ‘yellow’ and ‘banana’).
(E) : Excitatory synaptic conductances (top) and calcium levels and kernel dynamics (below) over learning. The specific conductances in the top panel corresponding to the kernel dynamics and calcium levels in the bottom panels are indicated with arrows. Note that kernel and calcium dynamics for the example ‘yellow’ synapse refer to the only ‘yellow’ synapse in the top panel which initially weakens and is later strengthened. MP-K, LTP-K stands for metaplasticity- and LTP kernel, while the solid line shows the midpoint of the LTP kernels. Max refers to the peak NMDA calcium during a single stimulus.

Performance analysis of learning using distributed synaptic inputs
(A) : Example illustration of synaptic distribution before (top) and after learning (bottom) of the 200 excitatory and 60 inhibitory inputs. S - strawberry, Y - yellow, B - banana, R - red, U - feature-unspecific inputs. Bordered markers (circles and diamonds) indicate synapses whose conductances are shown in (C).
(B) : Performance over training epochs with and without inhibitory plasticity. (Left) Mean performance of the setups with and without inhibitory synapses. (Right) Individual traces for the setup with inhibition. Each of these individual trials (31 in total) uses a unique random distribution of synapses.
(C) : Example of summed synaptic conductances (left) and voltage (right) in the soma, and four example dendrites (d1-d4) of the synaptic distribution in (A) in a trial where the NFBP is successfully solved. The sums of both excitatory (left) and inhibitory (right) inputs are shown.

Synaptic Plasticity Rules: Calcium and Dopamine Interactions in Synaptic Weight Modification
(A) : Synaptic weight updates following a dopamine peak. (Left) LTP kernel is a bell-shaped curve describing the amount of weight increase, which happens over a region of [Ca]NMDA. (Right) A wider bell-shaped kernel, i.e. the metaplasticity kernel, determines how the LTP kernel is shifted along the calcium level ([Ca]NMDA) axis following a peak in dopamine.
(B) : Synaptic weight updates following a dopamine pause. (Left) The LTD plasticity kernel. The LTD threshold is constant and set at 70 nM (Right). Metaplasticity describing how the LTP kernel shifts along the calcium axis following a dopamine pause.
(C) : A schematic of how the LTP kernel is shifted following a dopamine peak, and examples of NMDA calcium levels before and after the shift. The NMDA calcium levels increase following activation of the strengthened synapse due to the previous dopamine peak (illustrated with the red circle jumping to the blue circle). The LTP kernel also moves downwards, and in this example the strengthened synapse (blue circle) would stabilize, and thus would not increase in strength even though additional rewards arrive.
(D) : A schematic showing how the LTP kernel is shifted following a dopamine pause, and examples of NMDA calcium levels before and after the shift. The calcium level shortly before the dopamine pause is shown with the red circle. The dopamine pause causes a decrease in synaptic weight and a smaller calcium response results in the next synapse activation (illustrated with the red circle jumping to the blue one). Also, the LTP kernel moves up towards higher NMDA calcium levels, and thus the weakened synapse (blue circle) is more likely to stay out of the plastic region of the LTP kernel (unless it in the future is activated with dopamine peaks more regularly than with dopamine pauses).
(E) : Illustration of the inhibitory plasticity rule. (Left) Changes in synaptic weight for active (beige) and inactive (blue) synapses based on voltage dependent calcium levels in the dendritic shaft at the location of the inhibitory synapse. The dashed lines show the minimum/lower (CTL) and maximum/higher (CTH) threshold. (Right) Functions for updating the lower and higher thresholds with voltage dependent calcium level. The asterisk denotes the Ca level where the curves for active and inactive synapses meet, which is a point of zero update.


Excitatory Plasticity Parameters


Inhibitory Plasticity Parameters

Learning-induced synaptic plasticity with metaplasticity turned off
(A): Example voltage traces recorded at the soma and in two dendrites before and after learning. Without metaplasticity, no significant learning occurs—every stimulus triggers spikes (50% performance). (B, C): Evolution of synaptic conductances throughout learning (top row) and characteristic examples of peak calcium levels (dots), along with the constant kernels in single synapses (shaded areas, bottom rows). The same synaptic arrangement and initial synaptic weights as in Figure 3A is used (for both clustered and feature-unspecific synapses). (B) illustrates changes in clustered synapses on dendrite 1 (d1). (C) shows distributed, feature-unspecific synapses. The final outcome depends on the initial synaptic weights. Synapses whose calcium levels are within or above the LTP region cannot be weakened, because any LTD during irrelevant stimuli is counteracted by LTP occurring for the relevant stimuli. The calcium levels remain “trapped” within the LTP kernel, and the synaptic weights “zig-zag” around a steady state level (seen in the example clustered ‘red’, ‘yellow’ and ‘banana’ synapse, as well as for the example green and blue feature-unspecific synapse). The synapses with small initial weights, whose calcium levels are situated below the LTP kernel are weakened (example purple synapse). The MP-K and LTP-K denote metaplasticity and LTP kernels, respectively, with the solid line representing the midpoint of the kernels where LTP is strongest. In summary, without metaplasticity learning of NFBP is difficult if initial synaptic weights are not optimally tuned with regard to the postsynaptic calcium concentration. Metaplasticity allows the synapses to fall in and out of the LTP kernel range. They can be strengthened as long as their calcium levels are within the LTP kernel, are stabilized when they are above the LTP kernel, and can be weakened when they are below the LTP kernel.

Learning with and without inhibitory inputs in distal dendrites.
(A) The figure illustrates the impact of inhibitory inputs on learning dynamics in two distal dendrites (d1 and d2). The top panels depict somatic voltage, while the middle and bottom panels show dendritic voltages. The left panels show voltage in the middle of the learning process (approximately 480 training epochs) with and without inhibitory inputs. The right panels show voltage at the end of the learning process (approximately 960 training epochs). The dashed black line corresponds to the scenario without inhibitory inputs, and the solid red line indicates the presence of inhibitory inputs. In the middle of learning, inhibition prevents spiking for the irrelevant stimuli, thus leading to less LTD. This preserves synaptic strengths in both the clustered and feature-unspecific synapses, allowing the neuron to spike for the relevant stimuli at the end of learning. The overall decrease in unspecific input was 39% without inhibitory inputs and 35% with inhibitory inputs.
(B) The effect of inhibitory plasticity on learning in a dendrite receiving all four features. The left panel shows excitatory synaptic weights in d1 without inhibition, where all features weaken over time due to competition. The middle panel illustrates excitatory synapses in the presence of inhibition, where learning is better regulated. The right panel shows inhibitory synaptic weights.

Other learning paradigms
(A) : Learning of 5×5 and 3×3 feature combination tasks with randomly distributed synapses over 30 dendrites. The neuron is trained to spike for black-marked combinations and remain silent for white-marked ones. For the 25-feature (5×5) combination, the model used 500 excitatory synapses (50 per feature) with initial conductances of 1.125 ± 0.25 nS, along with 100 inhibitory synapses (10 per feature) with initial conductances of 0.1 ± 0.01 nS. For the 9-feature (3×3) combination, the model used 312 excitatory synapses (52 per feature) and 96 inhibitory synapses (16 per feature), with the same initial weights as in the 5×5 case.
(B) : Learning of a 2×2 linear feature combination task with randomly distributed synapses. Each feature is represented by 75 synapses distributed across 30 dendrites (300 total synapses), with an initial synaptic weight of 0.375 ± 0.25 nS. The top row shows the feature combinations used, followed by the performance across 30 different synaptic distributions. The bottom row illustrates somatic voltage traces before and after learning, where the neuron starts from a silent state (red trace) and learns to spike through metaplasticity and dopamine-driven weight updates. Additionally, in the subthreshold learning paradigm, dopamine was given without spiking for expected firings and suppressed when unexpected spiking occurred.
(C) : Evolution of synaptic conductances (left column) and calcium dynamics with metaplasticity activity (right column) during learning for each feature (X1, X2, Y1, Y2). The left panels display the weight adjustments of all synapses corresponding to each feature, with one synapse per feature bolded to highlight its individual trajectory. The right panels show the corresponding calcium traces and kernel dynamics for the bolded synapse. MP-K and LTP-K refer to metaplasticity and LTP kernels, respectively, with the solid line representing the midpoint of the kernels where LTP is strongest. ‘Max’ indicates the peak NMDA calcium level during a single stimulus.