Computational models of ion channels represent the building blocks of conductance-based, biologically inspired models of neurons and neural networks. Ion channels are still widely modelled by means of the formalism developed by the seminal work of Hodgkin and Huxley (HH), although the electrophysiological features of the channels are currently known to be better fitted by means of kinetic Markov-type models. The present study is aimed at showing why simplified Markov-type kinetic models are more suitable for ion channels modelling as compared to HH ones, and how a manual optimization process can be rationally carried out for both. Previously published experimental data of an illustrative ion channel (Na(V)1.5) are exploited to develop a step by step optimization of the two models in close comparison. A conflicting practical limitation is recognized for the HH model, which only supplies one parameter to model two distinct electrophysiological behaviours. In addition, a step by step procedure is provided to correctly optimize the kinetic Markov-type model. Simplified Markov-type kinetic models are currently the best option to closely approximate the known complexity of the macroscopic currents of ion channels. Their optimization can be achieved through a rationally guided procedure, and allows to obtain models with a computational burden that is comparable with HH models one.

The basal ganglia play an important role in decision making and selection of action primarily based on input from cortex, thalamus, and the dopamine system. Their main input structure, striatum, is central to this process. It consists of two types of projection neurons, together representing 95% of the neurons, and 5% of interneurons, among which are the cholinergic, fast-spiking, and low threshold-spiking subtypes. The membrane properties, somadendritic shape, and intrastriatal and extrastriatal synaptic interactions of these neurons are quite well described in the mouse, and therefore they can be simulated in sufficient detail to capture their intrinsic properties, as well as the connectivity. We focus on simulation at the striatal cellular/microcircuit level, in which the molecular/subcellular and systems levels meet. We present a nearly full-scale model of the mouse striatum using available data on synaptic connectivity, cellular morphology, and electrophysiological properties to create a microcircuit mimicking the real network. A striatal volume is populated with reconstructed neuronal morphologies with appropriate cell densities, and then we connect neurons together based on appositions between neurites as possible synapses and constrain them further with available connectivity data. Moreover, we simulate a subset of the striatum involving 10,000 neurons, with input from cortex, thalamus, and the dopamine system, as a proof of principle. Simulation at this biological scale should serve as an invaluable tool to understand the mode of operation of this complex structure. This platform will be updated with new data and expanded to simulate the entire striatum.

The striatum is the main input stage of the basal ganglia receiving extrinsic input from cortex and thalamus. The striatal projection neurons (SPN) constitute 95% of the neurons in the striatum in mice while the remaining 5% are cholinergic and GABAergic interneurons. The cholinergic (ChIN) and low-threshold spiking interneurons (LTS) are spontaneously active and form a striatal subnetwork involved in salience detection and goal-directed learning. Activation of ChINs has been shown to inhibit LTS via muscarinic receptor type 4 (M4R) and LTS in turn can modulate ChINs via nitric oxide (NO) causing a prolonged depolarization. Thalamic input prefentially excites ChINs, whereas input from motor cortex favours LTS, but can also excite ChINs. This varying extrinsic input with intrinsic reciprocal, yet opposing, effects raises the possibility of a slow input-dependent modulatory subnetwork. Here, we simulate this subnetwork using multicompartmental neuron models that incorporate data regarding known ion channels and detailed morphological reconstructions. The modelled connections replicate the experimental data on muscarinic (M4R) and nitric oxide modulation onto LTS and ChIN, respectively, and capture their physiological interaction. Finally, we show that the cortical and thalamic inputs triggering the opposing modulation within the network induce periods of increased and decreased spiking activity in ChINs and LTS. This could provide different temporal windows for selective modulation by acetylcholine and nitric oxide, and the possibility of interaction with the wider striatal microcircuit. 

The majority of excitatory synapses in the brain uses glutamate as neurotransmitter, and the synaptic transmission is primarily mediated by AMPA and NMDA receptors in postsynaptic neurons. Here, we present data-driven models of the postsynaptic currents of these receptors in excitatory synapses in mouse striatum. It is common to fit two decay time constants to the decay phases of the current profiles but then compute a single weighted mean time constant to describe them. We have shown that this approach does not lead to an improvement in the fitting, and, hence, we present a new model based on the use of both the fast and slow time constants and a numerical calculation of the peak time using Newton's method. Our framework allows for a more accurate description of the current profiles without needing extra data and without overburdening the comptuational costs. The user-friendliness of the method, here implemented in Python, makes it easily applicable to other data sets.

Striatum, the input stage of the basal ganglia, is important for sensory-motor integration, initiation and selection of behaviour, as well as reward learning. Striatum receives glutamatergic inputs from mainly cortex and thalamus. In rodents, the striatal projection neurons (SPNs), giving rise to the direct and the indirect pathway (dSPNs and iSPNs, respectively), account for 95% of the neurons and the remaining 5% are GABAergic and cholinergic interneurons. Interneuron axon terminals as well as local dSPN and iSPN axon collaterals form an intricate striatal network. Following chronic dopamine depletion as in Parkinson’s disease (PD), both morphological and electrophysiological striatal neuronal features are altered. Our goal with this \textit{in-silico} study is twofold: a) to predict and quantify how the intrastriatal network connectivity structure becomes altered as a consequence of the morphological changes reported at the single neuron level, and b) to investigate how the effective glutamatergic drive to the SPNs would need to be altered to account for the activity level seen in SPNs during PD. In summary we find that the richness of the connectivity motifs is significantly decreased during PD, while at the same time a substantial enhancement of the effective glutamatergic drive to striatum is present.