Background: Velocity and pattern of propagation of cardiac AP depends on structural andfunctional properties of the tissue, such as conductivity, dynamics of transmembrane ionicchannels and gap junctions (GJ). Gap junctions are clusters of channels that connect adjacent cells. A gap junction channel(GJC) is made of proteins named ­ connexins. Electrical behavior of GJCs depend on the typeand arrangement of their connexin composition. The dominating connexins in cardiacmyocytes are Cx43, Cx45 and Cx40. Methods and results: In current mathematical models, GJCs are considered to be passive.But, the experimental results, obtained by the dual­voltage clamp technique, show that GJCs display biophysical electrical properties such as voltage gating,i.e. a time and voltage dependence. Here we model Cx43 GJCs. We use the Hodgkin­Huxley formalism to describe GJCsconductance via one gating variable g j = g j (t, V j ). From our experimental results we obtain model parameters: the normalisedsteady state conductance and the time constant to reach the steady state, both voltagedependent. Once we have described the behavior of the single GJC, we write the mathematical model ofthe tissue, where we apply GJ current on specific parts of the cells’ membranes. Numerics and outlook: Some 3D numerical experiments are currently being performed on athin strip of cells, in order to compare the model’s results with the experimental ones. We use a monolayer of 50 × 3 cells, represented by cylinders of 100μm lengthand 10μm radius, with 2μm inter­cellular distance. We model GJCs on the cross sections of the cylinders. Finally, we apply an external stimulus on the border of the domain, and observe the propagation of theAP. Our goal is to make a mathematical model of the heterogeneous GJCs, including Cx45channels, as these have been shown to play a role in arrythmogenesis.Read less <