E. Hofer

Automatically Generated, Anatomically Accurate Meshes for Cardiac Electrophysiology Problems

Significant advancements in imaging technology and the dramatic increase in computer power over the last few years broke the ground for the construction of anatomically realistic models of the heart at an unprecedented level of detail. To effectively make use of high-resolution imaging datasets for modeling purposes, the imaged objects have to be discretized. This procedure is trivial for structured grids. However, to develop generally applicable heart models, unstructured grids are much preferable. In this study, a novel image-based unstructured mesh generation technique is proposed. It uses the dual mesh of an octree applied directly to segmented 3-D image stacks. The method produces conformal, boundary-fitted, and hexahedra-dominant meshes. The algorithm operates fully automatically with no requirements for interactivity and generates accurate volume-preserving representations of arbitrarily complex geometries with smooth surfaces. The method is very well suited for cardiac electrophysiological simulations. In the myocardium, the algorithm minimizes variations in element size, whereas in the surrounding medium, the element size is grown larger with the distance to the myocardial surfaces to reduce the computational burden. The numerical feasibility of the approach is demonstrated by discretizing and solving the monodomain and bidomain equations on the generated grids for two preparations of high experimental relevance, a left ventricular wedge preparation, and a papillary muscle.

Electroanatomical Characterization of Atrial Microfibrosis in a Histologically Detailed Computer Model

Fibrosis is thought to play an important role in the formation and maintenance of atrial fibrillation (AF). The propensity of fibrosis to increase AF vulnerability depends not only on its amount, its texture plays a crucial role as well. While the detection of fibrotic tissue patches in the atria with extracellular recordings is feasible based on the analysis of electrogram fractionation, as used in clinical practice to identify ablation targets, the classification of fibrotic texture is a more challenging problem. This study seeks to establish a method for the electroanatomical characterization of the fibrotic textures based on the analysis of electrogram fractionation. The proposed method exploits the dependence of fractionation patterns on the incidence direction of wavefronts which differs significantly as a function of texture. A histologically detailed computer model of the right atrial isthmus was developed for testing the method. A stimulation protocol was conceived which generated various incidence directions for any given recording site where electrograms were computed. A classification method is derived then for discriminating three types of fibrosis, no fibrosis (control), diffuse, and patchy fibrosis. Simulation results showed that electrogram fractionation and amplitudes and their dependence upon incidence direction allow a robust discrimination between different classes of fibrosis. Finally, to minimize the technical effort, sensitivity analysis was performed to identify a minimum number of incidence directions required for robust classification.

An Efficient Finite Element Approach for Modeling Fibrotic Clefts in the Heart

Advanced medical imaging technologies provide a wealth of information on cardiac anatomy and structure at a paracellular resolution, allowing to identify microstructural discontinuities which disrupt the intracellular matrix. Current state-of-the-art computer models built upon such datasets account for increasingly finer anatomical details, however, structural discontinuities at the paracellular level are typically discarded in the model generation process, owing to the significant costs which incur when using high resolutions for explicit representation. In this study, a novel discontinuous finite element (dFE) approach for discretizing the bidomain equations is presented, which accounts for fine-scale structures in a computer model without the need to increase spatial resolution. In the dFE method, this is achieved by imposing infinitely thin lines of electrical insulation along edges of finite elements which approximate the geometry of discontinuities in the intracellular matrix. Simulation results demonstrate that the dFE approach accounts for effects induced by microscopic size scale discontinuities, such as the formation of microscopic virtual electrodes, with vast computational savings as compared to high resolution continuous finite element models. Moreover, the method can be implemented in any standard continuous finite element code with minor effort.

Contributions of Purkinje-myocardial coupling to suppression and facilitation of early afterdepolarization-induced triggered activity

Electrical loading by ventricular myocardium modulates conduction system repolarization near Purkinje-ventricular junctions (PVJs). We investigated how that loading suppresses and facilitates early afterdepolarizations (EADs) under conditions where there is a high degree of functional coupling between tissue types, which is consistent with the anatomic arrangement at the peripheral conduction system-myocardial interface. Experiments were completed in eight rabbit right ventricular (RV) free wall preparations. Free-running Purkinje strands were locally superfused, and action potentials were recorded from strands. RV free walls were bathed in normal solution. Surface electrograms were recorded near strand insertions into downstream free wall myocardium. Detailed histology was performed to assemble a computer model with interspersed Purkinje and ventricular myocytes weakly coupled throughout the region. Delays from Purkinje upstrokes to downstream peripheral conduction system and myocardial activation were comparable between experiments and simulations, supporting model node-to-node electrical coupling, i.e., the functional coupling. Purkinje action potential duration (APD) prolongation with localized isoproterenol in experiments and calcium current enhancement in simulations failed to establish EADs. With myocardial APD prolongation by delayed rectifier potassium current inhibition or L-type calcium current enhancement accompanying Purkinje APD prolongation in simulations, however, EAD-induced triggered activity developed. Collectively, our findings suggest competing contributions of the myocardial sink when there is a high degree of functional coupling between tissue types, with the transition from suppression to facilitation of EAD-induced triggered activity depending critically upon myocardial APD prolongation.

Use of Cardiac Electric Near-Field Measurements to Determine Activation Times

AbstractIn a recent paper, we described the behavior of the cardiac electric near-field, E, parallel to the tissue surface during continuous conduction. We found that TE, the time at which the peak near-field,

Decomposition of fractionated local electrograms using an analytic signal model based on sigmoid functions.

\hat E

Oblique Propagation of Activation Allows the Detection of Uncoupling Microstructures from Cardiac Near Field Behavior

, occurs, is an accurate marker of local activation time. Examination of experimentally recorded E vector loops revealed a large variety of morphologies. We postulated that propagation around an obstacle could lead to the observed deviations in loop morphology. The purpose of this study was to determine if this was plausible, and if so, whether TE remains an accurate time marker of local activation under these conditions. We used a monodomain computer model of a sheet of cardiac tissue with a central conduction obstacle immersed in an unbounded volume conductor. Activation times