# Calcium Dynamics in Atrial Tissues

written in

## Cell+Cell+Cell+…+Cell$\neq$Tissue!

This is an ongoing project, and it consists of investigating the emergent phenomena resulting from coupling cells (from the model presented here) by coupling the membrane potential $V_m(r,t)$ of every cell and distributing alternating-non alternating domains of cells and see what happens.

Cells are coupled by junctions as:

With:

• $\Delta r=0.025[\mu m]$ ,
• $dt=0.01$ $[ms]$ (so $[D]=[\mu m]^2/[ms]$).
• $T_{stim}=400$ $[ms]$ (Alternating regime)

Early versions of these experiments showed wave blockage and re-entry waves in the $V_m$ field. Although I haven’t considered real human atria geometries I have dozens of cool results by setting up a variety of arrangements for domains of sick and healthy cells.

The one I show in this post is very interesting, it exhibits discordant alternans, induced only by the geometry. This means that the alternans within the sick cells domain occurs completely out of phase.

For a better understanding of what I am trying to communicate check the movies below.

$[Ca^{2+}]_{i}(x,y,t):$

And, $V_m(x,y,t):$

This is still an open project, although some results have been published, we are still looking for a reduction of some dynamical aspects ot the RyR dynamics, as well as to realistic homeomorphims of human atrial geometries.

# Movie Scripts.

To conclude with the post, I’m going to briefly discuss the scripts I use to make the movies. These are tweaked versions of the examples discussed in the matplotlib documentation, and basically are created reading data files labeled using time containing a “picture” of a particular field variable at that given instant.

It requires the initial and final times, and the encoder FFMpeg installed, and of course python and matplotlib. The movie is encoded in mp4 format, for web usage perhaps is better to use “webm” (using ffmpeg -i V2D.mp4 V2D.webm , for example).