Action potential of mouse urinary bladder smooth muscle (Mahapatra et al 2018)

 Download zip file   Auto-launch 
Help downloading and running models
Accession:243842
Urinary incontinence is associated with enhanced spontaneous phasic contractions of the detrusor smooth muscle (DSM). Although a complete understanding of the etiology of these spontaneous contractions is not yet established, it is suggested that the spontaneously evoked action potentials (sAPs) in DSM cells initiate and modulate the contractions. In order to further our understanding of the ionic mechanisms underlying sAP generation, we present here a biophysically detailed computational model of a single DSM cell. First, we constructed mathematical models for nine ion channels found in DSM cells based on published experimental data: two voltage-gated Ca2+ ion channels, an hyperpolarization-activated ion channel, two voltage-gated K+ ion channels, three Ca2+-activated K+ ion channels and a non-specific background leak ion channel. Incorporating these channels, our DSM model is capable of reproducing experimentally recorded spike-type sAPs of varying configurations, ranging from sAPs displaying after-hyperpolarizations to sAPs displaying after-depolarizations. Our model, constrained heavily by physiological data, provides a powerful tool to investigate the ionic mechanisms underlying the genesis of DSM electrical activity, which can further shed light on certain aspects of urinary bladder function and dysfunction.
Reference:
1 . Mahapatra C, Brain KL, Manchanda R (2018) A biophysically constrained computational model of the action potential of mouse urinary bladder smooth muscle. PLoS One 13:e0200712 [PubMed]
Citations  Citation Browser
Model Information (Click on a link to find other models with that property)
Model Type: Neuron or other electrically excitable cell;
Brain Region(s)/Organism:
Cell Type(s):
Channel(s): ATP-senstive potassium current; I Calcium; I h; I K,Ca; I K; I L high threshold; I T low threshold; IK Bkca; IK Skca;
Gap Junctions: Gap junctions;
Receptor(s):
Gene(s):
Transmitter(s):
Simulation Environment: NEURON;
Model Concept(s): Action Potentials; Action Potential Initiation; Calcium dynamics; Ion Channel Kinetics;
Implementer(s):
Search NeuronDB for information about:  I L high threshold; I T low threshold; I K; I h; I K,Ca; I Calcium; ATP-senstive potassium current; IK Bkca; IK Skca;
TITLE Voltage dependent Potassium channel (Kv) 
 
: Author: Chitaranjan Mahapatra (chitaranjan@iitb.ac.in)
: Computational Neurophysiology Lab
: Indian Institute of Technology Bombay, India 

: For details refer: 
: Mahapatra C, Brain KL, Manchanda R, A biophysically constrained computational model of the action potential 
: of mouse urinary bladder smooth muscle. PLOS One (2018) 



NEURON {
	SUFFIX KV2
	USEION k READ ek WRITE ik
	RANGE gkbar,ek, pt,qt,pw,vhi,sli ,vha,sla,ik
    GLOBAL finf,tauf ,ginf,taug, q10
	    
		}

UNITS {
	(molar) = (1/liter)
	(mM) = (millimolar)
	(mA) = (milliamp)
	(mV) = (millivolt)
}

PARAMETER {
	
    pt = 1
	qt = 1
	pw = 2
	vha = 1.1 (mV)
	sla= 11
	vhi = -58 (mV)
	sli= 15
    gkbar	= 0.006 (mho/cm2)
	ek	= -75	(mV)
	dt          (ms)
	v           (mV)
	celsius = 37(degC)
	q10 = 2.3
	
}

STATE {
	 f g   
}


ASSIGNED {
	
	     ik	   (mA/cm2)
         finf 
		 tauf
		  
		 ginf 
		 taug
		
        
         
}

INITIAL {
      
	  rates(v)
	  f=finf
	  g=ginf
	  	  
       }

BREAKPOINT {
	SOLVE states METHOD cnexp
	        ik = gkbar * f^pw*g* (v - ek)
			
			
  	        }
		
DERIVATIVE states {     
        rates(v)
        g' = (ginf - g)/taug
       f' =  (finf - f)/tauf		   

}

PROCEDURE rates(v(mV)) { LOCAL a,b,qtt
            qtt=q10^((celsius-37)/10)
	
        finf = 1/(1 + exp((vha-v)/sla)) 
        ginf = 1/(1 + exp((v-vhi)/sli))  
	    
          tauf = pt * qtt * (1/(1+(((v+15)/20)^2)))
	       taug = qt * (200+ (10/(1+(((v+54.18)/120)^5))))
}

FUNCTION trap0(v,th,a,q) {
	if (fabs(v-th) > 1e-6) {
	       trap0 = a * (v - th) / (1 - exp(-(v - th)/q))
	} else {
	        trap0 = a * q
 	}
}