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

Progress of Conductivity and Conduction Velocity Measured in Human and Animal Hearts

Cardiac conduction velocity (CV) is a critical electrophysiological characteristic of the myocardium, representing the speed at which electrical pulses propagate through cardiac tissue. It can be delineated into longitudinal, transverse, and normal components in the myocardium. The CV and its anisotropy ratio are crucial to both normal electrical conduction and myocardial contraction, as well as pathological conditions where it increases the risk of conduction block and reentry. This comprehensive review synthesizes longitudinal and transverse CV values from clinical and experimental studies of human infarct hearts, including findings from the isthmus and outer loop, alongside data derived from animal models. Additionally, we explore the anisotropic ratio of conductivities assessed through both animal and computational models. The review culminates with a synthesis of scientific evidence that guides the selection of CV and its corresponding conductivity in cardiac modeling, particularly emphasizing its application in patient-specific cardiac arrhythmia modeling.

Combining pharmacokinetic and electrophysiological models for early prediction of drug-induced arrhythmogenicity

BACKGROUND AND OBJECTIVE

In silico methods are gaining attention for predicting drug-induced Torsade de Pointes (TdP) in different stages of drug development. However, many computational models tended not to account for inter-individual response variability due to demographic covariates, such as sex, or physiologic covariates, such as renal function, which may be crucial when predicting TdP. This study aims to compare the effects of drugs in male and female populations with normal and impaired renal function using in silico methods.

METHODS

Pharmacokinetic models considering sex and renal function as covariates were implemented from data published in pharmacokinetic studies. Drug effects were simulated using an electrophysiologically calibrated population of cellular models of 300 males and 300 females. The population of models was built by modifying the endocardial action potential model published by O'Hara et al. (2011) according to the experimentally measured gene expression levels of 12 ion channels.

RESULTS

Fifteen pharmacokinetic models for CiPA drugs were implemented and validated in this study. Eight pharmacokinetic models included the effect of renal function and four the effect of sex. The mean difference in action potential duration (APD) between male and female populations was 24.9 ms (p<0.05). Our simulations indicated that women with impaired renal function were particularly susceptible to drug-induced arrhythmias, whereas healthy men were less prone to TdP. Differences between patient groups were more pronounced for high TdP-risk drugs. The proposed in silico tool also revealed that individuals with impaired renal function, electrophysiologically simulated with hyperkalemia (extracellular potassium concentration [K+]o = 7 mM) exhibited less pronounced APD prolongation than individuals with normal potassium levels. The pharmacokinetic/electrophysiological framework was used to determine the maximum safe dose of dofetilide in different patient groups. As a proof of concept, 3D simulations were also run for dofetilide obtaining QT prolongation in accordance with previously reported clinical values.

CONCLUSIONS

This study presents a novel methodology that combines pharmacokinetic and electrophysiological models to incorporate the effects of sex and renal function into in silico drug simulations and highlights their impact on TdP-risk assessment. Furthermore, it may also help inform maximum dose regimens that ensure TdP-related safety in a specific sub-population of patients.

Time-related factors predicting a positive response to cardiac resynchronisation therapy in patients with heart failure

This study aimed to identify time parameters predicting favourable CRT response. A total of 38 patients with ischemic cardiomyopathy, qualified for CRT implantation, were enrolled in the study. A 15% reduction in indexed end-systolic volume after 6 months was a criterion for a positive response to CRT. We evaluated QRS duration, measured from a standard ECG before and after CRT implantation and obtained from mapping with NOGA XP system (AEMM); and the delay, measured with the implanted device algorithm (DCD) and its change after 6 months (ΔDCD); and selected delay parameters between the left and right ventricles based on AEMM data. A total of 24 patients presented with a positive response to CRT versus 9 non-responders. After CRT implantation, we observed differences between responders and non-responders group in the reduction of QRS duration (31 ms vs. 16 ms), duration of paced QRS (123 ms vs. 142 ms), and the change of ΔDCDMaximum (4.9 ms vs. 0.44 ms) and ΔDCDMean (7.7 ms vs. 0.9 ms). The difference in selected parameters obtained during AEMM in both groups was related to interventricular delay (40.3 ms vs. 18.6 ms). Concerning local activation time and left ventricular activation time, we analysed the delays in individual left ventricular segments. Predominant activation delay of the posterior wall middle segment was associated with a better response to CRT. Some AEMM parameters, paced QRS time of less than 120 ms and reduction of QRS duration greater than 20 ms predict the response to CRT. ΔDCD is associated with favourable electrical and structural remodelling.Clinical trial registration: SUM No. KNW/0022/KB1/17/15.

The mechanisms of potassium loss in acute myocardial ischemia: New insights from computational simulations

Acute myocardial ischemia induces hyperkalemia (accumulation of extracellular potassium), a major perpetrator of lethal reentrant ventricular arrhythmias. Despite considerable experimental efforts to explain this pathology in the last decades, the intimate mechanisms behind hyperkalemia remain partially unknown. In order to investigate these mechanisms, we developed a novel computational model of acute myocardial ischemia which couples a) an electrophysiologically detailed human cardiomyocyte model that incorporates modifications to account for ischemia-induced changes in transmembrane currents, with b) a model of cardiac tissue and extracellular K ⁺ transport. The resulting model is able to reproduce and explain the triphasic time course of extracellular K ⁺ concentration within the ischemic zone, with values of [ K + ] o close to 14 mmol/L in the central ischemic zone after 30 min. In addition, the formation of a [ K + ] o border zone of approximately 1.2 cm 15 min after the onset of ischemia is predicted by the model. Our results indicate that the primary rising phase of [ K + ] o is mainly due to the imbalance between K ⁺ efflux, that increases slightly, and K ⁺ influx, that follows a reduction of the NaK pump activity by more than 50%. The onset of the plateau phase is caused by the appearance of electrical alternans (a novel mechanism identified by the model), which cause an abrupt reduction in the K ⁺ efflux. After the plateau, the secondary rising phase of [ K + ] o is caused by a subsequent imbalance between the K ⁺ influx, which continues to decrease slowly, and the K ⁺ efflux, which remains almost constant. Further, the study shows that the modulation of these mechanisms by the electrotonic coupling is the main responsible for the formation of the ischemic border zone in tissue, with K ⁺ transport playing only a minor role. Finally, the results of the model indicate that the injury current established between the healthy and the altered tissue is not sufficient to depolarize non-ischemic cells within the healthy tissue.

Cardiac Resynchronization Therapy for Adult Patients With a Failing Systemic Right Ventricle: A Multicenter Study

Background

The objective of this international multicenter study was to investigate both early and late outcomes of cardiac resynchronization therapy (CRT) in patients with a systemic right ventricle (SRV) and to identify predictors for congestive heart failure readmissions and mortality.

Methods and Results

This retrospective international multicenter study included 13 centers. The study population comprised 80 adult patients with SRV (48.9% women) with a mean age of 45±14 (range, 18–77) years at initiation of CRT. Median follow‐up time was 4.1 (25th–75th percentile, 1.3–8.3) years. Underlying congenital heart disease consisted of congenitally corrected transposition of the great arteries and dextro‐transposition of the great arteries in 63 (78.8%) and 17 (21.3%) patients, respectively. CRT resulted in significant improvement in functional class (before CRT: III, 25th–75th percentile, II–III; after CRT: II, 25th–75th percentile, II–III; P =0.005) and QRS duration (before CRT: 176±27; after CRT: 150±24 milliseconds; P =0.003) in patients with pre‐CRT ventricular pacing who underwent an upgrade to a CRT device (n=49). These improvements persisted during long‐term follow‐up with a marginal but significant increase in SRV function (before CRT; 30%, 25th–75th percentile, 25–35; after CRT: 31%, 25th–75th percentile, 21–38; P =0.049). In contrast, no beneficial change in the above‐mentioned variables was observed in patients who underwent de novo CRT (n=31). A quarter of all patients were readmitted for heart failure during follow‐up, and mortality at latest follow‐up was 21.3%.

Conclusions

This international experience with CRT in patients with an SRV demonstrated that CRT in selected patients with SRV dysfunction and pacing‐induced dyssynchrony yielded consistent improvement in QRS duration and New York Heart Association functional status, with a marginal increase in SRV function.

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

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

Optimization of Left Ventricle Pace Maker Location Using Echo-Based Fluid-Structure Interaction Models

INTRODUCTION

Cardiac pacing has been an effective treatment in the management of patients with bradyarrhythmia and tachyarrhythmia. Different pacemaker location has different responses, and pacemaker effectiveness to each individual can also be different. A novel image-based ventricle animal modeling approach was proposed to optimize ventricular pacemaker site for better cardiac outcome.

METHOD

One health female adult pig (weight 42.5 kg) was used to make a pacing animal model with different ventricle pacing locations. Ventricle surface electric signal, blood pressure and echo image were acquired 15 min after the pacemaker was implanted. Echo-based left ventricle fluid-structure interaction models were constructed to perform ventricle function analysis and investigate impact of pacemaker location on cardiac outcome. With the measured electric signal map from the pig associated with the actual pacemaker site, electric potential conduction of myocardium was modeled by material stiffening and softening in our model, with stiffening simulating contraction and softening simulating relaxation. Ventricle model without pacemaker (NP model) and three ventricle models with the following pacemaker locations were simulated: right ventricular apex (RVA model), posterior interventricular septum (PIVS model) and right ventricular outflow tract (RVOT model). Since higher peak flow velocity, flow shear stress (FSS), ventricle stress and strain are linked to better cardiac function, those data were collected for model comparisons.

RESULTS

At the peak of filling, velocity magnitude, FSS, stress and strain for RVOT and PIVS models were 13%, 45%, 18%, 13% and 5%, 30%, 10%, 5% higher than NP model, respectively. At the peak of ejection, velocity magnitude, FSS, stress and strain for RVOT and PIVS models were 50%, 44%, 54%, 59% and 23%, 36%, 39%, 53% higher than NP model, respectively. RVA model had lower velocity, FSS, stress and strain than NP model. RVOT model had higher peak flow velocity and stress/strain than PIVS model. It indicated RVOT pacemaker site may be the best location.

CONCLUSION

This preliminary study indicated that RVOT model had the best performance among the four models compared. This modeling approach could be used as "virtual surgery" to try various pacemaker locations and avoid risky and dangerous surgical experiments on real patients.

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.

Analysis of vulnerability to reentry in acute myocardial ischemia using a realistic human heart model

Electrophysiological alterations of the myocardium caused by acute ischemia constitute a pro-arrhythmic substrate for the generation of potentially lethal arrhythmias. Experimental evidence has shown that the main components of acute ischemia that induce these electrophysiological alterations are hyperkalemia, hypoxia (or anoxia in complete artery occlusion), and acidosis. However, the influence of each ischemic component on the likelihood of reentry is not completely established. Moreover, the role of the His-Purkinje system (HPS) in the initiation and maintenance of arrhythmias is not completely understood. In the present work, we investigate how the three components of ischemia affect the vulnerable window (VW) for reentry using computational simulations. In addition, we analyze the role of the HPS on arrhythmogenesis. A 3D biventricular/torso human model that includes a realistic geometry of the central and border ischemic zones with one of the most electrophysiologically detailed model of ischemia to date, as well as a realistic cardiac conduction system, were used to assess the VW for reentry. Four scenarios of ischemic severity corresponding to different minutes after coronary artery occlusion were simulated. Our results suggest that ischemic severity plays an important role in the generation of reentries. Indeed, this is the first 3D simulation study to show that ventricular arrhythmias could be generated under moderate ischemic conditions, but not in mild and severe ischemia. Moreover, our results show that anoxia is the ischemic component with the most significant effect on the width of the VW. Thus, a change in the level of anoxia from moderate to severe leads to a greater increment in the VW (40 ms), in comparison with the increment of 20 ms and 35 ms produced by the individual change in the level of hyperkalemia and acidosis, respectively. Finally, the HPS was a necessary element for the generation of approximately 17% of reentries obtained. The retrograde conduction from the myocardium to HPS in the ischemic region, conduction blocks in discrete sections of the HPS, and the degree of ischemia affecting Purkinje cells, are suggested as mechanisms that favor the generation of ventricular arrhythmias.

Prevention of non-response to cardiac resynchronization therapy: points to remember

Cardiac resynchronization therapy (CRT) is an important and effective therapy for end-stage heart failure. Non-response to CRT is one of the main obstacles to its application in clinical practice. There is no uniform consensus or definition of CRT “response.” Clinical symptoms, ventricular remodeling indices, and cardiovascular events have been reported to be associated with non-responders. To prevent non-response to CRT, three aspects should be thoroughly considered: preoperative patient selection, electrode implantation, and postoperative management. Preoperative selection of appropriate patients for CRT treatment is an important step in preventing non-response. Currently, the CRT inclusion criteria are mainly based on the morphology of QRS waves in deciding ventricular dyssynchrony. Echocardiography and cardiac magnetic resonance are being explored to predict nonresponse to CRT. The location of left ventricular electrode implantation is a current hot spot of research; it is important to identify the location of the latest exciting ventricular segment and avoid scars. Cardiac magnetic resonance and ultrasonic spot tracking are being progressively developed in this field. Some new techniques such as His Bundle pacing, endocardial electrodes, and novel sensors are also being investigated. Postoperative management of patients is another essential step towards preventing non-response; it mainly focuses on the treatment of the disease itself and CRT program control optimization. CRT treatment is just one part of the overall treatment of heart failure, and multidisciplinary efforts are needed to improve the overall outcome.