CA1 pyramidal neurons: effect of external electric field from power lines (Cavarretta et al. 2014)

 Download zip file   Auto-launch 
Help downloading and running models
Accession:151731
The paper discusses the effects induced by an electric field at power lines frequency.
Reference:
1 . Cavarretta F, Carnevale NT, Tegolo D, Migliore M (2014) Effects of low frequency electric fields on synaptic integration in hippocampal CA1 pyramidal neurons: implications for power line emissions. Front Cell Neurosci 8:310 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Dendrite; Extracellular;
Brain Region(s)/Organism: Hippocampus;
Cell Type(s): Hippocampus CA1 pyramidal cell;
Channel(s): I Na,t; I A; I K; I h;
Gap Junctions:
Receptor(s): AMPA;
Gene(s):
Transmitter(s): Glutamate;
Simulation Environment: NEURON;
Model Concept(s): Activity Patterns; Dendritic Action Potentials; Active Dendrites; Detailed Neuronal Models; Action Potentials; Synaptic Integration; Extracellular Fields;
Implementer(s): Carnevale, Ted [Ted.Carnevale at Yale.edu]; Migliore, Michele [Michele.Migliore at Yale.edu]; Cavarretta, Francesco [francescocavarretta at hotmail.it];
Search NeuronDB for information about:  Hippocampus CA1 pyramidal cell; AMPA; I Na,t; I A; I K; I h; Glutamate;
/
CavarrettaEtAl2014
readme.html
fzap.mod *
h.mod *
kadist.mod *
kaprox.mod *
kdrca1.mod *
na3n.mod *
naxn.mod *
netstimm.mod *
xtrau.mod *
anatscale.hoc *
biophys.hoc
calcd.hoc
calcrxcu.hoc
efheader.hoc *
geo5038804.hoc *
init.hoc
interpxyzu.hoc *
mosinit.hoc *
prova.ses
rigc.ses
screenshot.png
setnseg.hoc *
setpointersu.hoc *
synapses.hoc
voltage.ses
zapstimu.hoc
                            
celsius = 35

Rm = 28000
RmDend = Rm
RmSoma = Rm
RmAx = Rm

Cm    = 1
CmSoma = Cm
CmAx = Cm
CmDend = Cm

RaAll = 150
RaSoma =150  
RaAx = 50
Vrest =-65

forsec "axon" {
	insert pas e_pas=Vrest g_pas = 1/RmAx Ra=RaAx cm=CmAx
}

forsec "soma" {
	insert pas e_pas=Vrest g_pas = 1/RmSoma Ra=RaSoma cm=CmSoma
}

forsec "dendrite" {
	insert pas e_pas=Vrest g_pas = 1/RmDend Ra=RaAll cm=CmDend
}

forsec "user5" {
	insert pas e_pas=Vrest g_pas = 1/RmDend Ra=RaAll cm=CmDend
}


load_file("setnseg.hoc")

access soma
distance()

gna = 0.025
AXONM = 5
gkdr = 0.01
KMULT = 0.03
KMULTP =  0.03
ghd = 0.00005

forsec "axon" {   
	insert nax gbar_nax=gna*AXONM
	insert kdr gkdrbar_kdr=gkdr
	insert kap gkabar_kap=KMULTP
}

forsec "soma" {   
	insert hd ghdbar_hd=ghd vhalfl_hd=-73
	insert na3 gbar_na3=gna
	insert kdr gkdrbar_kdr=gkdr
	insert kap gkabar_kap=KMULTP				
}

for i = 0, numbasal - 1 dendrite[i] {
	insert hd ghdbar_hd=ghd vhalfl_hd=-73
	insert na3 gbar_na3=gna
	insert kdr gkdrbar_kdr=gkdr
	insert kap gkabar_kap=0
	insert kad gkabar_kad=0
	for(x) if(x > 0 && x < 1) {
		xdist = distance(x)
		ghdbar_hd(x) = ghd*(1+3*xdist/100)
		if(xdist > 100){
			vhalfl_hd=-81
			gkabar_kad(x) = KMULT*(1+xdist/100)
		} else {
			vhalfl_hd=-73
			gkabar_kap(x) = KMULTP*(1+xdist/100)
		}
	}
}

forsec "apical_dendrite" {
	insert hd ghdbar_hd=ghd
	insert na3 gbar_na3=gna
	insert kdr gkdrbar_kdr=gkdr
	insert kap gkabar_kap=0
	insert kad gkabar_kad=0
	for(x) if(x > 0 && x < 1) {
		xdist = distance(x)
		ghdbar_hd(x) = ghd*(1+3*xdist/100)
		if(xdist > 100) {
			vhalfl_hd=-81
			gkabar_kad(x) = KMULT*(1+xdist/100)
		} else {
			vhalfl_hd=-73
			gkabar_kap(x) = KMULTP*(1+xdist/100)
		}
	}
}

forsec "user5" {
	insert hd ghdbar_hd=ghd
	insert na3 gbar_na3=gna
	insert kdr gkdrbar_kdr=gkdr
	insert kap gkabar_kap=0
	insert kad gkabar_kad=0

	for(x) if(x > 0 && x < 1) {
		xdist = distance(x)
		ghdbar_hd(x) = ghd*(1+3*xdist/100)
		if(xdist > 100) {
			vhalfl_hd=-81
			gkabar_kad(x) = KMULT*(1+xdist/100)
		} else {
			vhalfl_hd=-73
			gkabar_kap(x) = KMULTP*(1+xdist/100)
		}
	}
}

proc init() {
	t = 0
	forall {
		v = Vrest
		if(ismembrane("nax") || ismembrane("na3")) {
			ena=55
		}
		if(ismembrane("kdr") || ismembrane("kap") || ismembrane("kad")) {
			ek=-90
		}
		if(ismembrane("hd")) {
			ehd_hd=-30
		}
	}
	finitialize(Vrest)
	fcurrent()
	forall {
		for(x) {
			if(ismembrane("na3") || ismembrane("nax")) {
				e_pas(x)=v(x)+(ina(x)+ik(x))/g_pas(x)
			}
			if(ismembrane("hd")) {
				e_pas(x)=e_pas(x)+i_hd(x)/g_pas(x)
			}
		}
	}
}

Loading data, please wait...