Effects of fibroblast-myocyte coupling on cardiac conduction and vulnerability to reentry: A computational study

The effect of fixed and functional remodelling on conduction velocity, wavefront propagation, and rotational activity formation in atrial fibrillation

AIMS

Pathophysiology of atrial fibrillation (AF) remains unclear. Interactions between scar and conduction velocity (CV) and their impact on wavefront propagation in sinus rhythm (SR) and rotational activity burden in AF were evaluated.

METHODS AND RESULTS

Local activation times (LATs) and voltage data were obtained from patients undergoing ablation for persistent AF. Omnipolar voltage (OV) and bipolar voltage (BV) data were obtained during AF and SR at pacing intervals of 600 and 250 ms. Local activation times were used to determine CV dynamics and their relationship to the underlying voltage and pivot points in SR. Computational modelling studies were performed to evaluate the impact of CVs and fibrosis on rotational activity burden in AF. Data from 60 patients with a total of 2 768 400 LAT and voltage points were analysed (46 140 ± 5689 points/patient). Voltage determined CV dynamics. Enhanced CV heterogeneity sites were predominantly mapped to low-voltage zones (LVZs) (0.2-0.49 mV) (128/168, 76.2%) rather than LVZs (<0.2 mV) and frequently co-located to pivot points (151/168, 89.9%). Atrial fibrillation OV maps correlated better with SR BV 250 ms than 600 ms maps, thereby representing fixed and functional remodelling. Sinus rhythm maps at 250 ms compared with 600 ms harboured a greater number of pivot points. Increased CV slowing and functional remodelling on computational models resulted in a greater rotational activity burden.

CONCLUSION

Conduction velocity dynamics are impacted by the degree of scar. Conduction velocity heterogeneity and functional remodelling impacts wavefront propagation in SR and rotational activity burden in AF. This study provides insight into the pathophysiology of AF and identifies potential novel ablation targets.

Slow left atrial conduction velocity in the anterior wall calculated by electroanatomic mapping predicts atrial fibrillation recurrence after catheter ablation-Systematic review and meta-analysis

Background

This study aimed to investigate and perform diagnostic test meta-analysis on whether slow left atrial conduction velocity (LACV) in the anterior wall calculated by electroanatomic mapping predicts atrial fibrillation (AF) recurrence after catheter ablation.

Methods

Extensive literature search was performed on PubMed, SCOPUS, and EuropePMC up to June 5, 2024. The exposure group included AF patients with slow LACV in the anterior wall, while the control group included AF patients without slow LACV in the anterior wall. Slow LACV in the anterior wall was defined as LACV below study-specific cut-off points in m/s, measured by invasive electroanatomic mapping. The primary outcome of this study was AF recurrence, defined as AF/Atrial Flutter/Atrial Tachyarrhythmias lasting over 30 s at least 3 months after the blanking period postablation.

Results

This systematic review and meta-analysis included seven studies, involving a sample size of 1428 patients with mean follow-up duration were 13 months. Patients with AF recurrence has slower LACV in the anterior wall (mean difference - 0.16 m/s [-0.18, -0.15], p < .001). Slow LACV in the anterior wall defined as LACV below 0.70-0.88 m/s was associated with increased AF (adjusted OR 3.41 [1.55, 7.50], p = .002). Slow LACV in the anterior wall has an AUROC of 0.80 [0.76-0.83], sensitivity of 70% [52, 84], specificity of 76% [67, 83], positive likelihood ratio of 2.9 [2.3, 3.6], negative likelihood ratio of 0.39 [0.25, 0.63] for predicting AF recurrence postablation.

Conclusion

Slow LACV in the anterior wall was associated with AF recurrence after catheter ablation.

Autonomic modulation impacts conduction velocity dynamics and wavefront propagation in the left atrium

AIMS

Atrial fibrosis and autonomic remodelling are proposed pathophysiological mechanisms in atrial fibrillation (AF). Their impact on conduction velocity (CV) dynamics and wavefront propagation was evaluated.

METHODS AND RESULTS

Local activation times (LATs), voltage, and geometry data were obtained from patients undergoing ablation for persistent AF. LATs were obtained at three pacing intervals (PIs) in sinus rhythm (SR). LATs were used to determine CV dynamics and their relationship to local voltage amplitude. The impact of autonomic modulation- pharmacologically and with ganglionated plexi (GP) stimulation, on CV dynamics, wavefront propagation, and pivot points (change in wavefront propagation of ≥90°) was determined in SR. Fifty-four patients were included. Voltage impacted CV dynamics whereby at non-low voltage zones (LVZs) (≥0.5 mV) the CV restitution curves are steeper [0.03 ± 0.03 m/s ΔCV PI 600-400 ms (PI1), 0.54 ± 0.09 m/s ΔCV PI 400-250 ms (PI2)], broader at LVZ (0.2-0.49 mV) (0.17 ± 0.09 m/s ΔCV PI1, 0.25 ± 0.11 m/s ΔCV PI2), and flat at very LVZ (<0.2 mV) (0.03 ± 0.01 m/s ΔCV PI1, 0.04 ± 0.02 m/s ΔCV PI2). Atropine did not change CV dynamics, while isoprenaline and GP stimulation resulted in greater CV slowing with rate. Isoprenaline (2.7 ± 1.1 increase/patient) and GP stimulation (2.8 ± 1.3 increase/patient) promoted CV heterogeneity, i.e. rate-dependent CV (RDCV) slowing sites. Most pivot points co-located to RDCV slowing sites (80.2%). Isoprenaline (1.3 ± 1.1 pivot increase/patient) and GP stimulation (1.5 ± 1.1 increase/patient) also enhanced the number of pivot points identified.

CONCLUSION

Atrial CV dynamics is affected by fibrosis burden and influenced by autonomic modulation which enhances CV heterogeneity and distribution of pivot points. This study provides further insight into the impact of autonomic remodelling in AF.

Sinoatrial node heterogeneity and fibroblasts increase atrial driving capability in a two-dimensional human computational model

Background: Cardiac pacemaking remains an unsolved matter from many perspectives. Extensive experimental and computational studies have been performed to describe the sinoatrial physiology across different scales, from the molecular to clinical levels. Nevertheless, the mechanism by which a heartbeat is generated inside the sinoatrial node and propagated to the working myocardium is not fully understood at present. This work aims to provide quantitative information about this fascinating phenomenon, especially regarding the contributions of cellular heterogeneity and fibroblasts to sinoatrial node automaticity and atrial driving. Methods: We developed a bidimensional computational model of the human right atrial tissue, including the sinoatrial node. State-of-the-art knowledge of the anatomical and physiological aspects was adopted during the design of the baseline tissue model. The novelty of this study is the consideration of cellular heterogeneity and fibroblasts inside the sinoatrial node for investigating the manner by which they tune the robustness of stimulus formation and conduction under different conditions (baseline, ionic current blocks, autonomic modulation, and external high-frequency pacing). Results: The simulations show that both heterogeneity and fibroblasts significantly increase the safety factor for conduction by more than 10% in almost all the conditions tested and shorten the sinus node recovery time after overdrive suppression by up to 60%. In the human model, especially under challenging conditions, the fibroblasts help the heterogeneous myocytes to synchronise their rate (e.g. -82% inσ C L under 25 nM of acetylcholine administration) and capture the atrium (with 25% L-type calcium current block). However, the anatomical and gap junctional coupling aspects remain the most important model parameters that allow effective atrial excitations. Conclusion: Despite the limitations to the proposed model, this work suggests a quantitative explanation to the astonishing overall heterogeneity shown by the sinoatrial node.

Slow conduction velocity predicts atrial fibrillation recurrence after radiofrequency ablation

OBJECTIVE

To evaluate the progression of electrophysiological phenomena in atrial fibrillation (AF) and elucidate the association between the left atrial conduction velocity (LACV) and AF recurrence after pulmonary vein isolation.

METHODS

A total of 188 AF patients (121 paroxysmal AF and 67 persistent AF) who underwent PVI for the first time were enrolled in this prospective study. The left atrium was mapped using a 20-pole electrode catheter combined with the CARTO3 system. The conduction distances and conduction times of the left atrium from the Bachmann bundle area to the mitral isthmus were calculated. Anterior, posterior, and septal LACV were calculated as conduction distance divided by conduction time.

RESULTS

The anterior, posterior, and septal LACVs in the AF recurrence group were slower than those in the nonrecurrence group (anterior: 0.807 [0.766, 0.848] and 1.048 [1.000, 1.093] m/s, p < .05; posterior: 1.037 [0.991, 1.084] vs. 1.315 [1.249, 1.380] m/s, p < .05; septal: 0.904 [0.862, 0.946] vs. 1.163 [1.107, 1.219] m/s, p < .05). The best cut-off value of anterior LACV for predicting AF recurrence was 0.887 m/s (sensitivity 73.9% and specificity 76.5%). Multivariate analysis showed slow anterior LACV <0.887 m/s was an independent predictor of AF recurrence with an adjusted odds ratio of 1.42 (1.04, 1.94).

CONCLUSIONS

Slowing conduction velocity is a predictor of AF recurrence after pulmonary vein isolation.

Myofibroblasts impair myocardial impulse propagation by heterocellular connexin43 gap-junctional coupling through micropores

Aim: Composite population of myofibroblasts (MFs) within myocardial tissue is known to alter impulse propagation, leading to arrhythmias. However, it remains unclear whether and how MFs alter their propagation patterns when contacting cardiomyocytes (CMs) without complex structural insertions in the myocardium. We attempted to unveil the effects of the one-sided, heterocellular CM-MF connection on the impulse propagation of CM monolayers without the spatial insertion of MFs as an electrical or mechanical obstacle. Methods and results: We evaluated fluo8-based spatiotemporal patterns in impulse propagation of neonatal rat CM monolayers cultured on the microporous membrane having 8-μm diameter pores with co-culture of MFs or CMs on the reverse membrane side (CM-MF model or CM-CM model, respectively). During consecutive pacing at 1 or 2 Hz, the CM monolayers exhibited forward impulse propagation from the pacing site with a slower conduction velocity (θ) and a larger coefficient of directional θ variation in the CM-MF model than that in the CM-CM model in a frequency-dependent manner (2 Hz >1 Hz). The localized placement of an MF cluster on the reverse side resulted in an abrupt segmental depression of the impulse propagation of the upper CM layer, causing a spatiotemporally non-uniform pattern. Dye transfer of the calcein loaded in the upper CM layer to the lower MF layer was attenuated by the gap-junction inhibitor heptanol. Immunocytochemistry identified definitive connexin 43 (Cx43) between the CMs and MFs in the membrane pores. MF-selective Cx43 knockdown in the MF layer improved both the velocity and uniformity of propagation in the CM monolayer. Conclusion: Heterocellular Cx43 gap junction coupling of CMs with MFs alters the spatiotemporal patterns of myocardial impulse propagation, even in the absence of spatially interjacent and mechanosensitive modulations by MFs. Moreover, MFs can promote pro-arrhythmogenic impulse propagation when in face-to-face contact with the myocardium that arises in the healing infarct border zone.

Fibroblasts in heart scar tissue directly regulate cardiac excitability and arrhythmogenesis

After heart injury, dead heart muscle is replaced by scar tissue. Fibroblasts can electrically couple with myocytes, and changes in fibroblast membrane potential can lead to myocyte excitability, which suggests that fibroblast-myocyte coupling in scar tissue may be responsible for arrhythmogenesis. However, the physiologic relevance of electrical coupling of myocytes and fibroblasts and its impact on cardiac excitability in vivo have never been demonstrated. We genetically engineered a mouse that expresses the optogenetic cationic channel ChR2 (H134R) exclusively in cardiac fibroblasts. After myocardial infarction, optical stimulation of scar tissue elicited organ-wide cardiac excitation and induced arrhythmias in these animals. Complementing computational modeling with experimental approaches, we showed that gap junctional and ephaptic coupling, in a synergistic yet functionally redundant manner, excited myocytes coupled to fibroblasts.

Nonmyocytes as electrophysiological contributors to cardiac excitation and conduction

Although cardiac action potential (AP) generation and propagation have traditionally been attributed exclusively to cardiomyocytes (CM), other cell types in the heart are also capable of forming electrically conducting junctions. Interactions between CM and nonmyocytes (NM) enable and modulate each other's activity. This review provides an overview of the current understanding of heterocellular electrical communication in the heart. Although cardiac fibroblasts were initially thought to be electrical insulators, recent studies have demonstrated that they form functional electrical connections with CM in situ. Other NM, such as macrophages, have also been recognized as contributing to cardiac electrophysiology and arrhythmogenesis. Novel experimental tools have enabled the investigation of cell-specific activity patterns in native cardiac tissue, which is expected to yield exciting new insights into the development of novel or improved diagnostic and therapeutic strategies.

Cardiomyocyte electrophysiology and its modulation: current views and future prospects

Normal and abnormal cardiac rhythms are of key physiological and clinical interest. This introductory article begins from Sylvio Weidmann's key historic 1950s microelectrode measurements of cardiac electrophysiological activity and Singh & Vaughan Williams's classification of cardiotropic targets. It then proceeds to introduce the insights into cardiomyocyte function and its regulation that subsequently emerged and their therapeutic implications. We recapitulate the resulting view that surface membrane electrophysiological events underlying cardiac excitation and its initiation, conduction and recovery constitute the final common path for the cellular mechanisms that impinge upon this normal or abnormal cardiac electrophysiological activity. We then consider progress in the more recently characterized successive regulatory hierarchies involving Ca2+ homeostasis, excitation-contraction coupling and autonomic G-protein signalling and their often reciprocal interactions with the surface membrane events, and their circadian rhythms. Then follow accounts of longer-term upstream modulation processes involving altered channel expression, cardiomyocyte energetics and hypertrophic and fibrotic cardiac remodelling. Consideration of these developments introduces each of the articles in this Phil. Trans. B theme issue. The findings contained in these articles translate naturally into recent classifications of cardiac electrophysiological targets and drug actions, thereby encouraging future iterations of experimental cardiac electrophysiological discovery, and testing directed towards clinical management. This article is part of the theme issue 'The heartbeat: its molecular basis and physiological mechanisms'.

Ephaptic Coupling as a Resolution to the Paradox of Action Potential Wave Speed and Discordant Alternans Spatial Scales in the Heart

Previous computer simulations have suggested that existing models of action potential wave propagation in the heart are not consistent with observed wave propagation behavior. Specifically, computer models cannot simultaneously reproduce the rapid wave speeds and small spatial scales of discordant alternans patterns measured experimentally in the same simulation. The discrepancy is important, because discordant alternans can be a key precursor to the development of abnormal and dangerous rapid rhythms in the heart. In this Letter, we show that this paradox can be resolved by allowing so-called ephaptic coupling to play a primary role in wave front propagation in place of conventional gap-junction coupling. With this modification, physiological wave speeds and small discordant alternans spatial scales both occur with gap-junction resistance values that are more in line with those observed in experiments. Our theory thus also provides support to the hypothesis that ephaptic coupling plays an important role in normal wave propagation.

Role of ephaptic coupling in discordant alternans domain sizes and action potential propagation in the heart

Discordant alternans, the spatially out-of-phase alternation of the durations of propagating action potentials in the heart, has been linked to the onset of fibrillation, a major cardiac rhythm disorder. The sizes of the regions, or domains, within which these alternations are synchronized are critical in this link. However, computer models employing standard gap junction-based coupling between cells have been unable to reproduce simultaneously the small domain sizes and rapid action potential propagation speeds seen in experiments. Here we use computational methods to show that rapid wave speeds and small domain sizes are possible when a more detailed model of intercellular coupling that accounts for so-called ephaptic effects is used. We provide evidence that the smaller domain sizes are possible, because different coupling strengths can exist on the wavefronts, for which both ephaptic and gap-junction coupling are involved, in contrast to the wavebacks, where only gap-junction coupling plays an active role. The differences in coupling strength are due to the high density of fast-inward (sodium) channels known to localize on the ends of cardiac cells, which are only active (and thus engage ephaptic coupling) during wavefront propagation. Thus, our results suggest that this distribution of fast-inward channels, as well as other factors responsible for the critical involvement of ephaptic coupling in wave propagation, including intercellular cleft spacing, play important roles in increasing the vulnerability of the heart to life-threatening tachyarrhythmias. Our results, combined with the absence of short-wavelength discordant alternans domains in standard gap-junction-dominated coupling models, also provide evidence that both gap-junction and ephaptic coupling are critical in wavefront propagation and waveback dynamics.

Abnormal Conduction Zone Detected by Isochronal Late Activation Mapping Accurately Identifies the Potential Atrial Substrate and Predicts the Atrial Fibrillation Ablation Outcome After Pulmonary Vein Isolation

BACKGROUND

The presence of abnormal substrate of left atrium is a predictor of atrial fibrillation (AF) recurrence after pulmonary vein isolation. We aimed to investigate the isochronal late activation mapping to access the abnormal conduction velocity for predicting AF ablation outcome.

METHODS

Forty-five paroxysmal AF patients (30 males, 57.8±8.7 years old) who underwent pulmonary vein isolation were enrolled. Isochronal late activation mapping was retrospectively constructed with 2 different windows of interest: from onset of P wave to onset of QRS wave on surface electrocardiography (W1) and 74 ms tracking back from the end of P wave (W2). Deceleration zone was defined as regions with 3 isochrones (DZa) or ≥4 isochrones (DZb) within a 1 cm radius on the isochronal late activation mapping, and the estimated conduction velocity (ECV) are 0.27 m/s and <0.20 m/s for DZa and DZb, respectively in W2. The distribution of deceleration zone was compared with the location of low-voltage zone (bipolar voltage ≤0.5 mV). Any recurrence of atrial arrhythmias was defined as the primary end point during follow ups after a 3-month blanking period.

RESULTS

Pulmonary vein isolation was performed in all patients, and there were 2 patients (4.4%) received additional extrapulmonary vein ablation. After a mean follow-up of 12.7±4.5 months, recurrence of AF occurred in 14 patients (31.1%). Patients with the presence of DZb in W2 had higher AF recurrence (Kaplan-Meier event rate estimates: HR, 9.41 [95% CI, 2.61-33.90]; log-rank P<0.0001). There were 52.6% of the DZb locations in W2 comparable to the distributions of low-voltage zone and 47.4% DZb were distributed in the area without low-voltage zone.

CONCLUSIONS

Deceleration zone detected by isochronal late activation mapping represents a critical AF substrate, it accurately predicts the AF recurrence following ablation in patients with paroxysmal AF.

Computational design of custom therapeutic cells to correct failing human cardiomyocytes

Background

Myocardial delivery of non-excitable cells-namely human mesenchymal stem cells (hMSCs) and c-kit⁺ cardiac interstitial cells (hCICs)-remains a promising approach for treating the failing heart. Recent empirical studies attempt to improve such therapies by genetically engineering cells to express specific ion channels, or by creating hybrid cells with combined channel expression. This study uses a computational modeling approach to test the hypothesis that custom hypothetical cells can be rationally designed to restore a healthy phenotype when coupled to human heart failure (HF) cardiomyocytes.

Methods

Candidate custom cells were simulated with a combination of ion channels from non-excitable cells and healthy human cardiomyocytes (hCMs). Using a genetic algorithm-based optimization approach, candidate cells were accepted if a root mean square error (RMSE) of less than 50% relative to healthy hCM was achieved for both action potential and calcium transient waveforms for the cell-treated HF cardiomyocyte, normalized to the untreated HF cardiomyocyte.

Results

Custom cells expressing only non-excitable ion channels were inadequate to restore a healthy cardiac phenotype when coupled to either fibrotic or non-fibrotic HF cardiomyocytes. In contrast, custom cells also expressing cardiac ion channels led to acceptable restoration of a healthy cardiomyocyte phenotype when coupled to fibrotic, but not non-fibrotic, HF cardiomyocytes. Incorporating the cardiomyocyte inward rectifier K⁺ channel was critical to accomplishing this phenotypic rescue while also improving single-cell action potential metrics associated with arrhythmias, namely resting membrane potential and action potential duration. The computational approach also provided insight into the rescue mechanisms, whereby heterocellular coupling enhanced cardiomyocyte L-type calcium current and promoted calcium-induced calcium release. Finally, as a therapeutically translatable strategy, we simulated delivery of hMSCs and hCICs genetically engineered to express the cardiomyocyte inward rectifier K⁺ channel, which decreased action potential and calcium transient RMSEs by at least 24% relative to control hMSCs and hCICs, with more favorable single-cell arrhythmia metrics.

Conclusion

Computational modeling facilitates exploration of customizable engineered cell therapies. Optimized cells expressing cardiac ion channels restored healthy action potential and calcium handling phenotypes in fibrotic HF cardiomyocytes and improved single-cell arrhythmia metrics, warranting further experimental validation studies of the proposed custom therapeutic cells.

Regional left atrial conduction velocity in the anterior wall is associated with clinical recurrence of atrial fibrillation after catheter ablation: efficacy in combination with the ipsilateral low voltage area

Background

Left atrial (LA) conduction velocity (CV) is an electrical remodeling parameter of atrial fibrillation (AF) substrate. However, the pathophysiological substrate of LA-CV and its impact on outcomes after catheter ablation for AF have not been well evaluated.

Methods

We retrospectively evaluated 119 patients with AF who underwent catheter ablation and electroanatomical mapping during sinus rhythm. To measure regional LA-CV, we took triplet sites (A, B, and C) on the activation map and calculated the magnitude of the matched orthogonal projection vector between vector-AB and vector-AC, indicating two-dimensional CV. The median of the LA-CVs from four triad sites in both the anterior and posterior walls was set as the ‘model LA-CV’. We evaluated the impact of the model LA-CV on recurrence after ablation and relationship between the model LA-CV and LA-low voltage area (LVA) of < 0.5 mV.

Results

During the 12-month follow-up, 29 patients experienced recurrence. The LA-CV model was significantly correlated with ipsilateral LVA. The lower anterior model LA-CV was significantly associated with recurrence, with the cut-off value of 0.80 m/s having a sensitivity of 72% and specificity of 67%. Multivariable analysis revealed that the anterior model LA-CV (hazard ratio, 0.09; 95% confidence interval, 0.01–0.94; p = 0.043) and anterior LA-LVA (hazard ratio, 1.06; 95% confidence interval, 1.00–1.11; p = 0.033) were independently associated with AF recurrence. The anterior LA-LVA was mildly correlated with the anterior model LA-CV (r = -0.358; p < 0.001), and patients with both lower LA-CV and greater anterior LA-LVA based on each cut-off value had the worst prognosis. However, decreased LA-CV was more likely to be affected by the distribution pattern of the LVA rather than the total size of the LVA.

Conclusion

Decreased anterior LA-CV was a significant predictor of AF recurrence and was a useful electrical parameter in addition to LA-LVA for estimating AF arrhythmogenicity.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12872-022-02881-6.

The physics of heart rhythm disorders

The global burden caused by cardiovascular disease is substantial, with heart disease representing the most common cause of death around the world. There remains a need to develop better mechanistic models of cardiac function in order to combat this health concern. Heart rhythm disorders, or arrhythmias, are one particular type of disease which has been amenable to quantitative investigation. Here we review the application of quantitative methodologies to explore dynamical questions pertaining to arrhythmias. We begin by describing single-cell models of cardiac myocytes, from which two and three dimensional models can be constructed. Special focus is placed on results relating to pattern formation across these spatially-distributed systems, especially the formation of spiral waves of activation. Next, we discuss mechanisms which can lead to the initiation of arrhythmias, focusing on the dynamical state of spatially discordant alternans, and outline proposed mechanisms perpetuating arrhythmias such as fibrillation. We then review experimental and clinical results related to the spatio-temporal mapping of heart rhythm disorders. Finally, we describe treatment options for heart rhythm disorders and demonstrate how statistical physics tools can provide insights into the dynamics of heart rhythm disorders.

Effects of fibroblast on electromechanical dynamics of human atrial tissue—insights from a 2D discrete element model

Roughly 75% of normal myocardial tissue volume is comprised of myocytes, however, fibroblasts by number are the most predominant cells in cardiac tissue. Previous studies have shown distinctive differences in cellular electrophysiology and excitability between myocytes and fibroblasts. However, it is still unclear how the electrical coupling between the two and the increased population of fibroblasts affects the electromechanical dynamics of cardiac tissue. This paper focuses on investigating effects of fibroblast-myocyte electrical coupling (FMEC) and fibroblast population on atrial electrical conduction and mechanical contractility by using a two-dimensional Discrete Element Method (DEM) model of cardiac tissue that is different to finite element method (FEM). In the model, the electro-mechanics of atrial cells are modelled by a biophysically detailed model for atrial electrical action potentials and myofilament kinetics, and the atrial fibroblasts are modelled by an active model that considers four active membrane ionic channel currents. Our simulation results show that the FMEC impairs myocytes’ electrical action potential and mechanical contractibility, manifested by reduced upstroke velocity, amplitude and duration of action potentials, as well as cell length shortening. At the tissue level, the FMEC slows down the conduction of excitation waves, and reduces strain of the tissue produced during a contraction course. These findings provide new insights into understandings of how FMEC impairs cardiac electrical and mechanical dynamics of the heart.

In silico Cell Therapy Model Restores Failing Human Myocyte Electrophysiology and Calcium Cycling in Fibrotic Myocardium

Myocardial delivery of human c-kit+ cardiac interstitial cells (hCICs) and human mesenchymal stem cells (hMSCs), an emerging approach for treating the failing heart, has been limited by an incomplete understanding of the effects on host myocardium. This computational study aims to model hCIC and hMSC effects on electrophysiology and calcium cycling of healthy and diseased human cardiomyocytes (hCM), and reveals a possible cardiotherapeutic benefit independent of putative regeneration processes. First, we developed an original hCIC mathematical model with an electrical profile comprised of distinct experimentally identified ion currents. Next, we verified the model by confirming it is representative of published experiments on hCIC whole-cell electrophysiology and on hCIC co-cultures with rodent cardiomyocytes. We then used our model to compare electrophysiological effects of hCICs to other non-excitable cells, as well as clinically relevant hCIC-hMSC combination therapies and fused hCIC-hMSC CardioChimeras. Simulation of direct coupling of hCICs to healthy or failing hCMs through gap junctions led to greater increases in calcium cycling with lesser reductions in action potential duration (APD) compared with hMSCs. Combined coupling of hCICs and hMSCs to healthy or diseased hCMs led to intermediate effects on electrophysiology and calcium cycling compared to individually coupled hCICs or hMSCs. Fused hCIC-hMSC CardioChimeras decreased healthy and diseased hCM APD and calcium transient amplitude compared to individual or combined cell treatments. Finally, to provide a theoretical basis for optimizing cell-based therapies, we randomized populations of 2,500 models incorporating variable hMSC and hCIC interventions and simulated their effects on restoring diseased cardiomyocyte electrophysiology and calcium handling. The permutation simulation predicted the ability to correct abnormal properties of heart failure hCMs in fibrotic, but not non-fibrotic, myocardium. This permutation experiment also predicted paracrine signaling to be a necessary and sufficient mechanism for this correction, counteracting the fibrotic effects while also restoring arrhythmia-related metrics such as upstroke velocity and resting membrane potential. Altogether, our in silico findings suggest anti-fibrotic effects of paracrine signaling are critical to abrogating pathological cardiomyocyte electrophysiology and calcium cycling in fibrotic heart failure, and support further investigation of delivering an optimized cellular secretome as a potential strategy for improving heart failure therapy.

Fibrotic Remodeling during Persistent Atrial Fibrillation: In Silico Investigation of the Role of Calcium for Human Atrial Myofibroblast Electrophysiology

During atrial fibrillation, cardiac tissue undergoes different remodeling processes at different scales from the molecular level to the tissue level. One central player that contributes to both electrical and structural remodeling is the myofibroblast. Based on recent experimental evidence on myofibroblasts' ability to contract, we extended a biophysical myofibroblast model with Ca2+ handling components and studied the effect on cellular and tissue electrophysiology. Using genetic algorithms, we fitted the myofibroblast model parameters to the existing in vitro data. In silico experiments showed that Ca2+ currents can explain the experimentally observed variability regarding the myofibroblast resting membrane potential. The presence of an L-type Ca2+ current can trigger automaticity in the myofibroblast with a cycle length of 799.9 ms. Myocyte action potentials were prolonged when coupled to myofibroblasts with Ca2+ handling machinery. Different spatial myofibroblast distribution patterns increased the vulnerable window to induce arrhythmia from 12 ms in non-fibrotic tissue to 22 ± 2.5 ms and altered the reentry dynamics. Our findings suggest that Ca2+ handling can considerably affect myofibroblast electrophysiology and alter the electrical propagation in atrial tissue composed of myocytes coupled with myofibroblasts. These findings can inform experimental validation experiments to further elucidate the role of myofibroblast Ca2+ handling in atrial arrhythmogenesis.

Micropattern platform promotes extracellular matrix remodeling by human PSC-derived cardiac fibroblasts and enhances contractility of co-cultured cardiomyocytes

In native heart tissue, cardiac fibroblasts provide the structural framework of extracellular matrix (ECM) while also influencing the electrical and mechanical properties of cardiomyocytes. Recent advances in the field of stem cell differentiation have led to the availability of human pluripotent stem cell-derived cardiac fibroblasts (iPSC-CFs) in addition to cardiomyocytes (iPSC-CMs). Here we use a novel 2D in vitro micropatterned platform that provides control over ECM geometry and substrate stiffness. When cultured alone on soft micropatterned substrates, iPSC-CFs are confined to the micropatterned features and remodel the ECM into anisotropic fibers. Similar remodeling and ECM production occurs when cultured with iPSC-CMs in a co-culture model. In addition to modifications in the ECM, our results show that iPSC-CFs influence iPSC-CM function with accelerated Ca2+ transient rise-up time and greater contractile strains in the co-culture conditions compared to when iPSC-CMs are cultured alone. These combined observations highlight the important role cardiac fibroblasts play in vivo and the need for co-culture models like the one presented here to provide more representative in vitro cardiac constructs.

Pulsed low-energy stimulation initiates electric turbulence in cardiac tissue

Interruptions in nonlinear wave propagation, commonly referred to as wave breaks, are typical of many complex excitable systems. In the heart they lead to lethal rhythm disorders, the so-called arrhythmias, which are one of the main causes of sudden death in the industrialized world. Progress in the treatment and therapy of cardiac arrhythmias requires a detailed understanding of the triggers and dynamics of these wave breaks. In particular, two very important questions are: 1) What determines the potential of a wave break to initiate re-entry? and 2) How do these breaks evolve such that the system is able to maintain spatiotemporally chaotic electrical activity? Here we approach these questions numerically using optogenetics in an in silico model of human atrial tissue that has undergone chronic atrial fibrillation (cAF) remodelling. In the lesser studied sub-threshold illumination régime, we discover a new mechanism of wave break initiation in cardiac tissue that occurs for gentle slopes of the restitution characteristics. This mechanism involves the creation of conduction blocks through a combination of wavefront-waveback interaction, reshaping of the wave profile and heterogeneous recovery from the excitation of the spatially extended medium, leading to the creation of re-excitable windows for sustained re-entry. This finding is an important contribution to cardiac arrhythmia research as it identifies scenarios in which low-energy perturbations to cardiac rhythm can be potentially life-threatening.

Are Interactions between Epicardial Adipose Tissue, Cardiac Fibroblasts and Cardiac Myocytes Instrumental in Atrial Fibrosis and Atrial Fibrillation?

Atrial fibrillation is very common among the elderly and/or obese. While myocardial fibrosis is associated with atrial fibrillation, the exact mechanisms within atrial myocytes and surrounding non-myocytes are not fully understood. This review considers the potential roles of myocardial fibroblasts and myofibroblasts in fibrosis and modulating myocyte electrophysiology through electrotonic interactions. Coupling with (myo)fibroblasts in vitro and in silico prolonged myocyte action potential duration and caused resting depolarization; an optogenetic study has verified in vivo that fibroblasts depolarized when coupled myocytes produced action potentials. This review also introduces another non-myocyte which may modulate both myocardial (myo)fibroblasts and myocytes: epicardial adipose tissue. Epicardial adipocytes are in intimate contact with myocytes and (myo)fibroblasts and may infiltrate the myocardium. Adipocytes secrete numerous adipokines which modulate (myo)fibroblast and myocyte physiology. These adipokines are protective in healthy hearts, preventing inflammation and fibrosis. However, adipokines secreted from adipocytes may switch to pro-inflammatory and pro-fibrotic, associated with reactive oxygen species generation. Pro-fibrotic adipokines stimulate myofibroblast differentiation, causing pronounced fibrosis in the epicardial adipose tissue and the myocardium. Adipose tissue also influences myocyte electrophysiology, via the adipokines and/or through electrotonic interactions. Deeper understanding of the interactions between myocytes and non-myocytes is important to understand and manage atrial fibrillation.

Increased fibroblast accumulation in the equine heart following persistent atrial fibrillation

BACKGROUND

Fibroblasts maintain the extracellular matrix homeostasis and may couple to cardiomyocytes through gap junctions and thereby increase the susceptibility to slow conduction and cardiac arrhythmias, such as atrial fibrillation (AF). In this study, we used an equine model of persistent AF to characterize structural changes and the role of fibroblasts in the development of an arrhythmogenic substrate for AF.

MATERIAL AND METHODS

Eleven horses were subjected to atrial tachypacing until self-sustained AF developed and were kept in AF for six weeks. Horses in sinus rhythm (SR) served as control. In terminal open-chest experiments conduction velocity (CV) was measured. Tissue was harvested and stained from selected sites. Automated image analysis was performed to assess fibrosis, fibroblasts, capillaries and various cardiomyocyte characteristics.

RESULTS

Horses in SR showed a rate-dependent slowing of CV, while in horses with persistent AF this rate-dependency was completely abolished (CV•basic cycle length relation p = 0.0295). Overall and interstitial amounts of fibrosis were unchanged, but an increased fibroblast count was found in left atrial appendage, Bachmann's bundle, intraatrial septum and pulmonary veins (p < 0.05 for all) in horses with persistent AF. The percentage of α-SMA expressing fibroblasts remained the same between the groups.

CONCLUSION

Persistent AF resulted in fibroblast accumulation in several regions, particularly in the left atrial appendage. The increased number of fibroblasts could be a mediator of altered electrophysiology during AF. Targeting the fibroblast proliferation and differentiation could potentially serve as a novel therapeutic target slowing down the structural remodeling associated with AF.

A new model of myofibroblast-cardiomyocyte interactions and their differences across species

Although coupling between cardiomyocytes and myofibroblasts is well known to affect the physiology and pathophysiology of cardiac tissues across species, relating these observations to humans is challenging because the effect of this coupling varies across species and because the sources of these effects are not known. To identify the sources of cross-species variation, we built upon previous mathematical models of myofibroblast electrophysiology and developed a mechanoelectrical model of cardiomyocyte-myofibroblast interactions as mediated by electrotonic coupling and transforming growth factor-β1. The model, as verified by experimental data from the literature, predicted that both electrotonic coupling and transforming growth factor-β1 interaction between myocytes and myofibroblast prolonged action potential in rat myocytes but shortened action potential in human myocytes. This variance could be explained by differences in the transient outward K⁺ current associated with differential Kv4.2 gene expression across species. Results are useful for efforts to extrapolate the results of animal models to the predicted effects in humans and point to potential therapeutic targets for fibrotic cardiomyopathy.

Intercellular Communication in the Heart: Therapeutic Opportunities for Cardiac Ischemia

The maintenance of tissue, organ, and organism homeostasis relies on an intricate network of players and mechanisms that assist in the different forms of cell-cell communication. Myocardial infarction, following heart ischemia and reperfusion, is associated with profound changes in key processes of intercellular communication, involving gap junctions, extracellular vesicles, and tunneling nanotubes, some of which have been implicated in communication defects associated with cardiac injury, namely arrhythmogenesis and progression into heart failure. Therefore, intercellular communication players have emerged as attractive powerful therapeutic targets aimed at preserving a fine-tuned crosstalk between the different cardiac cells in order to prevent or repair some of harmful consequences of heart ischemia and reperfusion, re-establishing myocardial function.

Action potential propagation and block in a model of atrial tissue with myocyte-fibroblast coupling

The electrical coupling between myocytes and fibroblasts and the spacial distribution of fibroblasts within myocardial tissues are significant factors in triggering and sustaining cardiac arrhythmias, but their roles are poorly understood. This article describes both direct numerical simulations and an asymptotic theory of propagation and block of electrical excitation in a model of atrial tissue with myocyte-fibroblast coupling. In particular, three idealized fibroblast distributions are introduced: uniform distribution, fibroblast barrier and myocyte strait-all believed to be constituent blocks of realistic fibroblast distributions. Primary action potential biomarkers including conduction velocity, peak potential and triangulation index are estimated from direct simulations in all cases. Propagation block is found to occur at certain critical values of the parameters defining each idealized fibroblast distribution, and these critical values are accurately determined. An asymptotic theory proposed earlier is extended and applied to the case of a uniform fibroblast distribution. Biomarker values are obtained from hybrid analytical-numerical solutions of coupled fast-time and slow-time periodic boundary value problems and compare well to direct numerical simulations. The boundary of absolute refractoriness is determined solely by the fast-time problem and is found to depend on the values of the myocyte potential and on the slow inactivation variable of the sodium current ahead of the propagating pulse. In turn, these quantities are estimated from the slow-time problem using a regular perturbation expansion to find the steady state of the coupled myocyte-fibroblast kinetics. The asymptotic theory gives a simple analytical expression that captures with remarkable accuracy the block of propagation in the presence of fibroblasts.

Slow whole left atrial conduction velocity after pulmonary vein isolation predicts atrial fibrillation recurrence

BACKGROUND

Atrial conduction velocity may represent atrial fibrillation (AF) substrate after pulmonary vein isolation (PVI). To elucidate the association between whole left atrial conduction velocity (LACV) and AF recurrence after PVI.

METHODS AND RESULTS

This observational study enrolled 279 patients (147 paroxysmal and 132 persistent AF) who underwent PVI alone as an initial AF ablation procedure. After PVI, the left atrium was mapped with a 20-pole multielectrode in conjunction with the CARTO3 system during 100-ppm right atrial pacing. Left atrial conduction distance and conduction time were calculated from the start to the end of the propagation wave front in the left atrium. LACVs on the anterior and posterior routes were calculated as conduction distance divided by conduction time. Anterior and posterior LACVs were slower in patients with AF recurrence than in those without (anterior, 0.79 [0.71, 0.86] vs 0.96 [0.90, 1.06] m/s, P < .001; posterior, 0.99 [0.89, 1.14] vs 1.10 [1.00, 1.29] m/s, P < .001). AF recurrence was best predicted by anterior LACV with a cut-off value of 0.87 m/s (sensitivity 87%, specificity 81%, and predictive accuracy 84%). Multivariate analysis demonstrated that a slow anterior LACV <0.87 m/s was an independent predictor of AF recurrence with an adjusted hazard ratio of 11.8 (6.36-22.0). Patients with anterior low-voltage areas demonstrated slower anterior LACV than those without low-voltage areas (0.89 [0.71, 1.00] vs 0.94 [0.87, 1.05] m/s, P < .001).

CONCLUSION

A slow anterior LACV was an excellent predictor of AF recurrence after PVI.

Molecular Signatures of Sinus Node Dysfunction Induce Structural Remodeling in the Right Atrial Tissue

The sinus node (SN) is located at the apex of the cardiac conduction system, and SN dysfunction (SND)-characterized by electrical remodeling-is generally attributed to idiopathic fibrosis or ischemic injuries in the SN. SND is associated with increased risk of cardiovascular disorders, including syncope, heart failure, and atrial arrhythmias, particularly atrial fibrillation. One of the histological SND hallmarks is degenerative atrial remodeling that is associated with conduction abnormalities and increased right atrial refractoriness. Although SND is frequently accompanied by increased fibrosis in the right atrium (RA), its molecular basis still remains elusive. Therefore, we investigated whether SND can induce significant molecular changes that account for the structural remodeling of RA. Towards this, we employed a rabbit model of experimental SND, and then compared the genome-wide RNA expression profiles in RA between SND-induced rabbits and sham-operated controls to identify the differentially expressed transcripts. The accompanying gene enrichment analysis revealed extensive pro-fibrotic changes within 7 days after the SN ablation, including activation of transforming growth factor-β (TGF-β) signaling and alterations in the levels of extracellular matrix components and their regulators. Importantly, our findings suggest that periostin, a matricellular factor that regulates the development of cardiac tissue, might play a key role in mediating TGF-β-signaling-induced aberrant atrial remodeling. In conclusion, the present study provides valuable information regarding the molecular signatures underlying SND-induced atrial remodeling, and indicates that periostin can be potentially used in the diagnosis of fibroproliferative cardiac dysfunctions.

Engineering microenvironment for human cardiac tissue assembly in heart-on-a-chip platform

Organ-on-a-chip systems have the potential to revolutionize drug screening and disease modeling through the use of human stem cell-derived cardiomyocytes. The predictive power of these tissue models critically depends on the functional assembly and maturation of human cells that are used as building blocks for organ-on-a-chip systems. To resemble a more adult-like phenotype on these heart-on-a-chip systems, the surrounding micro-environment of individual cardiomyocyte needs to be controlled. Herein, we investigated the impact of four microenvironmental cues: cell seeding density, types and percentages of non-myocyte populations, the types of hydrogels used for tissue inoculation and the electrical conditioning regimes on the structural and functional assembly of human pluripotent stem cell-derived cardiac tissues. Utilizing a novel, plastic and open-access heart-on-a-chip system that is capable of continuous non-invasive monitoring of tissue contractions, we were able to study how different micro-environmental cues affect the assembly of the cardiomyocytes into a functional cardiac tissue. We have defined conditions that resulted in tissues exhibiting hallmarks of the mature human myocardium, such as positive force-frequency relationship and post-rest potentiation.

Melatonin receptor activation protects against low potassium-induced ventricular fibrillation by preserving action potentials and connexin-43 topology in isolated rat hearts

Hypokalemia prolongs the QRS and QT intervals, deteriorates intercellular coupling, and increases the risk for arrhythmia. Melatonin preserves gap junctions and shortens action potential as potential antiarrhythmic mechanisms, but its properties under hypokalemia remain unknown. We hypothesized that melatonin protects against low potassium-induced arrhythmias through the activation of its receptors, resulting in action potential shortening and connexin-43 preservation. After stabilization in Krebs-Henseleit solution (4.5 mEq/L K⁺ ), isolated hearts from Wistar rats underwent perfusion with low-potassium (1 mEq/L) solution and melatonin (100 μmol/L), a melatonin receptor blocker (luzindole, 5 μmol/L), melatonin + luzindole or vehicle. The primary endpoint of the study was the prevention of ventricular fibrillation. Electrocardiography was used, and epicardial action potentials and heart function were measured and analyzed. The ventricular expression, dephosphorylation, and distribution of connexin-43 were examined. Melatonin reduced the incidence of low potassium-induced ventricular fibrillation from 100% to 59%, delayed the occurrence of ventricular fibrillation and induced a faster recovery of sinus rhythm during potassium restitution. Melatonin prevented QRS widening, action potential activation delay, and the prolongation of action potential duration at 50% of repolarization. Other ECG and action potential parameters, the left ventricular developed pressure, and nonsustained ventricular arrhythmias did not differ among groups. Melatonin prevented connexin-43 dephosphorylation and its abnormal topology (lateralization). Luzindole abrogated the protective effects of melatonin on electrophysiological properties and connexin-43 misdistribution. Our results indicate that melatonin receptor activation protects against low potassium-induced ventricular fibrillation, shortens action potential duration, preserves ventricular electrical activation, and prevents acute changes in connexin-43 distribution. All of these properties make melatonin a remarkable antifibrillatory agent.

Acetylcholine Delays Atrial Activation to Facilitate Atrial Fibrillation

Background

Acetylcholine (ACh) shortens action potential duration (APD) in human atria. APD shortening facilitates atrial fibrillation (AF) by reducing the wavelength for reentry. However, the influence of ACh on electrical conduction in human atria and its contribution to AF are unclear, particularly when combined with impaired conduction from interstitial fibrosis.

Objective

To investigate the effect of ACh on human atrial conduction and its role in AF with computational, experimental, and clinical approaches.

Methods

S1S2 pacing (S1 = 600 ms and S2 = variable cycle lengths) was applied to the following human AF computer models: a left atrial appendage (LAA) myocyte to quantify the effects of ACh on APD, maximum upstroke velocity (V _max ), and resting membrane potential (RMP); a monolayer of LAA myocytes to quantify the effects of ACh on conduction; and 3) an intact left atrium (LA) to determine the effects of ACh on arrhythmogenicity. Heterogeneous ACh and interstitial fibrosis were applied to the monolayer and LA models. To corroborate the simulations, APD and RMP from isolated human atrial myocytes were recorded before and after 0.1 μM ACh. At the tissue level, LAAs from AF patients were optically mapped ex vivo using Di-4-ANEPPS. The difference in total activation time (AT) was determined between AT initially recorded with S1 pacing, and AT recorded during subsequent S1 pacing without (n = 6) or with (n = 7) 100 μM ACh.

Results

In LAA myocyte simulations, S1 pacing with 0.1 μM ACh shortened APD by 41 ms, hyperpolarized RMP by 7 mV, and increased V _max by 27 mV/ms. In human atrial myocytes, 0.1 μM ACh shortened APD by 48 ms, hyperpolarized RMP by 3 mV, and increased V _max by 6 mV/ms. In LAA monolayer simulations, S1 pacing with ACh hyperpolarized RMP to delay total AT by 32 ms without and 35 ms with fibrosis. This led to unidirectional conduction block and sustained reentry in fibrotic LA with heterogeneous ACh during S2 pacing. In AF patient LAAs, S1 pacing with ACh increased total AT from 39.3 ± 26 ms to 71.4 ± 31.2 ms (p = 0.036) compared to no change without ACh (56.7 ± 29.3 ms to 50.0 ± 21.9 ms, p = 0.140).

Conclusion

In fibrotic atria with heterogeneous parasympathetic activation, ACh facilitates AF by shortening APD and slowing conduction to promote unidirectional conduction block and reentry.

Heterogeneous Effects of Fibroblast-Myocyte Coupling in Different Regions of the Human Atria Under Conditions of Atrial Fibrillation

Background: Atrial fibrillation (AF), the most common cardiac arrhythmia, is characterized by alteration of the action potential (AP) propagation. Under persistent AF, myocytes undergo electrophysiological and structural remodeling, which involves fibroblast proliferation and differentiation, modifying the substrate for AP propagation. The aim of this study was to analyze the effects on the AP of fibroblast-myocyte coupling during AF and its propagation in different regions of the atria.

Methods: Isolated myocytes were coupled to different numbers of fibroblasts using the established AP models and tissue simulations were performed by randomly distributing fibroblasts. Fibroblast formulations were updated to match recent experimental data. Major ion current conductances of the myocyte model were modified to simulate AP heterogeneity in four different atrial regions (right atrium posterior wall, crista terminalis, left atrium posterior wall, and pulmonary vein) according to experimental and computational studies.

Results: The results of the coupled myocyte-fibroblast simulations suggest that a more depolarized membrane potential and higher fibroblast membrane capacitance have a greater impact on AP duration and myocyte maximum depolarization velocity. The number of coupled fibroblasts and the stimulation frequency are determining factors in altering myocyte AP. Strand simulations show that conduction velocity tends to homogenize in all regions, while the left atrium is more likely to be affected by fibroblast and AP propagation block is more likely to occur. The pulmonary vein is the most affected region, even at low fibroblast densities. In 2D sheets with randomly placed fibroblasts, wavebreaks are observed in the low density (10%) central fibrotic zone and when fibroblast density increases (40%) propagation in the fibrotic region is practically blocked. At densities of 10 and 20% the width of the vulnerable window increases with respect to control but is decreased at 40%.

Conclusion: Myocyte-fibroblast coupling characteristics heterogeneously affect AP propagation and features in the different atrial zones, and myocytes from the left atria are more sensitive to fibroblast coupling.

Personalized Cardiac Computational Models: From Clinical Data to Simulation of Infarct-Related Ventricular Tachycardia

In the chronic stage of myocardial infarction, a significant number of patients develop life-threatening ventricular tachycardias (VT) due to the arrhythmogenic nature of the remodeled myocardium. Radiofrequency ablation (RFA) is a common procedure to isolate reentry pathways across the infarct scar that are responsible for VT. Unfortunately, this strategy show relatively low success rates; up to 50% of patients experience recurrent VT after the procedure. In the last decade, intensive research in the field of computational cardiac electrophysiology (EP) has demonstrated the ability of three-dimensional (3D) cardiac computational models to perform in-silico EP studies. However, the personalization and modeling of certain key components remain challenging, particularly in the case of the infarct border zone (BZ). In this study, we used a clinical dataset from a patient with a history of infarct-related VT to build an image-based 3D ventricular model aimed at computational simulation of cardiac EP, including detailed patient-specific cardiac anatomy and infarct scar geometry. We modeled the BZ in eight different ways by combining the presence or absence of electrical remodeling with four different levels of image-based patchy fibrosis (0, 10, 20, and 30%). A 3D torso model was also constructed to compute the ECG. Patient-specific sinus activation patterns were simulated and validated against the patient's ECG. Subsequently, the pacing protocol used to induce reentrant VTs in the EP laboratory was reproduced in-silico. The clinical VT was induced with different versions of the model and from different pacing points, thus identifying the slow conducting channel responsible for such VT. Finally, the real patient's ECG recorded during VT episodes was used to validate our simulation results and to assess different strategies to model the BZ. Our study showed that reduced conduction velocities and heterogeneity in action potential duration in the BZ are the main factors in promoting reentrant activity. Either electrical remodeling or fibrosis in a degree of at least 30% in the BZ were required to initiate VT. Moreover, this proof-of-concept study confirms the feasibility of developing 3D computational models for cardiac EP able to reproduce cardiac activation in sinus rhythm and during VT, using exclusively non-invasive clinical data.

Atrial Transcriptional Profiles of Molecular Targets Mediating Electrophysiological Function in Aging and Pgc-1β Deficient Murine Hearts

BACKGROUND

Deficiencies in the transcriptional co-activator, peroxisome proliferative activated receptor, gamma, coactivator-1β are implicated in deficient mitochondrial function. The latter accompanies clinical conditions including aging, physical inactivity, obesity, and diabetes. Recent electrophysiological studies reported that Pgc-1β^-/- mice recapitulate clinical age-dependent atrial pro-arrhythmic phenotypes. They implicated impaired chronotropic responses to adrenergic challenge, compromised action potential (AP) generation and conduction despite normal AP recovery timecourses and background resting potentials, altered intracellular Ca²⁺ homeostasis, and fibrotic change in the observed arrhythmogenicity.

OBJECTIVE

We explored the extent to which these age-dependent physiological changes correlated with alterations in gene transcription in murine Pgc-1β^-/- atria.

METHODS AND RESULTS

RNA isolated from murine atrial tissue samples from young (12-16 weeks) and aged (>52 weeks of age), wild type (WT) and Pgc-1β^-/- mice were studied by pre-probed quantitative PCR array cards. We examined genes encoding sixty ion channels and other strategic atrial electrophysiological proteins. Pgc-1β^-/- genotype independently reduced gene transcription underlying Na⁺-K⁺-ATPase, sarcoplasmic reticular Ca²⁺-ATPase, background K⁺ channel and cholinergic receptor function. Age independently decreased Na⁺-K⁺-ATPase and fibrotic markers. Both factors interacted to alter Hcn4 channel activity underlying atrial automaticity. However, neither factor, whether independently or interactively, affected transcription of cardiac Na⁺, voltage-dependent K⁺ channels, surface or intracellular Ca²⁺ channels. Nor were gap junction channels, β-adrenergic receptors or transforming growth factor-β affected.

CONCLUSION

These findings limit the possible roles of gene transcriptional changes in previously reported age-dependent pro-arrhythmic electrophysiologial changes observed in Pgc-1β^-/- atria to an altered Ca²⁺-ATPase (Atp2a2) expression. This directly parallels previously reported arrhythmic mechanism associated with p21-activated kinase type 1 deficiency. This could add to contributions from the direct physiological outcomes of mitochondrial dysfunction, whether through reactive oxygen species (ROS) production or altered Ca²⁺ homeostasis.

Reduced cardiomyocyte Na+ current in the age-dependent murine Pgc-1β-/- model of ventricular arrhythmia

Peroxisome proliferator-activated receptor-γ coactivator-1 deficient (Pgc-1β-/- ) murine hearts model the increased, age-dependent, ventricular arrhythmic risks attributed to clinical conditions associated with mitochondrial energetic dysfunction. These were accompanied by compromised action potential (AP) upstroke rates and impaired conduction velocities potentially producing arrhythmic substrate. We tested a hypothesis implicating compromised Na+ current in these electrophysiological phenotypes by applying loose patch-clamp techniques in intact young and aged, wild-type (WT) and Pgc-1β-/- , ventricular cardiomyocyte preparations for the first time. This allowed conservation of their in vivo extracellular and intracellular conditions. Depolarising steps elicited typical voltage-dependent activating and inactivating inward Na+ currents with peak amplitudes increasing or decreasing with their respective activating or preceding inactivating voltage steps. Two-way analysis of variance associated Pgc-1β-/- genotype with independent reductions in maximum peak ventricular Na+ currents from -36.63 ± 2.14 (n = 20) and -35.43 ± 1.96 (n = 18; young and aged WT, respectively), to -29.06 ± 1.65 (n = 23) and -27.93 ± 1.63 (n = 20; young and aged Pgc-1β-/- , respectively) pA/μm2 (p < 0.0001), without independent effects of, or interactions with age. Voltages at half-maximal current V*, and steepness factors k in plots of voltage dependences of both Na+ current activation and inactivation, and time constants for its postrepolarisation recovery from inactivation, remained indistinguishable through all experimental groups. So were the activation and rectification properties of delayed outward (K+ ) currents, demonstrated from tail currents reflecting current recoveries from respective varying or constant voltage steps. These current-voltage properties directly implicate decreases specifically in maximum available Na+ current with unchanged voltage dependences and unaltered K+ current properties, in proarrhythmic reductions in AP conduction velocity in Pgc-1β-/- ventricles.

Left atrial scarring and conduction velocity dynamics: Rate dependent conduction slowing predicts sites of localized reentrant atrial tachycardias

BACKGROUND

Low voltage zones (LVZs) are associated with conduction velocity (CV) slowing. Rate-dependent CV slowing may play a role in reentry mechanisms.

METHODS

Patients undergoing catheter ablation for AT were enrolled. Aim was to assess the relationship between rate-dependent CV slowing and sites of localized reentrant atrial tachycardias (AT). On a bipolar voltage map regions were defined as non-LVZs [≥0.5 mV], LVZs [0.2-0.5 mV] and very-LVZs [<0.2 mV]. Unipolar electrograms were recorded with a 64-pole basket catheter during uninterrupted atrial pacing at four pacing intervals (PIs) during sinus rhythm. CVs were measured between pole pairs along the wavefront path. Sites of rate-dependent CV slowing were defined as exhibiting a reduction in CV between PI = 600 ms and 250 ms of ≥20% more than the mean CV reduction seen between these PIs for that voltage zone. Rate-dependent CV slowing sites were correlated to sites of localized reentrant ATs as confirmed with conventional mapping, entrainment and response to ablation.

RESULTS

Eighteen patients were included (63 ± 10 years). Mean CV at 600 ms was 1.53 ± 0.19 m/s in non-LVZs, 1.14 ± 0.15 m/s in LVZs, and 0.73 ± 0.13 m/s in very-LVZs respectively (p < 0.001). Rate-dependent CV slowing sites were predominantly in LVZs [0.2-0.5 mV] (74.4 ± 10.3%; p < 0.001). Localized reentrant ATs were mapped to these sites in 81.8% of cases (sensitivity 81.8%, 95% CI 48.2-97.9% and specificity 83.9%, 95% CI 81.8-86.0%). Macro-reentrant or focal ATs were not mapped to sites of rate-dependent CV slowing.

CONCLUSIONS

Rate-dependent CV slowing sites are predominantly confined to LVZs [0.2-0.5 mV] and the resultant CV heterogeneity may promote reentry mechanisms. These may represent a novel adjunctive target for AT ablation.

Ca2+ Cycling Impairment in Heart Failure Is Exacerbated by Fibrosis: Insights Gained From Mechanistic Simulations

Heart failure (HF) is characterized by altered Ca²⁺ cycling, resulting in cardiac contractile dysfunction. Failing myocytes undergo electrophysiological remodeling, which is known to be the main cause of abnormal Ca²⁺ homeostasis. However, structural remodeling, specifically proliferating fibroblasts coupled to myocytes in the failing heart, could also contribute to Ca²⁺ cycling impairment. The goal of the present study was to systematically analyze the mechanisms by which myocyte–fibroblast coupling could affect Ca²⁺ dynamics in normal conditions and in HF. Simulations of healthy and failing human myocytes were performed using established mathematical models, and cells were either isolated or coupled to fibroblasts. Univariate and multivariate sensitivity analyses were performed to quantify effects of ion transport pathways on biomarkers computed from intracellular [Ca²⁺] waveforms. Variability in ion channels and pumps was imposed and populations of models were analyzed to determine effects on Ca²⁺ dynamics. Our results suggest that both univariate and multivariate sensitivity analyses are valuable methodologies to shed light into the ionic mechanisms underlying Ca²⁺ impairment in HF, although differences between the two methodologies are observed at high parameter variability. These can result from either the fact that multivariate analyses take into account ion channels or non-linear effects of ion transport pathways on Ca²⁺ dynamics. Coupling either healthy or failing myocytes to fibroblasts decreased Ca²⁺ transients due to an indirect sink effect on action potential (AP) and thus on Ca²⁺ related currents. Simulations that investigated restoration of normal physiology in failing myocytes showed that Ca²⁺ cycling can be normalized by increasing SERCA and L-type Ca²⁺ current activity while decreasing Na⁺–Ca²⁺ exchange and SR Ca²⁺ leak. Changes required to normalize APs in failing myocytes depended on whether myocytes were coupled to fibroblasts. In conclusion, univariate and multivariate sensitivity analyses are helpful tools to understand how Ca²⁺ cycling is impaired in HF and how this can be exacerbated by coupling of myocytes to fibroblasts. The design of pharmacological actions to restore normal activity should take into account the degree of fibrosis in the failing heart.

Structural remodeling and conduction velocity dynamics in the human left atrium: Relationship with reentrant mechanisms sustaining atrial fibrillation

Background

Rate-dependent conduction velocity (CV) slowing is associated with atrial fibrillation (AF) initiation and reentrant mechanisms.

Objective

The purpose of this study was to assess the relationship between bipolar voltage, CV dynamics, and AF drivers.

Methods

Patients undergoing catheter ablation for persistent AF (<24 months) were enrolled. Unipolar electrograms were recorded with a 64-pole basket catheter during atrial pacing at 4 pacing intervals (PIs) during sinus rhythm. CVs were measured between pole pairs along the wavefront path and correlated with underlying bipolar voltage. CV dynamics within low-voltage zones (LVZs <0.5 mV) were compared to those of non-LVZs (≥0.5 mV) and were correlated to driver sites mapped using CARTOFINDER (Biosense Webster).

Results

Eighteen patients were included (age 62 ± 10 years). Mean CV at 600 ms was 1.59 ± 0.13 m/s in non-LVZs vs 0.98 ± 0.23 m/s in LVZs (P <.001). CV decreased incrementally over all 4 PIs in LVZs, whereas in non-LVZs a substantial decrease in CV was only seen between PIs 300–250 ms (0.59 ± 0.09 m/s; P <.001). Rate-dependent CV slowing sites measurements, defined as exhibiting CV reduction ≥20% more than the mean CV reduction seen between PIs 600–250 ms for that voltage zone, were predominantly in LVZs (0.2–0.5 mV; 75.6% ± 15.5%; P <.001). Confirmed rotational drivers were mapped to these sites in 94.1% of cases (sensitivity 94.1%, 95% CI 71.3%–99.9%; specificity 77.9%, 95% CI 74.9%–80.7%).

Conclusion

CV dynamics are determined largely by the extent of remodeling. Rate-dependent CV slowing sites are predominantly confined to LVZs (0.2–0.5 mV), and the resultant CV heterogeneity may promote driver formation in AF.

Quantitative systems models illuminate arrhythmia mechanisms in heart failure: Role of the Na+ -Ca2+ -Ca2+ /calmodulin-dependent protein kinase II-reactive oxygen species feedback

Quantitative systems modeling aims to integrate knowledge in different research areas with models describing biological mechanisms and dynamics to gain a better understanding of complex clinical syndromes. Heart failure (HF) is a chronic complex cardiac disease that results from structural or functional disorders impairing the ability of the ventricle to fill with or eject blood. Highly interactive and dynamic changes in mechanical, structural, neurohumoral, metabolic, and electrophysiological properties collectively predispose the failing heart to cardiac arrhythmias, which are responsible for about a half of HF deaths. Multiscale cardiac modeling and simulation integrate structural and functional data from HF experimental models and patients to improve our mechanistic understanding of this complex arrhythmia syndrome. In particular, they allow investigating how disease-induced remodeling alters the coupling of electrophysiology, Ca²⁺ and Na⁺ handling, contraction, and energetics that lead to rhythm derangements. The Ca²⁺ /calmodulin-dependent protein kinase II, which expression and activity are enhanced in HF, emerges as a critical hub that modulates the feedbacks between these various subsystems and promotes arrhythmogenesis. This article is categorized under: Physiology > Mammalian Physiology in Health and Disease Models of Systems Properties and Processes > Mechanistic Models Models of Systems Properties and Processes > Cellular Models Models of Systems Properties and Processes > Organ, Tissue, and Physiological Models.

Ventricular pro-arrhythmic phenotype, arrhythmic substrate, ageing and mitochondrial dysfunction in peroxisome proliferator activated receptor-γ coactivator-1β deficient (Pgc-1β-/-) murine hearts

INTRODUCTION

Ageing and age-related bioenergetic conditions including obesity, diabetes mellitus and heart failure constitute clinical ventricular arrhythmic risk factors.

MATERIALS AND METHODS

Pro-arrhythmic properties in electrocardiographic and intracellular recordings were compared in young and aged, peroxisome proliferator-activated receptor-γ coactivator-1β knockout (Pgc-1β-/-) and wild type (WT), Langendorff-perfused murine hearts, during regular and programmed stimulation (PES), comparing results by two-way ANOVA.

RESULTS AND DISCUSSION

Young and aged Pgc-1β-/- showed higher frequencies and durations of arrhythmic episodes through wider PES coupling-interval ranges than WT. Both young and old, regularly-paced, Pgc-1β-/- hearts showed slowed maximum action potential (AP) upstrokes, (dV/dt)max (∼157 vs. 120-130 V s-1), prolonged AP latencies (by ∼20%) and shortened refractory periods (∼58 vs. 51 ms) but similar AP durations (∼50 ms at 90% recovery) compared to WT. However, Pgc-1β-/- genotype and age each influenced extrasystolic AP latencies during PES. Young and aged WT ventricles displayed distinct, but Pgc-1β-/- ventricles displayed similar dependences of AP latency upon (dV/dt)max resembling aged WT. They also independently increased myocardial fibrosis. AP wavelengths combining activation and recovery terms paralleled contrasting arrhythmic incidences in Pgc-1β-/- and WT hearts. Mitochondrial dysfunction thus causes pro-arrhythmic Pgc-1β-/- phenotypes by altering AP conduction through reducing (dV/dt)max and causing age-dependent fibrotic change.

Killing Many Birds With Two Stones: Hypoxia and Fibrosis Can Generate Ectopic Beats in a Human Ventricular Model

During cardiac diseases many types of anatomical and functional remodeling of cardiac tissue can occur. In this work, we focus on two conditions: hypoxia and fibrosis, which are part of complex pathological modifications that take place in many cardiac diseases (hypertrophic cardiomyopathy, hypertensive heart disease, and recurrent myocardial infarction) and respiratory diseases (obstructive pulmonary disease, obstructive sleep apnea, and cystic fibrosis). Using computational models of cardiac electrophysiology, we evaluate if the interplay between hypoxia and fibrosis is sufficient to trigger cardiac arrhythmia. We study the mechanisms behind the generation of ectopic beats, an arrhythmic trigger also known as premature ventricular contractions (PVCs), in regions with high hypoxia and fibrosis. First, we modify an electrophysiological model of myocytes of the human left ventricle to include the effects of hypoxia. Second, diffuse fibrosis is modeled by randomly replacing cardiac myocytes by non-excitable and non-conducting cells. The Monte Carlo method is used to evaluate the probability of a region to generate ectopic beats with respect to different levels of hypoxia and fibrosis. In addition, we evaluate the minimum size of three-dimensional slabs needed to sustain reentries for different stimulation protocols. The observed mechanism behind the initiation of ectopic beats is unidirectional block, giving rise to sustained micro-reentries inside the region with diffuse fibrosis and hypoxia. In summary, our results suggest that hypoxia and fibrosis are sufficient for the creation of a focal region in the heart that generates PVCs.

Directed fusion of cardiac spheroids into larger heterocellular microtissues enables investigation of cardiac action potential propagation via cardiac fibroblasts

Multicellular spheroids generated through cellular self-assembly provide cytoarchitectural complexities of native tissue including three-dimensionality, extensive cell-cell contacts, and appropriate cell-extracellular matrix interactions. They are increasingly suggested as building blocks for larger engineered tissues to achieve shapes, organization, heterogeneity, and other biomimetic complexities. Application of these tissue culture platforms is of particular importance in cardiac research as the myocardium is comprised of distinct but intermingled cell types. Here, we generated scaffold-free 3D cardiac microtissue spheroids comprised of cardiac myocytes (CMs) and/or cardiac fibroblasts (CFs) and used them as building blocks to form larger microtissues with different spatial distributions of CMs and CFs. Characterization of fusing homotypic and heterotypic spheroid pairs revealed an important influence of CFs on fusion kinetics, but most strikingly showed rapid fusion kinetics between heterotypic pairs consisting of one CF and one CM spheroid, indicating that CMs and CFs self-sort in vitro into the intermixed morphology found in the healthy myocardium. We then examined electrophysiological integration of fused homotypic and heterotypic microtissues by mapping action potential propagation. Heterocellular elongated microtissues which recapitulate the disproportionate CF spatial distribution seen in the infarcted myocardium showed that action potentials propagate through CF volumes albeit with significant delay. Complementary computational modeling revealed an important role of CF sodium currents and the spatial distribution of the CM-CF boundary in action potential conduction through CF volumes. Taken together, this study provides useful insights for the development of complex, heterocellular engineered 3D tissue constructs and their engraftment via tissue fusion and has implications for arrhythmogenesis in cardiac disease and repair.

Effects of fibroblast-myocyte coupling on the sinoatrial node activity: A computational study

While the sinoatrial node (SAN) is structurally heterogeneous, most computer simulations of electrical activity take into account SAN pacemaker cells only. Our aim was to investigate how fibroblasts affect the SAN activity. We simulated the rabbit sinoatrial node accounting for differences between central and peripheral pacemaker cells, and for fibroblast-myocyte electrical coupling. We have observed that only if fibroblast-myocyte coupling is taken into account, (1) action potential is initiated in the central part of the SAN (within 1.2 mm of the center of simulated tissue); otherwise, leading centers are located on the periphery; (2) few (1 to 6) leading centers initiate action potential in the SAN; otherwise, we observed more than 8 leading centers; (3) acetylcholine superfusion results in a shift of leading centers toward the SAN periphery; and (4) sinus pauses up to 1.9 second follow acetylcholine superfusion. We observed negligible effect of fibroblast-myocyte coupling on the period of SAN activation. We conclude that fibroblast-myocyte coupling may explain action potential initiation and propagation from the center of the SAN observed in experimental studies, while atrial load on the peripheral SAN fails to explain this fact.

Cardiac fibroblasts : Active players in (atrial) electrophysiology?

Fibrotic areas in cardiac muscle-be it in ventricular or atrial tissue-are considered as obstacles for conduction of the excitatory wave and can therefore facilitate re-entry, which may contribute to the sustenance of cardiac arrhythmias. Persistence of one of the most frequent arrhythmias, atrial fibrillation (AF), is accompanied by enhanced atrial fibrosis. Any kind of myocardial perturbation, whether via mechanical stress or ischemic damage, inflammation, or irregular and high-frequency electrical activity, activates fibroblasts. This leads to the secretion of paracrine factors and extracellular matrix proteins, especially collagen, and to the differentiation of fibroblasts into myofibroblasts. Excessive collagen production is the hallmark of fibrosis and impairs regular impulse propagation. In addition, direct electrical coupling between cardiomyocytes and nonmyocytes, such as fibroblasts and macrophages, via gap junctions affects conduction. Although fibroblasts are not electrically excitable, they express functional ion channels, in particular K⁺ channels and mechanosensitive channels, some of which could be involved in tissue remodeling. Here, we briefly review these aspects with special reference to AF.

Cardiac ischemia-insights from computational models

BACKGROUND

Complementary to clinical and experimental studies, computational cardiac modeling serves to obtain a comprehensive understanding of the cardiovascular system in order to analyze dysfunction, evaluate existing, and develop novel treatment strategies.

OBJECTIVES

We describe the basics of multiscale computational modeling of cardiac electrophysiology from the molecular ion channel to the whole body scale. By modeling cardiac ischemia, we illustrate how in silico experiments can contribute to our understanding of how the pathophysiological mechanisms translate into changes observed in diagnostic tools such as the electrocardiogram (ECG).

MATERIALS AND METHODS

Quantitative in silico modeling spans a wide range of scales from ion channel biophysics to ECG signals. For each of the scales, a set of mathematical equations describes electrophysiology in relation to the other scales. Integration of ischemia-induced changes is performed on the ion channel, single-cell, and tissue level. This approach allows us to study how effects simulated at molecular scales translate to changes in the ECG.

RESULTS

Ischemia induces action potential shortening and conduction slowing. Hence, ischemic myocardium has distinct and significant effects on propagation and repolarization of excitation, depending on the intramural extent of the ischemic region. For transmural and subendocardial ischemic regions, ST segment elevation and depression, respectively, were observed, whereas intermediate ischemic regions were found to be electrically silent (NSTEMI).

CONCLUSIONS

In silico modeling contributes quantitative and mechanistic insight into fundamental ischemia-related arrhythmogenic mechanisms. In addition, computational modeling can help to translate experimental findings at the (sub-)cellular level to the organ and body context (e. g., ECG), thereby providing a thorough understanding of this routinely used diagnostic tool that may translate into optimized applications.

Age-dependent electrocardiographic changes in Pgc-1β deficient murine hearts

Increasing evidence implicates chronic energetic dysfunction in human cardiac arrhythmias. Mitochondrial impairment through Pgc-1β knockout is known to produce a murine arrhythmic phenotype. However, the cumulative effect of this with advancing age and its electrocardiographic basis have not been previously studied. Young (12-16 weeks) and aged (>52 weeks), wild type (WT) (n = 5 and 8) and Pgc-1β^-/- (n = 9 and 6), mice were anaesthetised and used for electrocardiographic (ECG) recordings. Time intervals separating successive ECG deflections were analysed for differences between groups before and after β1-adrenergic (intraperitoneal dobutamine 3 mg/kg) challenge. Heart rates before dobutamine challenge were indistinguishable between groups. The Pgc-1β^-/- genotype however displayed compromised nodal function in response to adrenergic challenge. This manifested as an impaired heart rate response suggesting a functional defect at the level of the sino-atrial node, and a negative dromotropic response suggesting an atrioventricular conduction defect. Incidences of the latter were most pronounced in the aged Pgc-1β^-/- mice. Moreover, Pgc-1β^-/- mice displayed electrocardiographic features consistent with the existence of a pro-arrhythmic substrate. Firstly, ventricular activation was prolonged in these mice consistent with slowed action potential conduction and is reported here for the first time. Additionally, Pgc-1β^-/- mice had shorter repolarisation intervals. These were likely attributable to altered K⁺ conductance properties, ultimately resulting in a shortened QT_c interval, which is also known to be associated with increased arrhythmic risk. ECG analysis thus yielded electrophysiological findings bearing on potential arrhythmogenicity in intact Pgc-1β^-/- systems in widespread cardiac regions.

Densification of Type I Collagen Matrices as a Model for Cardiac Fibrosis

Cardiac fibrosis is a disease state characterized by excessive collagenous matrix accumulation within the myocardium that can lead to ventricular dilation and systolic failure. Current treatment options are severely lacking due in part to the poor understanding of the complexity of molecular pathways involved in cardiac fibrosis. To close this gap, in vitro model systems that recapitulate the defining features of the fibrotic cellular environment are in need. Type I collagen, a major cardiac extracellular matrix protein and the defining component of fibrotic depositions, is an attractive choice for a fibrosis model, but demonstrates poor mechanical strength due to solubility limits. However, plastic compression of collagen matrices is shown to significantly increase its mechanical properties. Here, confined compression of oligomeric, type I collagen matrices is utilized to resemble defining hallmarks seen in fibrotic tissue such as increased collagen content, fibril thickness, and bulk compressive modulus. Cardiomyocytes seeded on compressed matrices show a strong beating abrogation as observed in cardiac fibrosis. Gene expression analysis of selected fibrosis markers indicates fibrotic activation and cardiomyocyte maturation with regard to the existing literature. With these results, a promising first step toward a facile heart-on-chip model is presented to study cardiac fibrosis.

Age-dependent atrial arrhythmic phenotype secondary to mitochondrial dysfunction in Pgc-1β deficient murine hearts

INTRODUCTION

Ageing and several age-related chronic conditions including obesity, insulin resistance and hypertension are associated with mitochondrial dysfunction and represent independent risk factors for atrial fibrillation (AF).

MATERIALS AND METHODS

Atrial arrhythmogenesis was investigated in Langendorff-perfused young (3-4 month) and aged (>12 month), wild type (WT) and peroxisome proliferator activated receptor-γ coactivator-1β deficient (Pgc-1β-/-) murine hearts modeling age-dependent chronic mitochondrial dysfunction during regular pacing and programmed electrical stimulation (PES).

RESULTS AND DISCUSSION

The Pgc-1β-/- genotype was associated with a pro-arrhythmic phenotype progressing with age. Young and aged Pgc-1β-/- hearts showed compromised maximum action potential (AP) depolarization rates, (dV/dt)max, prolonged AP latencies reflecting slowed action potential (AP) conduction, similar effective refractory periods and baseline action potential durations (APD90) but shortened APD90 in APs in response to extrasystolic stimuli at short stimulation intervals. Electrical properties of APs triggering arrhythmia were similar in WT and Pgc-1β-/- hearts. Pgc-1β-/- hearts showed accelerated age-dependent fibrotic change relative to WT, with young Pgc-1β-/- hearts displaying similar fibrotic change as aged WT, and aged Pgc-1β-/- hearts the greatest fibrotic change. Mitochondrial deficits thus result in an arrhythmic substrate, through slowed AP conduction and altered repolarisation characteristics, arising from alterations in electrophysiological properties and accelerated structural change.

Spatiotemporally Non-Uniform Ca2+ Dynamics of Cardiac Purkinje Fibers in Mouse Myocardial Infarct

Surviving Purkinje fibers in myocardial infarct are regarded as an important substrate in arrhythmogenesis. However, poorly understood are functional properties of Purkinje fibers in the infarcted heart. We sought to visualize intracellular Ca²⁺ ([Ca²⁺]_i) dynamics of Purkinje fiber networks in the mouse myocardial infarct. Using 3- to 4-day-old or 7- to 9-day-old infarcted hearts after the left coronary-artery ligation corresponding, respectively, to acute or healing phase, we conducted rapid fluo4-fluorescence imaging on the endocardial surface of the left ventricular septum by macro-zoom fluorescence microscopy and rapid-scanning confocal microscopy. In contrast with the intact heart, where uniform Ca²⁺ transients propagated rapidly, the infarcted heart exhibited slow, non-uniform impulse propagations. On confocal microscopy, Purkinje fibers in the peri-infarct zone exhibited non-uniform [Ca²⁺]_i dynamics: beat-to-beat alternans of the Ca²⁺ transient amplitude in and among the individual fibers, whereas the intact fibers exhibited uniform Ca²⁺ transients. Such non-uniform [Ca²⁺]_i dynamics were more conspicuous in the acute infarcted hearts than in the healing ones. In accordance with [Ca²⁺]_i dynamics, fixed fluo4-loaded heart preparations exhibited definitive connexin-40 plaques in the peri-infarct Purkinje fibers, whereas the subjacent myocardium presented coagulative necrosis and granulation tissues, respectively. The surviving Purkinje fibers in the peri-infarct zone exhibited non-uniform [Ca²⁺]_i dynamics, which may lead to arrhythmogenesis.

The effects of ageing and adrenergic challenge on electrocardiographic phenotypes in a murine model of long QT syndrome type 3

Long QT Syndrome 3 (LQTS3) arises from gain-of-function Na_v1.5 mutations, prolonging action potential repolarisation and electrocardiographic (ECG) QT interval, associated with increased age-dependent risk for major arrhythmic events, and paradoxical responses to β-adrenergic agents. We investigated for independent and interacting effects of age and Scn5a+/ΔKPQ genotype in anaesthetised mice modelling LQTS3 on ECG phenotypes before and following β-agonist challenge, and upon fibrotic change. Prolonged ventricular recovery was independently associated with Scn5a+/ΔKPQ and age. Ventricular activation was prolonged in old Scn5a+/ΔKPQ mice (p = 0.03). We associated Scn5a+/ΔKPQ with increased atrial and ventricular fibrosis (both: p < 0.001). Ventricles also showed increased fibrosis with age (p < 0.001). Age and Scn5a+/ΔKPQ interacted in increasing incidences of repolarisation alternans (p = 0.02). Dobutamine increased ventricular rate (p < 0.001) and reduced both atrioventricular conduction (PR segment-p = 0.02; PR interval-p = 0.02) and incidences of repolarisation alternans (p < 0.001) in all mice. However, in Scn5a+/ΔKPQ mice, dobutamine delayed the changes in ventricular repolarisation following corresponding increases in ventricular rate. The present findings implicate interactions between age and Scn5a+/ΔKPQ in prolonging ventricular activation, correlating them with fibrotic change for the first time, adding activation abnormalities to established recovery abnormalities in LQTS3. These findings, together with dynamic electrophysiological responses to β-adrenergic challenge, have therapeutic implications for ageing LQTS patients.