Dyssynchronous heart failure (DHF) is a severe form of heart failure where conduction block in the left bundle branch causes delayed left ventricular electrical activation and discoordinated systolic contraction, dramatically reducing cardiac output. Cardiac resynchronization therapy (CRT) is a cost effective pacing treatment that has been shown to improve symptoms and survival, especially due to left ventricular reverse remodeling. However, approximately 50% of patients do not show objective evidence of reverse remodeling even after 6 months of CRT. A deeper understanding of the physiological mechanisms leading to positive long-term outcomes and identification of patients who are most likely to benefit are needed to maximize quality of care and minimize health risks and economic costs. The ability to predict the outcome and personalize CRT application for an individual patient from clinical measurements alone is challenging given the wide inter-patient variability of clinical features and pathophysiological complexity of DHF. In this work, we seek to answer questions regarding physiological mechanisms that are implicated in CRT response, baseline physiological features that are predictive of response, and personalized CRT application for an individual patient. For this purpose, we construct patient-specific computational models of DHF which integrate anatomical, electrophysiological, biomechanical, and hemodynamic clinical and empirical data to quantitatively characterize baseline and CRT physiology to understand how patients differ in response. The primary aims of this thesis will be to : 1) construct patient-specific computational models of DHF incorporating clinical and empirical measurements to test whether the models can recapitulate characteristics of DHF and predict measured acute effects of CRT; 2) test the hypothesis that CRT response physiologically depends on the severity of baseline heterogeneity of mechanical loading caused by electrical dyssynchrony and ventricular dilation; 3) test the hypothesis that CRT response can be predicted from novel model-derived biomarkers of electrical dyssynchrony. Through quantification and prediction of patient-specific cardiovascular physiology in disease and therapy, computational models have great potential to enhance the quality of medical care by providing novel diagnostic value to support clinical decisions regarding the best personalized approach to treat the individual patient

Computational modeling of tissue-scale cardiac electrophysiology requires numerically converged solutions to avoid spurious artifacts. The steep gradients inherent to cardiac action potential propagation necessitate fine spatial scales and therefore a substantial computational burden. The use of high-order interpolation methods has previously been proposed for these simulations due to their theoretical convergence advantage. In this study, we compare the convergence behavior of linear Lagrange, cubic Hermite, and the newly proposed cubic Hermite-style serendipity interpolation methods for finite element simulations of the cardiac monodomain equation. The high-order methods reach converged solutions with fewer degrees of freedom and longer element edge lengths than traditional linear elements. Additionally, we propose a dimensionless number, the cell Thiele modulus, as a more useful metric for determining solution convergence than element size alone. Finally, we use the cell Thiele modulus to examine convergence criteria for obtaining clinically useful activation patterns for applications such as patient-specific modeling where the total activation time is known a priori.

Patient-specific models of the ventricular myocardium, combined with the computational power to run rapid simulations, are approaching the level where they could be used for personalized cardiovascular medicine. A major remaining challenge is determining model parameters from available patient data, especially for models of the Purkinje-myocardial junctions (PMJs): the sites of initial ventricular electrical activation. There are no non-invasive methods for localizing PMJs in patients, and the relationship between the standard clinical ECG and PMJ model parameters is underexplored. Thus, this study aimed to determine the sensitivity of the QRS complex of the ECG to the anatomical location and regional number of PMJs. The QRS complex was simulated using an image-based human torso and biventricular model, and cardiac electrophysiology was simulated using Cardioid. The PMJs were modeled as discrete current injection stimuli, and the location and number of stimuli were varied within initial activation regions based on published experiments. Results indicate that the QRS complex features were most sensitive to the presence or absence of four "seed" stimuli, and adjusting locations of nearby "regional" stimuli provided finer tuning. Decreasing number of regional stimuli by an order of magnitude resulted in virtually no change in the QRS complex. Thus, a minimal 12-stimuli configuration was identified that resulted in physiological excitation, defined by QRS complex feature metrics and ventricular excitation pattern. Overall, the sensitivity results suggest that parameterizing PMJ location, rather than number, be given significantly higher priority in future studies creating personalized ventricular models from patient-derived ECGs.

In patients with dyssynchronous heart failure (DHF), cardiac conduction abnormalities cause the regional distribution of myocardial work to be non-homogeneous. Cardiac resynchronization therapy (CRT) using an implantable, programmed biventricular pacemaker/defibrillator, can improve the synchrony of contraction between the right and left ventricles in DHF, resulting in reduced morbidity and mortality and increased quality of life. Since regional work depends on wall stress, which cannot be measured in patients, we used computational methods to investigate regional work distributions and their changes after CRT. We used three-dimensional multi-scale patient-specific computational models parameterized by anatomic, functional, hemodynamic, and electrophysiological measurements in eight patients with heart failure and left bundle branch block (LBBB) who received CRT. To increase clinical translatability, we also explored whether streamlined computational methods provide accurate estimates of regional myocardial work. We found that CRT increased global myocardial work efficiency with significant improvements in non-responders. Reverse ventricular remodeling after CRT was greatest in patients with the highest heterogeneity of regional work at baseline, however the efficacy of CRT was not related to the decrease in overall work heterogeneity or to the reduction in late-activated regions of high myocardial work. Rather, decreases in early-activated regions of myocardium performing negative myocardial work following CRT best explained patient variations in reverse remodeling. These findings were also observed when regional myocardial work was estimated using ventricular pressure as a surrogate for myocardial stress and changes in endocardial surface area as a surrogate for strain. These new findings suggest that CRT promotes reverse ventricular remodeling in human dyssynchronous heart failure by increasing regional myocardial work in early-activated regions of the ventricles, where dyssynchrony is specifically associated with hypoperfusion, late systolic stretch, and altered metabolic activity and that measurement of these changes can be performed using streamlined approaches.

Background

Stereotactic ablative radiotherapy (SAbR) is an emerging therapy for refractory ventricular tachycardia (VT). However, the current workflow is complicated, and the precision and safety in patients with significant cardiorespiratory motion and VT targets near the stomach may be suboptimal.

Objective

We hypothesized that automated 12-lead electrocardiogram (ECG) mapping and respiratory-gated therapy may improve the ease and precision of SAbR planning and facilitate safe radiation delivery in patients with refractory VT.

Methods

Consecutive patients with refractory VT were studied at 2 hospitals. VT exit sites were localized using a 3-D computational ECG algorithm noninvasively and compared to available prior invasive mapping. Radiotherapy (25 Gy) was delivered at end-expiration when cardiac respiratory motion was ≥0.6 cm or targets were ≤2 cm from the stomach.

Results

In 6 patients (ejection fraction 29% ± 13%), 4.2 ± 2.3 VT morphologies per patient were mapped. Overall, 7 out of 7 computational ECG mappings (100%) colocalized to the identical cardiac segment when prior invasive electrophysiology study was available. Respiratory gating was associated with smaller planning target volumes compared to nongated volumes (71 ± 7 vs 153 ± 35 cc, P < .01). In 2 patients with inferior wall VT targets close to the stomach (6 mm proximity) or significant respiratory motion (22 mm excursion), no GI complications were observed at 9- and 12-month follow-up. Implantable cardioverter-defibrillator shocks decreased from 23 ± 12 shocks/patient to 0.67 ± 1.0 (P < .001) post-SAbR at 6.0 ± 4.9 months follow-up.

Conclusions

A workflow including computational ECG mapping and protocol-guided respiratory gating is feasible, is safe, and may improve the ease of SAbR planning. Studies to validate this workflow in larger populations are required.