Contractile and Genetic Characterization of Cardiac Constructs Engineered from Human Induced Pluripotent Stem Cells: Modeling of Tuberous Sclerosis Complex and the Effects of Rapamycin.

The implementation of three-dimensional tissue engineering concurrently with stem cell technology holds great promise for in vitro research in pharmacology and toxicology and modeling cardiac diseases, particularly for rare genetic and pediatric diseases for which animal models, immortal cell lines, and biopsy samples are unavailable. It also allows for a rapid assessment of phenotype-genotype relationships and tissue response to pharmacological manipulation. Mutations in the TSC1 and TSC2 genes lead to dysfunctional mTOR signaling and cause tuberous sclerosis complex (TSC), a genetic disorder that affects multiple organ systems, principally the brain, heart, skin, and kidneys. Here we differentiated healthy (CC3) and tuberous sclerosis (TSP8-15) human induced pluripotent stem cells (hiPSCs) into cardiomyocytes to create engineered cardiac tissue constructs (ECTCs). We investigated and compared their mechano-elastic properties and gene expression and assessed the effects of rapamycin, a potent inhibitor of the mechanistic target of rapamycin (mTOR). The TSP8-15 ECTCs had increased chronotropy compared to healthy ECTCs. Rapamycin induced positive inotropic and chronotropic effects (i.e., increased contractility and beating frequency, respectively) in the CC3 ECTCs but did not cause significant changes in the TSP8-15 ECTCs. A differential gene expression analysis revealed 926 up- and 439 down-regulated genes in the TSP8-15 ECTCs compared to their healthy counterparts. The application of rapamycin initiated the differential expression of 101 and 31 genes in the CC3 and TSP8-15 ECTCs, respectively. A gene ontology analysis showed that in the CC3 ECTCs, the positive inotropic and chronotropic effects of rapamycin correlated with positively regulated biological processes, which were primarily related to the metabolism of lipids and fatty and amino acids, and with negatively regulated processes, which were predominantly associated with cell proliferation and muscle and tissue development. In conclusion, this study describes for the first time an in vitro TSC cardiac tissue model, illustrates the response of normal and TSC ECTCs to rapamycin, and provides new insights into the mechanisms of TSC.

I-Wire Heart-on-a-Chip II: Biomechanical analysis of contractile, three-dimensional cardiomyocyte tissue constructs.

This companion study presents the biomechanical analysis of the "I-Wire" platform using a modified Hill model of muscle mechanics that allows for further characterization of construct function and response to perturbation. The I-Wire engineered cardiac tissue construct (ECTC) is a novel experimental platform to investigate cardiac cell mechanics during auxotonic contraction. Whereas passive biomaterials often exhibit nonlinear and dissipative behavior, active tissue equivalents, such as ECTCs, also expend metabolic energy to perform mechanical work that presents additional challenges in quantifying their properties. The I-Wire model uses the passive mechanical response to increasing applied tension to measure the inherent stress and resistance to stretch of the construct before, during, and after treatments. Both blebbistatin and isoproterenol reduced prestress and construct stiffness; however, blebbistatin treatment abolished subsequent force-generating potential while isoproterenol enhanced this property. We demonstrate that the described model can replicate the response of these constructs to intrinsic changes in force-generating potential in response to both increasing frequency of stimulation and decreasing starting length. This analysis provides a useful mathematical model of the I-Wire platform, increases the number of parameters that can be derived from the device, and serves as a demonstration of quantitative characterization of nonlinear, active biomaterials. We anticipate that this quantitative analysis of I-Wire constructs will prove useful for qualifying patient-specific cardiomyocytes and fibroblasts prior to their utilization for cardiac regenerative medicine. STATEMENT OF SIGNIFICANCE: Passive biomaterials may have non-linear elasticity and losses, but engineered muscle tissue also exhibits time- and force-dependent contractions. Historically, mathematical muscle models include series-elastic, parallel-elastic, contractile, and viscous elements. While hearts-on-a-chip can demonstrate in vitro the contractile properties of engineered cardiac constructs and their response to drugs, most of these use cellular monolayers that cannot be readily probed with controlled forces. The I-Wire platform described in the preceding paper by Sidorov et al. addresses these limitations with three-dimensional tissue constructs to which controlled forces can be applied. In this companion paper, we show how to characterize I-Wire constructs using a non-linear, active Hill model, which should be useful for qualifying cells prior to their use in cardiac regenerative medicine.

I-Wire Heart-on-a-Chip I: Three-dimensional cardiac tissue constructs for physiology and pharmacology.

Engineered 3D cardiac tissue constructs (ECTCs) can replicate complex cardiac physiology under normal and pathological conditions. Currently, most measurements of ECTC contractility are either made isometrically, with fixed length and without control of the applied force, or auxotonically against a variable force, with the length changing during the contraction. The "I-Wire" platform addresses the unmet need to control the force applied to ECTCs while interrogating their passive and active mechanical and electrical characteristics. A six-well plate with inserted PDMS casting molds containing neonatal rat cardiomyocytes cultured with fibrin for 13-15days is mounted on the motorized mechanical stage of an inverted microscope equipped with a fast sCMOS camera. A calibrated flexible probe provides strain load of the ECTC via lateral displacement, and the microscope detects the deflections of both the probe and the ECTC. The ECTCs exhibited longitudinally aligned cardiomyocytes with well-developed sarcomeric structure, recapitulated the Frank-Starling force-tension relationship, and demonstrated expected transmembrane action potentials, electrical and mechanical restitutions, and responses to both ß-adrenergic stimulation and blebbistatin. The I-Wire platform enables creation and mechanical and electrical characterization of ECTCs, and hence can be valuable in the study of cardiac diseases, drug screening, drug development, and the qualification of cells for tissue-engineered regenerative medicine. STATEMENT OF SIGNIFICANCE: There is a growing interest in creating engineered heart tissue constructs for basic cardiac research, applied research in cardiac pharmacology, and repair of damaged hearts. We address an unmet need to characterize fully the performance of these tissues with our simple "I-Wire" assay that allows application of controlled forces to three-dimensional cardiac fiber constructs and measurement of both the electrical and mechanical properties of the construct. The advantage of I-Wire over other approaches is that the constructs being measured are truly three-dimensional, rather than a single layer of cells grown within a microfluidic device. We anticipate that the I-Wire will be extremely useful for the evaluation of myocardial constructs created using cardiomyocytes derived from human induced pluripotent stem cells.

Glutamine and glutamate limit the shortening of action potential duration in anoxia-challenged rabbit hearts.

In clinical conditions, amino acid supplementation is applied to improve contractile function, minimize ischemia/reperfusion injury, and facilitate postoperative recovery. It has been shown that glutamine enhances myocardial ATP/APD (action potential duration) and glutathione/oxidized glutathione ratios, and can increase hexosamine biosynthesis pathway flux, which is believed to play a role in cardioprotection. Here, we studied the effect of glutamine and glutamate on electrical activity in Langendorff-perfused rabbit hearts. The hearts were supplied by Tyrode's media with or without 2.5 mmol/L glutamine and 150 µmol/L glutamate, and exposed to two 6-min anoxias with 20-min recovery in between. Change in APD was detected using a monophasic action potential probe. A nonlinear mixed-effects regression technique was used to evaluate the effect of amino acids on APD over the experiment. Typically, the dynamic of APD change encompasses three phases: short transient increase (more prominent in the first episode), slow decrease, and fast increase (starting with the beginning of recovery). The effect of both anoxic challenge and glutamine/glutamate was cumulative, being more pronounced in the second anoxia. The amino acids' protective effect became largest by the end of anoxia - 20.0% (18.9, 95% CI: [2.6 ms, 35.1 ms]), during the first anoxia and 36.6% (27.1, 95% CI: [7.7 ms, 46.6 ms]), during the second. Following the second anoxia, APD difference between control and supplemented hearts progressively increased, attaining 10.8% (13.6, 95% CI: [4.1 ms, 23.1 ms]) at the experiments' end. Our data reveal APD stabilizing and suggest an antiarrhythmic capacity of amino acid supplementation in anoxic/ischemic conditions.

Transmembrane current imaging in the heart during pacing and fibrillation.

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

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

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

Quantification of transmembrane currents during action potential propagation in the heart.

The measurement, quantitative analysis, theory, and mathematical modeling of transmembrane potential and currents have been an integral part of the field of electrophysiology since its inception. Biophysical modeling of action potential propagation begins with detailed ionic current models for a patch of membrane within a distributed cable model. Voltage-clamp techniques have revolutionized clinical electrophysiology via the characterization of the transmembrane current gating variables; however, this kinetic information alone is insufficient to accurately represent propagation. Other factors, including channel density, membrane area, surface/volume ratio, axial conductivities, etc., are also crucial determinants of transmembrane currents in multicellular tissue but are extremely difficult to measure. Here, we provide, to our knowledge, a novel analytical approach to compute transmembrane currents directly from experimental data, which involves high-temporal (200 kHz) recordings of intra- and extracellular potential with glass microelectrodes from the epicardial surface of isolated rabbit hearts during propagation. We show for the first time, to our knowledge, that during stable planar propagation the biphasic total transmembrane current (I(m)) dipole density during depolarization was ∼0.25 ms in duration and asymmetric in amplitude (peak outward current was ∼95 µA/cm(2) and peak inward current was ∼140 µA/cm(2)), and the peak inward ionic current (I(ion)) during depolarization was ∼260 µA/cm(2) with duration of ∼1.0 ms. Simulations of stable propagation using the ionic current versus transmembrane potential relationship fit from the experimental data reproduced these values better than traditional ionic models. During ventricular fibrillation, peak I(m) was decreased by 50% and peak I(ion) was decreased by 70%. Our results provide, to our knowledge, novel quantitative information that complements voltage- and patch-clamp data.

Identification and characterization of a compound that protects cardiac tissue from human Ether-à-go-go-related gene (hERG)-related drug-induced arrhythmias.

The human Ether-à-go-go-related gene (hERG)-encoded K(+) current, I(Kr) is essential for cardiac repolarization but is also a source of cardiotoxicity because unintended hERG inhibition by diverse pharmaceuticals can cause arrhythmias and sudden cardiac death. We hypothesized that a small molecule that diminishes I(Kr) block by a known hERG antagonist would constitute a first step toward preventing hERG-related arrhythmias and facilitating drug discovery. Using a high-throughput assay, we screened a library of compounds for agents that increase the IC(70) of dofetilide, a well characterized hERG blocker. One compound, VU0405601, with the desired activity was further characterized. In isolated, Langendorff-perfused rabbit hearts, optical mapping revealed that dofetilide-induced arrhythmias were reduced after pretreatment with VU0405601. Patch clamp analysis in stable hERG-HEK cells showed effects on current amplitude, inactivation, and deactivation. VU0405601 increased the IC(50) of dofetilide from 38.7 to 76.3 nM. VU0405601 mitigates the effects of hERG blockers from the extracellular aspect primarily by reducing inactivation, whereas most clinically relevant hERG inhibitors act at an inner pore site. Structure-activity relationships surrounding VU0405601 identified a 3-pyridiyl and a naphthyridine ring system as key structural components important for preventing hERG inhibition by multiple inhibitors. These findings indicate that small molecules can be designed to reduce the sensitivity of hERG to inhibitors.

Continuous-waveform constant-current isolated physiological stimulator.

We have developed an isolated continuous-waveform constant-current physiological stimulator that is powered and controlled by universal serial bus (USB) interface. The stimulator is composed of a custom printed circuit board (PCB), 16-MHz MSP430F2618 microcontroller with two integrated 12-bit digital to analog converters (DAC0, DAC1), high-speed H-Bridge, voltage-controlled current source (VCCS), isolated USB communication and power circuitry, two isolated transistor-transistor logic (TTL) inputs, and a serial 16 × 2 character liquid crystal display. The stimulators are designed to produce current stimuli in the range of ±15 mA indefinitely using a 20V source and to be used in ex vivo cardiac experiments, but they are suitable for use in a wide variety of research or student experiments that require precision control of continuous waveforms or synchronization with external events. The device was designed with customization in mind and has features that allow it to be integrated into current and future experimental setups. Dual TTL inputs allow replacement by two or more traditional stimulators in common experimental configurations. The MSP430 software is written in C++ and compiled with IAR Embedded Workbench 5.20.2. A control program written in C++ runs on a Windows personal computer and has a graphical user interface that allows the user to control all aspects of the device.

Amino acids as metabolic substrates during cardiac ischemia.

The heart is well known as a metabolic omnivore in that it is capable of consuming fatty acids, glucose, ketone bodies, pyruvate, lactate, amino acids and even its own constituent proteins, in order of decreasing preference. The energy from these substrates supports not only mechanical contraction, but also the various transmembrane pumps and transporters required for ionic homeostasis, electrical activity, metabolism and catabolism. Cardiac ischemia - for example, due to compromise of the coronary vasculature or end-stage heart failure - will alter both electrical and metabolic activity. While the effects of myocardial ischemia on electrical propagation and stability have been studied in depth, the effects of ischemia on metabolic substrate preference has not been fully appreciated: oxygen deprivation during ischemia will significantly alter the relative ability of the heart to utilize each of these substrates. Although changes in cardiac metabolism are understood to be an underlying component in almost all cardiac myopathies, the potential contribution of amino acids in maintaining cardiac electrical conductance and stability during ischemia is underappreciated. Despite clear evidence that amino acids exert cardioprotective effects in ischemia and other cardiac disorders, their role in the metabolism of the ischemic heart has yet to be fully elucidated. This review synthesizes the current literature of the metabolic contribution of amino acids during ischemia by analyzing relevant historical and recent research.

Regional increase of extracellular potassium leads to electrical instability and reentry occurrence through the spatial heterogeneity of APD restitution.

The heterogeneities of electrophysiological properties of cardiac tissue are the main factors that control both arrhythmia induction and maintenance. Although the local increase of extracellular potassium ([K(+)](o)) due to coronary occlusion is a well-established metabolic response to acute ischemia, the role of local [K(+)](o) heterogeneity in phase 1a arrhythmias has yet to be determined. In this work, we created local [K(+)](o) heterogeneity and investigated its role in fast pacing response and arrhythmia induction. The left marginal vein of a Langendorff-perfused rabbit heart was cannulated and perfused separately with solutions containing 4, 6, 8, 10, and 12 mM of K(+). The fluorescence dye was utilized to map the voltage distribution. We tested stimulation rates, starting from 400 ms down to 120 ms, with steps of 5-50 ms. We found that local [K(+)](o) heterogeneity causes action potential (AP) alternans, 2:1 conduction block, and wave breaks. The effect of [K(+)](o) heterogeneity on electrical stability and vulnerability to arrhythmia induction was largest during regional perfusion with 10 mM of K(+). We detected three concurrent dynamics: normally propagating activation when excitation waves spread over tissue perfused with normal K(+), alternating 2:2 rhythm near the border of [K(+)](o) heterogeneity, and 2:1 aperiodicity when propagation was within the high [K(+)](o) area. [K(+)](o) elevation changed the AP duration (APD) restitution and shifted the restitution curve toward longer diastolic intervals and shorter APD. We conclude that spatial heterogeneity of the APD restitution, created with regional elevation of [K(+)](o), can lead to AP instability, 2:1 block, and reentry induction.

The potential of dual camera systems for multimodal imaging of cardiac electrophysiology and metabolism.

Fluorescence imaging has become a common modality in cardiac electrodynamics. A single fluorescent parameter is typically measured. Given the growing emphasis on simultaneous imaging of more than one cardiac variable, we present an analysis of the potential of dual camera imaging, using as an example our straightforward dual camera system that allows simultaneous measurement of two dynamic quantities from the same region of the heart. The advantages of our system over others include an optional software camera calibration routine that eliminates the need for precise camera alignment. The system allows for rapid setup, dichroic image separation, dual-rate imaging, and high spatial resolution, and it is generally applicable to any two-camera measurement. This type of imaging system offers the potential for recording simultaneously not only transmembrane potential and intracellular calcium, two frequently measured quantities, but also other signals more directly related to myocardial metabolism, such as [K(+)](e), NADH, and reactive oxygen species, leading to the possibility of correlative multimodal cardiac imaging. We provide a compilation of dye and camera information critical to the design of dual camera systems and experiments.

Myofilament Ca2+ sensitization causes susceptibility to cardiac arrhythmia in mice.

In human cardiomyopathy, anatomical abnormalities such as hypertrophy and fibrosis contribute to the risk of ventricular arrhythmias and sudden death. Here we have shown that increased myofilament Ca2+ sensitivity, also a common feature in both inherited and acquired human cardiomyopathies, created arrhythmia susceptibility in mice, even in the absence of anatomical abnormalities. In mice expressing troponin T mutants that cause hypertrophic cardiomyopathy in humans, the risk of developing ventricular tachycardia was directly proportional to the degree of Ca2+ sensitization caused by the troponin T mutation. Arrhythmia susceptibility was reproduced with the Ca2+-sensitizing agent EMD 57033 and prevented by myofilament Ca2+ desensitization with blebbistatin. Ca2+ sensitization markedly changed the shape of ventricular action potentials, resulting in shorter effective refractory periods, greater beat-to-beat variability of action potential durations, and increased dispersion of ventricular conduction velocities at fast heart rates. Together these effects created an arrhythmogenic substrate. Thus, myofilament Ca2+ sensitization represents a heretofore unrecognized arrhythmia mechanism. The protective effect of blebbistatin provides what we believe to be the first direct evidence that reduction of Ca2+ sensitivity in myofilaments is antiarrhythmic and might be beneficial to individuals with hypertrophic cardiomyopathy.

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).

High-resolution high-speed panoramic cardiac imaging system.

A panoramic cardiac imaging system consisting of three high-speed CCD cameras has been developed to image the surface electrophysiology of a rabbit heart via fluorescence imaging using a voltage-sensitive fluorescent dye. A robust, unique mechanical system was designed to accommodate the three cameras and to adapt to the requirements of future experiments. A unified computer interface was created for this application - a single workstation controls all three CCD cameras, illumination, stimulation, and a stepping motor that rotates the heart. The geometric reconstruction algorithms were adapted from a previous cardiac imaging system. We demonstrate the system by imaging a polymorphic cardiac tachycardia.

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.

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.

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.

Spatiotemporal dynamics of damped propagation in excitable cardiac tissue.

Compared to steadily propagating waves (SPW), damped waves (DW), another solution to the nonlinear wave equation, are seldom studied. In cardiac tissue after electrical stimulation in an SPW wake, we observe DW with diminished amplitude and velocity that either gradually decrease as the DW dies, or exhibit a sharp amplitude increase after a delay to become an SPW. The cardiac DW-SPW transition is a key link in understanding defibrillation and stimulation close to the refractory period, and is ideal for a general study of DW dynamics.