Olga Solovyova

Mechanical interaction of heterogeneous cardiac muscle segments in silico: Effects on Ca2+ handling and action potential

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 Ca2+-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 Ca2+ 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 Ca2+ 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.

Activation sequence as a key factor in spatio-temporal optimization of myocardial function

Using one-dimensional models of myocardial tissue, implemented as chains of virtual ventricular muscle segments that are kinematically connected in series, we studied the role of the excitation sequence in spatio-temporal organization of cardiac function. Each model element was represented by a well-verified mathematical model of cardiac electro-mechanical activity. We found that homogeneous chains, consisting of identical elements, respond to non-simultaneous stimulation by generation of complex spatio-temporal heterogeneities in element deformation. These are accompanied by the establishment of marked gradients in local electro-mechanical properties of the elements (heterogeneity in action potential duration, Ca2+ transient characteristics and sarcoplasmic reticulum Ca2+ loading). In heterogeneous chains, composed of elements simulating fast and slow contracting cardiomyocytes from different transmural layers, we found that only activation sequences where stimulation of the slower elements preceded that of faster ones gave rise to optimization of the systems electro-mechanical function, which was confirmed experimentally. Based on the results obtained, we hypothesize that the sequence of activation of cardiomyocytes in different ventricular layers is one of the key factors of spatio-temporal organization of myocardium. Moreover, activation sequence and regional differences in intrinsic electro-mechanical properties of cardiac muscle must be matched in order to optimize myocardial function.

Mechanical inhomogeneity of myocardium studied in parallel and serial cardiac muscle duplexes: experiments and models

Abstract We investigate, both experimentally and theoretically, contribution of the myocardium mechanical inhomogeneity to the contractile function. We developed three approaches, named as Muscle Duplex Methods, to study the specific effects and mechanisms of interaction in the simplest myocardial system consisting of two muscular units connected either in parallel or in series. Our experimental approach is designed to study the interaction between two isolated mechanically inhomogeneous cardiac muscles. The virtual duplex approach is based on a mathematical model of the myocardium contraction. The hybrid duplex approach has been designed to support, in real time, interaction between a natural muscle and its virtual counterpart. Using these approaches we showed the existence of a fine alignment between mechanical characteristics of interacting inhomogeneous myocardial elements. Contractile properties of the elements together with particular sequences and time delays in their stimulation specifically determine this alignment. We term as “tuning effects” all the phenomena concerning the interaction between inhomogeneous systems elements. Within the framework of the mathematical model we showed that the key mechanism underlying tuning effects is a feedback between mechanical conditions and cooperative Ca2+ binding by troponin C. Thanks to the model analysis, we also hypothesize that mechanical inhomogeneity of myocardium is apt to produce its electrical inhomogeneity.

Mathematical model of the anatomy and fibre orientation field of the left ventricle of the heart

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.

Slow force response and auto-regulation of contractility in heterogeneous myocardium.

Classically, the slow force response (SFR) of myocardium refers to slowly developing changes in cardiac muscle contractility induced by external mechanical stimuli, e.g. sustained stretch. We present evidence for an intra-myocardial SFR (SFR(IM)), caused by the internal mechanical interactions of muscle segments in heterogeneous myocardium. Here we study isometric contractions of a pair of end-to-end connected functionally heterogeneous cardiac muscles (an in-series muscle duplex). Duplex elements can be either biological muscles (BM), virtual muscles (VM), or a hybrid combination of BM and VM. The VM implements an Ekaterinburg-Oxford mathematical model accounting for the ionic and myofilament mechanisms of excitation-contraction coupling in cardiomyocytes. SFR(IM) is expressed in gradual changes in the overall duplex force and in the individual contractility of each muscle, induced by cyclic auxotonic deformations of coupled muscles. The muscle that undergoes predominant cyclic shortening shows force enhancement upon return to its isometric state in isolation, whereas average cyclic lengthening may decrease the individual muscle contractility. The mechanical responses are accompanied with slow and opposite changes in the shape and duration of both the action potential and Ca²⁺ transient in the cardiomyocytes of interacting muscles. Using the mathematical model we found that the contractility changes in interacting muscles follow the alterations in the sarcoplasmic reticulum loading in cardiomyocytes which result from the length-dependent Ca²⁺ activation of myofilaments and intracellular mechano-electrical feedback. The SFR(IM) phenomena unravel an important mechanism of cardiac functional auto-regulation applicable to the heart in norm and pathology, especially to hearts with severe electrical and/or mechanical dyssynchrony.

Electrical Wave Propagation in an Anisotropic Model of the Left Ventricle Based on Analytical Description of Cardiac Architecture

We develop a numerical approach based on our recent analytical model of fiber structure in the left ventricle of the human heart. A special curvilinear coordinate system is proposed to analytically include realistic ventricular shape and myofiber directions. With this anatomical model, electrophysiological simulations can be performed on a rectangular coordinate grid. We apply our method to study the effect of fiber rotation and electrical anisotropy of cardiac tissue (i.e., the ratio of the conductivity coefficients along and across the myocardial fibers) on wave propagation using the ten Tusscher–Panfilov (2006) ionic model for human ventricular cells. We show that fiber rotation increases the speed of cardiac activation and attenuates the effects of anisotropy. Our results show that the fiber rotation in the heart is an important factor underlying cardiac excitation. We also study scroll wave dynamics in our model and show the drift of a scroll wave filament whose velocity depends non-monotonically on the fiber rotation angle; the period of scroll wave rotation decreases with an increase of the fiber rotation angle; an increase in anisotropy may cause the breakup of a scroll wave, similar to the mother rotor mechanism of ventricular fibrillation.

The cardiac muscle duplex as a method to study myocardial heterogeneity

This paper reviews the development and application of paired muscle preparations, called duplex, for the investigation of mechanisms and consequences of intra-myocardial electro-mechanical heterogeneity. We illustrate the utility of the underlying combined experimental and computational approach for conceptual development and integration of basic science insight with clinically relevant settings, using previously published and new data. Directions for further study are identified.

Contribution of mechanical factors to arrhythmogenesis in calcium overloaded cardiomyocytes: Model predictions and experiments

It is well-known that Ca²⁺ overload in cardiomyocytes may underlie arrhythmias. However, the possible contribution of mechanical factors to rhythm disturbances in Ca²⁺ overloaded myocytes has not been sufficiently investigated. We used a mathematical model of the electrical and mechanical activity of cardiomyocytes to reveal an essential role of the mechanisms of cardiac mechano-electric feedback in arrhythmogenesis in Ca²⁺ overloaded myocardium. In the model, the following mechanical factors increased Ca²⁺ overload in contracting cardiomyocytes and promoted rhythm disturbances: i) a decrease in the mechanical load for afterloaded contractions; and ii) a decrease in the initial length of sarcomeres for isometric twitches. In exact accordance with the model predictions, in experiments on papillary muscles from the right ventricle of guinea pigs with Ca²⁺ overloaded cardiomyocytes (using 0.5-1 μM of ouabain), we found that emergence of rhythm disturbances and extrasystoles depends on the mechanical conditions of muscle contraction.

Role of myocardial viscoelasticity in disturbances of electrical and mechanical activity in calcium overloaded cardiomyocytes: mathematical modeling.

Cardiomyocyte Ca(2+) overload is closely linked to cardiac arrhythmias. We have earlier shown in a mathematical model that myocardium mechanical activity may contribute to rhythm disturbances induced by Ca(2+) overload in cardiomyocytes with reduced Na(+)-K(+) pump work (Sulman et al., 2008). The same model is used here to address possible contribution of the passive mechanical properties of cardiac muscle (i.e. myocardial viscous and elastic properties) to the arrhythmogenesis. In a series of contractions at regular pacing rate of 75 beats/min a model with higher viscosity demonstrated essentially earlier appearance of extrasystoles due to a faster cardiomyocyte Ca(2+) loading up to a level triggering spontaneous Ca(2+) releases from the sarcoplasmic reticulum. The model predicts that myocardial elasticity also may affect arrhythmogenesis in cardiomyocytes overloaded with Ca(2+). Contribution of the mechanical properties of the myocardial tissue to the arrhythmia has been analyzed for wide ranges of both viscosity and elasticity coefficients. The results suggest that myocardial viscoelastic properties may be a factor affecting Ca(2+) handling in cardiomyocytes and contributing to cardiac mechano-electric feedback in arrhythmogenesis.

Mathematical modelling of mechano-electric feedback in cardiomyocytes

We earlier developed the mathematical model of electrical and mechanical activity in myocardium, which takes into account both direct coupling and feedback between excitation and contraction. In this paper, in the framework of the model we found conditions under which both the abrupt shortening and stretch of cardiac preparation can cause extra action potentials and hence anomalous deformations can be arrhythmia sources. In the framework of the model, we establish possible mechanisms underlying the Anrep phenomenon that reflects a relationship between myocardium contraction and vascular resistance in the intact heart.