Gap junction heterogeneity as mechanism for electrophysiologically distinct properties across the ventricular wall

Conductance heterogeneities induced by multistability in the dynamics of coupled cardiac gap junctions

In this paper, we study the propagation of the cardiac action potential in a one-dimensional fiber, where cells are electrically coupled through gap junctions (GJs). We consider gap junctional gate dynamics that depend on the intercellular potential. We find that different GJs in the tissue can end up in two different states: a low conducting state and a high conducting state. We first present evidence of the dynamical multistability that occurs by setting specific parameters of the GJ dynamics. Subsequently, we explain how the multistability is a direct consequence of the GJ stability problem by reducing the dynamical system's dimensions. The conductance dispersion usually occurs on a large time scale, i.e., thousands of heartbeats. The full cardiac model simulations are computationally demanding, and we derive a simplified model that allows for a reduction in the computational cost of four orders of magnitude. This simplified model reproduces nearly quantitatively the results provided by the original full model. We explain the discrepancies between the two models due to the simplified model's lack of spatial correlations. This simplified model provides a valuable tool to explore cardiac dynamics over very long time scales. That is highly relevant in studying diseases that develop on a large time scale compared to the basic heartbeat. As in the brain, plasticity and tissue remodeling are crucial parameters in determining the action potential wave propagation's stability.

Dispersion of ventricular repolarization: Temporal and spatial

Repolarization heterogeneity (RH) is an intrinsic property of ventricular myocardium and the reason for T-wave formation on electrocardiogram (ECG). Exceeding the physiologically based RH level is associated with appearance of life-threatening ventricular arrhythmias and sudden cardiac death. In this regard, an accurate and comprehensive evaluation of the degree of RH parameters is of importance for assessment of heart state and arrhythmic risk. This review is devoted to comprehensive consideration of RH phenomena in terms of electrophysiological processes underlying RH, cardiac electric field formation during ventricular repolarization, as well as clinical significance of RH and its reflection on ECG parameters. The formation of transmural, apicobasal, left-to-right and anterior-posterior gradients of action potential durations and end of repolarization times resulting from the heterogenous distribution of repolarizing ion currents and action potential morphology throughout the heart ventricles, and the different sensitivity of myocardial cells in different ventricular regions to the action of pharmacological agents, temperature, frequency of stimulation, etc., are being discussed. The review is focused on the fact that RH has different aspects – temporal and spatial, global and local; ECG reflection of various RH aspects and their clinical significance are being discussed. Strategies for comprehensive assessment of ventricular RH using different ECG indices reflecting various RH aspects are presented.

Mechanisms of Arrhythmias in the Brugada Syndrome

Arrhythmias in Brugada syndrome patients originate in the right ventricular outflow tract (RVOT). Over the past few decades, the characterization of the unique anatomy and electrophysiology of the RVOT has revealed the arrhythmogenic nature of this region. However, the mechanisms that drive arrhythmias in Brugada syndrome patients remain debated as well as the exact site of their occurrence in the RVOT. Identifying the site of origin and mechanism of Brugada syndrome would greatly benefit the development of mechanism-driven treatment strategies.

Effect of myocardial heterogeneity on ventricular electro-mechanical responses: a computational study

BACKGROUND

The heart wall exhibits three layers of different thicknesses: the outer epicardium, mid-myocardium, and inner endocardium. Among these layers, the mid-myocardium is typically the thickest. As indicated by preliminary studies, heart-wall layers exhibit various characteristics with regard to electrophysiology, pharmacology, and pathology. Construction of an accurate three-dimensional (3D) model of the heart is important for predicting physiological behaviors. However, the wide variability of myocardial shapes and the unclear edges between the epicardium and soft tissues are major challenges in the 3D model segmentation approach for identifying the boundaries of the epicardium, mid-myocardium, and endocardium. Therefore, this results in possible variations in the heterogeneity ratios between the epicardium, mid-myocardium, and endocardium. The objective of this study was to observe the effects of different thickness ratios of the epicardium, mid-myocardium, and endocardium on cardiac arrhythmogenesis, reentry instability, and mechanical responses during arrhythmia.

METHODS

We used a computational method and simulated three heterogeneous ventricular models: Model 1 had the thickest M cell layer and thinnest epicardium and endocardium. Model 2 had intermediate layer thicknesses. Model 3 exhibited the thinnest mid-myocardium and thickest epicardium and endocardium. Electrical and mechanical simulations of the three heterogeneous models were performed under normal sinus rhythm and reentry conditions.

RESULTS

Model 1 exhibited the highest probability of terminating reentrant waves, and Model 3 exhibited to experience greater cardiac arrhythmia. In the reentry simulation, at 8 s, Model 3 generated the largest number of rotors (eight), while Models 1 and 2 produced five and seven rotors, respectively. There was no significant difference in the cardiac output obtained during the sinus rhythm. Under the reentry condition, the highest cardiac output was generated by Model 1 (19 mL/s), followed by Model 2 (9 mL/s) and Model 3 (7 mL/s).

CONCLUSIONS

A thicker mid-myocardium led to improvements in the pumping efficacy and contractility and reduced the probability of cardiac arrhythmia. Conversely, thinner M cell layers generated more unstable reentrant spiral waves and hindered the ventricular pumping.

Trigger vs. Substrate: Multi-Dimensional Modulation of QT-Prolongation Associated Arrhythmic Dynamics by a hERG Channel Activator

Background: Prolongation of the QT interval of the electrocardiogram (ECG), underlain by prolongation of the action potential duration (APD) at the cellular level, is linked to increased vulnerability to cardiac arrhythmia. Pharmacological management of arrhythmia associated with QT prolongation is typically achieved through attempting to restore APD to control ranges, reversing the enhanced vulnerability to Ca²⁺-dependent afterdepolarisations (arrhythmia triggers) and increased transmural dispersion of repolarisation (arrhythmia substrate) associated with APD prolongation. However, such pharmacological modulation has been demonstrated to have limited effectiveness. Understanding the integrative functional impact of pharmacological modulation requires simultaneous investigation of both the trigger and substrate. Methods: We implemented a multi-scale (cell and tissue) in silico approach using a model of the human ventricular action potential, integrated with a model of stochastic 3D spatiotemporal Ca²⁺ dynamics, and parameter modification to mimic prolonged QT conditions. We used these models to examine the efficacy of the hERG activator MC-II-157c in restoring APD to control ranges, examined its effects on arrhythmia triggers and substrates, and the interaction of these arrhythmia triggers and substrates. Results: QT prolongation conditions promoted the development of spontaneous release events underlying afterdepolarisations during rapid pacing. MC-II-157c applied to prolonged QT conditions shortened the APD, inhibited the development of afterdepolarisations and reduced the probability of afterdepolarisations manifesting as triggered activity in single cells. In tissue, QT prolongation resulted in an increased transmural dispersion of repolarisation, which manifested as an increased vulnerable window for uni-directional conduction block. In some cases, MC-II-157c further increased the vulnerable window through its effects on I_Na. The combination of stochastic release event modulation and transmural dispersion of repolarisation modulation by MC-II-157c resulted in an integrative behavior wherein the arrhythmia trigger is reduced but the arrhythmia substrate is increased, leading to variable and non-linear overall vulnerability to arrhythmia. Conclusion: The relative balance of reduced trigger and increased substrate underlies a multi-dimensional role of MC-II-157c in modulation of cardiac arrhythmia vulnerability associated with prolonged QT interval.

Gq-activated fibroblasts induce cardiomyocyte action potential prolongation and automaticity in a three-dimensional microtissue environment

Cardiac fibroblasts (CFs) are known to regulate cardiomyocyte (CM) function in vivo and in two-dimensional in vitro cultures. This study examined the effect of CF activation on the regulation of CM electrical activity in a three-dimensional (3-D) microtissue environment. Using a scaffold-free 3-D platform with interspersed neonatal rat ventricular CMs and CFs, G_q-mediated signaling was selectively enhanced in CFs by Gα_q adenoviral infection before coseeding with CMs in nonadhesive hydrogels. After 3 days, the microtissues were analyzed by signaling assay, histological staining, quantitative PCR, Western blots, optical mapping with voltage- or Ca²⁺-sensitive dyes, and microelectrode recordings of CF resting membrane potential (RMP_CF). Enhanced G_q signaling in CFs increased microtissue size and profibrotic and prohypertrophic markers. Expression of constitutively active Gα_q in CFs prolonged CM action potential duration (by 33%) and rise time (by 31%), prolonged Ca²⁺ transient duration (by 98%) and rise time (by 65%), and caused abnormal electrical activity based on depolarization-induced automaticity. Constitutive G_q activation in CFs also depolarized RMP_CF from -33 to -20 mV and increased connexin 43 and connexin 45 expression. Computational modeling confers that elevated RMP_CF and increased cell-cell coupling between CMs and CFs in a 3-D environment could lead to automaticity. In conclusion, our data demonstrate that CF activation alone is capable of altering action potential and Ca²⁺ transient characteristics of CMs, leading to proarrhythmic electrical activity. Our results also emphasize the importance of a 3-D environment where cell-cell interactions are prevalent, underscoring that CF activation in 3-D tissue plays a significant role in modulating CM electrophysiology and arrhythmias.NEW & NOTEWORTHY In a three-dimensional microtissue model, which lowers baseline activation of cardiac fibroblasts but enables cell-cell, paracrine, and cell-extracellular matrix interactions, we demonstrate that selective cardiac fibroblast activation by enhanced G_q signaling, a pathophysiological trigger in the diseased heart, modulates cardiomyocyte electrical activity, leading to proarrhythmogenic automaticity.

Microinjection Technique for Assessment of Gap Junction Function

Gap junctions are essential for the proper function of many native mammalian tissues including neurons, cardiomyocytes, embryonic tissues, and muscle. Assessing these channels is therefore fundamental to understanding disease pathophysiology, developing therapies for a multitude of acquired and genetic conditions, and providing novel approaches to drug delivery and cellular communication. Microinjection is a robust, albeit difficult, technique, which provides considerable information that is superior to many of the simpler techniques due to its ability to isolate cells, quantify kinetics, and allow cross-comparison of multiple cell lines. Despite its user-dependent nature, the strengths of the technique are considerable and with the advent of new, automation technologies may improve further. This text describes the basic technique of microinjection and briefly discusses modern automation advances that can improve the success rates of this technique.

Passive ventricular remodeling in cardiac disease: focus on heterogeneity

Passive ventricular remodeling is defined by the process of molecular ventricular adaptation to different forms of cardiac pathophysiology. It includes changes in tissue architecture, such as hypertrophy, fiber disarray, alterations in cell size and fibrosis. Besides that, it also includes molecular remodeling of gap junctions, especially those composed by Connexin43 proteins (Cx43) in the ventricles that affect cell-to-cell propagation of the electrical impulse, and changes in the sodium channels that modify excitability. All those alterations appear mainly in a heterogeneous manner, creating irregular and inhomogeneous electrical and mechanical coupling throughout the heart. This can predispose to reentry arrhythmias and adds to a further deterioration into heart failure. In this review, passive ventricular remodeling is described in Hypertrophic Cardiomyopathy (HCM), Dilated Cardiomyopathy (DCM), Ischemic Cardiomyopathy (ICM), and Arrhythmogenic Cardiomyopathy (ACM), with a main focus on the heterogeneity of those alterations mentioned above.

Protective effect of Salvia miltiorrhiza aqueous extract on myocardium oxidative injury in ischemic-reperfusion rats

Salvia miltiorrhiza has strong antioxidative activity. They may have a strong potential as cardioprotective agents in ischemic-reperfusion injury. Experiments were carried out in Sprague-Dawley rats with myocardium ischemia reperfusion (IR). Myocardial injuries during IR were determined by changes in electrocardiogram analysis of arrhythmias, antioxidant enzyme activities, AST, CK-MB, lactate dehydrogenase (LDH) levels, and myocyte apoptosis. Results showed that S. miltiorrhiza aqueous extract (SAME) pre-treatment significantly decreased the ST-segment (ΣST120) and myocardium MDA, AST, CK-MB, lactate dehydrogenase (LDH) levels, increased myocardium antioxidant enzyme activities, and inhibit myocardium cell apoptosis. Furthermore, the SAME pre-treatment significantly upregulated p-JAK2 and p-STAT3 protein expression, decreased myocardium TNF-α and IL-6 concentrations in IR rats. The levels of TNF-α and IL-6 were positively correlated with the changes in myocardium p-JAK2 and p-STAT3 protein expression levels in IR rats. It can be concluded that the SAME pre-treatment has anti-ischemic and anti-apoptosis activity in heart IR rats. SAME pre-treatment protects heart against IR injury, at least in part, through its stimulating effects on injury-induced deactivation of JAK2/STAT3 signaling pathway.

High-throughput measurement of gap junctional intercellular communication

Gap junctional intercellular communication (GJIC) is a critical part of cellular activities and is necessary for electrical propagation among contacting cells. Disorders of gap junctions are a major cause for cardiac arrhythmias. Dye transfer through microinjection is a conventional technique for measuring GJIC. To overcome the limitations of manual microinjection and perform high-throughput GJIC measurement, here we present a new robotic microinjection system that is capable of injecting a large number of cells at a high speed. The highly automated system enables large-scale cell injection (thousands of cells vs. a few cells) without major operator training. GJIC of three cell lines of differing gap junction density, i.e., HeLa, HEK293, and HL-1, was evaluated. The effect of a GJIC inhibitor (18-α-glycyrrhetinic acid) was also quantified in the three cell lines. System operation speed, success rate, and cell viability rate were quantitatively evaluated based on robotic microinjection of over 4,000 cells. Injection speed was 22.7 cells per min, with 95% success for cell injection and >90% survival. Dye transfer cell counts and dye transfer distance correlated with the expected connexin expression of each cell type, and inhibition of dye transfer correlated with the concentration of GJIC inhibitor. Additionally, real-time monitoring of dye transfer enables the calculation of coefficients of molecular diffusion through gap junctions. This robotic microinjection dye transfer technique permits rapid assessment of gap junction function in confluent cell cultures.

Effects of deformation on transmural dispersion of repolarization using in silico models of human left ventricular wedge

Mechanical deformation affects the electrical activity of the heart through multiple feedback loops. The purpose of this work is to study the effect of deformation on transmural dispersion of repolarization and on surface electrograms using an in silico human ventricular wedge. To achieve this purpose, we developed a strongly coupled electromechanical cell model by coupling a human left ventricle electrophysiology model and an active contraction model reparameterized for human cells. This model was then embedded in tissue simulations on the basis of bidomain equations and nonlinear solid mechanics. The coupled model was used to evaluate effects of mechanical deformation on important features of repolarization and electrograms. Our results indicate an increase in the T-wave amplitude of the surface electrograms in simulations that account for the effects of cardiac deformation. This increased T-wave amplitude can be explained by changes to the coupling between neighboring myocytes, also known as electrotonic effect. The thickening of the ventricular wall during repolarization contributes to the decoupling of cells in the transmural direction, enhancing action potential heterogeneity and increasing both transmural repolarization dispersion and T-wave amplitude of surface electrograms. The simulations suggest that a considerable percentage of the T-wave amplitude (15%) may be related to cardiac deformation.

Repolarization gradients in the intact heart: transmural or apico-basal?

Controversies regarding the genesis of the T wave in the electrocardiogram and the role of midmural M cells in the intact heart include: In normal, intact canine and human hearts there is no significant transmural gradient in repolarization times. The T wave results primarily from apico-basal differences in repolarization times. Also, in the intact heart there is no midmural region of prolonged action potential duration. This contrasts with isolated preparations, such as the wedge preparation or myocardial slices or disaggregated myocytes in which M cells, with action potentials longer than those of endocardial and epicardial myocardium, can be found. This disparity in action potential duration probably results from partial uncoupling of myocardial cells in the regions where measurements are made, e.g., the cut surface of a wedge preparation. In regions of a wedge where cellular coupling is normal, or in isolated myocardial bundles or sheets, no evidence for M cells is detected. In some wedge preparations, a drug-induced large transmural repolarization gradient, involving M cells, can lead to Torsade de Pointes, possibly caused by so-called phase two reentry. In contrast, when a gradient of repolarization times of more than 100 ms was created in intact hearts, no evidence for reentry was found and no spontaneous arrhythmias occurred. In conclusion, in the intact heart, M cells appear not to contribute to repolarization gradients and arrhythmias. Furthermore, no significant repolarization gradients between endocardium and epicardium exist. The T wave in the body surface electrocardiogram is caused by apico-basal and anterior-posterior differences in repolarization times.

Selection of reference genes in different myocardial regions of an in vivo ischemia/reperfusion rat model for normalization of antioxidant gene expression

Background

Changes in cardiac gene expression due to myocardial injury are usually assessed in whole heart tissue. However, as the heart is a heterogeneous system, spatial and temporal heterogeneity is expected in gene expression.

Results

In an ischemia/reperfusion (I/R) rat model we evaluated gene expression of mitochondrial and cytoplasmatic superoxide dismutase (MnSod, Cu-ZnSod) and thioredoxin reductase (trxr1) upon short (4 h) and long (72 h) reperfusion times in the right ventricle (RV), and in the ischemic/reperfused (IRR) and the remote region (RR) of the left ventricle. Gene expression was assessed by Real-time reverse-transcription quantitative PCR (RT-qPCR). In order to select most stable reference genes suitable for normalization purposes, in each myocardial region we tested nine putative reference genes by geNorm analysis. The genes investigated were: Actin beta (actb), Glyceraldehyde-3-P-dehydrogenase (gapdh), Ribosomal protein L13A (rpl13a), Tyrosine 3-monooxygenase (ywhaz), Beta-glucuronidase (gusb), Hypoxanthine guanine Phosphoribosyltransferase 1 (hprt), TATA binding box protein (tbp), Hydroxymethylbilane synthase (hmbs), Polyadenylate-binding protein 1 (papbn1). According to our findings, most stable reference genes in the RV and RR were hmbs/hprt and hmbs/tbp/hprt respectively. In the IRR, six reference genes were recommended for normalization purposes; however, in view of experimental feasibility limitations, target gene expression could be normalized against the three most stable reference genes (ywhaz/pabp/hmbs) without loss of sensitivity. In all cases MnSod and Cu-ZnSod expression decreased upon long reperfusion, the former in all myocardial regions and the latter in IRR alone. trxr1 expression did not vary.

Conclusions

This study provides a validation of reference genes in the RV and in the anterior and posterior wall of the LV of cardiac ischemia/reperfusion model and shows that gene expression should be assessed separately in each region.

Cardiac electrical remodeling in health and disease

Electrical remodeling of the heart takes place in response to both functional (altered electrical activation) and structural (including heart failure and myocardial infarction) stressors. These electrophysiological changes produce a substrate that is prone to malignant ventricular arrhythmias. Understanding the cellular and molecular mechanisms of electrical remodeling is important in elucidating potential therapeutic targets designed to alter maladaptive electrical remodeling. For example, altered patterns of electrical activation lead primarily to electrical remodeling, without significant structural remodeling. By contrast, secondary remodeling arises in response to a structural insult. In this article we review cardiac electrical remodeling (predominantly in the ventricle) with an emphasis on the mechanisms causing these adaptations. These mechanisms suggest novel therapeutic targets for the management or prevention of the most devastating manifestation of heart disease, sudden cardiac death (SCD).

Transmural and apicobasal gradients in repolarization contribute to T-wave genesis in human surface ECG

The cellular basis of the T-wave morphology of surface ECG remains controversial in clinical cardiology. We examined the effect of action potential duration (APD) distribution on T-wave morphology using a realistic model of the human ventricle and torso. We developed a finite-element model of the ventricle consisting of ∼26 million elements, including the conduction system, each implemented with the ion current model of cardiomyocytes. This model was embedded in a torso model with distinct organ structures to obtain the standard ECG leads. The APD distribution was changed in the transmural direction by locating the M cells in either the endocardial or epicardial region. We also introduced apicobasal gradients by modifying the ion channel parameters. Both the transmural gradient (with M cells on the endocardial side) and the apicobasal gradient produced positive T waves, although a very large gradient was required for the apicobasal gradient. By contrast, T waves obtained with the transmural gradient were highly symmetric and, therefore, did not represent the true physiological state. Only combination of the transmural and the moderate apicobasal gradients produced physiological T waves in surface ECG. Positive T waves in surface ECG mainly originated from the transmural distribution of APD with M cells on the endocardial side, although the apicobasal gradient was also required to attain the physiological waveform.

Spatiotemporal induction of matrix metalloproteinase-9 transcription after discrete myocardial injury

Radiofrequency (RF) ablation of the myocardium causes discrete sites of injury. RF scars can expand, altering the extracellular matrix (ECM) structure and the continuity of the electrical syncytium of the adjacent myocardium. Matrix metalloproteinases (MMPs), such as MMP-9, contribute to ECM remodeling. However, whether and to what degree transcriptional induction of MMP-9 occurs after myocardial RF injury and the association with electrical conduction patterns after RF injury remains unexplored. This study examined MMP-9 gene promoter (M9PROM) activation after myocardial RF injury using mice in which the M9PROM was fused to a β-galactosidase (β-gal) reporter. RF lesions (0.5-mm probe, 80°C, 30 s) were created on the left ventricular (LV) epicardium of M9PROM mice (n=62) and terminally studied at 1 h, 1 d, 3 d, 7 d, 14 d, and 28 d after RF injury. M9PROM activation was localized through β-gal staining. The RF scar area and the area of β-gal staining were measured and normalized to LV area (planimetry). RF scar size increased from 1 h post-RF-injury values by 7 d and remained higher at 28 d. M9PROM activation became evident at 3 d and peaked at 7 d. Electrical conduction was measured (potentiometric dye mapping) at 7 d after RF injury. Heterogeneities in action potentials and electrical impulse propagation coincident with M9PROM activation were observed after RF injury. For example, conduction proximal to the RF site was slower than that in the remote myocardium (0.15±0.02 vs. 0.83±0.08 mm/ms, P<0.05). Thus, a unique spatiotemporal pattern of MMP-9 transcriptional activation occurred after discrete myocardial injury, which was associated with the development of electrical heterogeneity. Therefore, these findings suggest that changes in a key determinant of extracellular matrix remodeling, in addition to changes in myocardial structure, can contribute to arrhythmogenesis around the region of myocardial injury.