Computational and Mathematical Methods in Medicine
Volume 2018 (2018), Article ID 4253928, 12 pages
https://doi.org/10.1155/2018/4253928

Research Article

Computational Prediction of the Combined Effect of CRT and LVAD on Cardiac Electromechanical Delay in LBBB and RBBB

Aulia K. Heikhmakhtiar and Ki M. Lim

Department of IT Convergence Engineering, Kumoh National Institute of Technology, Gumi 39177, Republic of Korea

Correspondence should be addressed to Ki M. Lim; kmlim@kumoh.ac.kr

Received 17 May 2018; Revised 3 September 2018; Accepted 25 October 2018; Published 14 November 2018

Guest Editor: Ka L. Man

Copyright © 2018 Aulia K. Heikhmakhtiar and Ki M. Lim. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Two case reports showed that the combination of CRT and LVAD benefits the end-stage heart failure patients with prolonged QRS interval significantly. In one of the reports, the patient had the LVAD removed due to the recovery of the heart function. However, the quantification of the combined devices has yet to be conducted. This study aimed at computationally predicting the effects of CRT-only or combined with LVAD on electromechanical behaviour in the failing ventricle with left bundle branch blocked (LBBB) and right bundle branch blocked (RBBB) conditions. The subjects are normal sinus rhythm, LBBB, RBBB, LBBB with CRT-only, RBBB with CRT-only, LBBB with CRT + LVAD, and RBBB with CRT + LVAD. The results showed that the CRT-only shortened the total electrical activation time (EAT) in the LBBB and RBBB conditions by 20.2% and 17.1%, respectively. The CRT-only reduced the total mechanical activation time (MAT) and electromechanical delay (EMD) of the ventricle under LBBB by 21.3% and 10.1%, respectively. Furthermore, the CRT-only reduced the contractile adenosine triphosphate (ATP) consumption by 5%, increased left ventricular (LV) pressure by 6%, and enhanced cardiac output (CO) by 0.2 L/min under LBBB condition. However, CRT-only barely affects the ventricle under RBBB condition. Under the LBBB condition, CRT + LVAD increased LV pressure and CO by 10.5% and by 0.9 L/min, respectively. CRT + LVAD reduced ATP consumption by 15%, shortened the MAT by 23.4%, and shortened the EMD by 15.2%. In conclusion, we computationally predicted and quantified that the CRT + LVAD implementation is superior to CRT-only implementation particularly in HF with LBBB condition.

1. Introduction

Heart failure (HF) plays a major role in the number of death worldwide [1]. Thus, the study of heart diseases including cardiac arrhythmia, which progressively leads to HF condition, is very important. Two types of cardiac therapy devices are commonly used to treat patients with cardiac disease: cardiac resynchronization therapy (CRT) and left ventricular assist device (LVAD). CRT is considered as a valuable treatment for patients with dyssynchrony HF with QRS interval >120 ms and left ventricular ejection fraction (EF) <35% [2, 3]. A number of studies showed significant benefits of using CRT. CRT synchronizes systolic function [4–6], restores heart structure [7, 8], and improves symptoms and the quality of life of the patients identified by the improvement of exercise endurance [3, 8–12]. However, 30% of patients with HF failed to benefit from the CRT [9]. Hu et al. [13] described that three major responses could be obtained from the CRT treatment: resynchronization of the intraventricular contraction, efficient ventricular preloading by a properly timed atrial contraction, and reduction of mitral regurgitation. One of the study findings was that resynchronization of the intraventricular contraction itself did not necessarily lead to stroke work improvement. However, the synchronization of the atrioventricular firing time was essential.

LVAD supports the ventricular pumping via mechanical unloading for weakened heart. LVAD was initially used as a bridge to transplantation for patients with end-stage disease [14, 15]. However, because the availability of heart donors is very limited, LVAD is now used as destination therapy as it lasts for years [16–18]. LVAD reverses the damaged ventricle and recovers myocardial functionalities by repairing left ventricular (LV) mass [19, 20], heart chamber size [20, 21], improves mitral filling [22], and induces cardiomyocyte hypertrophy regression [21, 22]. The use of LVAD also increased the quality of life of the patients [23–25]. Previously, we observed the effect of LVAD on the electromechanical delay (EMD) under mild and severe HF conditions [26]. The results showed that the LVAD not only increased the EF but also shortened the EMD under the mild HF and even better under severe HF conditions. The comparison study of the symptomatic relief of CRT and LVAD by Delgado et al. concluded that the use of both devices could synergistically improve cardiac functions for severe HF treatment [27].

In 2011, a case report article stated that the CRT supported the restoration of end-stage HF in a 15-year-old patient who underwent LVAD implantation [28]. The CRT shortened the septal to posterior wall motion delay from 146 ms to 104 ms, and overall, it backed up the hemodynamic improvements. After such improvement, the patient successfully had the LVAD removed. In another case report article, the CRT and LVAD cooperatively restored the cardiac functions of a 62-year-old patient who had cardiogenic shock and left bundle branch block (LBBB) [29]. The report stated that the combination of CRT and LVAD performed a profound treatment on the weakening heart under LBBB condition. It recovered the heart hypertrophy, and the EF was increased. Based on these reports, we conducted a computational modelling which combined the CRT and LVAD to the failing ventricle with LBBB and RBBB in order to understand the mechanism of the two devices combined to provide such improvement to the heart.

To the best of our knowledge, a computational study that quantifies the effect of combined CRT and LVAD has yet to be conducted. This study aims at predicting, computationally with an electromechanical failing ventricles model, the effects of CRT alone and the combination of CRT and LVAD treatment in patients with LBBB or right bundle branch block (RBBB). We used a well-developed electromechanical ventricular model which had been used to observe different heart conditions, including the prolonged EMD in the dyssynchrony HF condition [30], mechanoelectrical feedback on scroll wave stability [31], and spontaneous arrhythmia in acute regional ischemia [32]. Recently, by using similar electromechanical coupling model, our group revealed the effect of KCNQ1 S140G mutation on arrhythmogenesis and pumping performance [33]. We used computational methods to overcome the measurement limitations and risks of the experimental study. We analyzed seven HF diseases and therapies: (i) normal sinus rhythm, (ii) LBBB, (iii) LBBB with CRT alone (LBBB + CRT), (iv) LBBB with CRT and LVAD combination (LBBB + CRT + LVAD), (v) RBBB, (vi) RBBB with CRT alone (RBBB + CRT), and (vii) RBBB with combined CRT and LVAD (RBBB + CRT + LVAD) conditions.

2. Materials and Methods

We followed an existing well-developed electromechanical failing ventricle model with the fibers and laminar sheet structures based on diffusion tensor magnetic resonance imaging [34–37]. In this study, our electromechanical model was coupled with circulatory systems and LVAD models similarly as Lim et al. with additional Purkinje networks compartment to simulate LBBB and RBBB conditions [38]. The electromechanical model consisted of electrophysiological and myofilament dynamics model coupled by Ca²⁺ transient. The LVAD function included the circulatory systems, which connected to the 3D ventricular mechanics. Figure 1 shows a full schematic of the system we used in this study.

Figure 1: Schematics of the electrical and mechanical elements coupled with transient calcium. Electrical element: it represents the currents, pumps, and exchanger of the ten Tusscher ionic model as explained in equation (1). Mechanical element: a schematic diagram of the finite-element ventricular mechanical model coupled with the circulatory and LVAD models. PRV, RV pressure; VRV, RV volume; PLV, LV pressure; VLV, LV volume; RPA, pulmonary artery resistance; CPA, pulmonary artery compliance; RPV, pulmonary vein resistance; CPV, pulmonary vein compliance; RMI, mitral valve resistance; CLA, left atrium compliance; RAO, aortic valve resistance; RSA, systemic artery resistance; CSA, systemic artery compliance; RSV, systemic vein resistance; CSV, systemic vein compliance; RTR tricuspid valve resistance; CRA, right atrium compliance; RPU, pulmonary valve resistance.

The electrical component was a failing ventricular mesh with 241,725 nodes and 1,298,751 tetrahedron elements. The electrical mesh has the characteristics of realistic heart compartment, which consists of endocardial, midmyocardial, and epicardial cells following the ten Tusscher et al. human ventricular tissue model [39]. The line mesh type representing Purkinje networks was mapped onto a 3D ventricle chamber as well at the endocardial region (Figure 2(c)).

(a)

(b)

(c)

Figure 2: (a) The heterogeneous mesh that has the characteristics of endocardial, midmyocardial, and epicardial cells. (b) The CRT pacing site for the LBBB heart was placed on the LV free-wall, while that for the RBBB heart was placed inside the RV at the bottom of the septum. (c) The electrical propagation by the Purkinje network in the sinus, LBBB, and RBBB conditions. The electrical activation by the Purkinje networks is indicated in red. In sinus pacing, the Purkinje network delivers the electrical signal to its terminals. In the LBBB condition, the Purkinje network delivered the signal only to the right network in the RV area. In the RBBB condition, the Purkinje network delivered the signal only to the left network in the LV area. CRT, cardiac resynchronization therapy; LBBB, left bundle branch block; LV, left ventricular; RBBB, right bundle branch block; RV, right ventricular.

The Purkinje networks induced the electrical activation sequences of sinus rhythm (normal), LBBB, and RBBB (Figure 2(c)). The CRT pacing site of the LBBB was placed at the LV free-wall, while the CRT pacing site of the RBBB was placed at the RV endocardial apex as shown in Figure 2(b). The electrical stimulation was first induced at the root node of Purkinje fiber model, which propagates to the terminals, hence stimulating the ventricular tissue. The electrical propagation in the Purkinje can be described by solving a one-dimensional wave propagation equation and triggered the ventricular activation [40].

The electrical propagation signal represents an ion exchange across the myocyte as described by ten Tusscher et al. [39]. The electrophysiological phenomenon for the single cell can be described as follows:

where represents the voltage difference of intracellular and extracellular, represents time, represents the sum of all transmembrane ionic currents, represents the current if an external stimulus is applied, and represents the cell capacitance per unit of surface area. consists of major ionic current as follows:

where is fast inward Na⁺ current, is inward rectifier current, is transient outward K⁺ current, is rapid delayed rectifier K⁺ current, is slow delayed rectifier K⁺ current, is L-type inward Ca²⁺ current, is Na⁺/Ca²⁺ exchanger current, is the Na⁺/K⁺ pump current, and are sarcoplasmic plateau Ca²⁺ and K⁺ currents, is background Ca²⁺ current, and is background Na⁺ current.

To represent the electrical propagation through the conduction in 3D spatial, the cardiac tissue in this case could be described by the combination of the Equation (1) with the cellular resistivity () and surface-to-volume ratio () in x, y, and z directions, respectively. This phenomenon can be described by the following partial differential equation:

The mathematical equation for calcium dynamics for coupling the electromechanical simulation was also described by ten Tusscher et al. [39]:

where represents the calcium released from the sarcoplasmic reticulum (SR) into the cell. represents the calcium pumping to restore the calcium again back to the SR. represents the calcium-induced calcium release current, is the activation gate of , which was the same with . is the total calcium inside the cell, which consists of , the calcium buffer inside the cell, and , the free calcium inside the cell. Accordingly, the is the sum of calcium in the SR, which includes , the calcium buffer, and , the free calcium inside the SR.

A Ca²⁺ transient serves as an input to the cell myofilament model representing the generation of active tension within each myocyte in which an ODE set and multiple algebraic equations describe Ca²⁺ binding to troponin C, cooperatively between regulatory proteins, and cross-bridge cycling. Ventricular contraction results from the active tension generation represented by the myofilament dynamics model described by Rice et al. [41]. Ventricular deformation is described by the equations of passive cardiac mechanics, with the myocardium being orthotropically hyperelastic, and nearly incompressible material with passive properties defined by an exponential strain-energy function. Simultaneous solutions of the myofilament model equations to those representing passive cardiac mechanics on the finite-element mesh constitute cardiac contraction. Considering the isometric contraction, we assumed that the sarcomere length (SL) is 0 at the initial value, dSL/dt = 0. To measure the isotonic contraction, the SL is solved using the following ordinary differential equation (ODE):

the viscosity and mass are described in Figure 1 at the mechanical compartment, and is the sum of the normalized force integrated toward time:

where the is the active force and the is the passive force. is the constant force at the resting length of the initial sarcomere length, and is the force during the isotonic or isometric contraction. Thus, is expressed as follows:

We calculated the ATP consumption by integrating the myofilament model in one cycle proposed by Rice et al. from each node [41]. In the myofilament model, the ATP consumption rate (E) is the outcome of cross-bridge detachment rate () and the single overlap fraction of thick filaments (SOVFThick):

To construct an integrated model of an LVAD-implanted cardiovascular system, we added a compartment of LVAD function between LV and systemic arteries in the circulatory system based on Kerckhoffs et al. [42] as described previously by Lim et al. [38]. Briefly, the LVAD component was modelled as a flow generator with a specific mean flow rate of 3 L/min. Constant-flow conditions were used to simulate the continuous LVAD with the inlet at the LV apex and the outlet at the ascending aorta.

For the simulation protocol, first we simulated the electrical model with the Purkinje delivering the signal representing sinus rhythm, LBBB, or RBBB until the steady state was achieved. We set the conduction velocity by 60 cm/s and the basic cycle length by 600 ms. Here, we incorporated the dyssynchrony HF conditions (LBBB and RBBB) using CRT-alone or combined CRT and LVAD. We then used Gaussian point for the interpolation of transient Ca²⁺ from the electrical simulation as the input to the mechanical simulation. To model the HF condition, we multiplied the constant of passive scaling in the strain-energy function by 5 times to stiffen the myocardium. The mechanical model was simulated for 20 seconds to reach the steady state. We compared the ATP consumption rates and tension activation during end-systolic volume (ESV) and the strain during end-diastolic volume (EDV) by integrating the information of them from each node. EMD was defined as the time interval between mechanical activation time (MAT) and electrical activation time (EAT). MAT was identified when the local strain was shortened to 10% before its maximum, while EAT was identified when the myocyte started to depolarize, which in our case was −50 mV [30].

3. Results and Discussion

3.1. Electrophysiological Simulation Results

Figure 3 shows the membrane potential propagation in the normal sinus rhythm (sinus), LBBB, LBBB + CRT, LBBB + CRT + LVAD, RBBB, RBBB + CRT, and RBBB + CRT + LVAD conditions. The electrical wavelength almost covered the whole ventricle tissue at 100 ms in the sinus condition, and the EAT in the sinus condition was 120 ms (Table 1). The EAT in the LBBB condition was 173 ms, which was the longest among other cases. CRT shortened the EAT in the LBBB condition to 138 ms, which was close to that in the sinus condition. In the RBBB heart, the EAT was 164 ms. CRT shortened the EAT in the RBBB heart to 136 ms, which was close to the sinus condition as well.

Figure 3: The membrane potential propagation of the sinus pacing, LBBB, LBBB + CRT, RBBB, and RBBB + CRT models. For the LBBB + CRT + LVAD and RBBB + CRT + LVAD models, we used the same electrical activation sequence as the LBBB + CRT and RBBB + CRT models, respectively, because the LVAD was incorporated into the mechanical computation. LBBB, left bundle branch block; CRT, cardiac resynchronization therapy; RBBB, right bundle branch block; LVAD, left ventricular assist device.

Table 1: Hemodynamic responses under normal, LBBB, LBBB + CRT, LBBB + CRT + LVAD, RBBB, RBBB + CRT, and RBBB + CRT + LVAD models.

Condition
EDV (mL)
ESV (mL)
CO (L/min)
EF (%)
Longest EAT (ms)
Average MAT (ms)
Average EMD (ms)

Normal sinus rhythm
88
55
3.4
38
120
157
78

LBBB
90
60
3
33.4
173
188
79

LBBB + CRT
89
57
3.2
36
138
148
71

LBBB + CRT + LVAD
—
—
3.9
—
138
144
67

RBBB
87
54
3.4
38
164
162
80

RBBB + CRT
88
55
3.4
38
136
157
81

RBBB + CRT + LVAD
—
—
4
—
136
155
79

Figure 4 shows the calcium activation sequences. The calcium activation sequence followed the electrical activation sequence for each case. Compared to the membrane potential, which deactivated after 450 ms, the Ca²⁺ deactivation occurred at 250 ms. CRT fastened the Ca²⁺ activation throughout the ventricle in the LBBB condition but insignificantly in the RBBB condition due to the CRT pacing site.

Figure 4: The Ca²⁺ propagation sequence following the membrane potential activation sequence in all cases. LBBB, left bundle branch block; CRT, cardiac resynchronization therapy; LVAD, left ventricular assist device; RBBB, right bundle branch block.

3.2. Cardiac Mechanics Simulation Results

Figure 5 shows the 3D ATP consumption rate, tension, and strain transmural distribution in the sinus, LBBB, LBBB + CRT, LBBB + CRT + LVAD, RBBB, RBBB + CRT, and RBBB + CRT + LVAD conditions. We pick one node at the LV free-wall as representative to compare the values of them. Overall, the ATP consumption rate and tension in the LBBB condition had the highest values, 3.3 and 3.5 times larger than those in the normal sinus condition, respectively. However, the CRT decreased the ATP consumption rate and tension in the LBBB condition to be 6% and 7% lower than those in the normal sinus condition. Furthermore, CRT and LVAD reduced the ATP consumption rate and tension in the LBBB condition to 10% and 12% lower than those in the normal sinus condition, respectively. The ATP consumption rate and tension in the RBBB condition were 8% and 9% higher than those in the control condition. The CRT reduced the ATP consumption rate and tension by 4% and 5%, respectively. CRT and LVAD reduced the ATP consumption rate and tension by 9% and 11% in the RBBB condition, respectively. The strain under the LBBB condition was 80 times larger than that in the normal condition (notified by major red region in the LV free-wall). The CRT reduced the strain in the LBBB 56 times lower than that in the normal sinus condition. CRT significantly restored the total strain under the LBBB condition with the pacing site at the LV free-wall. However, under LBBB + CRT + LVAD, the total strain was only eight times lower than that in the sinus condition. The continuous LVAD altered the overall strain in the ventricles. The total strain in the RBBB condition was 84 times larger than that in the normal condition. In the RBBB + CRT model, the strain activation was 122 times larger. In the RBBB + CRT + LVAD model, the strain activation was 23 times lower than that in the normal condition.

Figure 5: The 3D transmural distribution of ATP, tension, and strain in all cases. The snapshot for the ATP and tension was taken at the end-systolic volume, while the strain snapshot was taken at the end-diastolic volume time. ATP, adenosine triphosphate; LBBB, left bundle branch block; CRT, cardiac resynchronization therapy; LVAD, left ventricular assist device; RBBB, right bundle branch block.

Figure 6(a) shows the LV pressure-volume (PV) loop diagram in the sinus, LBBB, LBBB + CRT, LBBB + CRT + LVAD, RBBB, RBBB + CRT, and RBBB + CRT + LVAD conditions. The LV PV-loop in the normal sinus rhythm condition was the same as those in the RBBB and RBBB + CRT conditions. The EDV in the three conditions was 88 mL, ESV was 54.5 mL, and EF was 38%. The LBBB condition had the highest EDV and ESV of 90 mL and 60 mL, respectively. The stroke volume (SV) and EF were 30 mL and 33.4%, respectively. CRT increased the EF up to 36% in the LBBB condition. The EDV, ESV, and SV in the LBBB + CRT condition were 89 mL, 57 mL, and 32 mL, respectively. We did not quantify the EF in models in which LVAD model was incorporated because the PV-loop of the LV was altered due to the LVAD implementation. However, the CO was presented to compare the effects of CRT-alone and combined CRT and LVAD. The complete EDV, ESV, CO, EF, MAT, and EMD data in all conditions are provided in Table 1.

(a)

(b)

(c)

(d)

Figure 6: (a) LV pressure-volume loop, (b) LV pressure and aortic pressure, (c) ATP consumption rate, and (d) cardiac output of all cases. LV, left ventricular; LBBB, left bundle branch block; CRT, cardiac resynchronization therapy; LVAD, left ventricular assist device; RBBB, right bundle branch block; ATP, adenosine triphosphate.

Figure 6(b) shows the LV and systemic artery pressures coloured in black and red lines, respectively. In the normal sinus condition, the peak LV pressure was 148 mmHg. The LV pressure peaked at 358 ms, and the aortic valve was opened for 96 ms. In the LBBB condition, the LV peak pressure was 132.5 mmHg. The LV peak pressure time was prolonged to 407 ms, and the aortic valve was opened for 100 ms. The LV peak pressure in the LBBB + CRT condition was increased to 141 mmHg. The LV peak pressure time was shortened to 364 ms, and the aortic valve was opened for 96 ms. In the LBBB + CRT + LVAD condition, the peak LV pressure was 148 mmHg, the same as that in the normal condition, the LV peak pressure time was 363 ms, and the aortic valve was opened for 62 ms. The combination of CRT and LVAD restored the LV peak pressure and LV peak pressure time in the LBBB condition better than CRT alone. In the RBBB and RBBB + CRT conditions, the LV peak pressures and LV peak pressure time were the same at 148 mmHg and 353 ms, respectively, and the aortic valve was opened for 96 ms, the same as that in the normal condition. The reason was that the electrical activation of the LV was not altered in the RBBB condition. In the RBBB + CRT + LVAD condition, the peak LV pressure and LV peak pressure time were increased to 155 mmHg and 350 ms, respectively, and the aortic valve was opened for 63 ms. This finding shows that the combined CRT and LVAD increased the LV pressure and systemic artery more than the CRT-only implementation.

Figure 6(c) shows the overall ATP consumption rate from one cycle of steady state. In the normal sinus condition, the ATP consumption rate was 93 s⁻¹. The ATP consumption rate in LBBB condition was the highest at 100 s⁻¹. CRT reduced the ATP consumption rate in the LBBB condition by 5%, while CRT and LVAD reduced the ATP consumption rate by 15%. In the RBBB and RBBB + CRT conditions, the ATP consumption rates were the same as that in the normal condition, 93 s⁻¹. However, with CRT and LVAD support, the ATP consumption rate in the RBBB condition was decreased by 16% to 84 s⁻¹.

Figure 6(d) shows the CO of the seven subjects. CRT increased the CO slightly in the LBBB condition but not in the RBBB condition. With the combination of CRT and LVAD, the CO of the LBBB and RBBB conditions was significantly increased (Table 1).

Figure 7(a) shows the 3D contour of EAT, MAT, and EMD, while Figure 7(b) shows the MAT and EMD of all cases. As shown in Figure 7(a), the activation sequence of MAT and EMD was identical to the EAT. The MAT and EMD values in the normal condition were 157 and 78 ms, respectively. In the LBBB condition, the MAT and EMD values were the greatest at 188 ms and 79 ms, respectively. CRT shortened the MAT and EMD to 148 ms and 71 ms in the LBBB condition. Furthermore, the MAT and EMD in the LBBB condition were further shortened with the combined CRT and LVAD to 144 and 67 ms, respectively. This finding showed that combination of CRT and LVAD performed better than CRT alone despite having the same EAT. In the RBBB condition, the MAT and EMD values were 162 and 80 ms, respectively. In the RBBB + CRT condition, the mean MAT and EMD values were reduced slightly to 157 and 81 ms, respectively. In the RBBB + CRT + LVAD condition, the mean MAT and EMD values were reduced to 155 and 79 ms (close to the control condition), respectively. In the RBBB condition, CRT alone only slightly affected the electrical and mechanical responses. However, CRT + LVAD increased CO in the RBBB condition significantly as expected.

(a)

(b)

Figure 7: (a) Three-dimensional contour distribution of the EAT, MAT, and EMD and (b) the mean MAT and EMD values in all cases. EAT, electrical activation time; MAT, mechanical activation time; EMD, electromechanical delay; LBBB, left bundle branch block; CRT, cardiac resynchronization therapy; LVAD, left ventricular assist device; RBBB, right bundle branch block.

In general, we performed a simulation and analyzed seven cardiac diseases and therapy conditions by using an electromechanical ventricular model incorporated with a circulatory systems, CRT, and LVAD models. The models including normal sinus rhythm, LBBB, LBBB + CRT, LBBB + CRT + LVAD, RBBB, RBBB + CRT, and RBBB + CRT + LVAD. The major findings of this study are as follows:

(1)CRT shortened the longest EAT by 20.2% in the LBBB condition and 17.1% in the RBBB condition. CRT shortened EMD by 10.1% in the LBBB condition but did not shorten EMD in the RBBB condition.(2)Combination of CRT and LVAD treatment shortened EMD more than CRT alone in the LBBB condition (15.2%) and shortened EMD in the RBBB condition (1.3%).(3)CRT reduced the ATP consumption by 5% as well as the tension and strain in the LBBB condition. CRT also slightly increased the LV peak pressure by 6% and increased the CO by 0.2 L/min in the LBBB condition. However, CRT-alone did not affect these mechanical responses in the RBBB condition.(4)Combination of CRT and LVAD reduced the ATP consumption by 15% in the LBBB condition and by 16% in the RBBB condition. It also increased the LV pressure by 10.5% in the LBBB condition and 5.7% in the RBBB condition as well as the CO up to 4 L/min for both conditions, a degree greater than that in the control condition.

The pacing site at the tissue greatly affected synchronization. Placement at the LV free-wall showed faster activation throughout the chamber since the signal propagated evenly from the midway between the base and apex throughout the LV tissue. However, the pacing site at the RV endocardial apex showed a longer activation time for the electrical signal to propagate to the base.

In the RBBB condition, the electrical activation of the LV chamber was the same as normal sinus despite the electrical activation alteration in the RV. This is the major factor why the mechanical responses of the LV (LV PV-loop, LV pressure, and CO) were also the same as that in normal sinus rhythm. Even the CRT implementation in the RBBB condition did not affect the mechanical responses of the LV chamber. However, we observed significant improvement in CO in the RBBB + CRT + LVAD condition. CRT and LVAD also reduced the ATP consumption and tension and increased the LV and aortic pressure in the RBBB condition.

CRT-alone did not fully restore the mechanical responses under the LBBB condition to normal despite the resynchronization of the EAT, which shortened the MAT and EMD to back normal. On the contrary, the use of combined CRT and LVAD reduced the energy consumption (indicated by the ATP consumption rate), tension, increased the LV pressure, and produced CO by 18% more than that in the HF with sinus rhythm condition. In addition, the LVAD did not fully assist the blood distribution in the LBBB or RBBB condition.

The computational method allowed us to predict the electromechanical phenomenon in the dyssynchrony heart which underwent CRT and LVAD treatment. It is hardly possible to observe some parameters including ATP consumption, tension, MAT, and EMD from the patient’s heart of that condition in experimental procedure. This study quantified the electrical activation and hemodynamic responses in several dyssynchrony HF combined with CRT-only, and CRT and LVAD, which has never been conducted previously. As for the clinical impact, the results of this study can be used as reference to generally predict the effect of the combination of CRT and LVAD devices to the HF patients with LBBB. Though, there are some parameters need to be consider deeply by the cardiologist expert.

This computational study has several limitations that need to be addressed. The study did not follow the standard biventricular pacing method [43]. Instead, we used only one CRT pacing site: the RV endocardial apex (for the RBBB condition) or the LV lateral wall (for the LBBB condition). The LV pacing site was not placed at the optimal position as described before [44]. We did not observe or describe the RV mechanical responses in the RBBB condition; instead, we showed the LV responses only. We used LVAD instead of a right ventricular assist device to support the RBBB condition. Because of our limitations, we did not validate the results of our simulation with experimental data. In addition, we did not describe long-term effects such as recovery of cardiac functions by the combined CRT and LVAD, as previously described [28, 29].

4. Conclusions

In conclusion, although CRT-alone shortened the MAT and EMD to more than normal in the LBBB condition, the mechanical responses in the LBBB condition were not restored to normal. The combined CRT and LVAD shortened the MAT and EMD more than CRT-alone, restored the hemodynamic, and produced a greater CO than normal in the LBBB and RBBB conditions. Using the combined system, LVAD contributed to the MAT reduction by mechanical unloading, shortened the EMD, reduced ATP consumption, and reduced tension, which contributes to the recovery of the heart shape and function. In short, we computationally predicted and quantified that the CRT + LVAD implementation is superior to CRT-only implementation particularly in HF with the LBBB condition.

Data Availability

The methods and results data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research was partially supported by the MSIT (Ministry of Science, ICT), under the ITRC support program (IITP-2018-2014-0-00639) and and Global IT Talent Support Program (IITP-2017-0-01811) supervised by the IITP, and NRF under basic engineering research project (2016R1D1A1B0101440) and the EDISON (NRF-2011-0020576) Programs.

References

D. Lloyd-Jones, R. Adams, M. Carnethon et al., “Heart disease and stroke statistics—2009 update: a report from the American Heart Association Statistics Committee and Stroke Statistics Subcommittee,” Circulation, vol. 119, no. 3, pp. e21–e181, 2009.

C. Linde, C. Leclercq, S. Rex et al., “Long-term benefits of biventricular pacing in congestive heart failure: results from the MUltisite STimulation in cardiomyopathy (MUSTIC) study,” Journal of the American College of Cardiology, vol. 40, no. 1, pp. 111–118, 2002.

J. B. Young, W. T. Abraham, A. L. Smith et al., “Combined cardiac resynchronization and implantable cardioversion defibrillation in advanced chronic heart failure: the MIRACLE ICD Trial,” JAMA, vol. 289, no. 20, pp. 2685–2694, 2003.

M. Kawaguchi, T. Murabayashi, B. J. Fetics et al., “Quantitation of basal dyssynchrony and acute resynchronization from left or biventricular pacing by novel echo-contrast variability imaging,” Journal of the American College of Cardiology, vol. 39, no. 12, pp. 2052–2058, 2002.

G. S. Nelson, R. D. Berger, B. J. Fetics et al., “Left ventricular or biventricular pacing improves cardiac function at diminished energy cost in patients with dilated cardiomyopathy and left bundle-branch block,” Circulation, vol. 102, no. 25, pp. 3053–3059, 2000.

P. Søgaard, W. Y. Kim, H. K. Jensen et al., “Impact of acute biventricular pacing on left ventricular performance and volumes in patients with severe heart failure,” Cardiology, vol. 95, no. 4, pp. 173–182, 2001.

A. Auricchio, C. Stellbrink, S. Sack et al., “Long-term clinical effect of hemodynamically optimized cardiac resynchronization therapy in patients with heart failure and ventricular conduction delay,” Journal of the American College of Cardiology, vol. 39, no. 12, pp. 2026–2033, 2002.

M. G. S. J. Sutton, T. Plappert, K. E. Hilpisch et al., “Sustained reverse left ventricular structural remodeling with cardiac resynchronization at one year is a function of etiology quantitative doppler echocardiographic evidence from the Multicenter InSync Randomized Clinical Evaluation (MIRACLE),” Circulation, vol. 113, no. 2, pp. 266–272, 2006.

W. T. Abraham, W. G. Fisher, A. L. Smith et al., “Cardiac resynchronization in chronic heart failure,” New England Journal of Medicine, vol. 346, no. 24, pp. 1845–1853, 2002.

S. Cazeau, C. Leclercq, T. Lavergne et al., “Effects of multisite biventricular pacing in patients with heart failure and intraventricular conduction delay,” New England Journal of Medicine, vol. 344, no. 12, pp. 873–880, 2001.

J. G. Cleland, J.-C. Daubert, E. Erdmann et al., “The effect of cardiac resynchronization on morbidity and mortality in heart failure,” New England Journal of Medicine, vol. 352, no. 15, pp. 1539–1549, 2005.

A. J. Moss, W. J. Hall, D. S. Cannom et al., “Cardiac-resynchronization therapy for the prevention of heart-failure events,” New England Journal of Medicine, vol. 361, no. 14, pp. 1329–1338, 2009.

Y. Hu, V. Gurev, J. Constantino, and N. Trayanova, “Efficient preloading of the ventricles by a properly timed atrial contraction underlies stroke work improvement in the acute response to cardiac resynchronization therapy,” Heart Rhythm, vol. 10, no. 12, pp. 1800–1806, 2013.

R. John, F. Kamdar, K. Liao et al., “Improved survival and decreasing incidence of adverse events with the HeartMate II left ventricular assist device as bridge-to-transplant therapy,” Annals of Thoracic surgery, vol. 86, no. 4, pp. 1227–1235, 2008.

D. Casarotto, T. Bottio, A. Gambino, L. Testolin, and G. Gerosa, “The last to die is hope: prolonged mechanical circulatory support with a Novacor left ventricular assist device as a bridge to transplantation,” Journal of Thoracic and Cardiovascular Surgery, vol. 125, no. 2, pp. 417-418, 2003.

M. A. Daneshmand, K. Rajagopal, B. Lima et al., “Left ventricular assist device destination therapy versus extended criteria cardiac transplant,” Annals of Thoracic Surgery, vol. 89, no. 4, pp. 1205–1210, 2010.

E. J. Birks, “Left ventricular assist devices,” Heart, vol. 96, no. 1, pp. 63–71, 2010.

E. A. Rose, A. C. Gelijns, A. J. Moskowitz et al., “Long-term use of a left ventricular assist device for end-stage heart failure,” New England Journal of Medicine, vol. 345, no. 20, pp. 1435–1443, 2001.

G. T. Altemose, V. Gritsus, V. Jeevanandam, B. Goldman, and K. B. Margulies, “Altered myocardial phenotype after mechanical support in human beings with advanced cardiomyopathy,” Journal of heart and lung transplantation: The Official Publication of the International Society for Heart Transplantation, vol. 16, no. 7, pp. 765–773, 1997.

H. R. Levin, M. C. Oz, J. M. Chen et al., “Reversal of chronic ventricular dilation in patients with end-stage cardiomyopathy by prolonged mechanical unloading,” Circulation, vol. 91, no. 11, pp. 2717–2720, 1995.

A. Barbone, M. C. Oz, D. Burkhoff, and J. W. Holmes, “Normalized diastolic properties after left ventricular assist result from reverse remodeling of chamber geometry,” Circulation, vol. 104, no. 1, pp. I-229–I-232, 2001.

S. Nakatani, P. M. McCarthy, K. Kottke-Marchant et al., “Left ventricular echocardiographic and histologic changes: impact of chronic unloading by an implantable ventricular assist device,” Journal of the American College of Cardiology, vol. 27, no. 4, pp. 894–901, 1996.

K. L. Grady, P. Meyer, A. Mattea et al., “Improvement in quality of life outcomes 2 weeks after left ventricular assist device implantation,” Journal of Heart and Lung Transplantation, vol. 20, no. 6, pp. 657–669, 2001.

F. D. Pagani, L. W. Miller, S. D. Russell et al., “Extended mechanical circulatory support with a continuous-flow rotary left ventricular assist device,” Journal of the American College of Cardiology, vol. 54, no. 4, pp. 312–321, 2009.

M. S. Slaughter, J. G. Rogers, C. A. Milano et al., “Advanced heart failure treated with continuous-flow left ventricular assist device,” New England Journal of Medicine, vol. 361, no. 23, pp. 2241–2251, 2009.

A. K. Heikhmakhtiar, A. J. Ryu, E. B. Shim et al., “Influence of LVAD function on mechanical unloading and electromechanical delay: a simulation study,” Medical and Biological Engineering and Computing, vol. 56, no. 5, pp. 911–921, 2017.

R. M. Delgado and B. Radovancevic, “Symptomatic relief: left ventricular assist devices versus resynchronization therapy,” Heart Failure Clinics, vol. 3, no. 3, pp. 259–265, 2007.

J. Muratsu, M. Hara, I. Mizote et al., “The impact of cardiac resynchronization therapy in an end-stage heart failure patient with a left ventricular assist device as a bridge to recovery,” International Heart Journal, vol. 52, no. 4, pp. 246-247, 2011.

H. Keilegavlen, J. E. Nordrehaug, S. Faerestrand et al., “Treatment of cardiogenic shock with left ventricular assist device combined with cardiac resynchronization therapy: a case report,” Journal of Cardiothoracic Surgery, vol. 5, no. 1, p. 1, 2010.

J. Constantino, Y. Hu, A. C. Lardo, and N. A. Trayanova, “Mechanistic insight into prolonged electromechanical delay in dyssynchronous heart failure: a computational study,” American Journal of Physiology-Heart and Circulatory Physiology, vol. 305, no. 8, pp. H1265–H1273, 2013.

Y. Hu, V. Gurev, J. Constantino, J. D. Bayer, and N. A. Trayanova, “Effects of mechano-electric feedback on scroll wave stability in human ventricular fibrillation,” PLoS One, vol. 8, no. 4, Article ID e60287, 2013.

X. Jie, V. Gurev, and N. Trayanova, “Mechanisms of mechanically induced spontaneous arrhythmias in acute regional ischemia,” Circulation Research, vol. 106, no. 1, pp. 185–192, 2010.

K. M. Lim and D. U. Jeong, “Influence of the KCNQ1 S140G mutation on human ventricular arrhythmogenesis and pumping performance: simulation study,” Frontiers in Physiology, vol. 9, p. 926, 2018.

V. Gurev, T. Lee, J. Constantino, H. Arevalo, and N. A. Trayanova, “Models of cardiac electromechanics based on individual hearts imaging data,” Biomechanics and Modeling in Mechanobiology, vol. 10, no. 3, pp. 295–306, 2011.

J. Provost, V. Gurev, N. Trayanova, and E. E. Konofagou, “Mapping of cardiac electrical activation with electromechanical wave imaging: an in silico–in vivo reciprocity study,” Heart Rhythm, vol. 8, no. 5, pp. 752–759, 2011.

N. A. Trayanova, “Whole-heart modeling: applications to cardiac electrophysiology and electromechanics,” Circulation Research, vol. 108, no. 1, pp. 113–128, 2011.

N. A. Trayanova and J. J. Rice, “Cardiac electromechanical models: from cell to organ,” Frontiers in Physiology, vol. 2, p. 43, 2011.

K. M. Lim, J. Constantino, V. Gurev et al., “Comparison of the effects of continuous and pulsatile left ventricular-assist devices on ventricular unloading using a cardiac electromechanics model,” Journal of Physiological Sciences, vol. 62, no. 1, pp. 11–19, 2012.

K. ten Tusscher, D. Noble, P. Noble, and A. Panfilov, “A model for human ventricular tissue,” American Journal of Physiology-Heart and Circulatory Physiology, vol. 286, no. 4, pp. H1573–H1589, 2004.

O. Berenfeld and J. Jalife, “Purkinje-muscle reentry as a mechanism of polymorphic ventricular arrhythmias in a 3-dimensional model of the ventricles,” Circulation Research, vol. 82, no. 10, pp. 1063–1077, 1998.

J. J. Rice, F. Wang, D. M. Bers, and P. P. De Tombe, “Approximate model of cooperative activation and crossbridge cycling in cardiac muscle using ordinary differential equations,” Biophysical Journal, vol. 95, no. 5, pp. 2368–2390, 2008.

R. C. Kerckhoffs, M. L. Neal, Q. Gu et al., “Coupling of a 3D finite element model of cardiac ventricular mechanics to lumped systems models of the systemic and pulmonic circulation,” Annals of Biomedical Engineering, vol. 35, no. 1, pp. 1–18, 2007.

J. B. Shea and M. O. Sweeney, “Cardiac resynchronization therapy a patient’s guide,” Circulation, vol. 108, no. 9, pp. e64–e66, 2003.

Y. Hu, V. Gurev, J. Constantino, and N. Trayanova, “Optimizing cardiac resynchronization therapy to minimize ATP consumption heterogeneity throughout the left ventricle: a simulation analysis using a canine heart failure model,” Heart Rhythm, vol. 11, no. 6, pp. 1063–1069, 2014.