A two layers monodomain model of cardiac electrophysiology of the atria

Simulated P wave morphology in the presence of endo-epicardial activation delay

Aims

Evidences of asynchrony between epicardial and endocardial activation in the atrial wall have been reported. We used a computer model of the atria and torso to investigate the consequences of such activation delay on P wave morphology, while controlling for P wave duration.

Methods and results

We created 390 models of the atria based on the same geometry. These models differed by atrial wall thickness (from 2 to 3 mm), transmural coupling, and tissue conductivity in the endocardial and epicardial layers. Among them, 18 were in baseline, 186 had slower conduction in the epicardium layer and 186 in the endocardial layer. Conduction properties were adjusted in such a way that total activation time was the same in all models. P waves on a 16-lead system were simulated during sinus rhythm. Activation maps were similar in all cases. Endo-epicardial delay varied between -5.5 and 5.5 ms vs. 0 ± 0.5 ms in baseline. All P waves had the same duration but variability in their morphology was observed. With slower epicardial conduction, P wave amplitude was reduced by an average of 20% on leads V3-V5 and P wave area decreased by 50% on leads V1-V2 and by 40% on lead V3. Reversed, lower magnitude effects were observed with slower endocardial conduction.

Conclusion

An endo-epicardial delay of a few milliseconds is sufficient to significantly alter P wave morphology, even if the activation map remains the same.

Effects of Persistent Atrial Fibrillation-Induced Electrical Remodeling on Atrial Electro-Mechanics - Insights from a 3D Model of the Human Atria

Aims

Atrial stunning, a loss of atrial mechanical contraction, can occur following a successful cardioversion. It is hypothesized that persistent atrial fibrillation-induced electrical remodeling (AFER) on atrial electrophysiology may be responsible for such impaired atrial mechanics. This simulation study aimed to investigate the effects of AFER on atrial electro-mechanics.

Methods and Results

A 3D electromechanical model of the human atria was developed to investigate the effects of AFER on atrial electro-mechanics. Simulations were carried out in 3 conditions for 4 states: (i) the control condition, representing the normal tissue (state 1) and the tissue 2–3 months after cardioversion (state 2) when the atrial tissue recovers its electrophysiological properties after completion of reverse electrophysiological remodelling; (ii) AFER-SR condition for AF-remodeled tissue with normal sinus rhythm (SR) (state 3); and (iii) AFER-AF condition for AF-remodeled tissue with re-entrant excitation waves (state 4). Our results indicate that at the cellular level, AFER (states 3 & 4) abbreviated action potentials and reduced the Ca²⁺ content in the sarcoplasmic reticulum, resulting in a reduced amplitude of the intracellular Ca²⁺ transient leading to decreased cell active force and cell shortening as compared to the control condition (states 1 & 2). Consequently at the whole organ level, atrial contraction in AFER-SR condition (state 3) was dramatically reduced. In the AFER-AF condition (state 4) atrial contraction was almost abolished.

Conclusions

This study provides novel insights into understanding atrial electro-mechanics illustrating that AFER impairs atrial contraction due to reduced intracellular Ca²⁺ transients.

A two layers monodomain model of cardiac electrophysiology of the atria

Numerical simulations of the cardiac electrophysiology in the atria are often based on the standard bidomain or monodomain equations stated on a two-dimensional manifold. These simulations take advantage of the thinness of the atrial tissue, and their computational cost is reduced, as compared to three-dimensional simulations. However, these models do not take into account the heterogeneities located in the thickness of the tissue, like discontinuities of the fiber direction, although they can be a substrate for atrial arrhythmia (Hocini et al., Circulation 105(20):2442-2448, 2002; Ho et al., Cardiovasc Res 54(2):325-336, 2002; Nattel, Nature 415(6868):219-226, 2002). We investigate a two-dimensional model with two coupled, superimposed layers that allows to introduce three-dimensional heterogeneities, but retains a reasonable computational cost. We introduce the mathematical derivation of this model and error estimates with respect to the three-dimensional model. We give some numerical illustrations of its interest: we numerically show its convergence for vanishing thickness, introduce an optimization process of the coupling coefficient and assess its validity on physiologically relevant geometries. Our model would be an efficient tool to test the influence of three-dimensional fiber direction heterogeneities in reentries or atrial arrhythmia without using three-dimensional models.