A Markov model of human Cav2.3 channels and their modulation by Zn2+ (Neumaier et al 2020)

 Download zip file 
Help downloading and running models
Accession:261714
The Markov model for Cav2.3 channel gating in the absence of trace metals was developed based on channel structure, previous modeling studies and the ability to fit the data. Model parameters were optimized by fitting the model to macroscopic currents recorded with various electrophysiological protocols from HEK-293 cells stably transfected with human Cav2.3+ß3 channel subunits. The effects of Zn2+ were implemented by assuming that Zn2+ binding to a first site (KZn=0.003 mM) leads to electrostatic modification and mechanical slowing of one of the voltage-sensors while Zn2+-binding to a second, intra-pore site (KZn=0.1 mM) blocks the channel and modifies the opening and closing transitions.
Reference:
1 . Neumaier F, Apldogan S, Hescheler J and Schneider T (2020) Zn2+-induced changes in Cav2.3 channel function: An electrophysiological and modeling study Journal of General Physiology, accepted
Model Information (Click on a link to find other models with that property)
Model Type: Channel/Receptor;
Brain Region(s)/Organism: Human;
Cell Type(s):
Channel(s): I Calcium;
Gap Junctions:
Receptor(s):
Gene(s):
Transmitter(s):
Simulation Environment: NEURON;
Model Concept(s): Markov-type model; Ion Channel Kinetics; Simplified Models; Drug binding;
Implementer(s):
Search NeuronDB for information about:  I Calcium;
strdef cmd
objref volt_cl			  
objref vbox1, vbox2, vbox3		 
objref g_v, g_i	          
objref ms

y_min = -0.2

create soma
soma {
	diam = 21.851       
	L    = 21.851 	
	nseg = 1 	     
	cm   = 1         
	Ra   = 60        
	insert CaR 
	cai0_ca_ion = 0.000001
	cao0_ca_ion = 4
	celsius = 22
	volt_cl = new VClamp_plus(.5)
}


vbox3 = new VBox()
vbox3.intercept(1)
ms = new MechanismStandard("CaR")
ms.action("change_CaR()")
ms.in()
ms.panel()
vbox3.intercept(0)
vbox3.map("CaR Control", 1, 854, -1, -1)


proc change_CaR() {
	ms.out()
}

proc start() {
	volt_cl.amp(0) = hold_pot
	volt_cl.amp(2) = hold_pot
	
	g_v.erase()
	g_i.erase()
	g_v.size(0, tstop, hold_pot-10, end_cl+10)
	g_i.size(0, tstop, y_min, 0.1)

	for (i=st_cl; i<=end_cl; i=i+incr1) {
       volt_cl.amp[1]=i
       
       finitialize(v_init)
       g_v.beginline()
       g_i.beginline()
       
       while (t<tstop) {
       		
			dens=soma.ica(0.5) 
       		
       		g_i.line(t, dens)
       		g_v.line(t, soma.v(.5))
       		
       		fadvance()
       }
    }
}


proc Tstop() {
	tstop = volt_cl.dur[0]+volt_cl.dur[1]+volt_cl.dur[2]
}

access soma
execute(cmd)
	
{
xpanel("Run Control", 0)
v_init = -80
xvalue("Init","v_init", 1,"stdinit()", 1, 1 )
xbutton("Init & Run","run()")
xbutton("Stop","stoprun=1")
runStopAt = 5
xvalue("Continue til","runStopAt", 1,"{continuerun(runStopAt) stoprun=1}", 1, 1 )
runStopIn = 1
xvalue("Continue for","runStopIn", 1,"{continuerun(t + runStopIn) stoprun=1}", 1, 1 )
xbutton("Single Step","steprun()")
t = 18
xvalue("t","t", 2 )
tstop = 415
xvalue("Tstop","tstop", 1,"tstop_changed()", 0, 1 )
dt = 0.25
xvalue("dt","dt", 1,"setdt()", 0, 1 )
steps_per_ms = 40
xvalue("Points plotted/ms","steps_per_ms", 1,"setdt()", 0, 1 )
screen_update_invl = 0.05
xvalue("Scrn update invl","screen_update_invl", 1,"", 0, 1 )
realtime = 0.00999999
xvalue("Real Time","realtime", 0,"", 0, 1 )
xpanel(1,102)
}

{
xpanel("Clamp Control", 0)
xbutton("Start","start()")
st_cl=-80
xvalue("Start voltage", "st_cl")
end_cl=60
xvalue("End voltage", "end_cl")
incr1=10
xvalue("Increment", "incr1")
volt_cl.dur[0]=5
xvalue("Pre-pulse duration", "volt_cl.dur[0]")
volt_cl.dur[1]=400
xvalue("Test pulse duration", "volt_cl.dur[1]")
volt_cl.dur[2]=10
xvalue("Post-pulse duration", "volt_cl.dur[2]")
hold_pot=v_init
xvalue("Holding voltage", "hold_pot")
xpanel(1,537)
}


vbox1=new VBox()
vbox1.intercept(1)
g_v=new Graph()
g_v.size(0, tstop, hold_pot-10, end_cl+10)
vbox1.intercept(0)
vbox1.map("Voltage protocol", 319, 24, 625, 260)


vbox2=new VBox()
vbox2.intercept(1)
g_i=new Graph()
g_i.size(0, tstop, y_min, 0.1)
vbox2.intercept(0)
vbox2.map("IV currents", 320, 456, 625, 304)

Loading data, please wait...