tonic.mod 772 B

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748
  1. : tonic current with rectification
  2. : based on Pavlov et al (J Neuro 2009)
  3. NEURON{
  4. SUFFIX tonic
  5. NONSPECIFIC_CURRENT i
  6. RANGE i, v, a, b, g, e_gaba}
  7. PARAMETER{
  8. g = 0.001 (siemens/cm2)
  9. e_gaba = -80 (millivolt)
  10. }
  11. ASSIGNED{
  12. v (millivolt)
  13. i (milliampere)
  14. a (/ms)
  15. b (/ms)}
  16. STATE{o c}
  17. BREAKPOINT{
  18. SOLVE kin METHOD sparse
  19. i = g*o*(v-e_gaba)}
  20. INITIAL {SOLVE kin STEADYSTATE sparse}
  21. KINETIC kin{
  22. rates(v)
  23. ~ o<->c (b, a) : b is forward rate constant, a backward
  24. CONSERVE o+c=1}
  25. PROCEDURE rates(v(millivolt)) {
  26. LOCAL x, y
  27. UNITSOFF
  28. x = 0.1*(v+20)
  29. if (fabs(x)>1e-6){
  30. a = (50*x)/(1-exp(-x))
  31. } else{
  32. a=0.25*(v^2+(20*v)+200)
  33. }
  34. y = -0.08*(v-10)
  35. if(fabs(x)>1e-6){
  36. b = (20*y)/(1-exp(-y))
  37. } else{
  38. b = -0.064*(v^2-(45*v)+37.5)
  39. }
  40. UNITSON}