Modeling single neuron LFPs and extracellular potentials with LFPsim (Parasuram et al. 2016)

 Download zip file 
Help downloading and running models
Accession:190140
LFPsim - Simulation scripts to compute Local Field Potentials (LFP) from cable compartmental models of neurons and networks implemented in the NEURON simulation environment.
References:
1 . Parasuram H, Nair B, D'Angelo E, Hines M, Naldi G, Diwakar S (2016) Computational Modeling of Single Neuron Extracellular Electric Potentials and Network Local Field Potentials using LFPsim. Front Comput Neurosci 10:65 [PubMed]
2 . Diwakar S, Medini C, Nair M, Parasuram H, Vijayan A, Nair B (2017) Computational Neuroscience of Timing, Plasticity and Function in Cerebellum Microcircuits (Chapter 12) Computational Neurology and Psychiatry, Springer Series in Bio-/Neuroinformatics, √Črdi P:et al, ed. pp.343
Model Information (Click on a link to find other models with that property)
Model Type: Extracellular;
Brain Region(s)/Organism:
Cell Type(s): Cerebellum interneuron granule GLU cell; Hippocampus CA1 pyramidal GLU cell;
Channel(s):
Gap Junctions:
Receptor(s):
Gene(s):
Transmitter(s):
Simulation Environment: NEURON;
Model Concept(s): Extracellular Fields; Methods;
Implementer(s): Parasuram, Harilal [harilalp@am.amrita.edu]; Diwakar, Shyam [shyam at amrita.edu];
Search NeuronDB for information about:  Hippocampus CA1 pyramidal GLU cell; Cerebellum interneuron granule GLU cell;
: mea.mod

COMMENT
LFPsim - Simulation scripts to compute Local Field Potentials (LFP) from cable compartmental models of neurons and networks implemented in NEURON simulation environment.

LFPsim works reliably on biophysically detailed multi-compartmental neurons with ion channels in some or all compartments.

Last updated 12-March-2016
Developed by : Harilal Parasuram & Shyam Diwakar
Computational Neuroscience & Neurophysiology Lab, School of Biotechnology, Amrita University, India.
Email: harilalp@am.amrita.edu; shyam@amrita.edu
www.amrita.edu/compneuro 
ENDCOMMENT

NEURON {
	SUFFIX mea
	POINTER transmembrane_current_m
	RANGE mea_line0,mea_line1,mea_line2,mea_line3,mea_line4,mea_line5,mea_line6,mea_line7,mea_line8,mea_line9,mea_line10,mea_line11,mea_line12,mea_line13,mea_line14,mea_line15
	RANGE initial_part_line0,initial_part_line1,initial_part_line2,initial_part_line3,initial_part_line4,initial_part_line5,initial_part_line6,initial_part_line7,initial_part_line8,initial_part_line9,initial_part_line10,initial_part_line11,initial_part_line12,initial_part_line13,initial_part_line14,initial_part_line15	

}

PARAMETER {
	: default values put here

	}

ASSIGNED {
 
	transmembrane_current_m 
	initial_part_line0
	initial_part_line1
	initial_part_line2
	initial_part_line3
	initial_part_line4
	initial_part_line5
	initial_part_line6
	initial_part_line7
	initial_part_line8
	initial_part_line9
	initial_part_line10
	initial_part_line11
	initial_part_line12
	initial_part_line13
	initial_part_line14
	initial_part_line15

	mea_line0
	mea_line1
	mea_line2
	mea_line3
	mea_line4
	mea_line5
	mea_line6
	mea_line7
	mea_line8
	mea_line9
	mea_line10
	mea_line11
	mea_line12
	mea_line13
	mea_line14
	mea_line15


}

BREAKPOINT { 

	:Line Source Approximation
	mea_line0 = transmembrane_current_m * initial_part_line0 * 1e-1 	: 1e-1 (mA to uA) : calculated potential will be in uV
	mea_line1 = transmembrane_current_m * initial_part_line1 * 1e-1 
	mea_line2 = transmembrane_current_m * initial_part_line2 * 1e-1 
	mea_line3 = transmembrane_current_m * initial_part_line3 * 1e-1 
	mea_line4 = transmembrane_current_m * initial_part_line4 * 1e-1 
	mea_line5 = transmembrane_current_m * initial_part_line5 * 1e-1 
	mea_line6 = transmembrane_current_m * initial_part_line6 * 1e-1 
	mea_line7 = transmembrane_current_m * initial_part_line7 * 1e-1 
	mea_line8 = transmembrane_current_m * initial_part_line8 * 1e-1 
	mea_line9 = transmembrane_current_m * initial_part_line9 * 1e-1 
	mea_line10 = transmembrane_current_m * initial_part_line10 * 1e-1 
	mea_line11 = transmembrane_current_m * initial_part_line11 * 1e-1 
	mea_line12 = transmembrane_current_m * initial_part_line12 * 1e-1 
	mea_line13 = transmembrane_current_m * initial_part_line13 * 1e-1 
	mea_line14 = transmembrane_current_m * initial_part_line14 * 1e-1 
	mea_line15 = transmembrane_current_m * initial_part_line15 * 1e-1 

}


Loading data, please wait...