Efficient preloading of the ventricles by a properly timed atrial contraction underlies stroke work improvement in the acute response to cardiac resynchronization therapy

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.

Within-patient comparison of His-bundle pacing, right ventricular pacing, and right ventricular pacing avoidance algorithms in patients with PR prolongation: Acute hemodynamic study

AIMS

A prolonged PR interval may adversely affect ventricular filling and, therefore, cardiac function. AV delay can be corrected using right ventricular pacing (RVP), but this induces ventricular dyssynchrony, itself harmful. Therefore, in intermittent heart block, pacing avoidance algorithms are often implemented. We tested His-bundle pacing (HBP) as an alternative.

METHODS

Outpatients with a long PR interval (>200 ms) and intermittent need for ventricular pacing were recruited. We measured within-patient differences in high-precision hemodynamics between AV-optimized RVP and HBP, as well as a pacing avoidance algorithm (Managed Ventricular Pacing [MVP]).

RESULTS

We recruited 18 patients. Mean left ventricular ejection fraction was 44.3 ± 9%. Mean intrinsic PR interval was 266 ± 42 ms and QRS duration was 123 ± 29 ms. RVP lengthened QRS duration (+54 ms, 95% CI 42-67 ms, p < .0001) while HBP delivered a shorter QRS duration than RVP (-56 ms, 95% CI -67 to -46 ms, p < .0001). HBP did not increase QRS duration (-2 ms, 95% CI -8 to 13 ms, p = .6). HBP improved acute systolic blood pressure by mean of 5.0 mmHg (95% CI 2.8-7.1 mmHg, p < .0001) compared to RVP and by 3.5 mmHg (95% CI 1.9-5.0 mmHg, p = .0002) compared to the pacing avoidance algorithm. There was no significant difference in hemodynamics between RVP and ventricular pacing avoidance (p = .055).

CONCLUSIONS

HBP provides better acute cardiac function than pacing avoidance algorithms and RVP, in patients with prolonged PR intervals. HBP allows normalization of prolonged AV delays (unlike pacing avoidance) and does not cause ventricular dyssynchrony (unlike RVP). Clinical trials may be justified to assess whether these acute improvements translate into longer term clinical benefits in patients with bradycardia indications for pacing.

Left ventricular response after cardiac resynchronization therapy is related to early left atrial volume reduction

BACKGROUND/AIMS

The current study aimed to elucidate a time-course change in left atrial volume after cardiac resynchronization therapy (CRT) and to verify factors associated with left atrial volume reduction (LAVR) and its prognostic implications.

METHODS

The records of 97 patients were retrospectively reviewed after CRT. Echocardiographic data were analyzed at baseline before CRT, at early follow-up (FU) (≤ 1 year, median 6 months), and at late FU (median 30 months). Left ventricular volume response (LVVR) was defined as 15% reduction in left ventricular (LV) end-systolic volume (ESV). LAVR was classified into two groups by the median value at early FU: LAVR (≥ 7.5%) and no LAVR (< 7.5%).

RESULTS

LV ESV index continuously decreased from baseline to early FU and from early FU to late FU (106.1 ± 47.4 mL/m2 vs. 87.6 ± 51.6 mL/m2 vs. 72.5 ± 57.1 mL/m2). LA volume index decreased from baseline to early FU, but there were no reductions thereafter (51.8 ± 21.9 mL/m2 vs. 45.1 ± 19.6 mL/m2 vs. 44.9 ± 23.0 mL/m2). The only echocardiographic factor associated with LAVR was change in E velocity (odds ratio [OR], 1.04; p = 0.002). Early LAVR (OR, 10.05; p = 0.002) was an independent predictor for late LVVR.

CONCLUSION

LAVR was related to reduction in E velocity, suggesting its relation with optimization of LV filling pressure. Early LAVR was a predictor for LVVR to CRT in long-term FU.

Silencing miR-370-3p rescues funny current and sinus node function in heart failure

Bradyarrhythmias are an important cause of mortality in heart failure and previous studies indicate a mechanistic role for electrical remodelling of the key pacemaking ion channel HCN4 in this process. Here we show that, in a mouse model of heart failure in which there is sinus bradycardia, there is upregulation of a microRNA (miR-370-3p), downregulation of the pacemaker ion channel, HCN4, and downregulation of the corresponding ionic current, If, in the sinus node. In vitro, exogenous miR-370-3p inhibits HCN4 mRNA and causes downregulation of HCN4 protein, downregulation of If, and bradycardia in the isolated sinus node. In vivo, intraperitoneal injection of an antimiR to miR-370-3p into heart failure mice silences miR-370-3p and restores HCN4 mRNA and protein and If in the sinus node and blunts the sinus bradycardia. In addition, it partially restores ventricular function and reduces mortality. This represents a novel approach to heart failure treatment.

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.

4D cardiac electromechanical activation imaging

Cardiac abnormalities, a major cause of morbidity and mortality, affect millions of people worldwide. Despite the urgent clinical need for early diagnosis, there is currently no noninvasive technique that can infer to the electrical function of the whole heart in 3D and thereby localize abnormalities at the point of care. Here we present a new method for noninvasive 4D mapping of the cardiac electromechanical activity in a single heartbeat for heart disease characterization such as arrhythmia and infarction. Our novel technique captures the 3D activation wave of the heart in vivo using high volume-rate (500 volumes per second) ultrasound with a 32 × 32 matrix array. Electromechanical activation maps are first presented in a normal and infarcted cardiac model in silico and in canine heart during pacing and re-entrant ventricular tachycardia in vivo. Noninvasive 4D electromechanical activation mapping in a healthy volunteer and a heart failure patient are also determined. The technique described herein allows for direct, simultaneous and noninvasive visualization of electromechanical activation in 3D, which provides complementary information on myocardial viability and/or abnormality to clinical imaging.

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

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.

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.

Atrioventricular Conduction Delay Predicts Impaired Exercise Capacity in Patients with Heart Failure with Reduced Ejection Fraction

Background

Atrioventricular conduction delay (AVCD) impairs left ventricular (LV) filling and consequently leads to a reduction of cardiac output. We hypothesized that in patients with severely depressed LV function and coexisting intraventricular conduction disturbances (IVCD), AVCD can affect exercise performance. Therefore, we evaluated the association of AVCD and exercise capacity in patients with heart failure (HFREF) and coexisting IVCD.

Material/Methods

We included patients with stable, chronic HFREF, LVEF <35%, sinus rhythm, and QRS ≥120 ms. PR interval and peak oxygen consumption (VO_{2 peak}) were specifically investigated. Multiple regression analysis was used to adjust the association between PR interval and VO_{2 peak} for possible confounders.

Results

Most (57.5%) of the 40 included patients [20% female, aged 63±12, 47.5% of ischemic etiology (IHD)] were in NYHA class III. Mean PR interval was 196±38.1 ms. There were 26 (65%) patients with PR interval ≤200 ms and 14 (35%) with >200 ms. Groups were similar in clinical, laboratory, echocardiographic parameters, QRS morphology, and treatment regimens. VO_{2 peak} was lower in patients with longer PR interval group as compared to shorter PR interval group (12.3±4.1 vs. 17.06±4.4, p=0.002). In the regression model, PR interval, female sex, and IHD remained important predictors of VO_{2 peak} (_partial=−0.50, p=0.003; r_partial=−0.48, p=0.005; r_partial=−0.44, p=0.01; R²=0.61).

Conclusions

Delayed AV conduction contributes to decreased exercise capacity in patients with HFREF and coexisting IVCD.

Effect of Cardiac Resynchronization Therapy on Left Atrial Size and Function as Expressed by Speckle Tracking 2-Dimensional Strain

Changes in left atrial (LA) strain in patients treated with cardiac resynchronization therapy (CRT) remain not entirely explored. We prospectively evaluated long-term changes in LA size and function and their relation with left ventricular (LV) reverse remodeling and noninvasive hemodynamic variables in patients treated with CRT by 2-dimensional speckle tracking echocardiography. Thirty patients (62 ± 11 years, 63% men) underwent 2-dimensional speckle tracking echocardiography before implant and after 12 months. LA area, global and regional LA strains, LV ejection fraction (LVEF) and longitudinal strain, mitral regurgitation (MR), and diastolic variables were evaluated. At 12 months, CRT responders (60%) exhibited an increase in LA strain (11.4 ± 6.5% vs 16.5 ± 7.9%, p <0.001) and a reduction in LA area (p = 0.002), which were associated with an improvement in MR, E/E' ratio, LVEF, and LV longitudinal strain. In nonresponders, a worsening in LA strain (11.4 ± 6.8% vs 8.7 ± 4.6%, p = 0.017) and LA area (p = 0.002) occurred in parallel with an increase in E/E', whereas LVEF and LV longitudinal strain were unchanged. In conclusion, over long-term follow-up, LA size and strain improved in CRT responders, while worsening in nonresponders. Changes in LV function, filling pressures, and MR seem to be related to LA reverse remodeling, giving a feedback loop.

How computer simulations of the human heart can improve anti-arrhythmia therapy

Over the last decade, the state-of-the-art in cardiac computational modelling has progressed rapidly. The electrophysiological function of the heart can now be simulated with a high degree of detail and accuracy, opening the doors for simulation-guided approaches to anti-arrhythmic drug development and patient-specific therapeutic interventions. In this review, we outline the basic methodology for cardiac modelling, which has been developed and validated over decades of research. In addition, we present several recent examples of how computational models of the human heart have been used to address current clinical problems in cardiac electrophysiology. We will explore the use of simulations to improve anti-arrhythmic pacing and defibrillation interventions; to predict optimal sites for clinical ablation procedures; and to aid in the understanding and selection of arrhythmia risk markers. Together, these studies illustrate how the tremendous advances in cardiac modelling are poised to revolutionize medical treatment and prevention of arrhythmia.

A New MRI-Based Model of Heart Function with Coupled Hemodynamics and Application to Normal and Diseased Canine Left Ventricles

A methodology for the simulation of heart function that combines an MRI-based model of cardiac electromechanics (CE) with a Navier–Stokes-based hemodynamics model is presented. The CE model consists of two coupled components that simulate the electrical and the mechanical functions of the heart. Accurate representations of ventricular geometry and fiber orientations are constructed from the structural magnetic resonance and the diffusion tensor MR images, respectively. The deformation of the ventricle obtained from the electromechanical model serves as input to the hemodynamics model in this one-way coupled approach via imposed kinematic wall velocity boundary conditions and at the same time, governs the blood flow into and out of the ventricular volume. The time-dependent endocardial surfaces are registered using a diffeomorphic mapping algorithm, while the intraventricular blood flow patterns are simulated using a sharp-interface immersed boundary method-based flow solver. The utility of the combined heart-function model is demonstrated by comparing the hemodynamic characteristics of a normal canine heart beating in sinus rhythm against that of the dyssynchronously beating failing heart. We also discuss the potential of coupled CE and hemodynamics models for various clinical applications.

Patient-specific modelling of cardiac electrophysiology in heart-failure patients

Aims

Left-ventricular (LV) conduction disturbances are common in heart-failure patients and a left bundle-branch block (LBBB) electrocardiogram (ECG) type is often seen. The precise cause of this pattern is uncertain and is probably variable between patients, ranging from proximal interruption of the left bundle branch to diffuse distal conduction disease in the working myocardium. Using realistic numerical simulation methods and patient-tailored model anatomies, we investigated different hypotheses to explain the observed activation order on the LV endocardium, electrogram morphologies, and ECG features in two patients with heart failure and LBBB ECG.

Methods and results

Ventricular electrical activity was simulated using reaction–diffusion models with patient-specific anatomies. From the simulated action potentials, ECGs and cardiac electrograms were computed by solving the bidomain equation. Model parameters such as earliest activation sites, tissue conductivity, and densities of ionic currents were tuned to reproduce the measured signals. Electrocardiogram morphology and activation order could be matched simultaneously. Local electrograms matched well at some sites, but overall the measured waveforms had deeper S-waves than the simulated waveforms.

Conclusion

Tuning a reaction–diffusion model of the human heart to reproduce measured ECGs and electrograms is feasible and may provide insights in individual disease characteristics that cannot be obtained by other means.

Computational modeling of cardiac optogenetics: Methodology overview & review of findings from simulations

Cardiac optogenetics is emerging as an exciting new potential avenue to enable spatiotemporally precise control of excitable cells and tissue in the heart with low-energy optical stimuli. This approach involves the expression of exogenous light-sensitive proteins (opsins) in target heart tissue via viral gene or cell delivery. Preliminary experiments in optogenetically-modified cells, tissue, and organisms have made great strides towards demonstrating the feasibility of basic applications, including the use of light stimuli to pace or disrupt reentrant activity. However, it remains unknown whether techniques based on this intriguing technology could be scaled up and used in humans for novel clinical applications, such as pain-free optical defibrillation or dynamic modulation of action potential shape. A key step towards answering such questions is to explore potential optogenetics-based therapies using sophisticated computer simulation tools capable of realistically representing opsin delivery and light stimulation in biophysically detailed, patient-specific models of the human heart. This review provides (1) a detailed overview of the methodological developments necessary to represent optogenetics-based solutions in existing virtual heart platforms and (2) a survey of findings that have been derived from such simulations and a critical assessment of their significance with respect to the progress of the field.

Computational analysis of the effect of valvular regurgitation on ventricular mechanics using a 3D electromechanics model

Using a three-dimensional electromechanical model of the canine ventricles with dyssynchronous heart failure, we investigated the relationship between severity of valve regurgitation and ventricular mechanical responses. The results demonstrated that end-systolic tension in the septum and left ventricular free wall was significantly lower under the condition of mitral regurgitation (MR) than under aortic regurgitation (AR). Stroke work in AR was higher than that in MR. On the other hand, the difference in stroke volume between the two conditions was not significant, indicating that AR may cause worse pumping efficiency than MR in terms of consumed energy and performed work.

Methodology for image-based reconstruction of ventricular geometry for patient-specific modeling of cardiac electrophysiology

Patient-specific modeling of ventricular electrophysiology requires an interpolated reconstruction of the 3-dimensional (3D) geometry of the patient ventricles from the low-resolution (Lo-res) clinical images. The goal of this study was to implement a processing pipeline for obtaining the interpolated reconstruction, and thoroughly evaluate the efficacy of this pipeline in comparison with alternative methods. The pipeline implemented here involves contouring the epi- and endocardial boundaries in Lo-res images, interpolating the contours using the variational implicit functions method, and merging the interpolation results to obtain the ventricular reconstruction. Five alternative interpolation methods, namely linear, cubic spline, spherical harmonics, cylindrical harmonics, and shape-based interpolation were implemented for comparison. In the thorough evaluation of the processing pipeline, Hi-res magnetic resonance (MR), computed tomography (CT), and diffusion tensor (DT) MR images from numerous hearts were used. Reconstructions obtained from the Hi-res images were compared with the reconstructions computed by each of the interpolation methods from a sparse sample of the Hi-res contours, which mimicked Lo-res clinical images. Qualitative and quantitative comparison of these ventricular geometry reconstructions showed that the variational implicit functions approach performed better than others. Additionally, the outcomes of electrophysiological simulations (sinus rhythm activation maps and pseudo-ECGs) conducted using models based on the various reconstructions were compared. These electrophysiological simulations demonstrated that our implementation of the variational implicit functions-based method had the best accuracy.

Images as drivers of progress in cardiac computational modelling

Computational models have become a fundamental tool in cardiac research. Models are evolving to cover multiple scales and physical mechanisms. They are moving towards mechanistic descriptions of personalised structure and function, including effects of natural variability. These developments are underpinned to a large extent by advances in imaging technologies. This article reviews how novel imaging technologies, or the innovative use and extension of established ones, integrate with computational models and drive novel insights into cardiac biophysics. In terms of structural characterization, we discuss how imaging is allowing a wide range of scales to be considered, from cellular levels to whole organs. We analyse how the evolution from structural to functional imaging is opening new avenues for computational models, and in this respect we review methods for measurement of electrical activity, mechanics and flow. Finally, we consider ways in which combined imaging and modelling research is likely to continue advancing cardiac research, and identify some of the main challenges that remain to be solved.

Optimizing cardiac resynchronization therapy to minimize ATP consumption heterogeneity throughout the left ventricle: a simulation analysis using a canine heart failure model

BACKGROUND

Cardiac resynchronization therapy (CRT) has been demonstrated to lead to restoration of oxygen consumption homogeneity throughout the left ventricle (LV), which is important for long-term reverse remodeling of the ventricles. However, research has focused exclusively on identifying the LV pacing sites that led to acute hemodynamic improvements. It remains unclear whether there exist LV pacing sites that could both improve the hemodynamics and result in ATP consumption homogeneity throughout the LV, thus maximizing both CRT short-term and long-term benefits.

OBJECTIVE

The purpose of this study was to demonstrate the feasibility of optimizing CRT pacing locations to achieve maximal improvement in both ATPCTHI (an ATP consumption heterogeneity index) and stroke work.

METHODS

We used an magnetic resonance image-based electromechanical model of the dyssynchronous heart failure (DHF) canine ventricles. ATPCTHI and stroke work improvement were determined for each of 34 CRT pacing sites evenly spaced over the LV epicardium.

RESULTS

Results demonstrated the feasibility of determining the optimal LV pacing site that achieves simultaneous maximum improvements in ATPCTHI and stroke work. The optimal LV CRT pacing sites in the DHF canine ventricles were located midway between apex and base. The improvement in ATPCTHI decreased more rapidly with the distance from the optimal sites compared to stroke work improvement. CRT from the optimal sites homogenized ATP consumption by increasing septal ATP consumption and decreasing that of the lateral wall.

CONCLUSION

Simulation results using a canine heart failure model demonstrated that CRT can be optimized to achieve improvements in both ATPCTHI and stroke work.