A computational approach to understanding the cardiac electromechanical activation sequence in the normal and failing heart, with translation to the clinical practice of CRT

Use of acoustic cardiography to assess left ventricular electromechanical synchronization during left bundle branch pacing

Background

Left bundle branch pacing (LBBP) is a physiological pacing that captures the main left bundle or its proximal branch. Electromechanical activation time (EMAT) is an acoustic cardiographic metric that provides a simple method for evaluating left ventricular (LV) synchrony. Prolonged EMAT reflects impaired LV electromechanical coupling.

Objective

The purpose of this study was to explore whether EMAT can confirm that LBBP produces more satisfactory LV electromechanical synchronization than conventional right ventricular pacing modalities.

Methods

Patients with standard pacing indications and narrow QRS duration were recruited for this study. Unipolar pacing under 3 different modalities-right ventricular apical pacing (RVAP), right ventricular high septal pacing (RVHSP), and LBBP-were successively performed in each patient. Pacing parameters, echocardiographic characteristics, and acoustic cardiographic parameters at different pacing modalities and during normal rhythm were collected.

Results

A total of 55 patients were enrolled, and all had successful LBBP. Left ventricular activation time (LVAT) was significantly associated with EMAT, with LVAT vs EMAT correlation coefficient of 0.665 (P <.001). LVAT during LBBP was shorter than that during RVHSP (51.93 ± 2.732 ms vs 85.59 ± 2.240 ms; P <.001). EMAT of LBBP was significantly lower than either RVAP or RVHSP (95.44 ± 1.794 ms vs 143.32 ± 2.376 ms, and 132.22 ± 1.872 ms; both P <.001) but was similar to that of intrinsic rhythm (95.37 ± 2.271 ms; P = .862).

Conclusion

We found EMAT significantly prolonged in RVHSP and RVAP but not in the LBBP mode. This finding indicates superior electromechanical synchronization in patients having LBBP. EMAT measurement could be an additional method for identifying the ideal pacing position.

The basic heart anatomy and physiology from the cardiologist's perspective: Toward a better understanding of left ventricular mechanics, systolic, and diastolic function

A comprehensive understanding of the cardiac structure-function relationship is essential for proper clinical cardiac imaging. This review summarizes the basic heart anatomy and physiology from the perspective of a heart imager focused on myocardial mechanics. The main issues analyzed are the left ventricular (LV) architecture, the LV myocardial deformation through the cardiac cycle, the LV diastolic function basic parameters and the basic parameters of the LV deformation used in clinical practice for the LV function assessment.

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.

Desynchronization Strain Patterns and Contractility in Left Bundle Branch Block through Computer Model Simulation

Left bundle branch block (LBBB) is associated with specific septal-to-lateral wall activation patterns which are strongly influenced by the intrinsic left ventricular (LV) contractility and myocardial scar localization. The objective of this study was to propose a computational-model-based interpretation of the different patterns of LV contraction observed in the case of LBBB and preserved contractility or myocardial scarring. Two-dimensional transthoracic echocardiography was used to obtain LV volumes and deformation patterns in three patients with LBBB: (1) a patient with non-ischemic dilated cardiomyopathy, (2) a patient with antero-septal myocardial scar, and (3) a patient with lateral myocardial scar. Scar was confirmed by the distribution of late gadolinium enhancement with cardiac magnetic resonance imaging (cMRI). Model parameters were evaluated manually to reproduce patient-derived data such as strain curves obtained from echocardiographic apical views. The model was able to reproduce the specific strain patterns observed in patients. A typical septal flash with pre-ejection shortening, rebound stretch, and delayed lateral wall activation was observed in the case of non-ischemic cardiomyopathy. In the case of lateral scar, the contractility of the lateral wall was significantly impaired and septal flash was absent. In the case of septal scar, septal flash and rebound stretch were also present as previously described in the literature. Interestingly, the model was also able to simulate the specific contractile properties of the myocardium, providing an excellent localization of LV scar in ischemic patients. The model was able to simulate the electromechanical delay and specific contractility patterns observed in patients with LBBB of ischemic and non-ischemic etiology. With further improvement and validation, this technique might be a useful tool for the diagnosis and treatment planning of heart failure patients needing CRT.

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.

A rapid electromechanical model to predict reverse remodeling following cardiac resynchronization therapy

Cardiac resynchronization therapy (CRT) is an effective therapy for patients who suffer from heart failure and ventricular dyssynchrony such as left bundle branch block (LBBB). When it works, it reverses adverse left ventricular (LV) remodeling and the progression of heart failure. However, CRT response rate is currently as low as 50-65%. In theory, CRT outcome could be improved by allowing clinicians to tailor the therapy through patient-specific lead locations, timing, and/or pacing protocol. However, this also presents a dilemma: there are far too many possible strategies to test during the implantation surgery. Computational models could address this dilemma by predicting remodeling outcomes for each patient before the surgery takes place. Therefore, the goal of this study was to develop a rapid computational model to predict reverse LV remodeling following CRT. We adapted our recently developed computational model of LV remodeling to simulate the mechanics of ventricular dyssynchrony and added a rapid electrical model to predict electrical activation timing. The model was calibrated to quantitatively match changes in hemodynamics and global and local LV wall mass from a canine study of LBBB and CRT. The calibrated model was used to investigate the influence of LV lead location and ischemia on CRT remodeling outcome. Our model results suggest that remodeling outcome varies with both lead location and ischemia location, and does not always correlate with short-term improvement in QRS duration. The results and time frame required to customize and run this model suggest promise for this approach in a clinical setting.

Numerical evaluation of cardiac mechanical markers as estimators of the electrical activation time

Recent advances in the development of noninvasive cardiac imaging technologies have made it possible to measure longitudinal and circumferential strains at a high spatial resolution also at intramural level. Local mechanical activation times derived from these strains can be used as noninvasive estimates of electrical activation, in order to determine, eg, the origin of premature ectopic beats during focal arrhythmias or the pathway of reentrant circuits. The aim of this work is to assess the reliability of mechanical activation time markers derived from longitudinal and circumferential strains, denoted by AT_ell and AT_ecc , respectively, by means of three-dimensional cardiac electromechanical simulations. These markers are compared against the electrical activation time (AT_v ), computed from the action potential waveform, and the reference mechanical activation markers derived from the active tension and fiber strain waveforms, denoted by AT_ta and AT_eff , respectively. Our numerical simulations are based on a strongly coupled electromechanical model, including bidomain representation of the cardiac tissue, mechanoelectric (ie, stretch-activated channels) and geometric feedbacks, transversely isotropic strain energy function for the description of passive mechanics and detailed membrane and excitation-contraction coupling models. The results have shown that, during endocardial and epicardial ectopic stimulations, all the mechanical markers considered are highly correlated with AT_v , exhibiting correlation coefficients larger than 0.8. However, during multiple endocardial stimulations, mimicking the ventricular sinus rhythm, the mechanical markers are less correlated with the electrical activation time, because of the more complex resulting excitation sequence. Moreover, the inspection of the endocardial and epicardial isochrones has shown that the AT_ell and AT_ecc mechanical activation sequences reproduce only some qualitative features of the electrical activation sequence, such as the areas of early and late activation, but in some cases, they might yield wrong excitation sources and significantly different isochrones patterns.

Sensitivity Analysis of Ion Channel Conductance on Myocardial Electromechanical Delay: Computational Study

It is well known that cardiac electromechanical delay (EMD) can cause dyssynchronous heart failure (DHF), a prominent cardiovascular disease (CVD). This work computationally assesses the conductance variation of every ion channel on the cardiac cell to give rise to EMD prolongation. The electrical and mechanical models of human ventricular tissue were simulated, using a population approach with four conductance reductions for each ion channel. Then, EMD was calculated by determining the difference between the onset of action potential and the start of cell shortening. Finally, EMD data were put into the optimized conductance dimensional stacking to show which ion channel has the most influence in elongating the EMD. We found that major ion channels, such as L-type calcium (CaL), slow-delayed rectifier potassium (Ks), rapid-delayed rectifier potassium (Kr), and inward rectifier potassium (K1), can significantly extend the action potential duration (APD) up to 580 ms. Additionally, the maximum intracellular calcium (Cai) concentration is greatly affected by the reduction in channel CaL, Ks, background calcium, and Kr. However, among the aforementioned major ion channels, only the CaL channel can play a superior role in prolonging the EMD up to 83 ms. Furthermore, ventricular cells with long EMD have been shown to inherit insignificant mechanical response (in terms of how strong the tension can grow and how far length shortening can go) compared with that in normal cells. In conclusion, despite all variations in every ion channel conductance, only the CaL channel can play a significant role in extending EMD. In addition, cardiac cells with long EMD tend to have inferior mechanical responses due to a lack of Cai compared with normal conditions, which are highly likely to result in a compromised pump function of the heart.

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.

Optimal pacing sites in cardiac resynchronization by left ventricular activation front analysis

Cardiac resynchronization therapy (CRT) can substantially improve dyssynchronous heart failure and reduce mortality. However, about one-third of patients who are implanted, derive no measurable benefit from CRT. Non-response may partly be due to suboptimal activation of the left ventricle (LV) caused by electrophysiological heterogeneities. The goal of this study is to investigate the performance of a newly developed method used to analyze electrical wavefront propagation in a heart model including myocardial scar and compare this to clinical benchmark studies. We used computational models to measure the maximum activation front (MAF) in the LV during different pacing scenarios. Different heart geometries and scars were created based on cardiac MR images of three patients. The right ventricle (RV) was paced from the apex and the LV was paced from 12 different sites, single site, dual-site and triple site. Our results showed that for single LV site pacing, the pacing site with the largest MAF corresponded with the latest activated regions of the LV demonstrated during RV pacing, which also agrees with previous markers used for predicting optimal single-site pacing location. We then demonstrated the utility of MAF in predicting optimal electrode placements in more complex scenarios including scar and multi-site LV pacing. This study demonstrates the potential value of computational simulations in understanding and planning CRT.

Using cardiac ionic cell models to interpret clinical data

For over 100 years cardiac electrophysiology has been measured in the clinic. The electrical signals that can be measured span from noninvasive ECG and body surface potentials measurements through to detailed invasive measurements of local tissue electrophysiology. These electrophysiological measurements form a crucial component of patient diagnosis and monitoring; however, it remains challenging to quantitatively link changes in clinical electrophysiology measurements to biophysical cellular function. Multi-scale biophysical computational models represent one solution to this problem. These models provide a formal framework for linking cellular function through to emergent whole organ function and routine clinical diagnostic signals. In this review, we describe recent work on the use of computational models to interpret clinical electrophysiology signals. We review the simulation of human cardiac myocyte electrophysiology in the atria and the ventricles and how these models are being used to link organ scale function to patient disease mechanisms and therapy response in patients receiving implanted defibrillators, \cardiac resynchronisation therapy or suffering from atrial fibrillation and ventricular tachycardia. There is a growing use of multi-scale biophysical models to interpret clinical data. This allows cardiologists to link clinical observations with cellular mechanisms to better understand cardiopathophysiology and identify novel treatment strategies. This article is categorized under: Cardiovascular Diseases > Computational Models Cardiovascular Diseases > Biomedical Engineering Cardiovascular Diseases > Molecular and Cellular Physiology.

Patient-specific heart simulation can identify non-responders to cardiac resynchronization therapy

To identify non-responders to cardiac resynchronization therapy (CRT), various biomarkers have been proposed, but these attempts have not been successful to date. We tested the clinical applicability of computer simulation of CRT for the identification of non-responders. We used the multi-scale heart simulator “UT-Heart,” which can reproduce the electrophysiology and mechanics of the heart based on a molecular model of the excitation–contraction mechanism. Patient-specific heart models were created for eight heart failure patients who were treated with CRT, based on the clinical data recorded before treatment. Using these heart models, bi-ventricular pacing simulations were performed at multiple pacing sites adopted in clinical practice. Improvement in pumping function measured by the relative change of maximum positive derivative of left ventricular pressure (%ΔdP/dt_max) was compared with the clinical outcome. The operators of the simulation were blinded to the clinical outcome. In six patients, the relative reduction in end-systolic volume exceeded 15% in the follow-up echocardiogram at 3 months (responders) and the remaining two patients were judged as non-responders. The simulated %ΔdP/dt_max at the best lead position could identify responders and non-responders successfully. With further refinement of the model, patient-specific simulation could be a useful tool for identifying non-responders to CRT.

Comparison of Electromechanical Delay during Ventricular Tachycardia and Fibrillation under Different Conductivity Conditions Using Computational Modeling

Electromechanical delay (EMD) is the time interval between local myocyte depolarization and the onset of myofiber shortening. Previously, researchers measured EMD during sinus rhythm and ectopic pacing in normal and heart failure conditions. However, to our knowledge, there are no reports regarding EMD during another type of rhythms or arrhythmia. The goal of this study was to quantify EMD during sinus rhythm, tachycardia, and ventricular fibrillation conditions. We hypothesized that EMD under sinus rhythm is longer due to isovolumetric contraction which is imprecise during arrhythmia. We used a realistic model of 3D electromechanical ventricles. During sinus rhythm, EMD was measured in the last cycle of cardiac systole under steady conditions. EMD under tachycardia and fibrillation conditions was measured during the entire simulation, resulting in multiple EMD values. We assessed EMD for the following 3 conduction velocities (CVs): 31 cm/s, 51 cm/s, and 69 cm/s. The average EMD during fibrillation condition was the shortest corresponding to 53.45 ms, 55.07 ms, and 50.77 ms, for the CVs of 31 cm/s, 51 cm/s, and 69 cm/s, respectively. The average EMD during tachycardia was 58.61 ms, 58.33 ms, and 52.50 ms for the three CVs. Under sinus rhythm with action potential duration restitution (APDR) slope 0.7, the average EMD was 66.35 ms, 66.41 ms, and 66.60 ms in line with the three CVs. This result supports our hypothesis that EMD under sinus rhythm is longer than that under tachyarrhythmia conditions. In conclusion, this study observed and quantified EMD under tachycardia and ventricular fibrillation conditions. This simulation study has widened our understanding of EMD in 3D ventricles under chaotic conditions.

Cardiac Resynchronisation Therapy and Cellular Bioenergetics: Effects Beyond Chamber Mechanics

Cardiac resynchronisation therapy is a cornerstone in the treatment of advanced dyssynchronous heart failure. However, despite its widespread clinical application, precise mechanisms through which it exerts its beneficial effects remain elusive. Several studies have pointed to a metabolic component suggesting that, both in concert with alterations in chamber mechanics and independently of them, resynchronisation reverses detrimental changes to cellular metabolism, increasing energy efficiency and metabolic reserve. These actions could partially account for the existence of responders that improve functionally but not echocardiographically. This article will attempt to summarise key components of cardiomyocyte metabolism in health and heart failure, with a focus on the dyssynchronous variant. Both chamber mechanics-related and -unrelated pathways of resynchronisation effects on bioenergetics – stemming from the ultramicroscopic level – and a possible common underlying mechanism relating mechanosensing to metabolism through the cytoskeleton will be presented. Improved insights regarding the cellular and molecular effects of resynchronisation on bioenergetics will promote our understanding of non-response, optimal device programming and lead to better patient care.

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.

Electromechanical Model to Predict Cardiac Resynchronization Therapy

Cardiac resynchronization therapy (CRT) can substantially improve dyssynchronous heart failure and reduce mortality. However, one-third of the CRT patients derive no measurable benefit from CRT, due to suboptimal placement of the left ventricular (LV) lead. We introduce a pipeline for improved CRT-therapy by creating an electromechanical model using patient-specific geometric parameters allowing individualization of therapy. The model successfully mimics expected changes when variables for tension, stiffness, and conduction are entered. Changing LV pacing site had a notable effect on maximum pressure gradient (dP/dt_max) in the presence of cardiac scarring, causing non-uniform excitation propagation through the LV. Tailoring CRT to the individual requires simulations with patient-specific biventricular meshes including cardiac geometry and conductivity properties.

A review of inflammatory mechanism in airway diseases

BACKGROUND

Inflammation in the lung is the body's natural response to injury. It acts to remove harmful stimuli such as pathogens, irritants, and damaged cells and initiate the healing process. Acute and chronic pulmonary inflammation are seen in different respiratory diseases such as; acute respiratory distress syndrome, chronic obstructive pulmonary disease (COPD), asthma, and cystic fibrosis (CF).

FINDINGS

In this review, we found that inflammatory response in COPD is determined by the activation of epithelial cells and macrophages in the respiratory tract. Epithelial cells and macrophages discharge transforming growth factor-β (TGF-β), which trigger fibroblast proliferation and tissue remodeling. Asthma leads to airway hyper-responsiveness, obstruction, mucus hyper-production, and airway-wall remodeling. Cytokines, allergens, chemokines, and infectious agents are the main stimuli that activate signaling pathways in epithelial cells in asthma. Mutation of the CF transmembrane conductance regulator (CFTR) gene results in CF. Mutations in CFTR influence the lung epithelial innate immune function that leads to exaggerated and ineffective airway inflammation that fails to abolish pulmonary pathogens. We present mechanistic computational models (based on ordinary differential equations, partial differential equations and agent-based models) that have been applied in studying the complex physiological and pathological mechanisms of chronic inflammation in different airway diseases.

CONCLUSION

The scope of the present review is to explore the inflammatory mechanism in airway diseases and highlight the influence of aging on airways' inflammation mechanism. The main goal of this review is to encourage research collaborations between experimentalist and modelers to promote our understanding of the physiological and pathological mechanisms that control inflammation in different airway diseases.

Organ-level validation of a cross-bridge cycling descriptor in a left ventricular finite element model: effects of ventricular loading on myocardial strains

Although detailed cell‐based descriptors of cross‐bridge cycling have been applied in finite element (FE) heart models to describe ventricular mechanics, these multiscale models have never been tested rigorously to determine if these descriptors, when scaled up to the organ‐level, are able to reproduce well‐established organ‐level physiological behaviors. To address this void, we here validate a left ventricular (LV) FE model that is driven by a cell‐based cross‐bridge cycling descriptor against key organ‐level heart physiology. The LV FE model was coupled to a closed‐loop lumped parameter circulatory model to simulate different ventricular loading conditions (preload and afterload) and contractilities. We show that our model is able to reproduce a linear end‐systolic pressure volume relationship, a curvilinear end‐diastolic pressure volume relationship and a linear relationship between myocardial oxygen consumption and pressure–volume area. We also show that the validated model can predict realistic LV strain‐time profiles in the longitudinal, circumferential, and radial directions. The predicted strain‐time profiles display key features that are consistent with those measured in humans, such as having similar peak strains, time‐to‐peak‐strain, and a rapid change in strain during atrial contraction at late‐diastole. Our model shows that the myocardial strains are sensitive to not only LV contractility, but also to the LV loading conditions, especially to a change in afterload. This result suggests that caution must be exercised when associating changes in myocardial strain with changes in LV contractility. The methodically validated multiscale model will be used in future studies to understand human heart diseases.

A Numerical Study of Scalable Cardiac Electro-Mechanical Solvers on HPC Architectures

We introduce and study some scalable domain decomposition preconditioners for cardiac electro-mechanical 3D simulations on parallel HPC (High Performance Computing) architectures. The electro-mechanical model of the cardiac tissue is composed of four coupled sub-models: (1) the static finite elasticity equations for the transversely isotropic deformation of the cardiac tissue; (2) the active tension model describing the dynamics of the intracellular calcium, cross-bridge binding and myofilament tension; (3) the anisotropic Bidomain model describing the evolution of the intra- and extra-cellular potentials in the deforming cardiac tissue; and (4) the ionic membrane model describing the dynamics of ionic currents, gating variables, ionic concentrations and stretch-activated channels. This strongly coupled electro-mechanical model is discretized in time with a splitting semi-implicit technique and in space with isoparametric finite elements. The resulting scalable parallel solver is based on Multilevel Additive Schwarz preconditioners for the solution of the Bidomain system and on BDDC preconditioned Newton-Krylov solvers for the non-linear finite elasticity system. The results of several 3D parallel simulations show the scalability of both linear and non-linear solvers and their application to the study of both physiological excitation-contraction cardiac dynamics and re-entrant waves in the presence of different mechano-electrical feedbacks.

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.

The effect of left ventricular pacing on transmural activation delay in myopathic human hearts

Aims

Left ventricular (LV) epicardial pacing (LVEpiP) in human myopathic hearts does not decrease global epicardial activation delay compared with right ventricular (RV) endocardial pacing (RVEndoP); however, the effect on transmural activation delay has not been evaluated. To characterize the transmural electrical activation delay in human myopathic hearts during RVEndoP and LVEpiP compared with global epicardial activation delay.

Methods and results

Explanted hearts from seven patients (5 male, 46 ± 10 years) undergoing cardiac transplantation were Langendorff-perfused and mapped using an epicardial sock electrode array (112 electrodes) and 25 transmural plunge needles (four electrodes, 2 mm spacing), for a total of 100 unipolar transmural electrodes. Electrograms were recorded during LVEpiP and RVEndoP, and epicardial (sock) and transmural (needle) activation times, along with patterns of activation, were compared. There was no difference between the global epicardial activation times (LVEpiP 147 ± 8 ms vs. RVEndoP 156 ± 17 ms, P = 0.46). The mean LV transmural activation time during LVEpiP was significantly shorter than that during RVEndoP (125 ± 44 vs. 172 ± 43 ms, P < 0.001). During LVEpiP, of the transmural layers endo-, mid-myocardium and epicardium, LV endocardial layer was often the earliest compared with other transmural layers.

Conclusion

In myopathic human hearts, LVEpiP did not decrease global epicardial activation delays compared with RVEndoP. LV epicardial pacing led to early activation of the LV endocardium, revealing the importance of the LV endocardium even when pacing from the LV epicardium.

Influence of LVAD function on mechanical unloading and electromechanical delay: a simulation study

This study hypothesized that a left ventricular assist device (LVAD) shortens the electromechanical delay (EMD) by mechanical unloading. The goal of this study is to examine, by computational modeling, the influence of LVAD on EMD for four heart failure (HF) cases ranging from mild HF to severe HF. We constructed an integrated model of an LVAD-implanted cardiovascular system, then we altered the Ca²⁺ transient magnitude, with scaling factors 1, 0.9, 0.8, and 0.7 representing HF1, HF2, HF3, and HF4, respectively, in order of increasing HF severity. The four HF conditions are classified into two groups. Group one is the four HF conditions without LVAD, and group two is the conditions treated with continuous LVAD pump. The single-cell mechanical responses showed that EMD was prolonged with the higher load. The findings indicated that in group one, the HF-induced Ca2 + transient remodeling prolonged the mechanical activation time (MAT) and decreased the contractile tension, which reduced the left ventricle (LV) pressure, and increased the end-diastolic strain. In group two, LVAD shortened MAT of the ventricles. Furthermore, LVAD reduced the contractile tension, and end-diastolic strain, but increased the aortic pressure. The computational study demonstrated that LVAD shortens EMD by mechanical unloading of the ventricle.

A two dimensional electromechanical model of a cardiomyocyte to assess intra-cellular regional mechanical heterogeneities

Experimental studies on isolated cardiomyocytes from different animal species and human hearts have demonstrated that there are regional differences in the Ca2+ release, Ca2+ decay and sarcomere deformation. Local deformation heterogeneities can occur due to a combination of factors: regional/local differences in Ca2+ release and/or re-uptake, intra-cellular material properties, sarcomere proteins and distribution of the intracellular organelles. To investigate the possible causes of these heterogeneities, we developed a two-dimensional finite-element electromechanical model of a cardiomyocyte that takes into account the experimentally measured local deformation and cytosolic [Ca2+] to locally define the different variables of the constitutive equations describing the electro/mechanical behaviour of the cell. Then, the model was individualised to three different rat cardiac cells. The local [Ca2+] transients were used to define the [Ca2+]-dependent activation functions. The cell-specific local Young's moduli were estimated by solving an inverse problem, minimizing the error between the measured and simulated local deformations along the longitudinal axis of the cell. We found that heterogeneities in the deformation during contraction were determined mainly by the local elasticity rather than the local amount of Ca2+, while in the relaxation phase deformation was mainly influenced by Ca2+ re-uptake. Our electromechanical model was able to successfully estimate the local elasticity along the longitudinal direction in three different cells. In conclusion, our proposed model seems to be a good approximation to assess the heterogeneous intracellular mechanical properties to help in the understanding of the underlying mechanisms of cardiomyocyte dysfunction.

Joint influence of transmural heterogeneities and wall deformation on cardiac bioelectrical activity: A simulation study

The aim of this work is to investigate, by means of numerical simulations, the influence of myocardial deformation due to muscle contraction and relaxation on the cardiac repolarization process in presence of transmural intrinsic action potential duration (APD) heterogeneities. The three-dimensional electromechanical model considered consists of the following four coupled components: the quasi-static transversely isotropic finite elasticity equations for the deformation of the cardiac tissue; the active tension model for the intracellular calcium dynamics and cross-bridge binding; the anisotropic Bidomain model for the electrical current flow through the deforming cardiac tissue; the membrane model of ventricular myocytes, including stretch-activated channels. The numerical simulations are based on our finite element parallel solver, which employs Multilevel Additive Schwarz preconditioners for the solution of the discretized Bidomain equations and Newton-Krylov methods for the solution of the discretized non-linear finite elasticity equations. Our findings show that: (i) the presence of intrinsic transmural cellular APD heterogeneities is not fully masked by electrotonic current flow or by the presence of the mechanical deformation; (ii) despite the presence of transmural APD heterogeneities, the recovery process follows the activation sequence and there is no significant transmural repolarization gradient; (iii) with or without transmural APD heterogeneities, epicardial electrograms always display the same wave shape and discordance between the polarity of QRS complex and T-wave; (iv) the main effects of the mechanical deformation are an increase of the dispersion of repolarization time and APD, when computed over the total cardiac domain and over the endo- and epicardial surfaces, while there is a slight decrease along the transmural direction.

Tailor-made heart simulation predicts the effect of cardiac resynchronization therapy in a canine model of heart failure

Despite extensive studies on clinical indices for the selection of patient candidates for cardiac resynchronization therapy (CRT), approximately 30% of selected patients do not respond to this therapy. Herein, we examined whether CRT simulations based on individualized realistic three-dimensional heart models can predict the therapeutic effect of CRT in a canine model of heart failure with left bundle branch block. In four canine models of failing heart with dyssynchrony, individualized three-dimensional heart models reproducing the electromechanical activity of each animal were created based on the computer tomographic images. CRT simulations were performed for 25 patterns of three ventricular pacing lead positions. Lead positions producing the best and the worst therapeutic effects were selected in each model. The validity of predictions was tested in acute experiments in which hearts were paced from the sites identified by simulations. We found significant correlations between the experimentally observed improvement in ejection fraction (EF) and the predicted improvements in ejection fraction (P<0.01) or the maximum value of the derivative of left ventricular pressure (P<0.01). The optimal lead positions produced better outcomes compared with the worst positioning in all dogs studied, although there were significant variations in responses. Variations in ventricular wall thickness among the dogs may have contributed to these responses. Thus CRT simulations using the individualized three-dimensional heart models can predict acute hemodynamic improvement, and help determine the optimal positions of the pacing lead.

Lessons learned from multi-scale modeling of the failing heart

Heart failure constitutes a major public health problem worldwide. Affected patients experience a number of changes in the electrical function of the heart that predispose to potentially lethal cardiac arrhythmias. Due to the multitude of electrophysiological changes that may occur during heart failure, the scientific literature is complex and sometimes ambiguous, perhaps because these findings are highly dependent on the etiology, the stage of heart failure, and the experimental model used to study these changes. Nevertheless, a number of common features of failing hearts have been documented. Prolongation of the action potential (AP) involving ion channel remodeling and alterations in calcium handling have been established as the hallmark characteristics of myocytes isolated from failing hearts. Intercellular uncoupling and fibrosis are identified as major arrhythmogenic factors. Multi-scale computational simulations are a powerful tool that complements experimental and clinical research. The development of biophysically detailed computer models of single myocytes and cardiac tissues has contributed greatly to our understanding of processes underlying excitation and repolarization in the heart. The electrical, structural, and metabolic remodeling that arises in cardiac tissues during heart failure has been addressed from different computational perspectives to further understand the arrhythmogenic substrate. This review summarizes the contributions from computational modeling and simulation to predict the underlying mechanisms of heart failure phenotypes and their implications for arrhythmogenesis, ranging from the cellular level to whole-heart simulations. The main aspects of heart failure are presented in several related sections. An overview of the main electrophysiological and structural changes that have been observed experimentally in failing hearts is followed by the description and discussion of the simulation work in this field at the cellular level, and then in 2D and 3D cardiac structures. The implications for arrhythmogenesis in heart failure are also discussed including therapeutic measures, such as drug effects and cardiac resynchronization therapy. Finally, the future challenges in heart failure modeling and simulation will be discussed.

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.

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.

Relationship between mechanical dyssynchrony and intra-operative electrical delay times in patients undergoing cardiac resynchronization therapy

BACKGROUND

It is important to understand the relationship between electrical and mechanical ventricular activation in CRT patients. By measuring local electrical activation at multiple locations within the coronary veins and myocardial contraction at the same locations in the left ventricle, we determined the relationship between electrical and mechanical activation at potential left ventricular pacing locations.

METHODS

In this study, mechanical contraction times were computed using high temporal resolution cine cardiovascular magnetic resonance (CMR) data, while electrical activation times were derived from intra-procedural local electrograms.

RESULTS

In our cohort, there was a strong correlation between electrical and mechanical delay times within each patient (R2=0.78 ± 0.23). Additionally, the latest electrically activated location corresponded with the latest mechanically contracting location in 91% of patients.

CONCLUSIONS

This study provides initial evidence that our method of obtaining non-invasive mechanical activation patterns accurately reflects the underlying electromechanical substrate of intraventricular dyssynchrony.

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

BACKGROUND

The acute response to cardiac resynchronization therapy (CRT) has been shown to be due to 3 mechanisms: resynchronization of ventricular contraction, efficient preloading of the ventricles by a properly timed atrial contraction, and mitral regurgitation reduction. However, the contribution of each of the 3 mechanisms to the acute response to CRT, specifically stroke work improvement, has not been quantified.

OBJECTIVE

To use a magnetic resonance image-based anatomically accurate 3-dimensional model of failing canine ventricular electromechanics to quantify the contribution of each of the 3 mechanisms to stroke work improvement and identify the predominant mechanisms.

METHODS

An MRI-based electromechanical model of the failing canine ventricles assembled previously by our group was further developed and modified. Three different protocols were used to dissect the contribution of each of the 3 mechanisms to stroke work improvement.

RESULTS

Resynchronization of ventricular contraction did not lead to a significant stroke work improvement. Efficient preloading of the ventricles by a properly timed atrial contraction was the predominant mechanism underlying stroke work improvement. Stroke work improvement peaked at an intermediate atrioventricular delay, as it allowed ventricular filling by atrial contraction to occur at a low diastolic left ventricular pressure but also provided adequate time for ventricular filling before ventricular contraction. Reduction of mitral regurgitation by CRT led to stroke work worsening instead of improvement.

CONCLUSION

Efficient preloading of the ventricles by a properly timed atrial contraction is responsible for a significant stroke work improvement in the acute CRT response.

Stretch-activated BK channel and heart function

The heart is an organ that is exposed to extreme dynamic mechanical stimuli. From birth till death, the heart indefinitely repeats periodic contraction and dilation, i.e., shortening and elongation of cardiomyocytes. Mechanical stretch elicits a change in heart rate and may cause arrhythmia if it is excessive. Thus, mechanosensitivity is crucial to heart function. The molecule that is substantially involved in mechanosensitivity is a stretch-activated ion channel. Among several ion channels believed to be activated by stretch in the heart, the stretch-activated KCa (SAKCA) channel, a member of the group of large conductance (Big Potassium, BK) channels, shows a mechanosensitive (MS) response to membrane stretch. As BK channels respond to voltage and intracellular calcium concentration with large conductance, they are considered to be involved in repolarization after depolarization. Some BK channels are known to be activated by stretch and are expressed in a number of cells, including human osteoblasts and guinea pig intestinal neurons. The SAKCA channel was found to be sensitive to stretch in the chick heart. Given that the cardiomyocyte is unremittingly exposed to contraction and dilation and that it generates action potential and its contractility is modulated by intracellular calcium concentration, the SAKCA channel, which is dependent voltage and calcium, may be involved in action potential generation. It was recently reported that a BK channel is involved in the modulation of heart rate in the mouse. Further studies regarding the role of MS BK channels, including SAKCA, in the modulation of heart rate and contractility are expected.

Mechanistic insight into prolonged electromechanical delay in dyssynchronous heart failure: a computational study

In addition to the left bundle branch block type of electrical activation, there are further remodeling aspects associated with dyssynchronous heart failure (HF) that affect the electromechanical behavior of the heart. Among the most important are altered ventricular structure (both geometry and fiber/sheet orientation), abnormal Ca(2+) handling, slowed conduction, and reduced wall stiffness. In dyssynchronous HF, the electromechanical delay (EMD), the time interval between local myocyte depolarization and myofiber shortening onset, is prolonged. However, the contributions of the four major HF remodeling aspects in extending EMD in the dyssynchronous failing heart remain unknown. The goal of this study was to determine the individual and combined contributions of HF-induced remodeling aspects to EMD prolongation. We used MRI-based models of dyssynchronous nonfailing and HF canine electromechanics and constructed additional models in which varying combinations of the four remodeling aspects were represented. A left bundle branch block electrical activation sequence was simulated in all models. The simulation results revealed that deranged Ca(2+) handling is the primary culprit in extending EMD in dyssynchronous HF, with the other aspects of remodeling contributing insignificantly. Mechanistically, we found that abnormal Ca(2+) handling in dyssynchronous HF slows myofiber shortening velocity at the early-activated septum and depresses both myofiber shortening and stretch rate at the late-activated lateral wall. These changes in myofiber dynamics delay the onset of myofiber shortening, thus giving rise to prolonged EMD in dyssynchronous HF.

Computational cardiology: how computer simulations could be used to develop new therapies and advance existing ones

This article reviews the latest developments in computational cardiology. It focuses on the contribution of cardiac modelling to the development of new therapies as well as the advancement of existing ones for cardiac arrhythmias and pump dysfunction. Reviewed are cardiac modelling efforts aimed at advancing and optimizing existent therapies for cardiac disease (defibrillation, ablation of ventricular tachycardia, and cardiac resynchronization therapy) and at suggesting novel treatments, including novel molecular targets, as well as efforts to use cardiac models in stratification of patients likely to benefit from a given therapy, and the use of models in diagnostic procedures.

Effects of mechano-electric feedback on scroll wave stability in human ventricular fibrillation

Recruitment of stretch-activated channels, one of the mechanisms of mechano-electric feedback, has been shown to influence the stability of scroll waves, the waves that underlie reentrant arrhythmias. However, a comprehensive study to examine the effects of recruitment of stretch-activated channels with different reversal potentials and conductances on scroll wave stability has not been undertaken; the mechanisms by which stretch-activated channel opening alters scroll wave stability are also not well understood. The goals of this study were to test the hypothesis that recruitment of stretch-activated channels affects scroll wave stability differently depending on stretch-activated channel reversal potential and channel conductance, and to uncover the relevant mechanisms underlying the observed behaviors. We developed a strongly-coupled model of human ventricular electromechanics that incorporated human ventricular geometry and fiber and sheet orientation reconstructed from MR and diffusion tensor MR images. Since a wide variety of reversal potentials and channel conductances have been reported for stretch-activated channels, two reversal potentials, −60 mV and −10 mV, and a range of channel conductances (0 to 0.07 mS/µF) were implemented. Opening of stretch-activated channels with a reversal potential of −60 mV diminished scroll wave breakup for all values of conductances by flattening heterogeneously the action potential duration restitution curve. Opening of stretch-activated channels with a reversal potential of −10 mV inhibited partially scroll wave breakup at low conductance values (from 0.02 to 0.04 mS/µF) by flattening heterogeneously the conduction velocity restitution relation. For large conductance values (>0.05 mS/µF), recruitment of stretch-activated channels with a reversal potential of −10 mV did not reduce the likelihood of scroll wave breakup because Na channel inactivation in regions of large stretch led to conduction block, which counteracted the increased scroll wave stability due to an overall flatter conduction velocity restitution.

Understanding the mechanisms amenable to CRT response: from pre-operative multimodal image data to patient-specific computational models

This manuscript describes our recent developments towards better understanding of the mechanisms amenable to cardiac resynchronization therapy response. We report the results from a full multimodal dataset corresponding to eight patients from the euHeart project. The datasets include echocardiography, MRI and electrophysiological studies. We investigate two aspects. The first one focuses on pre-operative multimodal image data. From 2D echocardiography and 3D tagged MRI images, we compute atlas based dyssynchrony indices. We complement these indices with presence and extent of scar tissue and correlate them with CRT response. The second one focuses on computational models. We use pre-operative imaging to generate a patient-specific computational model. We show results of a fully automatic personalized electromechanical simulation. By case-per-case discussion of the results, we highlight the potential and key issues of this multimodal pipeline for the understanding of the mechanisms of CRT response and a better patient selection.

Computational cardiology: the heart of the matter

This paper reviews the newest developments in computational cardiology. It focuses on the contribution of cardiac modeling to the development of new therapies as well as the advancement of existing ones for cardiac arrhythmias and pump dysfunction. Reviewed are cardiac modeling efforts aimed at advancing and optimizing existent therapies for cardiac disease (defibrillation, ablation of ventricular tachycardia, and cardiac resynchronization therapy) and at suggesting novel treatments, including novel molecular targets, as well as efforts to use cardiac models in stratification of patients likely to benefit from a given therapy, and the use of models in diagnostic procedures.