# BUILDING NEUROCOMPUTATIONAL MODELS AT DIFFERENT LEVELS FOR BASAL GANGLIA CIRCUIT.

Abstract: The target of computational neuroscience studies can be considered in two-folds: understanding the connection between the physiology and functional aspects of the brain to develop new approaches for diagnosing and treatment of neurological disorders and behavioral deficits and understanding mind and consciousness to develop new intelligent technologies. The methods and approaches used in computational neuroscience have to overcome the complexity of the system in all aspects. So, different methods and approaches are developed for different scales not only for observing the phenomena, but also for modeling. In this paper, an approach is proposed to build a connection between different levels of modeling. A simple, linear system will be shown to give an understanding of the working principle of basal ganglia circuit which is modeled with a detailed spiking neural network approach. First, spiking neural network of basal ganglia circuit will be introduced and the role of dopamine on its functioning will be shown; then a simple linear system model will be given, and the relation between two models will be explained. The aim of this work is to show that even a simple model which is not sufficient for detailed understanding of the neuronal process, could give a coarse understanding of a complex phenomenon. Such simple models could be used as a starting point in building complex models and also can be benefited for implementing intelligent technologies.Keywords: Basal ganglia circuit, spiking neural network, mass model, dopamine.

1. Introduction

In computational neuroscience literature, there are numerous models of neurons and neural structures at different levels. One reason of this diversity is the collection of data at different levels using different measurement tools and methods. While based on single neuron measurements, it is possible to obtain data to build a detailed model of a neuronal behavior based on the role of ion channels [1,2], it is also possible to pinpoint the regions of the brain that are active during a task by fMRI (functional magnetic resonance imaging) and obtain neural field model which can mimic the collective activity of neurons at a specific region during a specific task [3]. The resolution of measurements also depend on the scale. While, temporal resolution of single neuron measurements and EEG (electroencephalogram) / LFP (local field potential) are better, spatial resolution of fMRI is superior to other techniques. Thus the models corresponding to different levels have to cope with all these different scales and the dynamics of the brain is either modeled by a set of nonlinear, ordinary differential equations or partial differential equations. Besides these, there are hybrid models, where some structures are modeled at neuron level; others are modeled as mass model [4,5]

Of course the role of different aspects on neurological disorders and diseases is another reason of this diversity. While mutual activity of neurons is responsible for some processes and the malfunctioning in their collective behavior give rise to deficits, the activity of ion channels and the concentration of ions and neurotransmitters are important in other cases. So models differ as they target all these different phenomena at different levels. This variety of models is needed since all provide information necessary to understand the complexity of cognitive processes and neurological disorders and diseases. In computational neuroscience, the level of the model built has to be decided considering the experimental results to be used and the cognitive process dealt with. In some cases, models at different levels should be considered together to have a better understanding.

There is another aspect which should be noted for the models in computational neuroscience other than building the connection between the physiology and functional aspects of the brain. As pointed out in abstract, computational neuroscience also focuses on understanding mind and consciousness to develop new intelligent technologies. For this aspect, the simplicity of the models is crucial, since implementation on a hardware and real time applications is possible only if the computational burden is manageable. Even though there are some attempts to develop special hardware for neural structures as SpiNNaker [6], neuromorphic hardware as Neurogrid [7], the scale of model should be kept small to use the models in mobile and robotics applications [4]. Thus, while modeling the behavior of a group of spiking neurons, the number of neurons in the model is not same as the number of neurons in the neuronal structure to be modeled, but scaled to a number that is capable to display the behavior. Also, the membrane potential of a neuron is modeled by either by first order dynamical systems as integrate and fire models [8,9] or second order dynamical systems as Izhikevich model [10], even though a detailed model could be obtained by adding ion channel dynamics to Hodgkin-Huxley model [11]. Thus, in neurorobotics applications and in developing new learning rules simple models are preferred, rather than detailed models.

In this paper, the objective is whether it is possible to foresee the behavior of a complicated computational model by a simple one. If this is fulfilled than it would be possible to build a connection between different levels of modeling and a tool can be developed to ease detailed modeling. Since models at each level are versatile as they point different aspects of the neural phenomenon, the aim is not to replace a detailed model by a simple one, but to use a coarse approach to understand a complex but detailed one. To show the possibility of such an approach basal ganglia network will be considered, and it will be modeled by spiking neural models and by a simple linear system.

Modeling basal ganglia network has been considered in computational neuroscience literature extensively [9,11-16] due to its role in voluntarily action selection, reward related learning and in neurological deficits and diseases as Parkinson's disease, Huntington's disease, and in behavioral deficits as addiction. Especially, models of basal ganglia network are developed to understand deep brain stimulation [2]. In recent years more attention is paid to the role of basal ganglia circuits in high level cognitive processes as decision making [17], substance dependence [18,19].

In the following section, a brief introduction to the basal ganglia circuits will be given, especially focusing on the role of dopamine in action selection. Then, the proposed spiking neuron model will be introduced and a simple linear mass model will be given. The simulation results obtained using BRIAN simulator and XPPAUT will be given and these results will be discussed. It will be shown that a connection between the the firing rate of the spiking neuron model and dynamic behavior of simple linear model can be drawn.

2. Basal Ganglia Circuit

Basal ganglia circuits proposed to have important role in motor activation and cognitive processes [20] especially their role in reward based learning and decision making pointed in various works [17,21-23]. Impairment of basal ganglia circuits manifest deficits in motor actions observed in neurodegenerative diseases such as Parkinson's and Huntington's disease, and also cause behavioral deficits observed in attention deficit hyperactivity disorder (ADHD), obsessive-compulsive disorder (OCD) and addiction [19,24-27]. These behavioral disorders and motor movement disorders are treated by deep brain stimulation (DBS), a well-known treatment of Parkinson's disease [28,29]. The role of basal ganglia in psychiatric disorders is considered more recently [30,31], and their treatment by DBS makes basal ganglia a target for functional and restorative neurosurgery [32,33].

In computational neuroscience, modeling basal ganglia circuits is an important subject and many models are proposed. While mostly focus more on the role of basal ganglia circuits in voluntary movement and action selection [9,11-16], there are also models for reinforcement learning [34-42]. Some models deal with the malfunctioning of basal ganglia circuits and how treatments can be developed [2,43-46]. Most of these work except [13,42,45] focus on simple spiking neuron networks. Here, both mass model simpler than the ones in [13,42,45] and spiking neuron networks will be considered.

Striatum is considered as the input structure of basal ganglia and it is together with Subthalamic nucleus (STN), Globus pallidus (internal(GPi), external(GPe) and ventral pallidum) and substantia nigra (pars compacta and pars reticulata) [47,48] form the direct, indirect and hyper-direct pathways of cortico- striatal circuit [20,24]. Normally, direct and indirect pathways are at an equilibrium state. Little perturbations on the output of basal ganglia circuit which correspond to the GPi/SNr (substantia nigra pars reticulata), result in the selection of an action. The role of the hyper-direct pathway is to perform the fine tuning between several possible output choices which are conducted by direct or indirect pathway [49]. Dopamine from substantia nigra pars compacta and ventral tegmental area modify the activation in basal ganglia circuit by acting on striatum.

The afferents of striatum are mainly cortical pyramidal neurons located in layer V and occasionally in layer III [50]. During motor activation, posterior putamen and the dorsolateral anterior putamen receive inputs from motor and motor association cortex [21,51].

Striatum is mostly composed of medium spiny neurons (MSNs) which comprise 80-95% [47] (90 - 95 % [48]) of striatum, remainder is mostly interneurons. Even though MSNs are structurally homogeneous they have different chemical properties and are classified according to their response to neurotransmitter dopamine (DA) [48]. The two most effective groups are MSN with D1 type and D2 type receptors. D2 type receptors are more abundant and they are considered to be promoting the selection of the latent reinforcers [52].

The stimulation of the D1 type DA receptors cause neuronal excitation while the stimulation of the D2 type receptors cause neuronal inhibition. Both D1 and D2 type receptors exist together on the membrane of a neuron and their net effect drives the neuronal output. While D1 type MSNs inhibit GPi neurons and form direct pathway, D2 type MSNs inhibit GPe neurons and form indirect pathway. The direct pathway promotes the action initiation and selection whereas indirect pathway prevents actions [21,53]. The role of DA on striatum behaviour is vital, many neurological diseases and disorders are due to malfunctioning of dopamine neurons in striatum [2,16,46].

The well-known and extensively studied basal ganglia action selection circuit [24] has two main pathways: direct and indirect. Both pathways start with the stimulus from cortex to the neostriatum (caudate and putamen - Str) and unite again at the output nucleus of the basal ganglia, GPi/SNr.

The direct pathway is responsible for action selection while the indirect pathway is responsible for the inhibition of the unwanted actions. Direct pathway starts with the glutamatergic projections from the cortex to the MSNs of striatum. Some of the striatal neurons have direct GABAergic (gamma aminobutyric acid) projections to the GPi. The connection between GPi and thalamus (THL) is inhibitory while the connection from thalamus to cortex is excitatory. The net result of the direct pathway is the inhibition of the inhibition on thalamus so the motor cortical areas are stimulated which is called disinhibition in the neuroscience literature [47,48]. In the indirect pathway different type of striatal neurons have inhibitory GABAergic connections to GPe. GPe inhibits STN while STN stimulates GPi via glutamatergic connection. The net result of the indirect pathway is the disinhibition on the GPi so thalamus is inhibited and motor cortical areas are less stimulated.

In addition to the direct and indirect pathways, there is another excitatory connection between the cortex and the STN and this is called the hyper-direct pathway [49]. Via hyper-direct pathway cortical activity is transferred to GPi over STN and the striatal pathways are short-circuited. By the help of the excitation from the STN, already tonically active GPi inhibits thalamus stronger. For action selection, ventral oralis nuclei of thalamus is active [54] and its efferent is layer IV of supplementary motor area of cortex [21,48,54,55].

3. Computational Models of Basal Ganglia Circuit at Different Levels

In this section, two different levels of modeling basal ganglia circuit will be considered and a spiking neuron model based on Izhikevich type neurons will be introduced first. Then, a simple continuous time differential equation set will be used to define a mass model of basal ganglia circuit. Such mass models for basal ganglia circuit have been given previously [13,42,56,57], but all these were composed of discrete time, nonlinear dynamical systems. These two model focus on the effect of dopamine on action selection and the role of dopamine is modeled by a parameter.

3.1. Spiking neural network model

In this work, a computational model of basal ganglia-thalamocortical circuit for action selection which is shown in Figure 1 is first built using point neurons. The point neuron model used in forming these groups is Izhikevich neuron [58] and the equations governing the neuron model with the reset condition are given in Equation 1 and Equation 2, respectively. While forming the model the neural activation pattern of each component of basal ganglia circuit is considered and to model different patterns different parameters of Izhikevich neuron model are used. These parameter values are given in the upper part of Table 1.

V' = 0.04[v.sup.2] + 5v +140 - u + ge - gi (1) U' = a(bv - u)

[mathematical expression not reproducible] (2)

In Equation 1, [g.sub.e] and [g.sub.i] represents the excitatory and inhibitory connections to the neuron. As a neuron is connected to numerous neurons sum of excitatory and inhibitory neurons effects the behavior of the neuron [59]. The dynamics of synapses and the parameter values are given in Equation 3, Equation 4 and Table 2, respectively.

[mathematical expression not reproducible] (3)

Here index e corresponds to excitatory, and i corresponds to inhibitory connections.

[mathematical expression not reproducible] (4)

As the connection increases when a presynaptic neuron fires, [W.sub.j--k,x] denotes this increase in the weight from neuron j to neuron k either excitatory or inhibitory.

While cortex in the model has 200 excitatory pyramidal neurons which have regular spiking activity. Striatum model consists of three groups of neuron populations: D1 type medium spiny neurons (MSND1), D2 type medium spiny neurons (MSND2) and interneurons (IN).

MSND1 and MSND2 neuron populations are composed of 100 neurons and IN group has 25 neurons. MSND1 and MSND2 are modeled as regular spiking neurons like cortex excitatory neurons. The neurons in GPi, GPe, IN and STN are modeled as fast spiking neurons with initial values [v.sub.i] =-75, [u.sub.i] =-16, while cortex inhibitory neurons have as initial condition [v.sub.i] =-65, [u.sub.i] =-15. The connection weights and probabilities for the model in Figure 1 are given in the lower part of Table 1. The number of point neurons considered for all structures are given in Table 3.

In order to model the role of dopamine on action selection, the synaptic connections have to be changed with dopamine and here, this is accomplished as in [60,61]. Thus, the Equations given in Equation 3, Equation 4 will be modulated with dopamine as in Equation 5.

[mathematical expression not reproducible] (5)

Here, with these equations implemented to the model, the effect of Dopamine can be investigated by changing DA parameter. In order to show different effect of dopamine on MSND1 and MSND2, the parameter [alpha] = DA is used for MSND1 group and [alpha] = 1/DA is used for MSND2 group.

3.2. Mass model

A mass model equations for basal ganglia circuits are formed by linear differential equations given in Equation 6. This model is inspired by the firing rate results of spiking neuron model given in previous section and knowledge of state space behavior of linear dynamical systems.

The behavior of cortex areas crtx are modeled by a Heaviside function. Striatum is represented by two state variables [msn.sub.D1] and [msn.sub.D2], which have afferent exitatory connection from crtx weighted by dopamine level (DA) denoted by [w.sub.DA1] and [w.sub.DA2], for [msn.sub.D1] and [msn.sub.D2], respectively. Each of neural structures other than striatum is modeled by a single state variable, thus [gp.sub.e], [gp.sub.i], stn and thl are represented by a single dynamical variable and they all have afferent and efferent connections, corresponding to the circuit given in Figure 1. The role of dopamine on action selection is investigated by changing dopamine level DA from low to high levels where the values of DA are taken as 0.25, 0.5 and 0.75 for low, normal, and high level respectively [62,63].

[mathematical expression not reproducible]

[mathematical expression not reproducible]

4. Results and Conclusion

The simulations are done for spiking neural network model and mass model proposed in Section 3. In both cases the role of dopamine on action selection is investigated similarly. The levels of dopamine are changed through time and its effect on the results are shown with the activity of THL population. This activity which will be called THL activity is depicted with raster plot and firing rate for spiking neuron model and with curves obtained by the dynamical system for mass model.

4.1. Spiking neural network results

The simulations for spiking neural network model are carried in Python based simulation environment BRIAN [59]. During the simulations the level of dopamine is modified by changing parameter [alpha] in Equation 5. As the effect of dopamine on D1 and D2 type receptors is different, and [alpha] =DA for MSND1 group and [alpha] =1/DA for MSND2 group, to change the value of [alpha], DA is taken as 1 for normal level, 0.9 for low level and 1.1 for high level. This modification is done through time and the level of dopamine is taken to be at normal level for the first 500ms and then switched to high level till 1500ms and then switched to normal level till 2500ms, to low level till 3500ms and again to normal level till the simulation ended at 4000ms. This switching is done to investigate the effect of dopamine level change on the action selection through time. The simulation results are given for a randomly chosen single neuron from MSND1 and MSND2 population and synaptic activity in Figure 2, where the effect of dopamine on single neuron activity can be followed. As it can be followed from Figure 2, the neuron chosen from MSND1 population is more active when the level of dopamine is high, and the neuron chosen from MSND2 population is more active when the level of dopamine is low. This change in activity is also projected to synaptic dynamics and synaptic activity in both populations show difference as the membrane potentials.

The outcome of the activity of neuronal populations in basal ganglia circuit determines the activity of THL which can be interpreted as the result of action selection. Thus eventhough the single neuron activity of striatum MSN population reflect the effect of dopamine, the overall result of all this neuronal activity is revealed in THL, so the population activity of MSND1, MSND2 and THL are given in Figure 3, with raster plots and firing rates of neurons in MSND1, MSND2 and THL populations.

In time course of the change of dopamine level can be followed from firing rates easily while raster plot gives the information of neuronal population at neuron level. The firing rates inspired the idea of proposing a simple model for the interaction of neuronal population, where each population activity is denoted by a single dynamical variable as proposed in mass model.

4.2. Mass model results

The simulation for the mass model are done using XPPAUT, a tool developed for dynamical system analysis. In order to demonstrate the analogy between the firing rate results of spiking neuronal network and mass model, the time course of dynamical behavior of the variables corresponding to MSND1, MSND2 and THL populations are given in Figure 4. Here, again the dopamine level is switched similarly from normal to high, then to normal level followed by low level and finally to normal level by changing the parameter DA in Equation 6. The simulation interval is different as the it is scaled for the mass model, but the results shown in Figure 4, resembles the change in firing rates given in Figure 3.

4.3. Discussion of the results

As the simulation results given in Subsections 4.1 and 4.2 reveal, the overall activity of neuronal population can be followed from the mass model proposed as a simple linear system. Of course, with this coarse linear system approach, it is not possible to investigate the synaptic activity, or the activity in a neuronal population, but a general idea about the effect of dopamine can be followed easily. So building a simple model which consolidate the afferent and efferent connections between neuronal populations and the effect of neurotransmitter, would be informative to grasp the dynamics behind the neuronal activity. Furthermore, such a simple model can be versatile for the implementation of biologically inspired approaches and developing new learning rules as in [4,5].

5. References

[1] S. A. Prescott, S. Ratte, Y. De Koninck, and T. J. Sejnowski, "Pyramidal neurons switch from integrators in vitro to resonators under in vivo-like conditions," Journal of neurophysiology, vol. 100, no. 6, pp. 3030-3042, 2008.

[2] D. Terman, J. E. Rubin, A. C. Yew, and C. J. Wilson, "Activity patterns in a model for the subthalamopallidal network of the basal ganglia," The Journal of Neuroscience, vol. 22, no. 7, pp. 2963-2976, 2002.

[3] G. Deco, V. K. Jirsa, P. A. Robinson, M. Breakspear, and K. Friston, "The dynamic brain: From spiking neurons to neural masses and cortical fields," PLoS Comput Biol, vol. 4, no. 8, p. e1000092, 08 2008.

[4] E. Erjelik, R. Elibol, and N. S. Sengor, "A model on building and modifying the stimulus action association in the brain," in 2015 23nd Signal Processing and Communications Applications Conference (SIU), pp. 2533-2536, May 2015.

[5] E. Erjelik and N. S. Sengor, "A neurocomputational model implemented on humanoid robot for learning action selection," in 2015 International Joint Conference on Neural Networks (IJCNN), pp.1-6, July 2015.

[6] S. B. Furber, F. Galluppi, S. Temple, and L. A. Plana, "The spinnaker project," Proceedings of the IEEE, vol. 102, no. 5, pp. 652-665, May 2014.

[7] B. V. Benjamin, P. Gao, E. McQuinn, S. Choudhary, A. R. Chandrasekaran, J. M. Bussat, R. Alvarez-Icaza, J. V. Arthur, P. A. Merolla, and K. Boahen, "Neurogrid: A mixed-analog-digital multichip system for large-scale neural simulations," Proceedings of the IEEE, vol. 102, no. 5, pp. 699-716, May 2014.

[8] M. Dipoppa and B. S. Gutkin, "Flexible frequency control of cortical oscillations enables computations required for working memory," Proceedings of the National Academy of Sciences, vol. 110, no. 31, pp. 12\,828-12\,833, 2013.

[9] K. Gurney, T. J. Prescott, and P. Redgrave, "A computational model of action selection in the basal ganglia. i. a new functional anatomy," Biological Cybernetics, vol. 84, no. 6, pp. 401-410, 2001.

[10] E. M. Izhikevich and G. M. Edelman, "Large-scale model of mammalian thalamocortical systems," Proceedings of the national academy of sciences, vol. 105, no. 9, pp. 3593-3598, 2008.

[11] C. Yucelgen, B. Denizdurduran, S. Metin, R. Elibol, and N. S. Sengor, "A Biophysical Network Model Displaying the Role of Basal Ganglia Pathways in Action Selection", in Artificial Neural Networks and Machine Learning--ICANN 2012: 22nd International Conference on Artificial Neural Networks, Lausanne, Switzerland, September 11-14, 2012, Proceedings, Part I, Berlin, Heidelberg: Springer Berlin Heidelberg, pp. 177-184, 2012. doi="10.1007/978-3-642-33269-2_23"

[12] K. Doya, "What are the computations of the cerebellum, the basal ganglia and the cerebral cortex?" Neural Networks, vol. 12, no. 7-8, pp. 961-974, 1999.

[13] J. G. Taylor and N. R. Taylor, "Analysis of recurrent corticobasal ganglia-thalamic loops for working memory," Biological Cybernetics, vol. 82, no. 5, pp. 415-432, 2000.

[14] K. Gurney, T. J. Prescott, and P. Redgrave, "A computational model of action selection in the basal ganglia. ii. analysis and simulation of behaviour," Biological Cybernetics, vol. 84, no. 6, pp. 411-423, 2001.

[15] M. D. Humphries, R. D. Stewart, and K. N. Gurney, "A physiologically plausible model of action selection and oscillatory activity in the basal ganglia," The Journal of Neuroscience, vol. 26, no. 50, pp. 12\,921-12\,942, 2006.

[16] M. Guthrie, C. Myers, and M. Gluck, "A neurocomputational model of tonic and phasic dopamine in action selection: A comparison with cognitive deficits in parkinson's disease," Behavioural Brain Research, vol. 200, no. 1, pp. 48-59, 2009.

[17] B. W. Balleine, M. R. Delgado, and O. Hikosaka, "The role of the dorsal striatum in reward and decision-making," Journal of Neuroscience vol. 27, no. 31, pp. 8161-8165, 2007.

[18] Q. J. Huys, A. Beck, P. Dayan, and A. Heinz, "Neurobiological understanding and computational models of addictive behavior," Phenomenological Neuropsychiatry: Bridging the Clinic and Clinical Neuroscience: Springer: Science}, 2011.

[19] D. Belin, S. Jonkman, A. Dickinson, T. W. Robbins, and B. J. Everitt, "Parallel and interactive learning processes within the basal ganglia: Relevance for the understanding of addiction," Behavioural Brain Research, vol. 199, no. 1, pp. 89-102, 2009.

[20] G. Alexander, M. Crutcher, and M. DeLong, "Basal ganglia-thalamocortical circuits: parallel substrates for motor, oculomotor, "prefrontal" and "limbic" functions," Progress in brain research}, vol. 85, pp. 119-146, 1990.

[21] S. N. Haber, J. L. Fudge, and N. R. McFarland, "Striatonigrostriatal pathways in primates form an ascending spiral from the shell to the dorsolateral striatum," The Journal of Neuroscience, vol. 20, no. 6, pp. 2369-2382, 2000.

[22] S. N. Haber and B. Knutson, "The reward circuit: Linking primate anatomy and human imaging," Neuropsychopharmacology, vol. 35, no. 1, pp. 4-26, 2010.

[23] G. S. Berns and T. J. Sejnowski, "How the Basal Ganglia Make Decisions".in Neurobiology of Decision-Making Berlin, Heidelberg: Springer Berlin Heidelberg, pp. 101-113, 1996.

[24] M. R. DeLong and T. Wichmann, "Circuits and circuit disorders of the basal ganglia," Archives of Neurology, vol. 64, no. 1, pp. 20-24, 2007.

[25] M. H. Teicher, C. M. Anderson, A. Polcari, C. A. Glod, L. C. Maas, and P. F. Renshaw, "Functional deficits in basal ganglia of children with attention-deficit/hyperactivity disorder shown with functional magnetic resonance imaging relaxometry," Nature Medicine, vol. 6, no. 4, pp. 470-473, 2000.

[26] R. C. Chacko, M. A. Corbin, and R. G. Harper, "Acquired obsessive-compulsive disorder associated with basal ganglia lesions," The Journal of Neuropsychiatry and Clinical Neurosciences, vol. 12, no. 2, pp. 269-272, 2000.

[27] M. G. Murer, K. Y. Tseng, F. Kasanetz, M. Belluscio, and L. A. Riquelme, "Brain oscillations, medium spiny neurons, and dopamine," Cellular and Molecular Neurobiology, vol. 22, no. 5, pp. 611-632, 2002.

[28] A. L. Benabid, "Deep brain stimulation for parkinson's disease," Current Opinion in Neurobiology, vol. 13, no. 6, pp. 696-706, 2003.

[29] C. D. Cunha, S. L. Boschen, A. Gomez-A, E. K. Ross, W. S. Gibson, H.-K. Min, K. H. Lee, and C. D. Blaha, "Toward sophisticated basal ganglia neuromodulation: Review on basal ganglia deep brain stimulation," Neuroscience \& Biobehavioral Reviews, vol. 58, pp. 186-210, 2015.

[30] J. Sui, G. D. Pearlson, Y. Du, Q. Yu, T. R. Jones, J. Chen, T. Jiang, J. Bustillo, and V. D. Calhoun, "In search of multimodal neuroimaging biomarkers of cognitive deficits in schizophrenia," Biological Psychiatry, vol. 78, no. 11, pp. 794-804, 2015.

[31] L. A. Gunaydin and A. C. Kreitzer, "Cortico-Basal Ganglia Circuit Function in Psychiatric Disease," in Annual Review of Physiology, vol. 78, pp. 327-350 2016.

[32] H. S. Mayberg, A. M. Lozano, V. Voon, H. E. McNeely, D. Seminowicz, C. Hamani, J. M. Schwalb, and S. H. Kennedy, "Deep brain stimulation for treatment-resistant depression," Neuron, vol. 45, no. 5, pp. 651-660, 2005.

[33] B. H. Kopell and B. D. Greenberg, "Anatomy and physiology of the basal ganglia: Implications for DBS in psychiatry," Neuroscience & Biobehavioral Reviews, vol. 32, no. 3, pp. 408-422, 2008.

[34] W. Schultz, P. Dayan, and P. R. Montague, "A neural substrate of prediction and reward," Science, vol. 275, no. 5306, pp. 1593-1599, 1997.

[35] R. E. Suri and W. Schultz, "Learning of sequential movements by neural network model with dopamine-like reinforcement signal, "Experimental Brain Research, vol. 121, no. 3, pp. 350-354, 1998.

[36] R. Suri, J. Bargas, and M. Arbib, "Modeling functions of striatal dopamine modulation in learning and planning," Neuroscience, vol. 103, no. 1, pp. 65-85, 2001.

[37] P. Dayan and B. W. Balleine, "Reward, motivation, and reinforcement learning," Neuron, vol. 36, no. 2, pp. 285-298, 2002.

[38] D. Joel, Y. Niv, and E. Ruppin, "Actor-critic models of the basal ganglia: new anatomical and computational perspectives," Neural Networks, vol. 15, no. 4-6, pp. 535-547, 2002.

[39] P. R. Montague, S. E. Hyman, and J. D. Cohen, "Computational roles for dopamine in behavioural control," Nature, vol. 431, pp. 760-767, 2004.

[40] M. Haruno and M. Kawato, "Heterarchical reinforcement-learning model for integration of multiple cortico-striatal loops: fMRI examination in stimulus-action-reward association learning," Neural Networks}, vol. 19, no. 8, pp. 1242-1254, 2006.

[41] G. S. Berns and T. J. Sejnowski, "A model of basal ganglia function unifying reinforcement learning and action selection," in Joint Symposium on Neural Computation}, 1994.

[42] N. Sengor, O. Karabacak, and U. Steinmetz, "A computational model of cortico-striato-thalamic circuits in goal-directed behaviour," in Artificial Neural Networks - ICANN 2008, ser. Lecture Notes in Computer Science. Springer Berlin Heidelberg, vol. 5164, pp. 328-337, 2008.

[43] B. S. Gutkin, S. Dehaene, and J.-P. Changeux, "A neurocomputational hypothesis for nicotine addiction," Proceedings of the National Academy of Sciences of the United States of America, vol. 103, no. 4, pp. 1106-1111, 2006.

[44] S. H. Ahmed, M. Graupner, and B. Gutkin, "Computational approaches to the neurobiology of drug addiction," Pharmacopsychiatry, vol. 42, no. S144, pp. 144-152, 2009.

[45] S. Metin and N. S. Sengor, "From occasional choices to inevitable musts: A computational model of nicotine addiction," Intell. Neuroscience, vol. 2012, pp. 18:18-18:18, Jan. 2012.

[46] M. M. McCarthy, C. Moore-Kochlacs, X. Gu, E. S. Boyden, X. Han, and N. Kopell, "Striatal origin of the pathologic beta oscillations in Parkinson's disease," Proceedings of the National Academy of Sciences, vol. 108, no. 28, pp. 11\,620-11\,625, 2011.

[47] L. Squire, D. Berg, F. Bloom, S. du Lac, A. Ghosh, and N. Spitzer, Fundamental Neuroscience, 3rd ed., Elsevier, 2008.

[48] E. R. Kandel, J. H. Schwartz, T. M. Jessell, and S. Mack, Principles of neural science, New York, Chicago, San Francisco: McGraw-Hill Medical, 2013.

[49] M. J. Frank, "Hold your horses: A dynamic computational role for the subthalamic nucleus in decision making," in Neural Networks, vol. 19, no. 8, pp. 1120-1136, 2006.

[50] A. Rosell and J. M. Gimenez-Amaya, "Anatomical re-evaluation of the corticostriatal projections to the caudate nucleus: a retrograde labeling study in the cat," in Neuroscience Research, vol. 34, no. 4, pp. 257-269, 1999.

[51] E. Y. Choi, B. T. T. Yeo, and R. L. Buckner, "The organization of the human striatum estimated by intrinsic functional connectivity," Journal of Neurophysiology, vol. 108, no. 8, pp. 2242-2263, 2012.

[52] R. N. Cardinal, C. A. Winstanley, T. W. Robbins, and B. J. Everitt, "Limbic corticostriatal systems and delayed reinforcement," Annals of the New York Academy of Sciences}, vol. 1021, no. 1, pp. 33-50, 2004.

[53] G. Alexander and M. Crutcher, "Functional architecture of basal ganglia circuits: neural substrates of parallel processing," Trends Neurosci., vol. 13, no. 7, pp. 266-271, 1990.

[54] S. Sakai, "Cerebellar thalamic and thalamocortical projections," in Handbook of the Cerebellum and Cerebellar Disorders, M. Manto, J. Schmahmann, F. Rossi, D. Gruol, and N. Koibuchi, Eds. Springer Netherlands, pp.529-547, 2013.

[55] P. van Donkelaar, J. F. Stein, R. E. Passingham, and R. C. Miall, "Temporary inactivation in the primate motor thalamus during visually triggered and internally generated limb movements," Journal of Neurophysiology, vol. 83, no. 5, pp. 2780-2790, 2000.

[56] B. Denizdurduran and N. S. Sengor, "A realization of goal-directed behavior-implementing a robot model based on cortico-striato-thalamic circuits." in ICAART, pp. 289-294, 2012.

[57] O. Karabacak and N. S. Sengor, "A Computational Model for the Effect of Dopamine on Action Selection During Stroop Test", In: Kollias S.D., Stafylopatis A., Duch W., Oja E. (eds) Artificial Neural Networks--ICANN 2006. ICANN 2006. Lecture Notes in Computer Science, vol 4131, pp. 485-494 Springer, Berlin, Heidelberg, 2006.

[58] E. Izhikevich, "Simple model of spiking neurons," Neural Networks, IEEE Transactions on}, vol. 14, no. 6, pp. 1569-1572, Nov 2003.

[59] D. Goodman and R. Brette, "Brian: a simulator for spiking neural networks in python," Frontiers in Neuroinformatics, vol. 2, p. 5, 2008.

[60] R. Elibol and N. S. Sengor, "A computational model to investigate the effect of dopamine on neural synchronization in striatum," 2015 International Joint Conference on Neural Networks (IJCNN), Killarney, pp. 1-5, 2015. doi: 10.1109/LTCNN.2015.7280762

[61] R. Elibol and N. S. Sengor, "A Computational Model of Neural Synchronization in Striatum.", in Computational Neurology and Psychiatry Volume 6 of the series Springer Series in Bio-/Neuroinformatics pp 373-403, 2017. doi: 10.1007/978-3-319-49959-8_13.

[62] R. Elibol and N. S. Sengor, "Looking at the role of direct and indirect pathways in basal ganglia networks at different levels," BMC Neuroscience, vol. 16, no. 1, p. P225, 2015. doi: 10.1186/1471-2202-16-S1-P225.

[63] R. Elibol and N. S. Sengor, "Modeling basal ganglia circuits with mass model equations," 2016 Medical Technologies National Congress (TIPTEKNO), Antalya, Turkey, 2016. doi: 10.1109/TIPTEKNO.2016.7863131.

Rahmi Elibol received B. Sc. degrees from Karadeniz Technical University, Turkey, in 2004, and M. Sc. degrees from Istanbul Technical University, Turkey, in 2013. He joined Istanbul Technical University as research assistant in 2011, where he is a PhD student now. His research interests include nonlinear dynamical systems, mathematical modeling and computational neuroscience

Neslihan Serap Sengor received B. Sc., M. Sc. And Ph.D. degrees from Istanbul Technical University, Turkey, in 1985, 1988, 1995 respectively. She joined Istanbul Technical University as research assistant in 1986, where she is Professor now. She worked as visiting scientist in Helsinki University of Technology 2000-2001 and summer of 2006. She also visited Lincoln University, UK as research associate in the summer of 2015. Her research interests include nonlinear circuits and systems, neural networks.

Rahmi Elibol (1), Neslihan Serap SENGOR (1)

(1) Istanbul Technical University, Electronics and Communication Engineering, Istanbul, Turkiye rahmielibol@itu.edu.tr, sengorn@itu.edu.tr

Received on: 31.01.2017

Accepted on: 13.03.2017

Table 1. Izhikevich model parameters and connection weights. [logical not] is used for inhibitory and [right arrow] for excitatory connections. Parameters RS FS a 0.02/ms 0.1/ms b 0.25/ms 0.2/ms c -65mV -65mV d 8mV/ms 2mV/ms Connections [w.sub.j--k,x] Probabilities CRTX [right arrow] MSND1 0.75V/s 0.5 IN [logical not] MSND1 1V/s 0.5 CRTX [right arrow] MSND2 0.75V/s 0.5 IN [logical not] MSND2 1V/s 0.5 MSND1 [logical not] IN 1V/s 0.25 MSND2 [logical not] IN 1V/s 0.25 MSND2 [logical not] GPe 1V/s 0.25 STN [right arrow] GPe 1V/s 0.25 GPe [logical not] STN 1V/s 0.25 CRTX [right arrow] STN 1V/s 0.25 MSND1 [logical not] GPi 0.75V/s 0.25 GPe [logical not] GPi 0.75V/s 0.25 GPi [logical not] THL 1V/s 0.25 Table 2. Synaptic time constant and frequencies of Poisson groups. [[tau].sub.syn] CRTX GPe, THL STN GPi 10ms 100Hz 150Hz 200Hz 250Hz Table 3. Number of neurons in each neural population given in Figure 1. Neural Population # of neurons Behaviour CRTX 200 RS MSND1 100 RS MSND2 100 RS IN 25 FS GPi 100 FS GPe 100 FS STN 100 FS THL 100 FS

Printer friendly Cite/link Email Feedback | |

Author: | Elibol, Rahmi; Sengor, Neslihan Serap |
---|---|

Publication: | Istanbul University - Journal of Electrical & Electronics Engineering |

Article Type: | Report |

Date: | Jan 1, 2017 |

Words: | 5917 |

Previous Article: | The Role of Sports Participation in Hemispheric Dominance: Assessment by Electrodermal Activity Signals. |

Next Article: | EMOTION RECOGNITION VIA GALVANIC SKIN RESPONSE: COMPARISON OF MACHINE LEARNING ALGORITHMS AND FEATURE EXTRACTION METHODS. |

Topics: |