|Computational quantification of the cardiac energy consumption during intra-aortic balloon pumping using a cardiac electromechanics model.|
|Jump to Full Text|
|PMID: 23341718 Owner: NLM Status: MEDLINE|
|To quantify the reduction in workload during intra-aortic balloon pump (IABP) therapy, indirect parameters are used, such as the mean arterial pressure during diastole, product of heart rate and peak systolic pressure, and pressure-volume area. Therefore, we investigated the cardiac energy consumption during IABP therapy using a cardiac electromechanics model. We incorporated an IABP function into a previously developed electromechanical model of the ventricle with a lumped model of the circulatory system and investigated the cardiac energy consumption at different IABP inflation volumes. When the IABP was used at inflation level 5, the cardiac output and stroke volume increased 11%, the ejection fraction increased 21%, the stroke work decreased 1%, the mean arterial pressure increased 10%, and the ATP consumption decreased 12%. These results show that although the ATP consumption is decreased significantly, stroke work is decreased only slightly, which indicates that the IABP helps the failed ventricle to pump blood efficiently.|
|Ki Moo Lim; Jeong Sang Lee; Min-Soo Gyeong; Jae-Sung Choi; Seong Wook Choi; Eun Bo Shim|
Related Documents :
|22967988 - Xanthine oxidase contributes to mitochondrial ros generation in an experimental model o...
23838208 - Systolic arterial blood pressure in small-breed dogs with degenerative mitral valve dis...
2499228 - Hypoxic hazards of traditional paper bag rebreathing in hyperventilating patients.
9627218 - Myocardial hemodynamic and metabolic changes during abdominal insufflation with carbon ...
2648468 - Duplex pulsed doppler us versus intracranial pressure in the neonate: clinical and expe...
11324978 - Ocular complications of latanoprost in uveitic glaucoma: three case reports.
|Type: Journal Article; Research Support, Non-U.S. Gov't Date: 2013-01-08|
|Title: Journal of Korean medical science Volume: 28 ISSN: 1598-6357 ISO Abbreviation: J. Korean Med. Sci. Publication Date: 2013 Jan|
|Created Date: 2013-01-23 Completed Date: 2013-07-01 Revised Date: 2013-07-11|
Medline Journal Info:
|Nlm Unique ID: 8703518 Medline TA: J Korean Med Sci Country: Korea (South)|
|Languages: eng Pagination: 93-9 Citation Subset: IM|
|Department of Medical IT Convergence Engineering, Kumoh National Institute of Technology, Gumi, Korea.|
|APA/MLA Format Download EndNote Download BibTex|
Heart Failure / pathology
Intra-Aortic Balloon Pumping*
Ventricular Function, Left
Journal ID (nlm-ta): J Korean Med Sci
Journal ID (iso-abbrev): J. Korean Med. Sci
Journal ID (publisher-id): JKMS
Publisher: The Korean Academy of Medical Sciences
© 2013 The Korean Academy of Medical Sciences.
Received Day: 13 Month: 8 Year: 2012
Accepted Day: 09 Month: 9 Year: 2012
Print publication date: Month: 1 Year: 2013
Electronic publication date: Day: 08 Month: 1 Year: 2013
Volume: 28 Issue: 1
First Page: 93 Last Page: 99
PubMed Id: 23341718
|Computational Quantification of the Cardiac Energy Consumption during Intra-Aortic Balloon Pumping Using a Cardiac Electromechanics Model|
|Ki Moo Lim1*|
|Jeong Sang Lee2*|
|Seong Wook Choi3|
|Eun Bo Shim3|
1Department of Medical IT Convergence Engineering, Kumoh National Institute of Technology, Gumi, Korea.
2Department of Thoracic and Cardiovascular Surgery, Seoul National University College of Medicine, & SMG-SNU Boramae Hospital, Seoul, Korea.
3Department of Mechanical and Biomedical Engineering, Kangwon National University, Chuncheon, Korea.
Address for Correspondence: Eun Bo Shim, PhD. Department of Mechanical & Biomedical Engineering, Kangwon National University, 1 Kangwondaehak-gil, Chuncheon, 200-701, Korea. Tel: +82.033-250-6318, Fax: +82.033-257-6595, firstname.lastname@example.org
[equal] *Ki Moo Lim and Jeong Sang Lee contributed equally to this work.
The intra-aortic balloon pump (IABP) is used to increase myocardial oxygen perfusion, while at the same time increasing cardiac output and decreasing the workload of the ventricle. This is realized via counter-pulsation of the IABP. Several parameters affect IABP treatment, such as the volume of the balloon, its location in the aorta, its rate of inflation and deflation, and the phase of IABP pumping. Several researchers have examined the effects of various parameters on IABP function in terms of increasing myocardial oxygen perfusion and decreasing the workload of the ventricle. When the volume of blood used to inflate and deflate the balloon increases, the IABP functions were improved (1). Inflation of the IABP at end-systole and deflation of the IABP at end-diastole were the optimal pumping phase in order to decrease the ventricle workload and increase aortic pressure during diastole, thereby increasing myocardial oxygen perfusion (2).
To quantify the reduction in workload, indirect parameters are used, such as the mean systolic pressure (3), product of heart rate and peak systolic pressure (3), and pressure-volume area (4). No study has quantified the reduction in workload with IABP directly because experimental methods for documenting and evaluating myocardial energy consumption throughout the ventricular volume are hampered by low spatiotemporal resolution. Computational modeling is an alternative approach that overcomes this limitation. However, the previous computational models of the cardiovascular system and IABP used a simple lumped model of the cardiovascular system, which cannot calculate the contractile cardiac energy consumption directly.
Previously, we developed a computational model of left ventricular assist device support using a three-dimensional electromechanical model of failing canine ventricles together with a lumped model of the circulatory system (5, 6). This model enabled us to quantify the changes in the local contractile energy consumption of the myocardium. In this study, we incorporated IABP function into the electromechanical model of a ventricle with a lumped model of the circulatory system. Then, we used the model to investigate the changes in the contractile energy consumption of the myocardium with an IABP using different inflation volumes.
The IABP-implanted cardiovascular system model was implemented by combining a three-dimensional (3D) finite element (FE) electromechanical model of the failed heart, a lumped model of the circulatory system, and a model of IABP function. Fig. 1 shows a schematic of the integrated model (A) and myofilament dynamics, which were adapted from Rice et al. (7). The models are described in detail below.
We used an anatomically accurate 3D-FE electromechanical model of failing canine ventricles, which incorporated a biophysical representation of cardiac myofilament dynamics (8). The electrical model based on the membrane dynamic model of Fox et al. (9) used magnetic resonance (MR) image-based canine ventricle geometry with approximately 1,400,000 nodes and 1,700,000 elements, which were developed using linear interpolation. The mechanical model based on the myofilament dynamic model of Rice et al. (7) used the same canine ventricle geometry with 356 nodes and 172 elements, which were developed using Hermite polynomials. Both meshes also have realistic fiber and laminar sheet structural information obtained from the diffusion tensor (DT) MR image. The passive mechanical properties of the myocardium were assumed to be transversely isotropic, hyperelastic, and nearly incompressible, defined by an exponential strain energy function. The local electrical activation times in the electrical mesh were mapped to the mechanical mesh. Fig. 2 shows the electrical activation time mapped to the mechanical mesh. Here, the local electrical activation time is defined as the instant at which transmembrane voltage exceeds 0 mV.
As an electrical wave propagates through the heart, the depolarization of each myocyte initiates the release of calcium (Ca) from intracellular stores. This is followed by the binding of Ca to troponin C and cross-bridge cycling. The cross-bridge cycling forms the basis for contractile protein movement and the development of active tension in the cell, resulting in deformation of the ventricles. These physiological contractile mechanisms are considered in the myofilament model (7).
To consider remodeling of the passive electromechanical properties associated with heart failure (HF), several changes were made to the electromechanical model (Table 1). First, we used the dilated geometry of a canine ventricle with heart failure (HF) (8). Then, the electrical conductivity was reduced by 30%, allowing for a total electrical activation time of 150 ms, according to the experimental results of Helm et al. (10). Third, to account for the increased stiffness of the failing myocardium, the passive scaling constant in the strain-energy function was increased fivefold (11). Finally, to introduce systolic dysfunction, the peak and time constant of the calcium transient function ([Ca](t) of Eq. 55 in the myofilament model of Rice et al. (7)), which served as an input, were reduced to 70% of the normal values. This resulted in a calcium transient with a reduced peak and prolonged relaxation rate, which are the key features of the remodeled Ca transient in HF (12).
We combined the FE heart model with a lumped model of the circulatory system and IABP to predict the cardiovascular response during IABP therapy. The model is formulated as an electrical analog model consisting of elements such as resistors, capacitors, and diodes. Using the hemodynamic parameters related to blood circulation, this model calculates the pressure, volume, and blood flow for each of the six compartments making up the vascular system and atria (Fig. 1): the systemic arteries and veins, pulmonary arteries and veins, left atrium, and right atrium. Individual elements in every compartment are modeled by considering the local resistance to blood flow and compliance of blood vessels. The atria are modeled as two-chamber pumps with time-varying compliance. The inertia of blood is ignored due to its relatively small effect on the canine cardiovascular response. The four heart valves that surround the ventricles and control the blood flow direction are modeled by diodes and resistances.
Briefly, the IABP component was modeled as the time-varying compliance of the systemic arteries. To generalize the patterns of inflation and deflation of an IABP, a harmonic waveform was used for the time-varying compliance of the systemic arteries. The harmonic waveforms for the compliance of the aorta were expressed as follows:
Here, CSA,IABP is the time-varying compliance of the systemic arteries with the IABP, CSA is the compliance of the systemic arteries without the IABP, SFIABP is a scale factor for the IABP, sf is the level of the scale factor, BCL indicates the basic cycle length of the ventricle, and ø indicates the time shift between the ventricular contraction cycle and IABP inflation cycle. The sf is major parameter which is proportional to the stroke volume of IABP and related to the pumping compliance of IABP. In this study, sf was given values of 0.05, 0.1, 0.15, 0.2, and 0.25, and ø was set to 3.66 radians (350 ms shift), instead of end-systole.
Fig. 2 shows the simulated pressure waveforms in the LV and systemic artery under HF (A), and HF with the IABP operating at level 1 (sf=5%, B), 2 (sf=10%, C), 3 (sf=15%, D), 4 (sf=20%, E), and 5 (sf=25%, F). The mean arterial pressure under HF was 125 mmHg, which matches that in congestive heart failure. As the level of IABP treatment was increased, the aortic pressure during diastole, which is related to coronary perfusion, increased and the left ventricular peak pressure during systole, which is related to the after-load, decreased.
Fig. 3 shows the calculated transmural distribution of the contractile ATP consumption rates at end systole and end diastole for HF without an IABP (Fig. 3A), and the model with the IABP operating at level 5 (Fig. 3B). The ATP consumption rate changed significantly in the left ventricle rather than the right during systole. The contractile ATP consumption decreased significantly in the presence of the IABP. Compared with the case involving HF without an IABP, the contractile ATP consumption decreased by 12% with IABP therapy (see the ATP consumption rate in Table 2).
Fig. 4 shows the calculated transmural distribution of the fiber strain at end systole and end diastole for HF without an IABP (Fig. 4A) and the model with the IABP operating at level 5 (Fig. 4B). The fiber strain decreased significantly in the presence of the IABP.
Fig. 5 shows the pressure-volume curves for the six cases studied: HF without IABP therapy, and HF with the IABP at levels 1 to 5. The pressure-volume curves shifted to the left and downward. The end-diastolic volume decreased by 8 mL (Fig. 6); this indicates that the IABP results in ventricular volume unloading by reducing the end-diastolic volume.
Table 2 shows the hemodynamic responses for HF and HF with the IABP at levels 1 to 5. While the cardiac output and stroke volume increased, the stroke work and ATP consumption decreased with IABP therapy due to the pressure unloading effect. The decreased left ventricular peak pressure also indicates the pressure unloading effect during IABP therapy.
This study investigated the changes in the contractile energy consumption of the myocardium with IABP with different inflation volumes using a cardiac electromechanical model with a lumped model of the circulatory system and IABP function. To generalize the patterns of inflation and deflation of the IABP, a harmonic waveform was used for the time-varying compliance of the systemic arteries. According to the IABP inflation volume, the arterial pressure increased during diastole, while the ventricular peak pressure decreased, which indicates improved coronary perfusion and ventricular unloading. To estimate the change in the ventricular workload directly, the contractile ATP consumption of the ventricle was computed (Fig. 3). With greater inflation, the ATP consumption decreased significantly. The transmural distribution of the mechanical strain was also compared between HF ventricles and HF ventricles with an IABP at level 5 in order to estimate the volume-unloading effect of the IABP (Fig. 4). The fiber strain also decreased significantly with an IABP. The pressure-volume diagrams in Fig. 5 were consistent with the fiber strain distributions.
The contractile ATP consumption and fiber strain varies transmurally. In general, the sarcomere length affects the density of effective cross-bridges and therefore the contractile force; this is known as the Frank-Starling law of the heart. A longer sarcomere length results in greater contractile ATP consumption because of the greater density of the attached cross bridges in myocytes, as expressed in Eq. 3 in the Appendix. Previously, we demonstrated that myofibers in the endocardium have greater end-diastolic lengths than myofibers in the mid-myocardium and epicardium (6). Fig. 4 also shows that the end-diastolic strain and ATP consumption rate of the endocardial sarcomere are greatest. This explains why the contractile ATP consumption by the endocardium was greater than that at the mid-wall or epicardium (Figs. 3 and 4). The cardiac output, stroke volume, stroke work, MAP, ATP consumption rate, and ejection fraction are all increased with IABP treatment (Table 2). Interesting point is that although the ATP consumption is decreased significantly, stroke work is decreased only slightly, which indicates that the IABP helps the failed ventricle to pump blood efficiently.
There are several limitations. We used an image-based electromechanical model of failing canine ventricles, in which the ventricular geometry, fiber architecture, and hemodynamics differ from those of human ventricles. We also considered only the contractile ATP consumption of the myocardium, although ATP is used in other ways, such as by sarcoplasmic/endoplasmic reticulum calcium ATPase (SERCA), plasma membrane Ca-ATPase (PMCA), and Na pumps. Their functions are changed at failing heart condition (13-15). Finally, we did not considered remodeled electrical activities of myocardial cell in failing heart. We assumed that the electrical conductivity changes in tissue level encapsulated those changes in cellular level. However, these limitations should not alter the main findings of this study greatly.
We conclude that IABP improves coronary perfusion as well as pressure unloading and volume unloading and also helps the failed ventricle to pump blood more efficiently.
This work was supported by the NRL (National Research Lab) program of National Research Foundation of Korea (No. ROA-2008-000-20127-0) and NRF grant (No. 2012-0008820).
The authors have no conflicts of interest to disclose.
The contractile energy consumption of the myocardium was quantified by calculating the contractile ATP consumption in the myofilament model of Rice et al. (7). The contractile ATP consumption rate, E, per unit volume was calculated as a function of the ATP-consuming cross-bridge detachment rate (gxbT) and the single overlap fraction of thick filaments (SOVFThick),
where gxbT, which indicates the ATP-consuming detachment transition rate, and SOVFThick, which indicates the single-overlap fraction of the thick filament, were functions taken from the original myofilament model of Rice et al. (7).
Then, gxbT can be derived as:
The term xbmodspecies scales all cross-bridge cycling rates to account for species-based differences (0.2 for canines). The temperature dependence of gxbT (QgxbT) is set to a default Q10 value of 6.25. The term gxbT had strain dependence in the rate modifier gxbmd, defined as:
where δp and δn set the effects of strain for positive and negative shortening velocities, respectively; xo is the distortion of the cross-bridge head; xXBPostR is the mean distortion of states xXBPostR, which is the probability of the transition to the post-rotated force-generating status. The ATP-consuming detachment transition rate gxbT is proportional to gxbmd, which depends on the strain of the cross bridges, and in turn, this strain depends on the velocity of contraction. See Rice et al. (7) for details. Finally, we calculated the ATP consumption distribution by integrating the ATP consumption rate with time for 600 ms and the ATP consumption of the entire ventricle by spatially integrating the local ATP consumption through the entire ventricular volume.
|1.||Lin CY,Galysh FT,Ho KJ,Patel AS. Response to single-segment intraaortic balloon pumping as related to aortic complianceAnn Thorac SurgYear: 1972134684765022424|
|2.||Clark JW Jr,Kane GR,Bourland HM. On the feasibility of closed-loop control of intra-aortic balloon pumpingIEEE Trans Biomed EngYear: 1973204044124754312|
|3.||Barnea O,Moore TW,Dubin SE,Jaron D. Cardiac energy considerations during intraaortic balloon pumpingIEEE Trans Biomed EngYear: 1990371701812312142|
|4.||Suga H,Hisano R,Goto Y,Yamada O,Igarashi Y. Effect of positive inotropic agents on the relation between oxygen consumption and systolic pressure volume area in canine left ventricleCirc ResYear: 1983533063186883652|
|5.||Lim KM,Kim IS,Choi SW,Min BG,Won YS,Kim HY,Shim EB. Computational analysis of the effect of the type of LVAD flow on coronary perfusion and ventricular afterloadJ Physiol SciYear: 20095930731619387788|
|6.||Lim KM,Constantino J,Gurev V,Zhu R,Shim EB,Trayanova NA. Comparison of the effects of continuous and pulsatile left ventricular-assist devices on ventricular unloading using a cardiac electromechanics modelJ Physiol SciYear: 201262111922076841|
|7.||Rice JJ,Wang F,Bers DM,de Tombe PP. Approximate model of cooperative activation and crossbridge cycling in cardiac muscle using ordinary differential equationsBiophys JYear: 2008952368239018234826|
|8.||Gurev V,Lee T,Constantino J,Arevalo H,Trayanova NA. Models of cardiac electromechanics based on individual hearts imaging data: image-based electromechanical models of the heartBiomech Model MechanobiolYear: 20111029530620589408|
|9.||Fox JJ,McHarg JL,Gilmour RF Jr. Ionic mechanism of electrical alternansAm J Physiol Heart Circ PhysiolYear: 2002282H516H53011788399|
|10.||Helm RH,Byrne M,Helm PA,Daya SK,Osman NF,Tunin R,Halperin HR,Berger RD,Kass DA,Lardo AC. Three-dimensional mapping of optimal left ventricular pacing site for cardiac resynchronizationCirculationYear: 200711595396117296857|
|11.||Wu Y,Bell SP,Trombitas K,Witt CC,Labeit S,LeWinter MM,Granzier H. Changes in titin isoform expression in pacing-induced cardiac failure give rise to increased passive muscle stiffnessCirculationYear: 20021061384138912221057|
|12.||O'Rourke B,Kass DA,Tomaselli GF,Kaab S,Tunin R,Marban E. Mechanisms of altered excitation-contraction coupling in canine tachycardia-induced heart failure, I: experimental studiesCirc ResYear: 19998456257010082478|
|13.||Dally S,Bredoux R,Corvazier E,Andersen JP,Clausen JD,Dode L,Fanchaouy M,Gelebart P,Monceau V,Del Monte F,et al. Ca2+-ATPases in non-failing and failing heart: evidence for a novel cardiac sarco/endoplasmic reticulum Ca2+-ATPase 2 isoform (SERCA2c)Biochem JYear: 200639524925816402920|
|14.||Bassani JW,Bassani RA. SERCA upregulation: breaking the positive feedback in heart failure?Cardiovasc ResYear: 20056758158216018992|
|15.||Muller-Ehmsen J,McDonough AA,Farley RA,Schwinger RH. Sodium pump isoform expression in heart failure: implication for treatmentBasic Res CardiolYear: 200297Suppl 1I25I3012479230|
Keywords: Intra-Aortic Balloon Pump, Cardiac Electromechanics Model, ATP Consumption, Stroke Work.
Previous Document: Peripheral artery disease in korean patients undergoing percutaneous coronary intervention: prevalen...
Next Document: The role of urinary liver-type fatty acid-binding protein in critically ill patients.