Effect of strength and timing of transmembrane current pulses on isolated ventricular myocytes.

INTRODUCTION: Little is known about how the amplitude and timing of transmembrane current pulses affect transmembrane potential (Vm) and action potential duration (APD) in isolated myocytes. METHODS AND RESULTS: Ten ventricular myocytes were isolated from five rabbit hearts. Each cell was paced at an S1 cycle length of 250 msec, and S2 pulses of 10-msec duration were delivered at various strengths and time intervals. For all S2 strengths (0.2 to 1.5 nA), the magnitude of changes in Vm did not depend on polarity during the plateau, but were larger for depolarizing pulses during phase 3 repolarization. However, the magnitude of changes in APD varied with polarity during the entire action potential for strengths ranging from 0.5 to 1.5 nA. Greater changes in APD occurred for hyperpolarizing pulses during the plateau and depolarizing pulses during phase 3. In addition, we used a cardiac phase variable to quantify the current threshold for regenerative depolarization and repolarization as a function of prestimulus Vm. Regenerative depolarization occurred during phase 3 repolarization, and its current threshold was less than that required for regenerative repolarization that occurred during the plateau. These data were compared to computer simulations in a patch of membrane represented by Luo-Rudy dynamic kinetics, and the results were qualitatively similar, including the higher threshold for regenerative repolarization compared to regenerative depolarization. CONCLUSION: This characterization of the nonlinear response of isolated cells to transmembrane current, including phase resetting, should aid in understanding the mechanisms of defibrillation because shock-induced changes in Vm and APD have been implicated as important factors in determining defibrillation success.

Computer simulation of stair falls to investigate scenarios in child abuse.

OBJECTIVES: To demonstrate the usefulness of computer simulation techniques in the investigation of pediatric stair falls. Since stair falls are a common falsely reported injury scenario in child abuse, our specific aim was to investigate the influence of stair characteristics on injury biomechanics of pediatric stair falls by using a computer simulation model. Our long-term goal is to use knowledge of biomechanics to aid in distinguishing between accidents and abuse. METHODS: A computer simulation model of a 3-year-old child falling down stairs was developed using commercially available simulation software. This model was used to investigate the influence that stair characteristics have on biomechanical measures associated with injury risk. Since femur fractures occur in unintentional and abuse scenarios, biomechanical measures were focused on the lower extremities. RESULTS: The number and slope of steps and stair surface friction and elasticity were found to affect biomechanical measures associated with injury risk. CONCLUSIONS: Computer simulation techniques are useful for investigating the biomechanics of stair falls. Using our simulation model, we determined that stair characteristics have an effect on potential for lower extremity injuries. Although absolute values of biomechanical measures should not be relied on in an unvalidated model such as this, relationships between accident-environment factors and biomechanical measures can be studied through simulation. Future efforts will focus on model validation.

Virtual electrode polarization leads to reentry in the far field.

INTRODUCTION: Our previous article examined cardiac vulnerability to reentry in the near field within the framework of the virtual electrode polarization (VEP) concept. The present study extends this examination to the far field and compares its predictions to the critical point hypothesis. METHODS AND RESULTS: We simulate the electrical behavior of a sheet of myocardium using a two-dimensional bidomain model. The fiber field is extrapolated from a set of rabbit heart fiber directions obtained experimentally. An S1 stimulus is applied along the top or left border. An extracellular line electrode on the top delivers a cathodal or anodal S2 stimulus. A VEP pattern matching that seen experimentally is observed and covers the entire sheet, thus constituting a far-field effect. Reentry arises from break excitation, make excitation, or a combination of both, and subsequent propagation through deexcited and recovered areas. Reentry occurs in cross-field, parallel-field, and uniform refractoriness protocols. For long coupling intervals (CIs) above CImake(min) (defined as the shortest CI at which make excitation can take place), rotors move away from the cathodal electrode and the S1 site for increases in S2 strength and CI, respectively. For cathodal S2 stimuli, findings are consistent with the critical point hypothesis. For CIs below CImake(min), reentry is initiated by break excitation only, and the resulting reentrant patterns are no longer consistent with those predicted by the critical point hypothesis. CONCLUSION: Shock-induced VEP can explain vulnerability in the far field. The VEP theory of vulnerability encompasses the critical point hypothesis for cathodal S2 shocks at long CIs.

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.

Roles of electric field and fiber structure in cardiac electric stimulation.

This study investigated roles of the variation of extracellular voltage gradient (VG) over space and cardiac fibers in production of transmembrane voltage changes (DeltaV(m)) during shocks. Eleven isolated rabbit hearts were arterially perfused with solution containing V(m)-sensitive fluorescent dye (di-4-ANEPPS). The epicardium received shocks from symmetrical or asymmetrical electrodes to produce nominally uniform or nonuniform VGs. Extracellular electric field and DeltaV(m) produced by shocks in the absolute refractory period were measured with electrodes and a laser scanner and were simulated with a bidomain computer model that incorporated the anterior left ventricular epicardial fiber field. Measurements and simulations showed that fibers distorted extracellular voltages and influenced the DeltaV(m). For both uniform and nonuniform shocks, DeltaV(m) depended primarily on second spatial derivatives of extracellular voltages, whereas the VGs played a smaller role. Thus, 1) fiber structure influences the extracellular electric field and the distribution of DeltaV(m); 2) the DeltaV(m) depend on second spatial derivatives of extracellular voltage.

Effects of electroporation on the transmembrane potential distribution in a two-dimensional bidomain model of cardiac tissue.

INTRODUCTION: Defibrillation shocks, when delivered through internal electrodes, establish transmembrane potentials (Vm) large enough to electroporate the membrane of cardiac cells. The effects of such shocks on the transmembrane potential distribution are investigated in a two-dimensional rectangular sheet of cardiac muscle modeled as a bidomain with unequal anisotropy ratios. METHODS AND RESULTS: The membrane is represented by a capacitance Cm, a leakage conductance g(l) and a variable electroporation conductance G, whose rate of growth depends exponentially on the square of Vm. The stimulating current Io, 0.05-20 A/m, is delivered through a pair of electrodes placed 2 cm apart for stimulation along fibers and 1 cm apart for stimulation across fibers. Computer simulations reveal three categories of response to Io: (1) Weak Io, below 0.2 A/m, cause essentially no electroporation, and Vm increases proportionally to Io. (2) Strong Io, between 0.2 and 2.5 A/m, electroporate tissue under the physical electrode. Vm is no longer proportional to Io; in the electroporated region, the growth of Vm is halted and in the region of reversed polarity (virtual electrode), the growth of Vm is accelerated. (3) Very strong Io, above 2.5 A/m, electroporate tissue under the physical and the virtual electrodes. The growth of Vm in all electroporated regions is halted, and a further increase of Io increases both the extent of the electroporated regions and the electroporation conductance G. CONCLUSION: These results indicate that electroporation of the cardiac membrane plays an important role in the distribution of Vm induced by defibrillation strength shocks.

Extension of refractoriness in a model of cardiac defibrillation.

This simulation study presents an inquiry into the mechanisms by which a strong electric shock halts life-threatening cardiac arrhythmias. It examines the "extension of refractoriness" hypothesis for defibrillation which postulates that the shock induces an extension of the refractory period of cardiac cells thus blocking propagating waves of arrhythmia and fibrillation. The present study uses a model of the defibrillation process that represents a sheet of myocardium as a biodomain with unequal anisotropy ratios. The tissue consists of curved fibers in which spiral wave reentry is initiated. The defibrillation shock is delivered via two line electrodes that occupy opposite tissue boundaries. Simulation results demonstrate that a large-scale region of depolarization is induced throughout most of the tissue. This depolarization extends the refractoriness of the cells in the region. In addition, new wavefronts are generated from the regions of induced hyperpolarization that further restrict the spiral wave pathway and cause its termination.

Impact of transvenous lead position on active-can ICD defibrillation: a computer simulation study.

Optimizing lead placement in transvenous defibrillation remains central to the clinical aspects of the defibrillation procedure. Studies involving superior vena cava (SVC) return electrodes have found that left ventricular (LV) leads or septal positioning of the right ventricular (RV) lead minimizes the voltage defibrillation threshold (VDFT) in endocardial lead-->SVC defibrillation systems. However, similar studies have not been conducted for active-can configurations. The goal of this study was to determine the optimal lead position to minimize the VDFT for systems incorporating an active can. This study used a high resolution finite element model of a human torso that includes the fiber architecture of the ventricular myocardium to find the role of lead positioning in a transvenous LEAD-->can defibrillation electrode system. It was found that, among single lead systems, posterior positioning of leads in the right ventricle lowers VDFTs appreciably. Furthermore, a septal location of leads resulted in lower VDFTs than free-wall positioning. Increasing the number of leads, and thus the effective lead surface area in the right ventricle also resulted in lower VDFTs. However, the lead configuration that resulted in the lowest VDFTs is a combination of mid-cavity right ventricle lead and a mid-cavity left ventricle lead. The addition of a left ventricular lead resulted in a reduction in the size of the low gradient regions and a change of its location from the left ventricular free wall to the septal wall.

Influence of anisotropy on local and global measures of potential gradient in computer models of defibrillation.

A heart-torso model including fiber orientation is used to calculate electric field strength in an active-can transvenous defibrillation system and estimate errors due to inadequate description of the anisotropy of the myocardium. Using a minimum potential gradient (5 V/cm) in a critical mass (95%) of the tissue, the estimated defibrillation voltage threshold for a right ventricular transvenous lead placement differs by only 4.5% when using isotropic myocardial conductivity compared to a model with realistic fiber architecture. In addition, pointwise comparisons of the two solutions reveal differences of 10.8% rms in potential gradient strength and 31.6% rms in current density magnitude in the myocardium, resulting in a change in the location of the low gradient regions. These results suggest that if a minimum potential gradient throughout the heart is necessary to avoid reinitiation of fibrillatory wave fronts, then isotropic models are adequate for modeling the electric field in the heart. Alternatively, the model demonstrates the use of physiologically based descriptions of anisotropy and fiber orientation, which will soon allow simulations of shock induced membrane polarization during defibrillation.