purkinje_reduced_PPR model_cvode.hoc

/*
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 */