Modeling dynamics in diseased cardiac tissue: Impact of model choice

Abstract

Computational models of cardiac conduction can play an important role in understanding the mechanisms of and developing treatments for fibrosis-induced cardiac arrhythmias. However, most studies investigating the arrhythmogenic role of cardiac fibrosis have modelled the myocardium as a locally homogenous continuum with slowed conduction and increased conduction anisotropy in the areas of fibrosis. While these continuum methods recreate normal conduction in healthy tissues and allow for increased computational efficiency, it remains unclear how effectively they are able to model complex conduction patterns such as re-entry. In this work, we show that continuous models incompletely represent re-entrant behavior because they are unable to capture changes in tip trajectory and restitution. We find that hybrid models, which incorporate discrete heterogeneity into the continuous tissue representation, are able to capture some aspects of complex re-entry and may allow for accurate simulation with preserved computational efficiency.

INTRODUCTION

Understanding the mechanisms of cardiac arrhythmia initiation and perpetuation has proven to be a challenging problem because of the difficulties in studying the dynamics of abnormal cardiac activity in vivo and recapitulating critical features of interest in vitro. Computational modeling of cardiac electrophysiology provides a useful framework for understanding the mechanisms of and developing treatments for arrhythmia because of the ability to carefully control parameters of interest. However, in order to accurately draw conclusions about the nature of pathological conduction, models must be constructed with sufficient fidelity to reproduce the structural components underlying human diseases.

Because cardiac tissue is a discontinuous network of highly connected individual cells, and because conduction is affected by the discrete cellular structure of the tissue,³³ cardiac electrical activity has often been modeled using microstructural model systems that can incorporate variation in the cell shape, size and orientation as well as distinct gap junctions between neighboring cells. This type of detailed structural model allows for the study of the effects of alterations in the tissue structure on conduction, including alterations in source-load balance, gap junctional distribution, and fibrosis.^17,18,32 However, because these discrete models require subcellular resolution, they are extremely computationally expensive. In order to speed up computation and simulate larger regions of conduction, continuous models with spatial resolutions at the cellular scale and higher are typically used.¹⁵ These models provide an idealized representation of the tissue, with homogenized electrical properties that aggregate the effects of tissue structure, connectivity and cellular orientation.¹⁴ Continuous models have proven capable of replicating the dynamics of conduction for models of healthy tissues;^5,35 however, recent evidence suggests that these models may be less predictive under diseased conditions such as ischemia³⁵ and reduced coupling,^5,11 or at the interfaces of regional heterogeneity.¹⁶

Recently, the effects of fibrosis on arrhythmogenesis and arrhythmia dynamics have been of growing interest. Cardiac fibrosis is an injury response that involves the remodeling of gap junctions that couple neighboring cells, proliferation of fibroblasts, and deposition of excess collagen that separates and replaces the functional cardiac tissue, decoupling the syncytium that is necessary for normal cardiac conduction. Fibrosis varies in density, distribution and size, and different types of fibrosis (dense, interstitial, diffuse, etc) have varying impacts on conduction.^19,20 Marrouche et al. have found that the degree of fibrosis (estimated via delayed enhancement MRI) can be predictive of the likelihood of recurrent arrhythmia after atrial fibrillation ablation.²² As such, fibrosis may be an important target to new biological therapies for arrhythmia.

The best way to represent fibrosis in tissue models is still an open question. Spach et al. introduced interstitial fibrosis in discrete models by representing non-conductive collagenous septa as the decoupling of transverse cellular connections,³² which maintains the electrical effects of fibrosis without occupying space in the tissue domain. In continuous computational models, fibrosis is traditionally incorporated by decreasing conductivity values to reproduce experimentally observed conduction slowing and conduction anisotropy. While this approach reproduces the macroscopic conduction behavior, it may not capture the microscopic effects of fibrosis that could be crucial in arrhythmogenesis. The more recent approach by Costa et al.^9,10 of decoupling elements in a coarse mesh to reproduce the effects of interstitial fibrosis is analogous to the approach of Spach et al. in the discrete model, and has shown promise in reproducing microscale conduction. However, it remains unclear how any of these approaches perform in the setting of complex behavior such as re-entrant spiral waves.

In this study, we compare spiral wave behavior in both discrete and two equivalent continuous models of interstitial fibrosis to understand whether continuous computational models that condense structural complexity into reduced conductivity are sufficient to capture the details of re-entry. We explore the potential source of divergence between discrete and uniform continuous models and consider the use of hybrid models that incorporate distinct decoupling septa into a coarse continuous model.

METHODS

Tissue models

Two-dimensional monodomain microstructural models of anisotropic cardiac monolayers that incorporate discrete, uniformly distributed gap junctions were randomly generated using previously described methods.¹⁸ Tissues used were 2 cm × 2 cm, with a spatial discretization of 20 μm. Individual myocytes were an average of 100 μm in length and 20 μm in width, arranged in a brick-wall form and aligned along their length. An intracellular resistivity of 150 Ω·cm and a gap junction conductance of 2.5 μS were used. Collagen septa of variable length (selected from a Poisson distribution with a mean length of 800 μm) were inserted, parallel to cellular orientation, to occupy between 0% and 30% of transverse interfaces between myocytes.¹⁸ [Fig. 1(a)]. All transverse junctions between cells separated by collagenous septa were removed to simulate the decoupling effect of collagen deposition. 70% of the remaining combined plicate transverse junctions were decoupled to replicate the adult phenotype of gap junctions concentrated at longitudinal cell junctions.¹⁷ Ten distinct microstructural models with random septa placement and transverse decoupling were constructed for each level of fibrosis (0%, 10%, 20%, and 30%).

Open in a separate window

FIG. 1.

Discrete model of fibrosis. (a) Fibrosis was modeled in the discrete tissue by inserting non-conductive collagen septa (solid) parallel to and between individual cardiac cells (gray). (b) Increasing the percent of transverse cell-to-cell boundaries interrupted by collagenous septa leads to a substantial decrease in the transverse conduction velocity (solid) and an increase in the anisotropy of conduction (dashed).

An equivalent homogenous continuous model with a spatial discretization of 100 μm was produced for each discrete model by adjusting the tissue conductivities to match the conduction velocities (CVs) of transverse and longitudinal conduction measured in the corresponding discrete model. In the presence of fibrosis, the continuous model's mean longitudinal conductivities decreased by 5.7% from 0% to 30% fibrosis (2.82 mS/cm to 2.66 mS/cm), while mean transverse conductivities decreased by 77.2% (0.445 mS/cm to 0.101 mS/cm). Hybrid continuous models were created by inserting longitudinally oriented decoupling septa of mean length of 800 microns into a homogenous continuous model to disrupt between 10% and 50% of transverse coupling between nodes. The electrical conductivities between the tissue nodes were then adjusted such that the combination of decoupling septa and tissue conductivity resulted in conduction velocities equivalent to each discrete model.

RESULTS

Propagation in microstructural models and equivalent continuous models

The degree of fibrosis, characterized by the density of collagenous septa, is expected to have an effect on the conduction velocity of wavefronts. Conduction velocity was measured in the discrete microstructural model after applying 2 Hz planar stimulation. As the degree of fibrosis was increased, conduction velocity in the longitudinal direction decreased slightly from 50.39 ± 0.01 cm/s with no fibrosis to 48.65 ± 0.01 cm/s with 30% fibrosis (n = 10 for each case). Transverse conduction velocity slowed substantially with increased fibrosis, from 19.23 ± 0.004 cm/s to 6.83 ± 0.90 cm/s as fibrosis was increased from 0% to 30% [Fig. 1(b), solid], resulting in an increase in the anisotropy ratio from 2.62:1 to 7.26:1 [Fig. 1(b), dashed].

The introduction of septa should affect the trajectory of spiral waves in the discrete tissue. Spiral waves were induced via cross-field stimulation. Increasing fibrosis leads to a prolongation of the mean cycle length from 76.27 ± 0.04 ms in the non-fibrotic tissue to 90.27 ± 3.13 ms in 30% fibrosis, an increase of 18.4% [Fig. 2(a), solid]. In addition, the variability in the cycle length was also observed to increase with increasing fibrosis, with minimal variability between non-fibrotic models and substantial variability between 30% fibrosis models. Unipolar electrograms recorded from point electrodes exhibited increasing complexity and fractionation in discrete tissues with 30% fibrosis. Representative electrograms are shown in Fig. 2(c) (left column).

Open in a separate window

FIG. 2.

Continuous model fails to capture the details of re-entry. Re-entry was induced in discrete models and their equivalent continuous representations by cross-field stimulation. (a) The discrete model exhibits a substantial increase in the cycle length of re-entry (solid), up to 18.4% in 30% fibrosis that is not captured by the equivalent continuous models (dashed). Cycle lengths are shown as mean ± standard error. (b) Transmembrane potential traces show a clear difference in cycle length between the discrete model (solid) and the homogenous continuous model (dashed) at high degrees of fibrosis (lower panel) that is not seen in the absence of fibrosis (upper panel). Scale bar represents 50 ms. Note that traces were time-shifted to align the first action potential. (c) Electrograms in discrete model simulations exhibit complexity and fractionation in the presence of fibrosis, which is incompletely captured in the equivalent continuous models.

An equivalent continuous model was created for each discrete microstructural model by matching the microstructural model's transverse and longitudinal conduction velocities following 2 Hz pacing, and spiral waves were induced in each continuous model. The non-fibrotic continuous model exhibited cycle lengths similar to those of the discrete non-fibrotic tissue (CL = 74.73 ms vs 76.27 ms). However, at higher degrees of fibrosis, the continuous tissues showed only a modest increase in cycle length, with a maximal increase to 82.24 ± 2.63 ms (10% increase) in 30%-fibrosis-equivalent tissue. As a result, the equivalent continuous model has an 8.9% error in the projected cycle length at 30% fibrosis, which is apparent from transmembrane potential traces from discrete and equivalent continuous simulations [Fig. 2(b)]. In addition, while unipolar electrograms from the equivalent continuous model are qualitatively similar to those of the discrete model in the absence of fibrosis, they fail to exhibit the fractionation seen in electrograms of fibrotic microstructural models [Fig. 2(c)].

Factors affecting the cycle length

To understand what factors lead to the longer cycle length in the discrete model that are not fully captured in the continuous models, the impact of tip trajectory and restitution behaviour was explored. Analysis of tip trajectory in the control (non-fibrotic) case showed that non-fibrotic discrete and continuous models exhibited similar tip trajectories that were hypotrochoid-like in shape, with a mean span of 3.4 mm in length and 0.96 mm in width [Fig. 3(a), left]. In the presence of 30% fibrosis, the tip trajectory in the continuous model remained unchanged in shape, but was significantly flattened, with a mean span of only 0.36 mm in width [Fig. 3(a), right]. However, in the discrete model, the tip trajectory was much longer and thinner (with a mean span of 8.2 mm in length and 0.42 mm in width) [Fig. 3(a), center], suggesting an anatomical rather than functional form of re-entry.

Open in a separate window

FIG. 3.

Factors contributing to continuous model cycle length error. (a) In the absence of fibrosis, both the discrete (left panel) and continuous (not shown) models exhibit a hypotrochoid-like spiral wave tip trajectory. In the presence of fibrosis, the discrete mode (center panel) transits from functional to anatomical re-entry with a long, linear tip trajectory, while the equivalent continuous model (right panel) continues to exhibit hypotrochoid-like tip trajectory, albeit significantly flattened due to slowed transverse conduction. Note that all figures show only a part of the domain. Scale bar = 2.5 mm. (b) Restitution behavior during longitudinal pacing is nearly equivalent in the discrete (solid) and homogenous continuous (dashed) models, both in the presence and absence of fibrosis. In the transverse direction, the continuous model restitution curve closely matches the discrete model in non-fibrotic tissues, with slight deviation (up to 14.0%) at intervals below 100 ms. In contrast, in fibrotic tissues, the continuous model exhibits significantly faster conduction (up to 43.4% faster) at intervals below 180 ms.

Restitution behavior in the discrete and continuous models with 0% and 20% fibrosis was examined using a standard S1-S2 pacing protocol, in both the longitudinal and transverse directions. In all cases, because the homogenous continuous models were constructed to replicate discrete model CVs at 2 Hz pacing, the CVs are equivalent at longer diastolic intervals (of 500 ms or greater). As the diastolic interval is decreased during longitudinal pacing, the continuous model behavior closely matches that of the discrete model [Fig. 3(b), top row]. Restitution in the continuous model deviates significantly from the discrete model only at intervals under 100 ms (p < 0.05). During transverse pacing of the non-fibrotic tissue, the restitution of the continuous model deviates from the discrete model at intervals below 140 ms, with a maximum error of 14.0% at an S1-S2 interval of 90 ms [Fig. 3(b), bottom left]. A larger deviation is noted in the case of fibrotic tissue [Fig. 3(b), bottom right], where the continuous model significantly underestimates CV slowing at S1-S2 intervals of 180 ms or below, with an error of 43.4% at an interval of 90 ms. As a result, the discrete model CV restitution curve has a maximum slope that is 1.58 times greater than that of the continuous model.

In order to understand how the character of fibrosis affects continuous model estimations of cycle length, the mean length of collagenous septa in the discrete model was adjusted between 0 and 800 μm, and the cycle length was examined in discrete and equivalent continuous models. Because the total amount of fibrosis was maintained at 30% of transverse node-to-node junctions, tissues with shorter fibrosis lengths contained more individual septa. In the absence of fibrosis, the continuous model underestimates the cycle length by 2.1% (1.53 ms). At a mean fibrosis length of 200 μm, there is minimal error between the discrete and equivalent continuous models (0.14% error); however, as fibrosis length is increased further, the continuous model predicts cycle lengths that are significantly shorter than that observed in the discrete model. (Fig. ​(Fig.44).

Open in a separate window

FIG. 4.

Effect of collagen septa size. Error between the discrete model (solid) and the equivalent continuous models (dashed) is dependent on the size of collagenous septa. Small collagenous septa result in a minimal difference between the discrete and continuous models. In all cases, total septa density was fixed at 30%.

Propagation in hybrid models

In order to capture the complexity of discrete fibrosis and maintain the computational efficiency of continuous models, we considered the development of hybrid models that incorporate coarse spatial discretization and distinct decoupling septa. Hybrid models were generated to replicate the behavior of 30% fibrotic discrete tissue. The continuous model was modified by disrupting 30% of transverse node-to-node connections with non-conductive septa (mean length 800 μm), and tissue conductivities were adjusted such that the combined effect of septa and conductivity tuning matched the longitudinal and transverse CVs observed in each discrete model. The resulting models exhibited spiral wave behavior with a mean cycle length of 83.9 ms, 1.6 ms longer than that in the homogenous continuous model, but still 6.37 ms shorter than the cycle length of re-entry in the discrete model (Fig. ​(Fig.5).5). The amount of connection-disrupting collagen in the hybrid models was then altered, and each model was re-tuned to match discrete planer conduction velocities. Increasing the fibrosis density above 30% in the hybrid model resulted in the prolongation of cycle lengths, and models with 50% transverse decoupling closely matched the cycle lengths of the discrete model (mean cycle length, 90.27 ms in 30% fibrotic discrete model vs 89.88 in 50% fibrotic hybrid model). Simulation of 1.5 s of spiral wave activity required substantially less computational time in the hybrid model than in the discrete model, and there was a minimal difference in the simulation time between the continuous and hybrid models [Fig. 5(b); a mean of 4.0 h for discrete model vs 0.19 h for the continuous model and 0.21 h for the hybrid model]. Three representative unipolar electrograms from the hybrid model are shown in Fig. 5(c). Electrograms from the hybrid tissues appear to be qualitatively similar in fractionation and complexity to those obtained from discrete tissues.

Open in a separate window

FIG. 5.

Hybrid models improve continuous model accuracy. The introduction of discrete heterogeneities (in the form of decoupling septa) in the continuous model allows for accurately replicating discrete model behavior. (a) The degree of fibrosis in the hybrid model was tuned to match discrete model cycle length. For each level of fibrosis, tissue conductivities were first adjusted to match discrete planar conduction velocities. The hybrid model with 50% fibrosis was selected for further use. (b) The hybrid model maintains 18× computational efficiency of the continuous model when each is run with equal computing resources. (c) The selected hybrid model produces electrograms that qualitatively display complexity and fractionation similar to the discrete model. (d) The tip trajectory in the hybrid model is similar to that of the discrete model with a long and narrow path. Note that this panel shows only a part of the domain. Scale bar = 2.5 mm. (e) The transverse restitution behavior of the hybrid tissue (dotted) remains the same as that of the continuous model (dashed) rather than that of the discrete model (solid).

In order to understand how the hybrid model re-creates discrete model behavior, we examined hybrid model restitution and tip trajectory. The tip trajectory profile of the spiral wave in the hybrid model follows a path very similar to the discrete model. The trajectory is thin and long (span of 7.3 mm in length and 0.35 mm in width) and appears to be anatomical in nature rather than functional [Fig. 5(e)]. However, the transverse restitution profile of the hybrid model closely matches that of the homogenous continuous model rather than the discrete model [Fig. 5(d)], indicating that the hybrid model's ability to replicate discrete cycle lengths is primarily due to tip trajectory rather than dynamic variation in conduction velocity.

DISCUSSION

In this study, we examined simulated re-entry in continuous models of cardiac propagation and considered whether continuous models of cardiac propagation are able to accurately capture the details of re-entry in fibrotic myocardium. We showed that (i) continuous models fail to fully capture the increases in the cycle length of re-entry observed in microstructural discrete models of interstitial fibrosis, (ii) differences between discrete and continuous models are due to differences in tip trajectory and restitution behavior, and (iii) hybrid models that incorporate distinct edge-decoupling on coarse discretization are able to, to some degree, recapitulate discrete model behavior without compromising the computational efficiency.

Using a discrete model of cardiac tissue, we first showed that increasing the interstitial fibrosis density in discrete, microstructural models of cardiac tissue causes transverse conduction velocity to decrease at a decreasing rate, and longitudinal velocity to be minimally decreased, as previously described.¹⁸ The resulting anisotropy ratios fall within that range previously described experimentally.^31,32 When spiral waves are induced in these tissues, we observed that as the degree of fibrosis is increased, the cycle length of re-entry increases, similar to previous findings in models of diffuse fibrosis.³⁸ In addition, as the degree of fibrosis was increased, the variability in the planar conduction velocity and the cycle length increased as well, suggesting that the spatial arrangement of fibrosis grows in importance at higher fibrosis densities.

Because several groups have modeled fibrosis in large tissue or organ models as a homogenous decrease in conductivity to match experimental CVs (for example, Refs. ²¹ and ⁴²), we used the same approach in developing equivalent continuous models for each discrete tissue model. Despite the fact that the continuous models each closely match the conduction behavior of discrete models during 2 Hz linear pacing, when re-entry was induced at higher densities of fibrosis, the resulting cycle lengths were significantly shorter than those in the discrete models [Fig. 2(a)]. This finding suggests that homogenized models are inadequate to fully capture the details of cardiac re-entry in the diseased tissue. The impact of this discrepancy in the context of whole atrial or ventricular models remains unclear. While the absolute error in the cycle length is relatively modest (8.9% of the discrete cycle length), an error on this magnitude could have a significant impact on the determination of spiral wave stability or instability.⁸ In addition, previous studies have used computational models to develop tools for interpretation of electrograms, but discrepancy in the complexity of electrograms in the discrete and continuous models suggests that continuous models do not produce electrograms that capture the complexity of the underlying substrate.

In order to understand this difference between continuous and discrete cycle lengths, we examined the differences in restitution behavior in these models. Even though the discrete fibrotic model and the equivalent continuous model have been calibrated to have an equal conduction velocity in response to 2 Hz planar pacing, propagation at different rates or of a premature stimulus may draw out differences between the models. Because spiral waves in the discrete tissues do not cycle at exactly 2 Hz, differences in restitution could affect the cycle length of re-entry. A comparison of non-fibrotic discrete and continuous tissues revealed a small difference in restitution behavior at very short coupling intervals [Fig. 3(b), left column]. These differences, possibly caused by source-load mismatches at cell-to-cell gap junctions due to decreased sodium excitability, may result in the small difference in the cycle length between discrete and continuous models in the absence of fibrosis. In comparison to non-fibrotic tissues, fibrotic tissues exhibited a larger difference in transverse restitution between discrete and continuous models [Fig. 3(b), bottom right]. The additional slowing seen in the discrete model is likely due to source-load mismatches at regions where multiple nearby septa result in tissue expansions, which have been shown to lead to local conduction slowing.³⁰ The resulting difference in restitution, previously described in models of diffuse fibrosis,¹² contributes to the prolongation of the cycle length in highly fibrotic tissues that is not captured by the continuous model. Restitution behavior plays a critical role in the initiation and maintenance of arrhythmia.^3,29 As such, it is vital for computational models to accurately capture restitution beyond just the effects on cycle length. Variation in the slope of the CV restitution curve has specifically been linked to the formation of alternans, a precursor to re-entry,⁷ and in the presence of fibrosis, the steepness of the discrete model's transverse CV restitution curve is poorly recapitulated by the continuous model. Taken together, these differences represent significant deficiencies in the continuous model under conditions where the cycle length is not constant, such as during alternans and the initiation of re-entry.

We also studied the path of the spiral wave tip during re-entry in continuous and discrete models and found a substantial prolongation in the path length of re-entry in the fibrotic tissue in the discrete model, but not in the continuous model. This near doubling of the longitudinal tip path length plays a critical role in the prolongation of the spiral wave cycle length. While adjusting conductivities in a continuous model may capture the impact of gap junctional remodeling during fibrosis, the effects of collagen deposition cannot be easily homogenized.

In recognition of potential limitations of using continuous models to replicate the effects of fibrosis, Costa et al. recently described a hybrid method for inserting non-conductive septa into a coarse finite element mesh,^9,10 analogous to the insertion of septa into our discrete model.¹⁸ This method has recently been used, in combination with simulated junctional remodeling and fibroblast proliferation, in several studies of atrial fibrosis by McDowell et al.^23,25,26 Several other methods of incorporating discrete fibrosis into a coarse model have been described. Percolation-based models include randomly distributed decoupled links or fully decoupled sites to reproduce the effects of fibrosis.^1,2,39 Studies using link-percolation are functionally equivalent to the introduction of non-conductive septa of a fixed size equal to the discretization of the domain, while site-percolation methods create space-occupying fibrosis that may reproduce functional changes in tissue geometry. The space-fractional model incorporates heterogeneity more abstractly as a deviation from the standard laws of diffusion.⁶ However, the fidelity of these methods, in comparison to discrete cardiac models, remains unclear. To investigate this, we developed a hybrid model by inserting non-conductive septa into a coarse finite difference mesh. The prevalence of septa was varied and each model was tuned to match planar conduction velocities of the 30%-fibrosis discrete models. We found that 30% transverse decoupling by septa was insufficient to match the discrete cycle length, likely because the coarse discretization of the hybrid model meant that the total length of inserted septa was substantially smaller than in the discrete model. By increasing the number of septa such that they decouple 50% of transverse node-node connections, we were able to closely match the cycle lengths of re-entry observed in the discrete tissue models. While this quantitative 20% difference in septa density is not directly translatable into human fibrotic models because of the inherent limitations of our 2D study, these simulations demonstrate that hybrid models are able to reproduce the macroscale behavior of the discrete model when the nature of the modeled fibrosis is appropriately tuned.

The hybrid model appears to recreate this behavior of the discrete model primarily by capturing the prolonged tip trajectory of the discrete model [Fig. 5(d)]. The presence of discrete heterogeneities in the hybrid model allows for a transition to anatomical re-entry. However, the hybrid model fails to capture the restitution effects of fibrosis in the discrete model. This is potentially because the source-load imbalance that leads to discrete CV slowing at shorter coupling intervals is lost at the hybrid models coarser spatial discretization. Despite this observed difference in restitution behavior, hybrid models show promise in reproducing the aspects of complex behaviors of the discrete model with a substantially reduced computational load, and may prove to be a valuable tool when carefully tuned for their intended purpose.

While homogenization is insufficient in modeling the effects of interstitial fibrosis, our findings suggest that other types of fibrosis, such as diffuse fibrosis, may be adequately represented by altering conductivities of continuous models because of the relatively small size of diffuse collagenous deposits. As such, hybrid models that incorporate both distinct heterogeneities and reduced overall conductivity may be ideally suited to recreate the physiology of conduction in a complex fibrotic substrate. The conductive changes in interstitial fibrosis do not occur in isolation, but rather in concert with a host of other fibrotic changes. Proliferation of myofibroblasts in the myocardium may increase the susceptibility to arrhythmia,²⁴ and the combination of myofibroblast proliferation and gap junctional remodeling is implicated in the initiation of re-entry.²⁶ Models seeking to study the mechanisms of and treatments for arrhythmia must therefore incorporate not only the effects of collagenous deposition, but also an accurate description of the impact of microscale myocyte-myocyte and myocyte-fibroblast coupling.

CONCLUSIONS

In this study, we examined the ability of continuous cardiac models to capture complex conduction behaviour such as re-entry in the setting of interstitial fibrosis. Our results indicate that the difference in restitution behaviour and spiral tip trajectory leads to modest error in continuous model estimations of cycle length. However, we demonstrate that hybrid models that incorporate non-conductive septa with coarse discretization are able to recapitulate some aspects of complex re-entry, while providing significant computational efficiency. We suggest that it is critical to carefully consider and validate methodologies of incorporating fibrosis in computational tissue models before these models can be used to gain a mechanistic understanding of and devise treatments for clinical arrhythmia.

ACKNOWLEDGMENTS

References