Model of arrhythmias in a cardiac cells network (Casaleggio et al. 2014)

 Download zip file   Auto-launch 
Help downloading and running models
" ... Here we explore the possible processes leading to the occasional onset and termination of the (usually) non-fatal arrhythmias widely observed in the heart. Using a computational model of a two-dimensional network of cardiac cells, we tested the hypothesis that an ischemia alters the properties of the gap junctions inside the ischemic area. ... In conclusion, our model strongly supports the hypothesis that non-fatal arrhythmias can develop from post-ischemic alteration of the electrical connectivity in a relatively small area of the cardiac cell network, and suggests experimentally testable predictions on their possible treatments."
1 . Casaleggio A, Hines ML, Migliore M (2014) Computational model of erratic arrhythmias in a cardiac cell network: the role of gap junctions. PLoS One 9:e100288 [PubMed]
Citations  Citation Browser
Model Information (Click on a link to find other models with that property)
Model Type: Realistic Network; Neuron or other electrically excitable cell;
Brain Region(s)/Organism:
Cell Type(s): Cardiac ventricular cell;
Channel(s): I K; I Sodium; I Calcium; I Potassium;
Gap Junctions: Gap junctions;
Simulation Environment: NEURON;
Model Concept(s): Spatio-temporal Activity Patterns; Detailed Neuronal Models; Action Potentials; Heart disease; Conductance distributions;
Implementer(s): Hines, Michael [Michael.Hines at]; Migliore, Michele [Michele.Migliore at];
Search NeuronDB for information about:  I K; I Sodium; I Calcium; I Potassium;

// network size parameters
default_var("Nx", 32) // number of cells in x-dimension
default_var("Ny", 64) // number of cells in y-dimension

// boundary conditions - 0 means cut, 1 means wrap, (2 means mirror?)
default_var("xwrap", 0)
default_var("ywrap", 0)

// gap connection parameter
default_var("ggap", 30) // nanosiemens
default_var("gsd", 0) // variance, nanosiemens
default_var("drift", 0) // drift

// ionic concentrations
default_var("ca_in", 1e-7) // used to set cai0_ca_ion (mM)
default_var("ca_out", 1.8) // used to set cao0_ca_ion (mM)
default_var("na_in", 20) // used to set nai0_na_ion (mM)
default_var("na_out", 140) // used to set nao0_na_ion (mM)
default_var("k_in", 141.59) // used to set ki0_k_ion (mM)
default_var("k_out", 5.4) // used to set ko0_k_ion (mM)

// stimulation parameters for NIClamp
default_var("msrow", 2) // beginning row for line of ms stim
default_var("mscol", 3) // beginning col for line of ms stim
default_var("nms", 4) // number of adjacent identical ms stim (in mscol)
default_var("msdel", 100) // (ms)
default_var("msdur", 1) // (ms)
default_var("msamp", 0) // (nA)
default_var("msdel1", 900) // (ms)
default_var("msn", 3) // (1)
default_var("ms2row", 8) // like msrow
default_var("ms2col", 8) // like mscol
default_var("ms2del", 100) // (ms)
default_var("ms2dur",5) // (msec)
default_var("ms2amp", 8) // (nA)
default_var("ms2del1", 800) // (msec)
default_var("ms2n", 200) // (1)

// simulation parameters
default_var("stoptime", 1000) // used to set tstop (ms)
default_var("temperature", 37) // used to set celsius (degC)
default_var("vinit", -75) // used to set v_init (mV)
default_var("ran_global", 0) // Random123_globalindex(ran_global)

// administrative parameters
default_var("use_halfgapspk", 0) // default HalfGap, otherwise, HalfGapSpk
default_var("checkpoint_interval", 10000) //interval for spike and trajec output
default_var("block_distrib", 0) // if 1 distribute in "square" blocks on ranks.
	// This reduces interprocessor communication but for greatest
	// effectiveness on a BlueGene one should also use an appropriate
	// rank to torus/core map

tstop = stoptime
celsius = temperature
v_init = vinit
cai0_ca_ion = ca_in
cao0_ca_ion = ca_out
nai0_na_ion = na_in
nao0_na_ion = na_out
ki0_k_ion = k_in
ko0_k_ion = k_out