Multisite pacing and myocardial scars: a computational study

Regional Left Ventricular Fiber Stress Analysis for Cardiac Resynchronization Therapy Response

Cardiac resynchronization therapy (CRT) is an effective treatment for a subgroup of heart failure (HF) patients, but more than 30% of those selected do not improve after CRT implantation. Imperfect pre-procedural criteria for patient selection and optimization are the main causes of the high non-response rate. In this study, we evaluated a novel measure for assessing CRT response. We used a computational modeling framework to calculate the regional stress of the left ventricular wall of seven CRT patients and seven healthy controls. The standard deviation of regional wall stress at the time of mitral valve closure (SD_MVC) was used to quantify dyssynchrony and compared between patients and controls and among the patients. The results show that SD_MVC is significantly lower in controls than patients and correlates with long-term response in patients, based on end-diastolic volume reduction. In contrast to our initial hypothesis, patients with lower SD_MVC respond better to therapy. The patient with the highest SD_MVC was the only non-responder in the patient cohort. The distribution of fiber stress at the beginning of the isovolumetric phase seems to correlate with the degree of response and the use of this measurement could potentially improve selection criteria for CRT implantation. Further studies with a larger cohort of patients are needed to validate these results.

Determinants of the time-to-peak left ventricular dP/dt (Td) and QRS duration with different fusion strategies in cardiac resynchronization therapy

Background

Cardiac resynchronization therapy (CRT) is helpful in selected patients; however, responder rates rarely exceed 70%. Optimization of CRT may therefore benefit a large number of patients. Time-to-peak dP/dt (Td) is a novel marker of myocardial synergy that reflects the degree of myocardial dyssynchrony with the potential to guide and optimize treatment with CRT. Optimal electrical activation is a prerequisite for CRT to be effective. Electrical activation can be altered by changing the electrical wave-front fusion resulting from pacing to optimize resynchronization. We designed this study to understand the acute effects of different electrical wave-front fusion strategies and LV pre-/postexcitation on Td and QRS duration (QRSd). A better understanding of measuring and optimizing resynchronization can help improve the benefits of CRT.

Methods

Td and QRSd were measured in 19 patients undergoing a CRT implantation. Two biventricular pacing groups were compared: pacing the left ventricle (LV) with fusion with intrinsic right ventricular activation (FUSION group) and pacing the LV and right ventricle (RV) at short atrioventricular delay (STANDARD group) to avoid fusion with intrinsic RV activation. A quadripolar LV lead enabled pacing from widely separated electrodes; distal (DIST), proximal (PROX) and both electrodes combined (multipoint pacing, MPP). The LV was stimulated relative in time to RV activation (either RV pace-onset or QRS-onset), with the LV stimulated prior to (PRE), simultaneous with (SIM) or after (POST) RV activation. In addition, we analyzed the interactions of the two groups (FUSION/STANDARD) with three different electrode configurations (DIST, PROX, MPP), each paced with three different degrees of LV pre-/postexcitation (PRE, SIM, POST) in a statistical model.

Results

We found that FUSION provided shorter Td and QRSd than STANDARD, MPP provided shorter Td and QRSd than DIST and PROX, and SIM provided both the shortest QRSd and Td compared to PRE and POST. The interaction analysis revealed that pacing MPP with fusion with intrinsic RV activation simultaneous with the onset of the QRS complex (MPP^*FUSION^*SIM) shortened QRSd and Td the most compared to all other modes and configurations. The difference in QRSd and Td from their respective references were significantly correlated (β = 1, R = 0.9, p < 0.01).

Conclusion

Pacing modes and electrode configurations designed to optimize electrical wave-front fusion (intrinsic RV activation, LV multipoint pacing and simultaneous RV and LV activation) shorten QRSd and Td the most. As demonstrated in this study, electrical and mechanical measures of resynchronization are highly correlated. Therefore, Td can potentially serve as a marker for CRT optimization.

Computational electrophysiology of the coronary sinus branches based on electro-anatomical mapping for the prediction of the latest activated region

This work dealt with the assessment of a computational tool to estimate the electrical activation in the left ventricle focusing on the latest electrically activated segment (LEAS) in patients with left bundle branch block and possible myocardial fibrosis. We considered the Eikonal-diffusion equation and to recover the electrical activation maps in the myocardium. The model was calibrated by using activation times acquired in the coronary sinus (CS) branches or in the CS solely with an electroanatomic mapping system (EAMS) during cardiac resynchronization therapy (CRT). We applied our computational tool to ten patients founding an excellent accordance with EAMS measures; in particular, the error for LEAS location was less than 4 mm. We also calibrated our model using only information in the CS, still obtaining an excellent agreement with the measured LEAS. The proposed tool was able to accurately reproduce the electrical activation maps and in particular LEAS location in the CS branches, with an almost real-time computational effort, regardless of the presence of myocardial fibrosis, even when information only at CS was used to calibrate the model. This could be useful in the clinical practice since LEAS is often used as a target site for the left lead placement during CRT.

Meshless Electrophysiological Modeling of Cardiac Resynchronization Therapy—Benchmark Analysis with Finite-Element Methods in Experimental Data

Computational models of cardiac electrophysiology are promising tools for reducing the rates of non-response patients suitable for cardiac resynchronization therapy (CRT) by optimizing electrode placement. The majority of computational models in the literature are mesh-based, primarily using the finite element method (FEM). The generation of patient-specific cardiac meshes has traditionally been a tedious task requiring manual intervention and hindering the modeling of a large number of cases. Meshless models can be a valid alternative due to their mesh quality independence. The organization of challenges such as the CRT-EPiggy19, providing unique experimental data as open access, enables benchmarking analysis of different cardiac computational modeling solutions with quantitative metrics. We present a benchmark analysis of a meshless-based method with finite-element methods for the prediction of cardiac electrical patterns in CRT, based on a subset of the CRT-EPiggy19 dataset. A data assimilation strategy was designed to personalize the most relevant parameters of the electrophysiological simulations and identify the optimal CRT lead configuration. The simulation results obtained with the meshless model were equivalent to FEM, with the most relevant aspect for accurate CRT predictions being the parameter personalization strategy (e.g., regional conduction velocity distribution, including the Purkinje system and CRT lead distribution).

A clinical-in silico study on the effectiveness of multipoint bicathodic and cathodic-anodal pacing in cardiac resynchronization therapy

Up to one-third of patients undergoing cardiac resynchronization therapy (CRT) are nonresponders. Multipoint bicathodic and cathodic-anodal left ventricle (LV) stimulations could overcome this clinical challenge, but their effectiveness remains controversial. Here we evaluate the performance of such stimulations through both in vivo and in silico experiments, the latter based on computer electromechanical modeling. Seven patients, all candidates for CRT, received a quadripolar LV lead. Four stimulations were tested: right ventricular (RVS); conventional single point biventricular (S-BS); multipoint biventricular bicathodic (CC-BS) and multipoint biventricular cathodic-anodal (CA-BS). The following parameters were processed: QRS duration; maximal time derivative of arterial pressure (dPdt_max); systolic arterial pressure (P_sys); and stroke volume (SV). Echocardiographic data of each patient were then obtained to create an LV geometric model. Numerical simulations were based on a strongly coupled Bidomain electromechanical coupling model. Considering the in vivo parameters, when comparing S-BS to RVS, there was no significant decrease in SV (from 45 ± 11 to 44 ± 20 ml) and 6% and 4% increases of dPdt_max and P_sys, respectively. Focusing on in silico parameters, with respect to RVS, S-BS exhibited a significant increase of SV, dPdt_max and P_sys. Neither the in vivo nor in silico results showed any significant hemodynamic and electrical difference among S-BS, CC-BS and CA-BS configurations. These results show that CC-BS and CA-BS yield a comparable CRT performance, but they do not always yield improvement in terms of hemodynamic parameters with respect to S-BS. The computational results confirmed the in vivo observations, thus providing theoretical support to the clinical experiments.

Determinants of LV dP/dtmax and QRS duration with different fusion strategies in cardiac resynchronisation therapy

BACKGROUND

We designed this study to assess the acute effects of different fusion strategies and left ventricular (LV) pre-excitation/post-excitation on LV dP/dtmax and QRS duration (QRSd).

METHODS

We measured LV dP/dtmax and QRSd in 19 patients having cardiac resynchronisation therapy (CRT). Two groups of biventricular pacing were compared: pacing the left ventricle (LV) with FUSION with intrinsic right ventricle (RV) activation (FUSION), and pacing the LV and RV with NO FUSION with intrinsic RV activation. In the NO FUSION group, the RV was paced before the expected QRS onset. A quadripolar LV lead enabled distal, proximal and multipoint pacing (MPP). The LV was stimulated relative in time to either RV pace or QRS-onset in four pre-excitation/post-excitation classes (PCs). We analysed the interactions of two groups (FUSION/NO FUSION) with three different electrode configurations, each paced with four different degrees of LV pre-excitation (PC1-4) in a statistical model.

RESULTS

LV dP/dtmax was higher with NO FUSION than with FUSION (769±46 mm Hg/s vs 746±46 mm Hg/s, p<0.01), while there was no difference in QRSd (NO FUSION 156±2 ms and FUSION 155±2 ms). LV dP/dtmax and QRSd increased with LV pre-excitation compared with pacing timed to QRS/RV pace-onset regardless of electrode configuration. Overall, pacing LV close to QRS-onset (FUSION) with MPP shortened QRSd the most, while LV dP/dtmax increased the most with LV pre-excitation.

CONCLUSION

We show how a beneficial change in QRSd dissociates from the haemodynamic change in LV dP/dtmax with different biventricular pacing strategies. In this study, LV pre-excitation was the main determinant of LV dP/dtmax, while QRSd shortens with optimal resynchronisation.