MicrocircuitDB: Direct recruitment of S1 pyramidal cells and interneurons via ICMS (Overstreet et al., 2013)

/*
CK Overstreet and JD Klein
Last Updated 12/20/2012
Major characteristics of cell bodies and axons reconstructed by hand from:
H.D. Schwark, “The distribution of intrinsic cortical axons in area 3b of cat primary somatosensory cortex”, Experimental Brain Research, vol. 78, no. 3, pp. 501-513, Aug. 1989.
*/

load_file("nrngui.hoc")

NAPICAL=42
NBOUTON=19

create soma, ap[NAPICAL], bt[NBOUTON]

access soma

connect ap[0](0), soma(0)
connect ap[1](0), ap[0](1)
connect ap[2](0), ap[0](1)
connect ap[3](0), ap[2](1)
connect ap[4](0), ap[3](1)
connect ap[5](0), ap[4](1)
connect ap[6](0), ap[4](1)
connect ap[7](0), ap[2](1)
connect ap[8](0), ap[2](1)
connect ap[9](0), ap[2](1)
connect ap[10](0), ap[9](1)
connect ap[11](0), ap[9](1)
connect ap[12](0), ap[11](1)
connect ap[13](0), ap[11](1)
connect ap[14](0), ap[13](1)
connect ap[15](0), ap[13](1)
connect ap[16](0), ap[9](1)
connect ap[17](0), ap[16](1)
connect ap[18](0), ap[16](1)
connect ap[19](0), ap[18](1)
connect ap[20](0), ap[19](1)
connect ap[21](0), ap[19](1)
connect ap[22](0), ap[18](1)
connect ap[23](0), ap[22](1)
connect ap[24](0), ap[23](1)
connect ap[25](0), ap[23](1)
connect ap[26](0), ap[22](1)
connect ap[27](0), ap[26](1)
connect ap[28](0), ap[26](1)
connect ap[29](0), ap[28](1)
connect ap[30](0), ap[29](1)
connect ap[31](0), ap[30](1)
connect ap[32](0), ap[29](1)
connect ap[33](0), ap[32](1)
connect ap[34](0), ap[28](1)
connect ap[35](0), ap[34](1)
connect ap[36](0), ap[28](1)
connect ap[37](0), ap[36](1)
connect ap[38](0), ap[37](1)
connect ap[39](0), ap[21](1)
connect ap[40](0), ap[39](1)
connect ap[41](0), ap[40](1)

connect bt[0](0), ap[1](1)
connect bt[1](0), ap[5](1)
connect bt[2](0), ap[6](1)
connect bt[3](0), ap[7](1)
connect bt[4](0), ap[8](1)
connect bt[5](0), ap[10](1)
connect bt[6](0), ap[12](1)
connect bt[7](0), ap[14](1)
connect bt[8](0), ap[15](1)
connect bt[9](0), ap[17](1)
connect bt[10](0), ap[20](1)
//connect bt[11](0), ap[21](1)
connect bt[11](0), ap[41](1)
connect bt[12](0), ap[24](1)
connect bt[13](0), ap[25](1)
connect bt[14](0), ap[27](1)
connect bt[15](0), ap[31](1)
connect bt[16](0), ap[33](1)
connect bt[17](0), ap[35](1)
connect bt[18](0), ap[38](1)


forall pt3dclear()

soma {
	pt3dadd(0, 0, 0, 10)
	pt3dadd(0, 10, 0, 10)
	insert axnode
	insert extracellular
	insert xtra
	xraxial=1e+09
	xg=1e+09
	xc=0
	e_extracellular=0

}

for i = 0, NAPICAL-1 ap[i] {
	insert axnode
	insert extracellular
	insert xtra
	xraxial=1e+09
	xg=1e+09
	xc=0
	e_extracellular=0
}

for i = 0, NBOUTON-1 bt[i] {
	insert axnode
	insert extracellular
	insert xtra
	xraxial=1e+09
	xg=1e+09
	xc=0
	e_extracellular=0

}

ap[0] {
	pt3dadd(0, 10, 0, 1)
	pt3dadd(25, -25, 0, 1)
}
ap[1] {
	pt3dadd(25, -25, 0, 0.5)
	pt3dadd(45, 25, 0, 0.25)
}
ap[2] {
	pt3dadd(25, -25, 0, 1)
	pt3dadd(25, -50, 0, 1)
}
ap[3] {
	pt3dadd(25, -50, 0, .5)
	pt3dadd(-50, -25, 0, .5)
}
ap[4] {
	pt3dadd(-50, -25, 0, .25)
	pt3dadd(-75, 25, 0, .25)
}
ap[5] {
	pt3dadd(-75, 25, 0, .25)
	pt3dadd(-125, 25, 0, .25)
}
ap[6] {
	pt3dadd(-50, -25, 0, .25)
	pt3dadd(-25, 75, 0, .25)
}
ap[7] {
	pt3dadd(25, -50, 0, .5)
	pt3dadd(100, -25, 0, .25)
}
ap[8] {
	pt3dadd(25, -50, 0, .5)
	pt3dadd(375, -75, 0, .25)
}
ap[9] {
	pt3dadd(25, -50, 0, 1)
	pt3dadd(25, -150, 0, 1)
}
ap[10] {
	pt3dadd(25, -150, 0, .5)
	pt3dadd(-50, -150, 0, .25)
}
ap[11] {
	pt3dadd(25, -150, 0, .5)
	pt3dadd(700, -100, 0, .5)
}
ap[12] {
	pt3dadd(700, -100, 0, .25)
	pt3dadd(750, -25, 0, .25)
}
ap[13] {
	pt3dadd(700, -100, 0, .5)
	pt3dadd(900, -125, 0, .25)
}
ap[14] {
	pt3dadd(900, -125, 0, .25)
	pt3dadd(975, -100, 0, .25)
}
ap[15] {
	pt3dadd(900, -125, 0, .25)
	pt3dadd(950, -150, 0, .25)
}
ap[16] {
	pt3dadd(25, -150, 0, 1)
	pt3dadd(25, -225, 0, 1)
}
ap[17] {
	pt3dadd(25, -225, 0, .5)
	pt3dadd(150, -275, 0, .25)
}
ap[18] {
	pt3dadd(25, -225, 0, 1)
	pt3dadd(25, -525, 0, 1)
}
ap[19] {
	pt3dadd(25, -525, 0, 1)
	pt3dadd(25, -700, 0, 1)
}
ap[20] {
	pt3dadd(25, -700, 0, .25)
	pt3dadd(100, -700, 0, .25)
}
ap[21] {
	pt3dadd(25, -700, 0, .5)
	pt3dadd(25, -775, 0, .5)
}
ap[22] {
	pt3dadd(25, -525, 0, 1)
	pt3dadd(950, -525, 0, 1)
}
ap[23] {
	pt3dadd(950, -525, 0, .5)
	pt3dadd(1000, -600, 0, .5)
}
ap[24] {
	pt3dadd(1000, -600, 0, .25)
	pt3dadd(975, -650, 0, .25)
}
ap[25] {
	pt3dadd(1000, -600, 0, .25)
	pt3dadd(1200, -700, 0, .25)
}
ap[26] {
	pt3dadd(950, -525, 0, 1)
	pt3dadd(1050, -525, 0, 1)
}
ap[27] {
	pt3dadd(1050, -525, 0, .5)
	pt3dadd(1300, -575, 0, .25)
}
ap[28] {
	pt3dadd(1050, -525, 0, .5)
	pt3dadd(1575, -525, 0, .5)
}
ap[29] {
	pt3dadd(1575, -525, 0, .25)
	pt3dadd(1575, -400, 0, .25)
}
ap[30] {
	pt3dadd(1575, -400, 0, .25)
	pt3dadd(1525, -350, 0, .25)
}
ap[31] {
	pt3dadd(1525, -350, 0, .25)
	pt3dadd(1550, -300, 0, .25)
}
ap[32] {
	pt3dadd(1575, -400, 0, .25)
	pt3dadd(1625, -350, 0, .25)
}
ap[33] {
	pt3dadd(1625, -350, 0, .25)
	pt3dadd(1650, -400, 0, .25)
}
ap[34] {
	pt3dadd(1575, -525, 0, .5)
	pt3dadd(1700, -525, 0, .5)
}
ap[35] {
	pt3dadd(1700, -525, 0, .25)
	pt3dadd(2000, -400, 0, .25)
}
ap[36] {
	pt3dadd(1575, -525, 0, .5)
	pt3dadd(1675, -600, 0, .5)
}
ap[37] {
	pt3dadd(1675, -600, 0, .25)
	pt3dadd(1825, -600, 0, .25)
}
ap[38] {
	pt3dadd(1825, -600, 0, .25)
	pt3dadd(1825, -350, 0, .25)
}
ap[39] {
	pt3dadd(25, -775, 0, .5)
	pt3dadd(75, -950, 0, .5)
}
ap[40] {
	pt3dadd(75, -950, 0, .5)
	pt3dadd(50, -1050, 0, .5)
}
ap[41] {
	pt3dadd(50, -1050, 0, .5)
	pt3dadd(0, -1500, 0, .25)
}


bt[0] {
	pt3dadd(45, 25, 0, 0.25)
	pt3dadd(46.9, 29.6, 0, 3)
}
bt[1] {
	pt3dadd(-125, 25, 0, 0.25)
	pt3dadd(-130.0, 25.0, 0, 3)
}
bt[2] {
	pt3dadd(-25, 75, 0, 0.25)
	pt3dadd(-23.8, 79.9, 0, 3)
}
bt[3] {
	pt3dadd(100, -25, 0, 0.25)
	pt3dadd(104.7, -23.4, 0, 3)
}
bt[4] {
	pt3dadd(375, -75, 0, 0.25)
	pt3dadd(380.0, -75.4, 0, 3)
}
bt[5] {
	pt3dadd(-50, -150, 0, 0.25)
	pt3dadd(-55.0, -150.0, 0, 3)
}
bt[6] {
	pt3dadd(750, -25, 0, 0.25)
	pt3dadd(752.8, -20.8, 0, 3)
}
bt[7] {
	pt3dadd(975, -100, 0, 0.25)
	pt3dadd(979.7, -98.4, 0, 3)
}
bt[8] {
	pt3dadd(950, -150, 0, 0.25)
	pt3dadd(954.5, -152.2, 0, 3)
}
bt[9] {
	pt3dadd(150, -275, 0, 0.25)
	pt3dadd(154.6, -276.9, 0, 3)
}
bt[10] {
	pt3dadd(100, -700, 0, 0.25)
	pt3dadd(105.0, -700.0, 0, 3)
}
bt[11] {
	pt3dadd(0, -1500, 0, 0.25)
	pt3dadd(-0.6, -1505.0, 0, 3)
}
bt[12] {
	pt3dadd(975, -650, 0, 0.25)
	pt3dadd(972.8, -654.5, 0, 3)
}
bt[13] {
	pt3dadd(1200, -700, 0, 0.25)
	pt3dadd(1204.5, -702.2, 0, 3)
}
bt[14] {
	pt3dadd(1300, -575, 0, 0.25)
	pt3dadd(1304.9, -576.0, 0, 3)
}
bt[15] {
	pt3dadd(1550, -300, 0, 0.25)
	pt3dadd(1552.2, -295.5, 0, 3)
}
bt[16] {
	pt3dadd(1650, -400, 0, 0.25)
	pt3dadd(1652.2, -404.5, 0, 3)
}
bt[17] {
	pt3dadd(2000, -400, 0, 0.25)
	pt3dadd(2004.6, -398.1, 0, 3)
}
bt[18] {
	pt3dadd(1825, -350, 0, 0.25)
	pt3dadd(1825.0, -345.0, 0, 3)
}


XORIGIN = 0
YORIGIN = 0
ZORIGIN = 0

access soma
pt3dchange(0, x3d(0)+XORIGIN, y3d(0)+YORIGIN, z3d(0)+ZORIGIN, diam3d(0))
pt3dchange(1, x3d(1)+XORIGIN, y3d(1)+YORIGIN, z3d(1)+ZORIGIN, diam3d(1))
nseg=21

for(i=0; i<=NAPICAL-1; i=i+1){
	access ap[i]
	pt3dchange(0, x3d(0)+XORIGIN, y3d(0)+YORIGIN, z3d(0)+ZORIGIN, diam3d(0))
	pt3dchange(1, x3d(1)+XORIGIN, y3d(1)+YORIGIN, z3d(1)+ZORIGIN, diam3d(1))
	nseg=101
}

for(i=0; i<=NBOUTON-1; i=i+1){
	access bt[i]
	pt3dchange(0, x3d(0)+XORIGIN, y3d(0)+YORIGIN, z3d(0)+ZORIGIN, diam3d(0))
	pt3dchange(1, x3d(1)+XORIGIN, y3d(1)+YORIGIN, z3d(1)+ZORIGIN, diam3d(1))
	nseg=5
}

forall Ra=70
access ap[0]