Synthesis of voltage-sensitive fluorescence signals from three-dimensional myocardial activation patterns

Extracellular Perinexal Separation Is a Principal Determinant of Cardiac Conduction

BACKGROUND

Cardiac conduction is understood to occur through gap junctions. Recent evidence supports ephaptic coupling as another mechanism of electrical communication in the heart. Conduction via gap junctions predicts a direct relationship between conduction velocity (CV) and bulk extracellular resistance. By contrast, ephaptic theory is premised on the existence of a biphasic relationship between CV and the volume of specialized extracellular clefts within intercalated discs such as the perinexus. Our objective was to determine the relationship between ventricular CV and structural changes to micro- and nanoscale extracellular spaces.

METHODS

Conduction and Cx43 (connexin43) protein expression were quantified from optically mapped guinea pig whole-heart preparations perfused with the osmotic agents albumin, mannitol, dextran 70 kDa, or dextran 2 MDa. Peak sodium current was quantified in isolated guinea pig ventricular myocytes. Extracellular resistance was quantified by impedance spectroscopy. Intercellular communication was assessed in a heterologous expression system with fluorescence recovery after photobleaching. Perinexal width was quantified from transmission electron micrographs.

RESULTS

CV primarily in the transverse direction of propagation was significantly reduced by mannitol and increased by albumin and both dextrans. The combination of albumin and dextran 70 kDa decreased CV relative to albumin alone. Extracellular resistance was reduced by mannitol, unchanged by albumin, and increased by both dextrans. Cx43 expression and conductance and peak sodium currents were not significantly altered by the osmotic agents. In response to osmotic agents, perinexal width, in order of narrowest to widest, was albumin with dextran 70 kDa; albumin or dextran 2 MDa; dextran 70 kDa or no osmotic agent, and mannitol. When compared in the same order, CV was biphasically related to perinexal width.

CONCLUSIONS

Cardiac conduction does not correlate with extracellular resistance but is biphasically related to perinexal separation, providing evidence that the relationship between CV and extracellular volume is determined by ephaptic mechanisms under conditions of normal gap junctional coupling.

Optical mapping of the pig heart in situ under artificial blood circulation

The emergence of optical imaging has revolutionized the investigation of cardiac electrical activity and associated disorders in various cardiac pathologies. The electrical signals of the heart and the propagation pathways are crucial for elucidating the mechanisms of various cardiac pathological conditions, including arrhythmia. The synthesis of near-infrared voltage-sensitive dyes and the voltage sensitivity of the FDA-approved dye Cardiogreen have increased the importance of optical mapping (OM) as a prospective tool in clinical practice. We aimed to develop a method for the high-spatiotemporal-resolution OM of the large animal hearts in situ using di-4-ANBDQBS and Cardiogreen under patho/physiological conditions. OM was adapted to monitor cardiac electrical behaviour in an open-chest pig heart model with physiological or artificial blood circulation. We detail the methods and display the OM data obtained using di-4-ANBDQBS and Cardiogreen. Activation time, action potential duration, repolarization time and conduction velocity maps were constructed. The technique was applied to track cardiac electrical activity during regional ischaemia and arrhythmia. Our study is the first to apply high-spatiotemporal-resolution OM in the pig heart in situ to record cardiac electrical activity qualitatively under artificial blood perfusion. The use of an FDA-approved voltage-sensitive dye and artificial blood perfusion in a swine model, which is generally accepted as a valuable pre-clinical model, demonstrates the promise of OM for clinical application.

Examination of the Effects of Conduction Slowing on the Upstroke of Optically Recorded Action Potentials

Introduction

The upstroke of optical action potentials (APs) recorded from intact hearts are generally recognized to be slower than those recorded with microelectrodes. This is thought to reflect spatial signal averaging within the volume of tissue that makes up the optical signal. However, to date, there has been no direct experimental study on the relationship between conduction velocity (CV) and optical AP upstroke morphology in the intact heart. Notably, it is known that sodium channel block and gap junction inhibition, which both slow CV, exert differential effects on the upstroke velocity of microelectrode-recorded APs. Whether such differences are evident in optical APs is not known. The present study sought to determine the relationship between tissue CV and optical AP upstroke velocity in intact mouse hearts.

Materials and Methods

Isolated, perfused mouse hearts were stained with the potentiometric dye Rh-237. Fluorescent signals were recorded from across the anterior surface of the left and right ventricles during constant pacing. Maximum rate of change in fluorescence (dF/dt_max) and tissue CV were assessed in control conditions, during an acute period of low-flow ischemia, and following perfusion of flecainide (1–3 μmol/L), a sodium channel blocker, or carbenoxolone (10–50 μmol/L), a gap junction inhibitor.

Results

During epicardial pacing, an anisotropic pattern was observed in both activation and dF/dt_max maps, with more rapid optical AP upstroke velocities orientated along the fastest conduction paths (and vice versa). Low-flow ischemia resulted in a time-dependent slowing of ventricular CV, which was accompanied by a concomitant reduction in optical AP upstroke velocity. All values returned to baseline on tissue reperfusion. Both flecainide and carbenoxolone were associated with a concentration-dependent reduction in CV and decrease in optical AP upstroke velocity, despite distinct mechanisms of action. Similar responses to carbenoxolone were observed for low- (156 μm pixel with) and high- (20 μm pixel width) magnification recordings. Comparison of data from all interventions revealed a linear relationship between CV and upstroke dF/dt.

Conclusion

In intact mouse hearts, slowing of optical AP upstroke velocity is directly proportional to the change in CV associated with low-flow ischemia, sodium channel block, and gap junction inhibition.

Acetylcholine Delays Atrial Activation to Facilitate Atrial Fibrillation

Background

Acetylcholine (ACh) shortens action potential duration (APD) in human atria. APD shortening facilitates atrial fibrillation (AF) by reducing the wavelength for reentry. However, the influence of ACh on electrical conduction in human atria and its contribution to AF are unclear, particularly when combined with impaired conduction from interstitial fibrosis.

Objective

To investigate the effect of ACh on human atrial conduction and its role in AF with computational, experimental, and clinical approaches.

Methods

S1S2 pacing (S1 = 600 ms and S2 = variable cycle lengths) was applied to the following human AF computer models: a left atrial appendage (LAA) myocyte to quantify the effects of ACh on APD, maximum upstroke velocity (V _max ), and resting membrane potential (RMP); a monolayer of LAA myocytes to quantify the effects of ACh on conduction; and 3) an intact left atrium (LA) to determine the effects of ACh on arrhythmogenicity. Heterogeneous ACh and interstitial fibrosis were applied to the monolayer and LA models. To corroborate the simulations, APD and RMP from isolated human atrial myocytes were recorded before and after 0.1 μM ACh. At the tissue level, LAAs from AF patients were optically mapped ex vivo using Di-4-ANEPPS. The difference in total activation time (AT) was determined between AT initially recorded with S1 pacing, and AT recorded during subsequent S1 pacing without (n = 6) or with (n = 7) 100 μM ACh.

Results

In LAA myocyte simulations, S1 pacing with 0.1 μM ACh shortened APD by 41 ms, hyperpolarized RMP by 7 mV, and increased V _max by 27 mV/ms. In human atrial myocytes, 0.1 μM ACh shortened APD by 48 ms, hyperpolarized RMP by 3 mV, and increased V _max by 6 mV/ms. In LAA monolayer simulations, S1 pacing with ACh hyperpolarized RMP to delay total AT by 32 ms without and 35 ms with fibrosis. This led to unidirectional conduction block and sustained reentry in fibrotic LA with heterogeneous ACh during S2 pacing. In AF patient LAAs, S1 pacing with ACh increased total AT from 39.3 ± 26 ms to 71.4 ± 31.2 ms (p = 0.036) compared to no change without ACh (56.7 ± 29.3 ms to 50.0 ± 21.9 ms, p = 0.140).

Conclusion

In fibrotic atria with heterogeneous parasympathetic activation, ACh facilitates AF by shortening APD and slowing conduction to promote unidirectional conduction block and reentry.

Isosbestic Point in Optical Mapping; Theoretical and Experimental Determination With Di-4-ANBDQPQ Transmembrane Voltage Sensitive Dye

Optical mapping methods utilize fluorescence dyes to measure dynamic response of cardiac tissue such as changes in transmembrane potential (V_m). For the commonly used V_m sensitive dyes, a dye absorption and emission spectra shift as V_m changes. Signals relevant to V_m are calculated as a relative fluorescence change with respect to the fluorescence baseline. The amplitude of the change depends on the long-pass (LP) filter cut-on wavelength, placed on the sensor side, and the excitation wavelength. An excitation wavelength near the absorption peak, termed the isosbestic point, results in minimal absorption coefficient change as absorption spectra shifts. Consequentially the fluorescence intensity virtually does not change, when fluorescence across the entire emission spectra is measured, irrelevant of V_m changes. In this study we experimentally determined the isosbestic point for a near infrared dye Di-4-ANBDQPQ. We then present a theoretical study examining the dye linear or non-linear response as the fractional fluorescence change of V_m change, due to emission spectra shift and amplitude change, over a range of excitation wavelengths and LP filters. Linear "optical" response is important to quantify certain aspects of cardiac dynamics such as the action potential (AP) shape and duration, especially when studying drug effects and dynamical substrates for arrhythmia development.

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

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

Validation and Trustworthiness of Multiscale Models of Cardiac Electrophysiology

Computational models of cardiac electrophysiology have a long history in basic science applications and device design and evaluation, but have significant potential for clinical applications in all areas of cardiovascular medicine, including functional imaging and mapping, drug safety evaluation, disease diagnosis, patient selection, and therapy optimisation or personalisation. For all stakeholders to be confident in model-based clinical decisions, cardiac electrophysiological (CEP) models must be demonstrated to be trustworthy and reliable. Credibility, that is, the belief in the predictive capability, of a computational model is primarily established by performing validation, in which model predictions are compared to experimental or clinical data. However, there are numerous challenges to performing validation for highly complex multi-scale physiological models such as CEP models. As a result, credibility of CEP model predictions is usually founded upon a wide range of distinct factors, including various types of validation results, underlying theory, evidence supporting model assumptions, evidence from model calibration, all at a variety of scales from ion channel to cell to organ. Consequently, it is often unclear, or a matter for debate, the extent to which a CEP model can be trusted for a given application. The aim of this article is to clarify potential rationale for the trustworthiness of CEP models by reviewing evidence that has been (or could be) presented to support their credibility. We specifically address the complexity and multi-scale nature of CEP models which makes traditional model evaluation difficult. In addition, we make explicit some of the credibility justification that we believe is implicitly embedded in the CEP modeling literature. Overall, we provide a fresh perspective to CEP model credibility, and build a depiction and categorisation of the wide-ranging body of credibility evidence for CEP models. This paper also represents a step toward the extension of model evaluation methodologies that are currently being developed by the medical device community, to physiological models.

Effects of mechanical feedback on the stability of cardiac scroll waves: A bidomain electro-mechanical simulation study

In this work, we investigate the influence of cardiac tissue deformation on re-entrant wave dynamics. We have developed a 3D strongly coupled electro-mechanical Bidomain model posed on an ideal monoventricular geometry, including fiber direction anisotropy and stretch-activated currents (SACs). The cardiac mechanical deformation influences the bioelectrical activity with two main mechanical feedback: (a) the geometric feedback (GEF) due to the presence of the deformation gradient in the diffusion coefficients and in a convective term depending on the deformation rate and (b) the mechano-electric feedback (MEF) due to SACs. Here, we investigate the relative contribution of these two factors with respect to scroll wave stability. We extend the previous works [Keldermann et al., Am. J. Physiol. Heart Circ. Physiol. 299, H134-H143 (2010) and Hu et al., PLoS One 8(4), e60287 (2013)] that were based on the Monodomain model and a simple non-selective linear SAC, while here we consider the full Bidomain model and both selective and non-selective components of SACs. Our simulation results show that the stability of cardiac scroll waves is influenced by MEF, which in case of low reversal potential of non-selective SACs might be responsible for the onset of ventricular fibrillation; GEF increases the scroll wave meandering but does not determine the scroll wave stability.

Determining the light scattering and absorption parameters from forward-directed flux measurements in cardiac tissue

We describe a method to accurately measure the light scattering model parameters from forward-directed flux (FDF) measurements carried out with a fiber-optic probe (optrode). Improved determination of light scattering parameters will, in turn, permit better modeling and interpretation of optical mapping in the heart using voltage-sensitive dyes. Using our optrode-based system, we carried out high spatial resolution measurements of FDF in intact and homogenized cardiac tissue, as well as in intralipid-based tissue phantoms. The samples were illuminated with a broad collimated beam at 660 and 532 nm. Measurements were performed with a plunge fiber-optic probe (NA=0.22) at a spatial resolution of up to 10  μm. In the vicinity of the illuminated surface, the FDF consistently manifested a fast decaying exponent with a space constant comparable with the decay rate of ballistic photons. Using a Monte Carlo model, we obtained a simple empirical formula linking the rate of the fast exponent to the scattering coefficient, the anisotropy parameter g, and the numerical aperture of the probe. The estimates of scattering coefficient based on this formula were validated in tissue phantoms. Potential applications of optical fiber-based FDF measurements for the evaluation of optical parameters in turbid media are discussed.

Voltage-Sensitive Fluorescence of Indocyanine Green in the Heart

So far, the optical mapping of cardiac electrical signals using voltage-sensitive fluorescent dyes has only been performed in experimental studies because these dyes are not yet approved for clinical use. It was recently reported that the well-known and widely used fluorescent dye indocyanine green (ICG), which has FDA approval, exhibits voltage sensitivity in various tissues, thus raising hopes that electrical activity could be optically mapped in the clinic. The aim of this study was to explore the possibility of using ICG to monitor cardiac electrical activity. Optical mapping experiments were performed on Langendorff rabbit hearts stained with ICG and perfused with electromechanical uncouplers. The residual contraction force and electrical action potentials were recorded simultaneously. Our research confirms that ICG is a voltage-sensitive dye with a dual-component (fast and slow) response to membrane potential changes. The fast component of the optical signal (OS) can have opposite polarities in different parts of the fluorescence spectrum. In contrast, the polarity of the slow component remains the same throughout the entire spectrum. Separating the OS into these components revealed two different voltage-sensitivity mechanisms for ICG. The fast component of the OS appears to be electrochromic in nature, whereas the slow component may arise from the redistribution of the dye molecules within or around the membrane. Both components quite accurately track the time of electrical signal propagation, but only the fast component is suitable for estimating the shape and duration of action potentials. Because ICG has voltage-sensitive properties in the entire heart, we suggest that it can be used to monitor cardiac electrical behavior in the clinic.

Optical Imaging of Cardiac Action Potential

This chapter reviews the major milestones and scientific achievements facilitated by optical imaging of the action potential in the heart over more than four decades since its introduction. We discuss the limitations of this technique, which sometimes are not fully recognized; the unresolved issues, such as motion artifacts, and the newest developments and future directions.

Evaluation of excitation propagation in the rabbit heart: optical mapping and transmural microelectrode recordings

Background

Because of the optical features of heart tissue, optical and electrical action potentials are only moderately associated, especially when near-infrared dyes are used in optical mapping (OM) studies.

Objective

By simultaneously recording transmural electrical action potentials (APs) and optical action potentials (OAPs), we aimed to evaluate the contributions of both electrical and optical influences to the shape of the OAP upstroke.

Methods and Results

A standard glass microelectrode and OM, using an near-infrared fluorescent dye (di-4-ANBDQBS), were used to simultaneously record transmural APs and OAPs in a Langendorff-perfused rabbit heart during atrial, endocardial, and epicardial pacing. The actual profile of the transmural AP upstroke across the LV wall, together with the OAP upstroke, allowed for calculations of the probing-depth constant (k ~2.1 mm, n = 24) of the fluorescence measurements. In addition, the transmural AP recordings aided the quantitative evaluation of the influences of depth-weighted and lateral-scattering components on the OAP upstroke. These components correspond to the components of the propagating electrical wave that are transmural and parallel to the epicardium. The calculated mean values for the depth-weighted and lateral-scattering components, whose sum comprises the OAP upstroke, were (in ms) 10.18 ± 0.62 and 0.0 ± 0.56 for atrial stimulation, 9.37 ± 1.12 and 3.01 ± 1.30 for endocardial stimulation, and 6.09 ± 0.79 and 8.16 ± 0.98 for epicardial stimulation; (n = 8 for each). For this dye, 90% of the collected fluorescence originated up to 4.83 ± 0.18 mm (n = 24) from the epicardium.

Conclusions

The co-registration of OM and transmural microelectrode APs enabled the probing depth of fluorescence measurements to be calculated and the OAP upstroke to be divided into two components (depth-weighted and lateral-scattering), and it also allowed the relative strengths of their effects on the shape of the OAP upstroke to be evaluated.

Evolution of action potential alternans in rabbit heart during acute regional ischemia

This study investigates the development of the spatiotemporal pattern of action potential alternans during acute regional ischemia. Experiments were carried out in isolated Langendorff-perfused rabbit heart using a combination of optical mapping and microelectrode recordings. The alternans pattern significantly changed over time and had a biphasic character reaching maximum at 6–9 min after occlusion. Phase I (3–11 minutes of ischemia) is characterized by rapid increase in the alternans magnitude and expansion of the alternans territory. Phase I is followed by gradual decline of alternans (Phase II) in both magnitude and territory. During both phases we observed significant beat-to-beat variations of the optical action potential amplitude (OAPA) alternans. Simultaneous microelectrode recordings from subepicardial and subendocardial layers showed that OAPA alternans coincided with intramural 2 : 1 conduction blocks. Our findings are consistent with the modeling studies predicting that during acute regional ischemia alternans can be driven by 2 : 1 conduction blocks in the ischemic region.

Simulating photon scattering effects in structurally detailed ventricular models using a Monte Carlo approach

Light scattering during optical imaging of electrical activation within the heart is known to significantly distort the optically-recorded action potential (AP) upstroke, as well as affecting the magnitude of the measured response of ventricular tissue to strong electric shocks. Modeling approaches based on the photon diffusion equation have recently been instrumental in quantifying and helping to understand the origin of the resulting distortion. However, they are unable to faithfully represent regions of non-scattering media, such as small cavities within the myocardium which are filled with perfusate during experiments. Stochastic Monte Carlo (MC) approaches allow simulation and tracking of individual photon “packets” as they propagate through tissue with differing scattering properties. Here, we present a novel application of the MC method of photon scattering simulation, applied for the first time to the simulation of cardiac optical mapping signals within unstructured, tetrahedral, finite element computational ventricular models. The method faithfully allows simulation of optical signals over highly-detailed, anatomically-complex MR-based models, including representations of fine-scale anatomy and intramural cavities. We show that optical action potential upstroke is prolonged close to large subepicardial vessels than further away from vessels, at times having a distinct “humped” morphology. Furthermore, we uncover a novel mechanism by which photon scattering effects around vessels cavities interact with “virtual-electrode” regions of strong de-/hyper-polarized tissue surrounding cavities during shocks, significantly reducing the apparent optically-measured epicardial polarization. We therefore demonstrate the importance of this novel optical mapping simulation approach along with highly anatomically-detailed models to fully investigate electrophysiological phenomena driven by fine-scale structural heterogeneity.

Toward microendoscopy-inspired cardiac optogenetics in vivo: technical overview and perspective

The ability to perform precise, spatially localized actuation and measurements of electrical activity in the heart is crucial in understanding cardiac electrophysiology and devising new therapeutic solutions for control of cardiac arrhythmias. Current cardiac imaging techniques (i.e. optical mapping) employ voltage- or calcium-sensitive fluorescent dyes to visualize the electrical signal propagation through cardiac syncytium in vitro or in situ with very high-spatiotemporal resolution. The extension of optogenetics into the cardiac field, where cardiac tissue is genetically altered to express light-sensitive ion channels allowing electrical activity to be elicited or suppressed in a precise cell-specific way, has opened the possibility for all-optical interrogation of cardiac electrophysiology. In vivo application of cardiac optogenetics faces multiple challenges and necessitates suitable optical systems employing fiber optics to actuate and sense electrical signals. In this technical perspective, we present a compendium of clinically relevant access routes to different parts of the cardiac electrical conduction system based on currently employed catheter imaging systems and determine the quantitative size constraints for endoscopic cardiac optogenetics. We discuss the relevant technical advancements in microendoscopy, cardiac imaging, and optogenetics and outline the strategies for combining them to create a portable, miniaturized fiber-based system for all-optical interrogation of cardiac electrophysiology in vivo.

Transmembrane current imaging in the heart during pacing and fibrillation

Recently, we described a method to quantify the time course of total transmembrane current (Im) and the relative role of its two components, a capacitive current (Ic) and a resistive current (Iion), corresponding to the cardiac action potential during stable propagation. That approach involved recording high-fidelity (200 kHz) transmembrane potential (Vm) signals with glass microelectrodes at one site using a spatiotemporal coordinate transformation via measured conduction velocity. Here we extend our method to compute these transmembrane currents during stable and unstable propagation from fluorescence signals of Vm at thousands of sites (3 kHz), thereby introducing transmembrane current imaging. In contrast to commonly used linear Laplacians of extracellular potential (Ve) to compute Im, we utilized nonlinear image processing to compute the required second spatial derivatives of Vm. We quantified the dynamic spatial patterns of current density of Im and Iion for both depolarization and repolarization during pacing (including nonplanar patterns) by calibrating data with the microelectrode signals. Compared to planar propagation, we found that the magnitude of Iion was significantly reduced at sites of wave collision during depolarization but not repolarization. Finally, we present uncalibrated dynamic patterns of Im during ventricular fibrillation and show that Im at singularity sites was monophasic and positive with a significant nonzero charge (Im integrated over 10 ms) in contrast with nonsingularity sites. Our approach should greatly enhance the understanding of the relative roles of functional (e.g., rate-dependent membrane dynamics and propagation patterns) and static spatial heterogeneities (e.g., spatial differences in tissue resistance) via recordings during normal and compromised propagation, including arrhythmias.

Diastolic field stimulation: the role of shock duration in epicardial activation and propagation

Detailed knowledge of tissue response to both systolic and diastolic shock is critical for understanding defibrillation. Diastolic field stimulation has been much less studied than systolic stimulation, particularly regarding transient virtual anodes. Here we investigated high-voltage-induced polarization and activation patterns in response to strong diastolic shocks of various durations and of both polarities, and tested the hypothesis that the activation versus shock duration curve contains a local minimum for moderate shock durations, and it grows for short and long durations. We found that 0.1-0.2-ms shocks produced slow and heterogeneous activation. During 0.8-1 ms shocks, the activation was very fast and homogeneous. Further shock extension to 8 ms delayed activation from 1.55 ± 0.27 ms and 1.63 ± 0.21 ms at 0.8 ms shock to 2.32 ± 0.41 ms and 2.37 ± 0.3 ms (N = 7) for normal and opposite polarities, respectively. The traces from hyperpolarized regions during 3-8 ms shocks exhibited four different phases: beginning negative polarization, fast depolarization, slow depolarization, and after-shock increase in upstroke velocity. Thus, the shocks of >3 ms in duration created strong hyperpolarization associated with significant delay (P < 0.05) in activation compared with moderate shocks of 0.8 and 1 ms. This effect appears as a dip in the activation-versus-shock-duration curve.

Combining wet and dry research: experience with model development for cardiac mechano-electric structure-function studies

Since the development of the first mathematical cardiac cell model 50 years ago, computational modelling has become an increasingly powerful tool for the analysis of data and for the integration of information related to complex cardiac behaviour. Current models build on decades of iteration between experiment and theory, representing a collective understanding of cardiac function. All models, whether computational, experimental, or conceptual, are simplified representations of reality and, like tools in a toolbox, suitable for specific applications. Their range of applicability can be explored (and expanded) by iterative combination of 'wet' and 'dry' investigation, where experimental or clinical data are used to first build and then validate computational models (allowing integration of previous findings, quantitative assessment of conceptual models, and projection across relevant spatial and temporal scales), while computational simulations are utilized for plausibility assessment, hypotheses-generation, and prediction (thereby defining further experimental research targets). When implemented effectively, this combined wet/dry research approach can support the development of a more complete and cohesive understanding of integrated biological function. This review illustrates the utility of such an approach, based on recent examples of multi-scale studies of cardiac structure and mechano-electric function.

Optical mapping at increased illumination intensities

Voltage-sensitive fluorescent dyes have become a major tool in cardiac and neuro-electrophysiology. Achieving high signal-to-noise ratios requires increased illumination intensities, which may cause photobleaching and phototoxicity. The optimal range of illumination intensities varies for different dyes and must be evaluated individually. We evaluate two dyes: di-4-ANBDQBS (excitation 660 nm) and di-4-ANEPPS (excitation 532 nm) in the guinea pig heart. The light intensity varies from 0.1 to 5  mW/mm2, with the upper limit at 5 to 10 times above values reported in the literature. The duration of illumination was 60 s, which in guinea pigs corresponds to 300 beats at a normal heart rate. Within the identified duration and intensity range, neither dye shows significant photobleaching or detectable phototoxic effects. However, light absorption at higher intensities causes noticeable tissue heating, which affects the electrophysiological parameters. The most pronounced effect is a shortening of the action potential duration, which, in the case of 532-nm excitation, can reach ∼30%. At 660-nm excitation, the effect is ∼10%. These findings may have important implications for the design of optical mapping protocols in biomedical applications.

Bridging experiments, models and simulations: an integrative approach to validation in computational cardiac electrophysiology

Computational models in physiology often integrate functional and structural information from a large range of spatiotemporal scales from the ionic to the whole organ level. Their sophistication raises both expectations and skepticism concerning how computational methods can improve our understanding of living organisms and also how they can reduce, replace, and refine animal experiments. A fundamental requirement to fulfill these expectations and achieve the full potential of computational physiology is a clear understanding of what models represent and how they can be validated. The present study aims at informing strategies for validation by elucidating the complex interrelations among experiments, models, and simulations in cardiac electrophysiology. We describe the processes, data, and knowledge involved in the construction of whole ventricular multiscale models of cardiac electrophysiology. Our analysis reveals that models, simulations, and experiments are intertwined, in an assemblage that is a system itself, namely the model-simulation-experiment (MSE) system. We argue that validation is part of the whole MSE system and is contingent upon 1) understanding and coping with sources of biovariability; 2) testing and developing robust techniques and tools as a prerequisite to conducting physiological investigations; 3) defining and adopting standards to facilitate the interoperability of experiments, models, and simulations; 4) and understanding physiological validation as an iterative process that contributes to defining the specific aspects of cardiac electrophysiology the MSE system targets, rather than being only an external test, and that this is driven by advances in experimental and computational methods and the combination of both.

Optical imaging of voltage and calcium in cardiac cells & tissues

Cardiac optical mapping has proven to be a powerful technology for studying cardiovascular function and disease. The development and scientific impact of this methodology are well-documented. Because of its relevance in cardiac research, this imaging technology advances at a rapid pace. Here, we review technological and scientific developments during the past several years and look toward the future. First, we explore key components of a modern optical mapping set-up, focusing on: (1) new camera technologies; (2) powerful light-emitting-diodes (from ultraviolet to red) for illumination; (3) improved optical filter technology; (4) new synthetic and optogenetic fluorescent probes; (5) optical mapping with motion and contraction; (6) new multiparametric optical mapping techniques; and (7) photon scattering effects in thick tissue preparations. We then look at recent optical mapping studies in single cells, cardiomyocyte monolayers, atria, and whole hearts. Finally, we briefly look into the possible future roles of optical mapping in the development of regenerative cardiac research, cardiac cell therapies, and molecular genetic advances.

Extracting surface activation time from the optically recorded action potential in three-dimensional myocardium

Optical mapping has become an indispensible tool for studying cardiac electrical activity. However, due to the three-dimensional nature of the optical signal, the optical upstroke is significantly longer than the electrical upstroke. This raises the issue of how to accurately determine the activation time on the epicardial surface. The purpose of this study was to establish a link between the optical upstroke and exact surface activation time using computer simulations, with subsequent validation by a combination of microelectrode recordings and optical mapping experiments. To simulate wave propagation and associated optical signals, we used a hybrid electro-optical model. We found that the time of the surface electrical activation (t(E)) within the accuracy of our simulations coincided with the maximal slope of the optical upstroke (t(F)*) for a broad range of optical attenuation lengths. This was not the case when the activation time was determined at 50% amplitude (t(F50)) of the optical upstroke. The validation experiments were conducted in isolated Langendorff-perfused rat hearts and coronary-perfused pig left ventricles stained with either di-4-ANEPPS or the near-infrared dye di-4-ANBDQBS. We found that t(F)* was a more accurate measure of t(E) than was t(F50) in all experimental settings tested (P = 0.0002). Using t(F)* instead of t(F50) produced the most significant improvement in measurements of the conduction anisotropy and the transmural conduction time in pig ventricles.

Single-detector simultaneous optical mapping of V(m) and [Ca(2+)](i) in cardiac monolayers

Simultaneous mapping of transmembrane voltage (V(m)) and intracellular Ca(2+) concentration (Ca(i)) has been used for studies of normal and abnormal impulse propagation in cardiac tissues. Existing dual mapping systems typically utilize one excitation and two emission bandwidths, requiring two photodetectors with precise pixel registration. In this study we describe a novel, single-detector mapping system that utilizes two excitation and one emission band for the simultaneous recording of action potentials and calcium transients in monolayers of neonatal rat cardiomyocytes. Cells stained with the Ca(2+)-sensitive dye X-Rhod-1 and the voltage-sensitive dye Di-4-ANEPPS were illuminated by a programmable, multicolor LED matrix. Blue and green LED pulses were flashed 180° out of phase at a rate of 488.3 Hz using a custom-built dual bandpass excitation filter that transmitted blue (482 ± 6 nm) and green (577 ± 31 nm) light. A long-pass emission filter (>605 nm) and a 504-channel photodiode array were used to record combined signals from cardiomyocytes. Green excitation yielded Ca(i) transients without significant crosstalk from V(m). Crosstalk present in V(m) signals obtained with blue excitation was removed by subtracting an appropriately scaled version of the Ca(i) transient. This method was applied to study delay between onsets of action potentials and Ca(i) transients in anisotropic cardiac monolayers.

Imaging electrical excitation inside the myocardial wall

Cardiac arrhythmias are often triggered by ectopic membrane depolarization originating deep inside the myocardial wall. Here we propose a new method utilizing a novel near-infrared voltage-sensitive fluorescent dye DI-4-ANBDQBS to determine the three-dimensional (3D) coordinates of the sources of such depolarization. We tested the method in live preparations of pig left and right ventricular myocardium (thickness 8-18 mm) and phantoms imitating the optical properties of myocardial tissue. The method utilizes an alternating transillumination approach that involves comparing pairs of simultaneously recorded broad-field epifluorescence and transillumination images produced at two alternating directions of illumination. Recordings were taken simultaneously by two CCD cameras facing the endocardial and epicardial surfaces of the heart at a frame rate up to 3 KHz. In live preparations, we were able to localize the origin of the depolarization wave with a precision of ±1.3mm in the transmural direction and 3 mm in the image plane. The accuracy of detection was independent of the depth of the source inside ventricular wall.

Construction and validation of anisotropic and orthotropic ventricular geometries for quantitative predictive cardiac electrophysiology

Reaction-diffusion computational models of cardiac electrophysiology require both dynamic excitation models that reconstruct the action potentials of myocytes as well as datasets of cardiac geometry and architecture that provide the electrical diffusion tensor D, which determines how excitation spreads through the tissue. We illustrate an experimental pipeline we have developed in our laboratories for constructing and validating such datasets. The tensor D changes with location in the myocardium, and is determined by tissue architecture. Diffusion tensor magnetic resonance imaging (DT-MRI) provides three eigenvectors e(i) and eigenvalues λ(i) at each voxel throughout the tissue that can be used to reconstruct this architecture. The primary eigenvector e(1) is a histologically validated measure of myocyte orientation (responsible for anisotropic propagation). The secondary and tertiary eigenvectors (e(2) and e(3)) specify the directions of any orthotropic structure if λ(2) is significantly greater than λ(3)-this orthotropy has been identified with sheets or cleavage planes. For simulations, the components of D are scaled in the fibre and cross-fibre directions for anisotropic simulations (or fibre, sheet and sheet normal directions for orthotropic tissues) so that simulated conduction velocities match values from optical imaging or plunge electrode experiments. The simulated pattern of propagation of action potentials in the models is partially validated by optical recordings of spatio-temporal activity on the surfaces of hearts. We also describe several techniques that enhance components of the pipeline, or that allow the pipeline to be applied to different areas of research: Q ball imaging provides evidence for multi-modal orientation distributions within a fraction of voxels, infarcts can be identified by changes in the anisotropic structure-irregularity in myocyte orientation and a decrease in fractional anisotropy, clinical imaging provides human ventricular geometry and can identify ischaemic and infarcted regions, and simulations in human geometries examine the roles of anisotropic and orthotropic architecture in the initiation of arrhythmias.

High-precision recording of the action potential in isolated cardiomyocytes using the near-infrared fluorescent dye di-4-ANBDQBS

The use of voltage-sensitive fluorescent dyes (VSD) for noninvasive measurement of the action potential (AP) in isolated cells has been hindered by low-photon yield of the preparation, dye toxicity, and photodynamic damage. Here we used a new red-shifted VSD, di-4-ANBDQBS, and a fast electron-multiplied charge-coupled device camera for optical AP (OAP) recording in guinea pig cardiac myocytes. Loading di-4-ANBDQBS did not alter APs recorded with micropipette. With short laser exposures (just enough to record one OAP every 1-5 min), di-4-ANBDQBS yielded fluorescent signals with very high signal-to-background ratios (change in fluorescence on depolarization/fluorescence at resting potential: 19.2 ± 4.1%) and signal-to-noise ratios (40 ± 13.2). Quantum chemical calculations comparing the ANBDQ chromophore to the conventional ANEP chromophore showed that the higher wavelength and the greater voltage sensitivity of the former have the same electro-optical origin: a longer path for electron redistribution in the excited state. OAP closely tracked simultaneously recorded electrical APs, permitting measurement of AP duration within 1% error. Prolonged laser exposure caused progressive AP duration prolongation and instability. However, these effects were alleviated or abolished by reducing the dye concentration and by perfusion with antioxidants. Thus the presented technique provides a unique opportunity for noninvasive AP recording in single cardiomyocytes.

Reconstructing subsurface electrical wave orientation from cardiac epi-fluorescence recordings: Monte Carlo versus diffusion approximation

The development of voltage-sensitive dyes has revolutionized cardiac electrophysiology and made optical imaging of cardiac electrical activity possible. Photon diffusion models coupled to electrical excitation models have been successful in qualitatively predicting the shape of the optical action potential and its dependence on subsurface electrical wave orientation. However, the accuracy of the diffusion equation in the visible range, especially for thin tissue preparations, remains unclear. Here, we compare diffusion and Monte Carlo (MC) based models and we investigate the role of tissue thickness. All computational results are compared to experimental data obtained from intact guinea pig hearts. We show that the subsurface volume contributing to the epi-fluorescence signal extends deeper in the tissue when using MC models, resulting in longer optical upstroke durations which are in better agreement with experiments. The optical upstroke morphology, however, strongly correlates to the subsurface propagation direction independent of the model and is consistent with our experimental observations.

Calculation of optical signal using three-dimensional bidomain/diffusion model reveals distortion of the transmembrane potential

Optical mapping experiments allow investigators to view the effects of electrical currents on the transmembrane potential, V(m), as a shock is applied to the heart. One important consideration is whether the optical signal accurately represents V(m). We have combined the bidomain equations along with the photon diffusion equation to study the excitation and emission of photons during optical mapping of cardiac tissue. Our results show that this bidomain/diffusion model predicts an optical signal that is much smaller than V(m) near a stimulating electrode, a result consistent with experimental observations. Yet, this model, which incorporates the effect of lateral averaging, also reveals an optical signal that overestimates V(m) at distances >1 mm away from the electrode. Although V(m) falls off with distance r from the electrode as exp(-r/lambda)/r, the optical signal decays as a simple exponential, exp(-r/lambda). Moreover, regions of hyperpolarization adjacent to a cathode are emphasized in the optical signal compared to the region of depolarization under the cathode. Imaging methods utilizing optical mapping techniques will need to account for these distortions to accurately reconstruct V(m).

Extracting intramural wavefront orientation from optical upstroke shapes in whole hearts

Information about intramural propagation of electrical excitation is crucial to understanding arrhythmia mechanisms in thick ventricular muscle. There is currently a controversy over whether it is possible to extract such information from the shape of the upstroke in optical mapping recordings. We show that even in the complex geometry of a whole guinea pig heart, optical upstroke morphology reveals the 3D wavefront orientation near the surface. To characterize the upstroke morphology, we use V(F)(*), the fractional level at which voltage-sensitive fluorescence, V(F), has maximal time derivative. Low values of V(F)(*)( approximately 0.2) indicate a wavefront moving away from the surface, high values of V(F)(*) ( approximately 0.6) a wavefront moving toward the surface, and intermediate values of V(F)(*) ( approximately 0.4) a wavefront moving parallel to the surface. We further performed computer simulations using Luo-Rudy II electrophysiology and a simplified 3D geometry. The simulated V(F)(*) maps for free wall and apical stimulations as well as for sinus rhythm are in good quantitative agreement with the averaged experimental results. Furthermore, computer simulations show that the effect of the curvature of the heart on wave propagation is negligible.

Intra-myocardial cusp waves and their manifestation in optical mapping signals

The rotating fiber orientation within the cardiac wall substantially affects the electrical propagation and can cause intra-myocardial cusp waves. Numerical simulations have shown that the cusps form in layers where propagation is perpendicular to the fiber orientation and lead to complex wave front morphologies. They can travel across layers and break through at the epi- or endocardial surfaces where they cause apparent accelerations of propagation. The validation of these results remains a major experimental challenge. In the present study, we investigate both computationally and experimentally how intramural cusp waves can be detected using optical imaging. Our simulations show that cusps alter the optical upstroke morphology and can be detected well before they reach the surface (up to 1 mm deep). Experiments in Langendorff-perfused guinea pig hearts are consistent with our numerical findings.

Modulation of shock-end virtual electrode polarisation as a direct result of 3D fluorescent photon scattering

Due to the large transmural variation in transmembrane potential following the application of strong electric shocks, it is thought that fluorescent photon scattering from depth plays a significant role in optical signal modulation at shock-end. For the first time, a model of photon scattering is used to accurately synthesize fluorescent signals over the irregular geometry of the rabbit ventricles following the application of such strong shocks. A bidomain representation of electrical activity is combined with finite element solutions to the photon diffusion equation, simulating both the excitation and emission processes, over an anatomically-based model of rabbit ventricular geometry and fiber orientation. Photon scattering from within a 3D volume beneath the epicardial optical recording site is shown to transduce differences in transmembrane potential within this volume through the myocardial wall. This leads directly to a significantly modulated optical signal response with respect to that predicted by the bidomain simulations, distorting epicardial virtual electrode polarization produced at shock-end. Furthermore, we show that this degree of distortion is very sensitive to the optical properties of the tissue, an important variable to consider during experimental mapping set-ups. These findings provide an essential first-step in aiding the interpretation of experimental optical mapping recordings following strong defibrillation shocks.

Photon scattering effects in optical mapping of propagation and arrhythmogenesis in the heart

BACKGROUND

Optical mapping is a widely used experimental tool providing high-resolution recordings of cardiac electrical activity. However, the technique is limited by signal distortion due to photon scattering in the tissue. Computational models of the fluorescence recording are capable of assessing these distortion effects, providing important insight to assist experimental data interpretation.

METHODS

We present results from a new panoramic optical mapping model, which is used to assess distortion in ventricular optical mapping signals during pacing and arrhythmogenesis arising from 3-dimensional photon scattering.

RESULTS/CONCLUSIONS

We demonstrate that accurate consideration of wavefront propagation within the complex ventricular structure, along with accurate representation of photon scattering in 3 dimensions, is essential to faithfully assess distortion effects arising during optical mapping. In this article, examined effects include (1) the specific morphology of the optical action potential upstroke during pacing and (2) the shift in the location of epicardial phase singularities obtained from fluorescent maps.

The role of photon scattering in optical signal distortion during arrhythmia and defibrillation

Optical mapping of arrhythmias and defibrillation provides important insights; however, a limitation of the technique is signal distortion due to photon scattering. The goal of this experimental/simulation study is to investigate the role of three-dimensional photon scattering in optical signal distortion during ventricular tachycardia (VT) and defibrillation. A three-dimensional realistic bidomain rabbit ventricular model was combined with a model of photon transport. Shocks were applied via external electrodes to induce sustained VT, and transmembrane potentials (V(m)) were compared with synthesized optical signals (V(opt)). Fluorescent recordings were conducted in isolated rabbit hearts to validate simulation results. Results demonstrate that shock-induced membrane polarization magnitude is smaller in V(opt) and in experimental signals as compared to V(m). This is due to transduction of potentials from weakly polarized midmyocardium to the epicardium. During shock-induced reentry and in sustained VT, photon scattering, combined with complex wavefront dynamics, results in optical action potentials near a filament exhibiting i), elevated resting potential, ii), reduced amplitude relative to pacing, and iii), dual-humped morphologies. A shift of up to 4 mm in the phase singularity location was observed in V(opt) maps when compared to V(m). This combined experimental/simulation study provides an interpretation of optical recordings during VT and defibrillation.

Depth-resolved optical imaging of transmural electrical propagation in perfused heart

We present a study of the 3-dimensional (3D) propagation of electrical waves in the heart wall using Laminar Optical Tomography (LOT). Optical imaging contrast is provided by a voltage sensitive dye whose fluorescence reports changes in membrane potential. We examined the transmural propagation dynamics of electrical waves in the right ventricle of Langendorf perfused rat hearts, initiated either by endo-cardial or epi-cardial pacing. 3D images were acquired at an effective frame rate of 667Hz. We compare our experimental results to a mathematical model of electrical transmural propagation. We demonstrate that LOT can clearly resolve the direction of propagation of electrical waves within the cardiac wall, and that the dynamics observed agree well with the model of electrical propagation in rat ventricular tissue.

Spatiotemporal relationship between intracellular Ca2+ dynamics and wave fragmentation during ventricular fibrillation in isolated blood-perfused pig hearts

Normal "master-slave" relationship between the action potential (AP) and intracellular Ca2+ transient (Ca(i)T) is sometimes altered during ventricular fibrillation (VF). The nature of AP/Ca(i)T dissociation during VF and its role in inducing wavebreaks (WBs) remain unclear. We simultaneously mapped AP (RH237) and Ca(i)T (Rhod-2) during VF in blood-perfused pig hearts. We computed AP and Ca(i)T dominant frequency (DF) and Ca(i)T delay in each AP cycle. We identified WBs as singularity points in AP phase movies and sites of conduction block (CB) as sites where an AP wavefront failed to propagate. We analyzed spatiotemporal relationship between abnormal AP/Ca(i)T sequences and CB sites. We used a calcium chelator (BAPTA-AM) to abolish Ca(i)T and test its involvement in WB formation. During VF, the DF difference between AP and Ca(i)T was <10% of the respective values in 95% of pixels, and 80% of all Ca(i)T upstrokes occurred during the initial 25% of the excitation cycle. Aberrant sequences of AP and Ca(i)T occurred almost exclusively near CB sites but could be traced to normal wavefront sequences away from CB sites. Thus, apparent AP/Ca(i)T dissociation was largely attributable to spatial uncertainty of the absolute position of block of each wave. BAPTA-AM reduced Ca(i)T amplitude to 30.5+/-12.9% of control and the DF of AP from 12.2+/-1.6 to 10.4+/-1.3 Hz (P<0.01), but did not significantly alter WB incidence (0.76+/-0.19 versus 0.72+/-0.19 SP/mm2). These results do not support presence of spontaneous, non-voltage-gated Ca(i)Ts during VF and suggest that AP/Ca(i)T dissociation is a consequence rather than a cause of wave fragmentation.

A guide to modelling cardiac electrical activity in anatomically detailed ventricles

One of the most recent trends in cardiac electrophysiology is the development of integrative anatomically accurate models of the heart, which include description of cardiac activity from sub-cellular and cellular level to the level of the whole organ. In order to construct this type of model, a researcher needs to collect a wide range of information from books and journal articles on various aspects of biology, physiology, electrophysiology, numerical mathematics and computer programming. The aim of this methodological article is to survey recent developments in integrative modelling of electrical activity in the ventricles of the heart, and to provide a practical guide to the resources and tools that are available for work in this exciting and challenging area.

Near-infrared voltage-sensitive fluorescent dyes optimized for optical mapping in blood-perfused myocardium

BACKGROUND

Styryl voltage-sensitive dyes (e.g., di-4-ANEPPS) have been used successfully for optical mapping in cardiac cells and tissues. However, their utility for probing electrical activity deep inside the myocardial wall and in blood-perfused myocardium has been limited because of light scattering and high absorption by endogenous chromophores and hemoglobin at blue-green excitation wavelengths.

OBJECTIVE

The purpose of this study was to characterize two new styryl dyes--di-4-ANBDQPQ (JPW-6003) and di-4-ANBDQBS (JPW-6033)--optimized for blood-perfused tissue and intramural optical mapping.

METHODS

Voltage-dependent spectra were recorded in a model lipid bilayer. Optical mapping experiments were conducted in four species (mouse, rat, guinea pig, and pig). Hearts were Langendorff perfused using Tyrode's solution and blood (pig). Dyes were loaded via bolus injection into perfusate. Transillumination experiments were conducted in isolated coronary-perfused pig right ventricular wall preparations.

RESULTS

The optimal excitation wavelength in cardiac tissues (650 nm) was >70 nm beyond the absorption maximum of hemoglobin. Voltage sensitivity of both dyes was approximately 10% to 20%. Signal decay half-life due to dye internalization was 80 to 210 minutes, which is 5 to 7 times slower than for di-4-ANEPPS. In transillumination mode, DeltaF/F was as high as 20%. In blood-perfused tissues, DeltaF/F reached 5.5% (1.8 times higher than for di-4-ANEPPS).

CONCLUSION

We have synthesized and characterized two new near-infrared dyes with excitation/emission wavelengths shifted >100 nm to the red. They provide both high voltage sensitivity and 5 to 7 times slower internalization rate compared to conventional dyes. The dyes are optimized for deeper tissue probing and optical mapping of blood-perfused tissue, but they also can be used for conventional applications.

Illumination and fluorescence collection volumes for fiber optic probes in tissue

Optical fibers can deliver light to, and collect it from, regions deep in tissue. However, reported illumination and fluorescence collection volumes adjacent to the fiber tip have been inconsistent, and systematic data on this topic are not available. Illumination and fluorescence collection profiles were characterized with high spatial resolution for different optical fibers in tissue and various fluids using two-photon flash photolysis and excitation. We confirm that illumination and fluorescence collection volumes for optical fibers are near identical. Collection volume is determined by the core dimensions and numerical aperture (NA) of the fiber and the scattering properties of the medium. For a multimode optical fiber with 100 microm core diam and NA=0.22, 80% of the total fluorescence is collected from a depth of 170 microm in tissue and 465 microm in nonscattering fluid. A semiempirical mathematical description of photon flux adjacent to the fiber tip was also developed and validated. This was used to quantify the extent of temporal blurring associated with propagation of a wavefront of altered fluorescence emission across the region addressed by fiber optic probes. We provide information that will facilitate the design of optical probes for tissue imaging or therapeutic applications.

Lenses and effective spatial resolution in macroscopic optical mapping

Optical mapping of excitation dynamically tracks electrical waves travelling through cardiac or brain tissue by the use of fluorescent dyes. There are several characteristics that set optical mapping apart from other imaging modalities: dynamically changing signals requiring short exposure times, dim fluorescence demanding sensitive sensors and wide fields of view (low magnification) resulting in poor optical performance. These conditions necessitate the use of optics with good light gathering ability, i.e. lenses having high numerical aperture. Previous optical mapping studies often used sensor resolution to estimate the minimum spatial feature resolvable, assuming perfect optics and infinite contrast. We examine here the influence of finite contrast and real optics on the effective spatial resolution in optical mapping under broad-field illumination for both lateral (in-plane) resolution and axial (depth) resolution of collected fluorescence signals.

Detection of intramyocardial scroll waves using absorptive transillumination imaging

Optical imaging using voltage-sensitive dyes has become an important tool for studying vortex-like electrical waves in the heart. Such waves, known as spiral or scroll waves, can spontaneously form in pathological ventricular myocardium, causing ventricular fibrillation and sudden death. Until recently, observations of scroll waves were limited to their surface manifestations, thus providing little information about the shape and location of their organizing center, the filament. We use computer modeling to assess the feasibility of visualizing filaments using dynamic transillumination imaging in conjunction with near-IR voltage-sensitive absorptive dyes (absorptive transillumination). We simulate transillumination signals produced by the intramural scroll waves in a realistic slab of ventricular tissue with trabeculated endocardial surface. The computations use a detailed ionic model of electrical excitation (LRd) coupled to a photon transport model for cardiac tissue. Our simulations show that dynamic absorptive transillumination data, with subsequent processing involving either amplitude maps, time-space plots, or power-of-the-dominant-frequency maps, can be used to reliably detect intramural scroll waves through the whole thickness (approximately 10 mm) of the ventricular wall. Neither variations in the thickness of the myocardial wall nor noise impeded the detection of intramural filaments.

What can we learn from the optically recorded epicardial action potential?

Optical mapping using voltage-sensitive fluorescent dyes has become a major tool for studying excitation propagation in the heart. Computational and experimental studies have indicated that the optical upstroke morphology reflects the orientation of the subsurface excitation front. In a recent whole heart computational study performed by Bishop et al. (Bishop, M. J., B. Rodriguez, J. Eason, J. P. Whiteley, N. Trayanova, and D. J. Gavaghan. 2006. Synthesis of voltage-sensitive optical signals: application to panoramic optical mapping. Biophys. J. 90:2938-2945), an example was provided of two different directions of propagation having nevertheless very similar epicardial optical upstrokes. The goal of this comment is to clarify the interpretation of optical upstroke morphologies and reconcile the results obtained by Bishop et al. with previous computational and experimental studies.

Method for the three-dimensional localization of intramyocardial excitation centers using optical imaging

This study explores the possibility of localizing the excitation centers of electrical waves inside the heart wall using voltage-sensitive dyes (fluorescent or absorptive). In the present study, we propose a method for the 3-D localization of excitation centers from pairs of 2-D images obtained in two modes of observation: reflection and transillumination. Such images can be obtained using high-speed charge-coupled device (CCD) cameras and photodiode arrays with time resolution up to 0.5 ms. To test the method, we simulate optical signals produced by point sources and propagating ellipsoidal waves in 1-cm-thick slabs of myocardial tissue. Solutions of the optical diffusion equation are constructed by employing the method of images with Robin boundary conditions. The coordinates of point sources as well as of the centers of expanding waves can be accurately determined using the proposed algorithm. The method can be extended to depth estimations of the outer boundaries of the expanding wave. The depth estimates are based on ratios of spatially integrated images. The method shows high tolerance to noise and can give accurate results even at relatively low signal-to-noise ratios. In conclusion, we propose a novel and efficient algorithm for the localization of excitation centers in 3-D cardiac tissue.

Characterization of the relationship between preshock state and virtual electrode polarization-induced propagated graded responses resulting in arrhythmia induction

BACKGROUND

Studies have demonstrated that failed defibrillation shocks often are followed by an electrically quiescent period (isoelectric window); however, the underlying mechanisms remain incompletely understood. We recently suggested a new mechanism termed "virtual electrode polarization-induced propagated graded responses" (VEPiPGRs) that might play a role in the origin of the global postshock activation following the isoelectric window.

OBJECTIVES

The purpose of this study to elucidate the circumstances under which VEPiPGR activations originate for shocks given to paced right ventricular preparations. Specifically, we examined the dependence of VEPiPGRs on coupling interval (CI) and shock polarity and whether VEPiPGRs emerge preferentially on the epicardium or the endocardium.

METHODS

Simultaneous endocardial and epicardial activity in isolated right ventricular preparations (n = 4) was imaged optically following shocks of strength +/-5A. All VEPiPGRs were analyzed, and the time T from shock end to activation onset was recorded (isoelectric window is the smallest T among activations that propagated globally).

RESULTS

VEPiPGR activations occurred for CIs in the range from 80 to 150 ms. Average duration of T was 64.5 +/- 18.15 ms, with T decreasing as CI increased (Tmax = 82 ms, Tmin = 46 ms, linear-fit slope = -0.675). The average earliest CI at which cathodal (+5A) shocks resulted in VEPiPGRs was 87 ms compared with 116 ms for anodal (-5A) shocks. All VEPiPGR activations emerged first on the epicardium in a focal pattern, and all induced ventricular fibrillation.

CONCLUSION

The global activation that terminates the isoelectric window could result from VEPiPGRs that find an exit pathway. VEPiPGRs originate at the sites of maximum action potential abbreviation by the shock, always on the epicardium for the preparation used here.

Synthesis of voltage-sensitive optical signals: application to panoramic optical mapping

Fluorescent photon scattering is known to distort optical recordings of cardiac transmembrane potentials; however, this process is not well quantified, hampering interpretation of experimental data. This study presents a novel model, which accurately synthesizes fluorescent recordings over the irregular geometry of the rabbit ventricles. Using the model, the study aims to provide quantification of fluorescent signal distortion for different optical characteristics of the preparation and of the surrounding medium. A bi-domain representation of electrical activity is combined with finite element solutions to the photon diffusion equation simulating both the excitation and emission processes, along with physically realistic boundary conditions at the epicardium, which allow simulation of different experimental setups. We demonstrate that distortion in the optical signal as a result of fluorescent photon scattering is truly a three-dimensional phenomenon and depends critically upon the geometry of the preparation, the scattering properties of the tissue, the direction of wavefront propagation, and the specifics of the experimental setup. Importantly, we show that in an anatomically accurate model of ventricular geometry and fiber orientation, the morphology of the optical signal does not provide reliable information regarding the intramural direction of wavefront propagation. These findings underscore the potential of the new model in interpreting experimental data.

Fluorescence imaging of cardiac propagation: spectral properties and filtering of optical action potentials

Fluorescence imaging using voltage-sensitive dyes is an important tool for studying electrical propagation in the heart. Yet, the low amplitude of the voltage-sensitive component in the fluorescence signal and high acquisition rates dictated by the rapid propagation of the excitation wave front make it difficult to achieve recordings with high signal-to-noise ratios. Although spatially and temporally filtering the acquired signals has become de facto one of the key elements of optical mapping, there is no consensus regarding their use. Here we characterize the spatiotemporal spectra of optically recorded action potentials and determine the distortion produced by conical filters of different sizes. On the basis of these findings, we formulate the criteria for rational selection of filter characteristics. We studied the evolution of the spatial spectra of the propagating wave front after epicardial point stimulation of the isolated, perfused right ventricular free wall of the pig heart stained with di-4-ANEPPS. We found that short-wavelength (<3 mm) spectral components represent primarily noise and surface features of the preparation (coronary vessels, fat, and connective tissue). The time domain of the optical action potential spectrum also lacks high-frequency components (>100 Hz). Both findings are consistent with the reported effect of intrinsic blurring caused by light scattering inside the myocardial wall. The absence of high-frequency spectral components allows the use of aggressive low-pass spatial and temporal filters without affecting the optical action potential morphology. We show examples where the signal-to-noise ratio increased up to 150 with <3% distortion. A generalization of our approach to the rational filter selection in various applications is discussed.