Increased interstitial loading reduces the effect of microstructural variations in cardiac tissue

Prokaryotic voltage-gated sodium channels are more effective than endogenous Nav1.5 channels in rescuing cardiac action potential conduction: an in silico study

Methods to augment Na+ current in cardiomyocytes hold potential for the treatment of various cardiac arrhythmias involving conduction slowing. Because the gene coding cardiac Na+ channel (Nav1.5) is too large to fit in a single adeno-associated virus (AAV) vector, new gene therapies are being developed to enhance endogenous Nav1.5 current (by overexpression of chaperon molecules or use of multiple AAV vectors) or to exogenously introduce prokaryotic voltage-gated Na+ channels (BacNav) whose gene size is significantly smaller than that of the Nav1.5. In this study, based on experimental measurements in heterologous expression systems, we developed an improved computational model of the BacNav channel, NavSheP D60A. We then compared in silico how NavSheP D60A expression vs. Nav1.5 augmentation affects the electrophysiology of cardiac tissue. We found that the incorporation of BacNav channels in both adult guinea pig and human cardiomyocyte models increased their excitability and reduced action potential duration. When compared with equivalent augmentation of Nav1.5 current in simulated settings of reduced tissue excitability, the addition of the BacNav current was superior in improving the safety of conduction under conditions of current source-load mismatch, reducing the vulnerability to unidirectional conduction block during premature pacing, preventing the instability and breakup of spiral waves, and normalizing the conduction and ECG in Brugada syndrome tissues with mutated Nav1.5. Overall, our studies show that compared with a potential enhancement of the endogenous Nav1.5 current, expression of the BacNav channels with their slower inactivation kinetics can provide greater anti-arrhythmic benefits in hearts with compromised action potential conduction.NEW & NOTEWORTHY Slow action potential conduction is a common cause of various cardiac arrhythmias; yet, current pharmacotherapies cannot augment cardiac conduction. This in silico study compared the efficacy of recently proposed antiarrhythmic gene therapy approaches that increase peak sodium current in cardiomyocytes. When compared with the augmentation of endogenous sodium current, expression of slower-inactivating bacterial sodium channels was superior in preventing conduction block and arrhythmia induction. These results further the promise of antiarrhythmic gene therapies targeting sodium channels.

Simulation of atrial fibrillation in a non-ohmic propagation model with dynamic gap junctions

Gap junctions exhibit nonlinear electrical properties that have been hypothesized to be relevant to arrhythmogenicity in a structurally remodeled tissue. Large-scale implementation of gap junction dynamics in 3D propagation models remains challenging. We aim to quantify the impact of nonlinear diffusion during episodes of arrhythmias simulated in a left atrial model. Homogenization of conduction properties in the presence of nonlinear gap junctions was performed by generalizing a previously developed mathematical framework. A monodomain model was solved in which conductivities were time-varying and depended on transjunctional potentials. Gap junction conductances were derived from a simplified Vogel-Weingart model with first-order gating and adjustable time constant. A bilayer interconnected cable model of the left atrium with 100  μm resolution was used. The diffusion matrix was recomputed at each time step according to the state of the gap junctions. Sinus rhythm and atrial fibrillation episodes were simulated in remodeled tissue substrates. Slow conduction was induced by reduced coupling and by diffuse or stringy fibrosis. Simulations starting from the same initial conditions were repeated with linear and nonlinear gap junctions. The discrepancy in activation times between the linear and nonlinear diffusion models was quantified. The results largely validated the linear approximation for conduction velocities >20 cm/s. In very slow conduction substrates, the discrepancy accumulated over time during atrial fibrillation, eventually leading to qualitative differences in propagation patterns, while keeping the descriptive statistics, such as cycle lengths, unchanged. The discrepancy growth rate was increased by impaired conduction, fibrosis, conduction heterogeneity, lateral uncoupling, fast gap junction time constant, and steeper action potential duration restitution.

Initiation and entrainment of multicellular automaticity via diffusion limited extracellular domains

Electrically excitable cells often spontaneously and synchronously depolarize in vitro and in vivo preparations. It remains unclear how cells entrain and autorhythmically activate above the intrinsic mean activation frequency of isolated cells with or without pacemaking mechanisms. Recent studies suggest that cyclic ion accumulation and depletion in diffusion-limited extracellular volumes modulate electrophysiology by ephaptic mechanisms (nongap junction or synaptic coupling). This report explores how potassium accumulation and depletion in a restricted extracellular domain induces spontaneous action potentials in two different computational models of excitable cells without gap junctional coupling: Hodgkin-Huxley and Luo-Rudy. Importantly, neither model will spontaneously activate on its own without external stimuli. Simulations demonstrate that cells sharing a diffusion-limited extracellular compartment can become autorhythmic and entrained despite intercellular electrical heterogeneity. Autorhythmic frequency is modulated by the cleft volume and potassium fluxes through the cleft. Additionally, inexcitable cells can suppress or induce autorhythmic activity in an excitable cell via a shared cleft. Diffusion-limited shared clefts can also entrain repolarization. Critically, this model predicts a mechanism by which diffusion-limited shared clefts can initiate, entrain, and modulate multicellular automaticity in the absence of gap junctions.

Multiparameter optimization of nonuniform passive diffusion properties for creating coarse-grained equivalent models of cardiac propagation

The arrhythmogenic role of discrete cardiac propagation may be assessed by comparing discrete (fine-grained) and equivalent continuous (coarse-grained) models. We aim to develop an optimization algorithm for estimating the smooth conductivity field that best reproduces the diffusion properties of a given discrete model. Our algorithm iteratively adjusts local conductivity of the coarse-grained continuous model by simulating passive diffusion from white noise initial conditions during 3-10 ms and computing the root mean square error with respect to the discrete model. The coarse-grained conductivity field was interpolated from up to 300 evenly spaced control points. We derived an approximate formula for the gradient of the cost function that required (in two dimensions) only two additional simulations per iteration regardless of the number of estimated parameters. Conjugate gradient solver facilitated simultaneous optimization of multiple conductivity parameters. The method was tested in rectangular anisotropic tissues with uniform and nonuniform conductivity (slow regions with sinusoidal profile) and random diffuse fibrosis, as well as in a monolayer interconnected cable model of the left atrium with spatially-varying fibrosis density. Comparison of activation maps served as validation. The results showed that after convergence the errors in activation time were < 1 ms for rectangular geometries and 1-3 ms in the atrial model. Our approach based on the comparison of passive properties (<10 ms simulation) avoids performing active propagation simulations (>100 ms) at each iteration while reproducing activation maps, with possible applications to investigating the impact of microstructure on arrhythmias.

Properties of cardiac conduction in a cell-based computational model

The conduction of electrical signals through cardiac tissue is essential for maintaining the function of the heart, and conduction abnormalities are known to potentially lead to life-threatening arrhythmias. The properties of cardiac conduction have therefore been the topic of intense study for decades, but a number of questions related to the mechanisms of conduction still remain unresolved. In this paper, we demonstrate how the so-called EMI model may be used to study some of these open questions. In the EMI model, the extracellular space, the cell membrane, the intracellular space and the cell connections are all represented as separate parts of the computational domain, and the model therefore allows for study of local properties that are hard to represent in the classical homogenized bidomain or monodomain models commonly used to study cardiac conduction. We conclude that a non-uniform sodium channel distribution increases the conduction velocity and decreases the time delays over gap junctions of reduced coupling in the EMI model simulations. We also present a theoretical optimal cell length with respect to conduction velocity and consider the possibility of ephaptic coupling (i.e. cell-to-cell coupling through the extracellular potential) acting as an alternative or supporting mechanism to gap junction coupling. We conclude that for a non-uniform distribution of sodium channels and a sufficiently small intercellular distance, ephaptic coupling can influence the dynamics of the sodium channels and potentially provide cell-to-cell coupling when the gap junction connection is absent.

The electrochemical wave traversing the heart during every beat is essential for cardiac pumping function and supply of blood to the body. Understanding the stability of this wave is crucial to understanding how lethal arrhythmias are generated. Despite this importance, our knowledge of the physical determinants of wave propagation are still evolving. One particular challenge has been the lack of accurate mathematical models of conduction at the cellular level. Because cardiac muscle is an electrical syncytium, in which direct charge transfer between cells drives wave propagation, classical bidomain and monodomain tissue models employ a homogenized approximation of this process. This approximation is not valid at the length scale of single cells, and prevents any analysis of how cellular structures impact cardiac conduction. Instead, so-called microdomain models must be used for these questions. Here we utilize a recently developed modelling framework that is well suited to represent small collections of cells. By applying this framework, we show that concentration of sodium channels at the longitudinal borders of myocytes accelerates cardiac conduction. We also demonstrate that when juxtaposed cells are sufficiently close, this non-uniform distribution induces large ephaptic currents, which contribute to intercellular coupling.

A microstructural model of reentry arising from focal breakthrough at sites of source-load mismatch in a central region of slow conduction

Regions of cardiac tissue that have a combination of focal activity and poor, heterogeneous gap junction coupling are often considered to be arrhythmogenic; however, the relationship between the properties of the cardiac microstructure and patterns of abnormal propagation is not well understood. The objective of this study was to investigate the effect of microstructure on the initiation of reentry from focal stimulation inside a poorly coupled region embedded in more well-coupled tissue. Two-dimensional discrete computer models of ventricular monolayers (1 × 1 cm) were randomly generated to represent heterogeneity in the cardiac microstructure. A small, central poorly coupled patch (0.40 × 0.40 cm) was introduced to represent the site of focal activity. Simulated unipolar electrogram recordings were computed at various points in the tissue. As the gap conductance of the patch decreased, conduction slowed and became increasingly complex, marked by fractionated electrograms with reduced amplitude. Near the limit of conduction block, isolated breakthrough sites occurred at single cells along the patch boundary and were marked by long cell-to-cell delays and negative deflections on electrogram recordings. The strongest determinant of the site of wavefront breakthrough was the connectivity of the brick wall architecture, which enabled current flow through small regions of overlapping cells to drive propagation into the well-coupled zone. In conclusion, breakthroughs at the size scale of a single cell can occur at the boundary of source-load mismatch allowing focal activations from slow conducting regions to produce reentry. These breakthrough regions, identifiable by distinct asymmetric, reduced amplitude electrograms, are sensitive to tissue architecture and may be targets for ablation.

Non-uniform dispersion of the source-sink relationship alters wavefront curvature

The distribution of cellular source-sink relationships plays an important role in cardiac propagation. It can lead to conduction slowing and block as well as wave fractionation. It is of great interest to unravel the mechanisms underlying evolution in wavefront geometry. Our goal is to investigate the role of the source-sink relationship on wavefront geometry using computer simulations. We analyzed the role of variability in the microscopic source-sink relationship in driving changes in wavefront geometry. The electrophysiological activity of a homogeneous isotropic tissue was simulated using the ten Tusscher and Panfilov 2006 action potential model and the source-sink relationship was characterized using an improved version of the Romero et al. safety factor formulation (SF_m2). Our simulations reveal that non-uniform dispersion of the cellular source-sink relationship (dispersion along the wavefront) leads to alterations in curvature. To better understand the role of the source-sink relationship in the process of wave formation, the electrophysiological activity at the initiation of excitation waves in a 1D strand was examined and the source-sink relationship was characterized using the two recently updated safety factor formulations: the SF_m2 and the Boyle-Vigmond (SF_VB) definitions. The electrophysiological activity at the initiation of excitation waves was intimately related to the SF_m2 profiles, while the SF_VB led to several counterintuitive observations. Importantly, with the SF_m2 characterization, a critical source-sink relationship for initiation of excitation waves was identified, which was independent of the size of the electrode of excitation, membrane excitability, or tissue conductivity. In conclusion, our work suggests that non-uniform dispersion of the source-sink relationship alters wavefront curvature and a critical source-sink relationship profile separates wave expansion from collapse. Our study reinforces the idea that the safety factor represents a powerful tool to study the mechanisms of cardiac propagation in silico, providing a better understanding of cardiac arrhythmias and their therapy.

Microscopic variations in interstitial and intracellular structure modulate the distribution of conduction delays and block in cardiac tissue with source-load mismatch

AIMS

Reentrant activity in the heart is often correlated with heterogeneity in both the intracellular structure and the interstitial structure surrounding cells; however, the combined effect of cardiac microstructure and interstitial resistivity in regions of source-load mismatch is largely unknown. The aim of this study was to investigate how microstructural variations in cell arrangement and increased interstitial resistivity influence the spatial distribution of conduction delays and block in poorly coupled regions of tissue.

METHODS AND RESULTS

Two-dimensional 0.6 cm × 0.6 cm computer models with idealized and realistic cellular structure were used to represent a monolayer of ventricular myocytes. Gap junction connections were distributed around the periphery of each cell at 10 μm intervals. Regions of source-load mismatch were added to the models by increasing the gap junction and interstitial resistivity in one-half of the tissue. Heterogeneity in cell shape and cell arrangement along the boundary between well-coupled and poorly coupled tissue increased variability in longitudinal conduction delays to as much as 10 ms before the onset of conduction block, resulting in wavefront breakthroughs with pronounced curvature at distinct points along the boundary. Increasing the effective interstitial resistivity reduced source-load mismatch at the transition boundary, which caused a decrease in longitudinal conduction delay and an increase in the number of wavefront breakthroughs.

CONCLUSION

Microstructural variations in cardiac tissue facilitate the formation of isolated sites of wavefront breakthrough that may enable abnormal electrical activity in small regions of diseased tissue to develop into more widespread reentrant activity.