//genesis /*************************** MS Model, Version 5.2 ********************* **************************** CaQ_channel.g ********************* Tom Sheehan tsheeha2@gmu.edu thsheeha@vt.edu 703-538-8361 ****************************************************************************** ******************************************************************************/ function create_CaQ str chanName = "CaQ_channel" str compPath = "/library" int c = 0 float increment = 0.00005 float x = -0.1 int xdivs = 3000 float xmin = -0.1 float xmax = 0.05 float mPower = 1.0 float hPower = 0.0 float mvHalfCaQ = -9.0e-3 float mkCaQ = -6.6e-3 float mInfCaQ= 0.0 float mTauCaQ = 0.0 float surf = 0.0 float qFactCaN = {3.0/{CaSlow}} float gMax = 6.0e-008 float theta = 0.0 float theta_exp = 0.0 pushe {compPath} create tabchannel {chanName} setfield {chanName} Xpower {mPower} Ypower {hPower} call {chanName} TABCREATE X {xdivs} {xmin} {xmax} for(c = 0; c < {xdivs} + 1; c = c + 1) /************************ Begin CaQ_mTau *********************/ //mTauCaQ = 0.377e-3*ones(vDiv+1,1); mTauCaQ = 3.77e-004 setfield {chanName} X_A->table[{c}] {mTauCaQ} /************************ End CaQ_mTau ***********************/ /************************ Begin CaQ_mInf *********************/ //mInfCaQ = 1./(1 + exp((vMemb - mvHalfCaQ)/mkCaQ)); theta = {{x} - {mvHalfCaQ}}/{mkCaQ} theta_exp = {exp {theta}} + 1.0 mInfCaQ = 1.0/{theta_exp} setfield {chanName} X_B->table[{c}] {mInfCaQ} /************************ End CaQ_mInf ***********************/ x = x + increment end tweaktau {chanName} X create ghk {chanName}GHK setfield {chanName}GHK Cout 2 // Carter & Sabatini 2004 uses 2mM, // Wolf 5mM setfield {chanName}GHK valency 2.0 setfield {chanName}GHK T {TEMPERATURE} setfield {chanName} Gbar {gMax*surf} addmsg {chanName} {chanName}GHK PERMEABILITY Gk pope end