Table of Contents Author Guidelines Submit a Manuscript
BioMed Research International
Volume 2018, Article ID 4798512, 10 pages
https://doi.org/10.1155/2018/4798512
Research Article

Effect of Transmural Differences in Excitation-Contraction Delay and Contraction Velocity on Left Ventricle Isovolumic Contraction: A Simulation Study

1Institute of Solid Mechanics, Mechatronics and Biomechanics, Faculty of Mechanical Engineering, University of Technology, Brno, Czech Republic
2Department of Cardiovascular Diseases, Faculty of Medicine, Masaryk University, Brno, Czech Republic
3Department of Physiology, Faculty of Medicine, Masaryk University, Brno, Czech Republic
4Institute of Thermomechanics, Czech Academy of Science, Prague, Czech Republic

Correspondence should be addressed to M. Pásek; zc.inum.dem@kesapm

Received 12 December 2017; Revised 1 March 2018; Accepted 13 March 2018; Published 10 May 2018

Academic Editor: Kazunori Uemura

Copyright © 2018 J. Vaverka et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Recent studies have shown that left ventricle (LV) exhibits considerable transmural differences in active mechanical properties induced by transmural differences in electrical activity, excitation-contraction coupling, and contractile properties of individual myocytes. It was shown that the time between electrical and mechanical activation of myocytes (electromechanical delay: EMD) decreases from subendocardium to subepicardium and, on the contrary, the myocyte shortening velocity (MSV) increases in the same direction. To investigate the physiological importance of this inhomogeneity, we developed a new finite element model of LV incorporating the observed transmural gradients in EMD and MSV. Comparative simulations with the model showed that when EMD or MSV or both were set constant across the LV wall, the LV contractility during isovolumic contraction (IVC) decreased significantly ( was reduced by 2 to 38% and IVC was prolonged by 18 to 73%). This was accompanied by an increase of transmural differences in wall stress. These results suggest that the transmural differences in EMD and MSV play an important role in physiological contractility of LV by synchronising the contraction of individual layers of ventricular wall during the systole. Reduction or enhancement of these differences may therefore impair the function of LV and contribute to heart failure.

1. Introduction

Myocardial infarction and heart failure belong to leading causes of death in the Western world. Both the clinical practice and laboratory research have shown that the mechanical performance of left ventricle (LV) is one of the most important factors that affect the pump function of the heart. Thus, understanding the processes underlying the physiological and pathophysiological function of LV is unavoidable for efficient treatment of reduced LV contractility. Computational modelling has been proposed and actively pursued as a tool for accelerating research in cardiac biomechanics. Up to now, many attempts were done to describe mechanical properties of the heart and to simulate their impact on cardiovascular hemodynamics under health and disease (for latest reviews see [13]).

One of the favourite approaches enabling simulation of the effect of electromechanical properties of cardiac muscle on development of tension in LV wall and intraventricular pressure during heart contraction is based on three-dimensional finite element (FE) models (e.g., [411]). To approach the real conditions of LV, studies published to date proposed the models with respect to (i) a real (asymmetric) geometry of human LV [58, 12]; (ii) electromechanical interaction with the right ventricle (RV) [6, 10]; (iii) transmural variation in fibre orientation [48, 1012]; (iv) sequential activation of different heart regions [10, 11]; (v) viscoelastic behaviour of the myocardial tissue [4, 13].

In 2004, Cordeiro et al. [17] showed that the time period between local electrical depolarization and onset of local shortening of myocytes (electromechanical delay: EMD) and the myocyte shortening velocity (MSV) varies transmurally in the ventricular myocardium. This reflects intrinsic differences in characteristics of membrane transport systems and excitation-contraction coupling in the myocytes from different layers of the ventricle [17, 18]. Later, Campbell et al. [19] embedded a cellular excitation-contraction coupling model with region-dependent parameters in a simple FE model of LV to study the effects of transmural electromechanical heterogeneity on local myocardial mechanics and global hemodynamics. Their study showed that a random assignment of parameters related to subendo-, midmyo-, and subepicardial cells to the LV model elements (i.e., a transmurally homogenous arrangement of the myocyte subtypes) affected significantly the transmural patterns of fibre and cross-fibre strains during early systole, but not in its later phases. Consequently, calculated parameters of the LV hemodynamic function (peak LV pressure, the maximum rate of its rise, and stroke volume) were almost identical to those in the control model in which the myocyte parameters were specified per elements according to myocardial region (transmurally heterogeneous arrangement of myocyte subtypes). However, the physiological importance of transmural differences in EMD and MSV during the systolic phases was not studied. Therefore, the aim of our work was to incorporate these wall depth dependent differences in the mechanical activity of cardiomyocytes into a newly developed LV model and to explore their role in the physiological function of LV during isovolumic contraction (IVC).

2. Methods

2.1. Geometry and FE Discretization

To simulate the ventricular contraction, a three-dimensional FE model was built using commercial FE software ANSYS®. Since the anatomical realism is not crucial for the problem addressed in this study, the model was based on a simplified geometry. The inner cavity and the outer surface of the ventricle are approximated by ellipsoids with the same axes of symmetry, truncated by a plane perpendicular to their major axis mimicking the LV basal plane. The dimensions were set according to literature (see Figure 1) to define the LV end-diastolic volume of 122 ml.

Figure 1: Geometry of the finite element model of left ventricle. The model consists of 7 layers (see details in the middle part), each of them containing differently oriented elements representing muscle fibres (bottom part).

The myocardial wall was discretized with 6048 quadratic hexahedral solid elements arranged in 7 layers. To mimic the active contraction of the myocytes, an additional reinforcing element in the form of a membrane parallel to both LV surfaces was embedded within each solid element. This membrane exhibits a unidirectional stiffness only to model the fibres aligned in a given direction. The angle of fibres against the circumferential direction (helix angle) varies linearly across the 7 layers from +60° in the inner (subendocardial) layer to −60° in the outer (subepicardial) layer [20, 21]. The thickness of the reinforcing membrane was set to occupy 80% of volume of the respective solid element which reflects the volume fraction of myocytes in ventricular myocardium [22]. The blood volume inside the cavity was modelled as an incompressible liquid using special hydrostatic fluid elements capable of modelling fluids fully enclosed by solids. The mesh density (Figure 1) was proven to be sufficient for the analysis.

2.2. Passive Properties of Myocardium

To describe the passive mechanical behaviour of myocardium, an incompressible transversely isotropic hyperelastic model was used. The preferred direction of the model representing the direction of muscle fibres was set for each solid element individually following the same rule as mentioned above for the reinforcing elements. The model was defined by a polynomial strain energy density function having the formwhere and are stress-like material parameters, is the first invariant of right Cauchy-Green deformation tensor describing the isotropic part of deformation, and is (pseudo) invariant of the right Cauchy-Green deformation tensor representing squared stretches in the fibre directions.

The material parameters were obtained as the best fit of different biaxial tests of myocardium published by Sommer et al. [23]; using our inhouse software Hyperfit (http://www.hyperfit.wz.cz) the following values were obtained: 0,347 kPa, 13,438 kPa, 48,846 kPa, 0,436 kPa, and 27,692 kPa. Since the passive response of myocardium is rate-dependent (effect of viscoelasticity), we utilized the results from equibiaxial tests performed at different testing speeds as well as from relaxation tests [23] and estimated that under strain rates typical for normally beating heart the myocardial stiffness is approximately 6 times higher compared with the normal testing speed. Therefore the above constants were multiplied by a factor of 6 in our model.

2.3. Active Fibre Contraction

As mentioned above, active contraction was incorporated into the model through reinforcing elements oriented in the local directions of myocytes in LV myocardium. To simulate heterogeneous time-dependent unidirectional contraction of myocytes, a simple approach based on prescribing time-dependent (negative) thermal strain in the fibre direction was used. A coefficient of thermal expansion was prescribed and the time course of contraction was governed by a gradual decrease in fictitious temperature (having no relation to real temperature or thermal properties of myocardium). The unidirectional stress-strain relation for activated myocytes was estimated from the observed linear relation between force and length in fully contracted isolated cardiomyocytes. This relationship was independent of the mode of contraction and initial sarcomere length [24] and conforms to the formulawhere , , and  kPa are Cauchy stress, logarithmic strain, and initial Young’s modulus, respectively. Due to limitations of ANSYS program, the above relation was incorporated into the model using a concept of multilinear elasticity with zero compression stiffness.

2.4. Boundary Conditions

In each simulation, the onset of contraction was calculated for each contractile element individually taking into account the time course of electrical activation of LV myocardium and the transmural differences in EMD. The electrical activation was assumed to start simultaneously on the whole endocardial surface and to spread towards the epicardial surface with the velocity of 47 cm/s [25, 26]. The time of electrical activation was calculated as the ratio of the shortest distance of centroid of a particular element from the cavity surface and the conduction velocity.

As exact measurements of EMD and MSV across the wall of human LV are missing we implemented these features in our model on the basis of measurements from dog hearts published by Cordeiro et al. [17]. They show a decrease of EMD in LV wall of 9 mm thickness from 47 ms (subendocardium) to 28 ms (subepicardium); thus we used a transmural EMD gradient of 2.1 ms/mm to simulate IVC under control conditions. In conjunction with the conduction velocity of 47 cm/s, this resulted in a simultaneous contraction of all contractile elements because the transmural depolarization gradient was completely counterbalanced by the EMD gradient.

The time course of contraction (controlled by temperature decrease; see Section 2.3) was defined on the basis of literature: a sequence of temperature values was prescribed to each contractile element to reproduce the experimental shortening traces of unloaded cells. As the initial part of shortening trace including IVC is typically fairly linear [2729], we assumed a linear shortening of myocytes. Consistent with the experimental traces from Cordeiro et al. [17, 30], the MSV in subepicardial layer was set two times higher than in subendocardial layer with linear decrease between them. The actual values were determined so that the intraventricular pressure increased from its end-diastolic value (1,3 kPa) to diastolic aortic pressure (10,7 kPa) in 60 ms (normal value of IVC time [31]). This normally contracting ventricle is denoted as the control model below.

The nonzero end-diastolic pressure causes the real ventricular wall to not be unloaded at the beginning of contraction. Thus the FE mesh in Figure 1 cannot represent the unloaded LV configuration and application of end-diastolic pressure on this geometry would lead to incorrect shape and volume of LV at the beginning of IVC. Therefore we employed an iterative algorithm based on the approach proposed by Bols et al. [32] to compute the unloaded (zero-pressure) geometry. As a result of the iterative process, the zero-pressure geometry was obtained with intraventricular volume of 105 ml. Upon the application of end-diastolic pressure to the hydrostatic fluid elements, the original ellipsoidal shape was restored and simulation of IVC begins. In all these simulations zero displacements were imposed on all nodes at the basal plane to model the firm attachment of LV myocardium to the stiff collagenous skeleton. As the end-diastolic pressure is balanced only by passive wall stresses with myocytes being relaxed, the contractile elements were inactivated during the initial “filling” phase until the contraction begins.

In addition to the control model, three artificial sets of boundary conditions were simulated to explain the role of heterogeneous EMD and MSV. First, we combined the layer-dependent MSV with EMD being uniform in the whole model. Thus the mechanical activation was asynchronous having the same pattern as the excitation front. Secondly, the subendocardial MSV was prescribed to all layers while the EMD time were kept the same as in the control simulation. Finally, both EMD and MSV were considered constant in the LV myocardium.

3. Results

3.1. Intraventricular Pressure and Wall Stress in Control Model

Behaviour of the control model is illustrated in Figure 2(a) showing a time course of intraventricular pressure rise from its end-diastolic value of 1.3 kPa (10 mm Hg) up to 10.7 kPa (80 mm Hg) within 60 ms. This is fully consistent with the physiological rise of intraventricular pressure measured in human LV during the preejection period (see digitalized experimental traces in Figure 2(a)). Regarding the maximum rate of intraventricular pressure rise , the model generates a value of 237.3 kPa/s (1780 mm Hg/s) that falls into the range of physiological values reported in the literature [33] and is attained in the end of IVC, just like in normal human heart [34]. These values show that the control model is capable of reliably simulating the physiological changes of intraventricular pressure during IVC in a real heart. Figure 2(b) shows a related time course of stress (tension) development in individual layers. It is apparent that during IVC the tension in the subepicardial layer is higher than in subendocardial layer. However, its transmural change is not monotonous and under control conditions is reduced in the inner parts of the LV wall (see Figure 4).

Figure 2: Time course of intraventricular pressure rise during IVC and the underlying increase of tension in LV wall under control conditions. (a) Simulated rise of intraventricular pressure and its comparison with digitalized experimental records published by Curtiss et al. [14], Manolas [15], and Kern et al. [16]. Note that the duration of IVK ( ms) and (1780 mm Hg/s) agrees with clinically measured values (61 ms, 63 ms, 61 ms, and 1790 mm Hg/s, 1730 mm Hg/s, and 2130 mm Hg/s, resp.). (b) Simulated rise of (active) stress along the direction of myocytes in the individual layers. The colours of individual lines (from blue to red) denote the layers of LV wall (from endocardium to epicardium - see Figure 1).
3.2. Effect of Transmural Differences in Electromechanical Delay and Myocyte Shortening Velocity on Intraventricular Pressure Rise during Isovolumic Contraction

The transmural differences in EMD and MSV belong to the important characteristics of LV closely related to its function. To explore the physiological importance of the transmural decrease of EMD (see Section 2.4) in the function of LV during IVC, we compared the simulated time course of intraventricular pressure rise in the control model with that in which EMD was set constant (47 ms) in all ventricular layers. The results presented in Figure 3(a) (full and dashed lines) clearly show that the transmural homogeneity of EMD leads to a prolongation of IVC from 60 to 71 ms. However, was only slightly changed (from 1780 mm Hg/s in control to 1740 mm Hg/s at constant EMD). An analysis of wall tension development during IVC under constant EMD revealed that this effect arose from a delay in the onset of contraction in the midmyocardial and subepicardial layers (Figure 3(b)) causing initial deceleration of the pressure rise.

Figure 3: Effect of transmural differences in the electromechanical delay (EMD) and in the myocyte shortening velocity (MSV) on the rate of ventricular pressure rise , duration of isovolumic contraction (IVC) and the underlying tension in LV wall. (a) Simulated increase of intraventricular pressure under control conditions and when either EMD, or MSV, or both were set constant in all ventricular layers. (b) Simulated rise of tension in individual layers when EMD was set constant across the LV wall. (c) Simulated rise of tension in individual layers when MSV was set constant across the LV wall. Colours of the individual lines in (b) and (c) denote the layers of LV wall from subendocardium (blue) to subepicardium (red).
Figure 4: Change of stress along the direction of myofibres across the LV wall (from endocardium: 0% to epicardium: 100%) in the end of IVC, under control conditions and when either EMD, or MSV, or both of them were set constant in all ventricular layers. The shaded area shows that the transmural change of stress in the inner parts of LV is reduced in control conditions.

To reveal the physiological effects of transmural differences in MSV on development of intraventricular pressure during IVC, we performed further simulation under condition of constant MSV in all ventricular layers. The results illustrated by the dotted line in Figure 3(a) show that this change would lead to a decreased during the whole IVC. This would cause a reduction of from 1780 mm Hg/s (control) to 1110 mm Hg/s and a prolongation of IVC from 60 ms (control) to 94 ms. Both changes are consistent with a decreased contractility and performance of LV. A change of tension development in individual ventricular layers underlying this effect is illustrated in Figure 3(c). As shown in the figure the rise of tension in individual layers is slowed down, which is the reason of the decreased dp/dt during IVC under this condition.

A combined effect of constant EMD and MSV on the time course of intraventricular pressure rise during IVC is shown by the dashed-dotted line in Figure 3(a). As apparent from the figure the model predicts that these conditions would cause a similar reduction of as in the case of constant MSV only (from 1780 mm Hg/s (control) to 1100 mm Hg/s). However, the prolongation of IVC duration would be more pronounced (from 60 to 104 ms). Hence the model shows that is almost exclusively affected by changes of MSV while the IVC duration is sensitive to changes of both MSV and EMD. These results suggest that while the homogenisation of both MSV and EMD across the LV wall contributes to a decreased contractility of LV, a reduction of transmural gradient of MSV plays a major role in this effect (see Table 1 for comparison).

Table 1: Values of and of IVC duration obtained from the model simulations under specified conditions.

4. Discussion

We present a new FE model of LV designed to simulate the effects of transmural differences in electromechanical characteristics of ventricular myocytes on development of tension in LV wall and intraventricular pressure during IVC. Here, a special attention is paid to exploration of the physiological role of transmural differences in EMD and MSV in these events.

4.1. Electrophysiological Contexts of Transmural Differences in Electromechanical Delay and Myocyte Shortening Velocity in Left Ventricle

The heterogeneity in electrophysiological properties of myocytes across the LV wall is commonly known. Substantial transmural distinctions in characteristics of membrane currents and action potential configurations have been demonstrated in myocytes from various species including humans [18, 3541]. Related transmural distinctions in the excitation-contraction coupling have been also documented [35, 4244]. However, transmural differences in mechanical activity of cardiomyocytes are still not conclusively explained and available experimental evidence is not fully consistent. In 2004, Cordeiro et al. [17] published their unique study examining the transmural heterogeneity of ionic currents, Ca2+ handling, and mechanical function in canine left ventricle. They identified several key distinctions in the time course of unloaded shortening among subepi-, midmyo-, and subendocardial cells. Subendocardial cells displayed a substantially longer EMD, while subepicardial cells showed the fastest MSV. Faster shortening of subepicardial cells than of those in the other myocardial layers was also observed in guinea pig hearts [45, 46] and in porcine multicellular preparations [47]. According to the coupled model of myocyte electromechanics designed by Campbell et al. [48], the distinctions in EMD in different cell types result from differences in the slope of the rising phase of intracellular Ca2+ transient. This interpretation is consistent with the observed transmural differences in a delay between onset of action potential (AP) and Ca2+ transient in human hearts [43]. As for the faster kinetics of subepicardial cells shortening, this phenomenon has been explained by the following: (i) higher activation of L-type calcium current during the early transient phase of AP repolarization induced by the presence of larger transient outward potassium current [17, 30, 48]; (ii) higher content of Ca2+ in the sarcoplasmic reticulum (SR) induced by the increased expression of SR Ca2+-ATPase [43, 44, 49]; and (iii) elevated expression of α myosin heavy chain (MHC) isoform in subepicardial myocytes which exhibit faster cross-bridge cycling kinetics [47, 50]. Besides, transmural differences in stiffness of cytoskeleton proteins [51] and in density of t-tubules as found in right atrium [52] might also play a potential role. Both shorter EMD and higher MSV in subepicardial myocytes are considered to compensate for the delayed arrival of AP into subepicardial layers, allowing synchronous contraction across the LV wall [17, 30]. This notion is highly consistent with numerous experimental and modelling studies showing that transmural fibre stress and strain are rather uniform in both diastole and systole (for review see Carruth et al. [53]). However, in vivo measurements of electrical activity and transmural myofibre mechanics in adult mongrel dogs during normal sinus rhythm or atrial pacing showed a surprisingly significant transmural gradient in the onset of myofibre shortening; the mean propagation velocity of myofibre shortening from endocardium to epicardium was 0.25 m/s [54]. In addition, the recent study on preparations from nonfailing human hearts indicated that shortening velocity of myocytes from different regions of LV does not exhibit significant differences [55]. Reasons of these contradictions are unclear. This would mean that a great part of energy produced in the heart is expended for mutual mechanical interaction between active and nonactive parts of the tissue which would decrease the pump efficiency of the heart.

4.2. Physiological Importance of Transmural Differences in Electromechanical Delay and Myocyte Shortening Velocity in Human Left Ventricle

Optimal function of the heart depends on ordered mechanical events that are orchestrated by timing of electrical events. This electromechanical coupling occurs at multiple anatomic levels: within atria, between atria and ventricles, between ventricles, and especially within LV. It was postulated that transmural uniformity of systolic fibre stress and strain likely helps to maximize the efficiency of conversion of regional contractile work into pumping function of the heart [56]. Synchronous onset of contraction in all model layers (see Section 2.2) resulting from the combination of transmural change of EMD (set on the basis of measurements done by Cordeiro et al. [17]) and average conduction velocity of 47 cm/s [25, 26] supports this view. The model simulations also showed that if one or both of these cellular parameters were set constant across the LV wall, the nonuniformity in fibre stress in inner parts of the wall increased remarkably (see the changes of Cauchy stress between 20 and 80% of wall depth in Figure 4). This was associated with a prolonged duration of IVC (by 18 to 73%) and reduced (by 2 to 38%, see Figure 3(a) and values in Table 1), indicating that a transmural homogenisation of electromechanical properties of ventricular myocytes impairs the function of the LV. Hence, our computations suggest that the transmural differences in the latency to onset of contraction and in the contraction velocity of myocytes play a physiologically important role in coordination of contraction of individual layers of ventricular myocardium, thereby ensuring a normal and most efficient LV function.

4.3. Clinical Implications

The results of our computational study indicate that a reduction of transmural gradients of EMD and MSV decreases the contractile function of LV because of discoordination of contraction in individual layers of LV. Such intraventricular dyssynchrony would lead to a decline in systolic performance, decreased mechanical efficiency, and a greater metabolic cost of LV contraction. Some changes initiated at cellular level can be identified to contribute to this pathological state. It has been documented that, in heart failure (HF), a greater reduction of in subepicardial myocytes diminishes the LV transmural gradient of early repolarization [57]. Given the relation among early repolarization, magnitude of , rate of rise of intracellular Ca2+ transient, and dynamics of cellular contraction [19, 58, 59] it is obvious that a reduction of subepicardial is an important factor underlying the attenuation of transmural gradient of EMD and MSV. It also appears that downregulation of α MHC isoform [60, 61] may contribute to the transmural reduction of both of these gradients because it exhibits higher expression in subepicardial myocytes [47]. Besides, transmural differences in disruption of excitation-contraction coupling induced by remodelling of t-tubules [62] may also play a role in this effect. On the other hand, an unphysiological increase of transmural gradients of EMD and MSV would lead to transmural desynchronization of contraction as well. It may happen for example under higher upregulation of Na-Ca exchanger current () or under higher downregulation of SR Ca2+ pump (SERCA) in subendocardial and midmyocardial cells compared to subepicardial cells as observed in HF [44, 63, 64]. Hence, all the above pathological changes typical for impaired function of the LV reduce its contractility in two ways: primarily because they decrease the contractile power of myocytes and secondarily because these changes in myocytes contractile power induce their nonsynchronous contraction across the LV wall. This complex effect of cellular impairment on LV function calls for therapeutic strategies that would not only restore the physiological concentrations of Ca2+ in cardiomyocytes (e.g., by inhibition of sodium-potassium pump or by stimulation or overexpression of SERCA) but also transmurally resynchronise their contraction. A supplementary pharmacological stimulation of subepicardial may be one of the promising ways to achieve this aim.

4.4. Limitations of the Model

The FE model used in this study succumbs to several limitations. Firstly, the model includes only the LV with its geometry approximated by two truncated ellipsoids. The omission of the mechanical interaction between the ventricles, together with the simplification of ventricular shape, affects the motion of the ventricle due to passive and active forces generated within myocardium and thus may slightly influence the calculated results. The kinematics of the ventricle is also affected by the simplified boundary condition at the base which prevents the basal movement of the ventricle. On the other hand, the method of reconstruction of the unloaded geometry is fully applicable to patient-specific geometries recorded under end-diastolic intraventricular pressure. Secondly, the activation pattern used in the model assumes simultaneous depolarization of the whole subendocardial layer and the subsequent spread of excitation towards the epicardium. This is a simplification of the real activation sequence which starts from a few early activated areas in subendocardium and proceeds around the cavity and towards the epicardium [65]. However, since the spread around the cavity is much faster than that in the transmural direction [65], the simultaneous activation of subendocardial layer represents a reasonable approximation. Thirdly, the mechanical response of passive myocardium is orthotropic rather than transversely isotropic because of the myocyte sheet structure [23, 66]; although the myocardium anisotropy is dominated by myofibre directions, a fully orthotropic hyperelastic model could provide even more accurate constitutive description of passive myocardium. Finally, the time course of contraction was derived on the basis of the assumption of linear shortening of unloaded myocytes. We believe that these limitations may change the results quantitatively but without a significant impact on the drawn conclusions.

4.5. Applicability of Model Results to Other Species

The results presented in this study hold true for humans and, to certain extent, for species with similar LV geometry, electromechanical properties of cardiomyocytes, and their similar interconnection with ventricular conduction system (Purkinje fibres). While it is clear that the changes in LV size and shape at the preserved transmural pattern of contraction would lead to substantial changes in the development of intraventricular pressure the impact of the two other factors is not clearly apparent. Nevertheless, considering the potential role of transmural gradient of in synchronisation of contraction across the human LV wall (see Section 4.2) it seems to be obvious that different mechanisms play the same role in other species that lack (e.g., guinea pig [67]) or have a slowly recovering (e.g., rabbit [68]). To unhide these mechanisms, special models accounting for these interspecies differences would be needed. Regarding the interconnection between conduction system and working myocytes, our model considered that the Purkinje fibres contact the endocardial surface at discrete points (resulting in near simultaneous activation) as has been found in human hearts [69]. However, beside endocardial also intramural Purkinje fibres have been demonstrated in ungulates including sheep [70], cow [71], pig [69], and horse [72], as well as in whales [73]. This implies that, unlike in our model, both endo- and midmyocardium are activated at the same time in these species. A tentative incorporation of this finding into our model led to transmural desynchronization of contraction starting in epicardial layers and propagating towards endocardium. The effect of this intervention on LV contractility was similar to that of constant EMD in our simulations (see Figure 3(a)) but less pronounced (IVC increased to 65 ms only and remained unchanged).

5. Conclusion

The presented results suggest that transmural differences in EMD and MSV play an important role in optimization of contractile function of LV. The control model with EMD and MSV distribution prescribed on the basis of experimental results showed the highest contraction efficiency and the increase of intraventricular pressure corresponding to physiological curves. In cases with the values of EMD, MSV, or both kept constant across the LV wall the IVC was prolonged and was reduced. This reduction of contractile power was accompanied by remarkable changes in distribution of LV wall stress. As the contractile properties of living myocytes are governed to a large extent by transmembrane ionic currents, their differential region-dependent pharmacological modulation (e.g., selective activation of in subepicardial myocytes) might be beneficial in clinical treatment of HF.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported through NETME Centre PLUS (LO1202) by financial means from the Ministry of Education, Youth, and Sports under the “National Sustainability Programme I” and through institutional support RVO: 61388998.

References

  1. R. Chabiniok, V. Y. Wang, M. Hadjicharalambous et al., “Multiphysics and multiscale modelling, data–model fusion and integration of organ physiology in the clinic: Ventricular cardiac mechanics,” Interface Focus, vol. 6, no. 2, Article ID 20150083, 2016. View at Publisher · View at Google Scholar · View at Scopus
  2. K. L. Sack, N. H. Davies, J. M. Guccione, and T. Franz, “Personalised computational cardiology: Patient-specific modelling in cardiac mechanics and biomaterial injection therapies for myocardial infarction,” Heart Failure Reviews, vol. 21, no. 6, pp. 815–826, 2016. View at Publisher · View at Google Scholar · View at Scopus
  3. A. Lopez-Perez, R. Sebastian, and J. M. Ferrero, “Three-dimensional cardiac computational modelling: methods, features and applications,” Biomedical Engineering Online, vol. 14, no. 1, article no. 35, 2015. View at Publisher · View at Google Scholar · View at Scopus
  4. J. J. Shen, F. Y. Xu, and W. A. Yang, “Finite element analysis of left ventricle during cardiac cycles in viscoelasticity,” Computers in Biology and Medicine, vol. 75, pp. 63–73, 2016. View at Publisher · View at Google Scholar · View at Scopus
  5. A. Meoli, E. Cutrì, A. Krishnamurthy et al., “A multiscale model for the study of cardiac biomechanics in single-ventricle surgeries: A clinical case,” Interface Focus, vol. 5, no. 2, Article ID 20140079, 2015. View at Publisher · View at Google Scholar · View at Scopus
  6. A. Palit, S. K. Bhudia, T. N. Arvanitis, G. A. Turley, and M. A. Williams, “Computational modelling of left-ventricular diastolic mechanics: Effect of fibre orientation and right-ventricle topology,” Journal of Biomechanics, vol. 48, no. 4, pp. 604–612, 2015. View at Publisher · View at Google Scholar · View at Scopus
  7. H. Gao, D. Carrick, C. Berry, B. E. Griffith, and X. Luo, “Dynamic finite-strain modelling of the human left ventricle in health and disease using an immersed boundary-finite element method,” IMA Journal of Applied Mathematics, vol. 79, no. 5, pp. 978–1010, 2014. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  8. M. Genet, L. C. Lee, R. Nguyen et al., “Distribution of normal human left ventricular myofiber stress at end diastole and end systole: A target for in silico design of heart failure treatments,” Journal of Applied Physiology, vol. 117, no. 2, pp. 142–152, 2014. View at Publisher · View at Google Scholar · View at Scopus
  9. H. Xia, K. Wong, and X. Zhao, “A fully coupled model for electromechanics of the heart,” Computational and Mathematical Methods in Medicine, Article ID 927279, 10 pages, 2012. View at Publisher · View at Google Scholar
  10. S. Göktepe and E. Kuhl, “Electromechanics of the heart: a unified approach to the strongly coupled excitation-contraction problem,” Computational Mechanics, vol. 45, no. 2-3, pp. 227–243, 2010. View at Publisher · View at Google Scholar · View at MathSciNet
  11. D. Nickerson, N. Smith, and P. Hunter, “New developments in a strongly coupled cardiac electromechanical model,” Europace, vol. 7, no. 2, pp. S118–S127, 2005. View at Publisher · View at Google Scholar · View at Scopus
  12. F. Dorri, P. F. Niederer, and P. P. Lunkenheimer, “A finite element model of the human left ventricular systole,” Computer Methods in Biomechanics and Biomedical Engineering, vol. 9, no. 5, pp. 319–341, 2006. View at Publisher · View at Google Scholar · View at Scopus
  13. O. Gültekin, G. Sommer, and G. A. Holzapfel, “An orthotropic viscoelastic model for the passive myocardium: continuum basis and numerical treatment,” Computer Methods in Biomechanics and Biomedical Engineering, vol. 19, no. 15, pp. 1647–1664, 2016. View at Publisher · View at Google Scholar · View at Scopus
  14. E. I. Curtiss, R. G. Matthews, and J. A. Shaver, “Mechanism of normal splitting of the second heart sound,” Circulation, vol. 51, no. 1, pp. 157–164, 1975. View at Publisher · View at Google Scholar · View at Scopus
  15. J. Manolas, “Invasive and noninvasive assessment of exercise-induced ischemic diastolic response using pressure transducers,” Current Cardiology Reviews, vol. 11, no. 1, pp. 90–99, 2015. View at Publisher · View at Google Scholar · View at Scopus
  16. M. J. Kern, P. Sorajja, and M. J. Lim, The cardiac catheterization handbook, Elsevier, Philadelphia, Pa, USA, 2015.
  17. J. M. Cordeiro, L. Greene, C. Heilmann, D. Antzelevitch, and C. Antzelevitch, “Transmural heterogeneity of calcium activity and mechanical function in the canine left ventricle,” American Journal of Physiology-Heart and Circulatory Physiology, vol. 286, no. 4, pp. H1471–H1479, 2004. View at Publisher · View at Google Scholar · View at Scopus
  18. C. Antzelevitch, S. Sicouri, S. H. Litovsky et al., “Heterogeneity within the ventricular wall. Electrophysiology and pharmacology of epicardial, endocardial, and M cells,” Circulation Research, vol. 69, no. 6, pp. 1427–1449, 1991. View at Publisher · View at Google Scholar · View at Scopus
  19. S. G. Campbell, E. Howard, J. Aguado-Sierra et al., “Effect of transmurally heterogeneous myocyte excitation-contraction coupling on canine left ventricular electromechanics,” Experimental Physiology, vol. 94, no. 5, pp. 541–552, 2009. View at Publisher · View at Google Scholar · View at Scopus
  20. C. Von Deuster, C. T. Stoeck, M. Genet, D. Atkinson, and S. Kozerke, “A reference dataset of in-vivo human left-ventricular fiber architecture in systole and diastole,” Journal of Cardiovascular Magnetic Resonance, pp. 1–4, 2015. View at Publisher · View at Google Scholar · View at Scopus
  21. D. D. Streeter Jr., H. M. Spotnitz, D. P. Patel, J. Ross Jr., and E. H. Sonnenblick, “Fiber orientation in the canine left ventricle during diastole and systole,” Circulation Research, vol. 24, no. 3, pp. 339–347, 1969. View at Publisher · View at Google Scholar · View at Scopus
  22. I. J. LeGrice, B. H. Smaill, L. Z. Chai, S. G. Edgar, J. B. Gavin, and P. J. Hunter, “Laminar structure of the heart: Ventricular myocyte arrangement and connective tissue architecture in the dog,” American Journal of Physiology-Heart and Circulatory Physiology, vol. 269, no. 2, pp. H571–H582, 1995. View at Publisher · View at Google Scholar · View at Scopus
  23. G. Sommer, A. J. Schriefl, M. Andrä et al., “Biomechanical properties and microstructure of human ventricular myocardium,” Acta Biomaterialia, vol. 24, article no. 3762, pp. 172–192, 2015. View at Publisher · View at Google Scholar · View at Scopus
  24. G. Iribe, M. Helmes, and P. Kohl, “Force-length relations in isolated intact cardiomyocytes subjected to dynamic changes in mechanical load,” American Journal of Physiology-Heart and Circulatory Physiology, vol. 292, no. 3, pp. H1487–H1497, 2007. View at Publisher · View at Google Scholar · View at Scopus
  25. P. Taggart, P. M. Sutton, T. Opthof et al., “Inhomogeneous transmural conduction during early ischaemia in patients with coronary artery disease,” Journal of Molecular and Cellular Cardiology, vol. 32, no. 4, pp. 621–630, 2000. View at Publisher · View at Google Scholar · View at Scopus
  26. A. V. Glukhov, V. V. Fedorov, P. W. Kalish et al., “Conduction remodeling in human end-stage nonischemic left ventricular cardiomyopathy,” Circulation, vol. 125, no. 15, pp. 1835–1847, 2012. View at Publisher · View at Google Scholar · View at Scopus
  27. B. Y. Poon, C. A. Ward, C. B. Cooper, W. R. Giles, A. R. Burns, and P. Kubes, “Alpha4-integrin mediates neutrophil-induced free radical injury to cardiac myocytes,” The Journal of Cell Biology, vol. 152, no. 5, pp. 857–866, 2001. View at Publisher · View at Google Scholar · View at Scopus
  28. S. U. Sys, G. W. De Keulenaer, and D. L. Brutsaert, “Reappraisal of the Multicellular Preparation for the In Vitro Physiopharmacological Evaluation of Myocardial Performance,” in Mechanisms of Work Production and Work Absorption in Muscle, G. H. Pollack and H. Sugi, Eds., vol. 453 of Advances in Experimental Medicine and Biology, pp. 441–451, Springer, Boston, MA, USA, 1998. View at Publisher · View at Google Scholar
  29. M. Boyett, J. Frampton, and M. Kirby, “The length, width and volume of isolated rat and ferret ventricular myocytes during twitch contractions and changes in osmotic strength,” Experimental Physiology, vol. 76, no. 2, pp. 259–270, 1991. View at Publisher · View at Google Scholar · View at Scopus
  30. J. M. Cordeiro, K. Calloe, R. Aschar-Sobbi et al., “Physiological roles of the transient outward current Ito in normal and diseased hearts,” Frontiers in Bioscience - Scholar, vol. 8, no. 1, pp. 143–159, 2016. View at Publisher · View at Google Scholar · View at Scopus
  31. R. Moreno, J. Zamorano, C. Almería et al., “Isovolumic contraction time by pulsed-wave Doppler tissue imaging in aortic stenosis,” European Heart Journal - Cardiovascular Imaging, vol. 4, no. 4, pp. 279–285, 2003. View at Publisher · View at Google Scholar · View at Scopus
  32. J. Bols, J. Degroote, B. Trachet, B. Verhegghe, P. Segers, and J. Vierendeels, “A computational method to assess the in vivo stresses and unloaded configuration of patient-specific blood vessels,” Journal of Computational and Applied Mathematics, vol. 246, pp. 10–17, 2013. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  33. C. F. Baicu, M. R. Zile, G. P. Aurigemma, and W. H. Gaasch, “Left ventricular systolic performance, function, and contractility in patients with diastolic heart failure,” Circulation, vol. 111, no. 18, pp. 2306–2312, 2005. View at Publisher · View at Google Scholar · View at Scopus
  34. R. L. Hamlin and C. del Rio, “dP/dt(max) - A measure of baroinometry,” Journal of Pharmacological and Toxicological Methods, vol. 66, no. 2, pp. 63–65, 2012. View at Publisher · View at Google Scholar · View at Scopus
  35. A. V. Glukhov, V. V. Fedorov, Q. Lou et al., “Transmural dispersion of repolarization in failing and nonfailing human ventricle,” Circulation Research, vol. 106, no. 5, pp. 981–991, 2010. View at Publisher · View at Google Scholar · View at Scopus
  36. T. Bányász, L. Fülöp, J. Magyar, N. Szentandrássy, A. Varró, and P. P. Nánási, “Endocardial versus epicardial differences in L-type calcium current in canine ventricular myocytes studied by action potential voltage clamp,” Cardiovascular Research, vol. 58, no. 1, pp. 66–75, 2003. View at Publisher · View at Google Scholar · View at Scopus
  37. M. A. McIntosh, S. M. Cobbe, and G. L. Smith, “Heterogeneous changes in action potential and intracellular Ca2+ in left ventricular myocyte sub-types from rabbits with heart failure,” Cardiovascular Research, vol. 45, no. 2, pp. 397–409, 2000. View at Publisher · View at Google Scholar · View at Scopus
  38. G. Li, J. Feng, L. Yue, and M. Carrier, “Transmural heterogeneity of action potentials and Ito1 in myocytes isolated from the human right ventricle,” American Journal of Physiology-Heart and Circulatory Physiology, vol. 275, no. 2, pp. H369–H377, 1998. View at Publisher · View at Google Scholar
  39. S. M. Bryant, X. Wan, S. J. Shipsey, and G. Hart, “Regional differences in the delayed rectifier current (IKr and IKs) contribute to the differences in action potential duration in basal left ventricular myocytes in guinea-pig,” Cardiovascular Research, vol. 40, no. 2, pp. 322–331, 1998. View at Publisher · View at Google Scholar · View at Scopus
  40. E. Wettwer, G. J. Amos, H. Posival, and U. Ravens, “Transient outward current in human ventricular myocytes of subepicardial and subendocardial origin,” Circulation Research, vol. 75, no. 3, pp. 473–482, 1994. View at Publisher · View at Google Scholar · View at Scopus
  41. R. B. Clark, R. A. Bouchard, E. Salinas-Stefanon, J. Sanchez-Chapula, and W. R. Giles, “Heterogeneity of action potential waveforms and potassium currents in rat ventricle,” Cardiovascular Research, vol. 27, no. 10, pp. 1795–1799, 1993. View at Publisher · View at Google Scholar · View at Scopus
  42. A. Khokhlova, N. Balakina-Vikulova, L. Katsnelson, G. Iribe, and O. Solovyova, “Transmural cellular heterogeneity in myocardial electromechanics,” The Journal of Physiological Sciences. View at Publisher · View at Google Scholar · View at Scopus
  43. Q. Lou, V. V. Fedorov, A. V. Glukhov, N. Moazami, V. G. Fast, and I. R. Efimov, “Transmural heterogeneity and remodeling of ventricular excitation-contraction coupling in human heart failure,” Circulation, vol. 123, no. 17, pp. 1881–1890, 2011. View at Publisher · View at Google Scholar · View at Scopus
  44. K. R. Laurita, R. Katra, B. Wible, X. Wan, and M. H. Koo, “Transmural heterogeneity of calcium handling in canine,” Circulation Research, vol. 92, no. 6, pp. 668–675, 2003. View at Publisher · View at Google Scholar · View at Scopus
  45. X. Wan, S. M. Bryant, and G. Hart, “A topographical study of mechanical and electrical properties of single myocytes isolated from normal guinea-pig ventricular muscle,” Journal of Anatomy, vol. 202, no. 6, pp. 525–536, 2003. View at Publisher · View at Google Scholar · View at Scopus
  46. S. M. Bryant, S. J. Shipsey, and G. Hart, “Regional differences in electrical and mechanical properties of myocytes from guinea-pig hearts with mild left ventricular hypertrophy,” Cardiovascular Research, vol. 35, no. 2, pp. 315–323, 1997. View at Publisher · View at Google Scholar · View at Scopus
  47. J. E. Stelzer, H. S. Norman, P. P. Chen, J. R. Patel, and R. L. Moss, “Transmural variation in myosin heavy chain isoform expression modulates the timing of myocardial force generation in porcine left ventricle,” The Journal of Physiology, vol. 586, no. 21, pp. 5203–5214, 2008. View at Publisher · View at Google Scholar · View at Scopus
  48. S. G. Campbell, S. N. Flaim, C. H. Leem, and A. D. McCulloch, “Mechanisms of transmurally varying myocyte electromechanics in an integrated computational model,” Philosophical Transactions of the Royal Society A: Mathematical, Physical & Engineering Sciences, vol. 366, no. 1879, pp. 3361–3380, 2008. View at Publisher · View at Google Scholar · View at Scopus
  49. T. Anderson, J. Wulfkuhle, E. Petricoin III, and R. L. Winslow, “High resolution mapping of the cardiac transmural proteome using reverse phase protein microarrays,” Molecular & Cellular Proteomics, vol. 10, no. 7, 2011. View at Publisher · View at Google Scholar · View at Scopus
  50. Y. Ait Mou, J.-Y. Le Guennec, E. Mosca, P. P. De Tombe, and O. Cazorla, “Differential contribution of cardiac sarcomeric proteins in the myofibrillar force response to stretch,” Pflügers Archiv - European Journal of Physiology, vol. 457, no. 1, pp. 25–36, 2008. View at Publisher · View at Google Scholar · View at Scopus
  51. O. Cazorla, A. Freiburg, M. Helmes et al., “Differential expression of cardiac titin isoforms and modulation of cellular stiffness,” Circulation Research, vol. 86, no. 1, pp. 59–67, 2000. View at Publisher · View at Google Scholar · View at Scopus
  52. M. Frisk, J. T. Koivumäki, P. A. Norseng, M. M. Maleckar, O. M. Sejersted, and W. E. Louch, “Variable t-tubule organization and Ca2+ homeostasis across the atria,” American Journal of Physiology-Heart and Circulatory Physiology, vol. 307, no. 4, pp. H609–H620, 2014. View at Publisher · View at Google Scholar · View at Scopus
  53. E. D. Carruth, A. D. McCulloch, and J. H. Omens, “Transmural gradients of myocardial structure and mechanics: Implications for fiber stress and strain in pressure overload,” Progress in Biophysics and Molecular Biology, vol. 122, no. 3, pp. 215–226, 2016. View at Publisher · View at Google Scholar · View at Scopus
  54. H. Ashikaga, B. A. Coppola, B. Hopenfeld, E. S. Leifer, E. R. McVeigh, and J. H. Omens, “Transmural Dispersion of Myofiber Mechanics. Implications for Electrical Heterogeneity In Vivo,” Journal of the American College of Cardiology, vol. 49, no. 8, pp. 909–916, 2007. View at Publisher · View at Google Scholar · View at Scopus
  55. P. Haynes, K. E. Nava, B. A. Lawson et al., “Transmural heterogeneity of cellular level power output is reduced in human heart failure,” Journal of Molecular and Cellular Cardiology, vol. 72, pp. 1–8, 2014. View at Publisher · View at Google Scholar · View at Scopus
  56. M. Vendelin, P. H. M. Bovendeerd, J. Engelbrecht, and T. Arts, “Optimizing ventricular fibers: Uniform strain or stress, but not ATP consumption, leads to high efficiency,” American Journal of Physiology-Heart and Circulatory Physiology, vol. 283, no. 3, pp. H1072–H1081, 2002. View at Publisher · View at Google Scholar · View at Scopus
  57. Q. He, Y. Feng, and Y. Wang, “Transient outward potassium channel: a heart failure mediator,” Heart Failure Reviews, vol. 20, no. 3, pp. 349–362, 2015. View at Publisher · View at Google Scholar · View at Scopus
  58. R. Sah, R. J. Ramirez, and P. H. Backx, “Modulation of Ca2+ release in cardiac myocytes by changes in repolarization rate: Role of phase-1 action potential repolarization in excitation-contraction coupling,” Circulation Research, vol. 90, no. 2, pp. 165–173, 2002. View at Publisher · View at Google Scholar · View at Scopus
  59. R. Sah, R. J. Ramirez, G. Y. Oudit et al., “Regulation of cardiac excitation-contraction coupling by action potential repolarization: Role of the transient outward potassium current (Ito),” The Journal of Physiology, vol. 546, no. 1, pp. 5–18, 2003. View at Publisher · View at Google Scholar · View at Scopus
  60. P. J. Reiser, M. A. Portman, X.-H. Ning, and C. S. Moravec, “Human cardiac myosin heavy chain isoforms in fetal and failing adult atria and ventricles,” American Journal of Physiology-Heart and Circulatory Physiology, vol. 280, no. 4, pp. H1814–H1820, 2001. View at Publisher · View at Google Scholar · View at Scopus
  61. S. Miyata, W. Minobe, M. R. Bristow, and L. A. Leinwand, “Myosin heavy chain isoform expression in the failing and nonfailing human heart,” Circulation Research, vol. 86, no. 4, pp. 386–390, 2000. View at Publisher · View at Google Scholar · View at Scopus
  62. A. Guo, C. Zhang, S. Wei, B. Chen, and L.-S. Song, “Emerging mechanisms of T-tubule remodelling in heart failure,” Cardiovascular Research, vol. 98, no. 2, pp. 204–215, 2013. View at Publisher · View at Google Scholar · View at Scopus
  63. W. Xiong, Y. Tian, D. DiSilvestre, and G. F. Tomaselli, “Transmural heterogeneity of Na+-Ca2+ exchange: Evidence for differential expression in normal and failing hearts,” Circulation Research, vol. 97, no. 3, pp. 207–209, 2005. View at Publisher · View at Google Scholar · View at Scopus
  64. J. Prestle, S. Dieterich, M. Preuss, U. Bieligk, and G. Hasenfuss, “Heterogeneous transmural gene expression of calcium-handling proteins and natriuretic peptides in the failing human heart,” Cardiovascular Research, vol. 43, no. 2, pp. 323–331, 1999. View at Publisher · View at Google Scholar · View at Scopus
  65. D. Durrer, R. T. van Dam, G. E. Freud, M. J. Janse, F. L. Meijler, and R. C. Arzbaecher, “Total excitation of the isolated human heart,” Circulation, vol. 41, no. 6, pp. 899–912, 1970. View at Publisher · View at Google Scholar · View at Scopus
  66. S. Dokos, B. H. Smaill, A. A. Young, and I. J. LeGrice, “Shear properties of passive ventricular myocardium,” American Journal of Physiology-Heart and Circulatory Physiology, vol. 283, no. 6, pp. H2650–H2659, 2002. View at Publisher · View at Google Scholar
  67. D. M. Nassal, X. Wan, H. Liu, and I. Deschênes, “Myocardial KChIP2 expression in Guinea pig resolves an expanded electrophysiologic role,” PLoS ONE, vol. 11, no. 1, Article ID e0146561, 2016. View at Publisher · View at Google Scholar · View at Scopus
  68. H. Cheng, M. B. Cannell, and J. C. Hancox, “Differential responses of rabbit ventricular and atrial transient outward current (Ito) to the Ito modulator NS5806,” Physiological Reports, vol. 5, no. 5, Article ID e13172, 2017. View at Publisher · View at Google Scholar · View at Scopus
  69. D. Sedmera and R. G. Gourdie, “Why do we have Purkinje fibers deep in our heart?” Physiological Research, vol. 63, pp. S9–S18, 2014. View at Google Scholar
  70. S. Ryu, S. Yamamoto, C. R. Andersen, K. Nakazawa, F. Miyake, and T. N. James, “Intramural Purkinje cell network of sheep ventricles as the terminal pathway of conduction system,” Anatomical Record, vol. 292, no. 1, pp. 12–22, 2009. View at Publisher · View at Google Scholar · View at Scopus
  71. P. W. Oosthoek, S. Viragh, W. H. Lamers, and A. F. M. Moorman, “Immunohistochemical delineation of the conduction system: II: The atrioventricular node and Purkinje fibers,” Circulation Research, vol. 73, no. 3, pp. 482–491, 1993. View at Publisher · View at Google Scholar · View at Scopus
  72. R. L. Hamlin and C. R. Smith, “Categorization of Common Domestic Mammals Based upon Their Ventricular Activation Process,” Annals of the New York Academy of Sciences, vol. 127, no. 1, pp. 195–203, 1965. View at Publisher · View at Google Scholar · View at Scopus
  73. N. Ono, T. Yamaguchi, H. Ishikawa et al., “Morphological varieties of the purkinje fiber network in mammalian hearts, as revealed by light and electron microscopy,” Archives of Histology and Cytology, vol. 72, no. 3, pp. 139–149, 2009. View at Publisher · View at Google Scholar · View at Scopus