Schematic of the simulation algorithm. At each global time step, there is sequential updating of (a) the stochastic dynamics of all couplons, (b) diffusion and buffering in all cytosolic voxels, and (c) membrane voltage and currents.