Cellular Size, Gap Junctions, and Sodium Channel Properties Govern Developmental Changes in Cardiac Conduction

Electrical conduction in cardiac ventricular tissue is regulated via sodium (Na⁺) channels and gap junctions (GJs). We and others have recently shown that Na⁺channels preferentially localize at the site of cell-cell junctions, the intercalated disc (ID), in adult cardiac tissue, facilitating coupling via the formation of intercellular Na⁺nanodomains, also termed ephaptic coupling (EpC). Several properties governing EpC vary with age, including Na⁺channel and GJ expression and distribution and cell size. Prior work has shown that neonatal cardiomyocytes have immature IDs with Na⁺channels and GJs diffusively distributed throughout the sarcolemma, while adult cells have mature IDs with preferentially localized Na⁺channels and GJs. In this study, we perform an in silico investigation of key age-dependent properties to determine developmental regulation of cardiac conduction. Simulations predict that conduction velocity (CV) biphasically depends on cell size, depending on the strength of GJ coupling. Total cell Na⁺channel conductance is predictive of CV in cardiac tissue with high GJ coupling, but not correlated with CV for low GJ coupling. We find that ephaptic effects are greatest for larger cells with low GJ coupling typically associated with intermediate developmental stages. Finally, simulations illustrate how variability in cellular properties during different developmental stages can result in a range of possible CV values, with a narrow range for both neonatal and adult myocardium but a much wider range for an intermediate developmental stage. Thus, we find that developmental changes predict associated changes in cardiac conduction.

1. Introduction

It is well-established that conduction in cardiac tissue is regulated by ionic currents and gap junction (GJ) coupling (Shaw and Rudy, 1997; Kucera et al., 2002). In ventricular tissue, the voltage-gated sodium (Na⁺) channel, Na_v1.5, is primarily responsible for generating the depolarizing Na⁺current (I_Na), and connexin43 (Cx43) is the primary GJ protein, facilitating the passive current flow between adjacent cells (Veeraraghavan et al., 2014), both of which mediate conduction. Altering either I_Naor GJ coupling can lead to changes in conduction and ultimately increase the risk of arrhythmias (Quan and Rudy, 1990; Shaw and Rudy, 1997; Rohr et al., 1998).

We (Veeraraghavan et al., 2015; Veeraraghavan and Gourdie, 2016; Mezache et al., 2020) and others (Kucera et al., 2002; Rhett et al., 2012; Agullo-Pascual et al., 2014; Leo-Macias et al., 2016) have shown that Na_v1.5channels preferentially localize at the intercalated disc (ID), the area of cell-cell junctions in cardiac tissue. Multiple in silico studies have hypothesized that I_Naat the ID can be altered via Na⁺nanodomain signaling at the intercellular cleft space (Kucera et al., 2002; Sperelakis, 2002; Mori et al., 2008; Lin and Keener, 2010; Wei et al., 2016; Tveito et al., 2017; Weinberg, 2017; Hichri et al., 2018; Jæger et al., 2019; Wei and Tolkacheva, 2020). This, in turn, can modulate cell-cell coupling through a mechanism known as ephaptic coupling (EpC). In this paper we consider two primary effects of EpC: electrical field effects and Na⁺depletion in the intercellular cleft (i.e., the narrow extracellular space between electrically coupled cells at the ID). To elaborate briefly, EpC is governed by the following: Na⁺influx during the action potential upstroke in an “upstream” or pre-junctional depolarizing cell during a propagating electrical wave decreases the electrical potential within the intercellular cleft. This reduction of the potential within the intercellular cleft then depolarizes the “downstream” or post-junctional apposing cell from the extracellular, rather than the intracellular, side of the cell membrane. Additionally, Na⁺influx reduces the Na⁺concentration within the intercellular cleft, which governs the flux of the Na⁺channels at the ID in both cells. The width or volume of the intercellular cleft space is one of the key properties governing the magnitude of these effects. When the intercellular cleft is narrow, both the elevated transmembrane potential (V_m) and locally depleted Na⁺concentration within the intercellular cleft reduce the Na⁺current driving force and, therefore, the Na⁺current. This reduction in I_Nahas been termed “self-attenuation” and has been shown to slow conduction velocity (CV) (Kucera et al., 2002; Sperelakis, 2002; George et al., 2016; Hichri et al., 2018).

Several key properties governing EpC and conduction overall are known to change during development: Cells in neonatal myocardium do not have fully formed IDs, and Na⁺channels and GJs are distributed diffusively throughout the sarcolemma (Fromaget et al., 1992; Vreeker et al., 2014). Consistent with reduced Na⁺channel expression, it has been shown that pediatric cardiomyocytes produce a reduced I_Na, compared to adult cardiomyocytes (Cai et al., 2011; Cordeiro et al., 2013). Cx43 is essentially undetectable until 23 weeks in utero and remain randomly distributed on the sarcolemma in neonatal cardiomyocytes (Peters et al., 1994; Hirschy et al., 2006; Vreeker et al., 2014; Swift et al., 2020). Vreeker et al. showed that Cx43 tends to relocate to the lateral membrane around 5 months postnatal and does not begin to preferentially localize at the ID until around 2.5–5 years old, with full preferential localization occurring at roughly 7 years of age (Vreeker et al., 2014). Na_v1.5channels, however, reside on the lateral membrane in neonatal cardiomyocytes and begin to begin to preferentially localize at the ID around 5 months postnatal, much earlier than Cx43 (Harrell et al., 2007; Vreeker et al., 2014). Additionally, studies have shown that adult cardiomyocytes are larger than neonatal cardiomyocytes (Cordeiro et al., 2013; Vreeker et al., 2014; Swift et al., 2020). This is especially important given that cell size broadly influences all electrical activity in the cell by altering surface area, cell volume, membrane capacitance, ion channel expression, etc. (Kato et al., 1996; Spach et al., 2000).

We hypothesize that the developmental-associated increase in both Na⁺and GJ current will increase CV from the neonatal stage to the adult and that this regulation will be influenced by the relative strength of EpC at different developmental stages. In the paper, we perform a wide parameter investigation, varying age-associated parameters including gap junctional conductance (f_gap), cell size (S), Na⁺channel density (ρ_Na), and Na⁺channel ID localization (ID_Na) and measure CV in simulated cardiac tissue. To our knowledge, no studies have investigated the interdependence of these parameters on conduction within a health myocardium during development from neonatal to adult tissue. While one study investigated the changes on impulse conduction in the canine myocardium, it was limited to a comparison of 8 week old postnatal to adult purkinje fibers (Rosen et al., 1981). Thus, in our study, we investigate conduction through ventricular tissue in a range of developmental stages and conditions. We find that CV has a biphasic dependence on cell size, in a manner that depends on GJ coupling. In addition, we find that ephaptic effects play a larger role in conduction for larger cells with low GJ coupling. Interestingly, simulations predict that variability in cellular properties in intermediate developmental stage between neonatal cardiomyocytes and adult cardiomyocytes can lead to a wide range of possible CV values, but this range narrows to a normal range with adult tissue-associated parameters.

2. Methods

Full details of the computational model are provided in Supplementary Material. Briefly, we simulate a 50-cell cable of guinea pig ventricular myocytes (Livshitz and Rudy, 2007) that incorporates a Markov chain model for the wild-type (WT) Na⁺channel (Clancy et al., 2002), shown in Figure 1. We note that a Markov chain formulation for I_Na was not necessary to reproduce our results, as prior studies have shown EpC effects are reproduced by Hodgkin-Huxley Na⁺channel models (Kucera et al., 2002; Weinberg, 2017; Moise et al., 2021), but rather is used to facilitate an appropriate comparison with our prior work simulating Na_v1.5gain-of-function mutations (Nowak et al., 2021), which were also modeled with a Markov chain. As in our previous work (Greer-Short et al., 2017; Weinberg, 2017; Nowak et al., 2020, 2021) and work performed by others (Kucera et al., 2002), we account for non-uniform Na⁺subcellular localization by spatially discretizing each cell into two ID membrane patches at the ends of the cell and axial membrane patches along the length of the cell. The number of axial patches varied with the size of the cell, as described below, with each axial patch fixed in length (L_p = 10 μm).

FIGURE 1

Figure 1. Schematic of the computational model. (A) Electric circuit representation of coupled myocytes. Intracellular nodes are coupled via a myoplasmic resistance (R_myo). End nodes are coupled via a gap junctional resistance (R_gap). Extracellular potentials at the disc and intercellular cleft (${\varphi }_e^{disc}$ and ${\varphi }_e^{cleft}$, respectively) are governed by a T-shaped network of two axial resistances in the intercellular cleft (R_cl) and one radial resistance (R_radial). (B) Na⁺concentration in diffusively coupled compartments, including intracellular Na⁺in the axial and disc compartments (${\left[N{a}^{+}\right]}_i^{ax}$ and ${\left[N{a}^{+}\right]}_i^{disc}$) and extracellular Na⁺in the intercellular cleft and bulk spaces (${\left[N{a}^{+}\right]}_e^{cleft}$ and ${\left[N{a}^{+}\right]}_e^{bulk}$). (C) Representation of age-associated change in model parameters, including changes in cell size (S), Na⁺ channel density (ρ_Na), and Na⁺channel ID localization (ID_Na).

Cells are coupled via GJs and EpC: gap junctional coupling is represented via gap junctional conductances coupling ID nodes of adjacent cells (g_gap, represented by resistor R_gap in the electrical circuit). EpC is represented by a T-shaped junction of two intercellular cleft resistances (R_cl) and a radial bulk (R_radial), which are proportional and inversely proportional to intercellular cleft width w, respectively. Intracellular nodes are coupled with a myoplasmic resistance (R_myo) (Figure 1A). The nominal value for g_gap is defined as $g_{gap}^0=1266nS$, and changes in GJ coupling are accounted for by adjusting the GJ scaling factor f_gap (between 0 and 1), such that $g_{gap}=f_{gap}{g}_{gap}^0$.

A wide range of macroscopic GJ conductance values have been measured experimentally, with values ranging from the low 10s of nanosiemens up to approximately 2,000 nS (Weingart, 1986; Wittenberg et al., 1986; White et al., 1990; Moreno et al., 1994; Kwak and Jongsma, 1996; Verheule et al., 1997; Kucera et al., 2002; Valiunas et al., 2002; Desplantez et al., 2007; McCain et al., 2012; Nielsen et al., 2012). While studies have shown increases in Cx43 expression with age, to our knowledge, no studies have quantified GJ conductances throughout development. Therefore, we performed simulations with GJ conductance spanning from the low to high end of physiological measurements and hypothesized that such a range is qualitatively similar to different developmental stages and consistent with Cx43 expression changes. Importantly, the values chosen were consistent with experimental measurements of conduction for neonatal and adult myocardium (Rosen et al., 1981; George et al., 2019; King et al., 2021). We note that age-associated changes in GJ localization are represented as changes in GJ coupling (i.e., changes in f_gap), as all GJs are located at cell ends in the one-dimensional tissue model.

We account for dynamic [Na⁺] in three spaces: (i) the ID (${\left[\text{N}{\text{a}}^{+}\right]}_i^{disc}$), mediated by $I_{Na}^{disc}$ and intracellular diffusion, (ii) the intercellular cleft (${\left[\text{N}{\text{a}}^{+}\right]}_e^{cleft}$), with volume proportional to w and mediated by ID Na⁺current $I_{Na}^{disc}$ and passive diffusion with the bulk extracellular space, and (iii) the axial intracellular space (${\left[\text{N}{\text{a}}^{+}\right]}_i^{ax}$), mediated by axial Na⁺current $I_{Na}^{ax}$ and intracellular diffusion. The cable was paced at one end with a specified basic cycle length (BCL). Unless otherwise stated, for all simulations, we utilize a BCL = 500 ms or a pacing rate of 2 Hz, which is normal pacing for the guinea pig model that is utilized for this study.

In addition to GJ coupling, we perform simulations in which we adjust several key age-associated properties: (1) We vary Na⁺channel localization at the ID (ID_Na) between 0.1 and 1, to account for Na⁺channel redistribution, where ID_Na = 1 represents 100% or all Na⁺ channels localized at the ID. (2) We define and vary a cell size scaling factor (S) between 0.2 and 1 to account for cell size growth, where 100% cell size represents an adult cell. The cell geometry is assumed to be cylindrical, with radius r and length L, defined as r = Sr₀ and L = n_pL_p, where the nominal adult radius r₀ = 11 μm, the number of axial patches $n_p=S{n}_p^0$, and the maximum axial membrane patches $n_p^0=10$. Note, we only consider values of Ssuch that n_p is a whole number, and that nominal adult cell length $L_0=n_p^0{L}_p=100\mu \text{m}$. Additionally, note that since S scales both cell length and radius, cell membrane surface area is scaled by S², e.g., S = 0.4 corresponds with surface area scaled by a factor of 0.16. (3) We also vary the Na⁺channel density (ρ_Na) between 0.2 and 1 to account for age-associated changes in Na⁺channel expression, where 100% ρ_Narepresents full expression of Na⁺channels on the cell membrane. The total cellular Na⁺conductance is proportional to both ρ_Naand total cell surface area, such that we can define a normalized total Na⁺channel conductance ($G_{Na}={\rho }_{Na}{S}^2$, also between 0 and 1), where the total Na⁺channel conductance (in physical units) is equal to G_Na, scaled by the nominal total Na⁺channel conductance ($G_{Na}^0=21.78\text{ mS}/\text{c}{\text{m}}^2$). Thus, for ρ_Naand Sof 100%, total Na⁺ channel conductance is $G_{Na}^0$. (4) Finally, we vary the intercellular cleft width (w) from 10 to 40 nm, consistent with intercellular cleft width ranges measured at the ID in our previous work (Veeraraghavan and Gourdie, 2016; Greer-Short et al., 2017; Nowak et al., 2020, 2021).

3. Results

3.1. Ephaptic Effects Are Enhanced for Larger Cells With Low Gap Junctional Coupling

Motivated by findings that both GJ coupling (Peters et al., 1994; Hirschy et al., 2006; Vreeker et al., 2014; Swift et al., 2020) and cell size (Kato et al., 1996; Spach et al., 2000; Cordeiro et al., 2013; Vreeker et al., 2014) increase with age, we first investigate how changes in cell size and GJ coupling influence conduction. We first consider moderate localization of Na⁺channels at the ID (ID_Na = 50%), high whole cell Na⁺channel density (ρ_Na = 100%), and a nominal intercellular cleft width (w = 20 nm). In Figure 2, the time series for transmembrane voltage (V_m), pre- and post-junctional Na⁺current ($I_{Na}^{pre}$and $I_{Na}^{post}$, respectively), GJ current (I_GJ), and the cleft voltage (V_cleft) are shown during the action potential upstroke at the same spatial location within the cardiac tissue for varying cell size and GJ coupling.

FIGURE 2

Figure 2. Conduction velocity (CV) depends on cell size and GJ coupling. (A) Transmembrane voltage of the post-junctional node of the ID (V_m), (B) Na⁺current at the pre-junctional node of the ID ($I_{Na}^{pre}$), (C) Na⁺current at the post-junctional node of the ID ($I_{Na}^{post}$), (D) GJ current at the ID (I_GJ), and (E) cleft voltage (V_cleft) are shown in tissue for low (50.6 nS), moderately low (101 nS), moderately high (253 nS), and high (1,266 nS) GJ coupling. For clarity, traces are shown for the same spatial point of 1.2 mm from the pacing site. Parameters: ID_Na= 50%, ρ_Na= 100%, w = 20 nm.

For low GJ coupling (Figure 2, left), the slowed action potential upstroke for the large cell size indicates that conduction slows as cell size increases. For tissue with small cell size, a large pre-junctional current $I_{Na}^{pre}$rapidly activates the post-junctional current $I_{Na}^{post}$ (Figures 2B,C, red), and the smaller membrane surface area results in a relatively larger GJ current density I_GJ (Figure 2D, red), compared with larger cell size. In contrast, for tissue with larger cell size, GJ current density is reduced (Figure 2D, purple), and further cleft hyperpolarization is enhanced, i.e., V_cleftis more negatively polarized (Figure 2E, purple). This enhanced EpC effect ultimately drives the self-attenuation mechanism, in which $I_{Na}^{post}$driving force is reduced. Collectively, these effects result in slower conduction for tissue with larger cell size.

For moderately low and moderately high GJ coupling (Figure 2, center left and center right), these EpC effects are reduced; however, cleft hyperpolarization is still enhanced for larger cell sizes (Figures 2B–E). For high GJ coupling (Figure 2, right), the faster action potential upstroke for large cell size indicates that conduction is enhanced as cell size increases. In this case, I_GJis sufficiently large, such that the self-attenuation effects on post-junctional Na⁺current are counterbalanced, in a manner that conduction is faster for tissue with larger cell size. Additionally, for larger cell size, there are fewer cell-cell junctions in a given length of tissue, which also results in overall faster conduction.

3.2. Conduction Velocity Depends on Key Cellular and Tissue Properties

We next more broadly investigated the interdependence between cell size, whole cell Na⁺channel density (ρ_Na), Na⁺channel ID localization (ID_Na), and GJ coupling on CV (Figure 3). We consider a wide range of parameters to assess the age-associated changes on CV. For all cases investigated, as expected, CV consistently increases with both increasing GJ coupling and ρ_Na. For low GJ coupling (Figure 3A), as illustrated above, CV decreases as cell size increases, due to both the reduced GJ current density and the self-attenuation mechanism previously stated. Increased Na⁺channel preferential localization at the ID tends to moderately increase CV for larger cell sizes due to increased EpC, which increases conduction for low GJ coupling.

FIGURE 3

Figure 3. Conduction velocity (CV) depends on key cellular and tissue properties. CV is shown as a function of cell size for different values of Na⁺channel densities (ρ_Na) for low (A), moderately low (B), moderately high (C), and high (D) GJ coupling and 10% (left), 50% (middle), and 90% (right) Na⁺channel ID localization (ID_Na). Parameters: Cleft width w = 20 nm. Parameter regimes associated with neonatal (gray boxes) and adult (black boxes) tissue are highlighted.

In contrast, for high GJ coupling (Figure 3D), CV increases as cell size increases due to fewer GJs per unit length, as described above. As ID_Naincreases, CV decreases moderately, more so for larger cell sizes, again due to increased EpC, which decreases conduction for high GJ coupling. From low to high GJ coupling, the CV-cell size relationship transitions from decreasing to increasing, such that for moderate GJ coupling strengths, CV exhibits a biphasic relationship with cell size (Figures 3B,C). However, for most conditions, CV only moderately varies across the wide range of cell sizes, which suggests GJ and EpC effects that depend on cell size are fairly balanced for these conditions. Increased Na⁺channel ID localization has only small effects, tending to increase CV for the moderately low GJ coupling cases.

We indicate two key regions of interest: the neonatal regime and the adult regime. The neonatal regime is indicated by gray boxes and represents the range of parameters associated with neonatal myocardium: small cell size, low Na⁺channel ID localization, and low to moderately low GJ coupling (Peters et al., 1994; Kato et al., 1996; Spach et al., 2000; Vreeker et al., 2014). The adult regime is indicated by black boxes and represents the range of parameters associated with adult myocardium: large cell size, high preferential Na⁺channel localization at the ID, and moderately high to high GJ coupling. In the neonatal regimes, CV is slow, consistent with the reduced GJ coupling, and further is sensitive to changes in cell size. In contrast, in the adult regime, CV is faster, consistent with the higher GJ coupling, and additionally does not depend on cell size. In both regimes, CV is highly sensitive to changes in Na⁺channel density. Collectively, this suggests that developmental changes result in conduction that is more robust to changes in cell size. Additionally, in Supplementary Material, we perform a similar broad investigation for a bradycardic pacing rate (1 Hz, or BCL of 1,000 ms), and we found that CV has a nearly identical relationship with cell size, Na⁺channel density and localization, and GJ coupling under these conditions (Supplementary Figure 1), as compared with the normal pacing rate.

3.3. Conduction Velocity Is Correlated to Total Cell Na⁺Conductance and Gap Junctional Coupling

The previous analysis shows CV had a consistent positive dependence on Na⁺current density ρ_Na, yet the relationship between CV and cell size depends on GJ coupling. As the total cell Na⁺conductance (G_Na) depends on both cell size and Na⁺channel density, we next investigate how CV depends on G_Nafor varying GJ coupling conditions (Figure 4). For all GJ coupling strengths, we find discrete vertical “columns” of points, which correspond to the same G_Navalue (i.e., same ρ_Naand S) and different values of Na⁺channel localization (ID_Na). From these data, we observe that increased ID_Natypically slows conduction, except for cases of moderately low GJ coupling and larger cell size. Additionally, we observe diagonal “bands” of points, which correspond with the same cell sizes and different ρ_Navalues, which illustrate that CV increases with increasing ρ_Na, as in Figure 2.

FIGURE 4

Figure 4. Correlation between conduction velocity (CV) and total cell Na⁺ conductance (G_Na) increases as gap junctional coupling increases. Conduction velocity (CV) is shown as a function of G_Na for different cell sizes and Na⁺ channel ID localization for (A) low, (B) moderately low, (C) moderately high, and (D) high GJ coupling. Pearson correlation coefficients r for low (r = −0.103), moderately low (r = 0.283), moderately high (r = 0.607), high (r = 0.865) GJ coupling. Parameters: Cleft width w = 20 nm.

With all of these dependencies within given “columns” or “bands,” we also consider if there is a distinct relationship between CV and G_Nafor a given GJ coupling strength. For low GJ coupling (Figure 4A), CV is overall very weakly correlated with G_Na (dashed black line), indicating that for these conditions, the total Na⁺channel conductance is not predictive of CV. However, as GJ coupling increases, the correlation between CV and G_Naincreases and approaches 1 for high GJ coupling (Figures 4B–D), demonstrating a closer relationship between conduction and total Na⁺current conductance.

Prior work from us and others has shown that the intercellular width (w) is a critical parameter governing the strength of EpC effects and ultimately conduction (Kucera et al., 2002; Lin and Keener, 2010; Greer-Short et al., 2017; Nowak et al., 2020), so we next investigate how the correlation between CV and G_Nadepends on w. In Figure 5, we plot the Pearson correlation coefficient r as a function of GJ coupling for different values of w. For all cleft widths, the correlation increases from near 0 to near 1 as GJ coupling increases, as in Figure 4. However, for narrow clefts (red), the correlation is less sensitive to changes in GJ coupling, such that the correlation coefficient is more positive for lower GJ coupling and less positive for higher GJ coupling, compared with wider clefts. These results are consistent with EpC playing a larger role governing conduction for lower GJ coupling, such that there is a stronger relationship between CV and G_Nafor these conditions.

FIGURE 5

Figure 5. Correlation between conduction velocity (CV) and total cell Na⁺ conductance depends on gap junctional coupling. The Pearson correlation coefficient r between CV and total cell Na⁺ conductance (G_Na) is shown as a function of GJ coupling conductance, for four different cleft width w values.

3.4. Conduction Velocity Depends on Developmental Stage

The above results illustrate the complex relationship between conduction and cellular/tissue properties known to alter with age and development. As previously discussed, neonatal tissues are associated with low Na⁺channel expression and ID localization and smaller cell size, while adult tissues have larger cells and high Na⁺channel expression and localization at the ID. We hypothesize that, during both disease and development, these properties remain variable in both time and between patients. For a final analysis, we investigate CV for conditions representing an age-associated progression, considering neonatal, intermediate developmental stages, and adult tissue (Figure 6). To account for variability, we consider a range of values for each of the key parameters investigated throughout this study: cell size, Na⁺channel density, and Na⁺channel localization at the ID. These parameters were varied over a 20% range. The minimum and maximum CV values over all parameter conditions for each stage were plotted (solid black lines), along with the average over the different parameter values (dashed black line). Due to lack of evidence of the precise order in which these parameters change throughout the developmental process, we consider different possible parameter conditions for the intermediate stages: uniform ρ_Naand ID_Naincrease, staged Sand ρ_Naincreases, staged Sand ID_Naincreases, and staged ρ_Naand ID_Naincreases ranging from neonatal to adult tissue parameter values. The parameter ranges for each developmental stage are listed in Supplementary Table 2.

FIGURE 6

Figure 6. Age-associated ranges of conduction velocity (CV). The range of CV values (see text for details) are shown as functions of age-dependent progression. Parameter ranges and markers for each stage are shown in Supplementary Table 2. The dashed line shown for each stage represents the average CV value over all conditions within the specified parameter ranges. Cleft width w = 20 nm.

For the neonatal stage, the range of CVs is small and conduction speeds are slow. For the first intermediate developmental stage, the CV range is moderately larger, but the conduction remains slow. However, for the second intermediate developmental stage, the CV range is dramatically increased, ranging from roughly 20 cm/s to over 50 cm/s, suggesting that conduction greatly increases and can approach adult myocardium values, but may be more variable during this intermediate stage. For the adult stage, the range of CVs decreases, such that conduction is consistently within the normal faster propagation regime, suggesting robust conduction regardless of small changes in cellular and tissue properties.

4. Discussion

In this study, we investigate the regulation of key age-dependent properties, specifically cell size, Na⁺channel density and localization, and GJ coupling, on cardiac conduction. We summarize our key findings: CV is consistently increased by increased Na⁺current density, across all conditions. However, simulations predict that CV biphasically depends on cell size, depending on the strength of GJ coupling. That is, CV increases with increasing cell size for high GJ coupling, yet decreases with increasing cell size for low GJ coupling. As a consequence, CV and total cell Na⁺channel conductance are well-correlated in cardiac tissue with high GJ coupling, but not correlated with CV for low GJ coupling. We predict that the role of EpC governing conduction changes during development, such that neonatal tissue is less sensitive to changes in EpC due to smaller Na⁺channel ID localization. We postulate that even though GJ coupling is low during this early stage, the small cell size and therefore small membrane capacitance is such that the lower Na⁺current density is still sufficient to maintain robust conduction in the myocardium. However, conduction is very slow during this stage due to the relatively high ratio of cell-cell junctions for a given length of tissue. These findings are consistent with Swift et al. (2020), who found slower atrioventricular conduction in neonatal and early postnatal rats, compared with adult rat myocardium.

As development progresses, model predictions are consistent with faster conduction. Further, simulations predict that a wide range of conduction velocities are possible during intermediate developmental stages due to variability in cellular/tissue properties and the relative timing of developmental changes. However, despite similar variability in parameter values, this variability in CV prediction narrows in adult tissue and is consistent with experimental measures (George et al., 2019; King et al., 2021). Interestingly, we see a larger influence of EpC effects in larger cell sizes with reduced GJ coupling, as would be the case in intermediate developmental stages and in adult myocardium for pathological conditions such heart failure (Smith et al., 1991; Peters et al., 1997; Yao et al., 2003; Akar et al., 2004; Poelzing and Rosenbaum, 2004), indicating a possible mechanism for maintained conduction during such transitional or diseased states. Additionally, this enhanced variability in conduction for intermediate developmental stages may be desirable, as this variability suggests an ability to adapt and modulate cardiac activity in response to developmental perturbations, which inherently vary significantly for different individuals.

Previous studies have investigated the properties of Na⁺channels and GJs and their roles in the developing heart. Work by Harrell et al. (2007) found that Na_v1.5mRNA was significantly more up-regulated in adult mouse hearts than in neonatal hearts. Similarly, Cordeiro et al. (2013) find that both peak and late I_Nais significantly smaller in neonatal canine cardiomyocytes compared to that of adults, and Cai et al. (2011) found the same in human atrial cardiomyocytes. Vreeker et al. (2014) show that the GJ protein Cx43 transitions from diffusely distributed on the cardiomyocyte membrane to highly associated with mechanical junctions at the ID during development in human samples, as similarly shown in rat myocardium (Angst et al., 1997). In contrast, Peters et al. (1994) found that GJs and adherens junctions have a highly correlated distribution over postnatal development of the ventricle in human samples. Incorporating developmental changes in GJ distribution, Spach et al. (2000) previously simulated a neonatal cardiomyocyte with diffuse GJs and smaller size and found that conduction was slower in the neonatal tissue, compared with adult myocardium with GJs primarily localized at the ID and larger size.

While these studies are crucial in understanding cardiac development, to our knowledge, our study is the first to combine all these key developmental changes occurring in conjunction with each other, specifically cellular size, gap junctional conductance, and Na⁺channel expression and distribution. Thus, we find that ephaptic effects are more pronounced in larger cells with low GJ coupling. We also find that CV biphasically depends on cell size in a manner dependent on GJ coupling: CV is relatively fast for both low GJ coupling and small cell size and high GJ coupling and large cell size. Interestingly, we find that CV is correlated total Na⁺conductance for high GJ coupling. Finally, by incorporating previous data on developmental changes (Peters et al., 1994; Spach et al., 2000; Harrell et al., 2007; Cai et al., 2011; Cordeiro et al., 2013; Vreeker et al., 2014), we predict how the variability in the ranges for conduction change in conjunction with developmental stages.

Recently, we investigated the age-dependent manifestation of a long QT type 3 (LQT3)-associated gain-of-function mutation in Na_v1.5 (Nowak et al., 2021). LQT3, while relatively rare, has a high mortality, reaching 49% (Vignati, 2007), and critically patients tend to remain asymptomatic until well after puberty (Beaufort-Krol et al., 2005; Wilde et al., 2016; Kutyifa et al., 2018). Interestingly, in pediatric patients with LQT syndromes, Na⁺channel blockers are commonly prescribed as a form of chronic management (Hanisch, 2001), and our results are consistent with the safety of this approach due to the weak dependence of conduction on overall total cell Na⁺channel conductance. In our study, we predicted that not only does the manifestation of arrhythmias depend on developmental stage, but that the sensitivity to changes in intercellular cleft width depended on these age-associated properties as well. Interestingly, Brugada syndrome (BrS), an inherited cardiac arrhythmia disorder caused a loss-of-function mutation in Na_v1.5, also often manifests later in life; the average first event occurs around 42 years old (Milman et al., 2017). Our study is consistent with this clinical manifestation, as simulations predict that in these early developmental stages with low GJ coupling, conduction is less sensitive to overall Na⁺conductance (Figure 4), due to the overall small cell sizes in neonatal tissue.

Finally, we acknowledge limitations of our study. To study developmental changes in Na⁺channel distribution, we incorporate this critical subcellular detail in the representation of cardiac tissue; however, our model is still a simplification of the complex cardiac tissue structure. Specifically, our model assumes a simplified cylindrical cell and uniform intercellular cleft, while the geometry of individual cells and the ID structure is known to be complex heterogeneous (Veeraraghavan et al., 2015). Future work will focus on investigating how ID structure changes in development and impacts conduction (Moise et al., 2021). Additionally, cardiac tissue is a three-dimensional structure, and our one-dimensional cable representation cannot fully represent all aspects of developmental changes, such as GJ localization along the lateral membrane. In a two- or three-dimensional tissue, GJ redistribution from the lateral membrane to the ID would also be expected to impact conduction heterogeneously dependent on the direction of wavefront propagation relative to the underlying tissue geometry, such that early developmental stages are associated with isotropic conduction while anisotropic conduction is associated with adult myocardium, consistent with previous work from Spach et al. (2000). Further, while we observe minimal differences in CV between normal and slow pacing rates, Entz et al. (2016) previously showed that CV slowing at faster pacing rates differed between longitudinal and transverse conduction, and further was regulated by extracellular ionic composition, suggesting that properties governing EpC may similarly modulate heart rate dependence. Additionally, expression of ion channels vary across individual patients, and these differences may impact conduction during development. Further, recent work from the Posnack lab has shown that several expression of several key ion channel and calcium handling proteins vary with developmental stage in rats (Swift et al., 2020). We are particularly interested in similar changes in human myocardium, and future work will incorporate such details as available and established in the literature.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author Contributions

MN performed computational studies and wrote the original manuscript draft. RV contributed to the manuscript draft. SP and SW designed simulations, provided funding, and contributed to the manuscript draft. All authors contributed to the article and approved the submitted version.

Funding

This study was supported by funding from the National Institutes of Health, grant numbers R01HL138003 (SW, SP), R01HL102298 (SP), and R01HL148736 (RV).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2021.731025/full#supplementary-material

References

Agullo-Pascual, E., Lin, X., Leo-Macias, A., Zhang, M., Liang, F.-X., Li, Z., et al. (2014). Super-resolution imaging reveals that loss of the c-terminus of connexin43 limits microtubule plus-end capture and nav1. 5 localization at the intercalated disc. Cardiovasc. Res. 104, 371–381. doi: 10.1093/cvr/cvu195

PubMed Abstract | CrossRef Full Text | Google Scholar

Akar, F. G., Spragg, D. D., Tunin, R. S., Kass, D. A., and Tomaselli, G. F. (2004). Mechanisms underlying conduction slowing and arrhythmogenesis in nonischemic dilated cardiomyopathy. Circ. Res. 95, 717–725. doi: 10.1161/01.RES.0000144125.61927.1c

PubMed Abstract | CrossRef Full Text | Google Scholar

Angst, B. D., Khan, L. U., Severs, N. J., Whitely, K., Rothery, S., Thompson, R. P., et al. (1997). Dissociated spatial patterning of gap junctions and cell adhesion junctions during postnatal differentiation of ventricular myocardium. Circ. Res. 80, 88–94. doi: 10.1161/01.RES.80.1.88

PubMed Abstract | CrossRef Full Text | Google Scholar

Beaufort-Krol, G. C., van den Berg, M. P., Wilde, A. A., van Tintelen, J. P., Viersma, J. W., Bezzina, C. R., et al. (2005). Developmental aspects of long QT syndrome type 3 and Brugada syndrome on the basis of a single SCN5A mutation in childhood. J. Am. Coll. Cardiol. 46, 331–337. doi: 10.1016/j.jacc.2005.03.066

PubMed Abstract | CrossRef Full Text | Google Scholar

Cai, B., Mu, X., Gong, D., Jiang, S., Li, J., Meng, Q., et al. (2011). Difference of sodium currents between pediatric and adult human atrial myocytes: evidence for developmental changes of sodium channels. Int. J. Biol. Sci. 7:708. doi: 10.7150/ijbs.7.708

PubMed Abstract | CrossRef Full Text | Google Scholar

Clancy, C. E., Tateyama, M., and Kass, R. S. (2002). Insights into the molecular mechanisms of bradycardia-triggered arrhythmias in long QT-3 syndrome. J. Clin. Investig. 110, 1251–1262. doi: 10.1172/JCI0215928

PubMed Abstract | CrossRef Full Text | Google Scholar

Cordeiro, J. M., Panama, B. K., Goodrow, R., Zygmunt, A. C., White, C., Treat, J. A., et al. (2013). Developmental changes in expression and biophysics of ion channels in the canine ventricle. J. Mol. Cell. Cardiol. 64, 79–89. doi: 10.1016/j.yjmcc.2013.09.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Desplantez, T., Dupont, E., Severs, N. J., and Weingart, R. (2007). Gap junction channels and cardiac impulse propagation. Journal of Membrane Biology 218, 13–28. doi: 10.1007/s00232-007-9046-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Entz, M. I., George, S. A., Zeitz, M. J., Raisch, T., Smyth, J. W., and Poelzing, S. (2016). Heart rate and extracellular sodium and potassium modulation of gap junction mediated conduction in guinea pigs. Front. Physiol. 7:16. doi: 10.3389/fphys.2016.00016

PubMed Abstract | CrossRef Full Text | Google Scholar

Fromaget, C., El Aoumari, A., and Gros, D. (1992). Distribution pattern of connexin 43, a gap junctional protein, during the differentiation of mouse heart myocytes. Differentiation 51, 9–20. doi: 10.1111/j.1432-0436.1992.tb00675.x

PubMed Abstract | CrossRef Full Text | Google Scholar

George, S. A., Bonakdar, M., Zeitz, M., Davalos, R. V., Smyth, J. W., and Poelzing, S. (2016). Extracellular sodium dependence of the conduction velocity-calcium relationship: evidence of ephaptic self-attenuation. Am. J. Physiol. Heart Circ. Physiol. 310, H1129–H1139. doi: 10.1152/ajpheart.00857.2015

PubMed Abstract | CrossRef Full Text | Google Scholar

George, S. A., Hoeker, G., Calhoun, P. J., Entz, M., Raisch, T. B., King, D. R., et al. (2019). Modulating cardiac conduction during metabolic ischemia with perfusate sodium and calcium in guinea pig hearts. Am. J. Physiol. Heart Circ. Physiol. 316, H849–H861. doi: 10.1152/ajpheart.00083.2018

PubMed Abstract | CrossRef Full Text | Google Scholar

Greer-Short, A., George, S. A., Poelzing, S., and Weinberg, S. H. (2017). Revealing the concealed nature of long-QT type 3 syndrome. Circ. Arrhythmia Electrophysiol. 10:e004400. doi: 10.1161/CIRCEP.116.004400

PubMed Abstract | CrossRef Full Text | Google Scholar

Hanisch, D. (2001). Pediatric arrhythmias. J. Pediatr. Nurs. 16, 351–362. doi: 10.1053/jpdn.2001.26571

PubMed Abstract | CrossRef Full Text | Google Scholar

Harrell, M. D., Harbi, S., Hoffman, J. F., Zavadil, J., and Coetzee, W. A. (2007). Large-scale analysis of ion channel gene expression in the mouse heart during perinatal development. Physiol. Genomics 28, 273–283. doi: 10.1152/physiolgenomics.00163.2006

PubMed Abstract | CrossRef Full Text | Google Scholar

Hichri, E., Abriel, H., and Kucera, J. P. (2018). Distribution of cardiac sodium channels in clusters potentiates ephaptic interactions in the intercalated disc. J. Physiol. 596, 563–589. doi: 10.1113/JP275351

PubMed Abstract | CrossRef Full Text | Google Scholar

Hirschy, A., Schatzmann, F., Ehler, E., and Perriard, J.-C. (2006). Establishment of cardiac cytoarchitecture in the developing mouse heart. Dev. Biol. 289, 430–441. doi: 10.1016/j.ydbio.2005.10.046

PubMed Abstract | CrossRef Full Text | Google Scholar

Jæger, K. H., Edwards, A. G., McCulloch, A., and Tveito, A. (2019). Properties of cardiac conduction in a cell-based computational model. PLoS Comput. Biol. 15:e1007042. doi: 10.1371/journal.pcbi.1007042

PubMed Abstract | CrossRef Full Text | Google Scholar

Kato, Y., Masumiya, H., Agata, N., Tanaka, H., and Shigenobu, K. (1996). Developmental changes in action potential and membrane currents in fetal, neonatal and adult guinea-pig ventricular myocytes. J. Mol. Cell. Cardiol. 28, 1515–1522. doi: 10.1006/jmcc.1996.0141

PubMed Abstract | CrossRef Full Text | Google Scholar

King, D. R., Entz, M., Blair, G. A., Crandell, I., Hanlon, A. L., Lin, J., et al. (2021). The conduction velocity-potassium relationship in the heart is modulated by sodium and calcium. Pflügers Arch. 473, 557–571. doi: 10.1007/s00424-021-02537-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Kucera, J. P., Rohr, S., and Rudy, Y. (2002). Localization of sodium channels in intercalated disks modulates cardiac conduction. Circ. Res. 91, 1176–1182. doi: 10.1161/01.RES.0000046237.54156.0A

PubMed Abstract | CrossRef Full Text | Google Scholar

Kutyifa, V., Daimee, U. A., McNitt, S., Polonsky, B., Lowenstein, C., Cutter, K., et al. (2018). Clinical aspects of the three major genetic forms of long QT syndrome (LQT 1, LQT 2, LQT 3). Ann. Noninvasive Electrocardiol. 23:e12537. doi: 10.1111/anec.12537

PubMed Abstract | CrossRef Full Text | Google Scholar

Kwak, B. R., and Jongsma, H. J. (1996). Regulation of cardiac gap junction channel permeability and conductance by several phosphorylating conditions. Mol. Cell. Biochem. 157, 93–99. doi: 10.1007/978-1-4613-1275-8_11

PubMed Abstract | CrossRef Full Text | Google Scholar

Leo-Macias, A., Agullo-Pascual, E., Sanchez-Alonso, J. L., Keegan, S., Lin, X., Arcos, T., et al. (2016). Nanoscale visualization of functional adhesion/excitability nodes at the intercalated disc. Nat. Commun. 7, 1–12. doi: 10.1038/ncomms10919

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, J., and Keener, J. P. (2010). Modeling electrical activity of myocardial cells incorporating the effects of ephaptic coupling. Proc. Natl. Acad. Sci. U.S.A. 107, 20935–20940. doi: 10.1073/pnas.1010154107

PubMed Abstract | CrossRef Full Text | Google Scholar

Livshitz, L. M., and Rudy, Y. (2007). Regulation of Ca2+ and electrical alternans in cardiac myocytes: role of camkii and repolarizing currents. Am. J. Physiol. Heart Circ. Physiol. 292, H2854–H2866. doi: 10.1152/ajpheart.01347.2006

PubMed Abstract | CrossRef Full Text | Google Scholar

McCain, M. L., Desplantez, T., Geisse, N. A., Rothen-Rutishauser, B., Oberer, H., Parker, K. K., et al. (2012). Cell-to-cell coupling in engineered pairs of rat ventricular cardiomyocytes: relation between cx43 immunofluorescence and intercellular electrical conductance. Am. J. Physiol. Heart Circ. Physiol. 302, H443–H450. doi: 10.1152/ajpheart.01218.2010

PubMed Abstract | CrossRef Full Text | Google Scholar

Mezache, L., Struckman, H. L., Greer-Short, A., Baine, S., Györke, S., Radwański, P. B., et al. (2020). Vascular endothelial growth factor promotes atrial arrhythmias by inducing acute intercalated disk remodeling. Sci. Rep. 10, 1–14. doi: 10.1038/s41598-020-77562-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Milman, A., Andorin, A., Gourraud, J.-B., Sacher, F., Mabo, P., Kim, S.-H., et al. (2017). Age of first arrhythmic event in brugada syndrome: data from the sabrus (survey on arrhythmic events in brugada syndrome) in 678 patients. Circ. Arrhythmia Electrophysiol. 10:e005222. doi: 10.1161/CIRCEP.117.005222

PubMed Abstract | CrossRef Full Text | Google Scholar

Moise, N., Struckman, H. L., Dagher, C., Veeraraghavan, R., and Weinberg, S. H. (2021). Intercalated disk nanoscale structure regulates cardiac conduction. J. Gen. Physiol. 153:e202112897. doi: 10.1085/jgp.202112897

PubMed Abstract | CrossRef Full Text | Google Scholar

Moreno, A., Rook, M., Fishman, G., and Spray, D. (1994). Gap junction channels: distinct voltage-sensitive and-insensitive conductance states. Biophys. J. 67, 113–119. doi: 10.1016/S0006-3495(94)80460-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Mori, Y., Fishman, G. I., and Peskin, C. S. (2008). Ephaptic conduction in a cardiac strand model with 3d electrodiffusion. Proc. Natl. Acad. Sci. U.S.A. 105, 6463–6468. doi: 10.1073/pnas.0801089105

PubMed Abstract | CrossRef Full Text | Google Scholar

Nielsen, M. S., Nygaard Axelsen, L., Sorgen, P. L., Verma, V., Delmar, M., and Holstein-Rathlou, N.-H. (2012). Gap junctions. Comprehens. Physiol. 2, 1981–2035. doi: 10.1002/cphy.c110051

PubMed Abstract | CrossRef Full Text | Google Scholar

Nowak, M. B., Greer-Short, A., Wan, X., Wu, X., Deschênes, I., Weinberg, S. H., et al. (2020). Intercellular sodium regulates repolarization in cardiac tissue with sodium channel gain-of-function. Biophys. J. 118:2829–2843. doi: 10.1016/j.bpj.2020.04.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Nowak, M. B., Poelzing, S., and Weinberg, S. H. (2021). Mechanisms underlying age-associated manifestation of cardiac sodium channel gain-of-function. J. Mol. Cell. Cardiol. 153, 60–71. doi: 10.1016/j.yjmcc.2020.12.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Peters, N. S., Coromilas, J., Severs, N. J., and Wit, A. L. (1997). Disturbed connexin43 gap junction distribution correlates with the location of reentrant circuits in the epicardial border zone of healing canine infarcts that cause ventricular tachycardia. Circulation 95, 988–996. doi: 10.1161/01.CIR.95.4.988

PubMed Abstract | CrossRef Full Text | Google Scholar

Peters, N. S., Severs, N. J., Rothery, S. M., Lincoln, C., Yacoub, M. H., and Green, C. R. (1994). Spatiotemporal relation between gap junctions and fascia adherens junctions during postnatal development of human ventricular myocardium. Circulation 90, 713–725. doi: 10.1161/01.CIR.90.2.713

PubMed Abstract | CrossRef Full Text | Google Scholar

Poelzing, S., and Rosenbaum, D. S. (2004). Altered connexin43 expression produces arrhythmia substrate in heart failure. Am. J. Physiol. Heart Circ. Physiol. 287, H1762–H1770. doi: 10.1152/ajpheart.00346.2004

PubMed Abstract | CrossRef Full Text | Google Scholar

Quan, W., and Rudy, Y. (1990). Unidirectional block and reentry of cardiac excitation: a model study. Circ. Res. 66, 367–382. doi: 10.1161/01.RES.66.2.367

PubMed Abstract | CrossRef Full Text | Google Scholar

Rhett, J. M., Ongstad, E. L., Jourdan, J., and Gourdie, R. G. (2012). Cx43 associates with na v 1.5 in the cardiomyocyte perinexus. J. Memb. Biol. 245, 411–422. doi: 10.1007/s00232-012-9465-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Rohr, S., Kucera, J. P., and Kléber, A. G. (1998). Slow conduction in cardiac tissue, I: effects of a reduction of excitability versus a reduction of electrical coupling on microconduction. Circ. Res. 83, 781–794. doi: 10.1161/01.RES.83.8.781

PubMed Abstract | CrossRef Full Text | Google Scholar

Rosen, M. R., Legato, M. J., and Weiss, R. M. (1981). Developmental changes in impulse conduction in the canine heart. Am. J. Physiol. Heart Circ. Physiol. 240, H546–H554. doi: 10.1152/ajpheart.1981.240.4.H546

PubMed Abstract | CrossRef Full Text | Google Scholar

Shaw, R. M., and Rudy, Y. (1997). Ionic mechanisms of propagation in cardiac tissue: roles of the sodium and l-type calcium currents during reduced excitability and decreased gap junction coupling. Circ. Res. 81, 727–741. doi: 10.1161/01.RES.81.5.727

PubMed Abstract | CrossRef Full Text | Google Scholar

Smith, J., Green, C., Peters, N., Rothery, S., and Severs, N. (1991). Altered patterns of gap junction distribution in ischemic heart disease. An immunohistochemical study of human myocardium using laser scanning confocal microscopy. Am. J. Pathol. 139:801.

PubMed Abstract | Google Scholar

Spach, M. S., Heidlage, J. F., Dolber, P. C., and Barr, R. C. (2000). Electrophysiological effects of remodeling cardiac gap junctions and cell size: experimental and model studies of normal cardiac growth. Circ. Res. 86, 302–311. doi: 10.1161/01.RES.86.3.302

PubMed Abstract | CrossRef Full Text | Google Scholar

Sperelakis, N. (2002). An electric field mechanism for transmission of excitation between myocardial cells. Circ. Res. 91, 985–987. doi: 10.1161/01.RES.0000045656.34731.6D

PubMed Abstract | CrossRef Full Text | Google Scholar

Swift, L. M., Burke, M., Guerrelli, D., Reilly, M., Ramadan, M., McCullough, D., et al. (2020). Age-dependent changes in electrophysiology and calcium handling: implications for pediatric cardiac research. Am. J. Physiol. Heart Circ. Physiol. 318, H354–H365. doi: 10.1152/ajpheart.00521.2019

PubMed Abstract | CrossRef Full Text | Google Scholar

Tveito, A., Jæger, K. H., Kuchta, M., Mardal, K.-A., and Rognes, M. E. (2017). A cell-based framework for numerical modeling of electrical conduction in cardiac tissue. Front. Phys. 5:48. doi: 10.3389/fphy.2017.00048

CrossRef Full Text | Google Scholar

Valiunas, V., Beyer, E. C., and Brink, P. R. (2002). Cardiac gap junction channels show quantitative differences in selectivity. Circ. Res. 91, 104–111. doi: 10.1161/01.RES.0000025638.24255.AA

PubMed Abstract | CrossRef Full Text | Google Scholar

Veeraraghavan, R., and Gourdie, R. G. (2016). Stochastic optical reconstruction microscopy-based relative localization analysis (STORM-RLA) for quantitative nanoscale assessment of spatial protein organization. Mol. Biol. Cell 27, 3583–3590. doi: 10.1091/mbc.e16-02-0125

PubMed Abstract | CrossRef Full Text | Google Scholar

Veeraraghavan, R., Gourdie, R. G., and Poelzing, S. (2014). Mechanisms of cardiac conduction: a history of revisions. Am. J. Physiol. Heart Circ. Physiol. 306, H619–H627. doi: 10.1152/ajpheart.00760.2013

PubMed Abstract | CrossRef Full Text | Google Scholar

Veeraraghavan, R., Lin, J., Hoeker, G. S., Keener, J. P., Gourdie, R. G., and Poelzing, S. (2015). Sodium channels in the cx43 gap junction perinexus may constitute a cardiac ephapse: an experimental and modeling study. Pflügers Arch. 467, 2093–2105. doi: 10.1007/s00424-014-1675-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Verheule, S., Van Kempen, M. J., Welscher, P. H., Kwak, B. R., and Jongsma, H. J. (1997). Characterization of gap junction channels in adult rabbit atrial and ventricular myocardium. Circ. Res. 80, 673–681. doi: 10.1161/01.RES.80.5.673

PubMed Abstract | CrossRef Full Text | Google Scholar

Vignati, G. (2007). Pediatric arrhythmias: which are the news? J. Cardiovasc. Med. 8, 62–66. doi: 10.2459/01.JCM.0000247438.12817.9e

PubMed Abstract | CrossRef Full Text | Google Scholar

Vreeker, A., Van Stuijvenberg, L., Hund, T. J., Mohler, P. J., Nikkels, P. G., and Van Veen, T. A. (2014). Assembly of the cardiac intercalated disk during pre-and postnatal development of the human heart. PLoS ONE 9:e94722. doi: 10.1371/journal.pone.0094722

PubMed Abstract | CrossRef Full Text | Google Scholar

Wei, N., Mori, Y., and Tolkacheva, E. G. (2016). The dual effect of ephaptic coupling on cardiac conduction with heterogeneous expression of connexin 43. J. Theoret. Biol. 397, 103–114. doi: 10.1016/j.jtbi.2016.02.029

PubMed Abstract | CrossRef Full Text | Google Scholar

Wei, N., and Tolkacheva, E. G. (2020). Interplay between ephaptic coupling and complex geometry of border zone during acute myocardial ischemia: effect on arrhythmogeneity. Chaos 30:033111. doi: 10.1063/1.5134447

PubMed Abstract | CrossRef Full Text | Google Scholar

Weinberg, S. (2017). Ephaptic coupling rescues conduction failure in weakly coupled cardiac tissue with voltage-gated gap junctions. Chaos 27:093908. doi: 10.1063/1.4999602

PubMed Abstract | CrossRef Full Text | Google Scholar

Weingart, R. (1986). Electrical properties of the nexal membrane studied in rat ventricular cell pairs. J. Physiol. 370, 267–284. doi: 10.1113/jphysiol.1986.sp015934

PubMed Abstract | CrossRef Full Text | Google Scholar

White, R. L., Doeller, J., Verselis, V., and Wittenberg, B. (1990). Gap junctional conductance between pairs of ventricular myocytes is modulated synergistically by H+ and Ca++. J. Gen. Physiol. 95, 1061–1075. doi: 10.1085/jgp.95.6.1061

PubMed Abstract | CrossRef Full Text | Google Scholar

Wilde, A. A., Moss, A. J., Kaufman, E. S., Shimizu, W., Peterson, D. R., Benhorin, J., et al. (2016). Clinical aspects of type 3 long-QT syndrome: an international multicenter study. Circulation 134, 872–882. doi: 10.1161/CIRCULATIONAHA.116.021823

PubMed Abstract | CrossRef Full Text | Google Scholar

Wittenberg, B. A., White, R., Ginzberg, R. D., and Spray, D. C. (1986). Effect of calcium on the dissociation of the mature rat heart into individual and paired myocytes: electrical properties of cell pairs. Circ. Res. 59, 143–150. doi: 10.1161/01.RES.59.2.143

PubMed Abstract | CrossRef Full Text | Google Scholar

Yao, J.-A., Hussain, W., Patel, P., Peters, N. S., Boyden, P. A., and Wit, A. L. (2003). Remodeling of gap junctional channel function in epicardial border zone of healing canine infarcts. Circ. Res. 92, 437–443. doi: 10.1161/01.RES.0000059301.81035.06

PubMed Abstract | CrossRef Full Text | Google Scholar