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

Improved algorithm for generating evenly-spaced streamlines from an orientation field on a triangulated surface

BACKGROUND

Vector fields such as cardiac fiber orientation can be visualized on a surface using streamlines. The application of evenly-spaced streamline generation to the construction of interconnected cable structure for cardiac propagation models has more stringent requirements imperfectly fulfilled by current algorithms.

METHOD

We developed an open-source C++/python package for the placement of evenly-spaced streamlines on a triangulated surface. The new algorithm improves upon previous works by more accurately handling streamline extremities, U-turns and limit cycles, by providing stronger geometrical guarantees on inter-streamline minimal distance, particularly when a high density of streamlines (up to 10μm spacing) is desired, and by making a more efficient parallel implementation available. The approach requires finding intersections between geometrical capsules and triangles to update an occupancy mask defined on the triangles. This enables fast streamline integration from thousands of seed points to identify optimal streamline placement.

RESULTS

The algorithm was assessed qualitatively on different left atrial models of fiber orientation with varying mesh resolutions (up to 375k triangles) and quantitatively by measuring streamline lengths and distribution of inter-streamline minimal distance. The complexity and the computational performance of the algorithm were studied as a function of streamline spacing in relation to triangular mesh resolution.

CONCLUSION

More accurate geometrical computations, attention to details and fine-tuning led to an algorithm more amenable to applications that require precise positioning of streamlines.

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 fibrillation: Insights from animal models, computational modeling, and clinical studies

Atrial fibrillation (AF) is the most common human arrhythmia, affecting millions of patients worldwide. A combination of risk factors and comorbidities results in complex atrial remodeling, which increases AF vulnerability and persistence. Insights from animal models, clinical studies, and computational modeling have advanced the understanding of the mechanisms and pathophysiology of AF. Areas of heterogeneous pathological remodeling, as well as altered electrophysiological properties, serve as a substrate for AF drivers and spontaneous activations. The complex and individualized presentation of this arrhythmia suggests that mechanisms-based personalized approaches will likely be needed to overcome current challenges in AF management. In this paper, we review the insights on the mechanisms of AF obtained from animal models and clinical studies and how computational models integrate this knowledge to advance AF clinical management. We also assess the challenges that need to be overcome to implement these mechanistic models in clinical practice.

Arrhythmogenic influence of mutations in a myocyte-based computational model of the pulmonary vein sleeve

In the heart, electrophysiological dysregulation arises from defects at many biological levels (from point mutations in ion channel proteins to gross structural abnormalities). These defects disrupt the normal pattern of electrical activation, producing ectopic activity and reentrant arrhythmia. To interrogate mechanisms that link these primary biological defects to macroscopic electrophysiologic dysregulation most prior computational studies have utilized either (i) detailed models of myocyte ion channel dynamics at limited spatial scales, or (ii) homogenized models of action potential conduction that reproduce arrhythmic activity at tissue and organ levels. Here we apply our recent model (EMI), which integrates electrical activation and propagation across these scales, to study human atrial arrhythmias originating in the pulmonary vein (PV) sleeves. These small structures initiate most supraventricular arrhythmias and include pronounced myocyte-to-myocyte heterogeneities in ion channel expression and intercellular coupling. To test EMI's cell-based architecture in this physiological context we asked whether ion channel mutations known to underlie atrial fibrillation are capable of initiating arrhythmogenic behavior via increased excitability or reentry in a schematic PV sleeve geometry. Our results illustrate that EMI's improved spatial resolution can directly interrogate how electrophysiological changes at the individual myocyte level manifest in tissue and as arrhythmia in the PV sleeve.

Understanding PITX2-Dependent Atrial Fibrillation Mechanisms through Computational Models

Atrial fibrillation (AF) is a common arrhythmia. Better prevention and treatment of AF are needed to reduce AF-associated morbidity and mortality. Several major mechanisms cause AF in patients, including genetic predispositions to AF development. Genome-wide association studies have identified a number of genetic variants in association with AF populations, with the strongest hits clustering on chromosome 4q25, close to the gene for the homeobox transcription PITX2. Because of the inherent complexity of the human heart, experimental and basic research is insufficient for understanding the functional impacts of PITX2 variants on AF. Linking PITX2 properties to ion channels, cells, tissues, atriums and the whole heart, computational models provide a supplementary tool for achieving a quantitative understanding of the functional role of PITX2 in remodelling atrial structure and function to predispose to AF. It is hoped that computational approaches incorporating all we know about PITX2-related structural and electrical remodelling would provide better understanding into its proarrhythmic effects leading to development of improved anti-AF therapies. In the present review, we discuss advances in atrial modelling and focus on the mechanistic links between PITX2 and AF. Challenges in applying models for improving patient health are described, as well as a summary of future perspectives.

Neural Tissue Degeneration in Rosenthal's Canal and Its Impact on Electrical Stimulation of the Auditory Nerve by Cochlear Implants: An Image-Based Modeling Study

Sensorineural deafness is caused by the loss of peripheral neural input to the auditory nerve, which may result from peripheral neural degeneration and/or a loss of inner hair cells. Provided spiral ganglion cells and their central processes are patent, cochlear implants can be used to electrically stimulate the auditory nerve to facilitate hearing in the deaf or severely hard-of-hearing. Neural degeneration is a crucial impediment to the functional success of a cochlear implant. The present, first-of-its-kind two-dimensional finite-element model investigates how the depletion of neural tissues might alter the electrically induced transmembrane potential of spiral ganglion neurons. The study suggests that even as little as 10% of neural tissue degeneration could lead to a disproportionate change in the stimulation profile of the auditory nerve. This result implies that apart from encapsulation layer formation around the cochlear implant electrode, tissue degeneration could also be an essential reason for the apparent inconsistencies in the functionality of cochlear implants.

Investigation of optimal sample preparation conditions with potassium triiodide and optimal imaging settings for microfocus computed tomography of excised cat hearts

OBJECTIVE

To determine optimal sample preparation conditions with potassium triiodide (I₂KI) and optimal imaging settings for microfocus CT (micro-CT) of excised cat hearts.

SAMPLE

7 excised hearts (weight range, 10 to 17.6 g) obtained from healthy adult cats after euthanasia by IV injection of pentobarbital sodium.

PROCEDURES

Following excision, the hearts were preserved in 10% formaldehyde solution. Six hearts were immersed in 1.25% I₂KI solution (n = 3) or 2.5% I₂KI solution (3) for a 12-day period. Micro-CT images were acquired at time 0 (prior to iodination) then approximately every 24 and 48 hours thereafter to determine optimal sample preparation conditions (ie, immersion time and concentration of I₂KI solution). Identified optimal conditions were then used to prepare the seventh heart for imaging; changes in voltage, current, exposure time, and gain on image quality were evaluated to determine optimal settings (ie, maximal signal-to-noise and contrast-to-noise ratios). Images were obtained at a voxel resolution of 30 μm. A detailed morphological assessment of the main cardiac structures of the seventh heart was then performed.

RESULTS

Immersion in 2.5% I₂KI solution for 48 hours was optimal for sample preparation. The optimal imaging conditions included a tube voltage of 100 kV, current of 150 μA, and exposure time of 354 milliseconds; scan duration was 12 minutes.

CONCLUSIONS AND CLINICAL RELEVANCE

Results provided an optimal micro-CT imaging protocol for excised cat hearts prepared with I₂KI solution that could serve as a basis for future studies of micro-CT for high resolution 3-D imaging of cat hearts.

Monitoring the dynamics of acute radiofrequency ablation lesion formation in thin-walled atria - a simultaneous optical and electrical mapping study

Background Radiofrequency ablation (RFA) is a common approach to treat cardiac arrhythmias. During this intervention, numerous strategies are applied to indirectly estimate lesion formation. However, the assessment of the spatial extent of these acute injuries needs to be improved in order to create well-defined and durable ablation lesions. Methods We investigated the electrophysiological characteristics of rat atrial myocardium during an ex vivo RFA procedure with fluorescence-optical and electrical mapping. By analyzing optical data, the temporal growth of punctiform ablation lesions was reconstructed after stepwise RFA sequences. Unipolar electrograms (EGMs) were simultaneously recorded by a multielectrode array (MEA) before and after each RFA sequence. Based on the optical results, we searched for electrical features to delineate these lesions from healthy myocardium. Results Several unipolar EGM parameters were monotonically decreasing when distances between the electrode and lesion boundary were smaller than 2 mm. The negative component of the unipolar EGM [negative peak amplitude (Aneg)] vanished for distances lesser than 0.4 mm to the lesion boundary. Median peak-to-peak amplitude (Vpp) was decreased by 75% compared to baseline. Conclusion Aneg and Vpp are excellent parameters to discriminate the growing lesion area from healthy myocardium. The experimental setup opens new opportunities to investigate EGM characteristics of more complex ablation lesions.

Challenges Associated with Interpreting Mechanisms of AF

Determining optimal treatment strategies for complex arrhythmogenesis in AF is confounded by the lack of consensus regarding the mechanisms causing AF. Studies report different mechanisms for AF, ranging from hierarchical drivers to anarchical multiple activation wavelets. Differences in the assessment of AF mechanisms are likely due to AF being recorded across diverse models using different investigational tools, spatial scales and clinical populations. The authors review different AF mechanisms, including anatomical and functional re-entry, hierarchical drivers and anarchical multiple wavelets. They then describe different cardiac mapping techniques and analysis tools, including activation mapping, phase mapping and fibrosis identification. They explain and review different data challenges, including differences between recording devices in spatial and temporal resolutions, spatial coverage and recording surface, and report clinical outcomes using different data modalities. They suggest future research directions for investigating the mechanisms underlying human AF.

Patient-specific simulations predict efficacy of ablation of interatrial connections for treatment of persistent atrial fibrillation

Aims

Treatments for persistent atrial fibrillation (AF) offer limited efficacy. One potential strategy aims to return the right atrium (RA) to sinus rhythm (SR) by ablating interatrial connections (IAC) to isolate the atria, but there is limited clinical data to evaluate this ablation approach. We aimed to use simulation to evaluate and predict patient-specific suitability for ablation of IAC to treat AF.

Methods and results

Persistent AF was simulated in 12 patient-specific geometries, incorporating electrophysiological heterogeneity and fibres, with IAC at Bachmann’s bundle, the coronary sinus, and fossa ovalis. Simulations were performed to test the effect of left atrial (LA)-to-RA frequency gradient and fibrotic remodelling on IAC ablation efficacy. During AF, we simulated ablation of one, two, or all three IAC, with or without pulmonary vein isolation and determined if this altered or terminated the arrhythmia. For models without structural remodelling, ablating all IAC terminated RA arrhythmia in 83% of cases. Models with the LA-to-RA frequency gradient removed had an increased success rate (100% success). Ablation of IACs is less effective in cases with fibrotic remodelling (interstitial fibrosis 50% success rate; combination remodelling 67%). Mean number of phase singularities in the RA was higher pre-ablation for IAC failure (success 0.6 ± 0.8 vs. failure 3.2 ± 2.5, P < 0.001).

Conclusion

This simulation study predicts that IAC ablation is effective in returning the RA to SR for many cases. Patient-specific modelling approaches have the potential to stratify patients prior to ablation by predicting if drivers are located in the LA or RA. We present a platform for predicting efficacy and informing patient selection for speculative treatments.

Understanding the Beat-to-Beat Variations of P-Waves Morphologies in AF Patients During Sinus Rhythm: A Scoping Review of the Atrial Simulation Studies

The remarkable advances in high-performance computing and the resulting increase of the computational power have the potential to leverage computational cardiology toward improving our understanding of the pathophysiological mechanisms of arrhythmias, such as Atrial Fibrillation (AF). In AF, a complex interaction between various triggers and the atrial substrate is considered to be the leading cause of AF initiation and perpetuation. In electrocardiography (ECG), P-wave is supposed to reflect atrial depolarization. It has been found that even during sinus rhythm (SR), multiple P-wave morphologies are present in AF patients with a history of AF, suggesting a higher dispersion of the conduction route in this population. In this scoping review, we focused on the mechanisms which modify the electrical substrate of the atria in AF patients, while investigating the existence of computational models that simulate the propagation of the electrical signal through different routes. The adopted review methodology is based on a structured analytical framework which includes the extraction of the keywords based on an initial limited bibliographic search, the extensive literature search and finally the identification of relevant articles based on the reference list of the studies. The leading mechanisms identified were classified according to their scale, spanning from mechanisms in the cell, tissue or organ level, and the produced outputs. The computational modeling approaches for each of the factors that influence the initiation and the perpetuation of AF are presented here to provide a clear overview of the existing literature. Several levels of categorization were adopted while the studies which aim to translate their findings to ECG phenotyping are highlighted. The results denote the availability of multiple models, which are appropriate under specific conditions. However, the consideration of complex scenarios taking into account multiple spatiotemporal scales, personalization of electrophysiological and anatomical models and the reproducibility in terms of ECG phenotyping has only partially been tackled so far.

Tissue-Specific Optical Mapping Models of Swine Atria Informed by Optical Coherence Tomography

Computational models and experimental optical mapping of cardiac electrophysiology serve as powerful tools to investigate the underlying mechanisms of arrhythmias. Modeling can also aid the interpretation of optical mapping signals, which may have different characteristics with respect to the underlying electrophysiological signals they represent. However, despite the prevalence of atrial arrhythmias such as atrial fibrillation, models of optical electrical mapping incorporating realistic structure of the atria are lacking. Therefore, we developed image-based models of atrial tissue using structural information extracted from optical coherence tomography (OCT), which can provide volumetric tissue characteristics in high resolution. OCT volumetric data of four swine atrial tissue samples were used to develop models incorporating tissue geometry, tissue-specific myofiber orientation, and ablation lesion regions. We demonstrated the use of these models through electrophysiology and photon scattering simulations. Changes in transmural electrical conduction were observed with the inclusion of OCT-derived, depth-resolved fiber orientation. Additionally, the amplitude of optical mapping signals were not found to correspond with lesion transmurality because of lesion geometry and electrical propagation occurring beyond excitation light penetration. This work established a framework for the development of tissue-specific models of atrial tissue derived from OCT imaging data, which can be useful in future investigations of electrophysiology and optical mapping signals with respect to realistic atrial tissue structure.

A technique for measuring anisotropy in atrial conduction to estimate conduction velocity and atrial fibre direction

BACKGROUND

Cardiac conduction properties exhibit large variability, and affect patient-specific arrhythmia mechanisms. However, it is challenging to clinically measure conduction velocity (CV), anisotropy and fibre direction. Our aim is to develop a technique to estimate conduction anisotropy and fibre direction from clinically available electrical recordings.

METHODS

We developed and validated automated algorithms for estimating cardiac CV anisotropy, from any distribution of recording locations on the atrial surface. The first algorithm is for elliptical wavefront fitting to a single activation map (method 1), which works well close to the pacing location, but decreases in accuracy further from the pacing location (due to spatial heterogeneity in the conductivity and fibre fields). As such, we developed a second methodology for measuring local conduction anisotropy, using data from two or three activation maps (method 2: ellipse fitting to wavefront propagation velocity vectors from multiple activation maps).

RESULTS

Ellipse fitting to CV vectors from two activation maps (method 2) leads to an improved estimation of longitudinal and transverse CV compared to method 1, but fibre direction estimation is still relatively poor. Using three activation maps with method 2 provides accurate estimation, with approximately 70% of atrial fibres estimated within 20∘. We applied the technique to clinical activation maps to demonstrate the presence of heterogeneous conduction anisotropy, and then tested the effects of this conduction anisotropy on predicted arrhythmia dynamics using computational simulation.

CONCLUSIONS

We have developed novel algorithms for calculating CV and measuring the direction dependency of atrial activation to estimate atrial fibre direction, without the need for specialised pacing protocols, using clinically available electrical recordings.

Variability in pulmonary vein electrophysiology and fibrosis determines arrhythmia susceptibility and dynamics

Success rates for catheter ablation of persistent atrial fibrillation patients are currently low; however, there is a subset of patients for whom electrical isolation of the pulmonary veins alone is a successful treatment strategy. It is difficult to identify these patients because there are a multitude of factors affecting arrhythmia susceptibility and maintenance, and the individual contributions of these factors are difficult to determine clinically. We hypothesised that the combination of pulmonary vein (PV) electrophysiology and atrial body fibrosis determine driver location and effectiveness of pulmonary vein isolation (PVI). We used bilayer biatrial computer models based on patient geometries to investigate the effects of PV properties and atrial fibrosis on arrhythmia inducibility, maintenance mechanisms, and the outcome of PVI. Short PV action potential duration (APD) increased arrhythmia susceptibility, while longer PV APD was found to be protective. Arrhythmia inducibility increased with slower conduction velocity (CV) at the LA/PV junction, but not for cases with homogeneous CV changes or slower CV at the distal PV. Phase singularity (PS) density in the PV region for cases with PV fibrosis was increased. Arrhythmia dynamics depend on both PV properties and fibrosis distribution, varying from meandering rotors to PV reentry (in cases with baseline or long APD), to stable rotors at regions of high fibrosis density. Measurement of fibrosis and PV properties may indicate patient specific susceptibility to AF initiation and maintenance. PV PS density before PVI was higher for cases in which AF terminated or converted to a macroreentry; thus, high PV PS density may indicate likelihood of PVI success.

Atrial fibrillation is the most commonly encountered cardiac arrhythmia, affecting a significant portion of the population. Currently, ablation is the most effective treatment but success rates are less than optimal, being 70% one-year post-treatment. There is a large effort to find better ablation strategies to permanently cure the condition. Pulmonary vein isolation by ablation is more or less the standard of care, but many questions remain since pulmonary vein ectopy by itself does not explain all of the clinical successes or failures. We used computer simulations to investigate how electrophysiological properties of the pulmonary veins can affect rotor formation and maintenance in patients suffering from atrial fibrillation. We used complex, biophysical representations of cellular electrophysiology in highly detailed geometries constructed from patient scans. We heterogeneously varied electrophysiological and structural properties to see their effects on rotor initiation and maintenance. Our study suggests a metric for indicating the likelihood of success of pulmonary vein isolation. Thus either measuring this clinically, or running patient-specific simulations to estimate this metric may suggest whether ablation in addition to pulmonary vein isolation should be performed. Our study provides motivation for a retrospective clinical study or experimental study into this metric.

Structural heterogeneity of the rat pulmonary vein myocardium: consequences on intracellular calcium dynamics and arrhythmogenic potential

Mechanisms underlying ectopic activity in the pulmonary vein (PV) which triggers paroxysmal atrial fibrillation are unknown. Although several studies have suggested that calcium signalling might be involved in these arrhythmias, little is known about calcium cycling in PV cardiomyocytes (CM). We found that individual PV CM showed a wide range of transverse tubular incidence and organization, going from their virtual absence, as described in atrial CM, to well transversally organised tubular systems, like in ventricular CM. These different types of CM were found in groups scattered throughout the tissue. The variability of the tubular system was associated with cell to cell heterogeneity of calcium channel (Ca_v1.2) localisation and, thereby, of Ca_v1.2-Ryanodine receptor coupling. This was responsible for multiple forms of PV CM calcium transient. Spontaneous calcium sparks and waves were not only more abundant in PV CM than in LA CM but also associated with a higher depolarising current. In conclusion, compared with either the atrium or the ventricle, PV myocardium presents marked structural and functional heterogeneity.

Tensile rupture of medial arterial tissue studied by X-ray micro-tomography on stained samples

Detailed characterization of damage and rupture mechanics of arteries is one the current challenges in vascular biomechanics, which requires developing suitable experimental approaches. This paper introduces an approach using in situ tensile tests in an X-ray micro-tomography setup to observe mechanisms of damage initiation and progression in medial layers of porcine aortic samples. The technique requires the use of sodium polytungstate as a contrast agent, of which the conditions for use are detailed in this paper. Immersion of the samples during 24h in a 15g/L concentrated solution provided the best compromise for viewing musculo-elastic units in this tissue. The process of damage initiation, delamination and rupture of medial tissue under tensile loading was observed and can be described as an elementary process repeating several times until complete failure. This elementary process initiates with a sudden mode I fracture of a group of musculo-elastic units, followed by an elastic recoil of these units, causing mode II separation of these, hence a delamination plane. The presented experimental approach constitutes a basis for observation of other constituents, or for investigations on other tissues and damage mechanisms.

In silico assessment of genetic variation in KCNA5 reveals multiple mechanisms of human atrial arrhythmogenesis

A recent experimental study investigating patients with lone atrial fibrillation identified six novel mutations in the KCNA5 gene. The mutants exhibited both gain- and loss-of-function of the atrial specific ultra-rapid delayed rectifier K⁺ current, I_Kur. The aim of this study is to elucidate and quantify the functional impact of these KCNA5 mutations on atrial electrical activity. A multi-scale model of the human atria was updated to incorporate detailed experimental data on I_Kur from both wild-type and mutants. The effects of the mutations on human atrial action potential and rate dependence were investigated at the cellular level. In tissue, we assessed the effects of the mutations on the vulnerability to unidirectional conduction patterns and dynamics of re-entrant excitation waves. Gain-of-function mutations shortened the action potential duration in single cells, and stabilised and accelerated re-entrant excitation in tissue. Loss-of-function mutations had heterogeneous effects on action potential duration and promoted early-after-depolarisations following beta-adrenergic stimulation. In the tissue model, loss-of-function mutations facilitated breakdown of excitation waves at more physiological excitation rates than the wild-type, and the generation of early-after-depolarisations promoted unidirectional patterns of excitation. Gain- and loss-of-function I_Kur mutations produced multiple mechanisms of atrial arrhythmogenesis, with significant differences between the two groups of mutations. This study provides new insights into understanding the mechanisms by which mutant I_Kur contributes to atrial arrhythmias. In addition, as I_Kur is an atrial-specific channel and a number of I_Kur-selective blockers have been developed as anti-AF agents, this study also helps to understand some contradictory results on both pro- and anti-arrhythmic effects of blocking I_Kur.

In a recent study, six mutations resulting in either gain-of-function or loss-of-function in the ultra-rapid delayed rectifier potassium current I_Kur, were identified to be associated with atrial fibrillation (AF). However, the causative link between the mutant I_Kur (either gain- or loss-of-function) and AF genesis, especially the difference and similarity between the two mutant groups, has not been elucidated. In our study, we used multiscale computational models to investigate the mechanism of arrhythmogenesis mediated by the two groups of mutations. The results suggest that the gain-of-function mutations shortened atrial action potential duration, stabilised and accelerated re-entrant excitation waves in tissue; the loss-of-function mutation promoted early-after-depolarisations following beta-adrenergic stimulation and thus wave breaks in tissue. We show these two groups of mutants carrying I_Kur produced multiple mechanisms of atrial arrhythmogenesis, with significant differences between the two groups. Our study provides new insights into understanding the mechanisms by which mutant I_Kur contributes to atrial arrhythmias. In addition, as I_Kur is an atrial-specific channel and a number of I_Kur-selective blockers have been developed as anti-AF agents, this study also helps to understand some contradictory results on both pro- and anti-arrhythmic effects of blocking I_Kur.

Anti-arrhythmic strategies for atrial fibrillation: The role of computational modeling in discovery, development, and optimization

Atrial fibrillation (AF), the most common cardiac arrhythmia, is associated with increased risk of cerebrovascular stroke, and with several other pathologies, including heart failure. Current therapies for AF are targeted at reducing risk of stroke (anticoagulation) and tachycardia-induced cardiomyopathy (rate or rhythm control). Rate control, typically achieved by atrioventricular nodal blocking drugs, is often insufficient to alleviate symptoms. Rhythm control approaches include antiarrhythmic drugs, electrical cardioversion, and ablation strategies. Here, we offer several examples of how computational modeling can provide a quantitative framework for integrating multiscale data to: (a) gain insight into multiscale mechanisms of AF; (b) identify and test pharmacological and electrical therapy and interventions; and (c) support clinical decisions. We review how modeling approaches have evolved and contributed to the research pipeline and preclinical development and discuss future directions and challenges in the field.

Electrophysiological Characteristics, Rhythm, Disturbances and Conduction Discontinuities Under Autonomic Stimulation in the Rat Pulmonary Vein Myocardium

INTRODUCTION

Despite the importance of neurogenic initiation of rapid firing from pulmonary veins (PVs), the mechanism of autonomic modulation of electrophysiological properties of the PV myocardium to form a substrate for atrial arrhythmia remains poorly understood.

METHODS AND RESULTS

A 2-microelectrode technique was used to characterize electrophysiological properties of rat PV myocardium and to explore PV arrhythmogenesis, at baseline, during electrical stimulation and/or under autonomic modulation. PV myocardium was characterized by prolonged action potential duration (APD), high degree of APD alternans, and spontaneous depolarizations. Autonomic stimulation resulted in significantly enhanced APD dispersion within the PV, which dynamically changed over time and was associated with intra-PV and atria-PV conduction blocks and could lead to spontaneous fibrillation-like high-frequency activity. In the distal part of the PV we found an unexcitable area that was characterized by depolarized resting potential (-50 ± 4 mV vs. -75 ± 2 mV vs. PV mouth, P < 0.01). This region could be activated during autonomic stimulation or fast pacing that led to multiple conduction discontinuities (uni- and bi-directional conduction blocks, Wenckebach periodicity, electrotonic modulation conduction block, echo phenomenon) in 17/23 preparations, including those occurring under norepinephrine superfusion (14/17) and during pacing frequency changes (3/17). PV echoes (unstable reentrant circuits) were found in 8/23 preparations. In some experiments, several types of conduction abnormalities were observed.

CONCLUSION

The PV myocardium demonstrates distinct electrophysiological characteristics, which could be considerably exaggerated by electrical stimulation and/or autonomic nervous system to dynamically form a functional substrate to support re-entry as well as focal activity.

Assessing the three-dimensional collagen network in soft tissues using contrast agents and high resolution micro-CT: Application to porcine iliac veins

The assessment of the three-dimensional architecture of collagen fibers inside vessel walls constitutes one of the bases for building structural models for the description of the mechanical behavior of these tissues. Multiphoton microscopy allows for such observations, but is limited to volumes of around a thousand of microns. In the present work, we propose to observe the collagenous network of vascular tissues using micro-CT. To get a contrast, three staining solutions (phosphotungstic acid, phosphomolybdic acid and iodine potassium iodide) were tested. Two of these stains were showed to lead to similar results and to a satisfactory contrast within the tissue. A detailed observation of a small porcine iliac vein sample allowed assessing the collagen fibers orientations within the medial and adventitial layers of the vein. The vasa vasorum network, which is present inside the adventitia of the vein, was also observed. Finally, the demonstrated micro-CT staining technique for the three-dimensional observation of thin soft tissues samples, like vein walls, contributes to the assessment of their structure at different scales while keeping a global overview of the tissue.

Optimization of catheter ablation of atrial fibrillation: insights gained from clinically-derived computer models

Atrial fibrillation (AF) is the most common heart rhythm disturbance, and its treatment is an increasing economic burden on the health care system. Despite recent intense clinical, experimental and basic research activity, the treatment of AF with current antiarrhythmic drugs and catheter/surgical therapies remains limited. Radiofrequency catheter ablation (RFCA) is widely used to treat patients with AF. Current clinical ablation strategies are largely based on atrial anatomy and/or substrate detected using different approaches, and they vary from one clinical center to another. The nature of clinical ablation leads to ambiguity regarding the optimal patient personalization of the therapy partly due to the fact that each empirical configuration of ablation lines made in a patient is irreversible during one ablation procedure. To investigate optimized ablation lesion line sets, in silico experimentation is an ideal solution. 3D computer models give us a unique advantage to plan and assess the effectiveness of different ablation strategies before and during RFCA. Reliability of in silico assessment is ensured by inclusion of accurate 3D atrial geometry, realistic fiber orientation, accurate fibrosis distribution and cellular kinetics; however, most of this detailed information in the current computer models is extrapolated from animal models and not from the human heart. The predictive power of computer models will increase as they are validated with human experimental and clinical data. To make the most from a computer model, one needs to develop 3D computer models based on the same functionally and structurally mapped intact human atria with high spatial resolution. The purpose of this review paper is to summarize recent developments in clinically-derived computer models and the clinical insights they provide for catheter ablation.

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.

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 (I_CaL, I_K1, I_to, and I_Kur), 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 (I_K1 and I_Kur 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.

Optimal iodine staining of cardiac tissue for X-ray computed tomography

X-ray computed tomography (XCT) has been shown to be an effective imaging technique for a variety of materials. Due to the relatively low differential attenuation of X-rays in biological tissue, a high density contrast agent is often required to obtain optimal contrast. The contrast agent, iodine potassium iodide (), has been used in several biological studies to augment the use of XCT scanning. Recently was used in XCT scans of animal hearts to study cardiac structure and to generate 3D anatomical computer models. However, to date there has been no thorough study into the optimal use of as a contrast agent in cardiac muscle with respect to the staining times required, which has been shown to impact significantly upon the quality of results. In this study we address this issue by systematically scanning samples at various stages of the staining process. To achieve this, mouse hearts were stained for up to 58 hours and scanned at regular intervals of 6–7 hours throughout this process. Optimal staining was found to depend upon the thickness of the tissue; a simple empirical exponential relationship was derived to allow calculation of the required staining time for cardiac samples of an arbitrary size.

A mathematical model of the unidirectional block caused by the pulmonary veins for anatomically induced atrial reentry

It is widely believed that the pulmonary veins (PVs) of the left atrium play the central role in the generation of anatomically induced atrial reentry but its mechanism has not been analytically explained. To understand this mechanism, a new analytic approach is proposed by adapting the geometric relative acceleration analysis from spacetime physics based on the hypothesis that a large relative acceleration can translate to a dramatic increase in the curvature of a wavefront and subsequently to conduction failure. By verifying the strong dependency of the propagational direction and the magnitude of anisotropy for conduction failure, this analytic method reveals that a unidirectional block can be generated by asymmetric propagation toward the PVs. This model is validated by computational tests in a T-shaped domain, computational simulations for three-dimensional atrial reentry and previous in-silico reports for anatomically induced atrial reentry.

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 (I_CaL, I_Kur, I_to, I_K1, I_Ks, I_NaCa) and on intracellular Ca²⁺ 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.