source: project/release/3/nemo/trunk/examples/AKP06/PotIhCa.scm @ 11996

;;  -*- Hen -*-

(nemo-model PotIhCa
            ((input v  cai cao)
             
             (const Vrest   = -68)
             (const diam    = 20)
             (const celsius = 24)
             
             (const temp_adj = (pow 3 (/ (- celsius 22) 10)))
             
             (defun ghk (v celsius ci co)
               (let ((F 96485.0) (R 8.3145))
                 (let ((zeta (/ (* 2e-3 F v) (* R (+ 273.19 celsius)))))
                   (if (< (abs (- 1.0 (exp (neg zeta)))) 1e-6)
                       (* 1e-6  (* 2 F) (- ci (* co (exp (neg zeta)))) (+ 1.0 (/ zeta 2)))
                       (/ (* 1e-6 (* 2 zeta F) (- ci (* co (exp (neg zeta))))) (- 1.0 (exp (neg zeta))))))))
             
             
             (component (type ion-channel) (name Kv1)
                        
                        (component (type gate)
                                   
                                   ;; rate functions
                                   (defun Kv1_amf (v) 
                                     (let ((cma 0.12889)
                                           (cka -33.90877)
                                           (cva 45))
                                       (* cma (exp (neg (/ (+ v cva) cka))))))
                                   
                                   (defun Kv1_bmf (v)
                                     (let ((cmb   0.12889)
                                           (ckb   12.42101)
                                           (cvb   45))
                                       (* cmb (exp (neg (/ (+ v cvb) ckb))))))
                                   
                                   (hh-ionic-conductance 
                                    (Kv1  ;; ion name: exported variables will be of the form {ion}_{id}
                                     (initial-m (/ (Kv1_amf Vrest) (+ (Kv1_amf Vrest) (Kv1_bmf Vrest))))
                                     (m-power   4)
                                     (h-power   0)
                                     (m-alpha   (* temp_adj (Kv1_amf v) ))
                                     (m-beta    (* temp_adj (Kv1_bmf v) )))
                                   
                                    ))


                        
                        (component (type pore)
                                   (const  gbar_Kv1  = 0.011)
                                   (output gbar_Kv1 ))
                        
                        (component (type permeating-substance) (name k)
                                   (const e_Kv1 = -85)
                                   (output e_Kv1 ))
                        
                        ) ;; end Kv1 current
             
             
             (component (type ion-channel) (name Kv4)

                        (component (type gate)
                                   
                                   ;; rate functions
                                   
                                   
                                   (defun Kv4_amf (v)  
                                     (let ((can    0.15743)
                                           (ckan   -32.19976)
                                           (cvan   57))
                                       (* can (exp (neg (/ (+ v cvan) ckan))))))
                                   

                                   (defun Kv4_bmf (v)  
                                     (let ((cbn   0.15743)
                                           (ckbn  37.51346)
                                           (cvbn  57))
                                       (* cbn (exp (neg (/ (+ v cvbn) ckbn))))))

                                   
                                   (defun Kv4_ahf (v)  
                                     (let ((cah   0.01342)
                                           (ckah  -7.86476)
                                           (cvah  60))
                                       (/ cah (+ 1.0 (exp (neg (/ (+ v cvah) ckah)))))))

                                   
                                   (defun Kv4_bhf (v)  
                                     (let ((cbh   0.04477)
                                           (ckbh  11.3615)
                                           (cvbh  54))
                                       (/ cbh (+ 1.0 (exp (neg (/ (+ v cvbh) ckbh)))))))

                                   (hh-ionic-conductance 
                                    (Kv4  ;; ion name: exported variables will be of the form {ion}_{id}
                                     (initial-m (/ (Kv4_amf Vrest) (+ (Kv4_amf Vrest) (Kv4_bmf Vrest))))
                                     (initial-h (/ (Kv4_ahf Vrest)  (+ (Kv4_ahf Vrest) (Kv4_bhf Vrest))) ) 
                                     (m-power   4)
                                     (h-power   1)
                                     (m-alpha   (* temp_adj (Kv4_amf v)))
                                     (m-beta    (* temp_adj (Kv4_bmf v)))
                                     (h-alpha   (* temp_adj (Kv4_ahf v)))
                                     (h-beta    (* temp_adj (Kv4_bhf v)))
                                     ))
                                   
                                   )

                        (component (type pore)
                                   (const  gbar_Kv4  = 0.0039)
                                   (output gbar_Kv4 ))
                        
                        (component (type permeating-substance) (name k)
                                   (const e_Kv4 = -85)
                                   (output e_Kv4 ))
                        
                        ) ;; end Kv4 current
             

             
             (component (type ion-channel) (name Ih)
                        
                        (component (type gate)
                                   
                                   ;; rate functions
                                   
                                   (defun Ih_inf (v)
                                     (let ((cvn   90.1)
                                           (ckn   -9.9))
                                       (/ 1.0 (+ 1.0 (exp (neg (/ (+ v cvn) ckn) ))))))
                                   
                                   (defun Ih_tau (v)   
                                     (let ((cct  190)
                                           (cat  720)
                                           (cvt  81.5)
                                           (ckt  11.9))
                                       (+ cct (* cat (exp (neg (pow (/ (+ v cvt) ckt) 2)))))))
                                   
                                   (hh-ionic-conductance 
                                    (Ih  ;; ion name: exported variables will be of the form {ion}_{id}
                                     (initial-m (Ih_inf Vrest))
                                     (m-power   1)
                                     (h-power   0)
                                     (m-inf     (Ih_inf v))
                                     (m-tau     (/ (Ih_tau v) temp_adj))
                                     ))
                                   
                                   )

                        (component (type pore)
                                   (const  gbar_Ih  = 0.0002)
                                   (output gbar_Ih ))
                        
                        (component (type permeating-substance) (name non-specific)
                                   (const e_Ih = -30)
                                   (output e_Ih ))
                        
                        ) ;; end Ih current

;;           (component (type ion-channel) (name CaP)
                        
;;                      (component (type gate)
                                   
;;                                 ;; rate functions
;;                                 (defun CaP_inf (v)  
;;                                   (let ((cv  19)
;;                                         (ck  5.5))
;;                                     (/ 1.0 (+ 1.0 (exp (neg (/ (+ v cv) ck)))))))
                                   
;;                                 (defun CaP_tau (v)  
;;                                   (if (> v -50) 
;;                                       (* 1e3 (+ 0.000191 (* 0.00376 (pow (exp (neg (/ (+ v 41.9) 27.8)))  2))))
;;                                       (* 1e3 (+ 0.00026367 (* 0.1278 (exp (* 0.10327 v)))))))

;;                                 (hh-ionic-conductance 
;;                                  (CaP  ;; ion name: exported variables will be of the form {ion}_{id}
;;                                   (initial-m  (CaP_inf Vrest))
;;                                   (m-power    1)
;;                                   (h-power    0)
;;                                   (m-inf      (CaP_inf v))
;;                                   (m-tau      (/ (CaP_tau v) temp_adj))))
                                   
;;                                 )

;;                      (component (type pore)
;;                                 (const gmax_CaP  = 0.01667)
;;                                 (gbar_CaP = (* gmax_CaP (ghk v celsius cai cao)))
;;                                 (output gbar_CaP ))
                        
;;                      (component (type accumulating-substance) (name ca) )
                        
                        
;;                      ) ;; end CaP current

             
             (component (type ion-channel) (name CaBK)
                        
                        (component (type gate)
                                   
                                   ;; rate functions
                                   
                                   (defun CaBK_zinf (ca)  
                                     (let ((zhalf 0.001))
                                       (/ 1 (+ 1 (/ zhalf ca)))))

                                   (const CaBK_ztau = 1.0)
                                   
                                   (defun CaBK_minf (v)
                                     (let ((cvm  28.9)
                                           (ckm  6.2))
                                       (/ 1.0 (+ 1.0 (exp (neg (/ (+ v 5.0 cvm) ckm)))))))
                                   
                                   (defun CaBK_mtau (v)
                                     (let ((ctm     0.000505)
                                           (cvtm1   86.4)
                                           (cktm1   -10.1)
                                           (cvtm2   -33.3)
                                           (cktm2   10))
                                       (+ ctm (/ 1.0 (+ (exp (neg (/ (+ v 5.0 cvtm1) cktm1))) 
                                                        (exp (neg (/ (+ v 5.0 cvtm2) cktm2))))))))
            
                                   (defun CaBK_hinf (v) 
                                     (let ((ch    0.085)
                                           (cvh   32)
                                           (ckh   -5.8))
                                       (+ ch (/ (- 1.0 ch) (+ 1.0 (exp (neg (/ (+ v 5.0 cvh) ckh))))))))
                                   
                                   (defun CaBK_htau (v)
                                     (let ((cth     0.0019)
                                           (cvth1   48.5)
                                           (ckth1   -5.2)
                                           (cvth2   -54.2)
                                           (ckth2   12.9))
                                       (+ cth (/ 1.0 ( + (exp (- (/ (+ v cvth1) ckth1))) 
                                                         (exp (- (/ (+ v cvth2) ckth2))))))))
                                   
                                   (state-complex
                                    (CaBK_z
                                     (transitions (-> zC zO  (/ (CaBK_zinf cai) CaBK_ztau))
                                                  (-> zO zC  (/ (- 1 (CaBK_zinf cai)) CaBK_ztau)))
                                     (initial   (CaBK_zinf 1e-4))
                                     (open zO)  (power 2)))
                                   
                                   (output CaBK_z )  

            
                                   (hh-ionic-conductance 
                                    (CaBK  ;; ion name: exported variables will be of the form {ion}_{id}
                                     (initial-m  (/ (CaBK_minf Vrest) temp_adj))
                                     (initial-h  (/ (CaBK_hinf Vrest) temp_adj ))
                                     (m-power    3)
                                     (h-power    1)
                                     (m-inf      (/ (CaBK_minf v) temp_adj ))
                                     (m-tau      (/ (CaBK_mtau v) temp_adj ))
                                     (h-inf      (/ (CaBK_hinf v) temp_adj) )
                                     (h-tau      (/ (CaBK_htau v) temp_adj) )))
                                   
                                   )
                        
                        (component (type pore)
                                   (const  gbar_CaBK  = 0.014)
                                   (output gbar_CaBK ))
                        
                        (component (type permeating-substance) (name k)
                                   (const e_CaBK = -85)
                                   (output e_CaBK ))
                        
                        ) ;; end BK current
             
             
             
             )) ;; end model