Slow conduction in cardiac tissue, II: effects of branching tissue geometry

Cardiac arrhythmogenesis: roles of ion channels and their functional modification

Cardiac arrhythmias cause significant morbidity and mortality and pose a major public health problem. They arise from disruptions in the normally orderly propagation of cardiac electrophysiological activation and recovery through successive cardiomyocytes in the heart. They reflect abnormalities in automaticity, initiation, conduction, or recovery in cardiomyocyte excitation. The latter properties are dependent on surface membrane electrophysiological mechanisms underlying the cardiac action potential. Their disruption results from spatial or temporal instabilities and heterogeneities in the generation and propagation of cellular excitation. These arise from abnormal function in their underlying surface membrane, ion channels, and transporters, as well as the interactions between them. The latter, in turn, form common regulatory targets for the hierarchical network of diverse signaling mechanisms reviewed here. In addition to direct molecular-level pharmacological or physiological actions on these surface membrane biomolecules, accessory, adhesion, signal transduction, and cytoskeletal anchoring proteins modify both their properties and localization. At the cellular level of excitation-contraction coupling processes, Ca2+ homeostatic and phosphorylation processes affect channel activity and membrane excitability directly or through intermediate signaling. Systems-level autonomic cellular signaling exerts both acute channel and longer-term actions on channel expression. Further upstream intermediaries from metabolic changes modulate the channels both themselves and through modifying Ca2+ homeostasis. Finally, longer-term organ-level inflammatory and structural changes, such as fibrotic and hypertrophic remodeling, similarly can influence all these physiological processes with potential pro-arrhythmic consequences. These normal physiological processes may target either individual or groups of ionic channel species and alter with particular pathological conditions. They are also potentially alterable by direct pharmacological action, or effects on longer-term targets modifying protein or cofactor structure, expression, or localization. Their participating specific biomolecules, often clarified in experimental genetically modified models, thus constitute potential therapeutic targets. The insights clarified by the physiological and pharmacological framework outlined here provide a basis for a recent modernized drug classification. Together, they offer a translational framework for current drug understanding. This would facilitate future mechanistically directed therapeutic advances, for which a number of examples are considered here. The latter are potentially useful for treating cardiac, in particular arrhythmic, disease.

Non-Invasive Electroanatomical Mapping: A State-Space Approach for Myocardial Current Density Estimation

Electroanatomical mapping is a method for creating a model of the electrophysiology of the human heart. Medical professionals routinely locate and ablate the site of origin of cardiac arrhythmias with invasive catheterization. Non-invasive localization takes the form of electrocardiographic (ECG) or magnetocardiographic (MCG) imaging, where the goal is to reconstruct the electrical activity of the human heart. Non-invasive alternatives to catheter electroanatomical mapping would reduce patients' risks and open new venues for treatment planning and prevention. This work introduces a new system state-based method for estimating the electrical activity of the human heart from MCG measurements. Our model enables arbitrary propagation paths and velocities. A Kalman filter optimally estimates the current densities under the given measurements and model parameters. In an outer optimization loop, these model parameters are then optimized via gradient descent. This paper aims to establish the foundation for future research by providing a detailed mathematical explanation of the algorithm. We demonstrate the feasibility of our method through a simplified one-layer simulation. Our results show that the algorithm can learn the propagation paths from the magnetic measurements. A threshold-based segmentation into healthy and pathological tissue yields a DICE score of 0.84, a recall of 0.77, and a precision of 0.93.

Lipomatous Metaplasia Facilitates Slow Conduction in Critical Ventricular Tachycardia Corridors Within Postinfarct Myocardium

BACKGROUND

Myocardial lipomatous metaplasia (LM) has been reported to be associated with post-infarct ventricular tachycardia (VT) circuitry.

OBJECTIVES

This study examined the association of scar versus LM composition with impulse conduction velocity (CV) in putative VT corridors that traverse the infarct zone in post-infarct patients.

METHODS

The cohort included 31 post-infarct patients from the prospective INFINITY (Intra-Myocardial Fat Deposition and Ventricular Tachycardia in Cardiomyopathy) study. Myocardial scar, border zone, and potential viable corridors were defined by late gadolinium enhancement cardiac magnetic resonance (LGE-CMR), and LM was defined by computed tomography. Images were registered to electroanatomic maps, and the CV at each electroanatomic map point was calculated as the mean CV between that point and 5 adjacent points along the activation wave front.

RESULTS

Regions with LM exhibited lower CV than scar (median = 11.9 vs 13.5 cm/s; P < 0.001). Of 94 corridors computed from LGE-CMR and electrophysiologically confirmed to participate in VT circuitry, 93 traversed through or near LM. These critical corridors displayed slower CV (median 8.8 [IQR: 5.9-15.7] cm/s vs 39.2 [IQR: 28.1-58.5]) cm/s; P < 0.001) than 115 noncritical corridors distant from LM. Additionally, critical corridors demonstrated low-peripheral, high-center (mountain shaped, 23.3%) or mean low-level (46.7%) CV patterns compared with 115 noncritical corridors distant from LM that displayed high-peripheral, low-center (valley shaped, 19.1%) or mean high-level (60.9%) CV patterns.

CONCLUSIONS

The association of myocardial LM with VT circuitry is at least partially mediated by slowing nearby corridor CV thus facilitating an excitable gap that enables circuit re-entry.

Enhanced optimization-based method for the generation of patient-specific models of Purkinje networks

Cardiac Purkinje networks are a fundamental part of the conduction system and are known to initiate a variety of cardiac arrhythmias. However, patient-specific modeling of Purkinje networks remains a challenge due to their high morphological complexity. This work presents a novel method based on optimization principles for the generation of Purkinje networks that combines geometric and activation accuracy in branch size, bifurcation angles, and Purkinje-ventricular-junction activation times. Three biventricular meshes with increasing levels of complexity are used to evaluate the performance of our approach. Purkinje-tissue coupled monodomain simulations are executed to evaluate the generated networks in a realistic scenario using the most recent Purkinje/ventricular human cellular models and physiological values for the Purkinje-ventricular-junction characteristic delay. The results demonstrate that the new method can generate patient-specific Purkinje networks with controlled morphological metrics and specified local activation times at the Purkinje-ventricular junctions.

Intracardiac Electrogram Targets for Ventricular Tachycardia Ablation

The pathogenesis of ventricular tachycardia (VT) in most patients with a prior myocardial scarring is reentry involving compartmentalized muscle fibers protected within the scar. Often the 12-lead ECG morphology of the VT itself is not available when treated with a defibrillator. Consequently, VT ablation takes on an interesting challenge of finding critical targets in sinus rhythm. High-density recordings are essential to evaluate a substrate based on whole electrogram voltage and activation delay, supplemented with substrate perturbation through alternate site pacing or introducing an extra stimulation. In this article, we discuss contemporary intracardiac electrogram targets for VT ablation, with explanation on each of their specific fundamental physiology.

Genome size drives morphological evolution in organ-specific ways

Morphogenesis is an emergent property of biochemical and cellular interactions during development. Genome size and the correlated trait of cell size can influence these interactions through effects on developmental rate and tissue geometry, ultimately driving the evolution of morphology. We tested whether variation in genome and body size is related to morphological variation in the heart and liver using nine species of the salamander genus Plethodon (genome sizes 29-67 gigabases). Our results show that overall organ size is a function of body size, whereas tissue structure changes dramatically with evolutionary increases in genome size. In the heart, increased genome size is correlated with a reduction of myocardia in the ventricle, yielding proportionally less force-producing mass and greater intertrabecular space. In the liver, increased genome size is correlated with fewer and larger vascular structures, positioning hepatocytes farther from the circulatory vessels that transport key metabolites. Although these structural changes should have obvious impacts on organ function, their effects on organismal performance and fitness may be negligible because low metabolic rates in salamanders relax selective pressure on function of key metabolic organs. Overall, this study suggests large genome and cell size influence the developmental systems involved in heart and liver morphogenesis.

Diversity and complexity of the cavotricuspid isthmus in rabbits: A novel scheme for classification and geometrical transformation of anatomical structures

The aim of this study was to describe the morphology of the cavotricuspid isthmus (CTI) in detail and introduce a comprehensive scheme to describe the topology of this region based on functional considerations. This may lead to a better understanding of isthmus-dependent flutter and fibrillation and to improved intervention strategies. We used images of the cavotricuspid isthmus from 52 rabbits of both sexes with a median weight of 3.40 ± 0.93 kg. The area of the CTI was 124.25 ± 42.14 mm² with 53.28 ± 21.13 mm² covered by pectinate muscles connecting the terminal crest and the vestibule. Isthmus length decreased from inferolateral (13.09 ±2.14 mm) to central (9.85 ± 2.14 mm) to paraseptal (4.88 ± 1.96 mm) resembling the overall human geometry. Ramification sites of pectinate muscles were identified and six levels dividing the CTI from posterior to anterior were introduced. This allowed the classification of pectinate muscle segments based on the connected ramification level. To account for the high inter-individual variations in size and shape, the CTI was projected onto a normalized reference frame using bilinear transformation. Furthermore, two measures of complexity were introduced: (i) the ramification index, which reflects the total number of muscle segments connected to a ramification site and (ii) the complexity index, which reflects the type of ramification (branching or merging site). Topological analysis showed that the complexity of the pectinate muscle network decreases from inferolateral to paraseptal and that the number of electrically uncoupled parallel pathways increases in the central section between the terminal crest and the vestibule which introduces potential reentry pathways.

Limitations and Pitfalls of Substrate Mapping for Ventricular Tachycardia

The fundamental hypothesis of substrate mapping for scar-mediated ventricular tachycardia is that surrogates of the isthmus can be identified and targeted with ablation during sinus rhythm. These surrogates include electrocardiographic indications for electric discontinuity such as fractionation, split, late, and long potentials, also evident as sites displaying activation slowing. However, ablation strategies targeting these surrogates during sinus rhythm have resulted in unacceptably high rates of clinical failures, promoting the idea that a more widespread ablation may be required. High-resolution mapping technologies provide an opportunity to examine the substrate at greater detail; however, their use has not yet translated into improved clinical outcomes. This may be related to ongoing efforts to examine the same surrogates at higher resolution instead of using high-resolution technologies for discovering new and potentially more specific surrogates. This article reviews common limitations and pitfalls of substrate mapping and discusses new opportunities for high-resolution mapping to increase the accuracy of substrate mapping: 1) multielectrode mapping catheters provide an opportunity to rapidly examine the substrate during electrophysiological conditions that more closely simulate ventricular tachycardia by means of activation from different directions and coupling intervals; 2) electrogram annotation methods based on the maximal negative derivative of the extracellular potential or maximal voltage are often inaccurate in nonuniform anisotropic tissue. The use of multielectrode catheters may improve the accuracy of electrogram annotation by using spatiotemporal dispersion of single-beat acquisitions and a localized indifferent reference; and 3) resetting and entrainment remain important methods for studying re-entry for and guiding ablation.

Simulation of diffuse and stringy fibrosis in a bilayer interconnected cable model of the left atrium

AIMS

The aim of this study is to design a computer model of the left atrium for investigating fibre-orientation-dependent microstructure such as stringy fibrosis.

METHODS AND RESULTS

We developed an approach for automatic construction of bilayer interconnected cable models from left atrial geometry and epi- and endocardial fibre orientation. The model consisted of two layers (epi- and endocardium) of longitudinal and transverse cables intertwined-like fabric threads, with a spatial discretization of 100 µm. Model validation was performed by comparison with cubic volumetric models in normal conditions. Then, diffuse (n = 2904), stringy (n = 3600), and mixed fibrosis patterns (n = 6840) were randomly generated by uncoupling longitudinal and transverse connections in the interconnected cable model. Fibrosis density was varied from 0% to 40% and mean stringy obstacle length from 0.1 to 2 mm. Total activation time, apparent anisotropy ratio, and local activation time jitter were computed during normal rhythm in each pattern. Non-linear regression formulas were identified for expressing measured propagation parameters as a function of fibrosis density and obstacle length (stringy and mixed patterns). Longer obstacles (even below tissue space constant) were independently associated with prolonged activation times, increased anisotropy, and local fluctuations in activation times. This effect was increased by endo-epicardial dissociation and mitigated when fibrosis was limited to the epicardium.

CONCLUSION

Interconnected cable models enable the study of microstructure in organ-size models despite limitations in the description of transmural structures.

Proarrhythmic Electrical Remodeling by Noncardiomyocytes at Interfaces With Cardiomyocytes Under Oxidative Stress

Life-threatening ventricular arrhythmias, typically arising from interfaces between fibrosis and surviving cardiomyocytes, are feared sequelae of structurally remodeled hearts under oxidative stress. Incomplete understanding of the proarrhythmic electrical remodeling by fibrosis limits the development of novel antiarrhythmic strategies. To define the mechanistic determinants of the proarrhythmia in electrical crosstalk between cardiomyocytes and noncardiomyocytes, we developed a novel in vitro model of interface between neonatal rat ventricular cardiomyocytes (NRVMs) and controls [NRVMs or connexin43 (Cx43)-deficient HeLa cells] vs. Cx43⁺ noncardiomyocytes [aged rat ventricular myofibroblasts (ARVFs) or HeLaCx43 cells]. We performed high-speed voltage-sensitive optical imaging at baseline and following acute H₂O₂ exposure. In NRVM-NRVM and NRVM-HeLa controls, no arrhythmias occurred under either experimental condition. In the NRVM-ARVF and NRVM-HeLaCx43 groups, Cx43⁺ noncardiomyocytes enabled passive decremental propagation of electrical impulses and impaired NRVM activation and repolarization, thereby slowing conduction and prolonging action potential duration. Following H₂O₂ exposure, arrhythmia triggers, automaticity, and non-reentrant and reentrant arrhythmias emerged. This study reveals that myofibroblasts (which generate cardiac fibrosis) and other noncardiomyocytes can induce not only structural remodeling but also electrical remodeling and that electrical remodeling by noncardiomyocytes can be particularly arrhythmogenic in the presence of an oxidative burst. Synergistic electrical remodeling between H₂O₂ and noncardiomyocytes may account for the clinical arrhythmogenicity of myofibroblasts at fibrotic interfaces with cardiomyocytes in ischemic/non-ischemic cardiomyopathies. Understanding the enhanced arrhythmogenicity of synergistic electrical remodeling by H₂O₂ and noncardiomyocytes may guide novel safe-by-design antiarrhythmic strategies for next-generation iatrogenic interfaces between surviving native cardiomyocytes and exogenous stem cells or engineered tissues in cardiac regenerative therapies.

Signature signal strategy: Electrogram-based ventricular tachycardia mapping

Multiple decades of work have recognized complexities of substrates responsible for ventricular tachycardia (VT). There is sufficient evidence that 3 critical components of a re-entrant VT circuit, namely, region of slow conduction, zone of unidirectional block, and exit site, are located in spatial vicinity to each other in the ventricular scar. Each of these components expresses characteristic electrograms in sinus rhythm, at initiation of VT, and during VT, respectively. Despite this, abnormal electrograms are widely targeted without appreciation of these signature electrograms during contemporary VT ablation. Our aim is to stimulate physiology-based VT mapping and a targeted ablation of VT. In this article, we focus on these 3 underappreciated aspects of the physiology of ischemic scar-related VT circuits that have practical applications during a VT ablation procedure. We explore the anatomic and functional elements underlying these distinctive bipolar electrograms, specifically the contribution of tissue branching, conduction restitution, and wave curvature to the substrate, as they pertain to initiation and maintenance of VT. We propose a VT ablation approach based on these 3 electrogram features that can be a potential practical means to recognize critical elements of a VT circuit and target ablation.

Factors Promoting Conduction Slowing as Substrates for Block and Reentry in Infarcted Hearts

The development of effective and safe therapies for scar-related ventricular tachycardias requires a detailed understanding of the mechanisms underlying the conduction block that initiates electrical re-entries associated with these arrhythmias. Conduction block has been often associated with electrophysiological changes that prolong action potential duration (APD) within the border zone (BZ) of chronically infarcted hearts. However, experimental evidence suggests that remodeling processes promoting conduction slowing as opposed to APD prolongation mark the chronic phase. In this context, the substrate for the initial block at the mouth of an isthmus/diastolic channel leading to ventricular tachycardia is unclear. The goal of this study was to determine whether electrophysiological parameters associated with conduction slowing can cause block and re-entry in the BZ. In silico experiments were conducted on two-dimensional idealized infarct tissue as well as on a cohort of postinfarction porcine left ventricular models constructed from ex vivo magnetic resonance imaging scans. Functional conduction slowing in the BZ was modeled by reducing sodium current density, whereas structural conduction slowing was represented by decreasing tissue conductivity and including fibrosis. The arrhythmogenic potential of APD prolongation was also tested as a basis for comparison. Within all models, the combination of reduced sodium current with structural remodeling more often degenerated into re-entry and, if so, was more likely to be sustained for more cycles. Although re-entries were also detected in experiments with prolonged APD, they were often not sustained because of the subsequent block caused by long-lasting repolarization. Functional and structural conditions associated with slow conduction rather than APD prolongation form a potent substrate for arrhythmogenesis at the isthmus/BZ of chronically infarcted hearts. Reduced excitability led to block while slow conduction shortened the wavelength of propagation, facilitating the sustenance of re-entries. These findings provide important insights for models of patient-specific risk stratification and therapy planning.

Dynamics of Pivoting Electrical Waves in a Cardiac Tissue Model

Through a detailed mathematical analysis we seek to advance our understanding of how cardiac tissue conductances govern pivoting (spiral, scroll, rotor, functional reentry) wave dynamics. This is an important problem in cardiology since pivoting waves likely underlie most reentrant tachycardias. The problem is complex, and to advance our methods of analysis we introduce two new tools: a ray tracing method and a moving-interface model. When used in combination with an ionic model, they permit us to elucidate the role played by tissue conductances on pivoting wave dynamics. Specifically we simulate traveling electrical waves with an ionic model that can reproduce the characteristics of plane and pivoting waves in small patches of cardiac tissue. Then ray tracing is applied to the simulated pivoting waves in a manner to expose their real displacement. In this exercise we find loci with special characteristics, as well as zones where a part of a pivoting wave quickly transitions from a regenerative to a non-regenerative propagation mode. The loci themselves and the monitoring of the ionic model state variables in this zone permit to elucidate several aspects of pivoting wave dynamics. We then formulate the moving-interface model based on the information gathered with the above-mentioned analysis. Equipped with a velocity profile v(s), s: distance along of the pivoting wave contour and the steady- state action potential duration (APD) of a plane wave during entrainment, APDss(T), at period T, this simple model can predict: shape, orbit of revolution, rotation period, whether a pivoting wave will break up or not, and whether the tissue will admit pivoting waves or not. Because v(s) and APDss(T) are linked to the ionic model, dynamical analysis with the moving-interface model conveys information on the role played by tissue conductances on pivoting wave dynamics. The analysis conducted here enables us to better understand previous results on the termination of pivoting waves. We surmise the method put forth here could become a means to discover how to alter tissue conductances in a manner to terminate pivoting waves at the origin of reentrant tachycardias.

Trabecular Architecture Determines Impulse Propagation Through the Early Embryonic Mouse Heart

Most embryonic ventricular cardiomyocytes are quite uniform, in contrast to the adult heart, where the specialized ventricular conduction system is molecularly and functionally distinct from the working myocardium. We thus hypothesized that the preferential conduction pathway within the embryonic ventricle could be dictated by trabecular geometry. Mouse embryonic hearts of the Nkx2.5:eGFP strain between ED9.5 and ED14.5 were cleared and imaged whole mount by confocal microscopy, and reconstructed in 3D at 3.4 μm isotropic voxel size. The local orientation of the trabeculae, responsible for the anisotropic spreading of the signal, was characterized using spatially homogenized tensors (3 × 3 matrices) calculated from the trabecular skeleton. Activation maps were simulated assuming constant speed of spreading along the trabeculae. The results were compared with experimentally obtained epicardial activation maps generated by optical mapping with a voltage-sensitive dye. Simulated impulse propagation starting from the top of interventricular septum revealed the first epicardial breakthrough at the interventricular grove, similar to experimentally obtained activation maps. Likewise, ectopic activation from the left ventricular base perpendicular to dominant trabecular orientation resulted in isotropic and slower impulse spreading on the ventricular surface in both simulated and experimental conditions. We conclude that in the embryonic pre-septation heart, the geometry of the A-V connections and trabecular network is sufficient to explain impulse propagation and ventricular activation patterns.

Information theory to tachycardia therapy: electrogram entropy predicts diastolic microstructure of reentrant ventricular tachycardia

There is no known strategy to differentiate which multicomponent electrograms in sinus rhythm maintain reentrant ventricular tachycardia (VT). Low entropy in the voltage breakdown of a multicomponent electrogram can localize conditions suitable for reentry but has not been validated against the classic VT activation mapping. We examined whether low entropy in a late and diversely activated ventricular scar region characterizes and differentiates the diastolic path of VT and represents protected tissue channels devoid of side branches. Intraoperative bipolar electrogram (BiEGM) activation and entropy maps were obtained during sinus rhythm in 17 patients with ischemic cardiomyopathy and compared with diastolic activation paths of VT (total of 39 VTs). Mathematical modeling of a zigzag main channel with side branches was also used to further validate structural representation of low entropy in the ventricular scar. A median of one region per patient (range: 1–2 regions) was identified in sinus rhythm, in which BiEGMwith the latest mean activation time and adjacent minimum entropy were assembled together in a high-activation dispersion region. These regions accurately recognized diastolic paths of 34 VTs, often to multiple inducible VTs within a single individual arrhythmogenic region. In mathematical modeling, side branching from the main channel had a strong influence on the BiEGMcomposition along the main channel. The BiEGMobtained from a long unbranched channel had the lowest entropy compared with those with multiple side branches. In conclusion, among a population of multicomponent sinus electrograms, those that demonstrate low entropy and are delayed colocalize to critical long-protected channels of VT. This information is pertinent for planning VT ablation in sinus rhythm.

NEW & NOTEWORTHY Entropy is a measure to quantify breakdown in information. Electrograms from a protected tissue channel can only possess a few states in their voltage and thus less information. In contrast, current-load interactions from side branches in unprotected channels introduce a number of dissimilar voltage deflections and thus high information. We compare here a mapping approach based on entropy against a rigorous reference standard of activation mapping during VT and entropy was assessed in sinus rhythm.

Quantitative spectral assessment of intracardiac electrogram characteristics associated with post infarct fibrosis and ventricular tachycardia

Background

Post-myocardial infarction (MI) remodeling contributes to increased electrophysiological and structural heterogeneity and arrhythmogenesis. Utilising the post-infarct ovine model our aim was to determine unipolar electrogram frequency characteristics consequent to this remodeling and the development of Ventricular Tachycardia (VT).

Methods and results

Mapping studies were performed on 14 sheep at >1 month post-MI induction. Sheep were divided into VT inducible (n = 7) and non-inducible (n = 7) groups. Multielectrode needles (n = 20) were deployed within and surrounding ventricular scar for electrophysiological assessment of electrogram amplitude and width. Spectral analysis of electrograms was undertaken using wavelet and fast fourier transformations (WFFT) to calculate root mean square (RMS) power intervals spanning 0-300Hz in 20Hz intervals. Quantitative assessment between electrophysiological and histological parameters including collagen density, and structural organization of the myocardium was performed.

Increasing myocardial scar density resulted in attenuation of electrogram amplitude and RMS values. (all p<0.01). Between groups there were no differences in electrogram amplitude (p = 0.37), however WFFT analysis revealed significantly higher RMS values in the VT group (p<0.05) in association with high frequency fractional components of the electrogram. As scar density increased, greater between-group differences in RMS were observed spanning this high frequency (200-280Hz) spectrum and which were proportionally dependent on the degree of structural disorganisation of the myocardium (p<0.001) and number of extrastimuli required to induce VT (p<0.05).

Conclusion

High frequency unipolar electrogram spectral characteristics were quantitatively co-influenced by the presence of fibrosis and degree of myocardial structural dissorganisation and were associated with the propensity for development of VT.

Microheterogeneity-induced conduction slowing and wavefront collisions govern macroscopic conduction behavior: A computational and experimental study

The incidence of cardiac arrhythmias is known to be associated with tissue heterogeneities including fibrosis. However, the impact of microscopic structural heterogeneities on conduction in excitable tissues remains poorly understood. In this study, we investigated how acellular microheterogeneities affect macroscopic conduction under conditions of normal and reduced excitability by utilizing a novel platform of paired in vitro and in silico studies to examine the mechanisms of conduction. Regular patterns of nonconductive micro-obstacles were created in confluent monolayers of the previously described engineered-excitable Ex293 cell line. Increasing the relative ratio of obstacle size to intra-obstacle strand width resulted in significant conduction slowing up to 23.6% and a significant increase in wavefront curvature anisotropy, a measure of spatial variation in wavefront shape. Changes in bulk electrical conductivity and in path tortuosity were insufficient to explain these observed macroscopic changes. Rather, microscale behaviors including local conduction slowing due to microscale branching, and conduction acceleration due to wavefront merging were shown to contribute to macroscopic phenomena. Conditions of reduced excitability led to further conduction slowing and a reversal of wavefront curvature anisotropy due to spatially non-uniform effects on microscopic slowing and acceleration. This unique experimental and computation platform provided critical mechanistic insights in the impact of microscopic heterogeneities on macroscopic conduction, pertinent to settings of fibrotic heart disease.

It is well known that perturbations in the heart structure are associated with the initiation and maintenance of clinically significant cardiac arrhythmia. While previous studies have examined how single structural perturbations affect local electrical conduction, our understanding of how numerous microscopic heterogeneities act in aggregate to alter macroscopic electrical behavior is limited. In this study, we utilized simplified engineered excitable cells that contain the minimal machinery of excitability and can be directly computationally modeled. By pairing experimental and computational studies, we showed that the microscopic branching and collisions of electrical waves slow and speed conduction, respectively, resulting in macroscopic changes in the speed and pattern of electrical activation. These microscale behaviors are significantly altered under reduced excitability, resulting in exaggerated collision effects. Overall, this study helps improve our understanding of how microscopic structural heterogeneities in excitable tissue lead to abnormal action potential propagation, conducive to arrhythmias.

Discretization-dependent model for weakly connected excitable media

Pattern formation has been widely observed in extended chemical and biological processes. Although the biochemical systems are highly heterogeneous, homogenized continuum approaches formed by partial differential equations have been employed frequently. Such approaches are usually justified by the difference of scales between the heterogeneities and the characteristic spatial size of the patterns. Under different conditions, for example, under weak coupling, discrete models are more adequate. However, discrete models may be less manageable, for instance, in terms of numerical implementation and mesh generation, than the associated continuum models. Here we study a model to approach discreteness which permits the computer implementation on general unstructured meshes. The model is cast as a partial differential equation but with a parameter that depends not only on heterogeneities sizes, as in the case of quasicontinuum models, but also on the discretization mesh. Therefore, we refer to it as a discretization-dependent model. We validate the approach in a generic excitable media that simulates three different phenomena: the propagation of action membrane potential in cardiac tissue, in myelinated axons of neurons, and concentration waves in chemical microemulsions.

Murine Electrophysiological Models of Cardiac Arrhythmogenesis

Cardiac arrhythmias can follow disruption of the normal cellular electrophysiological processes underlying excitable activity and their tissue propagation as coherent wavefronts from the primary sinoatrial node pacemaker, through the atria, conducting structures and ventricular myocardium. These physiological events are driven by interacting, voltage-dependent, processes of activation, inactivation, and recovery in the ion channels present in cardiomyocyte membranes. Generation and conduction of these events are further modulated by intracellular Ca2+ homeostasis, and metabolic and structural change. This review describes experimental studies on murine models for known clinical arrhythmic conditions in which these mechanisms were modified by genetic, physiological, or pharmacological manipulation. These exemplars yielded molecular, physiological, and structural phenotypes often directly translatable to their corresponding clinical conditions, which could be investigated at the molecular, cellular, tissue, organ, and whole animal levels. Arrhythmogenesis could be explored during normal pacing activity, regular stimulation, following imposed extra-stimuli, or during progressively incremented steady pacing frequencies. Arrhythmic substrate was identified with temporal and spatial functional heterogeneities predisposing to reentrant excitation phenomena. These could arise from abnormalities in cardiac pacing function, tissue electrical connectivity, and cellular excitation and recovery. Triggering events during or following recovery from action potential excitation could thereby lead to sustained arrhythmia. These surface membrane processes were modified by alterations in cellular Ca2+ homeostasis and energetics, as well as cellular and tissue structural change. Study of murine systems thus offers major insights into both our understanding of normal cardiac activity and its propagation, and their relationship to mechanisms generating clinical arrhythmias.

Myofibroblasts Electrotonically Coupled to Cardiomyocytes Alter Conduction: Insights at the Cellular Level from a Detailed In silico Tissue Structure Model

Fibrotic myocardial remodeling is typically accompanied by the appearance of myofibroblasts (MFBs). In vitro, MFBs were shown to slow conduction and precipitate ectopic activity following gap junctional coupling to cardiomyocytes (CMCs). To gain further mechanistic insights into this arrhythmogenic MFB-CMC crosstalk, we performed numerical simulations in cell-based high-resolution two-dimensional tissue models that replicated experimental conditions. Cell dimensions were determined using confocal microscopy of single and co-cultured neonatal rat ventricular CMCs and MFBs. Conduction was investigated as a function of MFB density in three distinct cellular tissue architectures: CMC strands with endogenous MFBs, CMC strands with coating MFBs of two different sizes, and CMC strands with MFB inserts. Simulations were performed to identify individual contributions of heterocellular gap junctional coupling and of the specific electrical phenotype of MFBs. With increasing MFB density, both endogenous and coating MFBs slowed conduction. At MFB densities of 5-30%, conduction slowing was most pronounced in strands with endogenous MFBs due to the MFB-dependent increase in axial resistance. At MFB densities >40%, very slow conduction and spontaneous activity was primarily due to MFB-induced CMC depolarization. Coating MFBs caused non-uniformities of resting membrane potential, which were more prominent with large than with small MFBs. In simulations of MFB inserts connecting two CMC strands, conduction delays increased with increasing insert lengths and block appeared for inserts >1.2 mm. Thus, electrophysiological properties of engineered CMC-MFB co-cultures depend on MFB density, MFB size and their specific positioning in respect to CMCs. These factors may influence conduction characteristics in the heterocellular myocardium.

Potassium channels in the Cx43 gap junction perinexus modulate ephaptic coupling: an experimental and modeling study

It was recently demonstrated that cardiac sodium channels (Nav1.5) localized at the perinexus, an intercalated disc (ID) nanodomain associated with gap junctions (GJ), may contribute to electrical coupling between cardiac myocytes via an ephaptic mechanism. Impairment of ephaptic coupling by acute interstitial edema (AIE)-induced swelling of the perinexus was associated with arrhythmogenic, anisotropic conduction slowing. Given that Kir2.1 has also recently been reported to localize at intercalated discs, we hypothesized that Kir2.1 channels may reside within the perinexus and that inhibiting them may mitigate arrhythmogenic conduction slowing observed during AIE. Using gated stimulated emission depletion (gSTED) and stochastic optical reconstruction microscopy (STORM) super-resolution microscopy, we indeed find that a significant proportion of Kir2.1 channels resides within the perinexus. Moreover, whereas Nav1.5 inhibition during AIE exacerbated arrhythmogenic conduction slowing, inhibiting Kir2.1 channels during AIE preferentially increased transverse conduction velocity-decreasing anisotropy and ameliorating arrhythmia risk compared to AIE alone. Comparison of our results with a nanodomain computer model identified enrichment of both Nav1.5 and Kir2.1 at intercalated discs as key factors underlying the experimental observations. We demonstrate that Kir2.1 channels are localized within the perinexus alongside Nav1.5 channels. Further, targeting Kir2.1 modulates intercellular coupling between cardiac myocytes, anisotropy of conduction, and arrhythmia propensity in a manner consistent with a role for ephaptic coupling in cardiac conduction. For over half a century, electrical excitation in the heart has been thought to occur exclusively via gap junction-mediated ionic current flow between cells. Further, excitation was thought to depend almost exclusively on sodium channels with potassium channels being involved mainly in returning the cell to rest. Here, we demonstrate that sodium and potassium channels co-reside within nanoscale domains at cell-to-cell contact sites. Experimental and computer modeling results suggest a role for these channels in electrical coupling between cardiac muscle cells via an ephaptic mechanism working in tandem with gap junctions. This new insight into the mechanism of cardiac electrical excitation could pave the way for novel therapies against cardiac rhythm disturbances.

Techniques for automated local activation time annotation and conduction velocity estimation in cardiac mapping

Measurements of cardiac conduction velocity provide valuable functional and structural insight into the initiation and perpetuation of cardiac arrhythmias, in both a clinical and laboratory context. The interpretation of activation wavefronts and their propagation can identify mechanistic properties of a broad range of electrophysiological pathologies. However, the sparsity, distribution and uncertainty of recorded data make accurate conduction velocity calculation difficult. A wide range of mathematical approaches have been proposed for addressing this challenge, often targeted towards specific data modalities, species or recording environments. Many of these algorithms require identification of activation times from electrogram recordings which themselves may have complex morphology or low signal-to-noise ratio. This paper surveys algorithms designed for identifying local activation times and computing conduction direction and speed. Their suitability for use in different recording contexts and applications is assessed.

Spatial characterization of electrogram morphology from transmural recordings in the intact normal heart

PURPOSE

Unipolar (UE) and bipolar electrograms (BE) are utilized to identify arrhythmogenic substrate. We quantified the effect of increasing distance from the source of propagation on local electrogram amplitude; and determined if transmural electrophysiological gradients exist with respect to propagation and stimulation depth.

METHODS

Mapping was performed on 5 sheep. Deployment of >50 quadripolar transmural needles in the LV were located in Cartesian space using Ensite. Contact electrograms from all needles were recorded during multisite bipolar pacing from epicardial then endocardial electrodes. Analysis was performed to determine stimulus distance to local activation time, peak negative amplitude (V-P), and peak-peak amplitude (VP-P) for (1) unfiltered UE, and (2) unfiltered and 30 Hz high-pass filtered BEs. Each sheep was analysed using repeated ANOVA.

RESULTS

Increasing distance from the pacing sites led to significant (p<0.01) attenuation of UEs (V-P = 7.0±0.5%; VP-P = 5.4±0.3% per cm). Attenuation of BE with distance was insignificant (Vp-p unfiltered = 2.2±0.5%; filtered = 1.7±1.4% per cm). Independent of pacing depth, significant (p<0.01) transmural electrophysiological gradients were observed, with highest amplitude occurring at epicardial layers for UE and endocardial layers for BE. Furthermore, during pacing, propagation was earlier at the epicardium than endocardial layer by 1.6±2.0 ms (UE) and 1.4±2.8 ms (BE) (all p>0.01) during endocardial stimulation, and 2.3±2.4 ms (UE) and 1.8±3.7 ms (BE) during epicardal stimulation (all p<0.01).

CONCLUSIONS

Electrogram amplitude is inversely proportional to propagation distance for unipolar modalities only, which affected V-P>VP-P. Conduction propagates preferentially via the epicardium during stimulation and is believed to contribute to a transmural amplitude gradient.

Studying dynamic events in the developing myocardium

Differentiation and conduction properties of the cardiomyocytes are critically dependent on physical conditioning both in vitro and in vivo. Historically, various techniques were introduced to study dynamic events such as electrical currents and changes in ionic concentrations in live cells, multicellular preparations, or entire hearts. Here we review this technological progress demonstrating how each improvement in spatial or temporal resolution provided answers to old and provoked new questions. We further demonstrate how high-speed optical mapping of voltage and calcium can uncover pacemaking potential within the outflow tract myocardium, providing a developmental explanation of ectopic beats originating from this region in the clinical settings.

Determinants of myocardial conduction velocity: implications for arrhythmogenesis

Slowed myocardial conduction velocity (θ) is associated with an increased risk of re-entrant excitation, predisposing to cardiac arrhythmia. θ is determined by the ion channel and physical properties of cardiac myocytes and by their interconnections. Thus, θ is closely related to the maximum rate of action potential (AP) depolarization [(dV/dt)_max], as determined by the fast Na⁺ current (I_Na); the axial resistance (r_a) to local circuit current flow between cells; their membrane capacitances (c_m); and to the geometrical relationship between successive myocytes within cardiac tissue. These determinants are altered by a wide range of pathophysiological conditions. Firstly, I_Na is reduced by the impaired Na⁺ channel function that arises clinically during heart failure, ischemia, tachycardia, and following treatment with class I antiarrhythmic drugs. Such reductions also arise as a consequence of mutations in SCN5A such as those occurring in Lenègre disease, Brugada syndrome (BrS), sick sinus syndrome, and atrial fibrillation (AF). Secondly, r_a, may be increased due to gap junction decoupling following ischemia, ventricular hypertrophy, and heart failure, or as a result of mutations in CJA5 found in idiopathic AF and atrial standstill. Finally, either r_a or c_m could potentially be altered by fibrotic change through the resultant decoupling of myocyte–myocyte connections and coupling of myocytes with fibroblasts. Such changes are observed in myocardial infarction and cardiomyopathy or following mutations in MHC403 and SCN5A resulting in hypertrophic cardiomyopathy (HCM) or Lenègre disease, respectively. This review defines and quantifies the determinants of θ and summarizes experimental evidence that links changes in these determinants with reduced myocardial θ and arrhythmogenesis. It thereby identifies the diverse pathophysiological conditions in which abnormal θ may contribute to arrhythmia.

Structure-function relationship in the sinus and atrioventricular nodes

Recently published optical mapping studies of larger mammals, including humans, have identified functionally discrete sinoatrial exit pathways of activation. This is in line with earlier mapping studies of the dog and the human but in contrast with findings in the mouse and the rabbit, wherein a propagation wave front pattern of activation has been described. It underpins the complex three-dimensional (3D) organization of the cardiac pacemaking and conduction system in larger species, wherein sinoatrial and atrioventricular nodal physiologies both demonstrate identifiable activation pathways, which coincide with anatomic landmarks and histologic architecture, so that in addition to muscle fiber orientation and cell coupling, these intrinsic factors act to determine excitation pathways. This complex 3D organization increases the effect of source-to-sink mismatch both by greater variability in the space constant of tissue and by the 3D projection of this effect in all directions. Mathematical modeling provides a means to study these interactions, and newer models should incorporate these additional factors and their effect into the 3D structure of large mammal physiology.

A finite element approach for modeling micro-structural discontinuities in the heart

The presence of connective tissue as well as interstitial clefts forms a natural barrier to the electrical propagation in the heart. At a microscopic scale, such uncoupling structures change the pattern of the electrical conduction from uniform towards complex and may play a role in the genesis of cardiac arrhythmias. The anatomical diversity of conduction structures and their topology at a microscopic size scale is overwhelming for experimental techniques. Mathematical models have been often employed to study the behavior of the electrical propagation at a sub-cellular level. However, very fine and computationally expensive meshes are required to capture all microscopic details found in the cardiac tissue. In this work, we present a numerical technique based on the finite element method which allows to reproduce the effects of microscopic conduction barriers caused by the presence of uncoupling structures without actually resolving these structures in a high resolution mesh, thereby reducing the computational costs significantly.

Ventricular fibrillation: dynamics and ion channel determinants

Ventricular fibrillation (VF) is the leading cause of sudden cardiac death. This brief review addresses issues relevant to the dynamics of the rotors responsible for functional reentry and VF. It also makes an attempt to summarize present-day knowledge of the manner in which the dynamic interplay between inward and outward transmembrane currents and the heterogeneous cardiac structure establish a substrate for the initiation and maintenance of rotors and VF. The fragmentary nature of our current understanding of ionic VF mechanisms does not even allow an approach toward a "Theory of VF". Yet some hope is provided by recently obtained insight into the roles played in VF by some of the sarcolemmal ion channels that control the excitation-recovery process. For example, strong evidence supports the idea that the interplay between the rapid-inward sodium current and the inward-rectifier potassium current controls rotor formation, as well as rotor stability and frequency. Solid evidence also exists for an involvement of L-type calcium current in the control of rotor frequency and in determining VF-to-ventricular tachycardia conversion. Less clear, however, is whether or not time dependent outward currents through voltage-gated potassium channels affect the fibrillatory process. Hopefully, taking advantage of currently available approaches of structural, molecular and cellular biology, together with computational and imaging techniques, will afford us the opportunity to further advance knowledge on VF mechanisms.

A computer study of the effects of branching dimension on safety factor distribution and propagation in a cardiac conduction network

Branching conduction networks are responsible for coordinated distribution of activation throughout the heart, but the effects of branching geometry on the efficiency of distribution and the safety of propagation remain unclear. We have developed a simplified computer model to investigate this issue in a systemic fashion. A simplified 2D model that reproduces key features of the right atrial pectinate muscle network has been developed. This consists of a large main strand that gives rise to regularly spaced perpendicular branches. Safety factor (SF) and activation time (AT) within the network are estimated for a range of different branch dimensions. The SF is reduced as branch dimension is increased due to the larger current load and is least at the proximal edge of the junction between the main strand and branches. On the other hand, activation efficiency depends on appropriate balance between the load imposed by branching and source which they subsequently provide. With repeated stimulation at progressively decreasing coupling intervals, the SF and conduction fall within the network, but the locations of block varied with branching dimension, again reflecting spatial variation in the balance between current source and current load. We hypothesize that the observed geometry of the branching network optimizes distribution efficiency and propagation safety.

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.

Novel micropatterned cardiac cell cultures with realistic ventricular microstructure

Systematic studies of cardiac structure-function relationships to date have been hindered by the intrinsic complexity and variability of in vivo and ex vivo model systems. Thus, we set out to develop a reproducible cell culture system that can accurately replicate the realistic microstructure of native cardiac tissues. Using cell micropatterning techniques, we aligned cultured cardiomyocytes at micro- and macroscopic spatial scales to follow local directions of cardiac fibers in murine ventricular cross sections, as measured by high-resolution diffusion tensor magnetic resonance imaging. To elucidate the roles of ventricular tissue microstructure in macroscopic impulse conduction, we optically mapped membrane potentials in micropatterned cardiac cultures with realistic tissue boundaries and natural cell orientation, cardiac cultures with realistic tissue boundaries but random cell orientation, and standard isotropic monolayers. At 2 Hz pacing, both microscopic changes in cell orientation and ventricular tissue boundaries independently and synergistically increased the spatial dispersion of conduction velocity, but not the action potential duration. The realistic variations in intramural microstructure created unique spatial signatures in micro- and macroscopic impulse propagation within ventricular cross-section cultures. This novel in vitro model system is expected to help bridge the existing gap between experimental structure-function studies in standard cardiac monolayers and intact heart tissues.

Decreased intercellular coupling improves the function of cardiac pacemakers derived from mouse embryonic stem cells

The aim of this study was to determine if embryonic stem cell derived cardiomyocyte aggregates (ESdCs) can act as pacemakers in spontaneously active cardiomyocyte preparations when their connexin isoform expression is tuned toward a more sinus nodal phenotype. Using microelectrode array recordings (MEAs), we demonstrate that mouse ESdCs establish electrical coupling with spontaneously active cardiomyocyte preparations (HL-1 monolayer) and obtain pacemaker dominance. WT- and Cx43(-/-)-ESdCs comparably established intercellular coupling with cardiac host tissue (Cx43(-/-): 86% vs. WT: 91%). Although both aggregates had a 100% success rate in pacing quiescent cardiac preparations, Cx43(-/-)-ESdCs had an increased likelihood of gaining pacemaker dominance (Cx43(-/-): 40% vs. WT: 13%) in spontaneously active preparations. No differences in size, beating frequency, V(m), or differentiation were detected between WT- and Cx43(-/-)-ESdCs but the intercellular coupling resistance in Cx43(-/-)-ESdCs was significantly increased (Cx43(-/-): 1.2nS vs. WT: 14.8nS). Lack of Cx43 prolonged the time until Cx43(-/-)-ESdCs established frequency synchronization with the host tissue. It further hampered the excitation spread from the cardiomyocyte preparation into the ESdC. However rectifying excitation spread in these co-cultures could not be unequivocally identified. In summary, ESdCs can function as dominant biological pacemakers and Cx43 expression is not a prerequisite for their electrical integration. Maintenance of pacemaker dominance depends critically on the pacemaker's gap junction expression benefiting those with increased intercellular coupling resistances. Our results provide important insight into the design of biological pacemakers that will benefit the use of cardiomyocytes for cell replacement therapy.

Subepicardial phase 0 block and discontinuous transmural conduction underlie right precordial ST-segment elevation by a SCN5A loss-of-function mutation

Two mechanisms are generally proposed to explain right precordial ST-segment elevation in Brugada syndrome: 1) right ventricular (RV) subepicardial action potential shortening and/or loss of dome causing transmural dispersion of repolarization; and 2) RV conduction delay. Here we report novel mechanistic insights into ST-segment elevation associated with a Na(+) current (I(Na)) loss-of-function mutation from studies in a Dutch kindred with the COOH-terminal SCN5A variant p.Phe2004Leu. The proband, a man, experienced syncope at age 22 yr and had coved-type ST-segment elevations in ECG leads V1 and V2 and negative T waves in V2. Peak and persistent mutant I(Na) were significantly decreased. I(Na) closed-state inactivation was increased, slow inactivation accelerated, and recovery from inactivation delayed. Computer-simulated I(Na)-dependent excitation was decremental from endo- to epicardium at cycle length 1,000 ms, not at cycle length 300 ms. Propagation was discontinuous across the midmyocardial to epicardial transition region, exhibiting a long local delay due to phase 0 block. Beyond this region, axial excitatory current was provided by phase 2 (dome) of the M-cell action potentials and depended on L-type Ca(2+) current ("phase 2 conduction"). These results explain right precordial ST-segment elevation on the basis of RV transmural gradients of membrane potentials during early repolarization caused by discontinuous conduction. The late slow-upstroke action potentials at the subepicardium produce T-wave inversion in the computed ECG waveform, in line with the clinical ECG.

Region of slowed conduction acts as core for spiral wave reentry in cardiac cell monolayers

Pathophysiological heterogeneity in cardiac tissue is related to the occurrence of arrhythmias. Of importance are regions of slowed conduction, which have been implicated in the formation of conduction block and reentry. Experimentally, it has been a challenge to produce local heterogeneity in a manner that is both reversible and well controlled. Consequently, we developed a dual-zone superfusion chamber that can dynamically create a small (5 mm) central island of heterogeneity in cultured cardiac cell monolayers. Three different conditions were studied to explore the effect of regionally slowed conduction on wave propagation and reentry: depolarization by elevated extracellular potassium, sodium channel inhibition with lidocaine, and cell-cell decoupling with palmitoleic acid. Using optical mapping of transmembrane voltage, we found that the central region of slowed conduction always served as the core region around which a spiral wave formed and then revolved following a period of rapid pacing. Because of the localized slowing in the core region, we observed experimentally for the first time an S shape of the spiral wave front near its tip. These results indicate that a small region of slowed conduction can play a crucial role in the formation, anchoring, and modulation of reentrant spiral waves.

Dynamic changes of cardiac conduction during rapid pacing

Slow conduction and unidirectional conduction block (UCB) are key mechanisms of reentry. Following abrupt changes in heart rate, dynamic changes of conduction velocity (CV) and structurally determined UCB may critically influence arrhythmogenesis. Using patterned cultures of neonatal rat ventricular myocytes grown on microelectrode arrays, we investigated the dynamics of CV in linear strands and the behavior of UCB in tissue expansions following an abrupt decrease in pacing cycle length (CL). Ionic mechanisms underlying rate-dependent conduction changes were investigated using the Pandit-Clark-Giles-Demir model. In linear strands, CV gradually decreased upon a reduction of CL from 500 ms to 230-300 ms. In contrast, at very short CLs (110-220 ms), CV first decreased before increasing again. The simulations suggested that the initial conduction slowing resulted from gradually increasing action potential duration (APD), decreasing diastolic intervals, and increasing postrepolarization refractoriness, which impaired Na(+) current (I(Na)) recovery. Only at very short CLs did APD subsequently shorten again due to increasing Na(+)/K(+) pump current secondary to intracellular Na(+) accumulation, which caused recovery of CV. Across tissue expansions, the degree of UCB gradually increased at CLs of 250-390 ms, whereas at CLs of 180-240 ms, it first increased and subsequently decreased. In the simulations, reduction of inward currents caused by increasing intracellular Na(+) and Ca(2+) concentrations contributed to UCB progression, which was reversed by increasing Na(+)/K(+) pump activity. In conclusion, CV and UCB follow intricate dynamics upon an abrupt decrease in CL that are determined by the interplay among I(Na) recovery, postrepolarization refractoriness, APD changes, ion accumulation, and Na(+)/K(+) pump function.

Optical imaging of arrhythmias in tissue culture

BACKGROUND

Cardiac cell cultures are becoming an important experimental system of minimal complexity that captures many of the salient features of myocardial tissue function and are simple enough that tissue parameters can be controlled systematically. Fundamental mechanisms that underlie normal and pathological electrophysiology at the tissue level can be studied. Of particular interest are spiral waves, which underlie many tachyarrhythmias and fibrillation.

METHODS

Methods of patterned growth were used to control tissue structure, and contact fluorescence imaging was used to visualize the spread of electrical waves in confluent monolayers of neonatal rat ventricular cells stained with voltage-sensitive dye.

RESULTS

Work is summarized regarding anisotropy, multiarmed spirals, cocultures of cardiac cells and skeletal myoblasts or mesenchymal stem cells, mechanical excitation, attachment of spiral waves to small anatomical obstacles, perturbation of spiral waves by external electric fields, and structure-based facilitation of spiral wave formation.

CONCLUSIONS

The cultured cell monolayer is a contemporary experimental model encompassing great versatility for basic studies of wavefront propagation and cardiac arrhythmias.

Connexin30.2 containing gap junction channels decelerate impulse propagation through the atrioventricular node

In the mammalian heart, gap junction channels between electrically coupled cardiomyocytes are necessary for impulse propagation and coordinated contraction of atria and ventricles. Recently, mouse connexin30.2 (Cx30.2) was shown to be expressed in the cardiac conduction system, predominantly in sinoatrial and atrioventricular (AV) nodes. The corresponding gap junctional channels expressed in HeLa cells exhibit the lowest unitary conductance (9 pS) of all connexin channels. Here we report that Cx30.2 slows down the propagation of excitation through the AV node. Mice expressing a LacZ reporter gene instead of the Cx30.2 coding region (Cx30.2(LacZ/LacZ)) exhibit a PQ interval that is approximately 25% shorter than in WT littermates. By recording atrial, His, and ventricular signals with intracardiac electrodes, we show that this decrease is attributed to significantly accelerated conduction above the His bundle (atrial-His interval: 27.9 +/- 5.1 ms in Cx30.2(LacZ/LacZ) versus 37.1 +/- 4.1 ms in Cx30.2(+/+) mice), whereas HV conduction is unaltered. Atrial stimulation revealed an elevated AV-nodal conduction capacity and faster ventricular response rates during induced episodes of atrial fibrillation in Cx30.2(LacZ/LacZ) mice. Our results show that Cx30.2 contributes to the slowdown of impulse propagation in the AV node and additionally limits the maximum number of beats conducted from atria to ventricles. Thus, it is likely to be involved in coordination of atrial and ventricular contraction and to fulfill a protective role toward pathophysiological states such as atrial tachyarrhythmias (e.g., atrial fibrillation) by preventing rapid conduction to the ventricles potentially associated with hemodynamic deterioration.