Lifetimes of epicardial rotors in panoramic optical maps of fibrillating swine ventricles

During ventricular fibrillation (VF), electrical activation waves are fragmented, and the heart cannot contract in synchrony. It has been proposed that VF waves emanate from stable periodic sources (often called "mother rotors"). The objective of the present study was to determine if stable rotors are consistently present on the epicardial surface of hearts comparable in size to human hearts. Using new optical mapping technology, we imaged VF from nearly the entire ventricular surface of six isolated swine hearts. Using newly developed pattern analysis algorithms, we identified and tracked VF wave fronts and phase singularities (PS; the pivot point of a reentrant wave front). We introduce the notion of a compound rotor in which the rotor's central PS can change and describe an algorithm for automatically identifying such patterns. This prevents rotor lifetimes from being inappropriately abbreviated by wave front fragmentation and collision events near the PS. We found that stable epicardial rotors were not consistently present during VF: only 1 of 17 VF episodes contained a compound rotor that lasted for the entire mapped interval of 4 s. However, shorter-lived rotors were common; 12.2 (SD 3.3) compound rotors with lifetime >200 ms were visible on the epicardium at any given instant. We conclude that epicardial mother rotors do not drive VF in this experimental model; if mother rotors do exist, they are intramural or septal. This paucity of persistent rotors suggests that individual rotors will eventually terminate by themselves and therefore that the continual formation of new rotors is critical for VF maintenance.     

The pinwheel experiment re-revisited.

Virtual electrode induced phase singularity hypothesis explains the origin of cardiac arrhythmias caused by artificial electrical induction of rotors, i.e. vortex-like self-sustained sources of activity. This mechanism is thought to underlie both stimulus-induced arrhythmias and shock defibrillation therapy. In this paper, we extend this hypothesis to three dimensions using the bidomain model of cardiac tissue. We predict that virtual electrode polarization can produce three topologically distinct types of rotors anchored to: (1) transmural I-shaped scroll wave filaments; (2) near-surface U-shaped scroll wave filaments; and (3) intramural O-shaped scroll wave filaments.     

Epicardial and intramural excitation during ventricular pacing: effect of myocardial structure

Published studies show that ventricular pacing in canine hearts produces three distinct patterns of epicardial excitation: elliptical isochrones near an epicardial pacing site, with asymmetric bulges; areas with high propagation velocity, up to 2 or 3 m/s and numerous breakthrough sites; and lower velocity areas (<1 m/s), where excitation moves across the epicardial projection of the septum. With increasing pacing depth, the magnitude of epicardial potential maxima becomes asymmetric. The electrophysiological mechanisms that generate the distinct patterns have not been fully elucidated. In this study, we investigated those mechanisms experimentally. Under pentobarbital anesthesia, epicardial and intramural excitation isochrone and potential maps have been recorded from 22 exposed or isolated dog hearts, by means of epicardial electrode arrays and transmural plunge electrodes. In five experiments, a ventricular cavity was perfused with diluted Lugol solution. The epicardial bulges result from electrotonic attraction from the helically shaped subepicardial portions of the wave front. The high-velocity patterns and the associated multiple breakthroughs are due to involvement of the Purkinje network. The low velocity at the septum crossing is due to the missing Purkinje involvement in that area. The asymmetric magnitude of the epicardial potential maxima and the shift of the breakthrough sites provoked by deep stimulation are a consequence of the epi-endocardial obliqueness of the intramural fibers. These results improve our understanding of intramural and epicardial propagation during premature ventricular contractions and paced beats. This can be useful for interpreting epicardial maps recorded at surgery or inversely computed from body surface ECGs.     

Model of aquaporin-4 supramolecular assembly in orthogonal arrays based on heterotetrameric association of M1-M23 isoforms

Tetramers of aquaporin-4 (AQP4) water channels form supramolecular assemblies in cell membranes called orthogonal arrays of particles (OAPs). We previously reported evidence that a short (M23) AQP4 isoform produced by alternative splicing forms OAPs by an intermolecular N-terminus interaction, whereas the full-length (M1) AQP4 isoform does not by itself form OAPs but can coassemble with M23 in OAPs as heterotetramers. Here, we developed a model to predict number distributions of OAP size, shape, and composition as a function M23:M1 molar ratio. Model specifications included: random tetrameric assembly of M1 with M23; intertetramer associations between M23 and M23, but not between M1 and M23 or M1; and a free energy constraint limiting OAP size. Model predictions were tested by total internal reflection fluorescence microscopy of AQP4-green-fluorescent protein chimeras and native gel electrophoresis of cells expressing different M23:M1 ratios. Experimentally validated model predictions included: 1), greatly increased OAP size with increasing M23:M1 ratio; 2), marked heterogeneity in OAP size at fixed M23:M1, with increased M23 fraction in larger OAPs; and 3), preferential M1 localization at the periphery of OAPs. The model was also applied to test predictions about binding to AQP4 OAPs of a pathogenic AQP4 autoantibody found in the neuroinflammatory demyelinating disease neuromyelitis optica. Our model of AQP4 OAPs links a molecular-level interaction of AQP4 with its supramolecular assembly in cell membranes.     

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.     

Noninvasive three-dimensional electrocardiographic imaging of ventricular activation sequence

Imaging the myocardial activation sequence is critical for improved diagnosis and treatment of life-threatening cardiac arrhythmias. It is desirable to reveal the underlying cardiac electrical activity throughout the three-dimensional (3-D) myocardium (rather than just the endocardial or epicardial surface) from noninvasive body surface potential measurements. A new 3-D electrocardiographic imaging technique (3-DEIT) based on the boundary element method (BEM) and multiobjective nonlinear optimization has been applied to reconstruct the cardiac activation sequences from body surface potential maps. Ultrafast computerized tomography scanning was performed for subsequent construction of the torso and heart models. Experimental studies were then conducted, during left and right ventricular pacing, in which noninvasive assessment of ventricular activation sequence by means of 3-DEIT was performed simultaneously with 3-D intracardiac mapping (up to 200 intramural sites) using specially designed plunge-needle electrodes in closed-chest rabbits. Estimated activation sequences from 3-DEIT were in good agreement with those constructed from simultaneously recorded intracardiac electrograms in the same animals. Averaged over 100 paced beats (from a total of 10 pacing sites), total activation times were comparable (53.3 +/- 8.1 vs. 49.8 +/- 5.2 ms), the localization error of site of initiation of activation was 5.73 +/- 1.77 mm, and the relative error between the estimated and measured activation sequences was 0.32 +/- 0.06. The present experimental results demonstrate that the 3-D paced ventricular activation sequence can be reconstructed by using noninvasive multisite body surface electrocardiographic measurements and imaging of heart-torso geometry. This new 3-D electrocardiographic imaging modality has the potential to guide catheter-based ablative interventions for the treatment of life-threatening cardiac arrhythmias.     

Validation of three-dimensional conduction models using experimental mapping: are we getting closer?

The anisotropic material properties, irregular geometry, and specialized conduction system of the heart all affect the three-dimensional (3D) spread of electrical activation. A limited number of research groups have tried accounting for these features in 3D conduction models to investigate more thoroughly their observations of cardiac electrical activity in 3D experimental preparations. The full potential of these large scale conduction models, however, has not been realized because of a lack of quantitative validation with experiment. Such validation is critical in order to use the models to predict the electrical response of the myocardium to drugs or electrical stimulation. In this paper, a quantitative, experimental validation of paced activation in a 3D conduction model of a 3 cm × 3 cm × 1 cm section of the ventricular wall is presented. Epicardial and intramural pacing stimuli were applied in the center of a 528 channel electrode plaque sutured to the left ventricle in dogs. Unipolar electrograms were recorded at 2 kHz during and after pacing. Fiber directions within the tissue below the electrodes were estimated histologically and from pace-mapping. Simulated epicardial electrograms were computed for surface paced beats using our 3D bidomain model of the mapped tissue volume incorporating the measured fiber directions. Extracellular potentials and isochronal maps resulting from paced activations in both model and experiment were directly compared. Preliminary results demonstrate that our 3D model reproduces qualitatively such key features of the experimental data as electrogram morphologies and epicardial conduction velocities. Though quantitative agreement between model and experiment was only moderate, the validation approach described herein is an essential first step in assessing the predictive capability of present day conduction models.     

Theoretical analysis of the magnetocardiographic pattern for reentry wave propagation in a three-dimensional human heart model

We present a computational study of reentry wave propagation using electrophysiological models of human cardiac cells and the associated magnetic field map of a human heart. We examined the details of magnetic field variation and related physiological parameters for reentry waves in two-dimensional (2-D) human atrial tissue and a three-dimensional (3-D) human ventricle model. A 3-D mesh system representing the human ventricle was reconstructed from the surface geometry of a human heart. We used existing human cardiac cell models to simulate action potential (AP) propagation in atrial tissue and 3-D ventricular geometry, and a finite element method and the Galerkin approximation to discretize the 3-D domain spatially. The reentry wave was generated using an S1-S2 protocol. The calculations of the magnetic field pattern assumed a horizontally layered conductor for reentry wave propagation in the 3-D ventricle. We also compared the AP and magnetocardiograph (MCG) magnitudes during reentry wave propagation to those during normal wave propagation. The temporal changes in the reentry wave motion and magnetic field map patterns were also analyzed using two well-known MCG parameters: the current dipole direction and strength. The current vector in a reentry wave forms a rotating spiral. We delineated the magnetic field using the changes in the vector angle during a reentry wave, demonstrating that the MCG pattern can be helpful for theoretical analysis of reentry waves.     

Predicting single spikes and spike patterns with the Hindmarsh–Rose model

Most simple neuron models are only able to model traditional spiking behavior. As physiologists discover and classify different electrical phenotypes, computational neuroscientists become interested in using simple phenomenological models that can exhibit these different types of spiking patterns. The Hindmarsh–Rose model is a three-dimensional relaxation oscillator which can show both spiking and bursting patterns and has a chaotic regime. We test the predictive powers of the Hindmarsh–Rose model on two different test databases. We show that the Hindmarsh–Rose model can predict the spiking response of rat layer 5 neocortical pyramidal neurons on a stochastic input signal with a precision comparable to the best known spiking models. We also show that the Hindmarsh–Rose model can capture qualitatively the electrical footprints in a database of different types of neocortical interneurons. When the model parameters are fit from sub-threshold measurements only, the model still captures well the electrical phenotype, which suggests that the sub-threshold signals contain information about the firing patterns of the different neurons.     

Imaging spatio-temporal patterns of long-term potentiation in mouse hippocampus

Spatio-temporal patterns of neuronal activity before and after the induction of long-term potentiation in mouse hippocampal slices were studied using a real-time high-resolution optical recording system. After staining the slices with voltage-sensitive dye, transmitted light images and extracellular field potentials were recorded in response to stimuli applied to CA1 stratum radiatum. Optical and electrical signals in response to single test pulses were enhanced for at least 30 minutes after brief high-frequency stimulation at the same site. In two-pathway experiments, potentiation was restricted to the tetanized pathway. The optical signals demonstrated that both the amplitude and area of the synaptic response were increased, in patterns not predictable from the initial, pretetanus, pattern of activation. Optical signals will be useful for investigating spatio-temporal patterns of synaptic enhancement underlying information storage in the brain.     

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.     

Unmanned Aircraft Systems complement biologging in spatial ecology studies

The knowledge about the spatial ecology and distribution of organisms is important for both basic and applied science. Biologging is one of the most popular methods for obtaining information about spatial distribution of animals, but requires capturing the animals and is often limited by costs and data retrieval. Unmanned Aircraft Systems (UAS) have proven their efficacy for wildlife surveillance and habitat monitoring, but their potential contribution to the prediction of animal distribution patterns and abundance has not been thoroughly evaluated. In this study, we assess the usefulness of UAS overflights to (1) get data to model the distribution of free‐ranging cattle for a comparison with results obtained from biologged (GPS‐GSM collared) cattle and (2) predict species densities for a comparison with actual density in a protected area. UAS and biologging derived data models provided similar distribution patterns. Predictions from the UAS model overestimated cattle densities, which may be associated with higher aggregated distributions of this species. Overall, while the particular researcher interests and species characteristics will influence the method of choice for each study, we demonstrate here that UAS constitute a noninvasive methodology able to provide accurate spatial data useful for ecological research, wildlife management and rangeland planning.     

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.     

2,6-Bis(3,4,5-trihydroxybenzylydene) derivatives of cyclohexanone: novel potent HIV-1 integrase inhibitors that prevent HIV-1 multiplication in cell-based assays

A number of 2,6-bisbenzylidenecyclohexane-1-one derivatives have been synthesized and tested as HIV-1 integrase (IN) inhibitors with the aim of obtaining compounds capable to elicit antiviral activity at non-cytotoxic concentrations in cell-based assays. 3,5-Bis(3,4,5-trihydroxybenzylidene)-4-oxocyclohexaneacetic acid (20d) resulted one of the most potent and selective derivatives in acutely infected MT-4 cells (EC(50) and CC(50) values of 2 and 40 microM, respectively). In enzyme assays with recombinant HIV-1 integrase (rIN), this compound proved able to inhibit both 3'-processing and disintegration with IC(50) values of 0.2 and 0.5 microM, respectively. In order to develop a model capable to predict the anti HIV-IN activity and useful to design novel derivatives, we performed a comparative molecular field analysis (CoMFA) like 3-D-QSAR. In our model the ligands were described quantitatively in the GRID program, and the model was optimized by selecting only the most informative variables in the GOLPE program. We found the predictive ability of the model to increase significantly when the number of variables was reduced from 20,925 to 1327. A Q(2) of 0.73 was obtained with the final model, confirming the predictive ability of the model. By studying the PLS coefficients in informative 3-D contour plots, ideas for the synthesis of new compounds could be generated.     

Virtual electrode polarization in the far field: implications for external defibrillation

We recently suggested that failure of implantable defibrillation therapy may be explained by the virtual electrode-induced phase singularity mechanism. The goal of this study was to identify possible mechanisms of vulnerability and defibrillation by externally applied shocks in vitro. We used bidomain simulations of realistic rabbit heart fibrous geometry to predict the passive polarization throughout the heart induced by external shocks. We also used optical mapping to assess anterior epicardium electrical activity during shocks in Langendorff-perfused rabbit hearts (n = 7). Monophasic shocks of either polarity (10-260 V, 8 ms, 150 microF) were applied during the T wave from a pair of mesh electrodes. Postshock epicardial virtual electrode polarization was observed after all 162 applied shocks, with positive polarization facing the cathode and negative polarization facing the anode, as predicted by the bidomain simulations. During arrhythmogenesis, a new wave front was induced at the boundary between the two regions near the apex but not at the base. It spread across the negatively polarized area toward the base of the heart and reentered on the other side while simultaneously spreading into the depth of the wall. Thus a scroll wave with a ribbon-shaped filament was formed during external shock-induced arrhythmia. Fluorescent imaging and passive bidomain simulations demonstrated that virtual electrode polarization-induced scroll waves underlie mechanisms of shock-induced vulnerability and failure of external defibrillation.     

Evaluating intramural virtual electrodes in the myocardial wedge preparation: simulations of experimental conditions

While defibrillation is the only means for prevention of sudden cardiac death, key aspects of the process, such as the intramural virtual electrodes (VEs), remain controversial. Experimental studies had attempted to assess intramural VEs by using wedge preparations and recording activity from the cut surface; however, applicability of this approach remains unclear. These studies found, surprisingly, that for strong shocks, the entire cut surface was negatively polarized, regardless of boundary conditions. The goal of this study is to examine, by means of bidomain simulations, whether VEs on the cut surface represent a good approximation to VEs in depth of the intact wall. Furthermore, we aim to explore mechanisms that could give rise to negative polarization on the cut surface. A model of wedge preparation was used, in which fiber orientation could be changed, and where the cut surface was subjected to permeable and impermeable boundary conditions. Small-scale mechanisms for polarization were also considered. To determine whether any distortions in the recorded VEs arise from averaging during optical mapping, a model of fluorescent recording was employed. The results indicate that, when an applied field is spatially uniform and impermeable boundary conditions are enforced, regardless of the fiber orientation VEs on the cut surface faithfully represent those intramurally, provided tissue properties are not altered by dissection. Results also demonstrate that VEs are sensitive to the conductive layer thickness above the cut surface. Finally, averaging during fluorescent recordings results in large negative VEs on the cut surface, but these do not arise from small-scale heterogeneities.     

Measurements of transmembrane potential and magnetic field at the apex of the heart

We studied the transmembrane potential and magnetic fields from electrical activity at the apex of the isolated rabbit heart experimentally using optical mapping and superconducting quantum interference device microscopy, and theoretically using monodomain and bidomain models. The cardiac apex has a complex spiral fiber architecture that plays an important role in the development and propagation of action currents during stimulation at the apex. This spiral fiber orientation contains both radial electric currents that contribute to the electrocardiogram and electrically silent circular currents that cannot be detected by the electrocardiogram but are detectable by their magnetic field, B_z. In our experiments, the transmembrane potential, V_m, was first measured optically and then B_z was measured with a superconducting quantum interference device microscope. Based on a simple model of the spiral structure of the apex, V_m was expected to exhibit circular wave front patterns and B_z to reflect the circular component of the action currents. Although the circular V_m wave fronts were detected, the B_z maps were not as simple as expected. However, we observed a pattern consistent with a tilted axis for the apex spiral fiber geometry. We were able to simulate similar patterns in both a monodomain model of a tilted stack of rings of dipole current and a bidomain model of a tilted stack of spiraled cardiac tissue that was stimulated at the apex. The fact that the spatial pattern of the magnetic data was more complex than the simple circles observed for V_m suggests that the magnetic data contain information that cannot be found electrically.     

Sequential pattern of the propagation of excitation in the heart of the sheep Ovis brachyurae

With the aid of intramural multipolar technique, the earliest focus of the depolarization is revealed in the myocard thickness of the cranioventral region of the right atrium of the sheep. From there the depolarization wave with a saw-edged front is spread along the thickness of the atria. The chronotopography of the intramural activation wave front shows a more complicated picture of the atrial excitation than a smiply radial one. Apparently it is connected with the presence of the atrial conduction system. The main mass of myocardium of free ventricle walls and the lower two thirds of the septum are activated by means of multifocal depolarization. The base third of the septum is the last to be activated. These features of excitation of the ovine ventricle myocardium can be explained by special character of distribution of the Purkinje fibers in myocardium ventricles.     

Methods in phenological mapping

Summary Phenological maps of flowering provide useful information about both spatial and temporal patterns of pollen emission, and their use could bring a substantial improvement of aerobiological forecasts. This paper presents a method for preparing flowering maps by computer on the basis of phenological data and topography. Data drawn from topographic maps are processed with the aid of an empirical model, based on the relationships between phenology and environment, for obtaining a phenological delay matrix. From this matrix it is possible to derive automatically various kinds of maps (chronological, synoptical and differential). An application is described, relating to the blooming patterns of a set of wild plants in a mountainous area of Northern Italy.     

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.