Abstract

Medium spiny neurons (MSNs) in the striatum, receiving input from cortex and projecting to the direct and indirect pathways express D1 and D2 type dopamine receptors, respectively. Recently, it has become clear that the two types of MSNs markedly differ in their mutual and recurrent connectivities as well as feedfonlvard inhibition from FSls. Therefore, to understand striatal function in action selection, it is of key importance to identify the role of the distinct connectivities within and between the two types of MSNs on the balance of their activity. Here, we used both a reduced firing rate model and numerical simulations of a spiking network model of the striatum to analyze the dynamic balance of spiking activities in D1 and D2 MSNs. We show that the asymmetric connectivity of the two types of MSNs renders the striatum into a threshold device, indicating the state of cortical input rates and correlations by the relative activity rates of D1 and D2 MSNs. Next, we describe how this striatal threshold can be effectively modulated by the activity of fast spiking interneurons, by the dopamine level, and by the activity of the GPe via pallidostriatal backprojections. We show that multiple mechanisms exist in the basal ganglia for biasing striatal output in favour of either the ‘Go' or the ‘No-Go' pathway. This new understanding of striatal network dynamics provides novel insights into the putative role of the striatum in various behavioral deficits in patients with Parkinson's disease, including increased reaction times, L-Dopa-induced dyskinesia, and deep brain stimulation-induced impulsivity.

Author Summary

BG dysfunction leads to brain disorders such as Parkinson’s disease. At the main input stage of the BG, the striatum, two competing pathways originate. Neurons projecting on these pathways either express D1 or D2 type dopamine receptors. Because activity of D1 or D2 neurons facilitate go or no-go type decision, it is important to study the balance of the D1 and D2 neuron activity. Contrary to the common assumption thus far, recent data shows an asymmetry in striatal circuit with D1 receiving higher inhibition from D2 and fast-spiking neurons. Here, we studied the functional implications of the asymmetric connectivity between D1 and D2 neurons. Our analysis and simulations show that the asymmetric connectivity between these neurons gives rise to a decision transition threshold (DTT), as a consequence D1 (D2) neurons have higher firing rate at lower (higher) average cortical firing rates. Importantly, DTT can be modulated by input correlations, local connectivity, feedforward inhibition and dopamine. Our results suggest that abnormal changes in the DTT could be a plausible mechanism underlying the cognitive and motor deficits associated with brain diseases involving BG malfunction.

Introduction

The striatum is the main input stage of the basal ganglia, receiving inputs from widely distributed areas in the corteX [1] , and projecting to the BG output nuclei Globus Pallidus Interna (GPi) and Substantia Nigra (SNr) via the so-called direct and indirect pathways, respectively [2]. The integration of multi-modal sensory signals with motor and/ or cognitive inputs in the striatum sets the stage for action selection. Therefore, to understand the computations performed by the basal ganglia, it is of key importance to characterize the dynamical properties of striatal network activity.

The remaining 5% are inhibitory interneurons, such as fast spiking interneurons (FSIs) and tonically active neurons (TANs) [3]. The MSNs are classified as D1 and D2 type neurons, depending on their dopamine receptor eXpression. Interestingly, D1 MSNs project to the GPi and SNr, forming the direct pathway, whereas D2 neurons project to Globus Pallidus Externa (GPe), forming the indirect pathway [2]. Consistent with basal ganglia anatomy and previous models [4—6], selective activation of D1 MSNs in the rat increases ambulation, whereas selective activation of D2 MSNs increases freezing behavior [7]. However, a complete shutdown of activity in either of the two neuron subpopulations might not occur in awake behaving animals [8]. In such a scenario,though, action selection could still be performed by a relative increase in the activity of one subpopulation compared to the other.

In such single population models, the striatal output is controlled by the strength of cortico-striatal synaptic weights. The recurrent inhibition is not strong enough to support win-ner-take-all dynamics as earlier speculated [10], but may be sufficient to allow for a winner-less competition [11] and may enhance the saliency of the cortical input representation in the stria-tum [12]. This, however, is a highly simplistic view of the recurrent inhibition within the stria-tum and, as will be described below, is inconsistent with experimental data, especially given the recent findings on the recurrent connectivity in the striatum.

Moreover, the striatal circuit also shows a highly specific connectivity in terms of the mutual inhibition between the MSN subpopulations [14, 15] and the feedforward inhibition from FSIs [16]. Paired neuron recordings showed that D2 MSNs make more and stronger connections to D1 MSNs, than vice versa. Furthermore, FSIs preferentially innervate D1 MSNs as compared to D2 MSNs (Fig 1A). The computational role of this specific connectivity within the striatum is not clear and cannot be inferred from previous models, which assumed a single homogeneously connected MSN population in the striatum. Specifically, it is important to identify the effect of the distinct D1 and D2 connectivities on the competition between the direct and indirect pathways.

We show that the firing rates of both D1 and D2 MSNs change in a non-monotonic manner in response to cortical input rates and correlations. Interestingly, higher output rates in D1 than in D2 MSNs, and vice versa, were observed for separate, non-overlapping ranges of cortical input rate and correlation. Correlations in the input can further change the range of cortical inputs for which either D1 or D2 MSNs have the higher firing rate. Thus, we argue that the striatum acts as a threshold device for cortical input rates by changing the magnitude of the difference between firing rates of D1 and D2 MSNs, depending on the level of cortical input rates and correlations. While the main determinant of the striatal threshold is the asymmetric connectivity among the various striatal elements, the threshold is not fixed and can be dynamically adjusted by the dopamine level, by the connectivity and firing rate of fast spiking interneurons (PSI), and by the GPe activity. That is, changes in the striatal threshold could reflect changes in the operating point of the striatum, behavioral context, and learning and reward history. These novel insights concerning the interactions between direct and indirect pathways suggest putative mechanistic explanations for the role of striatum in cognitive deficits such as L-Dopa-induced Dyskinesia (LID), deep brain stimulation (DBS)induced impulsivity and increased reaction times in Parkinson’s disease (PD) patients.

Results

These two pathways converge in the GPi/SNr, their relative activity balance lets the animal choose between a ‘Go’ or a ‘No-Go’ action [4, 5]. In this framework, the decision making process in the striatum is mediated by the selective activation of one of the two MSN (D1 or D2) subpopulations. However, recent eXperimental data suggest that such complete shutdown of activity in either one of the two neuron subpopulations may not occur in awake behaving animals [8]. Hence, we should consider the alternative possibility that action selection is performed by a relative increase in the activity of one subpopulation compared to the other. Therefore, to understand the effect of the recurrent, mutual connectivity and feedforward inhibition [14—16] on the relative balance of the activities in the direct and indirect pathways we studied the dynamics of the striatal network. Specifically, we evaluated the firing rates of the D1 and D2 MSNs, in response to cortical input rates and input correlations.

D1 MSNs require overall stronger input from cortex than D2 MSNs

In th6-GFP transgenic mice, FSIs preferentially target D1 MSNs as compared to D2 MSNs [16] (however cf. [15] ). This leads to the following two inequalities: 112 > [21, hp 2 IZF, with 112 denoting the connection from D2 MSN to D1 MSN, 121 the obverse, [IF the connection from FSI to D1 MSN and 12F likewise to D2 MSN (cf. Table 1). These inequalities imply that if the two MSN subpopulations receive the same amount of excitatory input, D2 MSNs will always have a higher firing rate.

In our spiking network simulations this corresponded to a lower mean firing rate of the D1 population compared to the D2 population. This bias might be considered functionally useful, because for similar input strengths (i.e. when the cortex does not impose any preference for either direct or indirect pathways), the indirect pathway will dominate the striatal network dynamics and the ‘No-Go’ would be the default state of the basal ganglia.

Moreover, pyramidal tract neurons have stronger connections to D1 MSNs than to D2 MSNs, while intratelencephalic tract neurons innervate D1 and D2 MSNs equally. [1]. Despite this evidence, current data are not sufficient to determine the exact strength and rate of excitatory inputs to D1 and D2 MSNs. To estimate how much additional excitation would be required for D1 MSNs to have their firing rates exceed over those of D2 MSNs, we systematically varied the drive of cortical inputs to D1 and D2 MSNs and calculated the response firing rates of the two subpopulations, for the firing rate model (Fig 2). In the following we discuss two different scenarios for systematically varying the cortical excitatory drive to D1 and D2 MSNs. Scenario 1: Increase the relative cortical input rate to D1 MSNs. One possibility is that the cortical input rate to D1 MSNs is higher than to D2 MSNs while the strength of the cor-tico-striatal synapses to the two types of MSNs is the same. This results in an ‘additive’ scenario, in which we kept the strength of the cortical inputs to D1 and D2 MSNs equal (IO;L = ICZ), and systematically increased the firing rate of the cortical inputs to both D1 and D2 MSNs:

ictx_d1 ictx_d2 ACTX

AMSN

This scenario is equivalent to a state when no learning has happened and the cortex has no preference for choosing between activating D1 or D2 MSNs. Here, we considered only those parameter ranges for which neither of the two subpopulations was completely shut down by mutual and/ or recurrent inhibition, because that would be a trivial solution of the network dynamics, which is also not supported by the experimental data [8]. When the two types of MSNs are driven with identical input rates (7Wle 2 law”), the output firing rate of D2 MSNs always exceeded that of D1 MSNs, irrespective of the magnitude of the input rate (Fig 2A: main diagonal, Fig 2B: points corresponding to ACTX = 0.0Hz). This shows the inherent bias of the striatal network towards higher firing rates of D2 MSNs [17] , for the reasons explained above. When the input rate to D1 MSNs exceeds that to D2 MSNs by a sufficiently large amount (here ACTX 2 2H2), D1 MSNs always fired at a higher rate than D2 MSNs, irrespective of the magnitude of the cortical input rate (Fig 2A: area below the dashed line, Fig 2B: area corresponding to ACTX > 2.0HZ).

In this case, the ‘winner’ of the competition between D1 and D2 MSNs depends on the size of the cortical input rates. For small input rates to the two MSNs, D1 MSNs have the relatively higher firing rates (i.e. AMSN > 0 HZ). However, as both input rates were increased, while maintaining the small difference between the two, the recurrent inhibition from D2 to D1 MSNs increased such that, beyond a certain input rate (here % 10 Hz), it exceeded the extra small excitatory drive to the D1 MSNs, resulting in the output rate of D2 MSNs surpassing that of D1 MSNs (i.e. AMSN < 0 HZ). Hence, in this transitional region with small differences in the excitatory inputs to both MSN types (OHZ < ACTX < 2H2), it was possible to have higher firing rate for either of D1 or D2 MSNs, depending on the size of the cortical input drive (Fig 2A, dashed line and 2B, set of points marked by the oval).

As a result, the difference between the two output rates AMSN reverses as the cortical input rate is increased. In other words: due to the asymmetric connectivity between D1 and D2 MSNs, the striatum can flexibly choose between the direct and the indirect pathway, depending on the level of the cortical input firing rates. By contrast, in a regime with a large difference between the excitatory drives to D1 and D2 MSNs, the striatum would lose this ability to switch between the two, depending on the cortical input. To understand the mechanism underlying the sign reversal we derived an eXpression for AMSN. In this input scenario the dynamics of the two MSN populations are given by:

It should be noted that absolute values of parameters are used for the analysis, and to calculate km and km (Eqs 4, 5) correct sign of the various weights must be used, i.e. I11, I12, I22, I21, I 1 F, IZF are negative and similarly, ACTX is positive. Where ACTX is again the extra excitatory input given to D1 MSNs. The response rate difference AMSN 2 km — km now reduces to:

The first term in Inpstr(term I) has an effective negative contribution to AMSN, due to the higher effective recurrent inhibition to D1 ((I11 + [12) > (I22 + I20). The contribution of FSI inhibition (term 11) could be effectively positive or negative, depending on the ratio of feedforward (I11: and IZF) and recurrent ((IuL + I12) and (I22 + I20) inhibition. In our case, it is negative because (IZFUH +

This approximative result of the linearized equation without leak was confirmed when we numerically calculated AMSN from the full Eqs. 2—3 (Fig 3C).

Therefore, to understand the contribution of the activity of fast spiking neurons (assuming their overall contribution to be negative), we estimated AMSN while independently varying XCTX and M31 (Fig 3A). We found that for a fixed value of XCTX the sign of AMSN could be reversed by an increase in M31. Similarly, for a constant value of XFSI the sign of AMSN could be reversed by an increase in XCTX. In fact, in this scenario, AMSN reverses its sign with respect to XCTX even in the absence of feedforward inhibition. Hence we conclude that the recurrent connectivity between the two MSN populations is the main determinant of the AMSN sign reversal. Note that in Fig 3C and 3B we showed the sign reversal of AMSN When both M31 and XCTX increased in proportion (cf. black lines in Fig 3A).

Instead of increasing the rate of cortical excitation, we could also increase the strength of the cortical synapses to the D1 MSNs to provide these neurons With more input. This gives rise to a second, in this case a ‘multiplicative’ scenario (Fig 2C and 2D):

As in scenario I, we estimated AMSN by systematically varying the input rate (XCTX) and the strength of the cortico-striatal projections (I01, ICZ). Notice that in this ‘multiplicative’ scenario the difference between the input drives of D1 and D2 MSNs increases proportionally with XCTX. Because the cortical input to D1 MSNs is now enhanced by the larger synaptic strength (IC1 > I02), the diagonal in Fig 2C is no longer in the regime where D2 MSNs dominate (AMSN < 0). Similarly to scenario I, also here also we found only a very small regime of input drive for which either D1 or D2 MSNs can have higher firing rate as cortical input rate is changed (the oval area marked in Fig 2D). Therefore, we argue that for a flexible switching between direct and indirect pathways the difference between the input rates/weights to the two MSN populations and their outputs rates should be small. As in the ‘additive’ scenario, we performed a simplified analysis of the dynamics of Eqs 21 and 22 without leak term and a linear transfer function and found the fixed point rates 7L D1 and AMSN now reduces to This expression can be split into the effective contributions of D1 (Dlefic) and D2 (Dzefi), respectively: Cortico-striatal projections are stronger than interneuronal projections to MSNs, therefore, IC1 > I11: and [CZ > IZF, and hence, Dlefic is positive and DZeff is negative. At low cortical rates, | Dlefl| is larger than | DZeficl due to the scaled cortical excitation for D1 (ICl > ICZ) and hence, the D1 type neurons dominate. However, for higher cortical driving rates, stronger feedforward inhibition to the D1 MSNs ensures that D2 MSNs exceed the firing rates of D1 MSNs, reversing the sign of AMSN (Fig 3D). The values are qualitatively similar to AMSN calculated numerically from the full Eqs 21 and 22 with the leak term and nonlinear transfer function (Fig 3F). In general, as is evident from the Eq 12, the sign of AMSN is determined by both XCTX and M31. However, in this input scenario the relation between AMSN and the firing rates XCTX and M31 is more complex than in the ‘additive’ input scenario. Similarly to the ‘additive’ input scenario, for a fixed value of XCTX the sign of AMSN could be reversed by an increase in M31 (Fig 3D). However, for a constant value of 21:31, an increase in XCTX changes the sign of AMSN from negative to positive. This is because a large value of XCTX is needed to overcome the effect of preferential connections of the FSIs to D1 MSNs. Indeed, the higher the firing rate of FSIs, the stronger cortical input would be needed to switch the sign of AMSN. Note that in Fig 3F and 3E we showed the sign reversal of AMSN when both 7L F51 and XCTX increased in proportion (cf. black lines in Fig 3D).

The ‘multiplicative’ input scenario requires preferential inhibition from FSis, whereas the ‘additive’ input scenario relies on asymmetrical projections between D1 and D2 MSNs (however, see also later in subsection “Ef-fect of symmetrical FSI projections on the DTT”). In all subsequent analyses, we used the scenario II, keeping ICl > [CZ and maintaining identical input rate XCTX to both D1 and D2 MSNs. The reason is that scenario I is not suitable to compare the effect of input correlations, as this requires the input rates for D1 and D2 MSNs to be equal. Decision transition threshold in the striatum

the ‘multiplicative scenario’ for detail). This scenario is shown for the firing rate model in Fig 1B and for the spiking neural network in Fig 4). In this state, we estimated the difference between D1 and D2 MSN output rates as a function of cortical input rate. For a constant difference in synaptic strength to D1 and D2 MSNs (AI = IC1 — [CZ = const), the difference between the output firing rates of D1 and D2 MSNs (AMSN) changes from positive to negative as the cortical input rate increases (Figs 1B and 4). That is, the response of the striatum to a range of cortical input rates can be divided into two regimes: a regime where D1 firing rate exceeds D2 firing rate and a regime where this situation is reversed. At the transition point 22%;, the D1 and D2 firing rates are equal. This behavior can be interpreted as a striatal bias to-tran wards the ‘Go’ pathway below a certain cortical input rate (XCTX < ACTX), switching to a bias

That is, the striatum acts as a threshold device, indicating the change in the cortical input rates by the sign of the difference between the firing rates of D1 and D2 MSNs—higher D1 activity reflecting lower cortical input rate and higher D2 activity reflecting higher cortical input rate. Therefore, we propose that the asymmetric striatal connectivity allows it to ‘sense’ small changes in the cortical activity reaching the striatum and to modulate the balance between ‘Go’ and ‘No-Go’ states of the basal ganglia accordingly.

It is possible to tune the strength of FSI inputs such that initially the striatum is biased towards the ‘NoGo’ pathway and a subsequent increase in the cortical input changes the bias to the ‘Go’ pathway (cf. Fig 3A and 3D).

However, in their model the striatum was considered a passive unit, the threshold of which could be tuned by the strength of cortico-striatal synapses. In our model, for the first time, we describe the striatum as an active participant in deciding the level of DTT. It should however be noted, that our model eXplores the DTT under steady state conditions only. In this regime, the striatal network settles into stable fixed points (cf. Materials and Methods) and, therefore, the DTT is a stable network state as well. A DTT might also emerge as a transient state in the network dynamics in the presence of additional mechanisms like short term and/ or long term plasticity of cortico-striatal and/ or striato-striatal synapses and intrinsic plasticity of striatal neurons. Further exploration of such mechanisms is, however, beyond the scope of the present work. The presence of such dynamic DTT suggests several possible control mechanisms, that may set the bias towards either of the two pathways along the basal ganglia downstream nuclei, depending on the animal behavioral state, motivation and acquired learning. In the following sections, we describe three possible mechanisms that can modify the balance of D1 and D2 activity by modulating the level of the DTT, thereby rendering it as a dynamic variable.

However, recent data shows that the basal ganglia network is more complicated than the simple feedforward model (cf. review by [19] ). In future it would be important to show whether GPi can change its activity according to the D1 and D2 MSNs firing rate differences, in a more updated model of the basal ganglia.

Robustness of DTT

Therefore, we tested the viability of the concept of DTT when the mean firing rate of the cortical inputs to the D1 and D2 MSNs was chosen from a low-pass filtered noise. In addition, the extra input to the D1 MSNs was chosen from a low-pass filtered noise. To quantify the DTT we measured AMSN in pre-DTT and post-DTT regions (Supplementary SI Fig). It is obvious that when ACTX is noisy the difference between the firing rate of the D1 and D2 MSNs around the DTT decreases, however AMSN in the pre-DTT and post-DTT remains large enough to be reliably measured (Supplementary SI Fig). When the standard deviation is increased upto 50% of the mean value of the ACTX, both pre-DTT and post-DTT areas decreased only by 20% (Supplementary Sl Fig). Correlations in the spiking is another way of generating transient fluctuations in the input to the D1 and D2 MSNs and the effect of correlation induced fluctuations is described below (cf. subsection “Effect of cortical spiking activity correlations on the DTT”)

Therefore, we tested the robustness of the DTT for variation in the chosen values of striatal network connectivity (e.g. I11, I12 etc.). There are up to eight network connectivity parameters (cortical input to the two MSNs populations, mutual and within MSNs connectivity and MSN and FSI connectivity) which could potentially affect the existence of the DTT. We systematically varied all the parameters (except the PSI connectivity to the MSNs) by i 100% around their mean values (as used throughout this manuscript). Visualisation of this high dimensional parameter space is difficult. Therefore, here we checked if the DTT exists when only one of the parameters is varied systematically by :100% of its mean value, while other parameters can take any value between i100% of their mean value. A DTT is said to exist if AMSN for increasing values of km and the given parameter combination changes sign from positive to negative values exactly once. This was detected algorithmically as well as by visual inspection. Indeed, the DTT is a robust to such changes in the network connectivity and for every parameter we tested here (cortical inputs to the MSNs and within and mutual connectivity of MSNs), we can find the DTT (cf. Supplementary S2A—S2F Fig). Only for very small values of the within connectivity of the D1 and D2 MSNs (I11 and I22),

In addition, we also checked the existence of the DTT for various pairs of the network connectivity parameters (e.g. 121 vs. I11). While there are some forbidden regions in which DTT does not exist, for a large variations in the striatal connectivity, the DTT can be observed (Supplementary 82 Fig). When the existence is visualised as a scatter diagram as a function of a pair of the network connectivity parameters, the mean of these clusters lie very close to the values we have used throughout this manuscript.

Effect of cortical spiking activity correlations on the DTT

Here we investigated how spike correlations in the cortical inputs to the striatum may affect the balance between D1 and D2 MSNs. Injection of correlated inputs to a population of MSNs requires the separation of input correlations into two categories [12]: correlations arising due to convergence of cortical projections onto individual neurons, i.e correlations within the pool of pre-synaptic neurons seen by individual MSNs (in the following referred to as ‘within’-correlations—W) and correlations arising due to divergence of cortical projections, i.e. correlations between the pools as seen by different MSNs, the so-called shared input (in the following referred to as ‘between’ correlations—B) (cf. Materials and Methods and Fig 5A). Note that because both B and W correlations arise from convergent and divergent projections of the same pre-synaptic population, between correlations are always smaller than or equal to within correlations(B 3 W [21, 22],cf. Materials and Methods).

These simulations were only performed for the spiking neural network model since modelling correlations in a mean field model is nontrivial, especially when post-synaptic neurons are recurrently connected. (cf. Materials and Methods, Fig 5A). As expected, within-correlations affected the individual firing rates in a non-monotonic manner [23] for both D1 and D2 MSNs. For weak input correlations, as D1 MSNs form stronger synapses with cortical afferents, (ICl > ICZ; scenario II), D1 MSNs spike at higher rate than D2 MSNs (Fig 5B). For high values of W the MSNs operated in a so called ‘spike-wasting’ regime, because the average count of coincident spikes (or the input variance) exceeded the spike threshold [23] and therefore, W caused no difference in the output firing rates of D1 and D2 MSNs (Fig 5C). That is, the difference in the firing rates of D1 and D2 MSNs (AMSN) changed non-mono-tonically as a function of W, reaching its maximal values for W = Wopt 7S 0 (Fig 6A and 6B).

For W < Wopt, increasing B resulted in a monotonic increase in A MSN. In this regime, B increased the correlation within D1 and D2 MSNs and, thereby, increased the effective inhibition of one population by the other, and vice versa. Because D1 MSNs received input with a slightly higher strength (IO;L > I02; scenario II), B induced more synchrony among D1MSNs and, hence, they were better able to inhibit the D2 MSNs (Fig 5E). However, this held only over a smaller parameter regime, because by design B cannot exceed W (cf. Materials and Methods). For W 2 Wopt, A MSN decreased monotonically with increasing B. In this regime, an increase in B enhanced the correlations within D1 and D2 MSNs, hence amounting to a ‘wasting of recurrent inhibition’ [12]. Increased correlations within D1 and D2 MSNs created strong but only short lasting inhibition leaving out longer lasting time windows without any recurrent inhibition (Fig 5D). In the absence of sufficiently long-lasting recurrent inhibition the difference A MSN monotonically decreased (Fig 6A and 6B).

Our results here now show that the ‘between’-correlations could play an even more important role in initiating the action-selection process in the striatum. However, whether B correlations improve or impair action selection strongly depends on the ‘within’-correlations. (Fig 6). With an increase in the cortical input rate, the amount of inhibition experienced by a MSN in the striatal network also increases due to the increase in MSN firing rates. This requires more coincident excitatory inputs (i.e. higher W) to overcome the increased inhibition and to drive the output neurons above firing threshold. Thus, Wop, increases with an increase in input rates (compare Fig 6A and 6B). That is, an increase in cortical input rates broadens the dynamic range of the interaction between B and W by increasing Wop, (Fig 6C, also discussed in [23]).

To illustrate this, we plotted AMSN as a function of the cortical input rates and within correlations (W), for a fixed value of B (Fig 6C, cf. Supplementary S4 Fig). Here, we can define three regions in the space of input firing rates and input correlations for which AMSN is positive, negative and zero respectively (Fig 6C), reflecting how the striatum will react to the first and second order statistics of the cortical inputs by activating the direct, the indirect or none of the two pathways. While the direct and indirect pathways reflect the striatal preference for a ‘Go’ and ‘No-Go’ action, respectively, AMSN % 0 could reflect a state of high conflict where it is useful to halt the decision making process until more information is available.

Effect of Dopamine on the DTT

Dopamine affects striatal function by modulating the intrinsic excitability of the MSNs and synaptic weights [24] and synaptic plasticity [25, 26] of the cortico-striatal projections. Because we are interested in how striatum processes cortical inputs, here we studied the effect of dopamine-induced modulation of the strength of cortico-striatal projections. Dopamine induces long term potentiation (LTP) in cortico-striatal synapses onto D1 MSNs and long term depression (LTD) in cortico-striatal synapses onto D2 MSNs [24]. This steady state effect of dopamine is often modelled by changing the weight of cortico-striatal synapses for D1 ([0) and for D2 (ICZ) [27]. Therefore, in our model we simulated the effect of dopamine depletion by decreasing ICI and increasing ICZ, whereas higher than normal dopamine was simulated by increasing I01 and decreasing I02.

This shifted the DTT towards lower cortical input rates (Fig 7A) reducing the regime under which the direct pathway could dominate the indirect pathway. Consistent with experimental observations and previous models [4, 5], our model suggests that dopamine depletion would (1) increase the firing rate of D2 MSNs, thereby, (2) introducing a preference for ‘No-Go’ type actions. By contrast, in a high dopamine state (e.g. as in PD patients with L-Dopa treatment) D1 MSNs received much stronger cortical input and were able to maintain higher output firing rate than D2 MSNs for a wider range of cortical inputs. That is, high dopamine shifted the DTT towards higher cortical inputs (Fig 7A) and increased the regime under which the direct pathway could dominate the indirect pathway. Consistent with experimental observations, our model suggests that an increase in steady state dopamine levels would (1) increase the firing rate of D1 MSNs [28] , thereby (2) introducing a preference for ‘Go’ type actions or dyskinesia. In summary, therefore, synaptic effects of dopamine can be seen as a change in the DTT of the striatum.

Effect of GPe induced disinhibition of FSI activity on the DTT

Despite being only a small portion of the neuronal population in the striatum, they can strongly modulate MSNs activity due to their high degree of convergent and divergent connections to the MSNs. Moreover, they form synaptic connections near or at MSN’s somata as opposed to glutamergic synapses or local MSN axon collaterals, which usually connect to spines or dendritic shafts [30]. Furthermore, recent data suggests that FSIs preferentially inhibit direct pathway MSNs (D1) as compared to indirect pathway MSNs (D2) [16].

For instance, Fig 4A showed that for small ACTX, D1 MSNs fire at higher rate than D2 MSNs at low inputs, however an increase in FSI firing rate would inhibit D1 MSNs more than D2 MSNs and, hence, reduce the DTT. By contrast, a reduction in the firing rate of FSIs will have an opposite effect. Thus, FSIs could play a vital role in maintaining the default state of the striatal circuit as a ‘No-Go’ state, which could be changed to a ‘Go’ state by exciting D1 MSNs more than D2 MSNs. Alternatively, a depletion of FSIs (e.g. as seen in patients with Tourette syndrome), a reduction in their firing rate or connectivity should reverse this scenario and enforce a ‘Go’ state as the default striatum state leading to im-pulsivity-related behavioral symptoms.

The effect of GPe on FSIs was modelled in our network model as a constant amount of inhibition on FSI firing rates. An increase in GPe activity (for instance due to an increase in STN activity) would effectively release D1 and D2 MSNs from the feedforward inhibition. However, the effect will be further amplified for D1 MSNs because these neurons receive more input from FSIs than D2 MSNs do. This can shift the bias from a ‘No-Go’ to ‘Go’ state for high cortical input rate (as shown in Fig 7B, dotted lines, for XCTX in the range (10—15 HZ)). The same effect of FSI activity can be observed in the Fig 3D. The same effect is shown for mean field model in Supplementary SS Fig. The temporal dynamics of such a modulation by GPe are shown in Supplementary S6 Fig for a mean field model. Thus, an increase in the activity of the FSIs would result in shifting the DTT to higher cortical rates. We argue that increased firing rates in STN might not only maintain a tight rein on GPi to prevent a premature response in a high conflict decision, as suggested by [5] , but also assist in resolving the conflict in the striatum via pallido-striatal back projections. In the dopamine depleted condition when the DTT is very small, the effect of the GPe back projection might not be able to modulate the DTT any further (Fig 7C). Similarly, tonically active interneuron (TANs) could exploit the FSI to D1 and D2 MSNs connectivities by modulating the strength of FSI synapses onto the respective MSNs [31], thereby, changing the DTT indirectly.

Effect of symmetrical FSI projections on the DTT

[14] and Gittis et al. [16]. Taverna et al. [14] showed that D2 MSNs inhibit D1 MSNs more strongly in terms of numbers of projections and synaptic strengths as compared to vice versa, whereas Gittis et al. [16] showed that in the LhX6-GFP transgenic mice, FSI has more projections to D1 than to D2, while the synaptic strengths of the projections to both MSN types are comparable. While it is unlikely that the preferential connectivity of FSIs to D1 MSNs is specific to the mouse strain (A. Gittis personal communication), another recent study by Planert et al. [15] , while supporting the results of Ta-verna et al. [14], did not report any significant preference in the connectivity of FSIs to D1 MSNs. Thus, even though Planert et al. drew their conclusion based on only a small number of neurons, the apparent inconsistency of the data in the two studies requires a reevaluation of the connectivity of FSIs to MSNs. We, therefore, analysed the striatal dynamics for equal innervation of MSNs by FSIs.

In this case, the inhibition from FSIs is considered same to be the same to both MSN subpopulations (I11: 2 I21: 2 I F). Therefore, the expression for AMSN in Eq 6 modifies to:

Given the nature of the MSNs interconnectivity (I12 + I 11 > 121 + I22) and the preference of FSI projections onto MSNs (IF < 1), Inpstr remains negative and increases (decreases) With M31 (XCTX). As a consequence, for a fixed value of M51, an increase in XCTX results in a sign reversal of AMSN, thereby, imposing a DTT. On the other hand, for a fixed value of XCTX, AMSN increases With increasing M51. However, note that if AMSN is positive, M31 enhances the magnitude of AMSN, Whereas When AMSN is negative, RPS] attenuates the magnitude of AMSN (Fig 8A). Thus, in a scenario Where XFSI and XCTX are comodulated a DTT can be observed (cf. Fig 8A, solid black line and Fig 8B). Symmetrical FSI projections in the multiplicative input scenario. In the multiplicative input scenario, When FSI projections to the D1 and D2 MSN subpopulations are equally strong (IIF 2 I21: 2 IF), the Eqs 13 and 14 reduce to: where commeff is the common FSI feedforward inhibition to both types of MSNs. In this scenario, an increase in the XCTX for a constant value of M31 enhances the activity of D2 MSNs and, thereby, reverses the sign of AMSN from positive to negative, imposing a DTT. By contrast, for a fixed value of XCTX, AMSN increases with increasing M31 reversing the sign of AMSN (Fig 8D) from negative to positive.

Fig 8D, solid black line and Fig 8E). This is because commefl is always positive due to higher effective inhibition to D1 (I12 + I11), and hence, the sum of Dlefic, D2617 and com-meficmonotonically increases as shown in Fig SD, SE and 8F.

In most cases, FSI activity is likely to be correlated with the cortical activity impinging onto the striatum. However, subtle changes in the correlation between the PSI and cortical activity can dramatically change striatal function. Moreover, these results also provide a possible eXperimental scenario to test predictions of our model. For instance, selective modulation and clamping of the activity of FSIs in a behavioural task could be used to determine whether the striatum operates in a multiplicative or additive mode, and whether FSIs have same the effective connectivity to both types of MSNs.

Discussion

Specifically, we focused on the firing rates of the D1 and D2 MSNs, because the relative balance of firing in these neurons is crucial for setting the stage for action-selection in the basal ganglia and reinforcement learning. Using both mean-field analysis and spiking network simulations of D1 and D2 MSNs we show that the experimentally observed asymmetry of recurrent and mutual connections among the two types of neurons creates a ‘decision transition thresh-old’(DTT) in the striatum to choose between the direct (‘Go’) and indirect (‘No-Go’) pathways. Based on this DTT, the striatum performs a threshold operation on the firing rates and spike correlations reaching the striatum from the neocortex and uses the outcome to bias the action-selection.

If ACTX is negative, the striatal dynamics are stuck in a state of constant bias towards D2 MSNs. By contrast, when ACTX is large and positive (more input to D1 MSNs), the striatal dynamics remain in a state of constant bias towards D1 MSNs. Thus, we predict that in order to have a flexible state in which both D1 and D2 MSNs could have higher firing rates and implement a preference towards either ‘Go’ or ‘No-Go’ actions, the D1 MSNs should receive only slightly higher cortical input than the D2 MSNs (cf. Fig 2 for details). A more exhaustive parameter search revealed many possibilities of implementing the thresholding operation. In an ‘additive’ input scenario, a DTT emerges for a constant firing rate of PSI and increasing cortical rates and vice versa (Fig 3A) together with a third, biologically more realistic possibility of correlated firing rates of cortical inputs and FSIs.

In addition, DTT can also be observed for increasing FSI firing rates and a constant cortical excitation (Fig 3D). This raises the issue of having an input-dependent DTT in the striatum, suggesting a multitude of mechanisms which could control both the existence as well as the value of the striatal DTT.

That is, in the space of input properties (firing rates and spike correlations) two regions exist for which either D1 MSNs or D2 MSNs have relatively higher output firing rates. These two regions are separated by the so called DTT. There is a third region, spanned by high input correlations (2 0.5) in the input parameter space, in which both neuron populations have equal firing rates. The notion of the DTT in the striatum naturally emerges as a consequence of the asymmetric connectivity between and among D1 and D2 MSNs (Figs 1 and 6). However, this does not imply that the DTT is fixed for the striatum. In fact, controlled modulation of the DTT is desired because it could provide a neuronal substrate for implementing notions like learning history, behavioral state, or motivation to perform a task. Furthermore, this view of striatal function can be extended to better understand neuronal mechanisms underlying various disorders of the basal ganglia. We identified several components of the BG network structure and activity that could potentially control and modulate the DTT.

We investigated this discrepancy between our observations and those described by Damodaran et al. [32] and tuned the leaky integrate and fire neurons to closely follow the input current and output firing rate (F-I) curves of the D1 and D2 MSNs by adjusting the passive properties of the neuron (i.e. membrane capacitance, resting membrane potential and time constant) according to the values described by Gertler et al. [13] (Supplementary 83A Fig). Using an asymmetric striatal network of these D1 and D2 MSNs, we show that a DTT exist for Scenario I (Supplementary 83B Fig). That is, different F-I curves and excitability of the D1 and D2 MSNs is not sufficient to abolish the DTT and create a balance of D1 and D2 output rates. Damodaran et al. [32] used a more complicated neuron model with some essential morphological features. To understand the reason for the discrepancy between our results and those described by Damodaran et al. [32], therefore, would require a more systematic analysis of the striatal connectivity and morphology and integrative properties of the MSNs, which is beyond the scope of the current manuscript.

Modulation of the DTT by cortical input correlations

To better understand the effect of input correlations we separately studied the role of correlations within and correlations between the pre-synaptic pools of the MSNs. Previously, we have shown that correlations between the pre-synaptic pools of MSNs (B) tend to reduce the signal-to-noise ratio of the striatal response to cortical inputs [12].

Here, we found that W correlations also affect the difference between D1 and D2 MSNs activity (AMSN) in a non-monotonic fashion (Fig 6). Interestingly, when input firing rates result in equal output f1r-ing rates of D1 and D2 MSNs, a modulation in W could break the symmetry and increase the firing rate of either one of the two MSN subpopulations. The value of within-correlations (Wopt) for which AMSN is maximal depends on the input firing rates and the value of between-correlations (B) (Fig 6).

In the regime W 2 Wopt, an increase in shared input (B) decreases the value of AMSN. However, in the regime W s Wopt, an increase in B increases the value of AMSN. Finally, for high values of input correlation lAMSNl is small and the two MSN populations have comparable firing rates, indicating the failure of the striatum to bias an action-selection process. Such an operating regime could halt the decision making process in high-conflict decision-making tasks. Indeed, we predict that in high conflict choice tasks, there should be an increase of correlation in the cortical activity reaching the striatum.

Modulation of the DTT by FSIs

For instance, an increase in FSI activity would reduce D1 MSNs activity and decrease the DTT whereas a decrease in FSI activity would tip the bias in favour of D1, and hence the ‘Go’ pathway. Such a decrease in FSI activity could be due to insufficient excitation from the corteX, or to high input from the pallido-striatal back backprojections. We suggest that selective modulation of the PSI activity, which could shift the DTT in the striatum from ‘No-Go’ to ‘Go’, could be a mechanism for decision making in a high-conflict situation. This functional concept may also help understand neuronal mechanisms underlying various behavioral phenomena. For instance, the arbitration of FSI via the STN-GPe pathway could eXplain longer reaction times, required to respond to a high-conflict decision. Similarly, increased impulsivity in PD patients with DBS could be explained by the fact that DBS increases GPe firing rates [33], which would decrease the activity of FSIs by pallido-striatal backprojections, resulting in lowering the DTT.

For example, we find that a preferential inhibition of D1 MSNs is required for a DTT to be present for ‘multiplicative’ input. By contrast, an ‘addi-tive’ scenario reveals a DTT, irrespective of the asymmetry of FSI projections onto MSNs. We showed that even in the ‘multiplicative’ input mode the striatum can exhibit a DTT, by clamping the FSIs rates, While increasing cortical input rates. We conclude that FSI activity plays a prominent role in determining the DTT in the striatum, underlining once more the functional importance of pallido-striatal backprojections to FSIs.

Modulation of the DTT by dopamine

This differential effect of dopamine provides a powerful mechanism to control the DTT and difference between the activities of the two neuronal subpopulations activities and, hence, the DTT. For instance, a dopamine induced increase in synaptic inputs to D1 MSNs shifts the DTT in the favour of D1 MSNs, and thereby increasing the cortical input regime for a ‘Go’ bias. This, in fact, could be a putative mechanism by which external administration of L-Dopa induces dyskinesia as a prominent side-effect. By contrast, a reduction in dopamine level reduces the regime for a ‘Go’ bias, while increasing the regime for a ‘No-Go’ bias. Thus, low dopamine results in higher firing rates in D2 MSNs, possibly manifesting in akinetic symptoms and inducing oscillations in the GPe-STN network [34].

At least at the population level, both D1 and D2 MSNs have been shown to be co-activated and neither of the two pathways were completely shutdown [8]. Consistent with these experimental findings, both D1 and D2 MSNs in our model are co-activated, however, the two populations differed in their firing rates. We note that here our goal is not to reproduce the experimental results of Cui et al. 2013 [8] which involved behaving mice performing self-paced lever pressing to retrieve the reward. Apart from the obvious simplicity of our network model, we think there might be an additional reason for this. In our model, we considered the striatum as unit representing a single action. Indeed, Cui et al. 2013 [8] suggest that a nonselective global deactivation of D2 MSNs in the striatum would abolish suppression of most unwanted motor programs and, hence, lead to hyperkinesia. It can be argued that our model represents such a global activation and deactivation of D1 and D2 in the striatum. The concept of a DTT, if extended to include multiple action channels, could be used to more closely address the results reported in [8].

Implications of asymmetry in the striatal circuit

Anatomical evidence, however, suggests that there is an inherent asymmetry in the striatal circuit in terms of mutual connections between D1 and D2 MSNs, as well as in feedforward projections from FSIs. In Fig 2, we showed that this asymmetry, for an appropriate activity regime, has the ability to reinforce the cortical bias (bias towards D1 at low cortical rates) or override it (bias towards D2 at high cortical rates). Taking cortical input correlations into consideration, these dynamics become even more complex, as is summarized in Fig 9. Low dopamine levels decrease the regime with D1 bias, while increasing the regime with D2 bias. Alternatively, an increase in dopamine levels increase the area for a D1 bias, while decreasing the area for a D2 bias. An increase in between-correlations in cortical input (B), however, pushes the contour towards the top left corner in Fig 9 (towards higher firing rates and lower values of within-correlation (W)), thereby enlarging the area of conflict (RHC). We therefore suggest, that the asymmetry in the striatal circuit is useful to cover different regimes of D1 and D2 bias, accommodating characteristic features of cortical input, including rates, input correlations and local dopamine levels.

In more realistic in vivo conditions, the connectivities of D1 and D2 MSNs may turn out to be different. In particular, if D1 and D2 MSNs would not differ in their connectivities and membrane properties, there is no reason to consider the striatum as a two population network and the role of action selection would need to be performed in the downstream nuclei. Such a single population striatum could initiate this process by ‘winner-less competition’ [11] or by enhancing the contrast of the input [12]. If the D1 and D2 MSNs would only differ in their membrane properties, then our analysis of asymmetric connectivity between the two neuron subpopulations would still be valid, because the effect of different neuron transfer functions could be reduced to a difference in effective connection strengths. In the extreme case, When D2 MSNs would receive more inhibition than D1 MSNs, our analysis remains valid, but we would have to reverse the interpretation.

Implications for the understanding brain disorders involving the basal ganglia

PD is thought to be a manifestation of dopamine depletion in the form of a domination of indirect and hyperdirect pathways, while weakening the direct pathway. Dyskinesia, on the other hand, is correlated with hyperactivity of the direct pathway and reduced activity of the indirect pathway [35]. The Tourette syndrome is correlated with decreased activity of FSIs [36]. This (dis) balance of striatal activity can be systematically modeled using the concept of the DTT. In fact, our model provides multiple neuronal mechanisms and a mechanistic understanding of several brain disorders by linking their specific pathophysiology to the DTT. Specifically, we explored three factors that may influence the DTT.

The DTT (basically, a threshold on the difference between the two) is a dynamic variable which may represent behavioral state, learning history and motivation level of the animal. Various mechanisms e.g. feedforward inhibition, dopamine, GPe backprojections, tonically active neurons exist that can modulate the DTT and thereby, provide the striatum with a rich computational repertoire.

Moreover, we included several factors that can modulate the DTT level by affecting chemical imbalances and changes in neuronal properties. Taken together, we provide a high-level functional model of the striatum, which can be easily linked to low-level properties. Typically, the striatum is included in large-scale functional models [5] , acting according to winner-take-all dynamics; an idea no longer supported by experiments [10]. Instead, we argue that our functional description of the striatum provides a more biologically realistic basis for large-scale functional models of basal ganglia function.

Model predictions and explanation of experimental data

The cognitive concept of conflict is described as a quantity that should increase with the absolute amount of activation and the number of competing representations [37]. In this condition, it is essential for the action selection system to stall the decision for further deliberation. A stalled decision in our model implies that the activity in the direct and indirect pathways are comparable i.e. AMSN % O (in our model AMSN % O for high cortical firing rates and correlations (Fig 6A)). Especially an increase in B expands the region in which AMSN % 0 (Fig 9). Therefore, we predict that high-conflict tasks should be associated with an increase in correlations in the activity of cortico-striatal projection neurons. Recent eXperimental data shows higher cortical firing rates in high-conflict situations [38]. According to our model this is not sufficient for decision stalling, because it would lead to a higher firing rate in D2 MSNs. Moreover, increased cortical rates should be accompanied with an increase in correlations. After stalling the decision, at some point a decision needs to be made in favor of either the direct or the indirect pathway. A rather simple way to resolve the conflict would be to change the cortical firing rates and correlations (cf. Fig 9). The pallido-striatal backprojections may also provide alternative mechanisms to resolve conflict and facilitate the decision making.

In our model, an increase in the steady state level of dopamine changes the DTT to higher cortical input rates, thereby increasing the region in which D1 MSNs can have their firing rate exceed those of D2 MSNs (Fig 7A). This observation could be used to explain the dyskinesia often reported in patients with PD on L-Dopa therapy (L-Dopa Induced Dyskinesia—LID). Undesired involuntary movements characteristic of LID can be attributed to the frequent undesired activation of D1 MSNs (direct pathway). This could be a direct consequence of a shift in DTT towards ‘Go’, which would result in a decrease in GPi activity during LID. Indeed, injection of dopamine agonists in MPTP treated nonhuman primates induced dyskinesia and showed a marked decrease in GPi firing rates [35].

According to our model, in a low dopamine state, there is a very small parameter regime in which D1 MSNs can have higher firing rate than D2 MSNs (Fig 7A). Therefore, we propose that in a low dopamine state even low-conflict tasks might require the arbitration by pallido-striatal backprojections via FSIs (as shown in Fig 7B and 85B), thereby increasing the average reaction time for voluntary actions in Parkinson’s disease patients.

In our model, a decrease in the steady state dopamine level results in an increase in the range of cortical input for which D2 MSNs have a higher firing rate than D1 MSNs. This might give rise to induced akinesia observed in dopamine-deplet-ed conditions. Such a bias towards ‘No-Go’ may have a twofold effect on the basal ganglia output. First, it would introduce an insufficient ‘Go’ bias in the GPi (firing rates do not decrease sufficiently) and second, in extreme dopamine depleted conditions it might introduce a strong inhibition on GPe, which can instigate beta band oscillations in the GPe-STN loop [34]. Furthermore, unlike in normal dopamine conditions, the STN-GPe loop is also insufficient to shift the balance of activity towards D1 MSNs (Figs 7C and 85C).

In several brain disorders, the number of FSIs and/or their connectivity to MSNs are altered. For instance, FSI count is reduced in Tourette syndrome and dystonia [39, 40]. Because FSIs preferentially connect to D1 MSNs, according to our model the loss of FSIs or their connectivity to MSNs would increase the activity in the direct pathway, compared to the indirect pathway (Fig 7B and 7C) and shift the DTT to higher values. This might underlie the symptoms of Tourette’s syndrome and dys-tonia. This explanation is partially supported by recent experimental findings which show that the symptoms of dystonia could be evoked in an animal model by a selective decrease of the firing rate of FSIs [36].

In our model, based on the experimental data, FSIs preferentially innervate D1 MSNs [16] in order to maintain the default state of bias as ‘No-Go’ under equal cortical drive, while allowing the cortical activity to switch the balance of D1 and D2 MSNs activity. Interestingly, in the mouse model of PD, three days after dopamine depletion, the FSIs doubled their connectivity to D2 MSNs as compared to D1 MSNs [41]. Our model suggests that dopamine depletion reduces the DTT such that the striatum mostly operates in a ‘No-Go’ state and, hence, cortical activity no longer suffices to switch the balance of D1 and D2 MSNs activity. In our model, an increase in the PSI to D2 MSNs connectivity will shift the DTT back in favour of D1 (Fig 4B). Therefore, we argue that the experimentally observed rewiring of the FSIs to D2 MSNs may be a compensatory mechanism to maintain striatal function in a low dopamine state. Increased impulsivity in PD patients with deep brain stimulation. Impulsivity could be considered as a too frequent activation of D1 MSNs or equivalently, a reduced DTT of the striatum. Deep brain stimulation (DBS) in the STN is known to increase GPe firing rates [33]. Such increased GPe firing rates would impose a sustained inhibition of FSIs, thereby shifting the DTT to higher cortical input rates and expanding the regime with a bias towards the ‘Go’ pathway, which may manifest itself as impulsivity at the behavioral level. Notice,that a similar trend was discussed in LID in Fig 7A. Therefore, we predict that striatal dynamics should show similar characteristics during impulsivity and LID. Indeed, there is evidence for the potential mechanistic overlap between behavioral ICDs (impulse control disorders) and motor (dyskine-sia) dopaminergic ramifications (cf. review in [42]).

Materials and Methods

The schematic of the striatal circuit is shown in Fig 1. The parameters considered for the rate model are included in Tables 1 and 2.

Mean field model

The mean population dynamics of the striatal circuit for scenario II were modelled using coupled nonlinear differential equations: The values of parameters used in Eqs 21 and 22 are listed in the Table 1. The nonlinear transfer function of the neuron 8(2) has a sigmoidal function of the form: The parameters In, 112, 121, 122 are calculated by considering the effective input received by a neuron in the population. This can be quantified as:

These parameters are listed in Table 2. This measure calculates the relative inhibitory strength of projections between D1 and D2: although the self connectivity of D1 is comparable to the number of projections it receives Table 3. Parameters for model neurons in network simulations. Parameter Stability of the striatal circuit. Linearizing around the fixed point With the Iacobian ma-triX: Where allows to analyze the stability by calculating the eigenvalues of the fixed point. We found the ei-genvalues to be real and negative for all values of the fixed points.

Network Simulations

The simulations were carried out in NEST [43] with networks of 4000 MSNs (2000 D1, 2000 D2) and 80 Fast Spiking Interneurons (FSI), consistent with the estimated ratio of MSN and F81 neurons [44, 45]. The neuron model parameters for FSIs and MSNs are listed in Table 3.

The cortical input parameters are summarized in Table 4. Connection probabilities for D1 MSNs, D2 MSNs and FSIs were taken from [14, 16] , also listed in Table 5. The leaky-integrate-and-fire (LIF) neuron model was used to simulate the neurons in the network with the subthreshold dynamics of the membrane potential Vx(t) described by: Table 4. Cortical input to striatal neurons.

recurrentIn the above equations, 15” "(1‘) describes the total cortical excitatory, re-current and feedforward inhibition to a neuron, C" is the membrane capacitance, Gx is the leak conductance, and Vrest is the resting membrane potential. When the membrane potential of the neuron reached Vth, a spike was elicited and the membrane potential was reset to Vrest for a refractory duration (trefz 2 ms.) An alpha function was used to model the excitatory synaptic input received by D 1, D2 and F81: where x 6 D1, D2, F81. The rise times for the excitatory inputs Tm was set to be identical for all neurons. In addition to excitatory input from the cortex, D1 and D2 MSNs received (self and mutual) recurrent inhibition as well as feedforward inhibition from the FSIs. The corresponding inhibitory conductance changes were modelled as:

The rise times rink for all inhibitory synaptic conductance transients were set to have identical values. Assuming the synaptic strengths as fixed, the total excitatory conductance G’e‘xc i(1‘) in a MSN 1' was given by by the inner sum, While the outer sum runs over all such excitatory synapses m in the set Kf‘ Ctx mn projecting onto neuron i. The set 1‘ represents the spike times of the excitatory neuron m. Viflh denote the reversal potentials of excitatory and inhibitory synaptic currents, respectively. Similarly, for a FSI neuron i, the excitatory conductance was: and the total synaptic current for the PSI neuron i was The parameter values for both MSNs and F813 in our network model are summarized in Table 3.

Generation of B and W correlations

The MIP model generates correlations by copying spikes from a spike train (the mother spike train) with a fixed probability (the copy probability, which determined the resulting correlation) to the individual spike trains. This process was implemented in the NEST simulator by introducing a synapse model that transmits spikes with a fixed probability: the lossy synapse. By making many convergent connections using the ‘lossy synapse’ we can mimic the random copying of spikes from the mother spike train to the children process. This way we controlled the correlations within (W) the pre-synaptic pool of

We started with a spike train M with a firing rate Rm. In the first layer of MIP we generated N spike trains (C1, C2, - - -, CN) with a fixed pairwise correlation (B’) and firing rate (Rm B’ ). In the second layer, each of these spike trains acted as a mother spike train to generate pairwise correlations in the pre-synaptic pools of individual striatal neurons. Each spike train C,- was used to generate correlated spike trains (Cl-0, Cil, - - -, Cim) that acted as the pre-synaptic pool of each striatal neuron with a pairwise correlation W and firing rate Rm B’ W (Fig 10). The effective correlation between the pre-synaptic pools of two neurons is then defined as B = B’ W. Thus, in this model we can control the correlations between and within the pre-synaptic pools by changing the copy probabilities in the first and second layer.

Relationship between B and W

B refers to the average correlation between the spiking activities of neurons in the pre-synaptic pools of two striatal neurons. The spike trains in each pre-synaptic pool are themselves correlated With a correlation coefficient W. For such pooled random variables, Bedenbaugh and Gerstein [22] derived the following relationship: where p12 is the correlation coefficient between the two pools of the pre-synaptic neurons, 11 is the size of each pre-synaptic pool. Because 0 3 p12 3 1, it follows that B s W.

In this case, the average correlations between the two pools will be small because each spike train has only copy in the other pool, while being uncorrelated with all others. However, when W > 0, more such pairs with nonzero correlation will occur, increasing the value of B. Because, B s W, the difference between the firing rates of the two MSN neuronal population in Fig 6 is estimated for the lower triangular region. This relationship is explicitly shown also in Fig 10B, where the measured input correlation between the pools (B) is shown to be less than W.

Model limitations

But ideally, there are multiple channels of competing motor programs which can be envisioned as mutually competing DTT’s. Thus, future work will include the extension of this model to incorporate multiple action channels. Secondly, the effect of input correlations is yet not formulated into the analytical framework. This is a nontrivial problem, due to the nonlinear transfer of input correlations into output firing rates, specifically if the output neurons are recurrently interconnected. Thirdly, our model explores the striatal dynamics at steady states, since the model showed stable fixed points. However, it is possible that in the presence of additional phenomena, such as plasticity of cortico-striatal or striato-striatal projections or intrinsic plasticity of the striatal neurons etc., the model may show interesting transient dynamics. Recent work proposed the existence of a ‘decision threshold’ in a different context [50]. They considered the ‘decision threshold’ as a function of the RT (reaction time) distribution of the rats in a decision task pertaining to the transient dynamics of the basal ganglia. A future extension to this work can explore the effect of DTT due to the asymmetry in the striatal connectivity on the distance to the ‘decision threshold’ in [50], i.e. the effect of DTT on the RTs of a decision. And lastly, though our model abides by many experimental observations, it contradicts one experimental result, to the best of our knowledge. In [51] , it was observed that FSIs increased their activity around a choice execution. In our model, a choice execution will correspond to the activation of the D1 pathway, which requires reduced activity in the FSIs, in contrast to the recordings in [51]. We think this contradiction is a consequence of our model of the striatum considering only one channel in the ‘Go’ and ‘No-Go’ pathways.

Supporting Information

(A) Cortical input rate applied to the striatal mean field model. The instantaneous input rate was choosen from a low pass filtered noise (time constant 2 5 msec) and standard deviation of :05 around the mean. The mean of the input rate was varied every 500 msec while the standard deviation remained fixed. (B) To mimic a noisy ACTX we provided additional noisy inputs to the D1 and D2 MSNs (Scenario 1). D1 MSNs received additional input with a mean of 1 HZ and standard deviation 0 — A — ctxl. (C) Similarly, the D2 MSNs also received low pass filtered noisy input with 0 HZ mean and standard deviation 0 — A — ctx2. Both a — A — ctxl and a — A — ctx2 were considered free variables here. (D) The instantaneous firing rates of D1 and D2 MSNs in response to the noisy cortical input and ACTX. A DTT can be observed in spite of noisy fluctuations. (E) To quantify the effect of noisy Actx, on the DTT we measured the areas around the crossover of the D1 and D2 MSNs firing rates. This is the area under the curve (yellow portion) for the time interval [t* — A,t*] (pre-DTT) and [t*,t* + A] (post-DTT), where t* is the time when DTT occurs. Because the gain of the D1 and D2 MSNs is different in pre-DTT (where km > km) and post-DTT regions (where km < km) we separately measured the area in the pre-DTT and post-DTT regimes. We normalised the pre-DTT and post-DTT areas with the areas measured for zero noise case. Here, we show the normalised pre-DTT and post-DTT areas for different values of a — A — ctxl and a — A — ctx2. As expected the increase in the input noise progressively decreases the pre-DTT area. Nevertheless, we can reliably measure DTT for large fluctuations compared to the mean. (TIFF)

(AF) Striatal parameters, A-ctxl, A-ctx2, I11, I22, I21, I12 were varied by up to i100% around their means as specified in the model. It is difficult to visualise this six dimensional space. We use two different Visualisations. First, we show Whether for i 100% variation in one of the parameters it is possible to find a DTT for any set of values of all other parameters (again Within :100% of their means). The “Solution” refers to When a DTT is observed for at least one parameter set, given a specific value of the parameter for Which robustness is tested. The values of parameters for Which no combination of other parameters yielded a DTT were marked as “No solution”. It can be observed that a DTT can be found for nearly all values of all parameters (AF). Only for very weak values of I 11 and I22 we could not find a combination of other parameters that would yield a DTT. (G—I) Next, we show the existence DTT for pairs of striatum network parameters. Each dot shows the existence of a DTT. Not only do these solutions follow the relational constraints as described in Taverna et al. (2010) and the model (e.g. I12 > I21), but also the centroids (marker With red asterisk) of these clusters lie very close to the values used in the model. This can be verified in the mean values shown in (AF). This indicates that these values are indeed a robust combination of striatal parameter values. (TIFF)

(A) F-I curves of D1 and D2 MSNs. Passive neuron properties were tuned to match the F-1 curves of the D1 and D2 MSNs to match With the experimental measurement of the F-1 curves of these neurons (shown in Gertler et al. 2008). (B) Firing rate of the D1 and D2 MSNs as a function of cortical input rate. D1 MSNs received extra input according to the scenario 1. D1 and D2 MSNs have different F-I curves shown in the panel A. These results show that a DTT exists in the striatum even When D1 and D2 MSNs have different F-I curves. The network connectivity is same as described in Table 5. (TIFF)

The blue and red regions mark the inputs for which D2 and D1 MSNs have higher firing rates. White regions show the inputs for which there is no difference between the firing rate of the two MSN populations. From these figures it is clear that B’ results in the increasing the region of high conflict.

(A) Fixed points for D1 and D2 MSNs plotted for different levels of dopamine. Darker shades of blue (red) correspond to D1 (D2) MSN activity for higher levels of dopamine. For lower than normal levels of dopamine, the DTT shifts to the left (from m 11 Hz to z 6 Hz). This decreases the regime with the bias towards D1 MSNs. For higher than normal levels of dopamine, the DTT shifts to right (from m 11 to % 19Hz). This in turn, increases the regime with a bias towards D1 MSNs. (B) and (C) same as Fig 7 in the main text but calculated for mean field equations. (TIFF)

(A) GPe activity is represented as an inhibition on FSI firing rates.

Dashed line represent the 95% confidence interval. (C,F) Inhibition on FSIs is unable to resolve the bias in favour of D1 in dopamine depletion conditions. (TIFF)

Acknowledgments

We thank Robert Schmidt, Alejandro Bujan, Ioannis Vlachos, Marko Filipovic, Jeanette Hellg-ren-Kotaleski and Stefan Rotter for their helpful comments on different versions of this manuscript.

Author Contributions

Performed the experiments: IB. Analyzed the data: IB AK. Contributed reagents/materials/ analysis tools: AK IB. Wrote the paper: AK IB AA.