Transmembrane voltage changes during unipolar stimulation of rabbit ventricle

A multi-scale simulation of retinal physiology

We present a detailed physiological model of the (human) retina that includes the biochemistry and electrophysiology of phototransduction, neuronal electrical coupling, and the spherical geometry of the eye. The model is a parabolic-elliptic system of partial differential equations based on the mathematical framework of the bi-domain equations, which we have generalized to account for multiple cell-types. We discretize in space with non-uniform finite differences and step through time with a custom adaptive time-stepper that employs a backward differentiation formula and an inexact Newton method. A refinement study confirms the accuracy and efficiency of our numerical method. Numerical simulations using the model compare favorably with experimental findings, such as desensitization to light stimuli and calcium buffering in photoreceptors. Other numerical simulations suggest an interplay between photoreceptor gap junctions and inner segment, but not outer segment, calcium concentration. Applications of this model and simulation include analysis of retinal calcium imaging experiments, the design of electroretinograms, the design of visual prosthetics, and studies of ephaptic coupling within the retina.

Living myocardial slices: Advancing arrhythmia research

Living myocardial slices (LMS) are ultrathin (150-400 µm) sections of intact myocardium that can be used as a comprehensive model for cardiac arrhythmia research. The recent introduction of biomimetic electromechanical cultivation chambers enables long-term cultivation and easy control of living myocardial slices culture conditions. The aim of this review is to present the potential of this biomimetic interface using living myocardial slices in electrophysiological studies outlining advantages, disadvantages and future perspectives of the model. Furthermore, different electrophysiological techniques and their application on living myocardial slices will be discussed. The developments of living myocardial slices in electrophysiology research will hopefully lead to future breakthroughs in the understanding of cardiac arrhythmia mechanisms and the development of novel therapeutic options.

Bidomain modeling of electrical and mechanical properties of cardiac tissue

Throughout the history of cardiac research, there has been a clear need to establish mathematical models to complement experimental studies. In an effort to create a more complete picture of cardiac phenomena, the bidomain model was established in the late 1970s to better understand pacing and defibrillation in the heart. This mathematical model has seen ongoing use in cardiac research, offering mechanistic insight that could not be obtained from experimental pursuits. Introduced from a historical perspective, the origins of the bidomain model are reviewed to provide a foundation for researchers new to the field and those conducting interdisciplinary research. The interplay of theory and experiment with the bidomain model is explored, and the contributions of this model to cardiac biophysics are critically evaluated. Also discussed is the mechanical bidomain model, which is employed to describe mechanotransduction. Current challenges and outstanding questions in the use of the bidomain model are addressed to give a forward-facing perspective of the model in future studies.

Reconstituting electrical conduction in soft tissue: the path to replace the ablationist

Cardiac arrhythmias are a leading cause of morbidity and mortality in the developed world. A common mechanism underlying many of these arrhythmias is re-entry, which may occur when native conduction pathways are disrupted, often by myocardial infarction. Presently, re-entrant arrhythmias are most commonly treated with antiarrhythmic drugs and myocardial ablation, although both treatment methods are associated with adverse side effects and limited efficacy. In recent years, significant advancements in the field of biomaterials science have spurred increased interest in the development of novel therapies that enable restoration of native conduction in damaged or diseased myocardium. In this review, we assess the current landscape of materials-based approaches to eliminating re-entrant arrhythmias. These approaches potentially pave the way for the eventual replacement of myocardial ablation as a preferred therapy for such pathologies.

Anisotropic Cardiac Conduction

Anisotropy is the property of directional dependence. In cardiac tissue, conduction velocity is anisotropic and its orientation is determined by myocyte direction. Cell shape and size, excitability, myocardial fibrosis, gap junction distribution and function are all considered to contribute to anisotropic conduction. In disease states, anisotropic conduction may be enhanced, and is implicated, in the genesis of pathological arrhythmias. The principal mechanism responsible for enhanced anisotropy in disease remains uncertain. Possible contributors include changes in cellular excitability, changes in gap junction distribution or function and cellular uncoupling through interstitial fibrosis. It has recently been demonstrated that myocyte orientation may be identified using diffusion tensor magnetic resonance imaging in explanted hearts, and multisite pacing protocols have been proposed to estimate myocyte orientation and anisotropic conduction in vivo. These tools have the potential to contribute to the understanding of the role of myocyte disarray and anisotropic conduction in arrhythmic states.

Optogenetic Hyperpolarization of Cardiomyocytes Terminates Ventricular Arrhythmia

Cardiac defibrillation to terminate lethal ventricular arrhythmia (VA) is currently performed by applying high energy electrical shocks. In cardiac tissue, electrical shocks induce simultaneously de- and hyperpolarized areas and only depolarized areas are considered to be responsible for VA termination. Because electrical shocks do not allow proper control over spatial extent and level of membrane potential changes, the effects of hyperpolarization have not been explored in the intact heart. In contrast, optogenetic methods allow cell type-selective induction of de- and hyperpolarization with unprecedented temporal and spatial control. To investigate effects of cardiomyocyte hyperpolarization on VA termination, we generated a mouse line with cardiomyocyte-specific expression of the light-driven proton pump ArchT. Isolated cardiomyocytes showed light-induced outward currents and hyperpolarization. Free-running VA were evoked by electrical stimulation of explanted hearts perfused with low K⁺ and the K_ATP channel opener Pinacidil. Optogenetic hyperpolarization was induced by epicardial illumination, which terminated VA with an average efficacy of ∼55%. This value was significantly higher compared to control hearts without illumination or ArchT expression (p = 0.0007). Intracellular recordings with sharp electrodes within the intact heart revealed hyperpolarization and faster action potential upstroke upon illumination, which should fasten conduction. However, conduction speed was lower during illumination suggesting enhanced electrical sink by hyperpolarization underlying VA termination. Thus, selective hyperpolarization in cardiomyocytes is able to terminate VA with a completely new mechanism of increased electrical sink. These novel insights could improve our mechanistic understanding and treatment strategies of VA termination.

Clinical deep brain stimulation strategies for orientation-selective pathway activation

OBJECTIVE

This study investigated stimulation strategies to increase the selectivity of activating axonal pathways within the brain based on their orientations relative to clinical deep brain stimulation (DBS) lead implants.

APPROACH

Previous work has shown how varying electrode shape and controlling the primary electric field direction through preclinical electrode arrays can produce orientation-selective axonal stimulation. Here, we significantly extend those results using computational models to evaluate the degree to which clinical DBS leads can direct stimulus-induced electric fields and generate orientation-selective activation of fiber pathways in the brain. Orientation-selective pulse paradigms were evaluated in conceptual models and in patient-specific models of subthalamic nucleus (STN)-DBS for treating Parkinson's disease.

MAIN RESULTS

Single-contact monopolar or two-contact bipolar stimulation through clinical DBS leads with cylindrical electrodes primarily activated axons orientated parallel to the lead. Conversely, multi-contact monopolar stimulation with a cathode-leading pulse waveform selectively activated axons perpendicular to the DBS lead. Clinical DBS leads with segmented rows of electrodes and a single current source provided additional angular resolution for activating axons oriented 0°, ±22.5°, ±45°, ±67.5°, or 90° relative to the lead shaft. Employing multiple independent current sources to deliver unequal amounts of current through these leads further increased the angular resolution of activation relative to the lead shaft. The patient-specific models indicated that multi-contact cathode configurations, which are rarely used in clinical practice, could increase activation of the hyperdirect pathway collaterals projecting into STN (a putative therapeutic target), while minimizing direct activation of the corticospinal tract of internal capsule, which can elicit sensorimotor side-effects when stimulated.

SIGNIFICANCE

When combined with patient-specific tissue anisotropy and patient-specific anatomical morphologies of neural pathways responsible for therapy and side effects, orientation-selective DBS approaches show potential to significantly improve clinical outcomes of DBS therapy for a range of existing and investigational clinical indications.

Incorporating inductances in tissue-scale models of cardiac electrophysiology

In standard models of cardiac electrophysiology, including the bidomain and monodomain models, local perturbations can propagate at infinite speed. We address this unrealistic property by developing a hyperbolic bidomain model that is based on a generalization of Ohm's law with a Cattaneo-type model for the fluxes. Further, we obtain a hyperbolic monodomain model in the case that the intracellular and extracellular conductivity tensors have the same anisotropy ratio. In one spatial dimension, the hyperbolic monodomain model is equivalent to a cable model that includes axial inductances, and the relaxation times of the Cattaneo fluxes are strictly related to these inductances. A purely linear analysis shows that the inductances are negligible, but models of cardiac electrophysiology are highly nonlinear, and linear predictions may not capture the fully nonlinear dynamics. In fact, contrary to the linear analysis, we show that for simple nonlinear ionic models, an increase in conduction velocity is obtained for small and moderate values of the relaxation time. A similar behavior is also demonstrated with biophysically detailed ionic models. Using the Fenton-Karma model along with a low-order finite element spatial discretization, we numerically analyze differences between the standard monodomain model and the hyperbolic monodomain model. In a simple benchmark test, we show that the propagation of the action potential is strongly influenced by the alignment of the fibers with respect to the mesh in both the parabolic and hyperbolic models when using relatively coarse spatial discretizations. Accurate predictions of the conduction velocity require computational mesh spacings on the order of a single cardiac cell. We also compare the two formulations in the case of spiral break up and atrial fibrillation in an anatomically detailed model of the left atrium, and we examine the effect of intracellular and extracellular inductances on the virtual electrode phenomenon.

Minor Variations in Electrode Pad Placement Impact Defibrillation Success

Defibrillation is essential for resuscitating patients with ventricular fibrillation (VF), but shocks often fail to defibrillate. We hypothesized that small variations in pad placement affect shock success, and that defibrillation waveform and shock dose could compensate for suboptimal pad placement. In 10 swine experiments, electrode pads were attached at 3 adjacent anterolateral positions, less than 3 centimeters apart. At each position, 24 episodes of VF were induced and shocked, 8 episodes for each of 3 defibrillation therapies. This resulted in 9 tested combinations of pad position and defibrillation therapy, with 80 episodes of VF for each combination. An episode consisted of 15 seconds of untreated VF, followed by a first shock and, if necessary, a repeat shock. Episodes were separated by four minutes of recovery. Both electrode pad position and therapy order were randomized by experiment. Primary outcome was defined as successful VF termination after the first shock; secondary outcome was the cumulative success of the first and second shocks. First shock efficacy varied widely across the 9 tested combinations of pad position and defibrillation therapy, ranging from 11.3% to 86.3%. When grouped by therapy, first shock efficacy varied significantly between the 3 pad positions: 38.3%, 48.3%, 36.7% (p = 0.02, ANOVA), and, when grouped by pad position, it varied significantly between therapies: 15.0%, 32.5%, 75.8% (p < 0.001, ANOVA). Cumulative 2-shock success varied significantly with therapy (p < 0.001, ANOVA) but not with pad position (p = 0.30, ANOVA). The lowest first shock success was at one position in 6 of 10 animals, at another position in 4 of 10 animals, and never at the third position. Small variations in pad placement can significantly affect defibrillation shock efficacy. However, anatomical variation between individuals and the challenging conditions of real-world resuscitations make optimal pad placement impractical. Suboptimal pad placement can be overcome with defibrillation waveform and shock dose.

Electrical Pacing of Cardiac Tissue Including Potassium Inward Rectification

In this study cardiac tissue is stimulated electrically through a small unipolar electrode. Numerical simulations predict that around an electrode are adjacent regions of depolarization and hyperpolarization. Experiments have shown that during pacing of resting cardiac tissue the hyperpolarization is often inhibited. Our goal is to determine if the inward rectifying potassium current (I_K1) causes the inhibition of hyperpolarization. Numerical simulations were carried out using the bidomain model with potassium dynamics specified to be inward rectifying. In the simulations, adjacent regions of depolarization and hyperpolarization were observed surrounding the electrode. For cathodal currents the virtual anode produces a hyperpolarization that decreases over time. For long duration pulses the current-voltage curve is non-linear, with very small hyperpolarization compared to depolarization. For short pulses, the hyperpolarization is more prominent. Without the inward potassium rectification, the current voltage curve is linear and the hyperpolarization is evident for both long and short pulses. In conclusion, the inward rectification of the potassium current explains the inhibition of hyperpolarization for long duration stimulus pulses, but not for short duration pulses.

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.

Electrophysiological changes correlated with temperature increases induced by high-intensity focused ultrasound ablation

To gain better understanding of the detailed mechanisms of high-intensity focused ultrasound (HIFU) ablation for cardiac arrhythmias, we investigated how the cellular electrophysiological (EP) changes were correlated with temperature increases and thermal dose (cumulative equivalent minutes [CEM43]) during HIFU application using Langendorff-perfused rabbit hearts. Employing voltage-sensitive dye di-4-ANEPPS, we measured the EP and temperature during HIFU using simultaneous optical mapping and infrared imaging. Both action potential amplitude (APA) and action potential duration at 50% repolarization (APD50) decreased with temperature increases, and APD50 was more thermally sensitive than APA. EP and tissue changes were irreversible when HIFU-induced temperature increased above 52.3 ± 1.4°C and log10(CEM43) above 2.16 ± 0.51 (n = 5), but were reversible when temperature was below 50.1 ± 0.8°C and log10(CEM43) below -0.9 ± 0.3 (n = 9). EP and temperature/thermal dose changes were spatially correlated with HIFU-induced tissue necrosis surrounded by a transition zone.

The strength-interval curve in cardiac tissue

The bidomain model describes the electrical properties of cardiac tissue and is often used to simulate the response of the heart to an electric shock. The strength-interval curve summarizes how refractory tissue is excited. This paper analyzes calculations of the strength-interval curve when a stimulus is applied through a unipolar electrode. In particular, the bidomain model is used to clarify why the cathodal and anodal strength-interval curves are different, and what the mechanism of the “dip” in the anodal strength-interval curve is.

Multiparametric optical mapping of the Langendorff-perfused rabbit heart

Optical imaging and fluorescent probes have significantly advanced research methodology in the field of cardiac electrophysiology in ways that could not have been accomplished by other approaches¹. With the use of the calcium- and voltage-sensitive dyes, optical mapping allows measurement of transmembrane action potentials and calcium transients with high spatial resolution without the physical contact with the tissue. This makes measurements of the cardiac electrical activity possible under many conditions where the use of electrodes is inconvenient or impossible¹. For example, optical recordings provide accurate morphological changes of membrane potential during and immediately after stimulation and defibrillation, while conventional electrode techniques suffer from stimulus-induced artifacts during and after stimuli due to electrode polarization¹.

The Langendorff-perfused rabbit heart is one of the most studied models of human heart physiology and pathophysiology. Many types of arrhythmias observed clinically could be recapitulated in the rabbit heart model. It was shown that wave patterns in the rabbit heart during ventricular arrhythmias, determined by effective size of the heart and the wavelength of reentry, are very similar to that in the human heart². It was also shown that critical aspects of excitation-contraction (EC) coupling in rabbit myocardium, such as the relative contribution of sarcoplasmic reticulum (SR), is very similar to human EC coupling³. Here we present the basic procedures of optical mapping experiments in Langendorff-perfused rabbit hearts, including the Langendorff perfusion system setup, the optical mapping systems setup, the isolation and cannulation of the heart, perfusion and dye-staining of the heart, excitation-contraction uncoupling, and collection of optical signals. These methods could be also applied to the heart from species other than rabbit with adjustments to flow rates, optics, solutions, etc.

Two optical mapping systems are described. The panoramic mapping system is used to map the entire epicardium of the rabbit heart^4-7. This system provides a global view of the evolution of reentrant circuits during arrhythmogenesis and defibrillation, and has been used to study the mechanisms of arrhythmias and antiarrhythmia therapy^8,9. The dual mapping system is used to map the action potential (AP) and calcium transient (CaT) simultaneously from the same field of view^10-13. This approach has enhanced our understanding of the important role of calcium in the electrical alternans and the induction of arrhythmia^14-16.

Effects of unipolar stimulation on voltage and calcium distributions in the isolated rabbit heart

BACKGROUND

The effect of electric stimulation on the polarization of cardiac tissue (virtual electrode effect) is well known; the corresponding response of intracellular calcium concentration ([Ca(2+)](i)) and its dependence on coupling interval between conditioning stimulus (S1) and test stimulus (S2) has yet to be elucidated.

OBJECTIVE

Because uncovering the transmembrane potential (V(m))-[Ca(2+)](i) relationship during an electric shock is imperative for understanding arrhythmia induction and defibrillation, we aimed to study simultaneous V(m) and [Ca(2+)](i) responses to strong unipolar stimulation.

METHODS

We used a dual-camera optical system to image concurrently V (m) and [Ca(2+)](i) responses to unipolar stimulation (20 ms +/- 20 mA) in Langendorff-perfused rabbit hearts. RH-237 and Rhod-2 fluorescent dyes were used to measure V(m) and [Ca(2+)](i), respectively. The S1-S2 interval ranged from 10 to 170 ms to examine stimulation during the action potential.

RESULTS

The [Ca(2+)](i) deflections were less pronounced than changes in V(m) for all S1-S2 intervals. For cathodal stimulation, [Ca(2+)](i) at the central virtual cathode region increased with prolongation of S1-S2 interval. For anodal stimulation, [Ca(2+)](i) at the central virtual anode area decreased with shortening of the S1-S2 interval. At very short S1-S2 intervals (10-20 ms), when S2 polarization was superimposed on the S1 action potential upstroke, the [Ca(2+)](i) distribution did not follow V(m) and produced a more complex pattern. After S2 termination [Ca(2+)](i) exhibited three outcomes in a manner similar to V(m): non-propagating response, break stimulation, and make stimulation.

CONCLUSIONS

Changes in the [Ca(2+)](i) distribution correlate with the behavior of the V (m) distribution for S1-S2 coupling intervals longer than 20 ms; at shorter intervals S2 creates more heterogeneous [Ca(2+)](i) distribution in comparison with V(m). Stimulation in diastole and at very short coupling intervals caused V(m)-[Ca(2+)](i) uncoupling at the regions of positive polarization (virtual cathode).

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.

The role of mechanoelectric feedback in vulnerability to electric shock

Experimental and clinical studies have shown that ventricular dilatation is associated with increased arrhythmogenesis and elevated defibrillation threshold; however, the underlying mechanisms remain poorly understood. The goal of the present study was to test the hypothesis that (1) stretch-activated channel (SAC) recruitment and (2) geometrical deformations in organ shape and fiber architecture lead to increased arrhythmogenesis by electric shocks following acute ventricular dilatation. To elucidate the contribution of these two factors, the study employed, for the first time, a combined electro-mechanical simulation approach. Acute dilatation was simulated in a model of rabbit ventricular mechanics by raising the LV end-diastolic pressure from 0.6 (control) to 4.2 kPa (dilated). The output of the mechanics model was used in the electrophysiological model. Vulnerability to shocks was examined in the control, the dilated ventricles, and in the dilated ventricles that also incorporated currents through SAC as a function of local strain, by constructing vulnerability grids. Results showed that dilatation-induced deformation alone decreased upper limit of vulnerability (ULV) slightly and did not result in increased vulnerability. With SAC recruitment in the dilated ventricles, the number of shock-induced arrhythmia episodes increased by 37% (from 41 to 56) and the lower limit of vulnerability (LLV) decreased from 9 to 7 V/cm, while ULV did not change. The heterogeneous activation of SAC caused by the heterogeneous fiber strain in the ventricular walls was the main reason for increased vulnerability to electric shocks since it caused dispersion of electrophysiological properties in the tissue, resulting in postshock unidirectional block and establishment of reentry.

Reverse Polarity Pacing:The Hemodynamic Benefit of Anodal Currents at Lead Tips forCardiac Resynchronization Therapy

BACKGROUND

Myocardial depolarization can be achieved with currents of either anodal or cathodal polarity. In contrast to conventional cathodal pacing, anodal pacing initially hyperpolarizes tissue and improves myocardial contractility in animal models.

METHODS AND RESULTS

In 13 patients undergoing cardiac resynchronization therapy (CRT) device implantation, we compared the mean left ventricular outflow velocity-time integral (LV-VTI) for anodal and cathodal polarities in three different pacing configurations. Intraoperative continuous-wave Doppler measurements were taken at a fixed interrogation angle, while polarities were switched during unipolar left ventricular, unipolar biventricular, and shared-coil biventricular pacing. Comparisons used identical pacing rates, intervals, and stimulus strengths. Patients had a mean ejection fraction of 0.18 +/- 0.08 and a mean QRS duration of 140 +/- 34 ms. All capture thresholds were less than 4.5 volts at a pulse width of 0.4 ms. Data were suitable for analysis in 37 of the 39 pairs of Doppler measurements. Anodal polarity significantly increased average LV-VTI in 36 of these 37 comparisons. The mean increase in LV-VTI for each configuration with anodal versus cathodal polarity was 2.8 +/- 2.6 cm (P < 0.001). The combined mean LV-VTI for all configurations was similarly higher for anodal polarity (24.4 +/- 11.7 cm) versus cathodal polarity (21.7 +/- 10.9 cm; P < 0.001).

CONCLUSION

Anodal pacing polarity significantly improves a measure of LV function compared to traditional cathodal currents. Anodal pacing, which can be achieved by a simple reversal of pacing circuit polarity, may represent an important therapeutic addition to future resynchronization devices.

Quantification of shock-induced microscopic virtual electrodes assessed by subcellular resolution optical potential mapping in guinea pig papillary muscle

INTRODUCTION

The primary objective of this study was the quantitative description of shock-induced, locally occurring virtual electrodes in natural cardiac tissue.

METHODS AND RESULTS

Multiscale optical potential mapping using 10x, 20x, and 40x magnifying objectives, achieving resolutions of 0.13, 0.065, and 0.033 mm, was performed when applying uniform shocks (+/-10 V/cm, 5 ms) during diastole and action potential plateau. A procedure was developed to identify local potential deviations as depolarizing or hyperpolarizing peaks and to quantify their occurrence and characteristic amplitudes, lateral extents, and dynamics. At shock onset, peaks of either polarity developed significantly faster (tau = 0.92 +/- 0.65 ms, N = 64) than the average bulk polarization (tau = 2.25 +/- 0.96 ms, P < 0.001) and appeared locally fixed, changing their polarity at shock reversal. The mean peak magnitude (21.2 +/- 12 mV) and the amplitude distribution were essentially independent from the magnification. The peak density continuously increased with decreasing peak extent (taken at 70% of the amplitude), reaching a maximum of approximately 3 peaks/mm2 in the range of approximately 30-65 microm. There was no correlation between peak amplitude and size throughout. Potentially exciting peaks were found with a density of 0.04-0.2 peaks/mm2 corresponding to estimated 1-5 peaks/mm3.

CONCLUSIONS

Our results suggest that microscopic inhomogeneities form a substantial substrate for far-field excitation in natural cardiac tissue. Here, we effectively bridged the gap between the extensively studied myocyte cultures and larger heart preparations.

Cathodal stimulation in the recovery phase of a propagating planar wave in the rabbit heart reveals four stimulation mechanisms

The stimulation of cardiac tissue in the recovery phase has significant importance in relation to reentry induction. In the theoretical experiment proposed by Winfree, termed the 'pinwheel' experiment, a point stimulus (S2) is applied in the wake of a freely propagating planar wave (S1). Reentry induced from this S1-S2 pinwheel protocol has been observed experimentally in heart preparations. However, in these experiments, which focused on activation outcomes, only mapping of extracellular voltages has been conducted. The lack of transmembrane potential (Vm) distribution data makes it impossible to analyse the underlying stimulation mechanisms which precede the reentry induction. In this work we sought to elucidate the stimulation mechanisms throughout the heart cycle using the pinwheel protocol. We examined the cardiac tissue responses during and immediately after cathodal stimulation in the refractory wake of a propagating planar wave. The voltage-sensitive dye di-4-ANEPPS was utilized to measure Vm directly from quasi two-dimensional preparations of cryoablated Langendorff-perfused rabbit hearts. Four stimulation mechanisms were observed that depended on the Vm magnitude during S2 cathodal stimulation. Make stimulation always occurred during diastolic stimulation. When stimulation was at the beginning of the relative refractory period (RRP), transitional make-break stimulation was detected. During the RRP the excitation was due to the break mechanism. While approaching the effective refractory period (ERP), the tissue response is characterized by a damped wave mediated response. These four stimulation mechanisms were observed in all hearts whether the S1 planar wave propagation was parallel or perpendicular to the fibre direction. This study is the first examination of Vm and the stimulation mechanisms throughout the cardiac cycle using the pinwheel protocol, and the results have implications in the development and improvement of pacing protocols for artificial cardiostimulators.

Simulations of optical mapping during electroporation

Experiments using optical mapping suggest that electroporation occurs in cardiac tissue when the transmembrane potential, Vm, is observed to be significantly less than +/- 400 mV. Our hypothesis, which we test by numerical simulation, is that Vm is greater than +/- 400 mV at the tissue surface, but optical mapping underestimates Vm because it averages over depth. Results indicate a significant underestimation of Vm. Experimental studies indicate a depolarization of the resting transmembrane potential, Vrest, after a strong shock. In a homogeneous model, electroporation only occurs near the tissue surface. Just as Vm during the stimulus is underestimated due to averaging, we hypothesize that the depolarization of Vrest is also underestimated.

Examination of depth-weighted optical signals during cardiac optical mapping: A simulation study

Optical mapping has become a powerful tool to explore complex cardiac propagation. Many experiments and studies claimed that the fluorescence obtained from tissue surface is the averaged response of the transmembrane potential upon probing depth rather than only on the surface. With the electrical propagation model and the photon transport model, the effects of depth-weighted optical signals are examined both during a normal excitation wave and a spiral wave. Our results indicate that depth-weighted optical signals may infer cardiac activation dynamics, such as the mode and the direction of the propagation, the spatial distribution of depolarization or repolarization.

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.

Average over depth during optical mapping of cardiac propagation

High-resolution optical mapping with voltagesensitive dyes has become a powerful tool to depict complex propagation patterns of cardiac transmembrane potentials. Many studied have proved that this optical signal obtained from tissue surface is the response of the transmembrane potential averaged upon depth rather than only surface. In order to investigate the differentia between the two transmembrane potentials, in this paper, we simulated cardiac propagation using Luo-Rudy (L-R) model and calculated the transmembrance potentials average over depth by an optical decay constant from 0.0 mm to 3.0mm. Our results suggest that the transmembrance potentials weighted average over depth is different from that on the tissue surface, the discrepancy between them depends on the depth of the fluorescence emission of the tissue. If only top layer of tissue (<1.0mm) contributes the fluorescence, optical mapping is an almost accurate representation of surface activation dynamics.

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.

Do Intramural Virtual Electrodes Facilitate Successful Defibrillation? Model-Based Analysis of Experimental Evidence

INTRODUCTION

Recent computer model and experimental studies have suggested that microscopic intramural collagenous planes may facilitate successful defibrillation through the generation of shock-induced virtual electrodes deep within the ventricular wall. Evidence supporting the existence of intramural virtual electrodes has been drawn from several recent studies, which map shock-induced membrane potential (Vm) over the cut transmural surface of dissected segments of porcine left ventricle (LV). The artificially created transmural boundary in these experiments is impermeable to intracellular current. It is not known how this constraint limits the interpretation of these experiments in terms of the shock response of the intact ventricle.

METHODS AND RESULTS

This study uses a realistic 3D computer model of LV myocardium to aid experimental interpretation. The model incorporates a microstructural description of intramural cleavage plane discontinuities measured by confocal microscopy of rat LV. Electrical shocks are applied across the model tissue, with and without introduced transmural boundaries. Shocks of varying strength (4-40 V/cm) are also applied to the model and the response analyzed. Results show that shock-induced Vm changes (deltaVm) on a transmural tissue boundary are significantly different to deltaVm of the intact ventricle, and the extent of difference depends on boundary orientation. However, the presence and qualitative behavior of intramural virtual electrodes is preserved irrespective of boundary placement. The model also confirms experimental observations that most rapid transmural activation occurs for shocks of strength 5-10 V/cm. Two distinct mechanisms suppress virtual electrode propagation, and hence slow tissue activation, outside of this optimal shock strength range.

CONCLUSIONS

This study supports the hypothesis that distributed microscopic intramural virtual electrodes contribute to rapid activation of the ventricular wall during defibrillation.

Use of light absorbers to alter optical interrogation with epi-illumination and transillumination in three-dimensional cardiac models

Cardiac optical mapping currently provides 2-D maps of transmembrane voltage-sensitive fluorescence localized near the tissue surface. Methods for interrogation at different depths are required for studies of arrhythmias and the effects of defibrillation shocks in 3-D cardiac tissue. We model the effects of coloading with a dye that absorbs excitation or fluorescence light on the radius and depth of the interrogated region with specific illumination and collection techniques. Results indicate radii and depths of interrogation are larger for transillumination versus epi-illumination, an effect that is more pronounced for broad-field excitation versus laser scanner. Coloading with a fluorescence absorber lessens interrogated depth for epi-illumination and increases it for transillumination, which is confirmed with measurements using transillumination of heart tissue slices. Coloading with an absorber of excitation light consistently decreases the interrogated depths. Transillumination and coloading also decrease the intensities of collected fluorescence. Thus, localization can be modified with wavelength-specific absorbers at the expense of a reduction in fluorescence intensity.

Cell and tissue responses to electric shocks

AIM

Existing models of myocardial membrane kinetics have not been able to reproduce the experimentally-observed negative bias in the asymmetry of transmembrane potential changes (DeltaV(m)) induced by strong electric shocks. The goals of this study are (1) to demonstrate that this negative bias could be reproduced by the addition, to the membrane model, of electroporation and an outward current, I(a), part of the K(+) flow through the L-type Ca(2+)-channel, and (2) to determine how such modifications in the membrane model affect shock-induced break excitation in a 2D preparation.

METHODS AND RESULTS

We conducted simulations of shocks in bidomain fibres and sheets with membrane dynamics represented by the Luo-Rudy dynamic model (LRd'2000), to which electroporation (LRd + EP model) and the outward current, I(a), activated upon strong shock-induced depolarization (aLRd model) was added. Assuming I(a) is a part of K(+) flow through the L-type Ca(2+)-channel enabled us to reproduce both the experimentally observed rectangularly-shaped positive DeltaV(m) and the value of near 2 of the negative-to-positive DeltaV(m) ratio. In the sheet, I(a) not only contributed to the negative bias in DeltaV(m) asymmetry at sites polarized by physical and virtual electrodes, but also restricted positive DeltaV(m). Electroporation, in its turn, was responsible for the decrease in cathode-break excitation threshold in the aLRd sheet, compared with the other two cases, as well as for the occurrence of the excitation after the shock-end rather than during the shock.

CONCLUSIONS

The incorporation of electroporation and I(a) in a membrane model ensures match between simulation results and experimental data. The use of the aLRd model results in a lower threshold for shock-induced break excitation.

Two-Photon Excitation of di-4-ANEPPS for Optical Recording of Action Potentials in Rabbit Heart

Cardiac action potentials have been measured with single-photon excitation (SPE) of transmembrane voltage-sensitive fluorescent dye. Two-photon excitation (TPE) may have advantages for localization and depth of the tissue region from which the action potential is measured. However measurements of action potentials with SPE have not been demonstrated. We sought to develop a method for TPE of di-4-ANEPPS and test whether the method yields voltage-dependent fluorescence in cardiac tissue. We modified our SPE and ratio-metric fluorescence recording system to use a femtosecond pulsed near-infrared laser. Modifications were made to enhance fluorescence collection efficiency and to block infrared laser light from entering the fluorescence collection system. Fluorescence was collected simultaneously in green (510-570 nm) and red (590-700 nm) wavelength bands. Action potentials were observed in the ratio of the green signal to the red signal, but were not observed above the noise level in either of the individual signals. Incorporation of a common-mode noise subtraction method revealed action potentials in green and red signals. We also found that the di-4-ANEPPS fluorescence emission spectrum for TPE at 930 nm was similar to the emission spectrum for SPE at 488 nm. The multiphoton method may be beneficial for highly localized cardiac optical measurements.

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.

Examination of stimulation mechanism and strength-interval curve in cardiac tissue

Understanding the basic mechanisms of excitability through the cardiac cycle is critical to both the development of new implantable cardiac stimulators and improvement of the pacing protocol. Although numerous works have examined excitability in different phases of the cardiac cycle, no systematic experimental research has been conducted to elucidate the correlation among the virtual electrode polarization pattern, stimulation mechanism, and excitability under unipolar cathodal and anodal stimulation. We used a high-resolution imaging system to study the spatial and temporal stimulation patterns in 20 Langendorff-perfused rabbit hearts. The potential-sensitive dye di-4-ANEPPS was utilized to record the electrical activity using epifluorescence. We delivered S1-S2 unipolar point stimuli with durations of 2-20 ms. The anodal S-I curves displayed a more complex shape in comparison with the cathodal curves. The descent from refractoriness for anodal stimulation was extremely steep, and a local minimum was clearly observed. The subsequent ascending limb had either a dome-shaped maximum or was flattened, appearing as a plateau. The cathodal S-I curves were smoother, closer to a hyperbolic shape. The transition of the stimulation mechanism from break to make always coincided with the final descending phase of both anodal and cathodal S-I curves. The transition is attributed to the bidomain properties of cardiac tissue. The effective refractory period was longer when negative stimuli were delivered than for positive stimulation. Our spatial and temporal analyses of the stimulation patterns near refractoriness show always an excitation mechanism mediated by damped wave propagation after S2 termination.

Use of translucent indium tin oxide to measure stimulatory effects of a passive conductor during field stimulation of rabbit hearts

Biomathematical models and experiments have indicated that passive extracellular conductors influence field stimulation. Because metallic conductors prevent optical mapping under the conductor, we have evaluated a passive translucent indium tin oxide (ITO) thin-film conductor to allow mapping of transmembrane potential (V(m)) and stimulatory current under the conductor. A 1-cm ITO disk was patterned photolithographically and positioned between 0.3-cm(2) mesh shock electrodes on the ventricular epicardium of isolated perfused rabbit hearts stained with 4-{2-[6-(dibutylamino)-2-naphthylenal]ethenyl}-1-(3-sulfopropyl)-, hydroxide, inner salt (di-4-ANEPPS). For a 1-A, 10-ms shock during the action potential plateau, optical maps from fluorescence collected using emission ratiometry (excitation at 488 nm and emissions at 510-570 and >590 nm) indicated that the disk altered V(m) by as much as the height of an action potential. DeltaV(m) became more positive near the edge of the disk, where the ITO conductance gradient was parallel to applied current, and more negative near the opposite edge, where the gradient was not parallel to current. For diastolic shocks, the disk expedited membrane excitation at the sites of positive DeltaV(m) in the heart and in a cardiac model with realistic ITO disk surface and interfacial conductances. Optical maps of ITO transmittance and the model indicated that the disk introduced anodal and cathodal stimulatory current at opposite edges of the disk. Thus ITO allows study of the stimulatory effects of a passive conductor in an electric field.

Mechanisms of superiority of ascending ramp waveforms: new insights into mechanisms of shock-induced vulnerability and defibrillation

Monophasic ascending ramp (AR) and descending ramp (DR) waveforms are known to have significantly different defibrillation thresholds. We hypothesized that this difference arises due to differences in mechanisms of arrhythmia induction for the two waveforms. Rabbit hearts (n = 10) were Langendorff perfused, and AR and DR waveforms (7, 20, and 40 ms) were randomly delivered from two line electrodes placed 10 mm apart on the anterior ventricular epicardium. We optically mapped cellular responses to shocks of various strengths (5, 10, and 20 V/cm) and coupling intervals (CIs; 120, 180, and 300 ms). Optical mapping revealed that maximum virtual electrode polarization (VEP) was reached at significantly different times for AR and DR of the same duration (P < 0.05) for all tested CIs. As a result, VEP for AR were stronger than for DR at the end of the shock. Postshock break excitation resulting from AR generated faster propagation and typically could not form reentry. In contrast, partially dissipated VEP resulting from DR generated slower propagation; the wavefront was able to propagate into deexcited tissue and thus formed a shock-induced reentry circuit. Therefore, for the same delivered energy, AR was less proarrhythmic compared with DR. An active bidomain model was used to confirm the electrophysiological results. The VEP hypothesis explains differences in vulnerability associated with monophasic AR and DR waveforms and, by extension, the superior defibrillation efficacy of the AR waveform compared with the DR waveform.

Asymmetry in membrane responses to electric shocks: insights from bidomain simulations

Models of myocardial membrane dynamics have not been able to reproduce the experimentally observed negative bias in the asymmetry of transmembrane potential changes (DeltaVm) induced by strong electric shocks delivered during the action potential plateau. The goal of this study is to determine what membrane model modifications can bridge this gap between simulation and experiment. We conducted simulations of shocks in bidomain fibers and sheets with membrane dynamics represented by the LRd'2000 model. We found that in the fiber, the negative bias in DeltaVm asymmetry could not be reproduced by addition of electroporation only, but by further addition of hypothetical outward current, Ia, activated upon strong shock-induced depolarization. Furthermore, the experimentally observed rectangularly shaped positive DeltaVm, negative-to-positive DeltaVm ratio (asymmetry ratio) = approximately 2, electroporation occurring at the anode only, and the increase in positive DeltaVm caused by L-type Ca2+-channel blockade were reproduced in the strand only if Ia was assumed to be a part of K+ flow through the L-type Ca2+-channel. In the sheet, Ia not only contributed to the negative bias in DeltaVm asymmetry at sites polarized by physical and virtual electrodes, but also restricted positive DeltaVm. Inclusion of Ia and electroporation is thus the bridge between experiment and simulation.

Characteristics of a charged-coupled-device-based optical mapping system for the study of cardiac arrhythmias

We develop an optical fluorescent mapping system that is able to record the action potential wavefront propagation within cardiac tissue samples with high spatial and temporal resolutions. The system's main component, the fluorescence acquisition device (customized CCD camera), offers a high spatial resolution of 128 x 128 pixels, with 12-bit digitization and a frame rate of 490 frames/s. The system is designed and implemented to image an area of approximately 20 x 20 mm at its minimum object distance of 140 mm, corresponding to a spatial resolution of approximately 3 line pairs/mm. Experiments using this system with di-4-ANEPPS-stained canine cardiac tissues with stimulated action potentials through external electrodes result in successful mappings of the distribution and propagation of the action potential wavefronts, showing the system's sensitivity to the change in fluorescence intensity in regions of action potentials. These data demonstrate this optical mapping system as a powerful device in the study of cardiac arrhythmia mechanisms.

High resolution magnetic images of planar wave fronts reveal bidomain properties of cardiac tissue

We magnetically imaged the magnetic action field and optically imaged the transmembrane potentials generated by planar wavefronts on the surface of the left ventricular wall of Langendorff-perfused isolated rabbit hearts. The magnetic action field images were used to produce a time series of two-dimensional action current maps. Overlaying epifluorescent images allowed us to identify a net current along the wavefront and perpendicular to gradients in the transmembrane potential. This is in contrast to a traditional uniform double-layer model where the net current flows along the gradient in the transmembrane potential. Our findings are supported by numerical simulations that treat cardiac tissue as a bidomain with unequal anisotropies in the intra- and extracellular spaces. Our measurements reveal the anisotropic bidomain nature of cardiac tissue during plane wave propagation. These bidomain effects play an important role in the generation of the whole-heart magnetocardiogram and cannot be ignored.

Examination of optical depth effects on fluorescence imaging of cardiac propagation

Optical mapping with voltage-sensitive dyes provides a high-resolution technique to observe cardiac electrodynamic behavior. Although most studies assume that the fluorescent signal is emitted from the surface layer of cells, the effects of signal attenuation with depth on signal interpretation are still unclear. This simulation study examines the effects of a depth-weighted signal on epicardial activation patterns and filament localization. We simulated filament behavior using a detailed cardiac model, and compared the signal obtained from the top (epicardial) layer of the spatial domain with the calculated weighted signal. General observations included a prolongation of the action upstroke duration, early upstroke initiation, and reduction in signal amplitude in the weighted signal. A shallow filament was found to produce a dual-humped action potential morphology consistent with previously reported observations. Simulated scroll wave breakup exhibited effects such as the false appearance of graded potentials, apparent supramaximal conduction velocities, and a spatially blurred signal with the local amplitude dependent upon the immediate subepicardial activity; the combination of these effects produced a corresponding change in the accuracy of filament localization. Our results indicate that the depth-dependent optical signal has significant consequences on the interpretation of epicardial activation dynamics.

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

Voltage-sensitive fluorescent dyes are commonly used to measure cardiac electrical activity. Recent studies indicate, however, that optical action potentials (OAPs) recorded from the myocardial surface originate from a widely distributed volume beneath the surface and may contain useful information regarding intramural activation. The first step toward obtaining this information is to predict OAPs from known patterns of three-dimensional (3-D) electrical activity. To achieve this goal, we developed a two-stage model in which the output of a 3-D ionic model of electrical excitation serves as the input to an optical model of light scattering and absorption inside heart tissue. The two-stage model permits unique optical signatures to be obtained for given 3-D patterns of electrical activity for direct comparison with experimental data, thus yielding information about intramural electrical activity. To illustrate applications of the model, we simulated surface fluorescence signals produced by 3-D electrical activity during epicardial and endocardial pacing. We discovered that OAP upstroke morphology was highly sensitive to the transmural component of wave front velocity and could be used to predict wave front orientation with respect to the surface. These findings demonstrate the potential of the model for obtaining useful 3-D information about intramural propagation.

Development of an optrode for intramural multisite optical recordings of Vm in the heart

INTRODUCTION

Optical mapping of transmembrane potential (Vm) is an important tool in the investigation of impulse propagation in the heart. It provides valuable information about spatiotemporal changes of Vm that cannot be obtained by other techniques, but it presently is limited to measurements from the heart surfaces. Therefore, the goal of this work was to develop a technique for intramural multisite optical measurements of Vm using fiberoptic technology.

METHODS AND RESULTS

An optrode, a bundle of thin optical fibers, was developed for measuring intramural optical signals at multiple sites in the heart. The optrode consisted of seven fibers with diameter of 225 microm arranged in a hexagonal pattern that were used to deliver excitation light to the myocardium, to collect the emitted fluorescence, and to project the light onto a 16 x 16 array of photodiode detectors. Rabbit hearts were stained with the Vm-sensitive dye RH-237. Fluorescence was excited using a 100-W Hg lamp. Intramural action potentials were recorded at multiple sites separated by 2 mm inside the left ventricle. Signal-to-noise (RMS) ratio was 21.2 +/- 12 (n = 7) without averaging or ratiometry and with negligible cross-talk (<1.9%) between the neighboring photodiodes. The size of the recording area for an individual fiber was estimated at approximately 0.8 mm.

CONCLUSION

These data demonstrate feasibility of multisite transmural measurements of Vm without signal averaging and ratiometry. This technique might become useful in studies of transmural impulse conduction during arrhythmias and defibrillation.

Fluorescence imaging of electrical activity in cardiac cells using an all-solid-state system

Tracking spatial and temporal determinants of cardiac arrhythmogenesis at the cellular level presents challenges to the optical mapping techniques employed. In this paper, we describe a compact system combining two nontraditional low-cost solutions for excitation light sources and emission filters in fluorescence measurements of transmembrane potentials, Vm, or intracellular calcium, [Ca2+]i in cardiac cell networks. This is the first reported use of high-power blue and green light emitting diodes (LEDs), to excite cell monolayers stained with Vm - (di-8-ANEPPS) or [Ca2+]i - (Fluo-3) sensitive dyes. In addition, we use simple techniques for fabrication of suitable thin emission filters with uniform properties, no auto-fluorescence, high durability and good flexibility for imaging Vm or [Ca2+]i. The battery-operated LEDs and the fabricated emission filters, integrated with a fiber-optic system for contact fluorescence imaging, were used as tools to characterize conduction velocity restitution at the macro-scale. The versatility of the LEDs for illumination is further emphasized through 1) demonstration of their usage for epi-illumination recordings at the single-cell level, and 2) demonstration of their unique high-frequency light modulation ability. The LEDs showed excellent stability as excitation light sources for fluorescence measurements; acceptable signal-to-noise ratio and negligible cell photodamage and indicator dye photobleaching were observed.

Field Stimulation of 2-D Sheets of Excitable Tissue

This paper develops equations for the transmembrane potentials (Vm) that occur in two-dimensional (2-D) sheets of tissue in response to field stimulation from an electrode near but not on the surface of the tissue. Comparison of results with those for one dimension shows that an additional term is present in the 2-D equations that influences the evolution of Vm in the interval between the end of the stimulus and the active propagation that may follow. The results provide an analytical framework for understanding Vm in response to field stimulation in two dimensions, both during the tissue's critical linear phase and thereafter.