New insights from a computational model on the relation between pacing site and CRT response

How to assess and treat right ventricular electromechanical dyssynchrony in post-repair tetralogy of Fallot: insights from imaging, invasive studies, and computational modelling

BACKGROUND AND AIMS

Right bundle branch block (RBBB) and resulting right ventricular (RV) electromechanical discoordination are thought to play a role in the disease process of subpulmonary RV dysfunction that frequently occur post-repair tetralogy of Fallot (ToF). We sought to describe this disease entity, the role of pulmonary re-valvulation, and the potential added value of RV cardiac resynchronization therapy (RV-CRT).

METHODS

Two patients with repaired ToF, complete RBBB, pulmonary regurgitation, and significantly decreased RV function underwent echocardiography, cardiac magnetic resonance, and an invasive study to evaluate the potential for RV-CRT as part of the management strategy. The data were used to personalize the CircAdapt model of the human heart and circulation. Resulting Digital Twins were analysed to quantify the relative effects of RV pressure and volume overload and to predict the effect of RV-CRT.

RESULTS

Echocardiography showed components of a classic RV dyssynchrony pattern which could be reversed by RV-CRT during invasive study and resulted in acute improvement in RV systolic function. The Digital Twins confirmed a contribution of electromechanical RV dyssynchrony to RV dysfunction and suggested improvement of RV contraction efficiency after RV-CRT. The one patient who underwent successful permanent RV-CRT as part of the pulmonary re-valvulation procedure carried improvements that were in line with the predictions based on his Digital Twin.

CONCLUSION

An integrative diagnostic approach to RV dysfunction, including the construction of Digital Twins may help to identify candidates for RV-CRT as part of the lifetime management of ToF and similar congenital heart lesions.

25 years of basic and translational science in EP Europace: novel insights into arrhythmia mechanisms and therapeutic strategies

In the last 25 years, EP Europace has published more than 300 basic and translational science articles covering different arrhythmia types (ranging from atrial fibrillation to ventricular tachyarrhythmias), different diseases predisposing to arrhythmia formation (such as genetic arrhythmia disorders and heart failure), and different interventional and pharmacological anti-arrhythmic treatment strategies (ranging from pacing and defibrillation to different ablation approaches and novel drug-therapies). These studies have been conducted in cellular models, small and large animal models, and in the last couple of years increasingly in silico using computational approaches. In sum, these articles have contributed substantially to our pathophysiological understanding of arrhythmia mechanisms and treatment options; many of which have made their way into clinical applications. This review discusses a representative selection of EP Europace manuscripts covering the topics of pacing and ablation, atrial fibrillation, heart failure and pro-arrhythmic ventricular remodelling, ion channel (dys)function and pharmacology, inherited arrhythmia syndromes, and arrhythmogenic cardiomyopathies, highlighting some of the advances of the past 25 years. Given the increasingly recognized complexity and multidisciplinary nature of arrhythmogenesis and continued technological developments, basic and translational electrophysiological research is key advancing the field. EP Europace aims to further increase its contribution to the discovery of arrhythmia mechanisms and the implementation of mechanism-based precision therapy approaches in arrhythmia management.

Comparison of novel ventricular pacing strategies using an electro-mechanical simulation platform

AIMS

Focus of pacemaker therapy is shifting from right ventricular (RV) apex pacing (RVAP) and biventricular pacing (BiVP) to conduction system pacing. Direct comparison between the different pacing modalities and their consequences to cardiac pump function is difficult, due to the practical implications and confounding variables. Computational modelling and simulation provide the opportunity to compare electrical, mechanical, and haemodynamic consequences in the same virtual heart.

METHODS AND RESULTS

Using the same single cardiac geometry, electrical activation maps following the different pacing strategies were calculated using an Eikonal model on a three-dimensional geometry, which were then used as input for a lumped mechanical and haemodynamic model (CircAdapt). We then compared simulated strain, regional myocardial work, and haemodynamic function for each pacing strategy. Selective His-bundle pacing (HBP) best replicated physiological electrical activation and led to the most homogeneous mechanical behaviour. Selective left bundle branch (LBB) pacing led to good left ventricular (LV) function but significantly increased RV load. RV activation times were reduced in non-selective LBB pacing (nsLBBP), reducing RV load but increasing heterogeneity in LV contraction. LV septal pacing led to a slower LV and more heterogeneous LV activation than nsLBBP, while RV activation was similar. BiVP led to a synchronous LV-RV, but resulted in a heterogeneous contraction. RVAP led to the slowest and most heterogeneous contraction. Haemodynamic differences were small compared to differences in local wall behaviour.

CONCLUSION

Using a computational modelling framework, we investigated the mechanical and haemodynamic outcome of the prevailing pacing strategies in hearts with normal electrical and mechanical function. For this class of patients, nsLBBP was the best compromise between LV and RV function if HBP is not possible.

A new Electromechanical Wave Imaging dispersion metric for the characterization of ventricular activation in different Cardiac Resynchronization Therapy pacing schemes

Conventional biventricular (BiV) pacing cardiac resynchronization therapy (CRT) is an established treatment for heart failure patients. Recently, multiple novel CRT delivering technologies such as His-Bundle pacing have been investigated as alternative pacing strategies for optimal treatment benefit. Electromechanical Wave Imaging (EWI), a high frame-rate echocardiography-based modality, is capable of visualizing the change from dyssynchronous activation to resynchronized BiV-paced ventricles in 3D. This proof-of-concept study introduces a new EWI-based dispersion metric to further characterize ventricular activation. Patients with His-Bundle device implantation (n=4), left-bundle branch block (n=10), right-ventricular (RV) pacing (n=10), or BiV pacing (n=15) were imaged, as well as four volunteers in normal sinus rhythm (NSR). EWI successfully mapped the ventricular activation resulting from His-Bundle pacing. Additionally, very similar activation patterns were obtained in the NSR subjects, confirming recovery of physiological activation with His pacing. The dispersion metric was the most sensitive EWI-based metric that identified His pacing as the most efficient treatment (lowest activation time spread), followed by BiV and RV pacing. More specifically, the dispersion metric significantly (p 0.005) distinguished His pacing from the other two pacing schemes as well as LBBB. The initial findings presented herein indicate that EWI and its new dispersion metric may provide a useful resynchronization evaluation clinical tool in CRT patients under both novel His-Bundle pacing and more conventional BiV pacing strategies.

Acute Hemodynamic Effects of Simultaneous and Sequential Multi-Point Pacing in Heart Failure Patients With an Expected Higher Rate of Sub-response to Cardiac Resynchronization Therapy: Results of Multicenter SYNSEQ Study

The aim of the SYNSEQ (Left Ventricular Synchronous vs. Sequential MultiSpot Pacing for CRT) study was to evaluate the acute hemodynamic response (AHR) of simultaneous (3P-MPP _syn) or sequential (3P-MPP _seq) multi-3-point-left-ventricular (LV) pacing vs. single point pacing (SPP) in a group of patients at risk of a suboptimal response to cardiac resynchronization therapy (CRT). Twenty five patients with myocardial scar or QRS ≤ 150 or the absence of LBBB (age: 66 ± 12 years, QRS: 159 ± 12 ms, NYHA class II/III, LVEF ≤ 35%) underwent acute hemodynamic assessment by LV + dP/dt_max with a variety of LV pacing configurations at an optimized AV delay. The change in LV + dP/dt _max (%ΔLV + dP/dt _max) with 3P-MPP _syn (15.6%, 95% CI: 8.8%-22.5%) was neither statistically significantly different to 3P-MPP _seq (11.8%, 95% CI: 7.6-16.0%) nor to SPP _basal (11.5%, 95% CI:7.1-15.9%) or SPP _mid (12.2%, 95% CI:7.9-16.5%), but higher than SPP _apical (10.6%, 95% CI:5.3-15.9%, p = 0.03). AHR (defined as a %ΔLV + dP/dt _max ≥ 10%) varied between pacing configurations: 36% (9/25) for SPP _apical, 44% (11/25) for SPP _basal, 54% (13/24) for SPP _mid, 56% (14/25) for 3P-MPP _syn and 48% (11/23) for 3P-MPP _seq.Fifteen patients (15/25, 60%) had an AHR in at least one pacing configuration. AHR was observed in 10/13 (77%) patients with a LBBB but only in 5/12 (42%) patients with a non-LBBB (p = 0.11). To conclude, simultaneous or sequential multipoint pacing compared to single point pacing did not improve the acute hemodynamic effect in a suboptimal CRT response population.

An accurate, robust, and efficient finite element framework with applications to anisotropic, nearly and fully incompressible elasticity

Fiber-reinforced soft biological tissues are typically modeled as hyperelastic, anisotropic, and nearly incompressible materials. To enforce incompressibility a multiplicative split of the deformation gradient into a volumetric and an isochoric part is a very common approach. However, the finite element analysis of such problems often suffers from severe volumetric locking effects and numerical instabilities. In this paper, we present novel methods to overcome volumetric locking phenomena for using stabilized P1-P1 elements. We introduce different stabilization techniques and demonstrate the high robustness and computational efficiency of the chosen methods. In two benchmark problems from the literature as well as an advanced application to cardiac electromechanics, we compare the approach to standard linear elements and show the accuracy and versatility of the methods to simulate anisotropic, nearly and fully incompressible materials. We demonstrate the potential of this numerical framework to accelerate accurate simulations of biological tissues to the extent of enabling patient-specific parameterization studies, where numerous forward simulations are required.

Optimization of cardiac resynchronization therapy based on a cardiac electromechanics-perfusion computational model

Cardiac resynchronization therapy (CRT) is an established treatment for left bundle branch block (LBBB) resulting in mechanical dyssynchrony. Approximately 1/3 of patients with CRT, however, are non-responders. To understand factors affecting CRT response, an electromechanics-perfusion computational model based on animal-specific left ventricular (LV) geometry and coronary vascular networks located in the septum and LV free wall is developed. The model considers contractility-flow and preload-activation time relationships, and is calibrated to simultaneously match the experimental measurements in terms of the LV pressure, volume waveforms and total coronary flow in the left anterior descending and left circumflex territories from 2 swine models under right atrium and right ventricular pacing. The model is then applied to investigate the responses of CRT indexed by peak LV pressure and (dP/dt)max at multiple pacing sites with different degrees of perfusion in the LV free wall. Without the presence of ischemia, the model predicts that basal-lateral endocardial region is the optimal pacing site that can best improve (dP/dt)max by 20%, and is associated with the shortest activation time. In the presence of ischemia, a non-ischemic region becomes the optimal pacing site when coronary flow in the ischemic region fell below 30% of its original value. Pacing at the ischemic region produces little response at that perfusion level. The optimal pacing site is associated with one that optimizes the LV activation time. These findings suggest that CRT response is affected by both pacing site and coronary perfusion, which may have clinical implication in improving CRT responder rates.

Computational models in cardiology

The treatment of individual patients in cardiology practice increasingly relies on advanced imaging, genetic screening and devices. As the amount of imaging and other diagnostic data increases, paralleled by the greater capacity to personalize treatment, the difficulty of using the full array of measurements of a patient to determine an optimal treatment seems also to be paradoxically increasing. Computational models are progressively addressing this issue by providing a common framework for integrating multiple data sets from individual patients. These models, which are based on physiology and physics rather than on population statistics, enable computational simulations to reveal diagnostic information that would have otherwise remained concealed and to predict treatment outcomes for individual patients. The inherent need for patient-specific models in cardiology is clear and is driving the rapid development of tools and techniques for creating personalized methods to guide pharmaceutical therapy, deployment of devices and surgical interventions.

Optimization of Lead Placement in the Right Ventricle During Cardiac Resynchronization Therapy. A Simulation Study

Patients suffering from heart failure and left bundle branch block show electrical ventricular dyssynchrony causing an abnormal blood pumping. Cardiac resynchronization therapy (CRT) is recommended for these patients. Patients with positive therapy response normally present QRS shortening and an increased left ventricle (LV) ejection fraction. However, around one third do not respond favorably. Therefore, optimal location of pacing leads, timing delays between leads and/or choosing related biomarkers is crucial to achieve the best possible degree of ventricular synchrony during CRT application. In this study, computational modeling is used to predict the optimal location and delay of pacing leads to improve CRT response. We use a 3D electrophysiological computational model of the heart and torso to get insight into the changes in the activation patterns obtained when the heart is paced from different regions and for different atrioventricular and interventricular delays. The model represents a heart with left bundle branch block and heart failure, and allows a detailed and accurate analysis of the electrical changes observed simultaneously in the myocardium and in the QRS complex computed in the precordial leads. Computational simulations were performed using a modified version of the O'Hara et al. action potential model, the most recent mathematical model developed for human ventricular electrophysiology. The optimal location for the pacing leads was determined by QRS maximal reduction. Additionally, the influence of Purkinje system on CRT response was assessed and correlation analysis between several parameters of the QRS was made. Simulation results showed that the right ventricle (RV) upper septum near the outflow tract is an alternative location to the RV apical lead. Furthermore, LV endocardial pacing provided better results as compared to epicardial stimulation. Finally, the time to reach the 90% of the QRS area was a good predictor of the instant at which 90% of the ventricular tissue was activated. Thus, the time to reach the 90% of the QRS area is suggested as an additional index to assess CRT effectiveness to improve biventricular synchrony.

Long-Term Outcomes of Cardiac Resynchronization Therapy Using Apical Versus Nonapical Left Ventricular Pacing

Background Experimental evidence indicates that left ventricular ( LV ) apical pacing is hemodynamically superior to nonapical LV pacing. Some studies have shown that an LV apical lead position is unfavorable in cardiac resynchronization therapy. We sought to determine whether an apical LV lead position influences cardiac mortality after cardiac resynchronization therapy. Methods and Results In this retrospective observational study, the primary end point of cardiac mortality was assessed in relation to longitudinal (basal, midventricular, or apical) and circumferential (anterior, lateral, or posterior) LV lead positions, as well as right ventricular (apical or septal), assigned using fluoroscopy. Lead positions were assessed in 1189 patients undergoing cardiac resynchronization therapy implantation over 15 years. After a median follow-up of 6.0 years (interquartile range: 4.4-7.7 years), an apical LV lead position was associated with lower cardiac mortality than a nonapical position (adjusted hazard ratio: 0.74; 95% confidence interval, 0.56-0.99) after covariate adjustment. There were no differences in total mortality or heart failure hospitalization. Death from pump failure was lower with apical than nonapical positions (adjusted hazard ratio: 0.69; 95% confidence interval, 0.51-0.94). Compared with a basal position, an apical LV position was also associated with lower risk of sudden cardiac death (adjusted hazard ratio: 0.34; 95% confidence interval, 0.13-0.93). No differences emerged between circumferential LV lead positions or right ventricular positions with respect to any end point. Conclusions In recipients of cardiac resynchronization therapy, an apical LV lead position was associated with better long-term cardiac survival than a nonapical position. This effect was due to a lower risk of pump failure and sudden cardiac death.

Computational Modeling for Cardiac Resynchronization Therapy

Cardiac resynchronization therapy (CRT) is an effective treatment for heart failure (HF) patients with an electrical substrate pathology causing ventricular dyssynchrony. However 40-50% of patients do not respond to treatment. Cardiac modeling of the electrophysiology, electromechanics, and hemodynamics of the heart has been used to study mechanisms behind HF pathology and CRT response. Recently, multi-scale dyssynchronous HF models have been used to study optimal device settings and optimal lead locations, investigate the underlying cardiac pathophysiology, as well as investigate emerging technologies proposed to treat cardiac dyssynchrony. However the breadth of patient and experimental data required to create and parameterize these models and the computational resources required currently limits the use of these models to small patient numbers. In the future, once these technical challenges are overcome, biophysically based models of the heart have the potential to become a clinical tool to aid in the diagnosis and treatment of HF.

Clinical Applications of Patient-Specific Models: The Case for a Simple Approach

Over the past several decades, increasingly sophisticated models of the heart have provided important insights into cardiac physiology and are increasingly used to predict the impact of diseases and therapies on the heart. In an era of personalized medicine, many envision patient-specific computational models as a powerful tool for personalizing therapy. Yet the complexity of current models poses important challenges, including identifying model parameters and completing calculations quickly enough for routine clinical use. We propose that early clinical successes are likely to arise from an alternative approach: relatively simple, fast, phenomenologic models with a small number of parameters that can be easily (and automatically) customized. We discuss examples of simple yet foundational models that have already made a tremendous impact on clinical education and practice, and make the case that reducing rather than increasing model complexity may be the key to realizing the promise of patient-specific modeling for clinical applications.