Activity patterns in a subthalamopallidal network of the basal ganglia model (Terman et al 2002)

 Download zip file 
Help downloading and running models
Accession:182758
"Based on recent experimental data, we have developed a conductance-based computational network model of the subthalamic nucleus and the external segment of the globus pallidus in the indirect pathway of the basal ganglia. Computer simulations and analysis of this model illuminate the roles of the coupling architecture of the network, and associated synaptic conductances, in modulating the activity patterns displayed by this network. Depending on the relationships of these coupling parameters, the network can support three general classes of sustained firing patterns: clustering, propagating waves, and repetitive spiking that may show little regularity or correlation. ...". Terman's XPP code and a partial implementation by Taylor Malone in NEURON and python are included.
Reference:
1 . Terman D, Rubin JE, Yew AC, Wilson CJ (2002) Activity patterns in a model for the subthalamopallidal network of the basal ganglia. J Neurosci 22:2963-76 [PubMed]
Citations  Citation Browser
Model Information (Click on a link to find other models with that property)
Model Type: Realistic Network;
Brain Region(s)/Organism: Basal ganglia;
Cell Type(s): Subthalamic nucleus principal GABA cell; Globus pallidus principal GABA cell; Subthalamus nucleus projection neuron; Globus pallidus neuron;
Channel(s): I K; I Na,t; I T low threshold; I Calcium; I_AHP;
Gap Junctions:
Receptor(s):
Gene(s):
Transmitter(s):
Simulation Environment: XPP; NEURON; Python;
Model Concept(s): Activity Patterns; Spatio-temporal Activity Patterns; Parkinson's;
Implementer(s): Terman, David [terman at math.ohio-state.edu]; Malone, Taylor ;
Search NeuronDB for information about:  Globus pallidus principal GABA cell; Subthalamic nucleus principal GABA cell; I Na,t; I T low threshold; I K; I Calcium; I_AHP;
## Set file for off_gtos.ode on Thu Nov 30 16:03:49 2000
96   Number of equations and auxiliaries
126   Number of parameters
# Numerical stuff
1    nout
40    nullcline mesh
8   Qual RK
999  total
0.2  DeltaT
0  T0
0  Transient
5000  Bound
1e-12  DtMin
1  DtMax
0.001  Tolerance
0.001  Abs. Tolerance
0  Max Delay
100   Eigenvector iterates
0.001  Eigenvector tolerance
0.001  Newton tolerance
0  Poincare plane
1e-05  Boundary value tolerance
1e-05  Boundary value epsilon
20   Boundary value iterates
0   Poincare None
1   Poincare variable
1   Poincare sign
0   Stop on Section
0   Delay flag
998.8000000000001  Current time
998.8000000000001  Last Time
1   MyStart
1   INFLAG
# Delays
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
# Bndry conds
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
# Old ICs
-77  V1
0.19  H1
0.15  N1
0.23  R1
0.06  CA1
0  S1
-77  V2
0.19  H2
0.15  N2
0.23  R2
0.06  CA2
0  S2
-53.2  V3
0.1  H3
0.45  N3
0.6  R3
0.12  CA3
0.44  S3
-53.2  V4
0.1  H4
0.45  N4
0.6  R4
0.12  CA4
0.44  S4
-77  V5
0.19  H5
0.15  N5
0.23  R5
0.06  CA5
0  S5
-77  V6
0.19  H6
0.15  N6
0.23  R6
0.06  CA6
0  S6
-53.2  V7
0.1  H7
0.45  N7
0.6  R7
0.12  CA7
0.44  S7
-53.2  V8
0.1  H8
0.45  N8
0.6  R8
0.12  CA8
0.44  S8
-95  VG1
-95  VG2
-77  VG3
-77  VG4
-95  VG5
-95  VG6
-77  VG7
-77  VG8
0.04  NG1
0.04  NG2
0.78  NG3
0.78  NG4
0.04  NG5
0.04  NG6
0.78  NG7
0.78  NG8
0.95  HG1
0.95  HG2
0.2  HG3
0.2  HG4
0.95  HG5
0.95  HG6
0.2  HG7
0.2  HG8
0.9  RG1
0.9  RG2
0.9  RG3
0.9  RG4
0.9  RG5
0.9  RG6
0.9  RG7
0.9  RG8
0  CAG1
0  CAG2
0.035  CAG3
0.035  CAG4
0  CAG5
0  CAG6
0.035  CAG7
0.035  CAG8
0.09  SG1
0.09  SG2
0.5  SG3
0.5  SG4
0.09  SG5
0.09  SG6
0.5  SG7
0.5  SG8
# Ending  ICs
-60.71791026671416  V1
0.1050817070569929  H1
0.3339047915553423  N1
0.5597894302066979  R1
0.2036376651616488  CA1
0.007004151871627591  S1
-60.71791026671416  V2
0.1050817070569929  H2
0.3339047915553423  N2
0.5597894302066979  R2
0.2036376651616488  CA2
0.007004151871627591  S2
-77.47069786413695  V3
0.2873203571604054  H3
0.04663827771005361  N3
0.5103140380255553  R3
0.3004438757198344  CA3
0.0009344040466171509  S3
-77.47069786413695  V4
0.2873203571604054  H4
0.04663827771005361  N4
0.5103140380255553  R4
0.3004438757198344  CA4
0.0009344040466171509  S4
-60.71791026671416  V5
0.1050817070569929  H5
0.3339047915553423  N5
0.5597894302066979  R5
0.2036376651616488  CA5
0.007004151871627591  S5
-60.71791026671416  V6
0.1050817070569929  H6
0.3339047915553423  N6
0.5597894302066979  R6
0.2036376651616488  CA6
0.007004151871627591  S6
-77.47069786413695  V7
0.2873203571604054  H7
0.04663827771005361  N7
0.5103140380255553  R7
0.3004438757198344  CA7
0.0009344040466171509  S7
-77.47069786413695  V8
0.2873203571604054  H8
0.04663827771005361  N8
0.5103140380255553  R8
0.3004438757198344  CA8
0.0009344040466171509  S8
-75.10391218638624  VG1
-75.10391218638624  VG2
-72.92625188400277  VG3
-72.92625188400277  VG4
-75.10391218638624  VG5
-75.10391218638624  VG6
-72.92625188400277  VG7
-72.92625188400277  VG8
0.2125076036533668  NG1
0.2125076036533668  NG2
0.1598693544220592  NG3
0.1598693544220592  NG4
0.2125076036533668  NG5
0.2125076036533668  NG6
0.1598693544220592  NG7
0.1598693544220592  NG8
0.7466959712543473  HG1
0.7466959712543473  HG2
0.7805846551106687  HG3
0.7805846551106687  HG4
0.7466959712543473  HG5
0.7466959712543473  HG6
0.7805846551106687  HG7
0.7805846551106687  HG8
0.7677417604636313  RG1
0.7677417604636313  RG2
0.7813970669429506  RG3
0.7813970669429506  RG4
0.7677417604636313  RG5
0.7677417604636313  RG6
0.7813970669429506  RG7
0.7813970669429506  RG8
0.05416816967535378  CAG1
0.05416816967535378  CAG2
0.06618847501824002  CAG3
0.06618847501824002  CAG4
0.05416816967535378  CAG5
0.05416816967535378  CAG6
0.06618847501824002  CAG7
0.06618847501824002  CAG8
0.5740762310666904  SG1
0.5740762310666904  SG2
0.005156071270063888  SG3
0.005156071270063888  SG4
0.5740762310666904  SG5
0.5740762310666904  SG6
0.005156071270063888  SG7
0.005156071270063888  SG8
# Parameters
2.25  gl
37.5  gna
45  gk
9  gahp
0.5  gca
0.5  gt
-60  vl
55  vna
-80  vk
140  vca
-30  thetam
-15  sigmam
-39  thetah
3.1  sigmah
-32  thetan
-8  sigman
-67  thetar
2  sigmar
-63  thetaa
-7.8  sigmaa
0.4  thetab
-0.1  sigmab
-39  thetas
-8  sigmas
1  tauh0
500  tauh1
-57  thh
3  sigmaht
0.75  phi
1  taun0
100  taun1
-80  thn
26  sigmant
40  taur0
17.5  taur1
68  thr
2.2  sigmart
0.2  phir
15  k1
5e-05  eps
22.5  kca
0  i1
0  i2
0  i3
0  i4
0  i5
0  i6
0  i7
0  i8
0  i9
0  i10
0  i11
0  i12
0  i13
0  i14
0  i15
0  i16
5  alpha
1  beta
30  thetag
4.5  gGtoS
-85  vGtoS
-39  thetH
-8  sigmH
120  gnag
30  gkg
30  gahpg
0.5  gtg
0.15  gcag
0.1  glg
55  vnag
-80  vkg
120  vcag
-55  vlg
-37  thmg
-10  sigmg
-58  thhg
12  sighg
-50  thng
-14  signg
-70  thrg
2  sigrg
30  taurg
-57  thag
-2  sigag
-35  thsg
-2  sigsg
0.05  tauhg0
0.27  tauhg1
-40  thhgt
12  shg
0.05  taung0
0.27  taung1
-40  thngt
12  sng
30  k1g
20  kcag
0.0001  epsg
1  phirg
0.05  phing
0.05  phihg
-0.6  ig1
-0.6  ig2
-0.6  ig3
-0.6  ig4
-0.6  ig5
-0.6  ig6
-0.6  ig7
-0.6  ig8
-1  ig9
-1  ig10
-1  ig11
-1  ig12
-1  ig13
-1  ig14
-1  ig15
-1  ig16
0.2  gStoG
0  vStoG
2  alphag
0.04  betag
20  thetagg
-85  vGtoG
0.1  gGtoG
-57  thetHg
-2  sigmHg
# Graphics
0.7071097251469876  rm
-0.4999999999826661  rm
0.4999958366025517  rm
0.7071038372138506  rm
0.5000041633974484  rm
-0.4999999999826661  rm
0  rm
0.7071038372138506  rm
0.7071097251469876  rm
0    
51    
2    
1    
0    
0    
13    
2    
1    
3    
0    
1    
2    
1    
0    
0    
1    
2    
1    
0    
0    
1    
2    
1    
0    
0    
1    
2    
1    
0    
0    
1    
2    
1    
0    
0    
1    
2    
1    
0    
0    
1    
2    
1    
0    
0    
1    
2    
1    
0    
-1000   
1000   
0    
0   3DFlag
1   Timeflag
0   Colorflag
0   Type
1  color scale
0   minscale
998.7999877929688   xmax
0   xmin
54.69213485717773   ymax
-79.85974884033203   ymin
12   zmax
-12   zmin
499.3999938964844   
0.002002402907932915   
-12.58380699157715   
0.01486415459256045   
0   
0.08333333333333333   
45   Theta
45   Phi
0    xshft
0    yshft
0    zshft
0   xlo
-79.85974884033203   ylo
0   
-80   
998.7999877929688   xhi
54.69213485717773   yhi
1000   
20   
# Transpose variables etc
V1
2   n columns
1   n rows
1   row skip
1   col skip
1   row 0
# Coupling stuff for H funs
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
# Array plot stuff
vg1
8   NCols
0   Row 1
5000   NRows
1   RowSkip
-90  Zmin
20  Zmax
# Torus information 
0    Torus flag 1=ON
6.283185308  Torus period
# Range information
gl
-1   eq-range stab col
0   shoot flag 1=on
10   eq-range steps
0  eq_range low
1  eq_range high
V1
V1
20   Range steps
0   Cycle color 1=on
1   Reset data 1=on
1   Use old I.C.s 1=yes
0  Par1 low
0  Par2 low
1  Par1 high
1  Par2 high
gl
0   BVP side
0   color cycle flag 1=on
10   BVP range steps
0  BVP range low
1  BVP range high
RHS etc ...
dV1/dT=-(CURR(V1,H1,N1,CA1,R1))-(ISYN(V1,SG3)+ISYN(V1,SG7))+I1  
dH1/dT=PHI*( HINF(V1)-H1 )/TAUH(V1)  
dN1/dT=PHI*( NINF(V1)-N1 )/TAUN(V1)  
dR1/dT=PHIR*(RINF(V1)-R1)/TAUR(V1)  
dCA1/dT=PHI*EPS*(-GCA*((SINF(V1))^2)*(V1-VCA)-IT(V1,R1) - KCA*CA1)  
dS1/dT=ALPHA*(1-S1)*HIN(V1-THETAG)-BETA*S1  
dV2/dT=-(CURR(V2,H2,N2,CA2,R2))-(ISYN(V2,SG4)+ISYN(V2,SG8))+I2  
dH2/dT=PHI*( HINF(V2)-H2)/TAUH(V2)  
dN2/dT=PHI*( NINF(V2)-N2)/TAUN(V2)  
dR2/dT=PHIR*(RINF(V2)-R2)/TAUR(V2)  
dCA2/dT=PHI*EPS*(-GCA*((SINF(V2))^2)*(V2-VCA)-IT(V2,R2) - KCA*CA2)  
dS2/dT=ALPHA*(1-S2)*HIN(V2-THETAG)-BETA*S2  
dV3/dT=-(CURR(V3,H3,N3,CA3,R3))-(ISYN(V3,SG5)+ISYN(V3,SG1))+I3  
dH3/dT=PHI*( HINF(V3)-H3)/TAUH(V3)  
dN3/dT=PHI*( NINF(V3)-N3)/TAUN(V3)  
dR3/dT=PHIR*(RINF(V3)-R3)/TAUR(V3)  
dCA3/dT=PHI*EPS*(-GCA*((SINF(V3))^2)*(V3-VCA)-IT(V3,R3) - KCA*CA3)  
dS3/dT=ALPHA*(1-S3)*HIN(V3-THETAG)-BETA*S3  
dV4/dT=-(CURR(V4,H4,N4,CA4,R4))-(ISYN(V4,SG6)+ISYN(V4,SG2))+I4  
dH4/dT=PHI*( HINF(V4)-H4)/TAUH(V4)  
dN4/dT=PHI*( NINF(V4)-N4)/TAUN(V4)  
dR4/dT=PHIR*(RINF(V4)-R4)/TAUR(V4)  
dCA4/dT=PHI*EPS*(-GCA*((SINF(V4))^2)*(V4-VCA)-IT(V4,R4) - KCA*CA4)  
dS4/dT=ALPHA*(1-S4)*HIN(V4-THETAG)-BETA*S4  
dV5/dT=-(CURR(V5,H5,N5,CA5,R5))-(ISYN(V5,SG7)+ISYN(V5,SG3))+I5  
dH5/dT=PHI*( HINF(V5)-H5)/TAUH(V5)  
dN5/dT=PHI*( NINF(V5)-N5)/TAUN(V5)  
dR5/dT=PHIR*(RINF(V5)-R5)/TAUR(V5)  
dCA5/dT=PHI*EPS*(-GCA*((SINF(V5))^2)*(V5-VCA)-IT(V5,R5) - KCA*CA5)  
dS5/dT=ALPHA*(1-S5)*HIN(V5-THETAG)-BETA*S5  
dV6/dT=-(CURR(V6,H6,N6,CA6,R6))-(ISYN(V6,SG8)+ISYN(V6,SG4))+I6  
dH6/dT=PHI*( HINF(V6)-H6)/TAUH(V6)  
dN6/dT=PHI*( NINF(V6)-N6)/TAUN(V6)  
dR6/dT=PHIR*(RINF(V6)-R6)/TAUR(V6)  
dCA6/dT=PHI*EPS*(-GCA*((SINF(V6))^2)*(V6-VCA)-IT(V6,R6) - KCA*CA6)  
dS6/dT=ALPHA*(1-S6)*HIN(V6-THETAG)-BETA*S6  
dV7/dT=-(CURR(V7,H7,N7,CA7,R7))-(ISYN(V7,SG1)+ISYN(V7,SG5))+I7  
dH7/dT=PHI*( HINF(V7)-H7)/TAUH(V7)  
dN7/dT=PHI*( NINF(V7)-N7)/TAUN(V7)  
dR7/dT=PHIR*(RINF(V7)-R7)/TAUR(V7)  
dCA7/dT=PHI*EPS*(-GCA*((SINF(V7))^2)*(V7-VCA)-IT(V7,R7) - KCA*CA7)  
dS7/dT=ALPHA*(1-S7)*HIN(V7-THETAG)-BETA*S7  
dV8/dT=-(CURR(V8,H8,N8,CA8,R8))-(ISYN(V8,SG2)+ISYN(V8,SG6))+I8  
dH8/dT=PHI*( HINF(V8)-H8)/TAUH(V8)  
dN8/dT=PHI*( NINF(V8)-N8)/TAUN(V8)  
dR8/dT=PHIR*(RINF(V8)-R8)/TAUR(V8)  
dCA8/dT=PHI*EPS*(-GCA*((SINF(V8))^2)*(V8-VCA)-IT(V8,R8) - KCA*CA8)  
dS8/dT=ALPHA*(1-S8)*HIN(V8-THETAG)-BETA*S8  
dVG1/dT=-CURRG(VG1,HG1,NG1,CAG1,RG1)-ISYNG(VG1,S1)+IG1-ISYNGG(VG1,SG2+SG8)  
dVG2/dT=-CURRG(VG2,HG2,NG2,CAG2,RG2)-ISYNG(VG2,S2)+IG2-ISYNGG(VG2,SG1+SG3)  
dVG3/dT=-CURRG(VG3,HG3,NG3,CAG3,RG3)-ISYNG(VG3,S3)+IG3-ISYNGG(VG3,SG2+SG4)  
dVG4/dT=-CURRG(VG4,HG4,NG4,CAG4,RG4)-ISYNG(VG4,S4)+IG4-ISYNGG(VG4,SG3+SG5)  
dVG5/dT=-CURRG(VG5,HG5,NG5,CAG5,RG5)-ISYNG(VG5,S5)+IG5-ISYNGG(VG5,SG4+SG6)  
dVG6/dT=-CURRG(VG6,HG6,NG6,CAG6,RG6)-ISYNG(VG6,S6)+IG6-ISYNGG(VG6,SG5+SG7)  
dVG7/dT=-CURRG(VG7,HG7,NG7,CAG7,RG7)-ISYNG(VG7,S7)+IG7-ISYNGG(VG7,SG6+SG8)  
dVG8/dT=-CURRG(VG8,HG8,NG8,CAG8,RG8)-ISYNG(VG8,S8)+IG8-ISYNGG(VG8,SG1+SG7)  
dNG1/dT= PHING*(NINFG(VG1)-NG1)/TAUNG(VG1)  
dNG2/dT= PHING*(NINFG(VG2)-NG2)/TAUNG(VG2)  
dNG3/dT= PHING*(NINFG(VG3)-NG3)/TAUNG(VG3)  
dNG4/dT= PHING*(NINFG(VG4)-NG4)/TAUNG(VG4)  
dNG5/dT= PHING*(NINFG(VG5)-NG5)/TAUNG(VG5)  
dNG6/dT= PHING*(NINFG(VG6)-NG6)/TAUNG(VG6)  
dNG7/dT= PHING*(NINFG(VG7)-NG7)/TAUNG(VG7)  
dNG8/dT= PHING*(NINFG(VG8)-NG8)/TAUNG(VG8)  
dHG1/dT= PHIHG*(HINFG(VG1)-HG1)/TAUHG(VG1)  
dHG2/dT= PHIHG*(HINFG(VG2)-HG2)/TAUHG(VG2)  
dHG3/dT= PHIHG*(HINFG(VG3)-HG3)/TAUHG(VG3)  
dHG4/dT= PHIHG*(HINFG(VG4)-HG4)/TAUHG(VG4)  
dHG5/dT= PHIHG*(HINFG(VG5)-HG5)/TAUHG(VG5)  
dHG6/dT= PHIHG*(HINFG(VG6)-HG6)/TAUHG(VG6)  
dHG7/dT= PHIHG*(HINFG(VG7)-HG7)/TAUHG(VG7)  
dHG8/dT= PHIHG*(HINFG(VG8)-HG8)/TAUHG(VG8)  
dRG1/dT=PHIRG*(RINFG(VG1)-RG1)/TAURG  
dRG2/dT=PHIRG*(RINFG(VG2)-RG2)/TAURG  
dRG3/dT=PHIRG*(RINFG(VG3)-RG3)/TAURG  
dRG4/dT=PHIRG*(RINFG(VG4)-RG4)/TAURG  
dRG5/dT=PHIRG*(RINFG(VG5)-RG5)/TAURG  
dRG6/dT=PHIRG*(RINFG(VG6)-RG6)/TAURG  
dRG7/dT=PHIRG*(RINFG(VG7)-RG7)/TAURG  
dRG8/dT=PHIRG*(RINFG(VG8)-RG8)/TAURG  
dCAG1/dT=EPSG*(-ICAG(VG1)-ITG(VG1,RG1) - KCAG*CAG1)  
dCAG2/dT=EPSG*(-ICAG(VG2)-ITG(VG2,RG2) - KCAG*CAG2)  
dCAG3/dT=EPSG*(-ICAG(VG3)-ITG(VG3,RG3) - KCAG*CAG3)  
dCAG4/dT=EPSG*(-ICAG(VG4)-ITG(VG4,RG4) - KCAG*CAG4)  
dCAG5/dT=EPSG*(-ICAG(VG5)-ITG(VG5,RG5) - KCAG*CAG5)  
dCAG6/dT=EPSG*(-ICAG(VG6)-ITG(VG6,RG6) - KCAG*CAG6)  
dCAG7/dT=EPSG*(-ICAG(VG7)-ITG(VG7,RG7) - KCAG*CAG7)  
dCAG8/dT=EPSG*(-ICAG(VG8)-ITG(VG8,RG8) - KCAG*CAG8)  
dSG1/dT=ALPHAG*(1-SG1)*HING(VG1-THETAGG)-BETAG*SG1  
dSG2/dT=ALPHAG*(1-SG2)*HING(VG2-THETAGG)-BETAG*SG2  
dSG3/dT=ALPHAG*(1-SG3)*HING(VG3-THETAGG)-BETAG*SG3  
dSG4/dT=ALPHAG*(1-SG4)*HING(VG4-THETAGG)-BETAG*SG4  
dSG5/dT=ALPHAG*(1-SG5)*HING(VG5-THETAGG)-BETAG*SG5  
dSG6/dT=ALPHAG*(1-SG6)*HING(VG6-THETAGG)-BETAG*SG6  
dSG7/dT=ALPHAG*(1-SG7)*HING(VG7-THETAGG)-BETAG*SG7  
dSG8/dT=ALPHAG*(1-SG8)*HING(VG8-THETAGG)-BETAG*SG8  

User-defined functions:
MINF(V) = 1./(1.+EXP((V-THETAM)/SIGMAM))  
HINF(V) = 1./(1.+EXP((V-THETAH)/SIGMAH))  
NINF(V) = 1./(1.+EXP((V-THETAN)/SIGMAN))  
RINF(V) = 1/(1+EXP((V-THETAR)/SIGMAR))  
AINF(V) = 1/(1+EXP((V-THETAA)/SIGMAA))  
BINF(R) = 1/(1+EXP((R-THETAB)/SIGMAB))-1/(1+EXP(-THETAB/SIGMAB))  
SINF(V) = 1./(1.+EXP((V-THETAS)/SIGMAS))  
TAUH(V) = TAUH0+TAUH1/(1+EXP((V-THH)/SIGMAHT))  
TAUN(V) = TAUN0+TAUN1/(1+EXP((V-THN)/SIGMANT))  
TAUR(V) = TAUR0+TAUR1/(1+EXP((V-THR)/SIGMART))  
HIN(V) = 1/(1+EXP((V-THETH)/SIGMH))  
MINFG(V) = 1./(1.+EXP((V-THMG)/SIGMG))  
HINFG(V) = 1/(1+EXP((V-THHG)/SIGHG))  
NINFG(V) = 1./(1.+EXP((V-THNG)/SIGNG))  
RINFG(V) = 1/(1+EXP((V-THRG)/SIGRG))  
AINFG(V) = 1/(1+EXP((V-THAG)/SIGAG))  
SINFG(V) = 1/(1+EXP((V-THSG)/SIGSG))  
TAUHG(V) = TAUHG0+TAUHG1/(1+EXP((V-THHGT)/SHG))  
TAUNG(V) = TAUNG0+TAUNG1/(1+EXP((V-THNGT)/SNG))  
HING(V) = 1/(1+EXP((V-THETHG)/SIGMHG))  
IL(V) = GL*(V-VL)  
INA(V,H) = GNA*(MINF(V))^3*H*(V-VNA)  
IK(V,N) = GK*N^4*(V-VK)  
IAHP(V,CA) = GAHP*(V-VK)*CA/(CA+K1)  
ICA(V) = GCA*((SINF(V))^2)*(V-VCA)  
IT(V,R) = GT*(AINF(V)**3)*(BINF(R)^2)*(V-VCA)  
ISYN(V,S) = GGTOS*S*(V-VGTOS)  
CURR(V,H,N,CA,R) = IL(V)+INA(V,H)+IK(V,N)+IAHP(V,CA)+ICA(V)+IT(V,R)  
ILG(V) = GLG*(V-VLG)  
INAG(V,H) = GNAG*(MINFG(V)^3)*H*(V-VNAG)  
IKG(V,N) = GKG*(N^4)*(V-VKG)  
IAHPG(V,CA) = GAHPG*(V-VKG)*CA/(CA+K1G)  
ICAG(V) = GCAG*((SINFG(V))^2)*(V-VCAG)  
ITG(V,R) = GTG*(AINFG(V)^3)*R*(V-VCAG)  
ISYNG(VG,S) = GSTOG*S*(VG-VSTOG)  
ISYNGG(VG,S) = GGTOG*S*(VG-VGTOG)  
CURRG(V,H,N,CA,R) = ITG(V,R)+INAG(V,H)+IKG(V,N)+ILG(V)+IAHPG(V,CA)+ICAG(V)