Therapeutic Hypothermia Reduces Peritoneal Dialysis Induced Myocardial Blood Flow Heterogeneity and Arrhythmia

Background: Moderate therapeutic hypothermia (TH) is a well-recognized cardio-protective strategy. The instillation of fluid into the peritoneum provides an opportunity to deliver moderate hypothermia as primary prevention against cardiovascular events. We aimed to to investigate both cardiac perfusion consequences (overall blood flow and detailed assessment of perfusion heterogeneity) and subsequently simulate the associated arrhythmic risk for patients undergoing peritoneal dialysis (PD) induced TH. Methods: Patients underwent high resolution myocardial perfusion scanning using high resolution 256 slice CT scanning, at rest and with adenosine stress. The first visit using the patient's usual PD regimen, on the second visit the same regime was utilized but with cooled peritoneal dialysate at 32°C. Myocardial blood flow (MBF) was quantified from generated perfusion maps, reconstructed in 3D. MBF heterogeneity was assessed by fractal dimension (FD) measurement on the 3D left ventricular reconstruction. Arrhythmogenicity was quantified from a sophisticated computational simulation using a multi-scale human 3D ventricle wedge electrophysiological computational model. Results: We studied 7 PD patients, mean age of 60 ± 7 and mean vintage dialysis of 23.6 ± 17.6 months. There were no significant different in overall segmental MBF between normothermic condition (NT) and TH. MBF heterogeneity was significantly decreased (-14%, p = 0.03) at rest and after stress (-14%, p = 0.03) when cooling was applied. Computational simulation showed that TH allowed a normalization of action potential, QT duration and T wave. Conclusion: TH-PD results in moderate hypothermia leading to a reduction in perfusion heterogeneity and simulated risk of non-terminating malignant ventricular arrhythmias.

Proarrhythmic Effects of Electrolyte Imbalance in Virtual Human Atrial and Ventricular Cardiomyocytes

Dialysis is prescribed to renal failure patients as a long-term chronic treatment. Whereas dialysis therapeutically normalizes serum electrolytes and removes small toxin molecules, it fails to alleviate fibroblast induced structural fibrosis, and unresponsive uremia. The simultaneous presence of altered electrolytes and fibrosis or uremia is thought to be pro-arrhythmogenic. This study explored potential arrhythmogenesis under pre-dialysis (high electrolyte levels) and post-dialysis (low physiological electrolyte levels) in the presence of fibrosis and uremia in human atrial and ventricular model cardiomyocytes.Two validated human cardiomyocyte models were used in this study that permitted simulation of cardiac atrial and ventricular detailed electrophysiology. Pathological conditions simulating active fibrosis and uremia were implemented in both models. Pre- and post-dialysis conditions were simulated using high and low electrolyte levels respectively. Arrythmogenesis was quantified by computing restitution curves that permitted identification of action potential duration and calcium transient alternans instabilities.In comparison to control conditions, fibrosis abbreviated action potential durations while uremia prolonged the same. Under pre-dialysis conditions, an elevation of serum electrolyte levels caused action potential durations to be abbreviated under both fibrosis and uremia. Alternans instability was observed in the ventricular cardiomyocyte model. Under post-dialysis conditions, lower levels of serum electrolytes promoted an abbreviated action potential duration under fibrosis but caused a large increase of the control and uremic action potential durations. Alternans instabilities were observed in the atrial cardiomyocyte model under post-dialysis conditions at physiological heart rates. The calcium transient restitution showed similar alternans instabilities.Co-existing conditions such as fibrosis and uremia in the presence of unphysiological electrolyte levels promote arrhythmogenesis and may require additional treatment to improve dialysis outcomes.Clinical Relevance. Knowledge of model response to clinically relevant conditions permits use of in silico modeling to better understand and dissect underlying arrhythmia mechanisms.

Elucidating the relationship between arrhythmia and ischemic heterogeneity: an in silico study

Dialysis causes blood flow defects in the heart that may augment electrophysiological heterogeneity in the form of increased number of ischemic zones in the human left ventricle. We computationally tested whether a larger number of ischemic zones aggravate arrhythmia using a 2D electrophysiological model of the human ventricle.A human ventricle cardiomyocyte model capable of simulating ischemic action potentials was adapted in this study. The cell model was incorporated into a spatial 2D model consisting of known number of ischemic zones. Inter-cellular gap junction coupling within ischemic zones was reduced to simulate slow conduction. Arrhythmia severity was assessed by inducing a re-entry, and quantifying the ensuing breakup and tissue pacing rates.Ischemia elevated the isolated cardiomyocyte's resting potential and reduced its action potential duration. In the absence of ischemic zones, the propensity in the 2D model to induce multiple re-entrant waves was low. The inclusion of ischemic zones provided the substrate for initiation of re-entrant waves leading to fibrillation. Dominant frequency, which measured the highest rate of pacing in the tissue, increased drastically with the inclusion of multiple ischemic zones. Re-entrant wave tip maximum numbers increased from 1 tip (no ischemic zone) to 34 tips when a large number (20) of ischemic zones were included. Computational limiting factors of our platform were identified using software profiling.Clinical significance. Dialysis may promote deleterious arrhythmias by increasing tissue level action potential dispersion.

Arterial Hypertension and Unusual Ascending Aortic Dilatation in a Neonate With Acute Kidney Injury: Mechanistic Computer Modeling

Background

Neonatal asphyxia caused kidney injury and severe hypertension in a newborn. An unusually dilatated ascending aorta developed. Dialysis and pharmacological treatment led to partial recovery of the ascending aortic diameters. It was hypothesized that the aortic dilatation may be associated with aortic stiffening, peripheral resistance, and cardiovascular changes. Mathematical modeling was used to better understand the potential causes of the hypertension, and to confirm our clinical treatment within the confines of the model's capabilities.

Methods

The patient's systolic arterial blood pressure showed hypertension. Echocardiographic exams showed ascending aorta dilatation during hypertension, which partially normalized upon antihypertensive treatment. To explore the underlying mechanisms of the aortic dilatation and hypertension, an existing lumped parameter hemodynamics model was deployed. Hypertension was simulated using realistic literature informed parameter values. It was also simulated using large parameter perturbations to demonstrate effects. Simulations were designed to permit examination of causal mechanisms. The hypertension inducing effects of aortic stiffnesses, vascular resistances, and cardiac hypertrophy on blood flow and pressure were simulated. Sensitivity analysis was used to stratify causes.

Results

In agreement with our clinical diagnosis, the model showed that an increase of aortic stiffness followed by augmentation of peripheral resistance are the prime causes of realistic hypertension. Increased left ventricular elastance may also cause hypertension. Ascending aortic pressure and flow increased in the simultaneous presence of left ventricle hypertrophy and augmented small vessel resistance, which indicate a plausible condition for ascending aorta dilatation. In case of realistic hypertension, sensitivity analysis showed that the treatment of both the large vessel stiffness and small vessel resistance are more important in comparison to cardiac hypertrophy.

Conclusion and Discussion

Large vessel stiffness was found to be the prime factor in arterial hypertension, which confirmed the clinical treatment. Treatment of cardiac hypertrophy appears to provide significant benefit but may be secondary to treatment of large vessel stiffness. The quantitative grading of pathophysiological mechanisms provided by the modeling may contribute to treatment recommendations. The model was limited due to a lack of data suitable to permit model identification.

3D ultrastructural organisation of calcium release units in the avian sarcoplasmic reticulum

Excitation-contraction coupling in vertebrate hearts is underpinned by calcium (Ca2+) release from Ca2+ release units (CRUs). CRUs are formed by clusters of channels called ryanodine receptors on the sarcoplasmic reticulum (SR) within the cardiomyocyte. Distances between CRUs influence the diffusion of Ca2+, thus influencing the rate and strength of excitation-contraction coupling. Avian myocytes lack T-tubules, thus Ca2+ from surface CRUs (peripheral couplings, PCs), must diffuse to internal CRU sites of the corbular SR (cSR) during centripetal propagation. Despite this, avian hearts achieve higher contractile rates and develop greater contractile strength than many mammalian hearts, which have T-tubules to provide simultaneous activation of the Ca2+ signal through the myocyte. We used 3D electron tomography to test the hypothesis that the intracellular distribution of CRUs in the avian heart permits faster and stronger contractions despite the absence T-tubules. Nearest edge-edge distances between PCs and cSR, and geometric information including surface area and volumes of individual cSR, were obtained for each cardiac chamber of the White Leghorn chicken. Computational modelling was then used to establish a relationship between CRUs distances and cell activation time in the avian heart. Our data suggest that cSR clustered close together along the Z-line is vital for rapid propagation of the Ca2+ signal from the cell periphery to the cell centre which would aid in the strong and fast contractions of the avian heart.

Computational Assessment of Blood Flow Heterogeneity in Peritoneal Dialysis Patients' Cardiac Ventricles

Dialysis prolongs life but augments cardiovascular mortality. Imaging data suggests that dialysis increases myocardial blood flow (BF) heterogeneity, but its causes remain poorly understood. A biophysical model of human coronary vasculature was used to explain the imaging observations, and highlight causes of coronary BF heterogeneity. Post-dialysis CT images from patients under control, pharmacological stress (adenosine), therapy (cooled dialysate), and adenosine and cooled dialysate conditions were obtained. The data presented disparate phenotypes. To dissect vascular mechanisms, a 3D human vasculature model based on known experimental coronary morphometry and a space filling algorithm was implemented. Steady state simulations were performed to investigate the effects of altered aortic pressure and blood vessel diameters on myocardial BF heterogeneity. Imaging showed that stress and therapy potentially increased mean and total BF, while reducing heterogeneity. BF histograms of one patient showed multi-modality. Using the model, it was found that total coronary BF increased as coronary perfusion pressure was increased. BF heterogeneity was differentially affected by large or small vessel blocking. BF heterogeneity was found to be inversely related to small blood vessel diameters. Simulation of large artery stenosis indicates that BF became heterogeneous (increase relative dispersion) and gave multi-modal histograms. The total transmural BF as well as transmural BF heterogeneity reduced due to large artery stenosis, generating large patches of very low BF regions downstream. Blocking of arteries at various orders showed that blocking larger arteries results in multi-modal BF histograms and large patches of low BF, whereas smaller artery blocking results in augmented relative dispersion and fractal dimension. Transmural heterogeneity was also affected. Finally, the effects of augmented aortic pressure in the presence of blood vessel blocking shows differential effects on BF heterogeneity as well as transmural BF. Improved aortic blood pressure may improve total BF. Stress and therapy may be effective if they dilate small vessels. A potential cause for the observed complex BF distributions (multi-modal BF histograms) may indicate existing large vessel stenosis. The intuitive BF heterogeneity methods used can be readily used in clinical studies. Further development of the model and methods will permit personalized assessment of patient BF status.

Computational assessment of the functional role of sinoatrial node exit pathways in the human heart

AIM

The human right atrium and sinoatrial node (SAN) anatomy is complex. Optical mapping experiments suggest that the SAN is functionally insulated from atrial tissue except at discrete SAN-atrial electrical junctions called SAN exit pathways, SEPs. Additionally, histological imaging suggests the presence of a secondary pacemaker close to the SAN. We hypothesise that a) an insulating border-SEP anatomical configuration is related to SAN arrhythmia; and b) a secondary pacemaker, the paranodal area, is an alternate pacemaker but accentuates tachycardia. A 3D electro-anatomical computational model was used to test these hypotheses.

METHODS

A detailed 3D human SAN electro-anatomical mathematical model was developed based on our previous anatomical reconstruction. Electrical activity was simulated using tissue specific variants of the Fenton-Karma action potential equations. Simulation experiments were designed to deploy this complex electro-anatomical system to assess the roles of border-SEPs and paranodal area by mimicking experimentally observed SAN arrhythmia. Robust and accurate numerical algorithms were implemented for solving the mono domain reaction-diffusion equation implicitly, calculating 3D filament traces, and computing dominant frequency among other quantitative measurements.

RESULTS

A centre to periphery gradient of increasing diffusion was sufficient to permit initiation of pacemaking at the centre of the 3D SAN. Re-entry within the SAN, micro re-entry, was possible by imposing significant SAN fibrosis in the presence of the insulating border. SEPs promoted the micro re-entry to generate more complex SAN-atrial tachycardia. Simulation of macro re-entry, i.e. re-entry around the SAN, was possible by inclusion of atrial fibrosis in the presence of the insulating border. The border shielded the SAN from atrial tachycardia. However, SAN micro-structure intercellular gap junctional coupling and the paranodal area contributed to prolonged atrial fibrillation. Finally, the micro-structure was found to be sufficient to explain shifts of leading pacemaker site location.

CONCLUSIONS

The simulations establish a relationship between anatomy and SAN electrical function. Microstructure, in the form of intercellular gap junction coupling, was found to regulate SAN function and arrhythmia.

3D anatomical reconstruction of human cardiac conduction system and simulation of bundle branch block after TAVI procedure

The cardiac conduction system (CCS) is responsible for the initiation and propagation of action potentials through the heart ensuring efficient pumping of blood. Understanding the anatomy of the CCS and its relationship with other major cardiac components is important to help understand arrhythmias and how certain procedures may increase the incidence of arrhythmias developing. We sectioned a whole human heart and performed Masson's trichrome histology in order to identify the components of the CCS. The histology images were used to construct a 3D anatomical model. We have shown that it is possible to create a 3D anatomical model of the human heart incorporating the CCS based on histological images, and that this model can be used to perform computer simulations of cardiac excitation. From the reconstruction we have been able to show the relative positions of the CCS components to each other. We have also shown how the close proximity of the CCS to the aortic valve can explain some of the conduction complications, such as left bundle branch block, that arise after aortic valve replacement procedures.

Ca(2+)-Clock-Dependent Pacemaking in the Sinus Node Is Impaired in Mice with a Cardiac Specific Reduction in SERCA2 Abundance

Background: The sarcoplasmic reticulum Ca²⁺-ATPase (SERCA2) pump is an important component of the Ca²⁺-clock pacemaker mechanism that provides robustness and flexibility to sinus node pacemaking. We have developed transgenic mice with reduced cardiac SERCA2 abundance (Serca2 KO) as a model for investigating SERCA2's role in sinus node pacemaking.

Methods and Results: In Serca2 KO mice, ventricular SERCA2a protein content measured by Western blotting was 75% (P < 0.05) lower than that in control mice (Serca2 FF) tissue. Immunofluorescent labeling of SERCA2a in ventricular, atrial, sinus node periphery and center tissue sections revealed 46, 45, 55, and 34% (all P < 0.05 vs. Serca2 FF) lower labeling, respectively and a mosaic pattern of expression. With telemetric ECG surveillance, we observed no difference in basal heart rate, but the PR-interval was prolonged in Serca2 KO mice: 49 ± 1 vs. 40 ± 1 ms (P < 0.001) in Serca2 FF. During exercise, heart rate in Serca2 KO mice was elevated to 667 ± 22 bpm, considerably less than 780 ± 17 bpm (P < 0.01) in Serca2 FF. In isolated sinus node preparations, 2 mM Cs⁺ caused bradycardia that was equally pronounced in Serca2 KO and Serca2 FF (32 ± 4% vs. 29 ± 5%), indicating no change in the pacemaker current, I_f. Disabling the Ca²⁺-clock with 2 μM ryanodine induced bradycardia that was less pronounced in Serca2 KO preparations (9 ± 1% vs. 20 ± 3% in Serca2 FF; P < 0.05), suggesting a disrupted Ca²⁺-clock. Mathematical modeling was used to dissect the effects of membrane- and Ca²⁺-clock components on Serca2 KO mouse heart rate and sinus node action potential. Computer modeling predicted a slowing of heart rate with SERCA2 downregulation and the heart rate slowing was pronounced at >70% reduction in SERCA2 activity.

Conclusions:Serca2 KO mice show a disrupted Ca²⁺-clock-dependent pacemaker mechanism contributing to impaired sinus node and atrioventricular node function.

Marine n-3 PUFAs modulate IKs gating, channel expression, and location in membrane microdomains

AIMS

Polyunsaturated fatty n-3 acids (PUFAs) have been reported to exhibit antiarrhythmic properties. However, the mechanisms of action remain unclear. We studied the electrophysiological effects of eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA) on IKs, and on the expression and location of Kv7.1 and KCNE1.

METHODS AND RESULTS

Experiments were performed using patch-clamp, western blot, and sucrose gradient techniques in COS7 cells transfected with Kv7.1/KCNE1 channels. Acute perfusion with both PUFAs increased Kv7.1/KCNE1 current, this effect being greater for DHA than for EPA. Similar results were found in guinea pig cardiomyocytes. Acute perfusion of either PUFA slowed the activation kinetics and EPA shifted the activation curve to the left. Conversely, chronic EPA did not modify Kv7.1/KCNE1 current magnitude and shifted the activation curve to the right. Chronic PUFAs decreased the expression of Kv7.1, but not of KCNE1, and induced spatial redistribution of Kv7.1 over the cell membrane. Cholesterol depletion with methyl-β-cyclodextrin increased Kv7.1/KCNE1 current magnitude. Under these conditions, acute EPA produced similar effects than those induced in non-cholesterol-depleted cells. A ventricular action potential computational model suggested antiarrhythmic efficacy of acute PUFA application under IKr block.

CONCLUSIONS

We provide evidence that acute application of PUFAs increases Kv7.1/KCNE1 through a probably direct effect, and shows antiarrhythmic efficacy under IKr block. Conversely, chronic EPA application modifies the channel activity through a change in the Kv7.1/KCNE1 voltage-dependence, correlated with a redistribution of Kv7.1 over the cell membrane. This loss of function may be pro-arrhythmic. This shed light on the controversial effects of PUFAs regarding arrhythmias.

Effects of human atrial ionic remodelling by β-blocker therapy on mechanisms of atrial fibrillation: a computer simulation

AIMS

Atrial anti-arrhythmic effects of β-adrenoceptor antagonists (β-blockers) may involve both a suppression of pro-arrhythmic effects of catecholamines, and an adaptational electrophysiological response to chronic β-blocker use; so-called 'pharmacological remodelling'. In human atrium, such remodelling decreases the transient outward (Ito) and inward rectifier (IK1) K(+) currents, and increases the cellular action potential duration (APD) and effective refractory period (ERP). However, the consequences of these changes on mechanisms of genesis and maintenance of atrial fibrillation (AF) are unknown. Using mathematical modelling, we tested the hypothesis that the long-term adaptational decrease in human atrial Ito and IK1 caused by chronic β-blocker therapy, i.e. independent of acute electrophysiological effects of β-blockers, in an otherwise un-remodelled atrium, could suppress AF.

METHODS AND RESULTS

Contemporarily, biophysically detailed human atrial cell and tissue models were used to investigate effects of the β-blocker-based pharmacological remodelling. Chronic β-blockade remodelling prolonged atrial cell APD and ERP. The incidence of small amplitude APD alternans in the CRN model was reduced. At the 1D tissue level, β-blocker remodelling decreased the maximum pacing rate at which APs could be conducted. At the three-dimensional organ level, β-blocker remodelling reduced the life span of re-entry scroll waves.

CONCLUSION

This study improves our understanding of the electrophysiological mechanisms of AF suppression by chronic β-blocker therapy. Atrial fibrillation suppression may involve a reduced propensity for maintenance of re-entrant excitation waves, as a consequence of increased APD and ERP.

Simulating the role of anisotropy in human atrial cardioversion

This computational study quantifies the effectiveness of feedback controlled low energy cardioversion in the anisotropic human atria. An established biophysical human cell model was adopted to reproduce Control and chronic atrial fibrillation (CAF) action potentials. The cell model was combined with a detailed human atrial geometry to construct a 3D realistic human atrial model. Scroll waves were simulated under Control and CAF conditions and the cardioversion parameters of stimulation strength and pacing duration were evaluated for scroll wave termination. Scroll waves were initiated at two locations in the atria to elicit the effects of scroll wave location. The role of anisotropy was highlighted by comparison to results from the isotropic case. Under Control conditions, scroll wave self-termination was rapid in the anisotropic case. Under CAF conditions, anisotropy caused the initiated scroll wave to degenerate into multiple scrolls with each evolving erratically or pinning to anatomical defects. The cardioversion successfully terminated scroll waves within 10 s, but the stimulus strength had a strong correlation to the location of the scroll wave. The low energy stimulation strength was always lower than the threshold stimulus. Anisotropy plays an important role in atrial electrical properties. Anisotropy aggravates CAF and leads to high frequency atrial pacing. The efficacy of cardioversion is significantly affected by anisotropy.