Non-Weak E-Fields Pyramidal Neurons (Reznik et. al.,2015)

 Download zip file 
Help downloading and running models
Accession:185512
Effect of Polarization Induced by Non-Weak Electric Fields on the Excitability of Elongated Neurons With Active Dendrite. In response to polarization, the active currents in the dendrites of pyramidal neurons play a pivotal role in the excitability of elongated neurons. Depending on a number of parameters either hyperpolarizing or depolarizing currents in the dendrite dominate as polarization is increased. Furthermore, the impact that these active dendrite channels (Ca, KAHP, etc) occur when only a small fraction of their channels are open.
Reference:
1 . Reznik RI, Barreto E, Sander E, So P (2015) Effects of polarization induced by non-weak electric fields on the excitability of elongated neurons with active dendrites. J Comput Neurosci [PubMed]
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): Hippocampus CA1 pyramidal cell; Hippocampus CA3 pyramidal cell; Neocortex V1 pyramidal corticothalamic L6 cell; Neocortex V1 pyramidal intratelencephalic L2-5 cell; Neocortex U1 pyramidal corticalthalamic L6 cell; Neocortex layer 4 pyramidal cell;
Channel(s): I h; I K,Ca;
Gap Junctions:
Receptor(s):
Gene(s):
Transmitter(s):
Simulation Environment: MATLAB;
Model Concept(s): Extracellular Fields;
Implementer(s): Barreto, Ernest ;
Search NeuronDB for information about:  Hippocampus CA1 pyramidal cell; Hippocampus CA3 pyramidal cell; Neocortex V1 pyramidal corticothalamic L6 cell; Neocortex V1 pyramidal intratelencephalic L2-5 cell; Neocortex U1 pyramidal corticalthalamic L6 cell; I h; I K,Ca;
  
%% MATLAB code to go with JCNS "Effect of Polarization Induced by Non-Weak Electric Fields on the Excitability of Elongated Neurons With Active Dendrite"
Robert I. Reznik  Ernest Barreto  Evelyn Sander  Paul So

RIR November 7, 2015

The paper explores the effects of polarization (VdsOut) on neural excitability. We use the two-compartment Pinsky-Rinzel model (1994)
In this work excitability is quantified by the TTFS of the neuron in response to stimulus.
The main stimulus used is ramp injected current (applied to soma). '
Later we use a model of synaptic AMPA applied to the dendrite instead of the ramp injection.
 Also we modify the PR neuron by adding an Ih current as modeled by Lippert and Booth (2009).

** I apologize for the extrememly long file names. I like to pack in a lot of metadata into the filenames themselves.

**************************************************************************

Main computational code

(1) Finding the equilibriums (if they exist)
	1.1 Use NumerEquilPR_db
		*input uses the PR structure
		*requires an initial guess but I found the guess can be quite off.
(2) Integrating the polarized PR neuron with ODE23(You are welcome to use other integrators I also used ODE45)
	2.1 In the folder Integration you will see a number of routines I would use
		2.1.1 SingIntegODE23PRWithInjCurr_db for use of the ramp protocol
		2.1.2 SingIntegODE23PRWithSyn_db for using the synaptic AMPA
		2.1.3 SingIntegODE23PRWithSynWithIntegParam_db is fairly versatile and can be used with a number of different protocols including syanptic and injected current,
		2.1.4 SingIntegODE23PRPlusHcurrWithSyn_db Synaptic AMPA and Ih current

	2.2 Also in the Integrate folder there is the initialize routine which uses a structure "aPR" and first populates it with default values. Later, you can set the values to whatever you want.
		2.2.1 IniPR_db

(3) The folder PR_WIh contains all of the additional routines needed to run code for the polarized PR with Ih current.	

***************************************************************	

General notes for the driver routines. All follow these general steps:

(1) Initialize parameters.
(2) Inside a loop over polarization (VdsOut)
	2.1 numerically search for an equilibrium (At this point there is no stimulus)
		(If it does not exist go to next increment in VdsOut loop)
	2.2 Intgrate ODE (with appropriate stimulus protocol)
		Default is to stop integration upon Vs passing up through 30 mV.
	2.3 optionally capture computational data such as currents and membrane potentials.
(3) There is the option, of course, to loop over additional parameters like Ek

************************************************
	
Some suggested Scripts and Figures to start with

(1) Fig4_ContIntmpt8 and Fig5_ContIntmpt3 reproduce figs 4 and 5 of paper.

(2) Plotting_ActiveDendCurr_AMPA_db
(3) TriipleLoopFixgAMPAVaryEkandgKAHP_db (This file shows a loop over 3 parameters for Polarized PR with AMPA
(4) PRwIhAMPALoopEkgKAHP.m (This one calls SingIntegODE23PRPlusHcurrWithSyn_db which uses both synaptic AMPA and the Ih currents) As currently set this might take awhile to run. The resulting dat is in computationaldata folder and figure 18 can be produced by running GridPlotLikeFig18PRwHJust3Grids.m
	

Reznik RI, Barreto E, Sander E, So P (2015) Effects of polarization induced by non-weak electric fields on the excitability of elongated neurons with active dendrites. J Comput Neurosci[PubMed]

References and models cited by this paper

References and models that cite this paper

Barreto E, Cressman JR (2011) Ion concentration dynamics as a mechanism for neuronal bursting Journal of Biological Physics 37:361-373 [Journal] [PubMed]

   Ion concentration dynamics as a mechanism for neuronal bursting (Barreto & Cressman 2011) [Model]

Berzhanskaya J, Gorchetchnikov A, Schiff SJ (2007) Switching between gamma and theta: dynamic network control using subthreshold electric fields Neurocomputing 70:10-12

Bikson M, Inoue M, Akiyama H, Deans JK, Fox JE, Miyakawa H, Jefferys JG (2004) Effects of uniform extracellular DC electric fields on excitability in rat hippocampal slices in vitro. J Physiol 557:175-90 [PubMed]

Bikson M, Name A, Rahman A (2013) Origins of specificity during tDCS: anatomical, activity-selective, and input-bias mechanisms. Front Hum Neurosci 7:688 [Journal] [PubMed]

Binder DK, Yao X, Zador Z, Sick TJ, Verkman AS, Manley GT (2006) Increased seizure duration and slowed potassium kinetics in mice lacking aquaporin-4 water channels. Glia 53:631-6 [Journal] [PubMed]

Bose A,Booth V (2004) Bursting in 2-compartment neurons: A case study of the Pinsky-Rinzel model

Cressman JR, Ullah G, Ziburkus J, Schiff SJ, Barreto E (2009) The influence of sodium and potassium dynamics on excitability, seizures, and the stability of persistent states: I. Single neuron dynamics. J Comput Neurosci 26:159-70 [Journal] [PubMed]

   Single neuron with dynamic ion concentrations (Cressman et al. 2009) [Model]

Cressman JR, Ullah G, Ziburkus J, Schiff SJ, Barreto E (2011) Erratum to: The influence of sodium and potassium dynamics on excitability, seizures, and the stability of persistent states: I. Single neuron dynamics. J Comput Neurosci 30:781

Csicsvari J, Jamieson B, Wise KD, Buzsaki G (2003) Mechanisms of gamma oscillations in the hippocampus of the behaving rat. Neuron 37:311-22 [PubMed]

Deans JK, Powell AD, Jefferys JG (2007) Sensitivity of coherent oscillations in rat hippocampus to AC electric fields. J Physiol 583:555-65 [Journal] [PubMed]

Dietzel I, Heinemann U, Lux HD (1989) Relations between slow extracellular potential changes, glial potassium buffering, and electrolyte and cellular volume changes during neuronal hyperactivity in cat brain. Glia 2:25-44 [PubMed]

Dyhrfjeld-Johnsen J, Morgan RJ, Soltesz I (2009) Double Trouble? Potential for Hyperexcitability Following Both Channelopathic up- and Downregulation of I(h) in Epilepsy. Front Neurosci 3:25-33 [PubMed]

Francis JT, Gluckman BJ, Schiff SJ (2003) Sensitivity of neurons to weak electric fields. J Neurosci 23:7255-61 [PubMed]

Gasparini S, DiFrancesco D (1999) Action of serotonin on the hyperpolarization-activated cation current (Ih) in rat CA1 hippocampal neurons. Eur J Neurosci 11:3093-100 [PubMed]

Ghai RS, Bikson M, Durand DM (2000) Effects of applied electric fields on low-calcium epileptiform activity in the CA1 region of rat hippocampal slices. J Neurophysiol 84:274-80 [PubMed]

Gluckman BJ, Neel EJ, Netoff TI, Ditto WL, Spano ML, Schiff SJ (1996) Electric field suppression of epileptiform activity in hippocampal slices. J Neurophysiol 76:4202-5 [PubMed]

Gluckman BJ, Nguyen H, Weinstein SL, Schiff SJ (2001) Adaptive electric field control of epileptic seizures. J Neurosci 21:590-600 [PubMed]

Gluckman BJ, So P, Netoff TI, Spano ML, Schiff SJ (1998) Stochastic resonance in mammalian neuronal networks. Chaos 8:588-598 [PubMed]

Golding NL, Mickus TJ, Katz Y, Kath WL, Spruston N (2005) Factors mediating powerful voltage attenuation along CA1 pyramidal neuron dendrites. J Physiol 568:69-82 [Journal] [PubMed]

   [22 reconstructed morphologies on NeuroMorpho.Org]
   Voltage attenuation in CA1 pyramidal neuron dendrites (Golding et al 2005) [Model]

Han CL, Hu W, Stead M, Zhang T, Zhang JG, Worrell GA, Meng FG (2014) Electrical stimulation of hippocampus for the treatment of refractory temporal lobe epilepsy. Brain Res Bull 109:13-21 [Journal] [PubMed]

Holt GR, Koch C (1999) Electrical interactions via the extracellular potential near cell bodies. J Comput Neurosci 6:169-84 [Journal] [PubMed]

   Extracellular Action Potential Simulations (Gold et al 2007) [Model]

Janigro D () Brain water and ion fluxes: a hard-to-die hypothesis to explain seizures. Epilepsy Curr 7:56-7 [Journal] [PubMed]

Lippert A, Booth V (2009) Understanding effects on excitability of simulated I (h) modulation in simple neuronal models. Biol Cybern 101:297-306 [Journal] [PubMed]

McNamara JO (1994) Cellular and molecular basis of epilepsy. J Neurosci 14:3413-25 [PubMed]

Mikkelsen R, Andreasen M, Nedergaard S (2013) Suppression of epileptiform activity by a single short-duration electric field in rat hippocampus in vitro. J Neurophysiol 109:2720-31 [Journal] [PubMed]

Miranda PC, Lomarev M, Hallett M (2006) Modeling the current distribution during transcranial direct current stimulation. Clin Neurophysiol 117:1623-9 [Journal] [PubMed]

Moody WJ, Futamachi KJ, Prince DA (1974) Extracellular potassium activity during epileptogenesis. Exp Neurol 42:248-63 [PubMed]

Park EH, Barreto E, Gluckman BJ, Schiff SJ, So P (2005) A Model of the Effects of Applied Electric Fields on Neuronal Synchronization J Comp Neurosci 19:53-70 [Journal]

Park EH, So P, Barreto E, Schiff SJ (2003) Electric field modulation of synchronization in neuronal networks Neurocomputing 52:169-175

Pinsky PF, Rinzel J (1994) Intrinsic and network rhythmogenesis in a reduced Traub model for CA3 neurons. J Comput Neurosci 1:39-60 [Journal] [PubMed]

   CA3 pyramidal cell: rhythmogenesis in a reduced Traub model (Pinsky, Rinzel 1994) [Model]

Pucihar G, Miklavcic D, Kotnik T (2009) A time-dependent numerical model of transmembrane voltage inducement and electroporation of irregularly shaped cells. IEEE Trans Biomed Eng 56:1491-501

Radman T, Ramos RL, Brumberg JC, Bikson M (2009) Role of cortical cell type and morphology in subthreshold and suprathreshold uniform electric field stimulation in vitro. Brain Stimul 2:215-28, 228.e1-3 [PubMed]

   [28 reconstructed morphologies on NeuroMorpho.Org]

Radman T, Su Y, An JH, Parra LC, Bikson M (2007) Spike timing amplifies the effect of electric fields on neurons: implications for endogenous field effects. J Neurosci 27:3030-6

Reato D, Rahman A, Bikson M, Parra LC (2010) Low-intensity electrical stimulation affects network dynamics by modulating population rate and spike timing. J Neurosci 30:15067-79

Richardson KA, Gluckman BJ, Weinstein SL, Glosch CE, Moon JB, Gwinn RP, Gale K, Schiff SJ (2003) In vivo modulation of hippocampal epileptiform activity with radial electric fields. Epilepsia 44:768-77

Richardson KA, Schiff SJ, Gluckman BJ (2005) Control of traveling waves in the Mammalian cortex. Phys Rev Lett 94:028103-92

Sunderam S, Chernyy N, Peixoto N, Mason JP, Weinstein SL, Schiff SJ, Gluckman BJ (2009) Seizure entrainment with polarizing low-frequency electric fields in a chronic animal epilepsy model. J Neural Eng 6:046009

Tranchina D, Nicholson C (1986) A model for the polarization of neurons by extrinsically applied electric fields. Biophys J 50:1139-56 [PubMed]

Traub RD, Dudek FE, Snow RW, Knowles WD (1985) Computer simulations indicate that electrical field effects contribute to the shape of the epileptiform field potential. Neuroscience 15:947-58

Traub RD, Dudek FE, Taylor CP, Knowles WD (1985) Simulation of hippocampal afterdischarges synchronized by electrical interactions. Neuroscience 14:1033-8 [PubMed]

Traub RD, Miles R (1991) Multiple modes of neuronal population activity emerge after modifying specific synapses in a model of the CA3 region of the hippocampus. Ann N Y Acad Sci 627:277-90 [PubMed]

Traub RD, Wong RK, Miles R, Michelson H (1991) A model of a CA3 hippocampal pyramidal neuron incorporating voltage-clamp data on intrinsic conductances. J Neurophysiol 66:635-50 [Journal] [PubMed]

Vigmond EJ, Perez Velazquez JL, Valiante TA, Bardakjian BL, Carlen PL (1997) Mechanisms of electrical coupling between pyramidal cells. J Neurophysiol 78:3107-16 [Journal] [PubMed]

Weiss SA, Faber DS (2010) Field effects in the CNS play functional roles. Front Neural Circuits 4:15 [Journal] [PubMed]

Yi GS, Wang J, Wei XL, Tsang KM, Chan WL, Deng B, Han CX (2014) Exploring how extracellular electric field modulates neuron activity through dynamical analysis of a two-compartment neuron model. J Comput Neurosci 36:383-99 [Journal] [PubMed]

(45 refs)