Initiation and propagation of ectopic waves: insights from an in vitro model of ischemia-reperfusion injury

Atrial Fibrillation and Chronic Coronary Ischemia: A Challenging Vicious Circle

Atrial fibrillation, the most frequent arrhythmia in clinical practice and chronic coronary syndrome, is one of the forms of coronary ischemia to have a strong dual relationship. Atrial fibrillation may accelerate atherosclerosis and may increase oxygen consumption in the myocardium, creating a mismatch between supply and demand, thus promoting the development or worsening of coronary ischemia. Chronic coronary syndrome alters the structure and function of gap junction proteins, affecting the conduction of action potential and leading to ischemic necrosis of cardiomyocytes and their replacement with fibrous tissue, in this way sustaining the focal ectopic activity in atrial myocardium. They have many risk factors in common, such as hypertension, obesity, type 2 diabetes mellitus, and dyslipidemia. It is vital for the prognosis of patients to break this vicious circle by controlling risk factors, drug therapies, of which antithrombotic therapy may sometimes be challenging in terms of prothrombotic and bleeding risk, and interventional therapies (revascularization and catheter ablation).

Assessing the arrhythmogenic risk of engineered heart tissue patches through in silico application on infarcted ventricle models

BACKGROUND

Post myocardial infarction (MI) ventricles contain fibrotic tissue and may have disrupted electrical properties, both of which predispose to an increased risk of life-threatening arrhythmias. Application of epicardial patches obtained from human-induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) are a potential long-term therapy to treat heart failure resulting from post MI remodelling. However, whether the introduction of these patches is anti- or pro-arrhythmic has not been studied.

METHODS

We studied arrhythmic risk using in silico engineered heart tissue (EHT) patch engraftment on human post-MI ventricular models. Two patient models were studied, including one with a large dense scar and one with an apparent channel of preserved viability bordered on both sides by scar. In each heart model a virtual EHT patch was introduced as a layer of viable tissue overlying the scarred area, with hiPSC-CMs electrophysiological properties. The incidence of re-entrant and sustained activation in simulations with and without EHT patches was assessed and the arrhythmia inducibility compared in the context of different EHT patch properties (conduction velocity (CV) and action potential duration (APD)). The impact of the EHT patch on the likelihood of focal ectopic impulse propagation was estimated by assessing the minimum stimulus strength and duration required to generate a propagating impulse in the scar border zone (BZ) with and without patch.

RESULTS

We uncovered two main mechanisms by which ventricular tachycardia (VT) risk could be either augmented or attenuated by the interaction of the patch with the tissue. In the case of isthmus-related VT, our simulations predict that EHT patches can prevent the induction of VT when the, generally longer, hiPSC-CMs APD is reduced towards more physiological values. In the case of large dense scar, we found that, an EHT patch with CV similar to the host myocardium does not promote VT, while EHT patches with lower CV increase the risk of VT, by promoting both non-sustained and sustained re-entry. Finally, our simulations indicate that electrically coupled EHT patches reduce the likelihood of propagation of focal ectopic impulses.

CONCLUSIONS

The introduction of EHT patches as a treatment for heart failure has the potential to augment or attenuate the risk of ventricular arrhythmias, and variations in the anatomic configuration of the substrate, the functional properties of the BZ and the electrophysiologic properties of the patch itself will determine the overall impact. Planning for delivery of this therapy will need to consider the possible impact on arrhythmia.

Critical Requirements for the Initiation of a Cardiac Arrhythmia in Rat Ventricle: How Many Myocytes?

Cardiovascular disease is the leading cause of death worldwide due in a large part to arrhythmia. In order to understand how calcium dynamics play a role in arrhythmogenesis, normal and dysfunctional Ca²⁺ signaling in a subcellular, cellular, and tissued level is examined using cardiac ventricular myocytes at a high temporal and spatial resolution using multiscale computational modeling. Ca²⁺ sparks underlie normal excitation-contraction coupling. However, under pathological conditions, Ca²⁺ sparks can combine to form Ca²⁺ waves. These propagating elevations of (Ca²⁺)_i can activate an inward Na⁺-Ca²⁺ exchanger current (I_NCX) that contributes to early after-depolarization (EADs) and delayed after-depolarizations (DADs). However, how cellular currents lead to full depolarization of the myocardium and how they initiate extra systoles is still not fully understood. This study explores how many myocytes must be entrained to initiate arrhythmogenic depolarizations in biophysically detailed computational models. The model presented here suggests that only a small number of myocytes must activate in order to trigger an arrhythmogenic propagating action potential. These conditions were examined in 1-D, 2-D, and 3-D considering heart geometry. The depolarization of only a few hundred ventricular myocytes is required to trigger an ectopic depolarization. The number decreases under disease conditions such as heart failure. Furthermore, in geometrically restricted parts of the heart such as the thin muscle strands found in the trabeculae and papillary muscle, the number of cells needed to trigger a propagating depolarization falls even further to less than ten myocytes.

Coronary heart disease and atrial fibrillation: a vicious cycle

The population suffering from coronary heart disease (CHD) complicated by atrial fibrillation (AF) is rising rapidly. A strong correlation between the two diseases has been reported, and the many common risk factors they share may play prominent roles in their development. In addition, CHD can directly promote the progression of AF by affecting reentry formation, focal ectopic activity, and neural remodeling. At the same time, AF also affects CHD through three aspects: 1) atherosclerosis, 2) the mismatch of blood supply and oxygen consumption, and 3) thrombosis. In conclusion, CHD and AF can aggravate each other and seem to form a vicious cycle. For patients with CHD complicated by AF, principal studies and guidelines have focused on antithrombotic treatment and rhythm control, which are paramount for these patients. Of note, our review sheds light on the strategies to break the cycle of the two diseases, which may be fundamental to treat these patients and optimize the benefit.

Variability in electrophysiological properties and conducting obstacles controls re-entry risk in heterogeneous ischaemic tissue

Ischaemia, in which inadequate blood supply compromises and eventually kills regions of cardiac tissue, can cause many types of arrhythmia, some life-threatening. A significant component of this is the effects of the resulting hypoxia, and concomitant hyperklaemia and acidosis, on the electrophysiological properties of myocytes. Clinical and experimental data have also shown that regions of structural heterogeneity (fibrosis, necrosis, fibro-fatty infiltration) can act as triggers for arrhythmias under acute ischaemic conditions. Mechanistic models have successfully captured these effects in silico. However, the relative significance of these separate facets of the condition, and how sensitive arrhythmic risk is to the extents of each, is far less explored. In this work, we use partitioned Gaussian process emulation and new metrics for source-sink mismatch that rely on simulations of bifurcating cardiac fibres to interrogate a model of heterogeneous ischaemic tissue. Re-entries were most sensitive to the level of hypoxia and the fraction of non-excitable tissue. In addition, our results reveal both protective and pro-arrhythmic effects of hyperklaemia, and present the levels of hyperklaemia, hypoxia and percentage of non-excitable tissue that pose the highest arrhythmic risks. This article is part of the theme issue 'Uncertainty quantification in cardiac and cardiovascular modelling and simulation'.

Multicellular In vitro Models of Cardiac Arrhythmias: Focus on Atrial Fibrillation

Atrial fibrillation (AF) is the most common cardiac arrhythmia in clinical practice with a large socioeconomic impact due to its associated morbidity, mortality, reduction in quality of life and health care costs. Currently, antiarrhythmic drug therapy is the first line of treatment for most symptomatic AF patients, despite its limited efficacy, the risk of inducing potentially life-threating ventricular tachyarrhythmias as well as other side effects. Alternative, in-hospital treatment modalities consisting of electrical cardioversion and invasive catheter ablation improve patients' symptoms, but often have to be repeated and are still associated with serious complications and only suitable for specific subgroups of AF patients. The development and progression of AF generally results from the interplay of multiple disease pathways and is accompanied by structural and functional (e.g., electrical) tissue remodeling. Rational development of novel treatment modalities for AF, with its many different etiologies, requires a comprehensive insight into the complex pathophysiological mechanisms. Monolayers of atrial cells represent a simplified surrogate of atrial tissue well-suited to investigate atrial arrhythmia mechanisms, since they can easily be used in a standardized, systematic and controllable manner to study the role of specific pathways and processes in the genesis, perpetuation and termination of atrial arrhythmias. In this review, we provide an overview of the currently available two- and three-dimensional multicellular in vitro systems for investigating the initiation, maintenance and termination of atrial arrhythmias and AF. This encompasses cultures of primary (animal-derived) atrial cardiomyocytes (CMs), pluripotent stem cell-derived atrial-like CMs and (conditionally) immortalized atrial CMs. The strengths and weaknesses of each of these model systems for studying atrial arrhythmias will be discussed as well as their implications for future studies.

Use of GelMA for 3D printing of cardiac myocytes and fibroblasts

Aim

To 3D print heart tissue, one must understand how the main two types of cardiac cells are affected by the printing process.

Materials & methods

Effects of gelatin methacryloyl (GelMA) concentration, extruder pressure and duration of UV exposure on survival of cardiac myocytes and fibroblasts were examined using lactate dehydrogenase and LIVE/DEAD assays, bioluminescence imaging and morphological assessment.

Results & conclusion

Cell survival within 3D printed cardiomyocyte-laden GelMA constructs was more sensitive to extruder pressure and GelMA concentrations than within 3D fibroblast-laden GelMA constructs. Cells within both types of constructs were adversely impacted by the UV curing step. Use of mixed cell populations and enrichment of bioink formulation with fibronectin led to an improvement of cardiomyocyte survival and spreading.

Ectopic beats arise from micro-reentries near infarct regions in simulations of a patient-specific heart model

Ectopic beats are known to be involved in the initiation of a variety of cardiac arrhythmias. Although their location may vary, ectopic excitations have been found to originate from infarct areas, regions of micro-fibrosis and other heterogeneous tissues. However, the underlying mechanisms that link ectopic foci to heterogeneous tissues have yet to be fully understood. In this work, we investigate the mechanism of micro-reentry that leads to the generation of ectopic beats near infarct areas using a patient-specific heart model. The patient-specific geometrical model of the heart, including scar and peri-infarct zones, is obtained through magnetic resonance imaging (MRI). The infarct region is composed of ischemic myocytes and non-conducting cells (fibrosis, for instance). Electrophysiology is captured using an established cardiac myocyte model of the human ventricle modified to describe ischemia. The simulation results clearly reveal that ectopic beats emerge from micro-reentries that are sustained by the heterogeneous structure of the infarct regions. Because microscopic information about the heterogeneous structure of the infarct regions is not available, Monte-Carlo simulations are used to identify the probabilities of an infarct region to behave as an ectopic focus for different levels of ischemia and different percentages of non-conducting cells. From the proposed model, it is observed that ectopic beats are generated when a percentage of non-conducting cells is near a topological metric known as the percolation threshold. Although the mechanism for micro-reentries was proposed half a century ago to be a source of ectopic beats or premature ventricular contractions during myocardial infarction, the present study is the first to reproduce this mechanism in-silico using patient-specific data.

Downregulation of beta-microglobulin to diminish T-lymphocyte lysis of non-syngeneic cell sources of engineered heart tissue constructs

The presence of non-autologous major histocompatibility complex class I (MHC-I) molecules on the surface of the grafted cells is one of the main reasons for their rejection in non-syngeneic hosts. We present a straightforward strategy to decrease the presence of MHC-I by shRNA inhibition of beta-2-microglobulin (B2M), a conservative light chain of MHC-I, on the surface of two main cell types that are used to engineer heart tissue constructs. Engineered heart tissue constructs can be generated by combining mouse WT19 fibroblasts and mouse embryonic stem cell-derived cardiac myocytes (mESC-CM). WT19 fibroblasts were stably transduced with an anti-B2M shRNA, which yielded a cell line with dramatically reduced B2M expression levels (16 ± 11% of mock treated control cell line). Interferon gamma treatment increased the levels of B2M expression by >3-fold in both control and transduced fibroblasts; yet, B2M expression levels still remained very low in the transduced cells. When compared with their unmodified counterparts, transduced fibroblasts caused 5.7-fold lesser activation of cognate T-cells. B2M depletion in mESC-CM was achieved by 72 h transduction with anti-B2M shRNA lentiviral particles. Transduced mESC-CM exhibited regular beating and expressed classical cardiac markers. When compared with their unmodified counterparts, transduced mESC-CM caused 2.5-fold lesser activation of cognate T-cells. In vivo assessment of B2M downregulation was performed by analyzing the preferential survival of B2M-downregulated cells in the intraperitoneal cavity of allogeneic mice. Both B2M-downregulated fibroblasts and B2M-downregulated myocytes survived significantly better when compared to their unmodified counterparts (2.01 ± 0.4 and 5.07 ± 1.6 fold increase in survival, respectively). In contrast, when modified WT19 fibroblasts were injected into the intraperitoneal cavity of syngeneic C57Bl/6 mice, no significant survival advantage was observed. Notably, the preferential survival of B2M-downregulated cells persisted in allogeneic hosts with normal levels of natural killer cells, although the effect was lesser in magnitude. Use of shRNA against beta-2-microglobulin offers a simple and effective approach to minimize immunogenicity of the main cellular components of cardiac tissue constructs in non-syngeneic recipients.

A microstructural model of reentry arising from focal breakthrough at sites of source-load mismatch in a central region of slow conduction

Regions of cardiac tissue that have a combination of focal activity and poor, heterogeneous gap junction coupling are often considered to be arrhythmogenic; however, the relationship between the properties of the cardiac microstructure and patterns of abnormal propagation is not well understood. The objective of this study was to investigate the effect of microstructure on the initiation of reentry from focal stimulation inside a poorly coupled region embedded in more well-coupled tissue. Two-dimensional discrete computer models of ventricular monolayers (1 × 1 cm) were randomly generated to represent heterogeneity in the cardiac microstructure. A small, central poorly coupled patch (0.40 × 0.40 cm) was introduced to represent the site of focal activity. Simulated unipolar electrogram recordings were computed at various points in the tissue. As the gap conductance of the patch decreased, conduction slowed and became increasingly complex, marked by fractionated electrograms with reduced amplitude. Near the limit of conduction block, isolated breakthrough sites occurred at single cells along the patch boundary and were marked by long cell-to-cell delays and negative deflections on electrogram recordings. The strongest determinant of the site of wavefront breakthrough was the connectivity of the brick wall architecture, which enabled current flow through small regions of overlapping cells to drive propagation into the well-coupled zone. In conclusion, breakthroughs at the size scale of a single cell can occur at the boundary of source-load mismatch allowing focal activations from slow conducting regions to produce reentry. These breakthrough regions, identifiable by distinct asymmetric, reduced amplitude electrograms, are sensitive to tissue architecture and may be targets for ablation.

Pacemaker interactions induce reentrant wave dynamics in engineered cardiac culture

Pacemaker interactions can lead to complex wave dynamics seen in certain types of cardiac arrhythmias. We use experimental and mathematical models of pacemakers in heterogeneous excitable media to investigate how pacemaker interactions can be a mechanism for wave break and reentrant wave dynamics. Embryonic chick ventricular cells are cultured in vitro so as to create a dominant central pacemaker site that entrains other pacemakers in the medium. Exposure of those cultures to a potassium channel blocker, E-4031, leads to emergence of peripheral pacemakers that compete with each other and with the central pacemaker. Waves emitted by faster pacemakers break up over the slower pacemaker to form reentrant waves. Similar dynamics are observed in a modified FitzHugh-Nagumo model of heterogeneous excitable media with two distinct sites of pacemaking. These findings elucidate a mechanism of pacemaker-induced reentry in excitable media.

Ratiometric imaging of calcium during ischemia-reperfusion injury in isolated mouse hearts using Fura-2

Background

We present an easily implementable method for measuring Fura-2 fluorescence from isolated mouse hearts using a commercially available switching light source and CCD camera. After calibration, it provides a good estimate of intracellular [Ca²⁺] with both high spatial and temporal resolutions, permitting study of changes in dispersion of diastolic [Ca²⁺], Ca²⁺ transient dynamics, and conduction velocities in mouse hearts. In a proof-of-principle study, we imaged isolated Langendorff-perfused mouse hearts with reversible regional myocardial infarctions.

Methods

Isolated mouse hearts were perfused in the Landendorff-mode and loaded with Fura-2. Hearts were then paced rapidly and subjected to 15 minutes of regional ischemia by ligation of the left anterior descending coronary artery, following which the ligation was removed to allow reperfusion for 15 minutes. Fura-2 fluorescence was recorded at regular intervals using a high-speed CCD camera. The two wavelengths of excitation light were interleaved at a rate of 1 KHz with a computer controlled switching light source to illuminate the heart.

Results

Fura-2 produced consistent Ca²⁺ transients from different hearts. Ligating the coronary artery rapidly generated a well defined region with a dramatic rise in diastolic Ca²⁺ without a significant change in transient amplitude; Ca²⁺ handling normalized during reperfusion. Conduction velocity was reduced by around 50% during ischemia, and did not recover significantly when monitored for 15 minutes following reperfusion.

Conclusions

Our method of imaging Fura-2 from isolated whole hearts is capable of detecting pathological changes in intracellular Ca²⁺ levels in cardiac tissue. The persistent change in the conduction velocities indicates that changes to tissue connectivity rather than altered intracellular Ca²⁺ handling may be underlying the electrical instabilities commonly seen in patients following a myocardial infarction.

Evolution of spiral and scroll waves of excitation in a mathematical model of ischaemic border zone

Abnormal electrical activity from the boundaries of ischemic cardiac tissue is recognized as one of the major causes in generation of ischemia-reperfusion arrhythmias. Here we present theoretical analysis of the waves of electrical activity that can rise on the boundary of cardiac cell network upon its recovery from ischaemia-like conditions. The main factors included in our analysis are macroscopic gradients of the cell-to-cell coupling and cell excitability and microscopic heterogeneity of individual cells. The interplay between these factors allows one to explain how spirals form, drift together with the moving boundary, get transiently pinned to local inhomogeneities, and finally penetrate into the bulk of the well-coupled tissue where they reach macroscopic scale. The asymptotic theory of the drift of spiral and scroll waves based on response functions provides explanation of the drifts involved in this mechanism, with the exception of effects due to the discreteness of cardiac tissue. In particular, this asymptotic theory allows an extrapolation of 2D events into 3D, which has shown that cells within the border zone can give rise to 3D analogues of spirals, the scroll waves. When and if such scroll waves escape into a better coupled tissue, they are likely to collapse due to the positive filament tension. However, our simulations have shown that such collapse of newly generated scrolls is not inevitable and that under certain conditions filament tension becomes negative, leading to scroll filaments to expand and multiply leading to a fibrillation-like state within small areas of cardiac tissue.

High resolution contrast ultrasound and NADH fluorescence imaging of myocardial perfusion in excised rat hearts

Simultaneous imaging of myocardial flow and hypoxia could be vital for identifying acute ischemic mechanisms that may trigger an arrhythmia. We have studied the distribution of flow and hypoxia in excised locally ischemic rat hearts using simultaneous contrast ultrasound imaging and beta-nicotinamide adenine dinucleotide (NADH) fluorescence imaging. Local ischemia was induced by controlling flow within a major coronary artery. Intra-myocardial flow was imaged using contrast high-resolution ultrasound (linear probe; 13-6 MHz). An ultrasound contrast agent (UCA) was used to highlight the ischemic border. We observed distinct borders between two perfusion beds. UCA images showed high contrast borders of flow. The progression of UCA through the tissue was clearly visible. Intramyocardial regions of flow overlap could be identified by superimposing images of UCA from two perfusion zones. Borders between hypoxic and normoxic tissue were clearly revealed by increased NADH fluorescence. Hypoxic borders were oriented along borders of flow. In summary, simultaneous ultrasound and NADH imaging of excised hearts from small animals provide high fidelity images for characterizing the distribution of flow and hypoxic tissue during acute localized ischemia.

Clinically relevant concentrations of di (2-ethylhexyl) phthalate (DEHP) uncouple cardiac syncytium

Di(2-ethylhexyl) phthalate (DEHP) is a widely used plasticizer found in a variety of polyvinyl chloride (PVC) medical products. The results of studies in experimental animals suggest that DEHP leached from flexible PVC tubing may cause health problems in some patient populations. While the cancerogenic and reproductive effects of DEHP are well recognized, little is known about the potential adverse impact of phthalates on the heart. This study examined the effects of clinically relevant concentrations of DEHP on neonatal rat cardiomyocytes. It was found that application of DEHP to a confluent, synchronously beating cardiac cell network, leads to a marked, concentration-dependent decrease in conduction velocity and asynchronous cell beating. The mechanism behind these changes was a loss of gap junctional connexin-43, documented using Western blot analysis, dye-transfer assay and immunofluorescence. In addition to its effect on electrical coupling, DEHP treatment also affected the mechanical movement of myocyte layers. The latter was linked to the decreased stiffness of the underlying fibroblasts, as the amount of triton-insoluble vimentin was significantly decreased in DEHP-treated samples. The data indicate that DEHP, in clinically relevant concentrations, can impair the electrical and mechanical behavior of a cardiac cell network. Applicability of these findings to human patients remains to be established.

Electrophysiological consequences of acute regional ischemia/reperfusion in neonatal rat ventricular myocyte monolayers

BACKGROUND

Electrophysiological changes promoting arrhythmias during acute regional ischemia/reperfusion are challenging to study in intact cardiac tissue because of complex 3-dimensional myocardial and vascular geometry. We characterized electrophysiological alterations and arrhythmias during regional ischemia/reperfusion in a simpler 2-dimensional geometry of cultured neonatal rat ventricular myocyte monolayers.

METHODS AND RESULTS

Optical mapping of intracellular Ca (Ca(i)) and voltage was performed with the use of Rhod 2-AM and Rh-237, respectively. Regional ischemia was mimicked by covering the central portion of monolayer with a glass coverslip, and reperfusion was mimicked by removing the coverslip. Monolayers were stained with fluorescent antibodies to detect total and dephosphorylated connexin-43 at various time points. During coverslip ischemia, action potential duration shortened, Ca(i) transient duration was prolonged, and local conduction velocity (CV) slowed progressively, with loss of excitability after 10.6 +/- 3.6 minutes. CV slowing was accompanied by connexin-43 dephosphorylation. During ischemia, spontaneous reentry occurred in 5 of 11 monolayers, initiated by extrasystoles arising from the border zone or unidirectional conduction block of paced beats. On reperfusion, excitability recovered within 1.0 +/- 0.8 minutes, but CV remained depressed for 9.0 +/- 3.0 minutes, promoting reentry in the reperfused zone. As connexin-43 phosphorylation recovered in the reperfused zone, CV normalized, and arrhythmias resolved.

CONCLUSIONS

Acute regional ischemia/reperfusion in neonatal rat ventricular myocyte monolayers recapitulates electrophysiological alterations and arrhythmias similar to those observed during acute coronary occlusion/reperfusion in intact hearts. During early reperfusion, slow recovery from connexin-43 dephosphorylation leads to persistent CV slowing, creating a highly arrhythmogenic substrate.

Locations of ectopic beats coincide with spatial gradients of NADH in a regional model of low-flow reperfusion

We studied the origins of ectopic beats during low-flow reperfusion after acute regional ischemia in excised rat hearts. The left anterior descending coronary artery was cannulated. Perfusate was delivered to the cannula using an high-performance liquid chromatography pump. This provided not only precise control of flow rate but also avoided mechanical artifacts associated with vessel occlusion and deocclusion. Optical mapping of epicardial transmembrane potential served to identify activation wavefronts. Imaging of NADH fluorescence was used to quantify local ischemia. Our experiments suggest that low-flow reperfusion of ischemic myocardium leads to a highly heterogeneous ischemic substrate and that the degree of ischemia between adjacent patches of tissue changes in time. In contrast to transient ectopic activity observed during full-flow reperfusion, persistent ectopic arrhythmias were observed during low-flow reperfusion. The origins of ectopic beats were traceable to areas of high spatial gradients of changes in NADH fluorescence caused by low-flow reperfusion.

Controlled regional hypoperfusion in Langendorff heart preparations

We describe a new approach that combines several techniques to allow abnormal electrical and calcium activity to be visualized within hypoperfused myocardial tissue. A flexible microcannula was inserted into the left anterior descending artery of Langendorff perfused rat hearts, an air-tight seal between the coronary artery and the cannula was created, and an HPLC pump was used to deliver a specified flowrate through the microcannula. High resolution optical mapping of NADH/calcium, NADH/voltage or calcium/voltage was then conducted using a dual camera system. The ECG was acquired using surface electrodes. This perfusion technique is superior to occluding a vessel by either a tie or a clamp because it allows precise control of the composition and amount of flow to a defined ischemic bed. Another advantage is that flow can be stopped and resumed remotely, without touching the heart. This allows ectopic beats, or other arrhythmogenic activity, such as alternans, to be recorded immediately after changes in flow are imposed. Altogether, the described method provides a powerful new tool to assess how coronary flow rate affects the degree of local ischemia by the ability to record abnormal patterns of electrical activity and associated intracellular calcium transients with high spatiotemporal resolution from epicardial areas as small as 100 x 100 microm.

Interaction between spiral and paced waves in cardiac tissue

For prevention of lethal arrhythmias, patients at risk receive implantable cardioverter-defibrillators, which use high-frequency antitachycardia pacing (ATP) to convert tachycardias to a normal rhythm. One of the suggested ATP mechanisms involves paced-induced drift of rotating waves followed by their collision with the boundary of excitable tissue. This study provides direct experimental evidence of this mechanism. In monolayers of neonatal rat cardiomyocytes in which rotating waves of activity were initiated by premature stimuli, we used the Ca(2+)-sensitive indicator fluo 4 to observe propagating wave patterns. The interaction of the spiral tip with a paced wave was then monitored at a high spatial resolution. In the course of the experiments, we observed spiral wave pinning to local heterogeneities within the myocyte layer. High-frequency pacing led, in a majority of cases, to successful termination of spiral activity. Our data show that 1) stable spiral waves in cardiac monolayers tend to be pinned to local heterogeneities or areas of altered conduction, 2) overdrive pacing can shift a rotating wave from its original site, and 3) the wave break, formed as a result of interaction between the spiral tip and a paced wave front, moves by a paced-induced drift mechanism to an area where it may become unstable or collide with a boundary. The data were complemented by numerical simulations, which was used to further analyze experimentally observed behavior.

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.

Optical mapping of V(m) and Ca(i)(2+) in a model of arrhythmias induced by local catecholamine application in patterned cell cultures

Catecholamines are known to provoke cardiac arrhythmias, but important aspects such as localization of the arrhythmia source in multicellular tissue and exact ionic mechanisms are not well-known. In this work, a multicellular model of arrhythmias caused by local epinephrine application was developed; V (m) and Ca(i)(2+) changes at the arrhythmia source were measured using fluorescent dyes and high-resolution optical mapping. Cultured strands of neonatal rat myocytes (width approximately 0.4 mm) were produced by patterned growth. Epinephrine (1 micromol/l) was applied over an area of 0.3-0.6 mm via two micropipettes, and strands were stimulated by burst pacing. Local epinephrine application caused triggered arrhythmias with cycle lengths of 202-379 ms and duration of >10 s in 9 out of 16 preparations. Optical V(m) mapping demonstrated that in 78% of cases, the source of arrhythmia was located at the boundary of the locally perfused area. Staining with Ca(i)(2+)-sensitive dye Fluo-4 prevented arrhythmia induction in most cases (85%) likely due to Ca(2+) buffering by the dye. Optical Ca(i)(2+) mapping revealed non-propagated Ca(i)(2+) oscillations at the boundary of the locally perfused area in 45% cases. In conclusion, we developed a new model of catecholamine-dependent arrhythmias allowing mapping of V(m) and Ca(i)(2+) at the arrhythmia source with microscopic resolution. The arrhythmias typically originated from the boundary of the epinephrine-perfused area. The location of the arrhythmia source correlated with localized Ca(i)(2+) oscillations suggesting that arrhythmias were caused by Ca(i)(2+) overload at these locations.

Protective role of protein kinase C epsilon activation in ischemia-reperfusion arrhythmia

PURPOSE

Ischemic heart disease carries an increased risk of malignant ventricular tachycardia (VT), fibrillation (VF), and sudden cardiac death. Protein kinase C (PKC) epsilon activation has been shown to improve the hemodynamics in hearts subjected to ischemia/reperfusion. However, very little is known about the role of epsilon PKC in reperfusion arrhythmias. Here we show that epsilon PKC activation is anti-arrhythmic and its inhibition is pro-arrhythmic.

METHOD

Langendorff-perfused isolated hearts from epsilonPKC agonist (epsilonPKC activation), antagonist (epsilonPKC inhibition) transgenic (TG), and wild-type control mice were subjected to 30 min stabilization period, 10 min global ischemia, and 30 min reperfusion. Action potentials (APs) and calcium transients (CaiT) were recorded simultaneously at 37 degrees C using optical mapping techniques. The incidence of VT and VF was assessed during reperfusion.

RESULTS

No VT/VF was seen in any group during the stabilization period in which hearts were perfused with Tyrode's solution. Upon reperfusion, 3 out of the 16 (19%) wild-type mice developed VT but no VF. In epsilonPKC antagonist group, in which epsilonPKC activity was downregulated, 10 out of 13 (76.9%) TG mice developed VT, of which six (46.2%) degenerated into sustained VF upon reperfusion. Interestingly, in epsilonPKC agonist mice, in which the activity of epsilonPKC was upregulated, no VF was observed and only 1 out of 12 mice showed only transient VT during reperfusion. During ischemia and reperfusion, CaiT decay was exceedingly slower in the antagonist mice compared to the other two groups.

CONCLUSION

Moderate in vivo activation of epsilonPKC exerts beneficial antiarrhythmic effect vis-a-vis the lethal reperfusion arrhythmias. Abnormal CaiT decay may, in part, contribute to the high incidence of reperfusion arrhythmias in the antagonist mice. These findings have important implications for the development of PKC isozyme targeted therapeutics and subsequently for the treatment of ischemic heart diseases.

Macroscopic optical mapping of excitation in cardiac cell networks with ultra-high spatiotemporal resolution

Optical mapping of cardiac excitation using voltage- and calcium-sensitive dyes has allowed a unique view into excitation wave dynamics, and facilitated scientific discovery in the cardiovascular field. At the same time, the structural complexity of the native heart has prompted the design of simplified experimental models of cardiac tissue using cultured cell networks. Such reduced experimental models form a natural bridge between single cells and tissue/organ level experimental systems to validate and advance theoretical concepts of cardiac propagation and arrhythmias. Macroscopic mapping (over >1cm(2) areas) of transmembrane potentials and intracellular calcium in these cultured cardiomyocyte networks is a relatively new development and lags behind whole heart imaging due to technical challenges. In this paper, we review the state-of-the-art technology in the field, examine specific aspects of such measurements and outline a rational system design approach. Particular attention is given to recent developments of sensitive detectors allowing mapping with ultra-high spatiotemporal resolution (>5 megapixels/s). Their interfacing with computer platforms to match the high data throughput, unique for this new generation of detectors, is discussed here. This critical review is intended to guide basic science researchers in assembling optical mapping systems for optimized macroscopic imaging with high resolution in a cultured cell setting. The tools and analysis are not limited to cardiac preparations, but are applicable for dynamic fluorescence imaging in networks of any excitable media.

Genesis of ectopic waves: role of coupling, automaticity, and heterogeneity

Many arrhythmias are believed to be triggered by ectopic sources arising from the border of the ischemic tissue. However, the development of ectopic activity from individual sources to a larger mass of cardiac tissue remains poorly understood. To address this critical issue, we used monolayers of neonatal rat cardiomyocytes to create conditions that promoted progression of ectopic activity from single cells to the network that consisted of hundreds of cells. To explain complex spatiotemporal patterns observed in these experiments we introduced a new theoretical framework. The framework's main feature is a parameter space diagram, which uses cell automaticity and coupling as two coordinates. The diagram allows one to depict network behavior, quantitatively address the heterogeneity factor, and evaluate transitions between different regimes. The well-organized wave trains were observed at moderate and high cell coupling values and network heterogeneity was found to be qualitatively unimportant for these regimes. In contrast, at lower values of coupling, spontaneous ectopic activity led to the appearance of fragmented ectopic waves. For these regimes, network heterogeneity played an essential role. The ectopic waves occasionally gave rise to spiral activity in two different regions within the parameter space via two distinct mechanisms. Together, our results suggest that localized ectopic waves represent an essential step in the progression of ectopic activity. These studies add to the understanding of initiation and progression of arrhythmias and can be applied to other phenomena that deal with assemblies of coupled oscillators.

The sensitivity of the heart to static magnetic fields

Static magnetic fields induce flow potentials in arterial flows in and around the heart, that have been detected as distortions in the ECG. The resultant currents flowing through the myocardium could alter the rate or rhythm of the heart. No such changes have been seen in animal experiments, or with humans, in static fields up to 8 T. The possible effects of such currents induced by fields larger than 8 T on cardiac pacemaker rate, and arrhythmogenesis are reviewed, using virtual cardiac tissues-computational models of cardiac electrophysiology. Arrhythmogenesis can be by the initiation of ectopic beats, or by re-entry, whose probability of occurrence is increased by any increase in the electrical heterogeneity, in particular, the action potential duration heterogeneity of the ventricle. Focal ectopic activity would be readily detectable, but since re-entrant arrhythmias are very rare events, even a large increase in their probability of occurrence still leaves them unlikely to be observed. Both of these two arrhythmogenic mechanisms would show a steep sigmoidal, or threshold dependence on induced current intensity, with the threshold for increasing the vulnerability to re-entry less than the threshold for initiating activity. Failure to observe them at fields less than 8 T provides only a lower bound for any threshold for arrhythmogenesis.

Multiple injury approach and its use for toxicity studies

We present an experimental approach that allows exposure of cells plated on a single coverslip to multiple distinct environments. The original chamber design created a small region of injury using geometrically defined flows of the control and ischemic solutions. Modifications of the original chamber design presented in this article produce a range of flow patterns that can be advantageous for a variety of imaging applications. These applications include: experiments that address effects of different treatments applied to a cell network, parallel testing of negative and positive controls using a single coverslip, border effect studies, evaluation of the treatment's reversibility, and simultaneous monitoring of a cell layer loaded with different fluorescent indicators. The method also can be used to reveal both micro- and macroscopic features of propagation, conduction, and cell coupling in a normal or altered cardiac cell network. These possibilities are illustrated in cultures of neonatal rat cardiomyocytes using oxidant- and calcium-sensitive fluorescent indicators.

Pinning force in active media

Pinning of vortices by defects plays an important role in various physical (superconductivity, superfluidity, etc.) or biological (propagation in cardiac muscle) situations. Which defects act as pinning centers? We propose a way to study this general problem by using an advection field to quantify the attraction between an obstacle and a vortex. A full solution is obtained for the real Ginzburg-Landau equation (RGLE). Two pinning mechanisms are found in excitable media. Our results suggest strong analogies with the RGLE when the heterogeneity is excitable. Unpinning from an unexcitable obstacle is qualitatively harder, resulting in a stronger pinning force. We discuss the implications of our results to control vortices and propose experiments in a chemical active medium and in cardiac tissue.

Scaffold topography alters intracellular calcium dynamics in cultured cardiomyocyte networks

Structural and functional changes ensue in cardiac cell networks when cells are guided by three-dimensional scaffold topography. We report enhanced synchronous pacemaking activity in association with slow diastolic rise in intracellular Ca2+concentration ([Ca2+]i) in cell networks grown on microgrooved scaffolds. Topography-driven changes in cardiac electromechanics were characterized by the frequency dependence of [Ca2+]iin syncytial structures formed of ventricular myocytes cultured on microgrooved elastic scaffolds (G). Cells were electrically paced at 0.5–5 Hz, and [Ca2+]iwas determined using microscale ratiometric (fura 2) fluorescence. Compared with flat (F) controls, the G networks exhibited elevated diastolic [Ca2+]iat higher frequencies, increased systolic [Ca2+]iacross the entire frequency range, and steeper restitution of Ca2+transient half-width ( n = 15 and 7 for G and F, respectively, P < 0.02). Significant differences in the frequency response of force-related parameters were also found, e.g., overall larger total area under the Ca2+transients and faster adaptation of relaxation time to pacing rate ( P < 0.02). Altered [Ca2+]idynamics were paralleled by higher occurrence of spontaneous Ca2+release and increased sarcoplasmic reticulum load ( P < 0.02), indirectly assessed by caffeine-triggered release. Electromechanical instabilities, i.e., Ca2+and voltage alternans, were more often observed in G samples. Taken together, these findings 1) represent some of the first functional electromechanical data for this in vitro system and 2) demonstrate direct influence of the microstructure on cardiac function and susceptibility to arrhythmias via Ca2+-dependent mechanisms. Overall, our results substantiate the idea of guiding cellular phenotype by cellular microenvironment, e.g., scaffold design in the context of tissue engineering.

Behavior of ectopic surface: effects of beta-adrenergic stimulation and uncoupling

By using both experimental and theoretical means, we have addressed the progression of ectopic activity from individual cardiac cells to a multicellular two-dimensional network. Experimental conditions that favor ectopic activity have been created by local perfusion of a small area of cardiomyocyte network (I-zone) with an isoproterenol-heptanol containing solution. The application of this solution initially slowed down and then fully blocked wave propagation inside the I-zone. After a brief lag period, ectopically active cells appeared in the I-zone, followed by evolution of the ectopic clusters into slowly propagating waves. The changing pattern of colliding and expanding ectopic waves confined to the I-zone persisted for as long as the isoproterenol-heptanol environment was present. On restoration of the control environment, the ectopic waves from the I-zone broke out into the surrounding network causing arrhythmias. The observed sequence of events was also modeled by FitzHugh-Nagumo equations and included a cell's arrangement of two adjacent square regions of 20 x 20 cells. The control zone consisted of well-connected, excitable cells, and the I-zone was made of weakly coupled cells (heptanol effect), which became spontaneously active as time evolved (isoproterenol effect). The dynamic events in the system have been studied numerically with the use of a finite difference method. Together, our experimental and computational data have revealed that the combination of low coupling, increased excitability, and spatial heterogeneity can lead to the development of ectopic waves confined to the injured network. This transient condition appears to serve as an essential step for the ectopic activity to "mature" before escaping into the surrounding control network.

Spontaneous Initiation and Termination of Complex Rhythms in Cardiac Cell Culture

INTRODUCTION

Complex cardiac arrhythmias often start and stop spontaneously. These poorly understood behaviors frequently are associated with pathologic modification of the structural heterogeneity and functional connectivity of the myocardium. To evaluate underlying mechanisms, we modify heterogeneity by varying the confluence of embryonic chick monolayer cultures that display complex bursting behaviors. A simple mathematical model was developed that reproduces the experimental behaviors and reveals possible generic mechanisms for bursting dynamics in heterogeneous excitable systems.

METHODS AND RESULTS

Wave propagation was mapped in embryonic chick myocytes monolayers using calcium-sensitive dyes. Monolayer confluence was varied by plating cultures with different cell densities and by varying times in culture. At high plating densities, waves propagate without breaks, whereas monolayers plated at low densities display spirals with frequent breaks and irregular activation fronts. Monolayers at intermediate densities display bursting rhythms in which there is paroxysmal starting and stopping of spiral waves of activity. Similar spatiotemporal patterns of activity were also observed as a function of the time in culture; irregular activity dominates the first 30 hours, followed by repetitive bursting dynamics until 54 hours, after which periodic target patterns or stable spirals prevail. In some quiescent cultures derived from older embryos, it was possible to trigger pacemaker activity following a single activation. We are able to reproduce all of these behaviors by introducing spatial heterogeneity and varying neighborhood size, equivalent to cell connectivity, in a spontaneous cellular automaton model containing a rate-dependent fatigue term.

CONCLUSION

We observe transitions from irregular propagating waves, to spiral waves that spontaneously start and stop, to target waves originating from localized pacemakers in cell culture and a simple theoretical model of heterogeneous excitable media. The results show how physiologic properties of spontaneous activity, heterogeneity, and fatigue can give rise to a wide range of different complex dynamic behaviors similar to clinically observed cardiac arrhythmias.

Reentrant waves in a ring of embryonic chick ventricular cells imaged with a Ca2+ sensitive dye

According to the classic model initially formulated by Mines, reentrant cardiac arrhythmias may be associated with waves circulating in a ring geometry. This study was designed to study the dynamics of reentry in a ring geometry of cardiac tissue culture. Reentrant calcium waves in rings of cultured embryonic chick cardiac myocytes were imaged using a macroscope to monitor the fluorescence of intracellular Calcium Green-1 dye. The rings displayed a variety of stable rhythms including pacemaker activity and spontaneous reentry. Waves originating from a localized pacemaker could lead to reentry as a consequence of unidirectional block. In addition, more complex patterns were observed due to the interactions between reentrant and pacemaker rhythms. These rhythms included instances in which pacemakers accelerated the reentrant rhythm, and instances in which the excitation was blocked in the vicinity of pacemakers. During reentrant activity an appropriately timed electrical stimulus could induce resetting of activity or cause complete annihilation of the propagating waves. This experimental preparation reveals many spontaneously occuring complex rhythms. These complex rhythms are hypothesized to reflect interactions between spontaneous pacemakers, wave propagation, refractory period, and overdrive suppression. This preparation may serve as a useful model system to further investigate complex dynamics arising during reentrant rhythms in cardiac tissue.