Interactions between calcium-induced arrhythmia triggers and the electrophysiological-anatomical substrate underlying the induction of atrial fibrillation.

Atrial fibrillation (AF) is the most common cardiac arrhythmia and is sustained by spontaneous focal excitations and re-entry. Spontaneous electrical firing in the pulmonary vein (PV) sleeves is implicated in AF generation. The aim of this simulation study was to identify the mechanisms determining the localisation of AF triggers in the PVs and their contribution to the genesis of AF. A novel biophysical model of the canine atria was used that integrates stochastic, spontaneous subcellular Ca2+ release events (SCRE) with regional electrophysiological heterogeneity in ionic properties and a detailed three-dimensional model of atrial anatomy, microarchitecture and patchy fibrosis. Simulations highlighted the importance of the smaller inward rectifier potassium current (IK1 ) in PV cells compared to the surrounding atria, which enabled SCRE more readily to result in delayed-afterdepolarisations that induced triggered activity. There was a leftward shift in the dependence of the probability of triggered activity on sarcoplasmic reticulum Ca2+ load. This feature was accentuated in 3D tissue compared to single cells (Δ half-maximal [Ca2+ ]SR  = 58 µM vs. 22 µM). In 3D atria incorporating electrical heterogeneity, excitations preferentially emerged from the PV region. These triggered focal excitations resulted in transient re-entry in the left atrium. Addition of fibrotic patches promoted localised emergence of focal excitations and wavebreaks that had a more substantial impact on generating AF-like patterns than the PVs. Thus, a reduced IK1 , less negative resting membrane potential, and fibrosis-induced changes of the electrotonic load all contribute to the emergence of complex excitation patterns from spontaneous focal triggers. KEY POINTS: Focal excitations in the atria are most commonly associated with the pulmonary veins, but the mechanisms for this localisation are yet to be elucidated. We applied a multi-scale computational modelling approach to elucidate the mechanisms underlying such localisations. Myocytes in the pulmonary vein region of the atria have a less negative resting membrane potential and reduced time-independent potassium current; we demonstrate that both of these factors promote triggered activity in single cells and tissues. The less negative resting membrane potential also contributes to heterogeneous inactivation of the fast sodium current, which can enable re-entrant-like excitation patterns to emerge without traditional conduction block.

Atrial Heterogeneity Generates Re-entrant Substrate during Atrial Fibrillation and Anti-arrhythmic Drug Action: Mechanistic Insights from Canine Atrial Models.

Anti-arrhythmic drug therapy is a frontline treatment for atrial fibrillation (AF), but its success rates are highly variable. This is due to incomplete understanding of the mechanisms of action of specific drugs on the atrial substrate at different stages of AF progression. We aimed to elucidate the role of cellular, tissue and organ level atrial heterogeneities in the generation of a re-entrant substrate during AF progression, and their modulation by the acute action of selected anti-arrhythmic drugs. To explore the complex cell-to-organ mechanisms, a detailed biophysical models of the entire 3D canine atria was developed. The model incorporated atrial geometry and fibre orientation from high-resolution micro-computed tomography, region-specific atrial cell electrophysiology and the effects of progressive AF-induced remodelling. The actions of multi-channel class III anti-arrhythmic agents vernakalant and amiodarone were introduced in the model by inhibiting appropriate ionic channel currents according to experimentally reported concentration-response relationships. AF was initiated by applied ectopic pacing in the pulmonary veins, which led to the generation of localized sustained re-entrant waves (rotors), followed by progressive wave breakdown and rotor multiplication in both atria. The simulated AF scenarios were in agreement with observations in canine models and patients. The 3D atrial simulations revealed that a re-entrant substrate was typically provided by tissue regions of high heterogeneity of action potential duration (APD). Amiodarone increased atrial APD and reduced APD heterogeneity and was more effective in terminating AF than vernakalant, which increased both APD and APD dispersion. In summary, the initiation and sustenance of rotors in AF is linked to atrial APD heterogeneity and APD reduction due to progressive remodelling. Our results suggest that anti-arrhythmic strategies that increase atrial APD without increasing its dispersion are effective in terminating AF.

Evolution and pharmacological modulation of the arrhythmogenic wave dynamics in canine pulmonary vein model.

AIMS: Atrial fibrillation (AF), the commonest cardiac arrhythmia, has been strongly linked with arrhythmogenic sources near the pulmonary veins (PVs), but underlying mechanisms are not fully understood. We aim to study the generation and sustenance of wave sources in a model of the PV tissue. METHODS AND RESULTS: A previously developed biophysically detailed three-dimensional canine atrial model is applied. Effects of AF-induced electrical remodelling are introduced based on published experimental data, as changes of ion channel currents (ICaL, IK1, Ito, and IKur), the action potential (AP) and cell-to-cell coupling levels. Pharmacological effects are introduced by blocking specific ion channel currents. A combination of electrical heterogeneity (AP tissue gradients of 5-12 ms) and anisotropy (conduction velocities of 0.75-1.25 and 0.21-0.31 m/s along and transverse to atrial fibres) can results in the generation of wave breaks in the PV region. However, a long wavelength (171 mm) prevents the wave breaks from developing into re-entry. Electrical remodelling leads to decreases in the AP duration, conduction velocity and wavelength (to 49 mm), such that re-entry becomes sustained. Pharmacological effects on the tissue heterogeneity and vulnerability (to wave breaks and re-entry) are quantified to show that drugs that increase the wavelength and stop re-entry (IK1 and IKur blockers) can also increase the heterogeneity (AP gradients of 26-27 ms) and the likelihood of wave breaks. CONCLUSION: Biophysical modelling reveals large conduction block areas near the PVs, which are due to discontinuous fibre arrangement enhanced by electrical heterogeneity. Vulnerability to re-entry in such areas can be modulated by pharmacological interventions.

Pro-arrhythmogenic effects of atrial fibrillation-induced electrical remodelling: insights from the three-dimensional virtual human atria.

Chronic atrial fibrillation (AF) is associated with structural and electrical remodelling in the atria, which are associated with a high recurrence of AF. Through biophysically detailed computer modelling, this study investigated mechanisms by which AF-induced electrical remodelling promotes and perpetuates AF. A family of Courtemanche-Ramirez-Nattel variant models of human atrial cell action potentials (APs), taking into account of intrinsic atrial electrophysiological properties, was modified to incorporate various experimental data sets on AF-induced changes of major ionic channel currents (ICaL, IKur, Ito, IK1, IKs, INaCa) and on intracellular Ca(2+) handling. The single cell models for control and AF-remodelled conditions were incorporated into multicellular three-dimensional (3D) atrial tissue models. Effects of the AF-induced electrical remodelling were quantified as the changes of AP profile, AP duration (APD) and its dispersion across the atria, and the vulnerability of atrial tissue to the initiation of re-entry. The dynamic behaviour of re-entrant excitation waves in the 3D models was characterised. In our simulations, AF-induced electrical remodelling abbreviated atrial APD non-uniformly across the atria; this resulted in relatively short APDs co-existing with marked regional differences in the APD at junctions of the crista terminalis/pectinate muscle, pulmonary veins/left atrium. As a result, the measured tissue vulnerability to re-entry initiation at these tissue junctions was increased. The AF-induced electrical remodelling also stabilized and accelerated re-entrant excitation waves, leading to rapid and sustained re-entry. Under the AF-remodelled condition, re-entrant scroll waves in the 3D model degenerated into persistent and erratic wavelets, leading to fibrillation. In conclusion, realistic 3D atrial tissue models indicate that AF-induced electrical remodelling produces regionally heterogeneous and shortened APD; these respectively facilitate initiation and maintenance of re-entrant excitation waves.

Patchy fibrosis promotes trigger-substrate interactions that both generate and maintain atrial fibrillation.

Fibrosis has been mechanistically linked to arrhythmogenesis in multiple cardiovascular conditions, including atrial fibrillation (AF). Previous studies have demonstrated that fibrosis can create functional barriers to conduction which may promote excitation wavebreak and the generation of re-entry, while also acting to pin re-entrant excitation in stable rotors during AF. However, few studies have investigated the role of fibrosis in the generation of AF triggers in detail. We apply our in-house computational framework to study the impact of fibrosis on the generation of AF triggers and trigger-substrate interactions in two- and three-dimensional atrial tissue models. Our models include a reduced and efficient description of stochastic, spontaneous cellular triggers as well as a simple model of heterogeneous inter-cellular coupling. Our results demonstrate that fibrosis promotes the emergence of focal excitations, primarily through reducing the electrotonic load on individual fibre strands. This enables excitation to robustly initiate within these single strands before spreading to neighbouring strands and inducing a full tissue focal excitation. Enhanced conduction block can allow trigger-substrate interactions that result in the emergence of complex, re-entrant excitation patterns. This study provides new insight into the mechanisms by which fibrosis promotes the triggers and substrate necessary to induce and sustain arrhythmia.

Postnatal development of transmural gradients in expression of ion channels and Ca²âº-handling proteins in the ventricle.

Transmural gradients in myocyte action potential duration (APD) and Ca(2+)-handling proteins are argued to be important for both the normal functioning of the ventricle and arrhythmogenesis. In rabbit, the transmural gradient in APD (left ventricular wedge preparation) is minimal in the neonate. During postnatal development, APD increases both in the epicardium and the endocardium, but the prolongation is more substantial in the endocardium leading to a significant transmural gradient. We have investigated changes in the expression of ion channels and also Ca(2+)-handling proteins in the subepicardial and subendocardial layers of the left ventricular free wall in neonatal (2-7 days of age) and adult male (~6 months of age) New Zealand White rabbits using quantitative PCR and also, when possible, in situ hybridisation and immunohistochemistry. In the adult, there were significant and substantial transmural gradients in Ca(v)1.2, KChIP2, ERG, K(v)LQT1, K(ir)2.1, NCX1, SERCA2a and RyR2 at the mRNA and, in some cases, protein level-in every case the mRNA or protein was more abundant in the epicardium than the endocardium. Of the eight transmural gradients seen in the adult, only three were observed in the neonate and, in two of these cases, the gradients were smaller than those in the adult. However, in the neonate there were also transmural gradients not observed in the adult: in HCN4, Na(v)1.5, minK, K(ir)3.1 and Cx40 mRNAs - in every case the mRNA was more abundant in the endocardium than the epicardium. If the postnatal changes in ion channel mRNAs are used to predict changes in ionic conductances, mathematical modelling predicts the changes in APD observed experimentally. It is concluded that many of the well known transmural gradients in the ventricle develop postnatally.

Mechanistic links between Na+ channel (SCN5A) mutations and impaired cardiac pacemaking in sick sinus syndrome.

RATIONALE: Familial sick sinus syndrome (SSS) has been linked to loss-of-function mutations of the SCN5A gene, which result in decreased inward Na(+) current, I(Na). However, the functional role of I(Na) in cardiac pacemaking is controversial, and mechanistic links between mutations and sinus node dysfunction in SSS are unclear. OBJECTIVE: To determine mechanisms by which the SCN5A mutations impair cardiac pacemaking. METHODS AND RESULTS: Action potential (AP) models for rabbit sinoatrial node (SAN) cells were modified to incorporate experimentally reported I(Na) changes induced by 2 groups of SCN5A gene mutations (affecting the activation and inactivation of I(Na), respectively). The cell models were incorporated into an anatomically detailed 2D model of the intact SAN-atrium. Effects of the mutations and vagal nerve activity on cardiac pacemaking at the single-cell and tissue levels were studied. Multielectrode extracellular potential recordings of activation pattern from intact SAN-atrium preparations were performed to test predictions of the models. At the single-cell level, the mutations slowed down pacemaking rates in peripheral, but not in central SAN cells that control the heart rhythm. However, in tissue simulations, the mutations not only slowed down pacemaking, but also compromised AP conduction across the SAN-atrium, leading to a possible SAN exit block or sinus arrest, the major features of SSS. Simulated vagal nerve activity amplified the bradycardiac effects of the mutations. Two groups of SCN5A mutations showed subtle differences in impairing the ability of the SAN to drive the surrounding atrium, primarily attributable to their differential effects on atrial excitability and conduction safety. Experimental data with tetrodotoxin and carbachol confirmed the simulation outcomes. CONCLUSIONS: Our study substantiates the causative link between SCN5A gene mutations and SSS and illustrates mechanisms by which the mutations impair the driving ability of the SAN.

Ionic mechanisms for electrical heterogeneity between rabbit Purkinje fiber and ventricular cells.

The intrinsic heterogeneity of electrical action potential (AP) properties between Purkinje fibers (PFs) and the ventricular wall, as well as within the wall, plays an important role in ensuring successful excitation of the ventricles. It can also be proarrhythmic due to nonuniform repolarization across the Purkinje-ventricular junction. However, the ionic mechanisms that underlie the marked AP differences between PFs and ventricular cells are not fully characterized. We studied such mechanisms by developing a new family of biophysically detailed AP models for rabbit PF cells and three transmural ventricular cell types. The models were based on and validated against experimental data recorded from rabbit at ionic channel, single cell, and tissue levels. They were then used to determine the functional roles of each individual ionic channel current in modulating the AP heterogeneity at the rabbit Purkinje-ventricular junction, and to identify specific currents responsible for the differential response of PFs and ventricular cells to pharmacological interventions.

Optimal velocity and safety of discontinuous conduction through the heterogeneous Purkinje-ventricular junction.

Slow and discontinuous wave conduction through nonuniform junctions in cardiac tissues is generally considered unsafe and proarrythmogenic. However, the relationships between tissue structure, wave conduction velocity, and safety at such junctions are unknown. We have developed a structurally and electrophysiologically detailed model of the canine Purkinje-ventricular junction (PVJ) and varied its heterogeneity parameters to determine such relationships. We show that neither very fast nor very slow conduction is safe, and there exists an optimal velocity that provides the maximum safety factor for conduction through the junction. The resultant conduction time delay across the PVJ is a natural consequence of the electrophysiological and morphological differences between the Purkinje fiber and ventricular tissue. The delay allows the PVJ to accumulate and pass sufficient charge to excite the adjacent ventricular tissue, but is not long enough for the source-to-load mismatch at the junction to be enhanced over time. The observed relationships between the conduction velocity and safety factor can provide new insights into optimal conditions for wave propagation through nonuniform junctions between various cardiac tissues.

Mechanisms of transition from normal to reentrant electrical activity in a model of rabbit atrial tissue: interaction of tissue heterogeneity and anisotropy.

Experimental evidence suggests that regional differences in action potential (AP) morphology can provide a substrate for initiation and maintenance of reentrant arrhythmias in the right atrium (RA), but the relationships between the complex electrophysiological and anatomical organization of the RA and the genesis of reentry are unclear. In this study, a biophysically detailed three-dimensional computer model of the right atrial tissue was constructed to study the role of tissue heterogeneity and anisotropy in arrhythmogenesis. The model of Lindblad et al. for a rabbit atrial cell was modified to incorporate experimental data on regional differences in several ionic currents (primarily, I(Na), I(CaL), I(K1), I(to), and I(sus)) between the crista terminalis and pectinate muscle cells. The modified model was validated by its ability to reproduce the AP properties measured experimentally. The anatomical model of the rabbit RA (including tissue geometry and fiber orientation) was based on a recent histological reconstruction. Simulations with the resultant electrophysiologically and anatomically detailed three-dimensional model show that complex organization of the RA tissue causes breakdown of regular AP conduction patterns at high pacing rates (>11.75 Hz): as the AP in the crista terminalis cells is longer, and electrotonic coupling transverse to fibers of the crista terminalis is weak, high-frequency pacing at the border between the crista terminalis and pectinate muscles results in a unidirectional conduction block toward the crista terminalis and generation of reentry. Contributions of the tissue heterogeneity and anisotropy to reentry initiation mechanisms are quantified by measuring action potential duration (APD) gradients at the border between the crista terminalis and pectinate muscles: the APD gradients are high in areas where both heterogeneity and anisotropy are high, such that intrinsic APD differences are not diminished by electrotonic interactions. Thus, our detailed computer model reconstructs complex electrical activity in the RA, and provides new insights into the mechanisms of transition from focal atrial tachycardia into reentry.

The canine virtual ventricular wall: a platform for dissecting pharmacological effects on propagation and arrhythmogenesis.

We have constructed computational models of canine ventricular cells and tissues, ultimately combining detailed tissue architecture and heterogeneous transmural electrophysiology. The heterogeneity is introduced by modifying the Hund-Rudy canine cell model in order to reproduce experimentally reported electrophysiological properties of endocardial, midmyocardial (M) and epicardial cells. These models are validated against experimental data for individual ionic current and action potential characteristics, and their rate dependencies. 1D and 3D heterogeneous virtual tissues are constructed, with detailed tissue architecture (anisotropy and orthotropy, due to fibre orientation and sheet structure) of the left ventricular wall wedge extracted from a diffusion tensor imaging data set. The models are used to study the effects of tissue heterogeneity and class III drugs on transmural propagation and tissue vulnerability to re-entry. We have determined relationships between the transmural dispersion of action potential duration (APD) and the vulnerable window in the 1D virtual ventricular wall, and demonstrated how changes in the transmural heterogeneity, and hence tissue vulnerability, can lead to generation of re-entry in the 3D ventricular wedge. Two class III drugs with opposite qualitative effects on transmural APD heterogeneity are considered: d-sotalol that increases transmural APD dispersion, and amiodarone that decreases it. Simulations with the 1D virtual ventricular wall show that under d-sotalol conditions the vulnerable window is substantially wider compared to amiodarone conditions, primarily in the epicardial region where unidirectional conduction block persists until the adjacent M cells are fully repolarised. Further simulations with the 3D ventricular wedge have shown that ectopic stimulation of the epicardial region results in generation of sustained re-entry under d-sotalol conditions, but not under amiodarone conditions or in control. Again, APD increase in M cells was identified as the major contributor to tissue vulnerability--re-entry was initiated primarily due to ectopic excitation propagating around the unidirectional conduction block in the M cell region. This suggests an electrophysiological mechanism for the anti- and proarrhythmic effects of the class III drugs: the relative safety of amiodarone in comparison to d-sotalol can be explained by relatively low transmural APD dispersion, and hence, a narrow vulnerable window and low probability of re-entry in the tissue.

Novel Computational Analysis of Left Atrial Anatomy Improves Prediction of Atrial Fibrillation Recurrence after Ablation.

The left atrium (LA) can change in size and shape due to atrial fibrillation (AF)-induced remodeling. These alterations can be linked to poorer outcomes of AF ablation. In this study, we propose a novel comprehensive computational analysis of LA anatomy to identify what features of LA shape can optimally predict post-ablation AF recurrence. To this end, we construct smooth 3D geometrical models from the segmentation of the LA blood pool captured in pre-procedural MR images. We first apply this methodology to characterize the LA anatomy of 144 AF patients and build a statistical shape model that includes the most salient variations in shape across this cohort. We then perform a discriminant analysis to optimally distinguish between recurrent and non-recurrent patients. From this analysis, we propose a new shape metric called vertical asymmetry, which measures the imbalance of size along the anterior to posterior direction between the superior and inferior left atrial hemispheres. Vertical asymmetry was found, in combination with LA sphericity, to be the best predictor of post-ablation recurrence at both 12 and 24 months (area under the ROC curve: 0.71 and 0.68, respectively) outperforming other shape markers and any of their combinations. We also found that model-derived shape metrics, such as the anterior-posterior radius, were better predictors than equivalent metrics taken directly from MRI or echocardiography, suggesting that the proposed approach leads to a reduction of the impact of data artifacts and noise. This novel methodology contributes to an improved characterization of LA organ remodeling and the reported findings have the potential to improve patient selection and risk stratification for catheter ablations in AF.

Novel MRI Technique Enables Non-Invasive Measurement of Atrial Wall Thickness.

Knowledge of atrial wall thickness (AWT) has the potential to provide important information for patient stratification and the planning of interventions in atrial arrhythmias. To date, information about AWT has only been acquired in post-mortem or poor-contrast computed tomography (CT) studies, providing limited coverage and highly variable estimates of AWT. We present a novel contrast agent-free MRI sequence for imaging AWT and use it to create personalized AWT maps and a biatrial atlas. A novel black-blood phase-sensitive inversion recovery protocol was used to image ten volunteers and, as proof of concept, two atrial fibrillation patients. Both atria were manually segmented to create subject-specific AWT maps using an average of nearest neighbors approach. These were then registered non-linearly to generate an AWT atlas. AWT was 2.4 ± 0.7 and 2.7 ± 0.7 mm in the left and right atria, respectively, in good agreement with post-mortem and CT data, where available. AWT was 2.6 ± 0.7 mm in the left atrium of a patient without structural heart disease, similar to that of volunteers. In a patient with structural heart disease, the AWT was increased to 3.1 ± 1.3 mm. We successfully designed an MRI protocol to non-invasively measure AWT and create the first whole-atria AWT atlas. The atlas can be used as a reference to study alterations in thickness caused by atrial pathology. The protocol can be used to acquire personalized AWT maps in a clinical setting and assist in the treatment of atrial arrhythmias.

Towards patient-specific modelling of lesion formation during radiofrequency catheter ablation for atrial fibrillation.

Radiofrequency catheter ablation procedures are a first-line method of clinical treatment for atrial fibrillation. However, they suffer from suboptimal success rates and are also prone to potentially serious adverse effects. These limitations can be at least partially attributed to the inter- and intra- patient variations in atrial wall thickness, and could be mitigated by patient-specific approaches to the procedure. In this study, a modelling approach to optimising ablation procedures in subject-specific 3D atrial geometries was applied. The approach enabled the evaluation of optimal ablation times to create lesions for a given wall thickness measured from MRI. A nonliner relationship was revealed between the thickness and catheter contact time required for fully transmural lesions. Hence, our approach based on MRI reconstruction of the atrial wall combined with subject-specific modelling of ablation can provide useful information for improving clinical procedures.

Slow Conduction in the Border Zones of Patchy Fibrosis Stabilizes the Drivers for Atrial Fibrillation: Insights from Multi-Scale Human Atrial Modeling.

Introduction: The genesis of atrial fibrillation (AF) and success of AF ablation therapy have been strongly linked with atrial fibrosis. Increasing evidence suggests that patient-specific distributions of fibrosis may determine the locations of electrical drivers (rotors) sustaining AF, but the underlying mechanisms are incompletely understood. This study aims to elucidate a missing mechanistic link between patient-specific fibrosis distributions and AF drivers. Methods: 3D atrial models integrated human atrial geometry, rule-based fiber orientation, region-specific electrophysiology, and AF-induced ionic remodeling. A novel detailed model for an atrial fibroblast was developed, and effects of myocyte-fibroblast (M-F) coupling were explored at single-cell, 1D tissue and 3D atria levels. Left atrial LGE MRI datasets from 3 chronic AF patients were segmented to provide the patient-specific distributions of fibrosis. The data was non-linearly registered and mapped to the 3D atria model. Six distinctive fibrosis levels (0-healthy tissue, 5-dense fibrosis) were identified based on LGE MRI intensity and modeled as progressively increasing M-F coupling and decreasing atrial tissue coupling. Uniform 3D atrial model with diffuse (level 2) fibrosis was considered for comparison. Results: In single cells and tissue, the largest effect of atrial M-F coupling was on the myocyte resting membrane potential, leading to partial inactivation of sodium current and reduction of conduction velocity (CV). In the 3D atria, further to the M-F coupling, effects of fibrosis on tissue coupling greatly reduce atrial CV. AF was initiated by fast pacing in each 3D model with either uniform or patient-specific fibrosis. High variation in fibrosis distributions between the models resulted in varying complexity of AF, with several drivers emerging. In the diffuse fibrosis models, waves randomly meandered through the atria, whereas in each the patient-specific models, rotors stabilized in fibrotic regions. The rotors propagated slowly around the border zones of patchy fibrosis (levels 3-4), failing to spread into inner areas of dense fibrosis. Conclusion: Rotors stabilize in the border zones of patchy fibrosis in 3D atria, where slow conduction enable the development of circuits within relatively small regions. Our results can provide a mechanistic explanation for the clinical efficacy of ablation around fibrotic regions.

Role of atrial tissue substrate and electrical activation pattern in fractionation of atrial electrograms: a computational study.

Complex fractionated atrial electrograms (CFAEs) are often used as a clinical marker for re-entrant drivers of atrial fibrillation. However, outcomes of clinical ablation procedures based on CFAEs are controversial and the mechanistic links between fractionation, re-entrant activity and the characteristics of the atrial substrate are not completely understood. We explore such links by simulating electrograms arising from both normal and re-entrant electrical activity in atrial tissue models. 2D and 3D tissue geometries with a range of conditions for intracellular coupling and myofiber orientation fields were studied. Electrograms were fractionated in the presence of complex atrial fiber fields and in 3D irregular geometries, due to far-field excitations. The complexity of the local electrical activity was not a strong determinant of the degree of fractionation. These results suggest that electrogram fractionation is more strongly linked to atrial substrate characteristics (including tissue geometry, fiber orientation and degree of intercelullar coupling) than to the electrical activation pattern sustaining atrial fibrillation.

A novel computational sheep atria model for the study of atrial fibrillation.

Sheep are often used as animal models for experimental studies into the underlying mechanisms of cardiac arrhythmias. Previous studies have shown that biophysically detailed computer models of the heart provide a powerful alternative to experimental animal models for underpinning such mechanisms. In this study, we have developed a family of mathematical models for the electrical action potentials of various sheep atrial cell types. The developed cell models were then incorporated into a three-dimensional anatomical model of the sheep atria, which was recently reconstructed and segmented based on anatomical features within different regions. This created a novel biophysically detailed computational model of the three-dimensional sheep atria. Using the model, we then investigated the mechanisms by which paroxysmal rapid focal activity in the pulmonary veins can transit to sustained atrial fibrillation. It was found that the anisotropic property of the atria arising from the fibre structure plays an important role in facilitating the development of fibrillatory atrial excitation waves, and the electrical heterogeneity plays an important role in its initiation.

Heterogeneous and anisotropic integrative model of pulmonary veins: computational study of arrhythmogenic substrate for atrial fibrillation.

Mechanisms underlying the genesis of re-entrant substrate for the most common cardiac arrhythmia, atrial fibrillation (AF), are not well understood. In this study, we develop a multi-scale three-dimensional computational model that integrates cellular electrophysiology of the left atrium (LA) and pulmonary veins (PVs) with the respective tissue geometry and fibre orientation. The latter is reconstructed in unique detail from high-resolution (approx. 70 µm) contrast micro-computed tomography data. The model is used to explore the mechanisms of re-entry initiation and sustenance in the PV region, regarded as the primary source of high-frequency electrical activity in AF. Simulations of the three-dimensional model demonstrate that an initial break-down of normal electrical excitation wave-fronts can be caused by the electrical heterogeneity between the PVs and LA. High tissue anisotropy is then responsible for the slow conduction and generation of a re-entrant circuit near the PVs. Evidence of such circuits has been seen clinically in AF patients. Our computational study suggests that primarily the combination of electrical heterogeneity and conduction anisotropy between the PVs and LA tissues leads to the generation of a high-frequency (approx. 10 Hz) re-entrant source near the PV sleeves, thus providing new insights into the arrhythmogenic mechanisms of excitation waves underlying AF.

Application of micro-computed tomography with iodine staining to cardiac imaging, segmentation, and computational model development.

Micro-computed tomography (micro-CT) has been widely used to generate high-resolution 3-D tissue images from small animals nondestructively, especially for mineralized skeletal tissues. However, its application to the analysis of soft cardiovascular tissues has been limited by poor inter-tissue contrast. Recent ex vivo studies have shown that contrast between muscular and connective tissue in micro-CT images can be enhanced by staining with iodine. In the present study, we apply this novel technique for imaging of cardiovascular structures in canine hearts. We optimize the method to obtain high-resolution X-ray micro-CT images of the canine atria and its distinctive regions-including the Bachmann's bundle, atrioventricular node, pulmonary arteries and veins-with clear inter-tissue contrast. The imaging results are used to reconstruct and segment the detailed 3-D geometry of the atria. Structure tensor analysis shows that the arrangement of atrial fibers can also be characterized using the enhanced micro-CT images, as iodine preferentially accumulates within the muscular fibers rather than in connective tissues. This novel technique can be particularly useful in nondestructive imaging of 3-D cardiac architectures from large animals and humans, due to the combination of relatively high speed ( ~ 1 h/per scan of the large canine heart) and high voxel resolution (36 µm) provided. In summary, contrast micro-CT facilitates fast and nondestructive imaging and segmenting of detailed 3-D cardiovascular geometries, as well as measuring fiber orientation, which are crucial in constructing biophysically detailed computational cardiac models.

Arrhythmogenic substrate for atrial fibrillation: insights from an integrative computational model of pulmonary veins.

Mechanisms underlying the genesis of re-entrant substrate for atrial fibrillation (AF) in the pulmonary veins (PVs) and left atrium (LA) are not well understood. We develop a biophysically detailed computational model for the PVs and surrounding LA tissue. The model integrates canine PV and LA single cell electrophysiology with the respective 3D tissue geometry and fiber orientation reconstructed from micro-CT data. The model simulations demonstrate that a combination of tissue anisotropy and electrical heterogeneity between the PVs and LA causes a break-down of normal electrical excitation wave-fronts. This leads to the generation of a high-frequency re-entrant source near the PV sleeves. Evidence of such sources have been seen clinically in AF patients. In summary, our modeling results provide new insights into the arrhythmogenic mechanisms of re-entrant excitation waves underlying AF.