Abstract
The atrioventricular node (AVN) is considered a “black box”, and the functioning of its dual pathways remains controversial and not fully understood. In contrast to numerous clinical studies, there are only a few mathematical models of the node. In this paper, we present a compact, computationally lightweight multi-functional rabbit AVN model based on the Aliev-Panfilov two-variable cardiac cell model. The one-dimensional AVN model includes fast (FP) and slow (SP) pathways, primary pacemaking in the sinoatrial node, and subsidiary pacemaking in the SP. To obtain the direction-dependent conduction properties of the AVN, together with gradients of intercellular coupling and cell refractoriness, we implemented the asymmetry of coupling between model cells. We hypothesized that the asymmetry can reflect some effects related to the complexity of the real 3D structure of AVN. In addition, the model is accompanied by a visualization of electrical conduction in the AVN, revealing the interaction between SP and FP in the form of ladder diagrams. The AVN model demonstrates broad functionality, including normal sinus rhythm, AVN automaticity, filtering of high-rate atrial rhythms during atrial fibrillation and atrial flutter with Wenckebach periodicity, direction-dependent properties, and realistic anterograde and retrograde conduction curves in the control case and the cases of FP and SP ablation. To show the validity of the proposed model, we compare the simulation results with the available experimental data. Despite its simplicity, the proposed model can be used both as a stand-alone module and as a part of complex three-dimensional atrial or whole heart simulation systems, and can help to understand some puzzling functions of AVN.
1 Introduction
The atrioventricular node (AVN) plays a key role in the cardiac electrical conduction system. It is located between the atria and ventricles, namely, in the base of the right atrium (Figure 1A). The AVN is the only site responsible for transmitting impulses originating in the sinoatrial node (SN) to the ventricles of the heart, coordinating the relationship between contraction periods of the heart chambers, and introducing a variable delay, allowing effective pumping of the blood in a wide range of cardiac rhythms (; ). Figure 1B demonstrates the detailed anatomical location of the rabbit AVN and surrounding tissues ().
FIGURE 1
The rabbit AVN contains two built-in functional pathways conducting impulses from the atria to His bundle (HB)—so-called fast pathway (FP) with a longer effective refractory period, and slow pathway (SP) with a shorter one (Figure 1B). The FP also has a smaller conduction latency than the SP (; ; ). The dual pathway AVN physiology provides protection of ventricles from atrial fibrillation (AFB) and supraventricular tachycardia by filtering excessive atrial electrical impulses (). On the other hand, the interplay of the FP and SP in the AVN can be a potential substrate for atrioventricular nodal re-entrant tachycardia (AVNRT) (; ). The AVN also acts as a subsidiary pacemaker and provides atrial and ventricular pacing in the case of SN failure. According to experimental studies, in such a case, the AVN pacemaker activity in the rabbit heart originates in the SP between the coronary sinus and the tricuspid valve in the inferior nodal extension ().
Despite numerous research efforts, many aspects of the AVN electrophysiological behavior remain unclear and debated (). In addition to animal experiments and clinical studies on humans, mathematical modeling at different levels of detalization and complexity may offer insight into the AVN “black box”. Currently, very few AVN models or atrial models incorporating AVN with dual conduction pathways exist (; ; ; ; ; ). The rabbit cardiac biophysically-detailed multicellular ion-channel model () is quite complex and requires reconfiguration, e.g., changing the number of model cells for simulation of different AVN phenomena. Moreover, some reproducibility issues were reported recently (). The model computes normal action potential propagation from the SN through the atrium to AVN, AVNRT, the effect of the FP and SP ablation, ICa,L block effect, the AVN filtering function during AFB, and AVN pacemaking.
There were also attempts to develop AVN models based on more computationally effective simplified non-linear two-variable models. In the work () the AVN model based on hundreds of Lienard-transformed modified Van der Pol () and FitzHugh-Nagumo (; ) equations, reproduces normal sinus rhythm and AVNRT. The three-dimensional anatomically-detailed model of the rabbit right atrium () includes the AVN with dual pathways and utilizes the cellular automaton method and modified FitzHugh-Nagumo two-variable models. The model simulates normal sinus rhythm and nodal behavior during AFB. In work () the human real-time hybrid automaton model was proposed, which can emulate sinus rhythm and the AVNRT. In the works (; ) the authors used a simple, functional mathematical model of the rabbit AVN with the inclusion of dual pathway physiology. This model consists of exponential mathematical equations and is based on rabbit heart preparation data, and is able to predict AVN conduction time and the interaction between the FP and SP wavefronts during regular and irregular atrial rhythms. The functions of the models include simulation of nodal recovery, concealed conduction, Wenckebach phenomena, AFB and atrial flutter (AFL), and FP and SP ablations. The network model of human AVN () consists of FP and SP interconnected at one end, which includes functions of the HB and Purkinje network. Both pathways are composed of interacting nodes, and the refractory periods and conduction delays of the latter are determined by the stimulation history of each node and described by exponential functions (). The model is able to represent RR interval series extracted from ECG data (both in the forms of histograms, Poincaré plots, and autocorrelation), and the latest version also incorporates autonomic tone ().
In our previous work (), we considered a simplified human AVN model with dual pathways, which consists of about 40 cells including excitable cells and pacemaking cells described by Aliev-Panfilov (A-P) () and modified Van der Pol models. The model reproduces normal sinus rhythm, AVN automaticity, reentry, and filtering function. However, we found that various variants of the modified Van der Pol equation are not suitable for the description of cells of subsidiary pacemakers.
In this work, we present a compact mathematical model of the rabbit AVN based on a single non-linear two-variable A-P model for the description of both quiescent excitable and pacemaking cells, which allows representing various parts of cardiac conduction system with different tissue types in a uniform and convenient way. The model consists of only 33 cells and requires relatively low computational resources. Adjusting to available experimental data for the rabbit heart made it possible to demonstrate a wide range of behavioral phenomena of AVN. This includes proper bidirectional conduction in the FP and SP both in the normal case and after the FP and SP ablations reflected in the calculated recovery curves, AVNRT, the AVN automaticity, filtering function during AFB and AFL, and Wenckebach periodicity. To achieve correct reproduction of the retrograde conduction characteristics, we assumed the presence of an asymmetry of the coupling between the AVN cells and implemented it in the model. In addition, we included visualization demonstrating the inner processes in the AVN “black box”, in particular, the laddergrams (Lewis ladder diagrams) () with the exact timing of each model cell in both FP and SP pathways during the excitation wave propagation. Finally, the comparison with published experimental findings provides validation of the simulated results.
2 Materials and methods
2.1 Structure and equations of the model
A schematic view of the one-dimensional rabbit cardiac conduction system model incorporating AVN with dual pathways is shown in Figure 2. The model consists of a two-row matrix, where the upper row includes the SN and peripheral sinus node (PS1–PS3) cells, atrial muscle (AM1–AM3), fast pathway (FP1–FP8), penetrating bundle (PB), and HB (HB1-HB6) cells, while the lower row includes only ten slow pathway (SP1–SP10) cells and additional intermediate atrial muscle cell (AM*). The latter is added to provide a smooth transformation of short rectangular-shaped excitation impulses into typical atrial action potential shapes in the case of AFB (with random pacing) and AFL (regular pacing) simulations. The arrows Antero and Retro denote the points of application of external S1 and S2 stimuli for anterograde and retrograde propagation study with S1S2 stimulation protocol, and the arrow AF indicates entrance for “internal” stimuli in the cases of AFB and AFL simulations. The premature S1S2 stimulation () refers to pacing that is performed at a constant basic cycle length (shorter than normal sinus rhythm) for 8–10 S1 impulses, followed by the introduction of an extra stimulus S2, and with stepwise reduction of the S1S2 interval until AVN refractoriness (total conduction block) is reached. Therefore, S1 and S2 mark the basic and test stimulus, respectively, and the same subscripts apply to atrial (A) and HB (H) stimulation or response.
FIGURE 2
Each model cell represents a group of real cardiac cells and is described by two-variable A-P model. The model for monodomain description of the cardiac cell consists of the following set of reaction-diffusion type non-linear ordinary differential equations:Where 0 ≤ u ≤ 1 and v are the dimensionless normalized transmembrane potential and slow recovery variables, respectively. Parameter k controls the magnitude of transmembrane current, and parameters μ1, μ2, a1, and a2 are adjusted to reproduce characteristics of cardiac tissue, ɛ sets the time scale of the recovery process determining the restitution properties of the action potential, ct is the time scaling coefficient introducing physical time in seconds into the system. In the general case, a1 = a2 > 0. Realistic transmembrane voltage values of the AVN cells in millivolts can be obtained with the rescaling of the u value by 120u − 80 mV for the AM, and 100u − 70 mV for HB (). Since the A-P model was originally developed for the description of canine ventricular myocyte, the parameters were adjusted for the rabbit AVN in terms of action potential duration and refractoriness (see Supplementary Tables S1, S2). Istim is the external stimulation current, Icoupl = ∇ ⋅ (D∇u) is the coupling (diffusion) current, where ∇ is a spatial gradient operator defined within the model tissue geometry, and D is a tensor of diffusion coefficients characterizing electrotonic interactions between neighboring cells via gap junctional coupling conductance. For the considered one-dimensional system, the coupling term in the discrete form reads as followswhere i is the model cell index, and d is the coupling (diffusion) coefficient in s−1 assigned for each pair of cells normalized on dimensionless distance. As will be discussed below, to enhance the direction-dependent conduction properties of the AVN model, we implemented a coupling asymmetry with the modified coupling term:where α > 0 is the coefficient of asymmetry applied to the right-hand side cell in each coupled pair of model cells. The values α = 1 correspond to symmetric coupling, α < 1 to the accelerated anterograde and slowed retrograde conduction, and α > 1 to the accelerated retrograde and slowed anterograde conduction.
The cells playing the role of natural primary (SN) and subsidiary (SP7 through SP10, FP8, and PB) pacemakers are gray-shaded on the scheme (Figure 2). They are described by the same A-P model with a negative coefficient a1. The oscillating (pacemaking) properties of the A-P model were considered recently (). Implementing the oscillating A-P model cells simplifies the AVN model description, providing easy control of intrinsic oscillation frequency and reliable overdrive suppression. That substitutes a significant advantage of the cardiac A-P model over the FitzHugh-Nagumo and modified Van der Pol models and their combinations, used in the modeling of the simplified cardiac conduction system (; ), due to the absence of the quiescent state in the modified Van der Pol model and the absence of well-defined threshold potential in both FitzHugh-Nagumo and modified Van der Pol models. Also, close values of the minimal diastolic potential of both non-pacemaking and pacemaking A-P model cells allow the driving of a relatively big number of non-pacemaking cells by a single pacemaking A-P model cell ().
2.2 Parameters
As a starting point for the search of optimal model cell parameters, we took the experimental data from (; ): minimal and maximal AVN conduction time for normal case and after the FP and SP ablations for anterograde conduction and retrograde conduction, the nodal effective refractory period (ERPN) for control case and after the SP ablation, and the exponental-like shapes of conduction curves for both anterograde and retrograde directions.
Taking into account prominent difference in antrograde and retrograde conduction curves observed experimentally, we assumed that the characteristics of FP and SP should include non-linear conduction properties (some degree of rectification, e.g., diode-like conduction) due to spatial gradients of refractory period, local intercellular coupling, coefficient of coupling asymmetry, and action potential duration. For the sake of simplicity, in the course of the model development, we paid attention to the differences between the FP and SP local conduction delays due to spatial distribution of coupling coefficient but not the macroscopic conduction velocities. Faster conduction velocity in the FP is not confirmed in experiments (). Anatomically shorter FP distance compared with the SP is considered a primary cause of a shorter conduction delay in the FP (). Thus, though in our phenomenological model, the SP is longer than the FP in terms of the number of model cells (10 vs. 8), this difference in the model cell numbers may not correspond to the real anatomical structure of the rabbit AVN. Proper pathway conduction delays were obtained by setting spatial distributions of the coupling coefficient d and coupling asymmetry α values, and refractory periods of the model cells (see Figure 3 and Supplementary Tables S1, S2). In this regard, in our model, the cells represent groups of real cells and are distinguished not by biophysical types (; ; ) but by a spatial gradual distribution of the parameters. Figure 3 shows the distribution of refractoriness in the uncoupled model cells, local coupling coefficient d, and coefficient of coupling asymmetry α. The refractory period values are in line with the available experimental data: 166 ± 30 ms for the rabbit SN and 68 ± 11 ms for atrial muscle (), 91 ± 10 ms for the inferior nodal extension (), and 141 ± 15 ms for the transitional zone (), which is an anatomical substrate of the FP. The profile of d approximately corresponds to the distribution of coupling conductance in the ion-channel multicellular model (). Action potential profiles of the uncoupled model cells in the FP and SP can be found in Supplementary Figure S1.
FIGURE 3
The FP and SP ablations were simulated by setting diffusion coefficient d = 0 in the middle of the pathways (between FP5 and FP6, and SP6 and SP7 cells).
2.3 External pacing
The S1S2 stimulation protocol consists of ten basic S1 impulses followed by a premature S2 impulse. We created the AVN conduction curves applying constant S1S1 basic cycles with 360 ms intervals (shorter than spontaneous sinus rhythm of 361 ms) followed by the test premature S1S2 cycle with stepwise reduction starting from the basic cycle length until the full AVN functional refractory block occurs (360–90 ms) (). The current stimulation impulses were of 1 ms duration and 1.3× greater than the threshold value (Istim = 280.0). The AFB and AFL were simulated by applying similar stimulation impulses through the intermediate cell AM* to eliminate the influence of the pacing impulse shape. The correct shape of AM* action potentials allowed more natural stimulation of the atrial part of the conduction system model. For the simulation of the AFB, a random sequence of the impulses with the same amplitude and duration was generated using MATLAB pseudorandom number generator with uniform distribution in the 75–125 ms range.
2.4 Visualization of laddergrams
An original algorithm was developed to visualize the laddergrams with the exact activation timing of each model cell. It allows tracing the wavefront propagation for each excitation sequence and in each pathway separately (with red color for the FP, and blue color for the SP) based on the activation time in each cell. The activation times in each model cell were recorded at the level of 0.25 during the action potential upstroke.
For anterograde propagation, the algorithm determines the leading pathway which provides the excitation to reach the PB first, thus, the traces coming down to the HB are marked in the same color. This is very useful since it allows visual differentiation between conduction through the SP and FP obtained in the simulations and easy comparison with the experimental method (; ), in which pronounced change in the signal amplitude in the inferior HB domain occurs when conduction switches from one pathway to another, and is especially important for visualization of chaotic rhythms ().
2.5 Numerical methods
The solution of the set of Eqs. (1) and (2) was performed with MATLAB (R2022a) standard ODE solver ode23 with relative tolerance 10–7 and absolute tolerance 10–10. The ode23 is a single-step solver and implements the explicit Runge-Kutta (2,3) pair (). Since ode23 uses an adaptive time stepping, for the reliable periodic and random pacing, MATLAB mod function was used to force the solver to accept input impulses at predefined time moments. The solver output with the variable time steps was evaluated with the finite discretization of 0.1 ms for further analysis. Open-end boundary conditions were used at both SN and HB6 ends of the model.
The MATLAB source code for this study is available in Figshare repository (https://doi.org/10.6084/m9.figshare.22015289.v1).
3 Results and their validation
3.1 Remarks on figures of excitation propagation
Figures 4–8 demonstrate different cases of simulated excitation propagation in the AVN. Each figure consists of a scheme and three panels. The scheme illustrates the propagation conduction sequence of interacting pathways in each case. In Figures 5–8 the scheme corresponds to the response on S2 stimulus. The top panel shows the laddergram with exact timing of each model cell excitation. The trace color in the HB part corresponds to the leading pathway. The middle and bottom panels demonstrate the sequences of action potentials passing between the SN and HB through the FP and SP, respectively.
FIGURE 4
FIGURE 5
FIGURE 6
FIGURE 7
FIGURE 8
3.2 Normal sinus rhythm
First, normal sinus rhythm with a cycle length of 361 ms (166 bpm) which is typical for rabbit heart (363 ± 21 ms, ()) was simulated, and the results are shown in Figure 4A. The numbers in the top panel indicate excitation wave latency with respect to the SN and the sinus rhythm in milliseconds. The excitation wave initiated at the SN propagates along the peripheral sinus node and AM cells and reaches the last atrium point (AM3) at 43.7 ms. Then the wave propagates through both the FP and SP. The wave in the FP arrives first at the PB cell at 92.4 ms, causing retrograde propagation through the SP. Both SP anterograde and retrograde excitation fronts meet in the inferior nodal extension at 99.8 ms and cancel each other. The leading FP conduction wave travels further down from the PB to the last HB cell of the model and arrives there at 139.6 ms. The obtained latency values are very close to the experimental measurements () shown in Figure 1B for the right atrium of the rabbit heart and computer simulations with complex 1D ion-channel based rabbit AVN model (). A similar scheme of AVN conduction via the FP and SP was also observed in the case of atrial pacing at long S1S2 intervals (see Figure 5A below).
3.3 AVN automaticity
In addition to its essential role in controlling the electrical impulse transmission between the atria and ventricles, the AVN can also serve as a subsidiary pacemaker when the SN fails (). In rabbits, the AVN leading pacemaker site is defined in the inferior nodal extension, generating electrical excitation with an average cycle length of about 600 ms (627 ± 111 ms, ()). In this case, the excitation spreads in two directions - upward to the atria and downward to the HB. Figure 4B shows the laddergram of AVN automaticity and the corresponding propagating action potentials in the absence of SN activity. The SP is the leading pathway for both anterograde and retrograde directions.
3.4 Atrial pacing
The most important properties of the AVN behavior are reflected in the so-called AV nodal conduction or recovery curve (). The characteristic steps of the AVN anterograde conduction curve calculation during atrial pacing are presented in Figures 5, 6. In the range of long S1S2 intervals, the conduction scheme is similar to that demonstrated in Figure 4A with the FP leading.
At the S1S2 interval shortened to approximately 140 ms, a transitional stage occurs when the FP and SP excitation waves reach the PB almost simultaneously with minimal difference in latencies (Figure 5A). In such a case, while our model determines the leading pathway mathematically, from a physiological point of view both pathways provide conduction down to the HB. Simultaneous input of both pathways was experimentally demonstrated in HB electrograms with intermediate height ().
For the S1S2 intervals shorter than 135 ms, only the SP becomes functional and allows passage of the excitation from the SN to the HB (Figure 5B).
By reducing the interval below 123 ms, we observed the onset of perpetual AVN reentrant tachycardia (AVNRT) (; ), when excitation repeatedly circulates within the AVN ring (Figure 6A). The existence of AVNRT in the rabbit AVN at short S1S2 intervals was reported in previous experiments (). An excitation wave propagating anterogradely along the SP, retrogradely travels back up the FP, and then activates the atrial muscle, after which subsequent reentry cycles occur. This phenomenon takes place because the FP has a significantly longer refractory period than the SP (Figure 3). When excitation travels via the SP, the upper part of the FP is blocked, but the lower FP part recovers its excitability and conducts in a retrograde direction (). According to the AVNRT classification (), the observed AVNRT belongs to the typical slow-fast type. AVNRT persisted with further shortening of S1S2 intervals until full functional block was achieved. In our simulations, the full functional block occurred at S1S2 intervals shorter than 95 ms when there was no conduction through both pathways (Figure 6B).
In all the above cases, the excitation wave from the SN pacemaker meets and annihilates with the retrograde wave from the atrial muscle caused by premature S1 impulses, and the SN activity is suppressed by S2 impulses.
Animation of S1S2 protocol simulation for anterograde atrial pacing is shown in Supplementary Video S1.
3.5 Effects of FP and SP ablations on anterograde conduction
Catheter ablation of the heart is used to destroy small areas in one of the AVN pathways, thereby stopping the propagation of unusual or undesirable electrical signals that travel through the AVN and cause irregular arrhythmias (). Figure 7 shows the simulation results of the atrial pacing after the SP and FP ablations. In the former case (Figure 7A), when S1S2 impulses are delivered at 110 ms intervals, the conduction of excitation is totally canceled since, at such short pacing intervals, functional conduction block in the FP takes place. Therefore, due to the broken AVN ring, SP ablation eliminates AVNRT shown in Figure 6A ().
FP ablation (Figure 7B) leads to increased conduction time (). Although the response to the S2 impulse looks similar to that shown in Figure 5B, where FP conduction is absent due to the functional conduction block associated with FP refractoriness, the increase in the conduction time is due to a different conduction propagation scenario caused by the preceding S1 impulse. Because ablation destroys pathway tissue, this conduction scenario extends over the entire range of S1S2 intervals (see complete anterograde conduction curves in Figure 9).
FIGURE 9
3.6 Retrograde conduction and FP and SP ablations
Retrograde conduction refers to conduction from the lower part of the cardiac conduction system to the atria through the AVN. Experimental and clinical studies have demonstrated the ability of the AV node to conduct impulses in the anterograde and retrograde directions (; ; ; ).
We simulated retrograde conduction applying S1S2 (H1H2) stimulation to the HB (at HB1 model cell, Retro arrow in Figure 2) and plotted pertinent recovery curves (see AVN conduction curves Subsection). Figure 8 demonstrates retrograde conduction in the control case (Figure 8A), after FP ablation (Figure 8B), and after SP ablation (Figure 8C) for S1S2 interval of 175 ms.
In the control and SP ablation cases, excitation propagates retrogradely via the FP, creating almost the same A1A2 intervals and with close H1A1 and H2A2 conduction times. In the former case, similar to the case of atrial pacing with long intervals, retrograde and anterograde waves in SP meet and annihilate, but in the upper part of the AVN model (Figure 8A). In the FP ablation case, the A1A2 interval was shorter (Figure 8B), while conduction passing through SP leads to longer H1A1 and H2A2 latencies (see complete retrograde conduction curves in Figure 9).
Animation of S1S2 protocol simulation for retrograde His bundle pacing in the control case is shown in Supplementary Video S2.
3.7 AVN conduction curves
Applying the S1S2 stimulation protocol to the atrial and HB cells, we constructed various conduction curves for control (pre-ablation) and post-ablation cases. The recovery curve represents the atrial-His conduction time A2H2 for anterograde conduction or His-atrial conduction time H2A2 for retrograde conduction plotted either versus the S1S2 interval (A1A2 for atrial pacing and H1H2 for HB pacing) or versus the recovery interval (H1A2 for atrial pacing and A1H2 for HB pacing) ().
According to the experimental findings, the AVN conduction curves demonstrate pronounced direction-dependent behavior (; ). Figure 9 shows the comparison of experimentally obtained (left panels) (modified from ()) and simulated (right panels) anterograde and retrograde recovery curves in the control and ablation cases plotted versus S1S2 intervals. The modification from the dependence on recovery time H1A2 to S1S2 (A1A2) interval was performed by shifting each conduction curve on its minimal conduction time (A1H1) value according to formula A1A2 = H1A2 + A1H1 (; ).
Simulated anterograde recovery curve (Figure 9A, right panel) in the control case demonstrates a typical exponential-like rise with reducing pacing interval (; ; ) with a pronounced bend and tilting at the point where conduction switches from the FP to SP with smooth change of the conduction curve from flat to steep without gap, which is typical for rabbit AVN (; ; ). At long S1S2 intervals, the SP conduction is overcome by the FP conduction (flat portion of the control curve, see also Figure 5A), while at short S1S2 intervals, excitation propagates through the SP only (Figures 5B, 6A) due to longer refractory period of the FP (Figure 3). In the S1S2 range of 95–122 ms, the perpetual AVNRT was observed (), as mentioned above (Figure 6A).
The ERPN is set as the longest S1S2 pacing interval that results in neither HB nor atrial response in the cases of atrial pacing (anterograde propagation) and HB pacing (retrograde propagation), respectively. FP ablation markedly increased conduction time in both anterograde (; ) and retrograde () conduction cases. The ERPN was 94 ms for anterograde conduction and did not change after FP ablation (91 ± 10 ms in ()). For the retrograde conduction, the ERPN was 163 ms and also did not change after FP ablation (143 ± 16 ms before and 149 ± 14 ms after FP ablation in ()).
During our initial simulations, the value of the upward shift of the retrograde control curve was relatively small (dashed line in Figure 9), while the experiments () demonstrated much slower retrograde conduction in the control case. The introduction of the asymmetry to the intercellular coupling in the FP, PB, and upper part of HB allowed additional slowing of the retrograde control conduction, better reflecting the AVN direction-dependent conduction properties ().
SP ablation blocked its conduction, while FP conduction remained unaffected. The ablation led to the prolongation of ERPN to 132 ms (141 ± 15 ms in ()) and the disappearance of the steep part of the anterograde conduction curve in the short range of S1–S2 intervals (Figure 9B). The ablation also eliminated the AVNRT, as observed in clinical practice () (compare Figures 6A, 7A).
In the retrograde case, SP ablation did not affect retrograde conduction and SP post-ablation curve virtually coincides with the control curve (Figure 9B). This, together with conduction scenarios shown in Figures 8A,C, confirms retrograde conduction through the FP in the control case. The ERPN was also prolonged to approximately the same value (163 ms) as in the retrograde pre-ablation case, matching with the experimental results (151 ± 12 ms in ()).
Figure 10 shows three main anterograde AVN functional curves (refractory curve H1H2, His-atrial curve H1A2, and recovery curve A2H2) obtained in experiments (left panels) and our simulations (right panels). The AVN functional refractory period (FRPN) is determined as the shortest interval between two consecutive impulses propagated from the atria to the HB (H1H2). In our simulations, FRPN was 180.5 ms, which is in agreement with experiments (175 ± 10 ms in ().
FIGURE 10
Each point of the refractory curve H1H2 obtained in our simulations equals the sum of A2H2 value on the recovery curve and H1A2 value on the His-atrial curve at the corresponding S1S2 interval (Figure 10A, right), which is in line with the previously proposed quantitative relationship between the three function curves () (Figure 10A, left). Figure 10B demonstrates another comparison of experimental () and simulated AVN refractory curves before and after SP ablation. SP ablation decreased the left steeply rising part of the refractory curve compared to the control curve, resulting in an increase in ERPN, but not FRPN.
Note that the experimental conduction curves shown in Figures 9, 10 were obtained for different rabbit preparations and slightly different experimental setups (; ; ). In particular, the AVN conduction measurements before and after FP and SP ablations were performed separately for each group of rabbits, thus making the experimental conduction curves rather different in Figures 9, 10. On the other hand, our model parameters were adjusted for the rabbit data corresponding to Figure 9A, thus our results show a greater quantitative difference with the experimental conduction curves demonstrated in Figures 9B, 10.
3.8 Filtering of atrial rhythm during atrial fibrillation
It has been demonstrated that dual pathways are involved in AVN conduction during AFB by limiting the number of atrial impulses transmitted to the ventricles (). To show the AVN filtering function and reveal the interaction between the pathways during AFB we set random atrial pacing with uniform distribution in the 75–125 ms range (Figure 11). In Figure 11, different conduction sequences for each atrial impulse are present, including complete passage of the excitation to HB through either pathway with the domination of the SP conduction over the FP (which may depend on mean random atrial rate), and simultaneous conduction block in both pathways. One can also see many episodes of retrograde propagation annihilating with an anterograde impulse resulting from a subsequent atrial stimulus. Besides, the sinus rhythm was completely suppressed by faster atrial pacing. This filtering behavior during AFB corresponds to the results observed in experiments (; ) and in simulations with complex 1D ion-channel based model (; ).
FIGURE 11
As an additional validation of the AVN model filtering function, we simulated the effect of SP and FP ablations on the AFB. Figure 12A shows laddergrams of SP ablation (top panel) and FP ablation (bottom panel) cases simulated with the same random atrial impulse pattern, as shown in Figure 11. The HB rate was visually lower for the SP ablation case. This is also confirmed by the statistical distributions obtained for all three cases simulated with random atrial stimulation during 350 s (Figure 12B). The reduction of the mean HB rate from 339 bpm in the control case to 292 bpm in the SP ablation case and to 318 bpm in the FP ablation case and prolongation of the mean HH interval from 177 ms to 205 ms and 189 ms, respectively, revealed SP ablation as a relatively more effective way of ventricular rate control than FP ablation, which is in agreement with experimental () and clinical () findings, and with mathematical modeling study ().
FIGURE 12
3.9 Wenckebach periodicity during atrial flutter
Dual pathway electrophysiology was shown to be also directly related to the occurrence of the Wenckebach phenomenon (; ). Figure 13 demonstrates an example of regular atrial pacing mimicking AFL that yields 5:4 Wenckebach periodicity corresponding to four out of five atrial impulses conducted to the HB, obtained at a regular atrial pacing interval of 125.2 ms (atrial/His rates 480/384 bpm). Figure 13A shows laddergram and action potential propagation through the FP and SP, while Figure 13B demonstrates AH conduction times (top panel) and HH intervals (bottom panel) of consecutive beats in the Wenckebach cycle. This “typical” Wenckebach pattern is characterized by the progressive increase of AVN conduction times (AH) and progressive shortening of HH intervals before the conduction block occurs (). Also, exactly in line with the experimentally obtained HB electrograms, the first atrial stimulus after the conduction block passed through the FP (depicted by red squares in Figure 13B), followed by a transition to the SP until its ultimate block (; ).
FIGURE 13
4 Discussion
4.1 Major achievements
In this work, we have developed a compact one-dimensional mathematical model of rabbit AVN extended structure, which includes SN and HB. The AVN model is based on the Aliev-Panfilov non-linear system of two-variable differential equations, describing both quiescent excitable and pacemaking cells. Despite its simplicity and a small number of model cells (33 elements), the model is multi-functional and does not require structural changes for reproduction of a wide variety of AVN dual-pathway electrophysiological phenomena. In particular, the model can correctly reproduce normal sinus rhythm, AVN subsidiary pacemaking, AVNRT, filtering AVN function during AFB and AFL with Wenckebach patterns, and situations after FP and SP ablations.
For the first time, we simulated the retrograde conduction in the AVN and implemented the asymmetry of the intercellular coupling in the FP, PB, and upper part of HB for accurate reproduction of AVN direction-dependent conduction properties (). Several conduction curves (Figures 9, 10) reflecting the general behavior of AVN were presented along with experimental data to demonstrate their good agreement.
Direction-dependent conduction can originate from specific 3D structural or functional characteristics (), such as the source-sink (current-to-load) mismatch (), or expression of connexins in heterotypic gap junctions (; ). The acceptable direction-dependent conduction behavior was achieved in the initial simulations with the spatial dispersion of action potential duration and gradients of local refractoriness and coupling coefficient (see Figure 3). However, the upward shift of the retrograde control curve (dashed line in Figure 9) was insufficient to reproduce the much slower retrograde conduction in the control case shown in experiments (). This led us to conclude that an additional factor is responsible for slowing down retrograde conduction in the FP. We hypothesized that the asymmetry of the intercellular coupling (Eq. (4)) plays the role of this factor of slow retrograde conduction. The gap junctions with asymmetric or rectifying electrical transmission properties are well known in neuroscience (). Rectifying coupling was also implemented in the simulation of the small intestine of mice (). The introduction of the asymmetry to the intercellular coupling in the FP, PB, and upper part of HB allowed additional slowing of the retrograde control conduction, better reflecting the AVN direction-dependent conduction properties (; ).
Our model is accompanied by a visualization tool, which discloses processes taking place inside the AVN “black box” by creating laddergrams and revealing the leading pathway for anterograde propagation. With the help of the laddergrams, it became possible to compare FP and SP functional state before and after ablation, including hidden SP activation over a broad range of atrial and HB pacing intervals, as well as the involvement of dual AVN pathways in the Wenckebach phenomena and during AFB.
We believe the proposed AVN model is suitable for cardiac research and educational purposes. The spatial dispersion of local refractoriness, coupling coefficient, and coupling asymmetry shown in Figure 3 can serve as a template for implementing ionic cellular models with different levels of complexity in the AVN modeling. Our model can be adjusted for other animal species and humans, taking into account specific features, for example, a characteristic gap (“jump”) in the conduction curve of human AVN (). The model can be integrated into the whole-heart models () and simulation frameworks such as openCARP () for near-realtime or even real-time simulations (). Another possible application of the model is utilization in the real-time testbed systems for implantable cardiac devices (; ). The proposed model can be also used for quick preliminary testing of hypotheses before accurate simulations with computationally expensive and biophysically detailed cellular models.
4.2 Comparison with previous AVN models
Similar to the complex multi-cellular ion-channel model () and functional mathematical model (), our model is based on experimental rabbit data. Due to the small number of elements and much simpler non-linear equation systems, our model allows faster computation and easier tuning than more biophysically detailed model (), keeping the model structure unchanged. On the other hand, functional mathematical models like () benefit from effortless fitting of exponential equations to particular rabbit cardiac preparation. The presented model is capable of reproducing almost all electrophysiological properties of AVN found in both types of models with the additional feature of retrograde conduction with realistic ERPN in control and FP and SP post-ablation cases.
4.3 Limitations
The concept of the presented model is based on monodomain formulation of electrical propagation in a 1D cable structure where FP and SP are electrically isolated from each other. It is possible that some macroscopic (anatomical) and microscopic (cytological and extracytologic matrix) properties () related to the complexity of real 3D AVN structure were not taken into account in the proposed AVN model. This can also explain the quantitative difference between experimental and simulated anterograde conduction curves in the range of short S1S2 pacing cycles and the retrograde conduction cases.
Normal conduction case shown in Figure 4A assumes that retrograde conduction exists in the SP. In Figure 4A, one can see that upon the excitation wave in the FP reaches the PB, the retrograde wave goes up through the SP, and the annihilation of the anterograde and retrograde waves in the SP occurs at about 100 ms. However, it was shown in the experiments () that in some rabbit heart preparations, retrograde conduction totally disappeared after the FP ablation, which means the absence of retrograde conduction in the intact SP. In such a case, normal conduction demonstrated in Figure 4A would look different for the SP. The SP trace will go straight down and end up at the PB due to the refractory state of the latter at about 110 ms. Similarly, in the case of spontaneous activity originating in the SP (Figure 4B), the retrograde part of the excitation wave in the SP will be absent, and retrograde conduction will be provided by the FP only.
Due to the simplicity of the two-variable A-P model used in this study, the proposed AVN model does not allow simulation of subtle effects on cardiac cell ion channels, such as the impact of drugs or change in ion concentrations as demonstrated in the more complex multicellular ion-channel based models ().
5 Conclusion
In this work, we have developed a compact multi-functional model of a rabbit atrioventricular node with dual pathways. The one-dimensional model is relatively simple, contains a small number of model cells, and is computationally efficient, allowing for fast simulations. Using the developed model, we successfully reproduced several electrophysiological phenomena accompanied by visualization of laddergrams. We believe the proposed AVN model has a high potential for utilization in research and education both as an independent module and as a part of complex three-dimensional whole-heart simulation systems and real-time testbed systems. The proposed model can also be used for quick preliminary testing of hypotheses before accurate and detailed simulations with computationally expensive and complex cellular models.
Statements
Data availability statement
The datasets presented in this study can be found in the article/Supplementary Material. The MATLAB source code is available in Figshare repository (https://doi.org/10.6084/m9.figshare.22015289.v1).
Author contributions
MR wrote the program code and performed computer simulations. ER analyzed the results. All authors contributed to the conception and design of the study and manuscript preparation and revisions, read and approved the submitted version.
Funding
This work was supported by Grant No. 20K12046, JSPS KAKENHI.
Acknowledgments
The authors would like to thank Katrina Armstrong for her assistance with the manuscript.
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.2023.1126648/full#supplementary-material
References
1
AiW.PatelN. D.RoopP. S.MalikA.AndalamS.YipE.et al (2018). A parametric computational model of the action potential of pacemaker cells. IEEE Trans. Biomed. Eng.65, 123–130. 10.1109/TBME.2017.2695537
2
AiW.PatelN. D.RoopP. S.MalikA.TrewM. L. (2020). Cardiac electrical modeling for closed-loop validation of implantable devices. IEEE Trans. Biomed. Eng.67, 536–544. 10.1109/TBME.2019.2917212
3
AlievR. R.PanfilovA. V. (1996). A simple two-variable model of cardiac excitation. Chaos Solit. Fractals7, 293–301. 10.1016/0960-0779(95)00089-5
4
BiasiN.SeghettiP.MercatiM.TognettiA. (2022). A reaction-diffusion heart model for the closed-loop evaluation of heart-pacemaker interaction. IEEE Access10, 121249–121260. 10.1109/access.2022.3222830
5
BilletteJ. (1987). Atrioventricular nodal activation during periodic premature stimulation of the atrium. Am. J. Physiol.252, H163–H177. 10.1152/ajpheart.1987.252.1.H163
6
BilletteJ.TadrosR. (2019). An integrated overview of AV node physiology. Pacing Clin. Electrophysiol.42, 805–820. 10.1111/pace.13734
7
CarmelietE. (2019). Conduction in cardiac tissue. Historical reflections. Physiol. Rep.7, e13860. 10.14814/phy2.13860
8
ClimentA. M.GuillemM. S.ZhangY.MilletJ.MazgalevT. N. (2011). Functional mathematical model of dual pathway AV nodal conduction. Am. J. Physiol. Heart Circ. Physiol.300, H1393–H1401. 10.1152/ajpheart.01175.2010
9
de CarvalhoA. P.de MelloW. C.HoffmanB. F. (1959). Electrophysiological evidence for specialized fiber types in rabbit atrium. Am. J. Physiol. Leg. Content196, 483–488. 10.1152/ajplegacy.1959.196.3.483
10
DenesP.WuD.DhingraR. C.ChuquimiaR.RosenK. M. (1973). Demonstration of dual A-V nodal pathways in patients with paroxysmal supraventricular tachycardia. Circulation48, 549–555. 10.1161/01.cir.48.3.549
11
DobrzynskiH.NikolskiV. P.SambelashviliA. T.GreenerI. D.YamamotoM.BoyettM. R.et al (2003). Site of origin and molecular substrate of atrioventricular junctional rhythm in the rabbit heart. Circ. Res.93, 1102–1110. 10.1161/01.RES.0000101913.95604.B9
12
FeldG. K.FleckR. P.FujimuraO.ProthroD. L.BahnsonT. D.IbarraM. (1994). Control of rapid ventricular response by radiofrequency catheter modification of the atrioventricular node in patients with medically refractory atrial fibrillation. Circulation90, 2299–2307. 10.1161/01.cir.90.5.2299
13
FitzHughR. (1961). Impulses and physiological states in theoretical models of nerve membrane. Biophys. J.1, 445–466. 10.1016/s0006-3495(61)86902-6
14
GarrigueS.TchouP. J.MazgalevT. N. (1999). Role of the differential bombardment of atrial inputs to the atrioventricular node as a factor influencing ventricular rate during high atrial rate. Cardiovasc. Res.44, 344–355. 10.1016/s0008-6363(99)00201-1
15
GeorgeS. A.FayeN. R.Murillo-BerliozA.LeeK. B.TrachiotisG. D.EfimovI. R. (2017). At the atrioventricular crossroads: Dual pathway electrophysiology in the atrioventricular node and its underlying heterogeneities. Arrhythmia Electrophysiol. Rev.6, 179–185. 10.15420/aer.2017.30.1
16
GilletteK.GsellM.StrocchiM.GranditsT.NeicA.ManningerM.et al (2022). A personalized real-time virtual model of whole heart electrophysiology. Front. Physiol.13, 907190. 10.3389/fphys.2022.907190
17
HansomS. P.GolianM.GreenM. S. (2021). The Wenckebach phenomenon. Curr. Cardiol. Rev.17, 10–16. 10.2174/1573403X16666200719022142
18
HongK. L.BorgesJ.GloverB. (2020). Catheter ablation for the management of atrial fibrillation: Current technical perspectives. Open Heart7, e001207. 10.1136/openhrt-2019–001207
19
HuckerW. J.FedorovV. V.FoyilK. V.MoazamiN.EfimovI. R. (2008). Images in cardiovascular medicine. Optical mapping of the human atrioventricular junction. Circulation117, 1474–1477. 10.1161/CIRCULATIONAHA.107.733147
20
HwangH. J.NgF. S.EfimovI. R. (2014). “Mechanisms of atrioventricular nodal excitability and propagation,” in Cardiac electrophysiology: From cell to bedside. Sixth Edition (Elsevier), 275–285. 10.1016/B978-1-4557-2856-5.00028-5
21
InadaS.HancoxJ. C.ZhangH.BoyettM. R. (2009). One-dimensional mathematical model of the atrioventricular node including atrio-nodal, nodal, and nodal-His cells. Biophys. J.97, 2117–2127. 10.1016/j.bpj.2009.06.056
22
InadaS.ShibataN.IwataM.HaraguchiR.AshiharaT.IkedaT.et al (2017). Simulation of ventricular rate control during atrial fibrillation using ionic channel blockers. J. Arrhythmia33, 302–309. 10.1016/j.joa.2016.12.002
23
JohnsonN.DenesP. (2008). The ladder diagram (a 100+ year history). Am. J. Cardiol.101, 1801–1804. 10.1016/j.amjcard.2008.02.085
24
KatritsisD. G.JosephsonM. E. (2013). Classification of electrophysiological types of atrioventricular nodal re-entrant tachycardia: A reappraisal. EP Eur.15, 1231–1240. 10.1093/europace/eut100
25
KatritsisD. G.MoradyF. (2022). “Basic intervals and atrial and ventricular conduction curves,” in Clinical cardiac electrophysiology. Editors KatritsisD. G.MoradyF. (Elsevier), 54–72. 10.1016/B978-0-323-79338-4.00015-7
26
KawashimaT.SatoF. (2021). First in situ 3D visualization of the human cardiac conduction system and its transformation associated with heart contour and inclination. Sci. Rep.11, 8636. 10.1038/s41598-021-88109-7
27
KerrC. R.PrystowskyE. N.BrowningD. J.StraussH. C. (1980). Characterization of refractoriness in the sinus node of the rabbit. Circ. Res.47, 742–756. 10.1161/01.res.47.5.742
28
KléberA. G.RudyY. (2004). Basic mechanisms of cardiac impulse propagation and associated arrhythmias. Physiol. Rev.84, 431–488. 10.1152/physrev.00025.2003
29
LiJ.InadaS.SchneiderJ. E.ZhangH.DobrzynskiH.BoyettM. R. (2014). Three-dimensional computer model of the right atrium including the sinoatrial and atrioventricular nodes predicts classical nodal behaviours. PLoS ONE9, e112547. 10.1371/journal.pone.0112547
30
LianJ.MüssigD.LangV. (2006). Computer modeling of ventricular rhythm during atrial fibrillation and ventricular pacing. IEEE Trans. Biomed. Eng.53, 1512–1520. 10.1109/TBME.2006.876627
31
LinL. J.BilletteJ.MedkourD.ReidM. C.TremblayM.KhalifeK. (2001). Properties and substrate of slow pathway exposed with a compact node targeted fast pathway ablation in rabbit atrioventricular node. J. Cardiovasc. Electrophysiol.12, 479–486. 10.1046/j.1540-8167.2001.00479.x
32
MasèM.GlassL.DisertoriM.RavelliF. (2012). Nodal recovery, dual pathway physiology, and concealed conduction determine complex AV dynamics in human atrial tachyarrhythmias. Am. J. Physiol. Heart Circ. Physiol.303, H1219–H1228. 10.1152/ajpheart.00228.2012
33
MazgalevT.MowreyK.EfimovI.FahyG. J.WagonerD. V.ChengY.et al (1997). Mechanism of atrioventricular nodal facilitation in rabbit heart: Role of proximal AV node. Am. J. Physiol. Heart Circ. Physiol.273, H1658–H1668. 10.1152/ajpheart.1997.273.4.H1658
34
MedkourD.BeckerA. E.KhalifeK.BilletteJ. (1998). Anatomic and functional characteristics of a slow posterior AV nodal pathway: Role in dual-pathway physiology and reentry. Circulation98, 164–174. 10.1161/01.cir.98.2.164
35
MitsuokaT.MazgalevT.DreifusL. S.MichelsonE. L. (1987). Differential vagal effects on antegrade vs. retrograde atrioventricular conduction. Am. J. Physiol.253, H1059–H1068. 10.1152/ajpheart.1987.253.5.H1059
36
MoeG. K.PrestonJ. B.BurlingtonH. (1956). Physiologic evidence for a dual A-V transmission system. Circ. Res.4, 357–375. 10.1161/01.res.4.4.357
37
NagumoJ.ArimotoS.YoshizawaS. (1962). An active pulse transmission line simulating nerve axon. Proc. IRE50, 2061–2070. 10.1109/jrproc.1962.288235
38
NarulaO. S. (1974). Retrograde pre-excitation. Comparison of antegrade and retrograde conduction intervals in man. Circulation50, 1129–1143. 10.1161/01.cir.50.6.1129
39
NikolskiV. P.JonesS. A.LancasterM. K.BoyettM. R.EfimovI. R. (2003). Cx43 and dual-pathway electrophysiology of the atrioventricular node and atrioventricular nodal reentry. Circ. Res.92, 469–475. 10.1161/01.RES.0000059304.97120.2F
40
O’HaraG.GendreauR.BilletteJ.AmellalF.NayebpourM.TalajicM.et al (1993). Rate-dependent functional properties of retrograde atrioventricular nodal conduction in experimental animals. Am. J. Physiol.265, H1257–H1264. 10.1152/ajpheart.1993.265.4.H1257
41
Palacios-PradoN.HuetterothW.PeredaA. E. (2014). Hemichannel composition and electrical synaptic transmission: Molecular diversity and its implications for electrical rectification. Front. Cell. Neurosci.8, art. 324. 10.3389/fncel.2014.00324
42
ParsonsS. P.HuizingaJ. D. (2020). A myogenic motor pattern in mice lacking myenteric interstitial cells of Cajal explained by a second coupled oscillator network. Am. J. Physiol. Gastrointest. Liver. Physiol.318, G243. 10.1152/ajpgi.00311.2019
43
PlankG.LoeweA.NeicA.AugustinC.HuangY. L.GsellM. A. F.et al (2021). The openCARP simulation environment for cardiac electrophysiology. Comput. Methods Programs Biomed.208, 106223. 10.1016/j.cmpb
44
PlappertF.WallmanM.AbdollahpurM.PlatonovP. G.ÖstensonS.SandbergF. (2022). An atrioventricular node model incorporating autonomic tone. Front. Physiol.13, 976468. 10.3389/fphys.2022.976468
45
PodziemskiP.ŻebrowskiJ. J. (2013). Liénard-type models for the simulation of the action potential of cardiac nodal cells. Phys. D.261, 52–61. 10.1016/j.physd.2013.06.007
46
ReidM. C.BilletteJ.KhalifeK.TadrosR. (2003). Role of compact node and posterior extension in direction-dependent changes in atrioventricular nodal function in rabbit. J. Cardiovasc. Electrophysiol.14, 1342–1350. 10.1046/j.1540-8167.2003.03382.x
47
RyzhiiE.RyzhiiM. (2014). A heterogeneous coupled oscillator model for simulation of ECG signals. Comput. Methods Programs Biomed.117, 40–49. 10.1016/j.cmpb.2014.04.009
48
RyzhiiM.RyzhiiE. (2019). Optimization of dual pathway AV nodal conduction model. J. Phys. Conf. Ser.1372, 012078. 10.1088/1742-6596/1372/1/012078
49
RyzhiiM.RyzhiiE. (2022). Pacemaking function of two simplified cell models. PLoS ONE17, e0257935. 10.1371/journal.pone.0257935
50
SchölzelC.BlesiusV.ErnstG.GoesmannA.DominikA. (2021). Countering reproducibility issues in mathematical models with software engineering techniques: A case study using a one-dimensional mathematical model of the atrioventricular node. PLoS ONE16, e0254749. 10.1371/journal.pone.0254749
51
ShampineL. F.ReicheltM. W. (1997). The MATLAB ODE suite. SIAM J. Sci. Comput.18, 1–22. 10.1137/s1064827594276424
52
TadrosR.BilletteJ. (2012). Rate-dependent AV nodal function: Closely bound conduction and refractory properties. J. Cardiovasc. Electrophysiol.23, 302–308. 10.1111/j.1540-8167.2011.02180.x
53
TalajicM.PapadatosD.VillemaireC.GlassL.NattelS. (1991). A unified model of atrioventricular nodal conduction predicts dynamic changes in Wenckebach periodicity. Circ. Res.68, 1280–1293. 10.1161/01.res.68.5.1280
54
TempleI. P.InadaS.DobrzynskiH.BoyettM. R. (2013). Connexins and the atrioventricular node. Heart rhythm.10, 297–304. 10.1016/j.hrthm.2012.10.020
55
TrayanovaN. A. (2011). Whole-heart modeling: Applications to cardiac electrophysiology and electromechanics. Circ. Res.108, 113–128. 10.1161/CIRCRESAHA.110.223610
56
van der PolB.van der MarkJ. (1928). LXXII. The heartbeat considered as a relaxation oscillation, and an electrical model of the heart. Lond. Edinb. Dublin Physiological Mag. J. Sci.6, 763–775. 10.1080/14786441108564652
57
WallmanM.SandbergF. (2018). Characterisation of human AV-nodal properties using a network model. Med. Biol. Eng. Comput.56, 247–259. 10.1007/s11517-017-1684-0
58
XuB.BilletteJ.LavalléeM. (2006). Concealed conduction in nodal dual pathways: Depressed conduction, prolonged refractoriness, or reset excitability cycle?Heart rhythm.3, 212–221. 10.1016/j.hrthm.2005.11.011
59
YeW. G.YueB.AoyamaH.KimN. K.CameronJ. A.ChenH.et al (2017). Junctional delay, frequency, and direction-dependent uncoupling of human heterotypic Cx45/Cx43 gap junction channels. J. Mol. Cell. Cardiol.111, 17–26. 10.1016/j.yjmcc.2017.07.117
60
ZhangY.BharatiS.MowreyK. A.MazgalevT. N. (2003). His electrogram alternans reveal dual atrioventricular nodal pathway conduction during atrial fibrillation: The role of slow-pathway modification. Circulation107, 1059–1065. 10.1161/01.cir.0000051464.52601.f4
61
ZhangY.BharatiS.SulaymanR.MowreyK. A.TchouP. J.MazgalevT. N. (2004). Atrioventricular nodal fast pathway modification: Mechanism for lack of ventricular rate slowing in atrial fibrillation. Cardiovasc. Res.61, 45–55. 10.1016/j.cardiores.2003.10.023
62
ZhangY.MazgalevT. N. (2011). AV nodal dual pathway electrophysiology and Wenckebach periodicity. J. Cardiovasc. Electrophysiol.22, 1256–1262. 10.1111/j.1540-8167.2011.02068.x
Summary
Keywords
atrioventricular node, rabbit heart model, aliev-panfilov model, dual pathway, conduction curve, coupling asymmetry, retrograde conduction, laddergram
Citation
Ryzhii M and Ryzhii E (2023) A compact multi-functional model of the rabbit atrioventricular node with dual pathways. Front. Physiol. 14:1126648. doi: 10.3389/fphys.2023.1126648
Received
18 December 2022
Accepted
22 February 2023
Published
10 March 2023
Volume
14 - 2023
Edited by
Elena Tolkacheva, University of Minnesota Twin Cities, United States
Reviewed by
Peter Michael Van Dam, University Medical Center Utrecht, Netherlands
Sanjay Ram Kharche, Western University, Canada
Updates
Copyright
© 2023 Ryzhii and Ryzhii.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Maxim Ryzhii , m-ryzhii@u-aizu.ac.jp
This article was submitted to Cardiac Electrophysiology, a section of the journal Frontiers in Physiology
Disclaimer
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.