A simulation study of the reaction of human heart to biphasic electrical shocks
© Popp et al; licensee BioMed Central Ltd. 2004
Received: 19 January 2004
Accepted: 22 June 2004
Published: 22 June 2004
This article presents a study, which examines the effects of biphasic electrical shocks on human ventricular tissue. The effects of this type of shock are not yet fully understood. Animal experiments showed the superiority of biphasic shocks over monophasic ones in defibrillation. A mathematical computer simulation can increase the knowledge of human heart behavior.
The research presented in this article was done with different models representing a three-dimensional wedge of ventricular myocardium. The electrophysiology was described with Priebe-Beuckelmann model. The realistic fiber twist, which is specific to human myocardium was included. Planar electrodes were placed at the ends of the longest side of the virtual cardiac wedge, in a bath medium. They were sources of electrical shocks, which varied in magnitude from 0.1 to 5 V. In a second arrangement ring electrodes were placed directly on myocardium for getting a better view on secondary electrical sources. The electrical reaction of the tissue was generated with a bidomain model.
The reaction of the tissue to the electrical shock was specific to the initial imposed characteristics. Depolarization appeared in the first 5 ms in different locations. A further study of the cardiac tissue behavior revealed, which features influence the response of the considered muscle. It was shown that the time needed by the tissue to be totally depolarized is much shorter when a biphasic shock is applied. Each simulation ended only after complete repolarization was achieved. This created the possibility of gathering information from all states corresponding to one cycle of the cardiac rhythm.
The differences between the reaction of the homogeneous tissue and a tissue, which contains cleavage planes, reveals important aspects of superiority of biphasic pulses. ...
Ventricular fibrillation is a major health problem. One of its characteristics is that death can occur in few minutes if no therapeutical intervention is applied. Up to now, electrical shocks are the only effective known therapy against ventricular fibrillation. It is important to understand the related phenomena in order to identify the optimal parameters of the applied defibrillation shocks.
Recently, studies were done to evaluate cardiac rhythms following the first defibrillation shock, comparing biphasic truncated exponential, monophasic damped sinusoidal and monophasic truncated exponential waveforms in patients experiencing out-of-hospital ventricular fibrillation cardiac arrest . The results of this research illustrated the defibrillation chances through a statistical method. The impacts of a specific waveform of electrical shocks can also be studied with computer models, which could actually bring a more detailed look over the phenomena. The theoretically obtained data should be in the end compared to the experimental ones.
The exogenous electric current- applied for defibrillating the heart- is redistributed in the cardiac tissue according to the conductivity variation. This is given by the anisotropic electric properties: the conductivity along the cellular axis being greater than the average conductivity across it. Experimental measurements revealed that the orientation of ventricular myocytes depends on the intramural position . This dependence can be approximated by a helix. In the human left ventricular epicardium the angle of the myocytes fiber orientation is equal to -75°, it is 0° in the midwall and 70° in the endocardial area.
Another factor, which modifies the path of the electric current, is the discontinuous character of the cardiac muscle. Cleavage planes are the most important source of discontinuity. They are parting the myocardium in distinct layers. The muscle sheets have a predominantly radial orientation in apicobasal transmural sections. Nevertheless, these layers are not perfectly parallel. They intersect each other forming an interconnected system .
Studies of electrical response of the heart to external stimuli include inhomogenity and discontinuity of the cardiac tissue. It is acknowledged that defibrillation shocks would not be successful if the myocardium would be homogeneous and continuous .
An interesting feature that can be observed after an electrical shock is applied, is the volume ratio of the tissue, which is immediately affected. The areas from which the depolarization could spread forming depolarization fronts are therefore identified. Muscle discontinuities like blood vessels, collagenous septa and cleavage planes induce the appearance of secondary electrical sources. This means that in the neighboring regions depolarization fronts are produced .
This article reports on the influence of biphasic electric shocks and the corresponding waveforms on human ventricular myocardium. For this investigation a three-dimensional computer model of human cardiac tissue was used. It included realistic fiber orientation and cleavage planes.
Many aspects of physiologic and pathophysiologic behavior of the myocardium could be explained with basic electrophysiologic principles . This motivated the experimenters to gather data of intra-, extra- and intercellular electrophysiologic quantities from specific functional regions and from the whole heart. Flow and concentrations of ions, as well as voltages across membranes were of special interest to them. The components which are considered to be relevant for the cell membrane, were: its conductance, ionic channels, pumps, exchangers and some intracellular structures, e.g. the sarcoplasmic reticulum with the storage and release of Ca 2+. The measurement results were used for constructing mathematical models of different levels of abstraction.
For the simulations presented in this article the Priebe-Beuckelmann model of human ventricular tissue had been used . It was basically constructed from a Luo-Rudy ventricular model . The modifications adopted were five currents based on experimental data obtained from human myocytes: I Kr , I Ks , I Ca , I to and I K1. The remaining currents were scaled to fit human cell data. Through these a more refined description of acquired arrhythmogenic substrates was brought.
As all the other electrophysiological models, the used one provided the possibility of calculating the values of its specific ionic currents of different type passing through the membrane and the transmembrane voltage. The transport of ions was correlated to the resistance of the ion channels, which is nonlinear depending on time and voltage. This type of relationship was given by the existence of both passive and active transport mechanisms of the cell membrane channels.
where C m is the membrane capacity, I mem represents the total transmembrane current and I inter is the intercellular source current. The transmembrane current consists in this model of:
The excitation propagation in the cardiac tissue is described by a bidomain reaction diffusion model . In this representation the myocardium was reproduced as a net of cells embedded in an interstitial medium. The myocytes are electrically coupled by an anisotropic current flow through the intra- and extracellular domain and through the gap junctions. In the two domains, the current flow is calculated using Poisson's equations for stationary electrical fields. This model presents a realistic description of the electrical coupling of human ventricular myocytes. The anisotropic electrical properties are incorporated using realistic fiber orientation of the cardiac muscle. The bidomain equations, which were giving us the variation of I inter and Φ e (extracellular potential) as function of the given V m , were implemented with finite difference method and solved with the Gauss-Seidel technique.
∇ ((σ i + σ e )∇Φ e ) + ∇ (σ i ∇V m ) = 0 (3)
∇ (σ i ∇V m ) + ∇ (σ i ∇Φ e ) = -I inter (4)
where σ i and σ e are conductivity tensors of the intra- and extracellular domain, respectively. The blood is a conductive medium (σ b ) and does not contain excitable cells. According to theses characteristics its mathematical description is given by the following equation:
∇ ((σ b )∇Φ e ) = 0 (5)
The current transition from the cardiac tissue to the monodomain bath region is done with the fulfillment of specific boundary conditions. At the frontier the normal component of current between a bidomain and monodomain tissue is continuous in the extracellular space, the normal component of current between a bidomain and monodomain tissue is zero in the intracellular space, and the extracellular potential and the monodomain potential are continuous .
Two types of electrodes, characterized by different geometric forms, were the source of the electric signal. In both cases they were fixed at the ends of the longest axis (z = 0, z = 2 cm). In case of quadratic planar electrodes (1 × 1 cm2 in size), they were placed in bath medium. Because of the chosen conditions, the applied electric signal had to cross 1 mm of the bath medium before entering the myocardial area. In this way the extracellular potential was homogeneously distributed to the myocytes placed on the boundary surfaces. The second type of electrodes had a ring shape. Their outside diameters were 1 cm and their inside diameters were 0.98 cm. The ring electrodes were directly attached to the myocardium. Through them both monophasic and biphasic electric impulses were delivered. The calculation time step was 10 μs. The simulations were done with a Power Mac G5. The machine had 2 GHz CPU speed. The average time needed for each computation was 10 hours.
The first objective was to study the cardiac tissue reaction to different types of biphasic signals given by planar electrodes. Past studies revealed that cleavage planes are drastically reducing the time needed for the entire tissue to be depolarized after an electrical shock . Firstly a cardiac wedge, which was not fragmentized by cleavage planes was used. In this way the differences between the obtained results were enhanced. A continuous medium needs more time for passing from one phase to another, especially when the depolarization front is planar.
The second objective of the study was to understand the phenomena related to tissue discontinuities. For this, a model, which included interlaminar clefts between layers of cardiomyocytes was used. The virtual wedge had a tissue with realistic fiber orientation that was twice intersected by isolating half filled grids (see fig. 3). While the extracellular domain remains continuous, the intracellular domain contains two planes with a chessboard pattern of voxels with zero intracellular conductivity. The used geometrical pattern is showing that cardiac cleavages are intersecting each other. In this way we obtained a schematic description of the structure of the heart.
The simulation showed that- when the electrodes were rectangular- the depolarization front was continuous (movie 1 [see additional file 1]). For this case the intramural myocardium was depolarized by secondary sources created near to the cleavage planes. The side area closest to the anode was depolarized, while the opposite one was hyperpolarized. When the extracellular potential was strong enough to produce a depolarization front from the depolarized region that was near to the cathode, also the cleavage planes were inducing depolarization fronts.
The presentation of our work started with a comparison between different types of biphasic electric shock, which are used in practice for treating a fibrillating human heart. All along the study the reaction of resting tissue to an electrical shock was shown. This reflects the way healthy tissue reacts, so it is helpful for getting a general fundamental understanding of the phenomena. Information gathered in this work will be used in future studies, which will present the effects of the electrical shocks on fibrillating tissue. The data obtained from our simulations showed that the sinusoidal and the exponentially decaying signals were more efficient in exciting a tissue than the rectangular ones. It was observed that a very strong hyperpolarization induced through the second phase of the electrical shock may lead to the failure of the therapy. This result indicates that the magnitude of the second phase of the electrical shock has to be smaller than the magnitude of its first phase. Nevertheless this magnitude has to be high enough to produce depolarization fronts.
The present article shows the importance of biphasic signals over monophasic ones for ring and planar electrodes. The given comparisons with monophasic impulses is reflecting some of the causes that can lead to the failure of complete depolarization. For example, monophasic shocks of same temporal length and magnitude as biphasic ones, applied through ring electrodes can produce a non excitable area which will stop the depolarization front. Comparing the results it was concluded that it is recommended to create depolarization fronts uniformly distributed in space for producing secondary electrical sources in the intramural region. A planar depolarizating front is activating the areas near to the cleavage planes much easier, so high voltages are not requested. In this way the clefts induce the appearance of secondary electric sources, which will start depolarizing the cardiac inner wall simultaneously with the primary sources. Another reason which sustains the bigger efficiency of the biphasic defibrillating shocks, compared to the monophasic ones, is the fact that a fibrillating tissue has already some refractory regions, so a depolarization produced at both ends (as was the presented reaction of the myocardium to the biphasic electrical impulse) would increase the chance of the entire tissue to be depolarized.
The simulations were done in three-dimensional human ventricular myocardium, which included a realistic fiber orientation and also cleavage planes. All these features, put together form a very complex representation of the cardiac muscle. The clefts constitute one of the most important features of the model, since they have been poorly studied from the electrical point of view in the past . Considering this it proves to be of high importance to understand their basic role. If a fibrillation would be ongoing in the virtual myocardium an electrically chaotic behavior would be present. In such circumstances it would be very difficult to identify the exact influence of the cleavage planes. This information is needed in order to take a good decision if the clefts should or not be included in the defibrillation simulation of the fibrillating cardiac human heart.
The chances of the defibrillation success depend on many factors. This is why the article does not exclusively address the total excitation time of the modeled tissue. On contrary, it is also focused on identifying the depolarization sources. These two combined information give a measure of how fast the myocardium can be brought in an electrically equivalent state. The affirmations related to the activation of the-cleavage planes, were the results of the statistics done with many computer experiments. At present we are developing a mathematical theory in order to explain the statistics.
- Carpenter J, Rea T, Murray JA, Kudenchuk PJ, Eisenberg MS: Deflbrillation waveform and post-shock rhythm in out-of-hospital ventricular fibrillation cardiac arrest. Resuscitation. 2003, 59: 189-96. 10.1016/S0300-9572(03)00183-7.View ArticlePubMedGoogle Scholar
- Streeter DD: Gross Morphology and fiber geometry of the Heart. In Handbook of Physiology: The Cardiovascular System. Edited by: Bethesda B. 1979, American Physiology Society, I: 61-112.Google Scholar
- Scollan DF, Holmes A, Zhang J, Winslow RL: Reconstruction of cardiac ventricular geometry and fiber orienta tion using magnetic resonance imaging. Ann Biomed Eng. 2000, 28 (8): 934-944. 10.1114/1.1312188.View ArticlePubMedPubMed CentralGoogle Scholar
- Trayanova N: Far-field stimulation of cardiac tissue. Herzschrittmacher-Therapie und Elektrophysiologie. 1999, 10: 137-148. 10.1007/s003990050055.View ArticleGoogle Scholar
- Hooks D, Tomlinson K, Marden S, LeGrice I, Smaill B, Pullan A, Hunter P: Cardiac microstructure: implications for electrical propagation and defibrillation in the heart. Circulation Research. 2002, 91 (4): 331-338. 10.1161/01.RES.0000031957.70034.89.View ArticlePubMedGoogle Scholar
- Bett GCL, Rasmusson RL: Computer Models of Ion Channels. In Quantitative Cardiac Electrophysiology. Edited by: Cabo C, Rosenbaum DS. 2002, New York: Marcel Dekker, Inc, 1-60.View ArticleGoogle Scholar
- Priebe L, Beuckelmann DI: Simulation Study of Cellular Electric Properties in Heart Failure. Circ Res. 1998, 82: 1206-1223.View ArticlePubMedGoogle Scholar
- Luo CH, Rudy Y: A Model of the Ventricular Cardiac Action Potential. Circ Res. 1991, 68 (6): 1501-1526.View ArticlePubMedGoogle Scholar
- Henriquez CS, Muzikant AL, Smoak CK: Anisotropy, Fiber Curvature and Bath Loading Effects on Activation in Thin and Thick Cardiac Tissue Preparations: Simulations in a Three-Dimensional Bidomain Model. J Cardiovasc Electrophysiol. 1996, 7: 424-444.View ArticlePubMedGoogle Scholar
- Krassowska W, Neu JC: Effective boundary conditions for syncytial tissues. IEEE Trans Biomed Eng. 1994, 41: 143-150. 10.1109/10.284925.View ArticlePubMedGoogle Scholar
- Popp IM, Seemann G, Dössel O: Investigation of the Influence of Electric Fields on Human Ventricular My ocardium Including Realistic Fibre Orientation: A Simulation Study. In Computers in Cardiology. 2003, 30 (IEEE 2003): 213-216.Google Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2261/4/9/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article: verbatim copying and redistribution of this article are permitted in all media for any purpose, provided this notice is preserved along with the article's original URL.