Myocardial Fibrosis in a 3D Model: Effect of Texture on Wave Propagation

: Non-linear electrical waves propagate through the heart and control cardiac contraction. Abnormal wave propagation causes various forms of the heart disease and can be lethal. One of the main causes of abnormality is a condition of cardiac ﬁbrosis, which, from mathematical point of view, is the presence of multiple non-conducting obstacles for wave propagation. The ﬁbrosis can have different texture which varies from diffuse (e.g., small randomly distributed obstacles), patchy (e.g., elongated interstitional stria), and focal (e.g., post-infarct scars) forms. Recently, Nezlobinsky et al. (2020) used 2D biophysical models to quantify the effects of elongation of obstacles (ﬁbrosis texture) and showed that longitudinal and transversal propagation differently depends on the obstacle length resulting in anisotropy for wave propagation. In this paper, we extend these studies to 3D tissue models. We show that 3D consideration brings essential new effects; for the same obstacle length in 3D systems, anisotropy is about two times smaller compared to 2D, however, wave propagation is more stable with percolation threshold of about 60% (compared to 35% in 2D). The percolation threshold increases with the obstacle length for the longitudinal propagation, while it decreases for the transversal propagation. Further, in 3D, the dependency of velocity on the obstacle length for the transversal propagation disappears.


Introduction
Non-linear waves of electrical excitation propagate through the heart and initiate cardiac contraction. Myocardial tissue consists of inter-connected excitable cells-cardiomyocytes forming so-called myocardial fibers ensuring prevailing conducting pathways in the tissue. In normal myocardium, the electrical wave propagation is about 2-3 times faster in the longitudinal direction of the myocardial fibers than in the transversal directions. From a mathematical point of view, the electrical waves belong to a large class of waves in the Reaction-Diffusion systems which are widely studied in applied mathematics [1]. The structural properties of the myocardium are accounted in the models as electro-diffusion anisotropy of the medium with higher diffusion coefficient along the myofiber direction. In many forms of heart disease, normal wave propagation is disturbed by an excessive growth of the fraction of connective tissue within myocardium including inexcitable cells-fibroblasts and myofibroblasts. This pathological condition is called cardiac fibrosis.
Cardiac fibrosis is a common attribute of myocardial diseases of different etiology such as ischaemic cardiomyopathy and myocardial infarction, dilated cardiomyopathy, aortic stenosis, hypertrophic cardiomyopathy, and myocarditis [2]. Fibrosis is shown to be a predictor of adverse outcomes, including heart failure, ventricular arrhythmias, sudden cardiac death and all-cause Effects of more complex structures of fibrosis formed of interstitial stria surrounding myofibers either randomly or being organized in patchy areas that may replace a large fraction of myocardial tissue have not been sufficiently studied neither in experiment nor in simulations. Role of the texture and type of fibrosis in the rhythm disturbances and the risk of their complications is not quantitatively evaluated. One of the first works demonstrating the importance of the texture and the characteristic dimensions of unexcitable obstacles in the tissue for the onset of spiral waves was the work of A. Pertsov's group (see [28]). In our recent study, [29], we made a first attempt to study the effects of fibrotic texture representing some features of interstitial fibrosis with primarily linear accumulations of collagen separating bundles of myocytes with little alteration in the alignment of myocytes or bundles [30]. One of main features of this pattern is that inexcitable inclusions to cardiac tissue here have some directional dependency laying between myofibers. In idealised representation they can be seen as thin elongated obstacles directed along cardiac fibers. In [29], we considered case of parallel fibers, thus texture of linear inexcitable elements which all have the same orientation and we studied propagation in such 2D models of myocardium for wave along this direction (longitudinal) and perpendicular to it (transversal) propagation. Main findings of this study revealed an opposite change in the conduction velocity (CV) with the obstacle length depending on the direction of the wavefront propagation. The CV increases with obstacle length if the excitation propagates along the obstacle strips, while the CV decreases under transverse wavefront propagation. This causes increasing tissue anisotropy with obstacle elongation [28]. This effect was explained by zig-zag propagation in the tissue with more short excitation path along the obstacles than across. The 2D model we used has natural limitations not allowing the wave to propagate in the third spacial direction transverse to the longitudinal coordinate of the obstacles. In this article, we extend our study of structured fibrosis to 3D myocardial tissue models and analyze if the third spatial dimension and the depth of the 3D tissue volume affect the characteristics of the excitation wave and the arrhythmogenic power of the fibrosis.
In particular, the aim of this paper is to study in detail effects of randomly generated textures of fibrosis with various lengths and percentages of obstacles on wave propagation in cardiac tissue slab of various depth.

Electrophysiology Model
To describe propagation of the excitation wave in the myocardial tissue we used a 3-dimensional monodomain formulation: where V is transmembrane potential, D -electro-diffusion matrix for anisotropic tissue, in this study we considered D constant, I ion -sum of all transmembrane ionic currents, described with biophysically detailed cardiac action potential model TP06 [31]. Initial conditions were set as the rest potential V = V rest for the cardiac tissue. Boundary conditions were formulated as the no flux through the boundaries: where n-is the normal to the boundary. The problem was solved in a 3D cuboid simulating myocardial slab with inclusions of structured fibrosis. Fibrosis elements of the tissue were simulated as non-conducting inexcitable obstacles and considered as the boundaries (no flux) for the myocardial elements.

Numerical Methods
To solve the problem (1) and (2) we used a finite-difference method with 18-point stencil discretization scheme as described in [32] with 0.25 mm for the spatial step and 0.02 ms for the time step.

Fibrosis Pattern Generation
We considered 3D finite-difference mesh as matrix X × Y × Z (where X, Y, Z are the numbers of elements in x, y and z direction) with each element to be either excitable myocardial tissue or inexcitable fibrosis. In general we used matrices 200 × 200 × 10 elements (or 50 mm × 50 mm × 2.5 mm) in size, but we also considered matrices with 200 × 200 × Z elements, where Z might be 1, 2, 3 or 5.
Similarly to Nezlobinsky et al. [29], we uniformly distributed fibrosis elements and varied length of elements (in x-axis direction) and density of elements distribution. To distribute fibrosis elements through the mesh we implemented two approaches. Fibrosis pattern was generated in the same manner as described in [29], namely, we subdivided each X × Y × Z matrix in x or y direction into Y · Z or X · Z rows with size 1 × 200 × 1 or 200 × 1 × 1 respectively. Then we subdivided each row into blocks of length n and assign it as fibrotic with a probability p. This approach allowed us to set fibrosis percentage precisely close to desired value.
Examples of wave propagation in a myocardial slab with fibrosis of different texture are shown in Figure 3. Upper panels A-B show an example of elongated fibrosis texture with linear 3 × 1 × 1 elements of 3 x-node length occupying 60% of the tissue volume. Panel A shows the wave propagation in the longitudinal direction of the fibrotic obstacles, panel B shows the wave in the transversal direction. Bottom panel C shows an example of the diffuse fibrosis texture of the same density. Details on the wave dynamics are described in the Results section.

Shortest Path Calculation
The shortest possible path of excitation propagation between the opposite sides of the myocardial cuboid was calculated to demonstrate the zig-zag propagation through the fibrotic texture with elongated obstacles. We used iterative multi-dimensional binary dilation routine. Figure 1 shows and 2D example of our algorithm. Firstly, we mark an initial side of the mesh (all cells in left lateral surface except fibrosis) as "activated" and applied 3D binary dilatation operation until at least one cell on the opposite side of mesh becomes "activated" (see Figure 1A). Then, starting from this early approached cell on the opposite side (if there were several points we randomly picked one single point), we repeated this procedure in backward direction. The number of dilation iterations was taken as the length of shortest possible path for excitation propagation (see Figure 1B). Very first activated opposite cells for both steps were then used to calculate the trace of the shortest path. To trace the shortest possible path we build an undirected graph connecting neighbouring "active" cells in the fibrosis matrix and then breadth-first search algorithm with cells obtained on previous step assigned to be source and target nodes was used . An example of the shortest possible path trace is shown in Figure 2.

Results
We evaluated properties of the excitation wave propagating either in the longitudinal direction of the fibrosis stria (along the long face of the strips in the direction of horizontal axis x) and across them (transversal, along axis y). We varied several parameters: (1) the depth of the myocardial slab from monolayer, which is equivalent 2D tissue to a 3D multilayer slab (2, 3, 5 and 10 layers); (2) the length of fibrosis stria from random 0.25 mm cubic elements simulating diffuse fibrosis (number of nodes n = 1 on axis x) to 1.0 × 0.25 × 0.25 mm 3 strips (n = 4) elongated along axis x; (3) the density of the fibrosis from 10 to 60%.

Wavefront Velocity Depends on Fibrosis Fraction and Propagation Direction
Representative simulations of the wave propagation in the longitudinal and transverse directions in the 3D myocardial volume with structural fibrosis of 60% density formed of linear 3 × 1 × 1 elements of 3 x-node length are shown in Figure 3. We compare propagation of the wave initiated by the stimulation on the left lateral surface and spreading along the obstacle elongation (see panel A) or on the upper lateral surface and travelling across the obstacles (panel B). We see that the wavefront has a complex shape due to interaction with obstacles. We also see that the wave propagates along the obstacles about two times faster than across them. For example, 50 ms after stimulation the wave propagated at a distance of approximately 25 mm for the longitudinal direction. However, for the transversal propagation the wavefront covered about 16 mm of the tissue. Panel C shows wave propagation in the model with diffuse fibrosis of the same density initiated on the upper surface of the volume. Here, average dynamics of the wave does not depend on the direction as the fibrotic elements are cubic. The conduction velocity at diffuse fibrosis is much slower than that at the longitudinal wavefront direction for the elongated obstacles, and is near to that for the transversal wavefront direction (CV = 0.29 mm/ms in panel C versus 0.48 mm/ms in panel A, and 0.24 mm/ms in panel B). The red color shows activated cells, the gray color indicates non-activated cells and transparent cells within the gray zone indicates fibrosis. Fibrosis density is 60%. Snapshots of excitation propagation are shown at 50 (left panels) and 100 ms (right panels) after excitation onset. The figure illustrates faster excitation propagation for the wavefront co-directed along the elongation of the fibrosis elements. In this case (panel A, longitudinal wavefront), about two times larger part of the myocardial volume is activated at the certain time moments as compared to the wave propagation across the fibrosis elongation (panel B, transversal wavefront). Panel C shows wave propagation in the model with 60% diffuse fibrosis. Here, excitation wave is initiated at the upper surface, and is independent of the wavefront direction. The conduction velocity at diffuse fibrosis is much slower than that at the longitudinal wavefront direction for the elongated obstacles, and is near to that for the transversal wavefront direction (CV = 0.29 mm/ms in panel C versus 0.48 mm/ms in panel A, and 0.24 mm/ms in panel B). Figure 4 shows the dependence of the conduction velocity (CV) on the percentage of the fibrosis for textures with 0.25 to 1 mm long elements (n = 1, 2, 3, 4 x-nodes). Upper panels A, B show the results we obtained in a 3D monolayer of the tissue (0.25 mm depth, 1 z-layer). As this case is equivalent to 2D model, we use it as a reference for comparison of 2D and 3D models, in particular with the results reported earlier [29]. We see that the velocity significantly decreases with an increase in the fibrosis percentage. For the diffuse fibrosis (n = 1) the velocity does not depend on the direction of activation (see red lines in Figure 4) and is used as the reference to compare with results for structural fibrosis. In the monolayer model, elongation of obstacles has a pronounced opposite effect on the longitudinal and transversal CV with increasing longitudinal velocity with obstacle length but decreasing transversal velocity.
Middle C, D panels in Figure 4 show the results obtained in the multilayer tissue of the 2.5 mm depth (10 z-layers). It is seen that qualitative effects of fibrosis density and obstacle length on the CV are similar to that of the monolayer model. At the same time, the tissue depth allowing excitation to propagate around obstacles in the additional transversal direction (z-axis). It differently affect the longitudinal and transversal propagation. We see a significant effect of the obstacle length on the longitudinal CV similar to monolayer case, while the transversal CV in the ten-layer tissue model is almost independent of the obstacle length in contrast to the monolayer model.  Therefore, the texture of fibrosis produces anisotropic propagation with increasing longitudinal/transversal anisotropy ratio with the density of fibrosis and obstacle length. This is illustrated in Figure 5. The anisotropy ratio in multilayer model reaches 2:1 with increase in the fibrosis density, the dependency is steeper at longer obstacles and more steep in the monolayer vs. multilayer models. The dependency is non-linear at higher than 25% fibrosis percentage.

Wave Propagation Stopping Depends on the Tissue Depth, Fibrosis Percentage and Wavefront Direction
The CV dependencies in Figure 4 are spanned for different ranges of fibrosis percentage for monolayer models (panels A,B) and multilayer models (panels C,D), indicating different percolation threshold for 2D and 3D cases. The fibrosis percentage in the monolayer models is limited by 35%, after which the wave cannot propagate anymore, while in the multilayer models the wave propagates for the fibrosis density up to 55%. Close to percolation threshold we observed a steep increase in the variability of CV for both models (see increasing STD whiskers on the curves in Figure 4). For example, in the monolayer models with 35% fibrosis and three-node long obstacles, coefficient of variability in CV is of 21.8% for the longitudinal direction and 11.5% for the transverse direction. In the multilayer models with 60% fibrosis and three-node long obstacles the coefficients of variability in CV are 1.8% and 3.6% respectively. Figure 6 characterizes the percolation threshold statistically: it shows the fraction of models with failed wave propagation depending on the fibrosis percentage in the mono-(upper panels) and multilayer (lower panels) model. Here, we used the binary dilation method to analyze if the wave is able to approach the opposite side of the tissue volume (see Methods section for detail) and counted a number of models with propagation failing. The dependencies are approximated by logistic Hill curves.
It is seen that the curves are right-shifted for longitudinal wavefront (left panels) compared to transversal wavefront (right panels) and in the multilayer models (lower panels) compared to monolayer models (upper panels). This means that the probability of excitation failure is higher for for the transversal versus longitudinal wavefront propagation and in the monolayer against multilayer tissue. Moreover, the curves shift in the opposite directions for the longitudinal and transversal propagation when the obstacle length is increasing.
We see that for the longitudinal propagation the probability of excitation failure decreases while for the transversal propagation the probability of excitation failure increases with increase in the obstacle length.
To clarify these results further we show the dependencies of the percolation threshold FP 95 (i.e., fibrosis percentage where propagation failure occurs for 95% of models) on the obstacle length for the longitudinal and transversal propagation in the mono-and multilayer models (Figure 7). It is seen that FP 95 is much higher in the multilayer (upper lines) compared to the monolayer (lower lines) models, and FP 95 is higher for the longitudinal wavefront direction (solid lines) compared to the transversal (dotted lines). For longitudinal wave direction the minimal FP 95 is observed for the diffuse fibrosis and FP 95 is about 43% for the monolayer models and FP 95 and is about 65% for the multilayer models. For transversal propagation the minimal FP 95 is lower than the minimal FP 95 for the longitudinal propagation. It is observed for the fibrosis of 4-node long and FP 95 is about 35% for the monolayer models and FP 95 and is about 60% for the multilayer models. Note, that in Figures 4 and  5 we considered the latter fibrosis percentages as upper limits for the mono-and multilayer models.

Effects of Fibrosis on the Conduction Velocity Depends on the Tissue Depth
In Figure 8 we compare effects of obstacle length on the CV in the ten-layer models (red lines) with that in the mono-layer models at 35% fibrosis. We see that the length of the obstacle substantially affects longitudinal propagation and the dependency in 3D is similar to that in 2D. However, for the transversal propagation the velocity is almost independent on the obstacle length in 3D. To quantify further effects of the tissue depth on the degree of CV change with the length of fibrosis obstacles, we demonstrate the ratio of the average CV for the models with 1 mm (n = 4) obstacle length to the model of diffuse fibrosis of the same percentage ( Figure 9) at different thickness of the tissue (number of layers is 1, 2, 3, 5, and 10). We see big change in velocity between 1 and 2 layers. However starting from 3-layer thick tissue, the curves become independent of the number of layers. So, for both wavefront directions the thickness when 3D effects become dominating is about 0.5 mm.

Shortest Path for Wave Propagation Depends on the Obstacle 0 and the Tissue Depth
In paper [29], we explained an opposite effect of increasing obstacle length on the longitudinal and transverse propagation velocity by a zig-zig propagation [8]. It was characterized by estimating the shortest path length between opposite boundaries of the medium. Here, we use the similar approach and calculate the shortest path from one to the opposite face of the 3D myocardial slab in the longitudinal and transverse direction (see Methods section in detail). As shown in Figure 2 the path is more straight and short in the longitudinal direction and it is more complex and long in the transverse direction. The distribution of the shortest path length for the 3D mono-(upper panels) and multilayer models (lower panels) for the same fibrosis density of 30% for obstacle length n = 4 in comparison with diffuse fibrosis with n = 1 are shown in Figure 10. It is seen that the average path is much shorter for the longitudinal wavefront direction than that for the diffuse fibrosis and the latter is shorter than the path for the transversal wavefront direction. Interestingly, for a particular texture shown in Figure 2 the ratio of conduction velocities in the longitudinal and transversal directions is 2.47, which is reasonably close to the ratio of average path lengths (2.23). Note that the average path length almost does not change with the tissue depth for the longitudinal wavefront direction. In contrast, the average path for the transverse direction essentially shortens in the multilayer tissue (116.71 for transverse fibrosis of n = 4 node length and 30% density. see Figure 10) coming closer to the value for the average value for the diffuse fibrosis models.

Discussion
In this paper, we extended our previous study on the effects of fibrosis texture on wave propagation in cardiac tissue by transition from 2D to 3D models with varying depth of the myocardial slab. First, we checked if the effects we observed in the 2D models are reproduced in the 3D formulations of a monolayer model with 1 node depth of the tissue. Then, we increased the number of layers in the slab to evaluate effects of the third spatial dimension on the properties of the wave in the fibrosis textures of various densities and obstacle elongation. The simple model of straight elongated obstacle stria in the myocardial tissue were considered as a simplified representation of interstitial fibrosis [8] localized between myocardial fibers replacing part of cardiomyocytes. To evaluate effects of obstacle elongation in comparison with more rounded small elements of the diffuse fibrosis, we varied the obstacle length considering the texture with fibrosis strips of equal length. Then, we varied fibrosis percentage in a rather wide range to approach such density where excitation wave stopped due to obstacle texture blocking the propagation at a certain distance from the initially stimulated area.
We show that the elongation of obstacles has a substantial effect on the waves. It results in anisotropic wave propagation in which the velocity of the longitudinal propagation increases with increasing elongation, while the velocity of the transversal direction decreases (Figure 4). The extent of the effect increases with the fibrosis percentage and with the the obstacle length.
We found that qualitative effects of the fibrosis texture are similar in the 3D and 2D models of the tissue, which validates main predictions based on the 2D models and allows one to translate them into the more adequate 3D settings. At the same time we showed that the effects of fibrosis in the 3D tissue are less expressed due to additional spacial dimension transversal to the obstacle elongation that allows the wave to find more ways to bypass the fibrotic obstacles. The effects of the third tissue dimension are more pronounced for the transversal wavefront direction increasing the conduction velocity in the 3D multilayer models compared to monolayer models and making almost negligible difference between the CV at various obstacle lengths at the same fibrosis percentage (Figure 4).
We showed that in 3D tissue the dependencies of the tissue anisotropy on the fibrosis density are much less steep ( Figure 5), so the higher fibrosis percentage is needed to approach the same extent of anisotropy as produced by the 2D models. For instance, an anisotropy of 2:1 was approached at about 30% fibrosis in the monolayer models while the same anisotropy was demonstrated in the models of 60% fibrosis in ten-layer deep models.
We showed that threshold fibrosis density for conduction block is also about twice higher in the 3D multilayer models as compared to the monolayer models (see Figure 7). Of note, the threshold fibrosis density increases with obstacle elongation for the longitudinal wavefront direction, while decreases for the transversal direction that makes much less percentage of fibrosis to be more dangerous.
The geometric approach we used here for determining the shortest path of the wave is the main approach which was used in most of the papers on fibrosis [12,13,19,34]. In real myocardial tissue, the path of the wave is governed by the dynamic electrotonic load (source-sink relationship), which can affect the propagation up to producing local blocks in spite of the geometrical connectivity [3]. In our previous work [29] we specially analyzed if the geometry consideration only is sufficient to estimate the path and to discriminate the cases of different elongation of the fibrosis elements. We performed additional simulations initiating a propagation wave by point stimulation at one boundary and calculating its shortest path to the opposite boundary. We found that the source-sink relationship slightly increases the path of zig-zag propagation as compared to percolation path for both the longitudinal and transversal wavefront direction, and an increase in the transversal path was about 15%, while for the longitudinal only about 3%. We also found that the ratio of the average path lengths and the ratio of the conduction velocities (anisotropy) in the longitudinal and transversal directions were almost the same for a particular fibrosis texture. Thus, we showed that the percolation path provided a good estimate for the real path of the wave and the difference in the geometry-based paths at fibrosis of different textures reflects adequately the real difference rather underestimating it quantitatively.
These simulations suggest that the severity of the fibrosis effects on the myocardial tissue is essentially dependent on the local density of the fibrosis and its orientation relative to the myocardial fibers that can vary depending on the tissue depth and orientation of the myofibers within the heart chambers.
Our research is a computation prediction on some 3D effects of fibrosis texture on wave propagation. It would be interesting to compare it with available experimental data. There is a lot of experimental data which qualitatively confirms the results of our study. In many papers and in various preparations it was shown that fibrosis decreases the conduction velocity of the wave propagation (see reviews [3,4]). Further in [35] it was shown that increased interstitial fibrosis affects transverse conduction much more than the longitudinal conduction. Further, one of our results is a complex shape of the wave in the presence fibrosis ( Figure 3) is also qualitatively confirmed in [36], where it was shown that increased fibrosis has also been associated with heterogeneioty of conduction. It was shown in patients with dilated cardiomyopathy that the amount of myocardial fibrosis correlated with the severity of abnormal propagation [30].
However, currently quantitative comparison is impossible. This is because it is challenging to create cardiac tissues with controlled textures of fibrosis in which length interstitial stria gradually changes. However, we think that using modern technologies of cardiac cell co-cultures [37] and/or optogenetics [38] it will be possible to mimic multiple elongated obstacles with controlled length and study wave propagation there. There are also recent developments which allow construction of such 3D preparations [39]. It would be interesting to use these technologies to study wave propagation in such experimental systems and compare it with the predictions of our modelling.
Our approach has several limitations. The major limitation is that we studied models of idealised fibrosis formed by linear inexcitable obstacles of the same length. The geometry and texture of fibrosis in real tissue is much more complex on many aspects. Fibrosis can be formed not only by fibroblasts, but also by myofibroblasts, which are much larger in size. Substantial amounts of collagen are present in fibrotic cardiac tissue and have different space scale and organisation than fibroblasts. Further, direction, size and shape of fibrotic inclusions are not so regular as we assumed in our paper. It would be interesting to perform studies on realistic textures of cardiac fibrosis in the human heart similar to those from [40] obtained for 3D tissue [41] using computational approaches presented in study [42,43]. Further, a novel approach of modelling which takes into account shapes of cardiac cells can be an interesting tool to better represent fibrosis textures [26]. In this paper we studied effects of fibrosis texture on the wave propagation. It would be also interesting to study effects of fibrosis texture on the onset and properties of cardiac arrhythmias which either occur due to re-entry mechanism of abnormal propagation related to the structural abnormality of the tissue [44] or due to ectopic activity related to the abnormal depolarization-repolarization of cardiac cells [45].
We considered 3D isotropic tissue without taking into account pre-excisting anysotrophy of myocardium due to myofiber structure with different electro-diffusion along and across the fibers. This is also the task for our future studies on realistic textures of cardiac fibrosis.

Numerical Implications for Fibrosis Modelling
Currently two main approaches are used to represent fibrotic remodeling: (1) modification of the diffusion tensor in the regions where fibrosis is present [20] and (2) introduction of real obstacles or internal boundary conditions (this study, [19]). Although, in our view approach (2) can potentially better represent texture of the tissue and better reproduce mechanisms of initiation of arrhythmias, its numerical implementation is much more challenging. In our paper we used approach (2), however, our approach can also be applied for simulations performed using approach (1) in combination with patient specific data in the following way. For some patients with cardiac myopathy a biopsy of endocardial tissue is often taken. The bioptic tissue can be studied using histological methods and a patient specific texture of fibrosis for this given patient is available. This texture can be input to computer and fibrotic tissue can be modelled as obstacles, as we did in our paper. Then we can compute the velocities of propagation for wave in three orthogonal directions and based on that calculate diffusion coefficients as well as a rule for their modification. We illustrate it with a specific example. We generated a fibrotic texture in a small 3D patch of cardiac tissue (Supplementary Materials Figure S1). The preexisting anisotropy was 1:0.5:0.5 with respect to velocity. From studies of wave propagation in 3 orthogonal direction we found that for given fibrosis texture (50% fibrosis produced by obstacles of 0.75 mm × 0.25 mm) decrease in velocities can be reproduced by decrease in diffusion coefficient by 1.75, 2.08, 2.08 for the longitudinal and 2 transversal directions. Thus if one wants to introduce effect of such fibrosis to the large scale model, the diffusion coefficients in fibrotic regions should be reduced by such factors, e.g., instead of 1:0.25 :0.25 one should use 0.57:0.12:0.12.
In conclusion, in 3D idealised models for human myocardial tissue, we showed that inexcitable fibrotic elements produce anisotropy of cardiac tissue that depends on the fibrosis percentage, obstacle elongation and tissue depth. The anisotropy results from zig-zag propagation depending on the elongation of fibrosis elements. We showed that a threshold fibrosis density to stop excitation propagation depends on the obstacle length and tissue depth.
Supplementary Materials: The following are available online at http://www.mdpi.com/2227-7390/8/8/1352/s1, Figure S1: Electrical wave propagation in a model with structural fibrosis and a reduced model simulating propagation timing in the fibrosis tissue. Table S1: Electro-diffusion coefficients (D), anisotropy (longitudinal to transverse ratio) and conduction velocity (CV) in the models with structural fibrosis and corresponding reduced model.