Effects of cardiac mechanical heterogeneity on the electrical function of the heart are difficult to assess experimentally, yet they pose a serious (patho-)physiological challenge. Here, we present an in silico study of the effects of mechanical heterogeneity on action potential duration (APD) in mechanically interacting muscle regions and consequent effects on the dispersion of repolarization, a well-established determinant of cardiac arrhythmogenesis. Using a novel mathematical description of ventricular electromechanical activity (virtual muscle), we first assessed how differences in intrinsic contractile properties affect the electrical behavior of cardiac muscle representations. In spite of identical electrophysiological model descriptions in virtual muscle samples, faster muscle models show shorter APD than their slower counterparts. This is a consequence of Ca 2+-mediated feedback from mechanical to electrical activity in the individual muscle models. This mechano-electric feedback (MEF) is, of course, significantly more complex in native cardiac tissue, as the heterogeneous muscle elements interact both mechanically and electrically. Cardiac mechanical heterogeneity, in its most reduced form, can be represented by a duplex consisting of two mechanically interacting muscle segments. Our in silico model of heterogeneous myocardium therefore consists of two individual virtual muscles that are mechanically interconnected in-series to form a virtual heterogeneous duplex. During isometric contraction of the duplex (i.e. at constant external length), internal mechanical interactions affect Ca 2+ handling and APD of muscle elements, resulting in an increased dispersion of repolarization beyond the intrinsic APD differences. Duplex electromechanical activity is strongly affected by the activation sequence of its elements. Late activation of the faster (subepicardial type) duplex element, postponed by time-lags that correspond to normal transmural activation delays, optimizes duplex contractility and smoothes out intrinsic APD differences, thereby reducing dispersion in repolarization. This smoothing effect is not observed upon delayed activation of the slower (subendocardial type) duplex element. In both settings, changes in repolarization timing follow a nonlinear dependence of APD on activation delay. Furthermore, asynchronous activation of identical elements in a homogeneous duplex causes an impairment of contractile function and increases dispersion of repolarization. This suggests that the normal electrical activation sequence in the heart requires matching mechanical and electrical heterogeneity for optimal cardiac performance. On the subcellular level, our results suggest that mechanical modulation of Ca 2+ handling is a key mechanism of MEF in heterogeneous myocardium, which contributes to the matching of local mechanical and/or electrical activity to global hemodynamic demand.

The heart is structurally and functionally a highly non-homogenous organ, yet its main function as a pump can only be achieved by the co-ordinated contraction of millions of ventricular cells. This apparent contradiction gives rise to the hypothesis that 'well-organised' inhomogeneity may be a pre-requisite for normal cardiac function. Here, we present a set of novel experimental and theoretical tools for the study of this concept. Heterogeneity, in its most condensed form, can be simulated using two individually controlled, mechanically interacting elements (duplex). We have developed and characterised three different types of duplexes: (i) biological duplex, consisting of two individually perfused biological samples (like thin papillary muscles or a trabeculae), (ii) virtual duplex, made-up of two interacting mathematical models of cardiac muscle, and (iii) hybrid duplex, containing a biological sample that interacts in real-time with a virtual muscle. In all three duplex types, in-series or in-parallel mechanical interaction of elements can be studied during externally isotonic, externally isometric, and auxotonic modes of contraction and relaxation. Duplex models, therefore, mimic (patho-)physiological mechano-electric interactions in heterogeneous myocardium at the multicellular level, and in an environment that allows one to control mechanical, electrical and pharmacological parameters. Results obtained using the duplex method show that: (i) contractile elements in heterogeneous myocardium are not 'independent' generators of tension/shortening, as their ino- and lusitropic characteristics change dynamically during mechanical interaction-potentially matching microscopic contractility to macroscopic demand, (ii) mechanical heterogeneity contributes differently to action potential duration (APD) changes, depending on whether mechanical coupling of elements is in-parallel or in-series, which may play a role in mechanical tuning of distant tissue regions, (iii) electro-mechanical activity of mechanically interacting contractile elements is affected by their activation sequence, which may optimise myocardial performance by smoothing intrinsic differences in APD. In conclusion, we present a novel set of tools for the experimental and theoretical investigation of cardiac mechano-electric interactions in healthy and/or diseased heterogeneous myocardium, which allows for the testing of previously inaccessible concepts.

A mathematical model of the cardiomyocyte electromechanical function is used to study contribution of mechanical factors to rhythm disturbances in the case of the cardiomyocyte calcium overload. Particular attention is paid to the overload caused by diminished activity of the sodium-potassium pump. It is shown in the framework of the model, where mechano-calcium feedback is accounted for that myocardium mechanics may significantly enhance arrhythmogenicity of the calcium overload. Specifically, a role of cross-bridge attachment/detachment processes, a role of mechanical conditions of myocardium contractions (length, load), and a role of myocardium viscosity in the case of simulated calcium overload have been revealed. Underlying mechanisms are analyzed. Several approaches are designed in the model and compared to each other for recovery of the valid myocardium electrical and mechanical performance in the case of the partially suppressed sodium-potassium pump.

A mathematical model for the regulation of mechanical activity in cardiac muscle has been developed based on a three-element rheological model of this muscle. The contractile element has been modeled taking into account the results of extensive mechanical tests that involved the recording of length-force and force-velocity relations and muscle responses to short-time deformations during various phases of the contraction-relaxation cycle. The best agreement between the experimental and the mathematical modeling results was obtained when a postulate stating two types of cooperativity to regulate the calcium binding by troponin was introduced into the model. Cooperativity of the first type is due to the dependence of the affinity of troponin C for Ca2+ on the concentration of myosin crossbridges in the vicinity of a given troponin C. Cooperativity of the second type assumes an increase in the affinity of a given troponin C for Ca2+ when the latter is bound by molecules neighboring troponin.

BackgroundOne of the main factors affecting propagation of electrical waves and contraction in ventricles of the heart is anisotropy of cardiac tissue. Anisotropy is determined by orientation of myocardial fibres. Determining fibre orientation field and shape of the heart is important for anatomically accurate modelling of electrical and mechanical function of the heart. The aim of this paper is to introduce a theoretical rule-based model for anatomy and fibre orientation of the left ventricle (LV) of the heart and to compare it with experimental data. We suggest explicit analytical formulae that allow us to obtain the left ventricle form and its fibre direction field. The ventricle band concept of cardiac architecture given by Torrent-Guasp is chosen as the model postulate.MethodsIn our approach, anisotropy of the heart is derived from some general principles. The LV is considered as a set of identical spiral surfaces, each of which can be produced from the other by rotation around one vertical axis. Each spiral surface is filled with non-intersecting curves which represent myocardial fibres.For model verification, we use experimental data on fibre orientation in human and canine hearts.ResultsLV shape and anisotropy are represented by explicit analytical expressions in a curvilinear 3-D coordinate system. The derived fibre orientation field shows good qualitative agreement with experimental data. The model reveals the most thorough quantitative simulation of fibre angles at the LV middle zone.ConclusionsOur analysis shows that the band concept can generate realistic anisotropy of the LV. Our model shows good qualitative agreement between the simulated fibre orientation field and the experimental data on LV anisotropy, and the model can be used for various numerical simulations to study the effects of anisotropy on cardiac excitation and mechanical function.