Macrophages Facilitate Electrical Conduction in the Heart.

Organ-specific functions of tissue-resident macrophages in the steady-state heart are unknown. Here, we show that cardiac macrophages facilitate electrical conduction through the distal atrioventricular node, where conducting cells densely intersperse with elongated macrophages expressing connexin 43. When coupled to spontaneously beating cardiomyocytes via connexin-43-containing gap junctions, cardiac macrophages have a negative resting membrane potential and depolarize in synchrony with cardiomyocytes. Conversely, macrophages render the resting membrane potential of cardiomyocytes more positive and, according to computational modeling, accelerate their repolarization. Photostimulation of channelrhodopsin-2-expressing macrophages improves atrioventricular conduction, whereas conditional deletion of connexin 43 in macrophages and congenital lack of macrophages delay atrioventricular conduction. In the Cd11bDTR mouse, macrophage ablation induces progressive atrioventricular block. These observations implicate macrophages in normal and aberrant cardiac conduction.

Mechano-electrical interactions and heterogeneities in wild-type and drug-induced long QT syndrome rabbits.

Electromechanical reciprocity - comprising electro-mechanical (EMC) and mechano-electric coupling (MEC) - provides cardiac adaptation to changing physiological demands. Understanding electromechanical reciprocity and its impact on function and heterogeneity in pathological conditions - such as (drug-induced) acquired long QT syndrome (aLQTS) - might lead to novel insights in arrhythmogenesis. Our aim is to investigate how electrical changes impact on mechanical function (EMC) and vice versa (MEC) under physiological conditions and in aLQTS. To measure regional differences in EMC and MEC in vivo, we used tissue phase mapping cardiac MRI and a 24-lead ECG vest in healthy (control) and IKr -blocker E-4031-induced aLQTS rabbit hearts. MEC was studied in vivo by acutely increasing cardiac preload, and ex vivo by using voltage optical mapping (OM) in beating hearts at different preloads. In aLQTS, electrical repolarization (heart rate corrected RT-interval, RTn370) was prolonged compared to control (P < 0.0001) with increased spatial and temporal RT heterogeneity (P < 0.01). Changing electrical function (in aLQTS) resulted in significantly reduced diastolic mechanical function and prolonged contraction duration (EMC), causing increased apico-basal mechanical heterogeneity. Increased preload acutely prolonged RTn370 in both control and aLQTS hearts (MEC). This effect was more pronounced in aLQTS (P < 0.0001). Additionally, regional RT-dispersion increased in aLQTS. Motion-correction allowed us to determine APD-prolongation in beating aLQTS hearts, but limited motion correction accuracy upon preload-changes prevented a clear analysis of MEC ex vivo. Mechano-induced RT-prolongation and increased heterogeneity were more pronounced in aLQTS than in healthy hearts. Acute MEC effects may play an additional role in LQT-related arrhythmogenesis, warranting further mechanistic investigations. KEY POINTS: Electromechanical reciprocity comprising excitation-contraction coupling (EMC) and mechano-electric feedback loops (MEC) is essential for physiological cardiac function. Alterations in electrical and/or mechanical heterogeneity are known to have potentially pro-arrhythmic effects. In this study, we aimed to investigate how electrical changes impact on the mechanical function (EMC) and vice versa (MEC) both under physiological conditions (control) and in acquired long QT syndrome (aLQTS). We show that changing the electrical function (in aLQTS) results in significantly altered mechanical heterogeneity via EMC and, vice versa, that increasing the preload acutely prolongs repolarization duration and increases electrical heterogeneity, particularly in aLQTS as compared to control. Our results substantiate the hypothesis that LQTS is an ?electro-mechanical', rather than a 'purely electrical', disease and suggest that acute MEC effects may play an additional role in LQT-related arrhythmogenesis.

Hyperinsulinemic Hypoglycemia Associated with a CaV1.2 Variant with Mixed Gain- and Loss-of-Function Effects.

The voltage-dependent L-type calcium channel isoform CaV1.2 is critically involved in many physiological processes, e.g., in cardiac action potential formation, electromechanical coupling and regulation of insulin secretion by beta cells. Gain-of-function mutations in the calcium voltage-gated channel subunit alpha 1 C (CACNA1C) gene, encoding the CaV1.2 α1-subunit, cause Timothy syndrome (TS), a multisystemic disorder that includes autism spectrum disorders and long QT (LQT) syndrome. Strikingly, TS patients frequently suffer from hypoglycemia of yet unproven origin. Using next-generation sequencing, we identified a novel heterozygous CACNA1C mutation in a patient with congenital hyperinsulinism (CHI) and associated hypoglycemic episodes. We characterized the electrophysiological phenotype of the mutated channel using voltage-clamp recordings and in silico action potential modeling experiments. The identified CaV1.2L566P mutation causes a mixed electrophysiological phenotype of gain- and loss-of-function effects. In silico action potential modeling supports that this mixed electrophysiological phenotype leads to a tissue-specific impact on beta cells compared to cardiomyocytes. Thus, CACNA1C variants may be associated with non-syndromic hyperinsulinemic hypoglycemia without long-QT syndrome, explained by very specific electrophysiological properties of the mutated channel. We discuss different biochemical characteristics and clinical impacts of hypoglycemia in the context of CACNA1C variants and show that these may be associated with significant morbidity for Timothy Syndrome patients. Our findings underline that the potential of hypoglycemia warrants careful attention in patients with CACNA1C variants, and such variants should be included in the differential diagnosis of non-syndromic congenital hyperinsulinism.

Cycle length statistics during human atrial fibrillation reveal refractory properties of the underlying substrate: a combined in silico and clinical test of concept study.

AIMS: The treatment of atrial fibrillation beyond pulmonary vein isolation has remained an unsolved challenge. Targeting regions identified by different substrate mapping approaches for ablation resulted in ambiguous outcomes. With the effective refractory period being a fundamental prerequisite for the maintenance of fibrillatory conduction, this study aims at estimating the effective refractory period with clinically available measurements. METHODS AND RESULTS: A set of 240 simulations in a spherical model of the left atrium with varying model initialization, combination of cellular refractory properties, and size of a region of lowered effective refractory period was implemented to analyse the capabilities and limitations of cycle length mapping. The minimum observed cycle length and the 25% quantile were compared to the underlying effective refractory period. The density of phase singularities was used as a measure for the complexity of the excitation pattern. Finally, we employed the method in a clinical test of concept including five patients. Areas of lowered effective refractory period could be distinguished from their surroundings in simulated scenarios with successfully induced multi-wavelet re-entry. Larger areas and higher gradients in effective refractory period as well as complex activation patterns favour the method. The 25% quantile of cycle lengths in patients with persistent atrial fibrillation was found to range from 85 to 190 ms. CONCLUSION: Cycle length mapping is capable of highlighting regions of pathologic refractory properties. In combination with complementary substrate mapping approaches, the method fosters confidence to enhance the treatment of atrial fibrillation beyond pulmonary vein isolation particularly in patients with complex activation patterns.

Molecular Mechanism of Autosomal Recessive Long QT-Syndrome 1 without Deafness.

KCNQ1 encodes the voltage-gated potassium (Kv) channel KCNQ1, also known as KvLQT1 or Kv7.1. Together with its ß-subunit KCNE1, also denoted as minK, this channel generates the slowly activating cardiac delayed rectifier current IKs, which is a key regulator of the heart rate dependent adaptation of the cardiac action potential duration (APD). Loss-of-function mutations in KCNQ1 cause congenital long QT1 (LQT1) syndrome, characterized by a delayed cardiac repolarization and a prolonged QT interval in the surface electrocardiogram. Autosomal dominant loss-of-function mutations in KCNQ1 result in long QT syndrome, called Romano-Ward Syndrome (RWS), while autosomal recessive mutations lead to Jervell and Lange-Nielsen syndrome (JLNS), associated with deafness. Here, we identified a homozygous KCNQ1 mutation, c.1892_1893insC (p.P631fs*20), in a patient with an isolated LQT syndrome (LQTS) without hearing loss. Nevertheless, the inheritance trait is autosomal recessive, with heterozygous family members being asymptomatic. The results of the electrophysiological characterization of the mutant, using voltage-clamp recordings in Xenopus laevis oocytes, are in agreement with an autosomal recessive disorder, since the IKs reduction was only observed in homomeric mutants, but not in heteromeric IKs channel complexes containing wild-type channel subunits. We found that KCNE1 rescues the KCNQ1 loss-of-function in mutant IKs channel complexes when they contain wild-type KCNQ1 subunits, as found in the heterozygous state. Action potential modellings confirmed that the recessive c.1892_1893insC LQT1 mutation only affects the APD of homozygous mutation carriers. Thus, our study provides the molecular mechanism for an atypical autosomal recessive LQT trait that lacks hearing impairment.

Modeling effects of voltage dependent properties of the cardiac muscarinic receptor on human sinus node function.

The cardiac muscarinic receptor (M2R) regulates heart rate, in part, by modulating the acetylcholine (ACh) activated K+ current IK,ACh through dissociation of G-proteins, that in turn activate KACh channels. Recently, M2Rs were noted to exhibit intrinsic voltage sensitivity, i.e. their affinity for ligands varies in a voltage dependent manner. The voltage sensitivity of M2R implies that the affinity for ACh (and thus the ACh effect) varies throughout the time course of a cardiac electrical cycle. The aim of this study was to investigate the contribution of M2R voltage sensitivity to the rate and shape of the human sinus node action potentials in physiological and pathophysiological conditions. We developed a Markovian model of the IK,ACh modulation by voltage and integrated it into a computational model of human sinus node. We performed simulations with the integrated model varying ACh concentration and voltage sensitivity. Low ACh exerted a larger effect on IK,ACh at hyperpolarized versus depolarized membrane voltages. This led to a slowing of the pacemaker rate due to an attenuated slope of phase 4 depolarization with only marginal effect on action potential duration and amplitude. We also simulated the theoretical effects of genetic variants that alter the voltage sensitivity of M2R. Modest negative shifts in voltage sensitivity, predicted to increase the affinity of the receptor for ACh, slowed the rate of phase 4 depolarization and slowed heart rate, while modest positive shifts increased heart rate. These simulations support our hypothesis that altered M2R voltage sensitivity contributes to disease and provide a novel mechanistic foundation to study clinical disorders such as atrial fibrillation and inappropriate sinus tachycardia.

Postpartum hormones oxytocin and prolactin cause pro-arrhythmic prolongation of cardiac repolarization in long QT syndrome type 2.

AIMS: Women with long QT syndrome 2 (LQT2) have a particularly high postpartal risk for lethal arrhythmias. We aimed at investigating whether oxytocin and prolactin contribute to this risk by affecting repolarization. METHODS AND RESULTS: In female transgenic LQT2 rabbits (HERG-G628S, loss of IKr), hormone effects on QT/action potential duration (APD) were assessed (0.2-200 ng/L). Hormone effects (200 ng/L) on ion currents and cellular APD were determined in transfected cells and LQT2 cardiomyocytes. Hormone effects on ion channels were assessed with qPCR and western blot. Experimental data were incorporated into in silico models to determine the pro-arrhythmic potential. Oxytocin prolonged QTc and steepened QT/RR-slope in vivo and prolonged ex vivo APD75 in LQT2 hearts. Prolactin prolonged APD75 at high concentrations. As underlying mechanisms, we identified an oxytocin- and prolactin-induced acute reduction of IKs-tail and IKs-steady (-25.5%, oxytocin; -13.3%, prolactin, P < 0.05) in CHO-cells and LQT2-cardiomyocytes. IKr currents were not altered. This oxytocin-/prolactin-induced IKs reduction caused APD90 prolongation (+11.9%/+13%, P < 0.05) in the context of reduced/absent IKr in LQT2 cardiomyocytes. Hormones had no effect on IK1 and ICa,L in cardiomyocytes. Protein and mRNA levels of CACNA1C/Cav1.2 and RyR2 were enhanced by oxytocin and prolactin. Incorporating these hormone effects into computational models resulted in reduced repolarization reserve and increased propensity to pro-arrhythmic permanent depolarization, lack of capture and early afterdepolarizations formation. CONCLUSIONS: Postpartum hormones oxytocin and prolactin prolong QT/APD in LQT2 by reducing IKs and by increasing Cav1.2 and RyR2 expression/transcription, thereby contributing to the increased postpartal arrhythmic risk in LQT2.

Influence of left atrial size on P-wave morphology: differential effects of dilation and hypertrophy.

AIMS: Chronic left atrial enlargement (LAE) increases the risk of atrial fibrillation. Electrocardiogram (ECG) criteria might provide a means to diagnose LAE and identify patients at risk; however, current criteria perform poorly. We seek to characterize the potentially differential effects of atrial dilation vs. hypertrophy on the ECG P-wave. METHODS AND RESULTS: We predict effects on the P-wave of (i) left atrial dilation (LAD), i.e. an increase of LA cavity volume without an increase in myocardial volume, (ii) left atrial concentric hypertrophy (LACH), i.e. a thickened myocardial wall, and (iii) a combination of the two. We performed a computational study in a cohort of 72 anatomical variants, derived from four human atrial anatomies. To model LAD, pressure was applied to the LA endocardium increasing cavity volume by up to 100%. For LACH, the LA wall was thickened by up to 3.3 mm. P-waves were derived by simulating atrial excitation propagation and computing the body surface ECG. The sensitivity regarding changes beyond purely anatomical effects was analysed by altering conduction velocity by 25% in 96 additional model variants. Left atrial dilation prolonged P-wave duration (PWd) in two of four subjects; in one subject a shortening, and in the other a variable change were seen. Left atrial concentric hypertrophy, in contrast, consistently increased P-wave terminal force in lead V1 (PTF-V1) in all subjects through an enlarged amplitude while PWd was unaffected. Combined hypertrophy and dilation generally enhanced the effect of hypertrophy on PTF-V1. CONCLUSION: Isolated LAD has moderate effects on the currently used P-wave criteria, explaining the limited utility of PWd and PTF-V1 in detecting LAE in clinical practice. In contrast, PTF-V1 may be a more sensitive indicator of LA myocardial hypertrophy.

Anatomical and spiral wave reentry in a simplified model for atrial electrophysiology.

For modeling the propagation of action potentials in the human atria, various models have been developed in the past, which take into account in detail the influence of the numerous ionic currents flowing through the cell membrane. Aiming at a simplified description, the Bueno-Orovio-Cherry-Fenton (BOCF) model for electric wave propagation in the ventricle has been adapted recently to atrial physiology. Here, we study this adapted BOCF (aBOCF) model with respect to its capability to accurately generate spatio-temporal excitation patterns found in anatomical and spiral wave reentry. To this end, we compare results of the aBOCF model with the more detailed one proposed by Courtemanche, Ramirez and Nattel (CRN model). We find that characteristic features of the reentrant excitation patterns seen in the CRN model are well captured by the aBOCF model. This opens the possibility to study origins of atrial fibrillation based on a simplified but still reliable description.

Electrophysiological characterization of a large set of novel variants in the SCN5A-gene: identification of novel LQTS3 and BrS mutations.

SCN5A encodes for the α-subunit of the cardiac voltage-gated sodium channel Nav1.5. Gain-of-function mutations in SCN5A are related to congenital long QT syndrome (LQTS3) characterized by delayed cardiac repolarization, leading to a prolonged QT interval in the ECG. Loss-of-function mutations in SCN5A are related to Brugada syndrome (BrS), characterized by an ST-segment elevation in the right precordial leads (V1-V3). The aim of this study was the characterization of a large set of novel SCN5A variants found in patients with different cardiac phenotypes, mainly LQTS and BrS. SCN5A variants of 13 families were functionally characterized in Xenopus laevis oocytes using the two-electrode voltage-clamp technique. We found in most of the cases, but not all, that the electrophysiology of the variants correlated with the clinically diagnosed phenotype. A susceptibility to develop LQTS can be suggested in patients carrying the variants S216L, K480N, A572D, F816Y, and G983D. However, taking the phenotype into account, the presence of the variants in genomic data bases, the mutational segregation, combined with our in vitro and in silico experiments, the variants S216L, S262G, K480N, A572D, F816Y, G983D, and T1526P remain as variants of unknown significance. However, the SCN5A variants R568H and A993T can be classified as pathogenic LQTS3 causing mutations, while R222stop and R2012H are novel BrS causing mutations.

Influence of the earliest right atrial activation site and its proximity to interatrial connections on P-wave morphology.

AIMS: P-wave morphology correlates with the risk for atrial fibrillation (AF). Left atrial (LA) enlargement could explain both the higher risk for AF and higher P-wave terminal force (PTF) in lead V1. However, PTF-V1 has been shown to correlate poorly with LA size. We hypothesize that PTF-V1 is also affected by the earliest activated site (EAS) in the right atrium and its proximity to inter-atrial connections (IAC), which both show tremendous variability. METHODS AND RESULTS: Atrial excitation was triggered from seven different EAS in a cohort of eight anatomically personalized computational models. The posterior IACs were non-conductive in a second set of simulations. Body surface ECGs were computed and separated by left and right atrial contributions. Mid-septal EAS yielded the highest PTF-V1. More anterior/superior and more inferior EAS yielded lower absolute PTF-V1 values deviating by a factor of up to 2.0 for adjacent EAS. Earliest right-to-left activation was conducted via Bachmann's Bundle (BB) for anterior/superior EAS and shifted towards posterior IACs for more inferior EAS. Non-conducting posterior IACs increased PTF-V1 by up to 150% compared to intact posterior IACs for inferior EAS. LA contribution to the P-wave integral was 24% on average. CONCLUSION: The electrical contributor's site of earliest activation and intactness of posterior IACs affect PTF-V1 significantly by changing LA breakthrough sites independent from LA size. This should be considered for interpretation of electrocardiographical signs of LA abnormality and LA enlargement.

In silico screening of the key cellular remodeling targets in chronic atrial fibrillation.

Chronic atrial fibrillation (AF) is a complex disease with underlying changes in electrophysiology, calcium signaling and the structure of atrial myocytes. How these individual remodeling targets and their emergent interactions contribute to cell physiology in chronic AF is not well understood. To approach this problem, we performed in silico experiments in a computational model of the human atrial myocyte. The remodeled function of cellular components was based on a broad literature review of in vitro findings in chronic AF, and these were integrated into the model to define a cohort of virtual cells. Simulation results indicate that while the altered function of calcium and potassium ion channels alone causes a pronounced decrease in action potential duration, remodeling of intracellular calcium handling also has a substantial impact on the chronic AF phenotype. We additionally found that the reduction in amplitude of the calcium transient in chronic AF as compared to normal sinus rhythm is primarily due to the remodeling of calcium channel function, calcium handling and cellular geometry. Finally, we found that decreased electrical resistance of the membrane together with remodeled calcium handling synergistically decreased cellular excitability and the subsequent inducibility of repolarization abnormalities in the human atrial myocyte in chronic AF. We conclude that the presented results highlight the complexity of both intrinsic cellular interactions and emergent properties of human atrial myocytes in chronic AF. Therefore, reversing remodeling for a single remodeled component does little to restore the normal sinus rhythm phenotype. These findings may have important implications for developing novel therapeutic approaches for chronic AF.

Arrhythmic potency of human ether-a-go-go-related gene mutations L532P and N588K in a computational model of human atrial myocytes.

AIMS: Human ether-à-go-go-related gene (hERG) missense mutations N588K and L532P are both associated with atrial fibrillation (AF). However, the underlying gain-of-function mechanism is different. The aim of this computational study is to assess and understand the arrhythmogenic mechanisms of these genetic disorders on the cellular and tissue level as a basis for the improvement of therapeutic strategies. METHODS AND RESULTS: The IKr formulation of an established model of human atrial myocytes was adapted by using the measurement data of wild-type and mutant hERG channels. Restitution curves of the action potential duration and its slope, effective refractory period (ERP), conduction velocity, reentry wavelength (WL), and the vulnerable window (VW) were determined in a one-dimensional (1D) tissue strand. Moreover, spiral wave inducibility and rotor lifetime in a 2D tissue patch were evaluated. The two mutations caused an increase in IKr regarding both peak amplitude and current integral, whereas the duration during which IKr is active was decreased. The WL was reduced due to a shorter ERP. Spiral waves could be initiated by using mutation models as opposed to the control case. The frequency dependency of the VW was reversed. CONCLUSION: Both mutations showed an increased arrhythmogenicity due to decreased refractory time in combination with a more linear repolarization phase. The effects were more pronounced for mutation L532P than for N588K. Furthermore, spiral waves presented higher stability and a more regular pattern for L532P. These in silico investigations unveiling differences of mutations affecting the same ion channel may help to advance genotype-guided AF prevention and therapy strategies.

In-silico assessment of the dynamic effects of amiodarone and dronedarone on human atrial patho-electrophysiology.

AIMS: The clinical efficacy in preventing the recurrence of atrial fibrillation (AF) is higher for amiodarone than for dronedarone. Moreover, pharmacotherapy with these drugs is less successful in patients with remodelled substrate induced by chronic AF (cAF) and patients suffering from familial AF. To date, the reasons for these phenomena are only incompletely understood. We analyse the effects of the drugs in a computational model of atrial electrophysiology. METHODS AND RESULTS: The Courtemanche-Ramirez-Nattel model was adapted to represent cAF remodelled tissue and hERG mutations N588K and L532P. The pharmacodynamics of amiodarone and dronedarone were investigated with respect to their dose and heart rate dependence by evaluating 10 descriptors of action potential morphology and conduction properties. An arrhythmia score was computed based on a subset of these biomarkers and analysed regarding circadian variation of drug concentration and heart rate. Action potential alternans at high frequencies was observed over the whole dronedarone concentration range at high frequencies, while amiodarone caused alternans only in a narrow range. The total score of dronedarone reached critical values in most of the investigated dynamic scenarios, while amiodarone caused only minor score oscillations. Compared with the other substrates, cAF showed significantly different characteristics resulting in a lower amiodarone but higher dronedarone concentration yielding the lowest score. CONCLUSION: Significant differences exist in the frequency and concentration-dependent effects between amiodarone and dronedarone and between different atrial substrates. Our results provide possible explanations for the superior efficacy of amiodarone and may aid in the design of substrate-specific pharmacotherapy for AF.

Patient-specific modeling of atrial fibrosis increases the accuracy of sinus rhythm simulations and may explain maintenance of atrial fibrillation.

Left atrial fibrosis is thought to contribute to the manifestation of atrial fibrillation (AF). Late Gadolinium enhancement (LGE) MRI has the potential to image regions of low perfusion, which can be related to fibrosis. We show that a simulation with a patient-specific model including left atrial regional fibrosis derived from LGE-MRI reproduces local activation in the left atrium more precisely than the regular simulation without fibrosis. AF simulations showed a spontaneous termination of the arrhythmia in the absence of fibrosis and a stable rotor center in the presence of fibrosis. The methodology may provide a tool for a deeper understanding of the mechanisms maintaining AF and eventually also for the planning of substrate-guided ablation procedures in the future.

Whole-heart computational modelling provides further mechanistic insights into ST-elevation in Brugada syndrome.

Background: Brugada syndrome (BrS) is characterized by dynamic ST-elevations in right precordial leads and increased risk of ventricular fibrillation and sudden cardiac death. As the mechanism underlying ST-elevation and malignant arrhythmias is controversial computational modeling can aid in exploring the disease mechanism. Thus we aim to test the main competing hypotheses ('delayed depolarization' vs. 'early repolarization') of BrS in a whole-heart computational model. Methods: In a 3D whole-heart computational model, delayed epicardial RVOT activation with local conduction delay was simulated by reducing conductivity in the epicardial RVOT. Early repolarization was simulated by instead increasing the transient outward potassium current (Ito) in the same region. Additionally, a reduction in the fast sodium current (INa) was incorporated in both models. Results: Delayed depolarization with local conduction delay in the computational model resulted in coved-type ST-elevation with negative T-waves in the precordial surface ECG leads. 'Saddleback'-shaped ST-elevation was obtained with reduced substrate extent or thickness. Increased Ito simulations showed early repolarization in the RVOT with a descending but not coved-type ST-elevation. Reduced INa did not show a significant effect on ECG morphology. Conclusions: In this whole-heart BrS computational model of both major hypotheses, realistic coved-type ECG resulted only from delayed epicardial RVOT depolarization with local conduction delay but not early repolarizing ion channel modifications. These simulations provide further support for the depolarization hypothesis as electrophysiological mechanism underlying BrS.

Impact of amiodarone and cisapride on simulated human ventricular electrophysiology and electrocardiograms.

AIMS: Amiodarone and cisapride are both known to prolong the QT interval, yet the two drugs have different effects on arrhythmia. Cisapride can cause torsades de pointes while amiodarone is found to be anti-arrhythmic. A computational model was used to investigate the action of these two drugs. METHODS AND RESULTS: In a biophysically detailed model, the ion current conductivities affected by both drugs were reduced in order to simulate the pharmacological effects in healthy and ischaemic cells. Furthermore, restitution curves of the action potential duration (APD), effective refractory period, conduction velocity, wavelength, and the vulnerable window were determined in a one-dimensional (1D) tissue strand. Moreover, cardiac excitation propagation was computed in a 3D model of healthy ventricles. The corresponding body surface potentials were calculated and standard 12-lead electrocardiograms were derived. Both cisapride and amiodarone caused a prolongation of the QT interval and the refractory period. However, cisapride did not significantly alter the conduction-related properties, such as e.g. the wavelength or vulnerable window, whereas amiodarone had a larger impact on them. It slightly flattened the APD restitution slope and furthermore reduced the conduction velocity and wavelength. CONCLUSION: Both drugs show similar prolongation of the QT interval, although they present different electrophysiological properties in the single-cell as well as in tissue simulations of cardiac excitation propagation. These computer simulations help to better understand the underlying mechanisms responsible for the initiation or termination of arrhythmias caused by amiodarone and cisapride.

A computational model of rabbit geometry and ECG: Optimizing ventricular activation sequence and APD distribution.

Computational modeling of electrophysiological properties of the rabbit heart is a commonly used way to enhance and/or complement findings from classic lab work on single cell or tissue levels. Yet, thus far, there was no possibility to extend the scope to include the resulting body surface potentials as a way of validation or to investigate the effect of certain pathologies. Based on CT imaging, we developed the first openly available computational geometrical model not only of the whole heart but also the complete torso of the rabbit. Additionally, we fabricated a 32-lead ECG-vest to record body surface potential signals of the aforementioned rabbit. Based on the developed geometrical model and the measured signals, we then optimized the activation sequence of the ventricles, recreating the functionality of the Purkinje network, and we investigated different apico-basal and transmural gradients in action potential duration. Optimization of the activation sequence resulted in an average root mean square error between measured and simulated signal of 0.074 mV/ms for all leads. The best-fit T-Wave, compared to measured data (0.038 mV/ms), resulted from incorporating an action potential duration gradient from base to apex with a respective shortening of 20 ms and a transmural gradient with a shortening of 15 ms from endocardium to epicardium. By making our model and measured data openly available, we hope to give other researchers the opportunity to verify their research, as well as to create the possibility to investigate the impact of electrophysiological alterations on body surface signals for translational research.

Novel insights into the electrophysiology of murine cardiac macrophages: relevance of voltage-gated potassium channels.

AIMS: Macrophages (MΦ), known for immunological roles, such as phagocytosis and antigen presentation, have been found to electrotonically couple to cardiomyocytes (CM) of the atrioventricular node via Cx43, affecting cardiac conduction in isolated mouse hearts. Here, we characterize passive and active electrophysiological properties of murine cardiac resident MΦ, and model their potential electrophysiological relevance for CM. METHODS AND RESULTS: We combined classic electrophysiological approaches with 3D florescence imaging, RNA-sequencing, pharmacological interventions, and computer simulations. We used Cx3cr1eYFP/+ mice wherein cardiac MΦ are fluorescently labelled. FACS-purified fluorescent MΦ from mouse hearts were studied by whole-cell patch-clamp. MΦ electrophysiological properties include: membrane resistance 2.2±0.1 GΩ (all data mean±SEM), capacitance 18.3±0.1 pF, resting membrane potential -39.6±0.3 mV, and several voltage-activated, outward or inwardly rectifying potassium currents. Using ion channel blockers (barium, TEA, 4-AP, margatoxin, XEN-D0103, and DIDS), flow cytometry, immuno-staining, and RNA-sequencing, we identified Kv1.3, Kv1.5, and Kir2.1 as channels contributing to observed ion currents. MΦ displayed four patterns for outward and two for inward-rectifier potassium currents. Additionally, MΦ showed surface expression of Cx43, a prerequisite for homo- and/or heterotypic electrotonic coupling. Experimental results fed into development of an original computational model to describe cardiac MΦ electrophysiology. Computer simulations to quantitatively assess plausible effects of MΦ on electrotonically coupled CM showed that MΦ can depolarize resting CM, shorten early and prolong late action potential duration, with effects depending on coupling strength and individual MΦ electrophysiological properties, in particular resting membrane potential and presence/absence of Kir2.1. CONCLUSION: Our results provide a first electrophysiological characterization of cardiac resident MΦ, and a computational model to quantitatively explore their relevance in the heterocellular heart. Future work will be focussed at distinguishing electrophysiological effects of MΦ-CM coupling on both cell types during steady-state and in patho-physiological remodelling, when immune cells change their phenotype, proliferate, and/or invade from external sources.