/*
Purkinje_reduced_PPR model_cvode.hoc
A hoc file for a cerebellar Purkinje neuron.
Adapted by SHERRY-ANN BROWN in "Virtual NEURON: a Strategy For Merged Biochemical and Electrophysiological Modeling"(Brown et al., 2010)
from original code written by MIYASHO Tsugumichi and MIYAKAWA Hiroyoshi.
This model runs on NEURON (version 3.2.3).
Description of the original version of this model can be found in
"Low-threshold potassium channels and a low-threshold calcium channel regulate
Ca2+ spike firing in the dendrites of cerebellar Purkinje neurons : a modeling study"(T. Miyasho et al., Brain Res 891, 2001, 106-115)
Morphology of the original cell used for the original model was taken from a paper by D.P.Shelton (Neuroscience 14, 1985, 111-131).
Channel descriptions are taken from a paper by E.De.Shutter and
J.M.Bower. (J. Neurophysiol., 71, 1994, 375-400)
*/
xopen("$(NEURONHOME)/lib/hoc/noload.hoc")
{load_proc("nrnmainmenu")}
/*
{
load_template("Inserter")
load_template("Electrode")
load_template("PointProcessManager")
load_template("PointProcessLoCaTcvodeor")
load_template("ShowMechanism")
}
*/
NumMainDendrite=1
NumSmoothDistalDendriteshort=1
NumSmoothDistalDendritelong=1
NumSpinyDistalDendrite=1
NumAdjacentDendrite=1
NumSpineNeck=1
NumSpine=1
NumIO1sm=1
NumIO1sp=1
NumIO2sm=1
NumIO2sp=1
NumIO3sm=1
NumIO3sp=1
NumIO4sm=1
NumIO4sp=1
NumIO5sp=1
NumIO6sp=1
create Soma, MainDendrite[NumMainDendrite], SmoothDistalDendriteshort[NumSmoothDistalDendriteshort], SmoothDistalDendritelong[NumSmoothDistalDendritelong], SpinyDistalDendrite[NumSpinyDistalDendrite], AdjacentDendrite[NumAdjacentDendrite], SpineNeck[NumSpineNeck], Spine[NumSpine], IO1sm[NumIO1sm], IO1sp[NumIO1sp], IO2sm[NumIO2sm], IO2sp[NumIO2sp], IO3sm[NumIO3sm], IO3sp[NumIO3sp], IO4sm[NumIO4sm], IO4sp[NumIO4sp], IO5sp[NumIO5sp], IO6sp[NumIO6sp]
access Soma
Soma connect MainDendrite(0),1
Soma {nseg=1 L=22 diam=22 }
MainDendrite {nseg=1 L=18 diam=4 }
MainDendrite connect IO1sm (0),1
IO1sm {nseg=1 L=140.5 diam=7.74 }
IO1sm connect IO1sp (0),1
IO1sp {nseg=1 L=95 diam=5.63 }
IO1sm connect IO2sm (0),1
IO2sm {nseg=1 L=82.5 diam=4.85 }
IO2sm connect IO2sp (0),1
IO2sp {nseg=1 L=159 diam=5.34 }
MainDendrite connect SmoothDistalDendriteshort (0),1
SmoothDistalDendriteshort {nseg=3 L=44 diam=5 }
SmoothDistalDendriteshort connect IO3sm (0),0.5
IO3sm {nseg=1 L=33 diam=2.63 }
IO3sm connect IO3sp (0),1
IO3sp {nseg=1 L=69 diam=3.67 }
SmoothDistalDendriteshort connect SmoothDistalDendritelong (0),1
SmoothDistalDendritelong {nseg=3 L=74 diam=5.3 }
SmoothDistalDendritelong connect IO4sm (0),0.5
IO4sm {nseg=1 L=134 diam=4.75 }
IO4sm connect IO4sp (0),1
IO4sp {nseg=1 L=204 diam=5.271633548 }
SmoothDistalDendritelong connect SpinyDistalDendrite (0),1
SpinyDistalDendrite {nseg=3 L=9.02 diam=2 }
SpinyDistalDendrite connect IO5sp (0),0.5
IO5sp {nseg=1 L=21 diam=1.918877022 }
SpinyDistalDendrite connect AdjacentDendrite (0),1
AdjacentDendrite {nseg=15 L=27.98 diam=1.4 }
AdjacentDendrite connect IO6sp (0),0.5
IO6sp {nseg=1 L=36 diam=1.86 }
AdjacentDendrite connect SpineNeck (0),1
SpineNeck {nseg=1 L=0.66 diam=0.2 }
SpineNeck connect Spine (0),1
Spine {nseg=1 L=0.57 diam=0.57 }
Soma {cm=0.8
insert Leak gl_Leak = 0.00006 el_Leak = -80
insert NaFcvode gnabar_NaFcvode = 10
insert NaPcvode
insert CaP2cvode cai = 0.00012723632 cao = 2.4 gcabar_CaP2cvode = 0.0005
insert CaTcvode cai = 0.00012723632 cao = 2.4 gcabar_CaTcvode = 0
insert CaEcvode cai = 0.00012723632 cao = 2.4 gcabar_CaEcvode = 0
insert Khhcvode
insert KMcvode
insert KAcvode
insert KDcvode gkbar_KDcvode = 0
insert Khcvode gkbar_Khcvode = 0.0005
insert cad taur_cad = 2 cainf_cad = 4e-5
kt_cad = 4e-5 kd_cad = 4e-5
}
MainDendrite {cm=0.8
insert Leak gl_Leak = 0.00021 el_Leak = -80
insert CaP2cvode cai = 0.0040600975 cao = 2.4 gcabar_CaP2cvode = 0.004
insert CaTcvode cai = 0.0040600975 cao = 2.4 gcabar_CaTcvode = 0.0015
insert CaEcvode cai = 0.0040600975 cao = 2.4 gcabar_CaEcvode = 0.008
insert Khhcvode gkbar_Khhcvode = 0.0006
insert KMcvode gkbar_KMcvode = 0.00001
insert KAcvode gkbar_KAcvode = 0.08
insert KDcvode gkbar_KDcvode = 0.09
insert KC3cvode gkbar_KC3cvode = 0.06
insert K23cvode
insert cad taur_cad = 2 cainf_cad = 4e-5
kt_cad = 4e-5 kd_cad = 4e-5
}
SmoothDistalDendriteshort {cm=0.8
insert Leak gl_Leak = 0.00021 el_Leak = -80
insert CaP2cvode cai = 0.0040600975 cao = 2.4 gcabar_CaP2cvode = 0.004
insert CaTcvode cai = 0.0040600975 cao = 2.4 gcabar_CaTcvode = 0.0015
insert CaEcvode cai = 0.0040600975 cao = 2.4 gcabar_CaEcvode = 0.008
insert Khhcvode gkbar_Khhcvode = 0.0006
insert KMcvode gkbar_KMcvode = 0.00001
insert KAcvode gkbar_KAcvode = 0.08
insert KDcvode gkbar_KDcvode = 0.09
insert KC3cvode gkbar_KC3cvode = 0.06
insert K23cvode
insert cad taur_cad = 2 cainf_cad = 4e-5
kt_cad = 4e-5 kd_cad = 4e-5
}
SmoothDistalDendritelong {cm=0.8
insert Leak gl_Leak = 0.00021 el_Leak = -80
insert CaP2cvode cai = 0.0040600975 cao = 2.4 gcabar_CaP2cvode = 0.004
insert CaTcvode cai = 0.0040600975 cao = 2.4 gcabar_CaTcvode = 0.0015
insert CaEcvode cai = 0.0040600975 cao = 2.4 gcabar_CaEcvode = 0.008
insert Khhcvode gkbar_Khhcvode = 0.0006
insert KMcvode gkbar_KMcvode = 0.00001
insert KAcvode gkbar_KAcvode = 0.08
insert KDcvode gkbar_KDcvode = 0.09
insert KC3cvode gkbar_KC3cvode = 0.06
insert K23cvode
insert cad taur_cad = 2 cainf_cad = 4e-5
kt_cad = 4e-5 kd_cad = 4e-5
}
SpinyDistalDendrite {cm=1.5
insert Leak gl_Leak = 0.00021 el_Leak = -80
insert CaP2cvode cai = 0.0040600975 cao = 2.4 gcabar_CaP2cvode = 0.004
insert CaTcvode cai = 0.0040600975 cao = 2.4 gcabar_CaTcvode = 0.0015
insert CaEcvode cai = 0.0040600975 cao = 2.4 gcabar_CaEcvode = 0.008
insert Khhcvode gkbar_Khhcvode = 0.0006
insert KMcvode gkbar_KMcvode = 0.00001
insert KAcvode gkbar_KAcvode = 0.08
insert KDcvode gkbar_KDcvode = 0.09
insert KC3cvode gkbar_KC3cvode = 0.06
insert K23cvode
insert cad taur_cad = 2 cainf_cad = 4e-5
kt_cad = 4e-5 kd_cad = 4e-5
}
AdjacentDendrite {cm=1.5
insert Leak gl_Leak = 0.00021 el_Leak = -80
insert CaP2cvode cai = 0.0040600975 cao = 2.4 gcabar_CaP2cvode = 0.004
insert CaTcvode cai = 0.0040600975 cao = 2.4 gcabar_CaTcvode = 0.0015
insert CaEcvode cai = 0.0040600975 cao = 2.4 gcabar_CaEcvode = 0.008
insert Khhcvode gkbar_Khhcvode = 0.0006
insert KMcvode gkbar_KMcvode = 0.00001
insert KAcvode gkbar_KAcvode = 0.08
insert KDcvode gkbar_KDcvode = 0.09
insert KC3cvode gkbar_KC3cvode = 0.06
insert K23cvode
insert cad taur_cad = 2 cainf_cad = 4e-5
kt_cad = 4e-5 kd_cad = 4e-5
}
Spine {cm=0.8
insert Leak gl_Leak = 0.00021 el_Leak = -80
insert CaP2cvode cai = 0.0040600975 cao = 2.4 gcabar_CaP2cvode = 0.004
insert CaTcvode cai = 0.0040600975 cao = 2.4 gcabar_CaTcvode = 0.0015
insert CaEcvode cai = 0.0040600975 cao = 2.4 gcabar_CaEcvode = 0.008
insert Khhcvode gkbar_Khhcvode = 0.0006
insert KMcvode gkbar_KMcvode = 0.00001
insert KAcvode gkbar_KAcvode = 0.08
insert KDcvode gkbar_KDcvode = 0.09
insert KC3cvode gkbar_KC3cvode = 0.06
insert K23cvode
insert cad taur_cad = 2 cainf_cad = 4e-5
kt_cad = 4e-5 kd_cad = 4e-5
}
IO1sm {cm= 0.84197833
insert Leak gl_Leak = 0.000221019 el_Leak = -80
insert CaP2cvode cai = 0.0040600975 cao = 2.4 gcabar_CaP2cvode = 0.004209892
insert CaTcvode cai = 0.0040600975 cao = 2.4 gcabar_CaTcvode = 0.001578709
insert CaEcvode cai = 0.0040600975 cao = 2.4 gcabar_CaEcvode = 0.008419783
insert Khhcvode gkbar_Khhcvode = 0.000631484
insert KMcvode gkbar_KMcvode = 1.05247E-05
insert KAcvode gkbar_KAcvode = 0.084197833
insert KDcvode gkbar_KDcvode = 0.094722562
insert KC3cvode gkbar_KC3cvode = 0.063148375
insert K23cvode gkbar_K23cvode = 0.000410464
insert cad taur_cad = 2 cainf_cad = 4e-5 cainit_cad = 0.0040600975
kt_cad = 4e-5 kd_cad = 4e-5 icaadjust_cad= 0.950143215
}
IO1sp {cm= 7.924229275
insert Leak gl_Leak = 0.001109392 el_Leak = -80
insert CaP2cvode cai = 0.0040600975 cao = 2.4 gcabar_CaP2cvode = 0.021131278
insert CaTcvode cai = 0.0040600975 cao = 2.4 gcabar_CaTcvode = 0.007924229
insert CaEcvode cai = 0.0040600975 cao = 2.4 gcabar_CaEcvode = 0.042262556
insert Khhcvode gkbar_Khhcvode = 0.003169692
insert KMcvode gkbar_KMcvode = 5.28282E-05
insert KAcvode gkbar_KAcvode = 0.422625561
insert KDcvode gkbar_KDcvode = 0.475453756
insert KC3cvode gkbar_KC3cvode = 0.316969171
insert K23cvode gkbar_K23cvode = 0.0020603
insert cad taur_cad = 2 cainf_cad = 4e-5 cainit_cad = 0.0040600975
kt_cad = 4e-5 kd_cad = 4e-5 icaadjust_cad= 0.189292857
}
IO2sm {cm= 0.731103206
insert Leak gl_Leak = 0.000191915 el_Leak = -80
insert CaP2cvode cai = 0.0040600975 cao = 2.4 gcabar_CaP2cvode = 0.003655516
insert CaTcvode cai = 0.0040600975 cao = 2.4 gcabar_CaTcvode = 0.001370819
insert CaEcvode cai = 0.0040600975 cao = 2.4 gcabar_CaEcvode = 0.007311032
insert Khhcvode gkbar_Khhcvode = 0.000548327
insert KMcvode gkbar_KMcvode = 9.13879E-06
insert KAcvode gkbar_KAcvode = 0.073110321
insert KDcvode gkbar_KDcvode = 0.082249111
insert KC3cvode gkbar_KC3cvode = 0.05483274
insert K23cvode gkbar_K23cvode = 0.000356413
insert cad taur_cad = 2 cainf_cad = 4e-5 cainit_cad = 0.0040600975
kt_cad = 4e-5 kd_cad = 4e-5 icaadjust_cad= 1.094236756
}
IO2sp {cm= 6.509497645
insert Leak gl_Leak = 0.00091133 el_Leak = -80
insert CaP2cvode cai = 0.0040600975 cao = 2.4 gcabar_CaP2cvode = 0.01735866
insert CaTcvode cai = 0.0040600975 cao = 2.4 gcabar_CaTcvode = 0.006509498
insert CaEcvode cai = 0.0040600975 cao = 2.4 gcabar_CaEcvode = 0.034717321
insert Khhcvode gkbar_Khhcvode = 0.002603799
insert KMcvode gkbar_KMcvode = 4.33967E-05
insert KAcvode gkbar_KAcvode = 0.347173208
insert KDcvode gkbar_KDcvode = 0.390569859
insert KC3cvode gkbar_KC3cvode = 0.260379906
insert K23cvode gkbar_K23cvode = 0.001692469
insert cad taur_cad = 2 cainf_cad = 4e-5 cainit_cad = 0.0040600975
kt_cad = 4e-5 kd_cad = 4e-5 icaadjust_cad= 0.230432528
}
IO3sm {cm= 0.515767253
insert Leak gl_Leak = 0.000135389 el_Leak = -80
insert CaP2cvode cai = 0.0040600975 cao = 2.4 gcabar_CaP2cvode = 0.002578836
insert CaTcvode cai = 0.0040600975 cao = 2.4 gcabar_CaTcvode = 0.000967064
insert CaEcvode cai = 0.0040600975 cao = 2.4 gcabar_CaEcvode = 0.005157673
insert Khhcvode gkbar_Khhcvode = 0.000386825
insert KMcvode gkbar_KMcvode = 6.44709E-06
insert KAcvode gkbar_KAcvode = 0.051576725
insert KDcvode gkbar_KDcvode = 0.058023816
insert KC3cvode gkbar_KC3cvode = 0.038682544
insert K23cvode gkbar_K23cvode = 0.000251437
insert cad taur_cad = 2 cainf_cad = 4e-5 cainit_cad = 0.0040600975
kt_cad = 4e-5 kd_cad = 4e-5 icaadjust_cad= 1.55108723
}
IO3sp {cm= 3.16498077
insert Leak gl_Leak = 0.000443097 el_Leak = -80
insert CaP2cvode cai = 0.0040600975 cao = 2.4 gcabar_CaP2cvode = 0.008439949
insert CaTcvode cai = 0.0040600975 cao = 2.4 gcabar_CaTcvode = 0.003164981
insert CaEcvode cai = 0.0040600975 cao = 2.4 gcabar_CaEcvode = 0.016879897
insert Khhcvode gkbar_Khhcvode = 0.001265992
insert KMcvode gkbar_KMcvode = 2.10999E-05
insert KAcvode gkbar_KAcvode = 0.168798974
insert KDcvode gkbar_KDcvode = 0.189898846
insert KC3cvode gkbar_KC3cvode = 0.126599231
insert K23cvode gkbar_K23cvode = 0.000822895
insert cad taur_cad = 2 cainf_cad = 4e-5 cainit_cad = 0.0040600975
kt_cad = 4e-5 kd_cad = 4e-5 icaadjust_cad= 0.473936529
}
IO4sm {cm= 0.675232434
insert Leak gl_Leak = 0.000177249 el_Leak = -80
insert CaP2cvode cai = 0.0040600975 cao = 2.4 gcabar_CaP2cvode = 0.003376162
insert CaTcvode cai = 0.0040600975 cao = 2.4 gcabar_CaTcvode = 0.001266061
insert CaEcvode cai = 0.0040600975 cao = 2.4 gcabar_CaEcvode = 0.006752324
insert Khhcvode gkbar_Khhcvode = 0.000506424
insert KMcvode gkbar_KMcvode = 8.44041E-06
insert KAcvode gkbar_KAcvode = 0.067523243
insert KDcvode gkbar_KDcvode = 0.075963649
insert KC3cvode gkbar_KC3cvode = 0.050642433
insert K23cvode gkbar_K23cvode = 0.000329176
insert cad taur_cad = 2 cainf_cad = 4e-5 cainit_cad = 0.0040600975
kt_cad = 4e-5 kd_cad = 4e-5 icaadjust_cad= 1.184777212
}
IO4sp {cm= 4.232977441
insert Leak gl_Leak = 0.000592617 el_Leak = -80
insert CaP2cvode cai = 0.0040600975 cao = 2.4 gcabar_CaP2cvode = 0.01128794
insert CaTcvode cai = 0.0040600975 cao = 2.4 gcabar_CaTcvode = 0.004232977
insert CaEcvode cai = 0.0040600975 cao = 2.4 gcabar_CaEcvode = 0.02257588
insert Khhcvode gkbar_Khhcvode = 0.001693191
insert KMcvode gkbar_KMcvode = 2.82198E-05
insert KAcvode gkbar_KAcvode = 0.225758797
insert KDcvode gkbar_KDcvode = 0.253978646
insert KC3cvode gkbar_KC3cvode = 0.169319098
insert K23cvode gkbar_K23cvode = 0.001100574
insert cad taur_cad = 2 cainf_cad = 4e-5 cainit_cad = 0.0040600975
kt_cad = 4e-5 kd_cad = 4e-5 icaadjust_cad= 0.354360499
}
IO5sp {cm= 1.838873147
insert Leak gl_Leak = 0.000257442 el_Leak = -80
insert CaP2cvode cai = 0.0040600975 cao = 2.4 gcabar_CaP2cvode = 0.004903662
insert CaTcvode cai = 0.0040600975 cao = 2.4 gcabar_CaTcvode = 0.001838873
insert CaEcvode cai = 0.0040600975 cao = 2.4 gcabar_CaEcvode = 0.009807323
insert Khhcvode gkbar_Khhcvode = 0.000735549
insert KMcvode gkbar_KMcvode = 1.22592E-05
insert KAcvode gkbar_KAcvode = 0.098073235
insert KDcvode gkbar_KDcvode = 0.110332389
insert KC3cvode gkbar_KC3cvode = 0.073554926
insert K23cvode gkbar_K23cvode = 0.000478107
insert cad taur_cad = 2 cainf_cad = 4e-5 cainit_cad = 0.0040600975
kt_cad = 4e-5 kd_cad = 4e-5 icaadjust_cad= 0.815716953
}
IO6sp {cm= 1.93409695
insert Leak gl_Leak = 0.000270774 el_Leak = -80
insert CaP2cvode cai = 0.0040600975 cao = 2.4 gcabar_CaP2cvode = 0.005157592
insert CaTcvode cai = 0.0040600975 cao = 2.4 gcabar_CaTcvode = 0.001934097
insert CaEcvode cai = 0.0040600975 cao = 2.4 gcabar_CaEcvode = 0.010315184
insert Khhcvode gkbar_Khhcvode = 0.000773639
insert KMcvode gkbar_KMcvode = 1.2894E-05
insert KAcvode gkbar_KAcvode = 0.103151837
insert KDcvode gkbar_KDcvode = 0.116045817
insert KC3cvode gkbar_KC3cvode = 0.077363878
insert K23cvode gkbar_K23cvode = 0.000502865
insert cad taur_cad = 2 cainf_cad = 4e-5 cainit_cad = 0.0040600975
kt_cad = 4e-5 kd_cad = 4e-5 icaadjust_cad= 0.775555744
}
global_ra = 250
set_ra()
celsius = 37
nrnmainmenu() /* shows the main panel window on the screen */