Synchronization by D4 dopamine receptor-mediated phospholipid methylation (Kuznetsova, Deth 2008)

 Download zip file 
Help downloading and running models
Accession:112968
"We describe a new molecular mechanism of dopamine-induced membrane protein modulation that can tune neuronal oscillation frequency to attention related gamma rhythm. This mechanism is based on the unique ability of D4 dopamine receptors (D4R) to carry out phospholipid methylation (PLM) that may affect the kinetics of ion channels. We show that by deceasing the inertia of the delayed rectifier potassium channel, a transition to 40 Hz oscillations can be achieved. ..."
Reference:
1 . Kuznetsova AY, Deth RC (2008) A model for modulation of neuronal synchronization by D4 dopamine receptor-mediated phospholipid methylation. J Comput Neurosci 24:314-29 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Realistic Network; Neuron or other electrically excitable cell; Synapse;
Brain Region(s)/Organism:
Cell Type(s): Neocortex V1 L6 pyramidal corticothalamic GLU cell; Neocortex V1 L2/6 pyramidal intratelencephalic GLU cell;
Channel(s): I Na,t; I T low threshold; I K; I K,Ca;
Gap Junctions:
Receptor(s):
Gene(s):
Transmitter(s): Dopamine;
Simulation Environment: XPP;
Model Concept(s): Activity Patterns; Oscillations; Synchronization; Simplified Models; Signaling pathways;
Implementer(s): Kuznetsova, Anna [anna.kuznetsova at utsa.edu];
Search NeuronDB for information about:  Neocortex V1 L6 pyramidal corticothalamic GLU cell; Neocortex V1 L2/6 pyramidal intratelencephalic GLU cell; I Na,t; I T low threshold; I K; I K,Ca; Dopamine;
# Kuznetsova and Deth, JCN, 2007
# 100-cell model:  80 E cells (Ve1, Ve2, ...Ve8) and 20 I (Vi1, Vi2) cells 
# all to all reciprocally connected 
# variables: V-voltage, R-K channel, C-Ca channel, H- Ca dependent K chanel, S,F - synapse 

init Ve[1...5]=-0.75
init Ve[6...10]=-0.71
init Ve[11...15]=-0.73
init Ve[16...25]=-0.72
init Ve[26...30]=-0.735
init Ve[31...40]=-0.715
init Ve[41...45]=-0.737
init Ve[46...50]=-0.718
init Ve[51...60]=-0.728
init Ve[61...70]=-0.74
init Ve[71...75]=-0.726
init Ve[76...80]=-0.713
init Vi[1...6]=-0.742
init Vi[7...12]=-0.714
init Vi[13...20]=-0.728
init Re[1...80]=0.26
init Ce[1...80]=0.1
init He[1...80]=0.1
init Ri[1...20]=0.26
init Ci[1...20]=0.1

par Ie=0.6
#continuos input in all E cells, canceled in I cells
par GCe=0.1, GCi=0.25
par GHe=4.0

par TRe[1...2]=2.1
par TRe[3...4]=2.12
par TRe[5...6]=2.3
par TRe[7...8]=2.14
par TRe[9...10]=2.2
par TRe[11...12]=2.25
par TRe[13...14]=2.29
par TRe[15...16]=2.15
par TRe[17...18]=2.18
par TRe[19...20]=2.2
par TRe[21...22]=2.24
par TRe[23...24]=2.15
par TRe[25...26]=2.25
par TRe[27...28]=2.11
par TRe[29...30]=2.28
par TRe[31...32]=2.18
par TRe[33...34]=2.22
par TRe[35...36]=2.28
par TRe[37...38]=2.3
par TRe[39...40]=2.23
par TRe[41...42]=2.13
par TRe[43...44]=2.22
par TRe[45...46]=2.14
par TRe[47...48]=2.26
par TRe[49...50]=2.14
par TRe[51...52]=2.24
par TRe[53...54]=2.28
par TRe[55...56]=2.3
par TRe[57...58]=2.27
par TRe[59...60]=2.17
par TRe[61...62]=2.11
par TRe[63...64]=2.14
par TRe[65...66]=2.28
par TRe[67...68]=2.205
par TRe[69...70]=2.26
par TRe[71...72]=2.16
par TRe[73...74]=2.29
par TRe[75...76]=2.12
par TRe[77...78]=2.26
par TRe[79...80]=2.18

# this is time constant for K channel in E cells, it is varied from 6 to 2

par GSee=0.07595
#6/79
par GSei=0.15
#6/40
par GSi=0.15789
#3/19=0.15789
#from I to E and I to I

par TRi[1..20]=1.5
# the same for in I cells, usualy it is not varied 

par TSe=2., TSi=8.
#par ESe=0, ESi=-0.75
#excit=0  inhib=-0.75
par W=-0.1

AlE=Se1+Se2+Se3+Se4+Se5+Se6+Se7+Se8+Se9+Se10+Se11+Se12+Se13+Se14+Se15+Se16+Se17+Se18+Se19+Se20 \\
+Se21+Se22+Se23+Se24+Se25+Se26+Se27+Se28+Se29+Se30+Se31+Se32+Se33+Se34+Se35+Se36+Se37+Se38+Se39+Se40 \\
+Se41+Se42+Se43+Se44+Se45+Se46+Se47+Se48+Se49+Se50+Se51+Se52+Se53+Se54+Se55+Se56+Se57+Se58+Se59+Se60 \\
+Se61+Se62+Se63+Se64+Se65+Se66+Se67+Se68+Se69+Se70+Se71+Se72+Se73+Se74+Se75+Se76+Se77+Se78+Se79+Se80

AlI=Si1+Si2+Si3+Si4+Si5+Si6+Si7+Si8+Si9+Si10+Si11+Si12+Si13+Si14+Si15+Si16+Si17+Si18+Si19+Si20

Fe[1..80]'=(1./TSe)*(-Fe[j]+heav(Ve[j]-W))
Se[1..80]'=(1./TSe)*(-Se[j]+Fe[j])

Fi[1..20]'=(1./TSi)*(-Fi[j]+heav(Vi[j]-W))
Si[1..20]'=(1./TSi)*(-Si[j]+Fi[j])

Ve[1..80]'=-Minf(Ve[j])*(Ve[j]-0.5)-26.*Re[j]*(Ve[j]+0.95)-GCe*Ce[j]*(Ve[j]-1.2)-GHe*He[j]*(Ve[j]+0.95)-GSi*(Ve[j]+0.75)*AlI-GSee*(Ve[j]+0.0)*(AlE-Se[j])+Ie
Re[1..80]'=(1./TRe[j])*(-Re[j]+Rinf(Ve[j]))
Ce[1..80]'=(1./14.)*(-Ce[j]+Cinf(Ve[j]))
He[1..80]'=(1./45.)*(-He[j]+3.*Ce[j])

Vi[1..20]'=-Minf(Vi[j])*(Vi[j]-0.5)-26.*Ri[j]*(Vi[j]+0.95)-GCi*Ci[j]*(Vi[j]-1.2)-GSi*(Vi[j]+0.75)*(AlI-Si[j])-GSei*(Vi[j]+0.0)*AlE
Ri[1..20]'=(1./TRi[j])*(-Ri[j]+Rinf(Vi[j]))
Ci[1..20]'=(1./14.)*(-Ci[j]+Cinf(Vi[j]))

Cinf(V)=8.*(V+0.725)^2
Minf(V)=17.8+47.6*V+33.8*V*V
Rinf(V)=1.24+3.7*V+3.2*V*V

aux n=(Ve1+Ve2+Ve3+Ve4+Ve5+Ve6+Ve7+Ve8+Ve9+Ve10+Ve11+Ve12+Ve13+Ve14+Ve15+Ve16+Ve17+Ve18+Ve19+Ve20\\
+Ve21+Ve22+Ve23+Ve24+Ve25+Ve26+Ve27+Ve28+Ve29+Ve30+Ve31+Ve32+Ve33+Ve34+Ve35+Ve36+Ve37+Ve38+Ve39+Ve40\\
+Ve41+Ve42+Ve43+Ve44+Ve45+Ve46+Ve47+Ve48+Ve49+Ve50+Ve51+Ve52+Ve53+Ve54+Ve55+Ve56+Ve57+Ve58+Ve59+Ve60\\
+Ve61+Ve62+Ve63+Ve64+Ve65+Ve66+Ve67+Ve68+Ve69+Ve70+Ve71+Ve72+Ve73+Ve74+Ve75+Ve76+Ve77+Ve78+Ve79+Ve80)/80.

#if(Ve1-W>0)then(Q(Ve1-W)=1)else(Q(Ve1-W)=0)

@ METHOD=stiff, TOLERANCE=.00001
@ MAXSTOR=400000, TOTAL=1000, XP=t,YP=n, BELL=0
@ xmin=0.0,xmax=1000,ymin=-1,ymax=0.5
@ DT=0.01, xlo=0.0,ylo=-1.0,xhi=1000,yhi=0.5,bound=30000

done

Loading data, please wait...