Modeling the effects of dopamine on network synchronization (Komek et al. 2012)

 Download zip file 
Help downloading and running models
Dopamine modulates cortical circuit activity in part through its actions on GABAergic interneurons, including increasing the excitability of fast-spiking interneurons. Though such effects have been demonstrated in single cells, there are no studies that examine how such mechanisms may lead to the effects of dopamine at a neural network level. In this study, we investigated the effects of dopamine on synchronization in two simulated neural networks; one biophysical model composed of Wang-Buzsaki neurons and a reduced model with theta neurons. In both models, we show that parametrically varying the levels of dopamine, modeled through the changes in the excitability of interneurons, reveals an inverted-U shaped relationship, with low gamma band power at both low and high dopamine levels and optimal synchronization at intermediate levels. Moreover, such a relationship holds when the external input is both tonic and periodic at gamma band range. Together, our results indicate that dopamine can modulate cortical gamma band synchrony in an inverted-U fashion and that the physiologic effects of dopamine on single fast-spiking interneurons can give rise to such non-monotonic effects at the network level.
1 . K├Âmek K, Bard Ermentrout G, Walker CP, Cho RY (2012) Dopamine and gamma band synchrony in schizophrenia--insights from computational and empirical studies. Eur J Neurosci 36:2146-55 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Realistic Network;
Brain Region(s)/Organism:
Cell Type(s): Abstract Wang-Buzsaki neuron; Abstract theta neuron;
Gap Junctions:
Transmitter(s): Dopamine;
Simulation Environment: XPP;
Model Concept(s): Synchronization;
Implementer(s): Ermentrout, Bard []; Komek, Kubra [kkomek at];
Search NeuronDB for information about:  Dopamine;
# Theta neuron network with 50e and 20i cells with all-to-all
# connectivity and tonic input to both the e- and i-cells
# Parameters
par I_e=1, I_i=0
par gei=.4,gii=0.15,gee=0.1, gie=2
par tauz=50,gz=.2
par sige=0.5, sigi=.02
par taue=1, taui=3
#ODEs for e-cells, i-cells and adaptation
# Synapses
global 1 x[0..49]-pi {x[j]=-pi}
global 1 y[0..19]-pi {y[j]=-pi}
# Initial Conditions
wiener we[0..49]
wiener wi[0..19]
par b=100
aux swgt=.8*se+.2*si
@ meth=euler,
@ total=1400,trans=400,maxstor=1000000
@ dt=0.02,nOut=10
@ xp=t,yp=se,xlo=600,xhi=1200,ylo=0,yhi=.5

Loading data, please wait...