Mechanisms of mechanically induced spontaneous arrhythmias in acute regional ischemia

A discrete electromechanical model for human cardiac tissue: effects of stretch-activated currents and stretch conditions on restitution properties and spiral wave dynamics

We introduce an electromechanical model for human cardiac tissue which couples a biophysical model of cardiac excitation (Tusscher, Noble, Noble, Panfilov, 2006) and tension development (adjusted Niederer, Hunter, Smith, 2006 model) with a discrete elastic mass-lattice model. The equations for the excitation processes are solved with a finite difference approach, and the equations of the mass-lattice model are solved using Verlet integration. This allows the coupled problem to be solved with high numerical resolution. Passive mechanical properties of the mass-lattice model are described by a generalized Hooke's law for finite deformations (Seth material). Active mechanical contraction is initiated by changes of the intracellular calcium concentration, which is a variable of the electrical model. Mechanical deformation feeds back on the electrophysiology via stretch-activated ion channels whose conductivity is controlled by the local stretch of the medium. We apply the model to study how stretch-activated currents affect the action potential shape, restitution properties, and dynamics of spiral waves, under constant stretch, and dynamic stretch caused by active mechanical contraction. We find that stretch conditions substantially affect these properties via stretch-activated currents. In constantly stretched medium, we observe a substantial decrease in conduction velocity, and an increase of action potential duration; whereas, with dynamic stretch, action potential duration is increased only slightly, and the conduction velocity restitution curve becomes biphasic. Moreover, in constantly stretched medium, we find an increase of the core size and period of a spiral wave, but no change in rotation dynamics; in contrast, in the dynamically stretching medium, we observe spiral drift. Our results may be important to understand how altered stretch conditions affect the heart's functioning.

Mechanobiology in cardiac physiology and diseases

Mechanosensitivity is essential for heart function just as for all other cells and organs in the body, and it is involved in both normal physiology and diseases processes of the cardiovascular system. In this review, we have outlined the relationship between mechanosensitivity and heart physiology, including the Frank-Starling law of the heart and mechanoelectric feedback. We then focused on molecules involved in mechanotransduction, particularly mechanosensitive ion channels. We have also discussed the involvement of mechanosensitivity in heart diseases, such as arrhythmias, hypertrophy and ischaemic heart disease. Finally, mechanobiology in cardiogenesis is described with regard to regenerative medicine.

[Molecular and genetic background of sudden cardiac death]

Despite recent findings on the functional, structural and genetic background of sudden cardiac death, the incidence is still relatively high in the entire population. A thorough knowledge on susceptibility, as well as pathophysiology behind the development of malignant arrhythmias will help us to identify individuals at risk and prevent sudden cardiac death. This article presents a review of the current literature on the role of altered intracellular Ca2+ handling, acute myocardial ischaemia, cardiac autonomic innervation, renin-angiotensin-aldosterone system, monogenic and complex heritability in the pathogenesis of sudden cardiac death.

Computational cardiology: the heart of the matter

This paper reviews the newest developments in computational cardiology. It focuses on the contribution of cardiac modeling to the development of new therapies as well as the advancement of existing ones for cardiac arrhythmias and pump dysfunction. Reviewed are cardiac modeling efforts aimed at advancing and optimizing existent therapies for cardiac disease (defibrillation, ablation of ventricular tachycardia, and cardiac resynchronization therapy) and at suggesting novel treatments, including novel molecular targets, as well as efforts to use cardiac models in stratification of patients likely to benefit from a given therapy, and the use of models in diagnostic procedures.

Distension of the ischemic region predicts increased ventricular fibrillation inducibility following coronary occlusion in swine

INTRODUCTION AND OBJECTIVES

Distension of the ischemic region has been related to an increased incidence of spontaneous ventricular arrhythmias following coronary occlusion. This study analyzed whether regional ischemic distension predicts increased ventricular fibrillation inducibility after coronary occlusion in swine.

METHODS

In 18 anesthetized, open-chest pigs, the left anterior descending coronary artery was ligated for 60 min. Myocardial segment length in the ischemic region was monitored by means of ultrasonic crystals. Programmed stimulation was applied at baseline and then continuously between 10 and 60 min after coronary occlusion.

RESULTS

Coronary occlusion induced a rapid increase in end-diastolic length in the ischemic region, which reached 109.4% (0.9%) of baseline values 10 min after occlusion (P<.001). On average, 6.6 (0.5) stimulation protocols were completed and 5.4 (0.6) ventricular fibrillation episodes induced between 10 and 60 min of coronary occlusion. Neither baseline serum potassium levels nor the size of the ischemic region were significantly related to ventricular fibrillation inducibility. In contrast, the increase in end-diastolic length 10 min after coronary occlusion was associated directly (r=0.67; P=.002) with the number of induced ventricular fibrillation episodes and inversely (r=-0.55; P=.018) with the number of extrastimuli needed for ventricular fibrillation induction.

CONCLUSIONS

Regional ischemic expansion predicts increased ventricular fibrillation inducibility following coronary occlusion. These results highlight the potential influence of mechanical factors, acting not only on the triggers but also on the substrate, in the genesis of malignant ventricular arrhythmias during acute ischemia.

Mechano-regulation of the beating heart at the cellular level--mechanosensitive channels in normal and diseased heart

The heart as a contractile hollow organ finely tunes mechanical parameters such as stroke volume, stroke pressure and cardiac output according to filling volumes, filling pressures via intrinsic and neuronal routes. At the cellular level, cardiomyocytes in beating hearts are exposed to large mechanical stress during successive heart beats. Although the mechanisms of excitation-contraction coupling are well established in mammalian heart cells, the putative contribution of mechanosensitive channels to Ca²⁺ homeostasis, Ca²⁺ signaling and force generation has been primarily investigated in relation to heart disease states. For instance, transient receptor potential channels (TRPs) are up-regulated in animal models of congestive heart failure or hypertension models and seem to play a vital role in pathological Ca²⁺ overload to cardiomyocytes, thus aggravating the pathology of disease at the cellular level. Apart from that, the contribution of mechanosensitive channels (MsC) in the normal beating heart to the downstream force activation cascade has not been addressed. We present an overview of the current literature and concepts of mechanosensitive channel involvement in failing hearts and cardiomyopathies and novel data showing a likely contribution of Ca²⁺ influx via mechanosensitive channels in beating normal cardiomyocytes during systolic shortening.

Computational approaches to understand cardiac electrophysiology and arrhythmias

Cardiac rhythms arise from electrical activity generated by precisely timed opening and closing of ion channels in individual cardiac myocytes. These impulses spread throughout the cardiac muscle to manifest as electrical waves in the whole heart. Regularity of electrical waves is critically important since they signal the heart muscle to contract, driving the primary function of the heart to act as a pump and deliver blood to the brain and vital organs. When electrical activity goes awry during a cardiac arrhythmia, the pump does not function, the brain does not receive oxygenated blood, and death ensues. For more than 50 years, mathematically based models of cardiac electrical activity have been used to improve understanding of basic mechanisms of normal and abnormal cardiac electrical function. Computer-based modeling approaches to understand cardiac activity are uniquely helpful because they allow for distillation of complex emergent behaviors into the key contributing components underlying them. Here we review the latest advances and novel concepts in the field as they relate to understanding the complex interplay between electrical, mechanical, structural, and genetic mechanisms during arrhythmia development at the level of ion channels, cells, and tissues. We also discuss the latest computational approaches to guiding arrhythmia therapy.

Toward an integrative computational model of the Guinea pig cardiac myocyte

The local control theory of excitation-contraction (EC) coupling asserts that regulation of calcium (Ca²⁺) release occurs at the nanodomain level, where openings of single L-type Ca²⁺ channels (LCCs) trigger openings of small clusters of ryanodine receptors (RyRs) co-localized within the dyad. A consequence of local control is that the whole-cell Ca²⁺ transient is a smooth continuous function of influx of Ca²⁺ through LCCs. While this so-called graded release property has been known for some time, its functional importance to the integrated behavior of the cardiac ventricular myocyte has not been fully appreciated. We previously formulated a biophysically based model, in which LCCs and RyRs interact via a coarse-grained representation of the dyadic space. The model captures key features of local control using a low-dimensional system of ordinary differential equations. Voltage-dependent gain and graded Ca²⁺ release are emergent properties of this model by virtue of the fact that model formulation is closely based on the sub-cellular basis of local control. In this current work, we have incorporated this graded release model into a prior model of guinea pig ventricular myocyte electrophysiology, metabolism, and isometric force production. The resulting integrative model predicts the experimentally observed causal relationship between action potential (AP) shape and timing of Ca²⁺ and force transients, a relationship that is not explained by models lacking the graded release property. Model results suggest that even relatively subtle changes in AP morphology that may result, for example, from remodeling of membrane transporter expression in disease or spatial variation in cell properties, may have major impact on the temporal waveform of Ca²⁺ transients, thus influencing tissue level electromechanical function.

Emergence of spiral wave activity in a mechanically heterogeneous reaction-diffusion-mechanics system

We perform a numerical study of emergent spiral wave activity in a two-dimensional reaction-diffusion-mechanics medium with a regional inhomogeneity in active and passive mechanical properties. We find that self-sustaining spiral wave activity emerges for a wide range of mechanical parameters of the inhomogeneity via five mechanisms. We classify these mechanisms, relate them to parameters of the inhomogeneity, and discuss how these results can be applied to understand the onset of cardiac arrhythmias due to regional mechanical heterogeneity.

Modeling cardiac electromechanics and mechanoelectrical coupling in dyssynchronous and failing hearts: insight from adaptive computer models

Computer models have become more and more a research tool to obtain mechanistic insight in the effects of dyssynchrony and heart failure. Increasing computational power in combination with increasing amounts of experimental and clinical data enables the development of mathematical models that describe electrical and mechanical behavior of the heart. By combining models based on data at the molecular and cellular level with models that describe organ function, so-called multi-scale models are created that describe heart function at different length and time scales. In this review, we describe basic modules that can be identified in multi-scale models of cardiac electromechanics. These modules simulate ionic membrane currents, calcium handling, excitation-contraction coupling, action potential propagation, and cardiac mechanics and hemodynamics. In addition, we discuss adaptive modeling approaches that aim to address long-term effects of diseases and therapy on growth, changes in fiber orientation, ionic membrane currents, and calcium handling. Finally, we discuss the first developments in patient-specific modeling. While current models still have shortcomings, well-chosen applications show promising results on some ultimate goals: understanding mechanisms of dyssynchronous heart failure and tuning pacing strategy to a particular patient, even before starting the therapy.

Regional increase in extracellular potassium can be arrhythmogenic due to nonuniform muscle contraction in rat ventricular muscle

In the ischemic myocardium, extracellular potassium ([K(+)](o)) increases to ≥20 mmol/l. To determine how lethal arrhythmias occur during ischemia, we investigated whether the increased spatial pattern of [K(+)](o), i.e., a regional or a global increase, affects the incidence of arrhythmias. Force, sarcomere length, membrane potential, and nonuniform intracellular Ca(2+) ([Ca(2+)](i)) were measured in rat ventricular trabeculae. A "regional" or "global" increase in [K(+)](o) was produced by exposing a restricted region of muscle to a jet of 30 mmol/l KCl or by superfusing trabeculae with a solution containing 30 mmol/l KCl, respectively. The increase in [Ca(2+)](i) (Ca(CW)) during Ca(2+) waves was measured (24°C, 3.0 mmol/l [Ca(2+)](o)). A regional increase in [K(+)](o) caused nonuniform [Ca(2+)](i) and contraction. In the presence of isoproterenol, the regional increase in [K(+)](o) induced sustained arrhythmias in 10 of 14 trabeculae, whereas the global increase did not induce such arrhythmias. During sustained arrhythmias, Ca(2+) surged within the jet-exposed region. In the absence of isoproterenol, the regional increase in [K(+)](o) increased Ca(CW), whereas the global increase decreased it. This increase in Ca(CW) with the regional increase in [K(+)](o) was not suppressed by 100 μmol/l streptomycin, whereas it was suppressed by 1) a combination of 10 μmol/l cilnidipine and 3 μmol/l SEA0400; 2) 20 mmol/l 2,3-butanedione monoxime; and 3) 10 μmol/l blebbistatin. A regional but not a global increase in [K(+)](o) induces sustained arrhythmias, probably due to nonuniform excitation-contraction coupling. The same mechanism may underlie arrhythmias during ischemia.

Electromechanical wave imaging for noninvasive mapping of the 3D electrical activation sequence in canines and humans in vivo

Cardiovascular diseases rank as America's primary killer, claiming the lives of over 41% of more than 2.4 million Americans. One of the main reasons for this high death toll is the severe lack of effective imaging techniques for screening, early detection and localization of an abnormality detected on the electrocardiogram (ECG). The two most widely used imaging techniques in the clinic are CT angiography and echocardiography with limitations in speed of application and reliability, respectively. It has been established that the mechanical and electrical properties of the myocardium change dramatically as a result of ischemia, infarction or arrhythmia; both at their onset and after survival. Despite these findings, no imaging technique currently exists that is routinely used in the clinic and can provide reliable, non-invasive, quantitative mapping of the regional, mechanical, and electrical function of the myocardium. Electromechanical Wave Imaging (EWI) is an ultrasound-based technique that utilizes the electromechanical coupling and its associated resulting strain to infer to the underlying electrical function of the myocardium. The methodology of EWI is first described and its fundamental performance is presented. Subsequent in vivo canine and human applications are provided that demonstrate the applicability of Electromechanical Wave Imaging in differentiating between sinus rhythm and induced pacing schemes as well as mapping arrhythmias. Preliminary validation with catheter mapping is also provided and transthoracic electromechanical mapping in all four chambers of the human heart is also presented demonstrating the potential of this novel methodology to noninvasively infer to both the normal and pathological electrical conduction of the heart.

At the heart of computational modelling

The link between experimental data and biophysically based mathematical models is key to computational simulation meeting its potential to provide physiological insight. However, despite the importance of this link, scrutiny and analysis of the processes by which models are parameterised from data are currently lacking. While this situation is common to many areas of physiological modelling, to provide a concrete context, we use examples drawn from detailed models of cardiac electro-mechanics. Using this biophysically detailed cohort of models we highlight the specific issues of model parameterization and propose this process can be separated into three stages: observation, fitting and validation. Finally, future research challenges and directions in this area are discussed.

Comparison of the effects of continuous and pulsatile left ventricular-assist devices on ventricular unloading using a cardiac electromechanics model

Left ventricular-assist devices (LVADs) are used to supply blood to the body of patients with heart failure. Pressure unloading is greater for counter-pulsating LVADs than for continuous LVADs. However, several clinical trials have demonstrated that myocardial recovery is similar for both types of LVAD. This study examined the contractile energy consumption of the myocardium with continuous and counter-pulsating LVAD support to ascertain the effect of the different LVADs on myocardial recovery. We used a three-dimensional electromechanical model of canine ventricles, with models of the circulatory system and an LVAD. We compared the left ventricular peak pressure (LVPP) and contractile ATP consumption between pulsatile and continuous LVADs. With the continuous and counter-pulsating LVAD, the LVPP decreased to 46 and 10%, respectively, and contractile ATP consumption decreased to 60 and 50%. The small difference between the contractile ATP consumption of these two types of LVAD may explain the comparable effects of the two types on myocardial recovery.

Electromechanical feedback with reduced cellular connectivity alters electrical activity in an infarct injured left ventricle: a finite element model study

Myocardial infarction (MI) significantly alters the structure and function of the heart. As abnormal strain may drive heart failure and the generation of arrhythmias, we used computational methods to simulate a left ventricle with an MI over the course of a heartbeat to investigate strains and their potential implications to electrophysiology. We created a fully coupled finite element model of myocardial electromechanics consisting of a cellular physiological model, a bidomain electrical diffusion solver, and a nonlinear mechanics solver. A geometric mesh built from magnetic resonance imaging (MRI) measurements of an ovine left ventricle suffering from a surgically induced anteroapical infarct was used in the model, cycled through the cardiac loop of inflation, isovolumic contraction, ejection, and isovolumic relaxation. Stretch-activated currents were added as a mechanism of mechanoelectric feedback. Elevated fiber and cross fiber strains were observed in the area immediately adjacent to the aneurysm throughout the cardiac cycle, with a more dramatic increase in cross fiber strain than fiber strain. Stretch-activated channels decreased action potential (AP) dispersion in the remote myocardium while increasing it in the border zone. Decreases in electrical connectivity dramatically increased the changes in AP dispersion. The role of cross fiber strain in MI-injured hearts should be investigated more closely, since results indicate that these are more highly elevated than fiber strain in the border of the infarct. Decreases in connectivity may play an important role in the development of altered electrophysiology in the high-stretch regions of the heart.

Cardiac electromechanical models: from cell to organ

The heart is a multiphysics and multiscale system that has driven the development of the most sophisticated mathematical models at the frontiers of computational physiology and medicine. This review focuses on electromechanical (EM) models of the heart from the molecular level of myofilaments to anatomical models of the organ. Because of the coupling in terms of function and emergent behaviors at each level of biological hierarchy, separation of behaviors at a given scale is difficult. Here, a separation is drawn at the cell level so that the first half addresses subcellular/single-cell models and the second half addresses organ models. At the subcellular level, myofilament models represent actin–myosin interaction and Ca-based activation. The discussion of specific models emphasizes the roles of cooperative mechanisms and sarcomere length dependence of contraction force, considered to be the cellular basis of the Frank–Starling law. A model of electrophysiology and Ca handling can be coupled to a myofilament model to produce an EM cell model, and representative examples are summarized to provide an overview of the progression of the field. The second half of the review covers organ-level models that require solution of the electrical component as a reaction–diffusion system and the mechanical component, in which active tension generated by the myocytes produces deformation of the organ as described by the equations of continuum mechanics. As outlined in the review, different organ-level models have chosen to use different ionic and myofilament models depending on the specific application; this choice has been largely dictated by compromises between model complexity and computational tractability. The review also addresses application areas of EM models such as cardiac resynchronization therapy and the role of mechano-electric coupling in arrhythmias and defibrillation.

Electromechanical models of the ventricles

Computational modeling has traditionally played an important role in dissecting the mechanisms for cardiac dysfunction. Ventricular electromechanical models, likely the most sophisticated virtual organs to date, integrate detailed information across the spatial scales of cardiac electrophysiology and mechanics and are capable of capturing the emergent behavior and the interaction between electrical activation and mechanical contraction of the heart. The goal of this review is to provide an overview of the latest advancements in multiscale electromechanical modeling of the ventricles. We first detail the general framework of multiscale ventricular electromechanical modeling and describe the state of the art in computational techniques and experimental validation approaches. The powerful utility of ventricular electromechanical models in providing a better understanding of cardiac function is then demonstrated by reviewing the latest insights obtained by these models, focusing primarily on the mechanisms by which mechanoelectric coupling contributes to ventricular arrythmogenesis, the relationship between electrical activation and mechanical contraction in the normal heart, and the mechanisms of mechanical dyssynchrony and resynchronization in the failing heart. Computational modeling of cardiac electromechanics will continue to complement basic science research and clinical cardiology and holds promise to become an important clinical tool aiding the diagnosis and treatment of cardiac disease.

Mapping of cardiac electrical activation with electromechanical wave imaging: an in silico-in vivo reciprocity study

BACKGROUND

Electromechanical wave imaging (EWI) is an entirely noninvasive, ultrasound-based imaging method capable of mapping the electromechanical activation sequence of the ventricles in vivo. Given the broad accessibility of ultrasound scanners in the clinic, the application of EWI could constitute a flexible surrogate for the 3-dimensional electrical activation.

OBJECTIVE

The purpose of this report is to reproduce the electromechanical wave (EW) using an anatomically realistic electromechanical model, and establish the capability of EWI to map the electrical activation sequence in vivo when pacing from different locations.

METHODS

EWI was performed in 1 canine during pacing from 3 different sites. A high-resolution dynamic model of coupled cardiac electromechanics of the canine heart was used to predict the experimentally recorded electromechanical wave. The simulated 3-dimensional electrical activation sequence was then compared with the experimental EW.

RESULTS

The electrical activation sequence and the EW were highly correlated for all pacing sites. The relationship between the electrical activation and the EW onset was found to be linear, with a slope of 1.01 to 1.17 for different pacing schemes and imaging angles.

CONCLUSION

The accurate reproduction of the EW in simulations indicates that the model framework is capable of accurately representing the cardiac electromechanics and thus testing new hypotheses. The one-to-one correspondence between the electrical activation and the EW sequences indicates that EWI could be used to map the cardiac electrical activity. This opens the door for further exploration of the technique in assisting in the early detection, diagnosis, and treatment monitoring of rhythm dysfunction.

Whole-heart modeling: applications to cardiac electrophysiology and electromechanics

Recent developments in cardiac simulation have rendered the heart the most highly integrated example of a virtual organ. We are on the brink of a revolution in cardiac research, one in which computational modeling of proteins, cells, tissues, and the organ permit linking genomic and proteomic information to the integrated organ behavior, in the quest for a quantitative understanding of the functioning of the heart in health and disease. The goal of this review is to assess the existing state-of-the-art in whole-heart modeling and the plethora of its applications in cardiac research. General whole-heart modeling approaches are presented, and the applications of whole-heart models in cardiac electrophysiology and electromechanics research are reviewed. The article showcases the contributions that whole-heart modeling and simulation have made to our understanding of the functioning of the heart. A summary of the future developments envisioned for the field of cardiac simulation and modeling is also presented. Biophysically based computational modeling of the heart, applied to human heart physiology and the diagnosis and treatment of cardiac disease, has the potential to dramatically change 21st century cardiac research and the field of cardiology.

Construction and validation of anisotropic and orthotropic ventricular geometries for quantitative predictive cardiac electrophysiology

Reaction-diffusion computational models of cardiac electrophysiology require both dynamic excitation models that reconstruct the action potentials of myocytes as well as datasets of cardiac geometry and architecture that provide the electrical diffusion tensor D, which determines how excitation spreads through the tissue. We illustrate an experimental pipeline we have developed in our laboratories for constructing and validating such datasets. The tensor D changes with location in the myocardium, and is determined by tissue architecture. Diffusion tensor magnetic resonance imaging (DT-MRI) provides three eigenvectors e(i) and eigenvalues λ(i) at each voxel throughout the tissue that can be used to reconstruct this architecture. The primary eigenvector e(1) is a histologically validated measure of myocyte orientation (responsible for anisotropic propagation). The secondary and tertiary eigenvectors (e(2) and e(3)) specify the directions of any orthotropic structure if λ(2) is significantly greater than λ(3)-this orthotropy has been identified with sheets or cleavage planes. For simulations, the components of D are scaled in the fibre and cross-fibre directions for anisotropic simulations (or fibre, sheet and sheet normal directions for orthotropic tissues) so that simulated conduction velocities match values from optical imaging or plunge electrode experiments. The simulated pattern of propagation of action potentials in the models is partially validated by optical recordings of spatio-temporal activity on the surfaces of hearts. We also describe several techniques that enhance components of the pipeline, or that allow the pipeline to be applied to different areas of research: Q ball imaging provides evidence for multi-modal orientation distributions within a fraction of voxels, infarcts can be identified by changes in the anisotropic structure-irregularity in myocyte orientation and a decrease in fractional anisotropy, clinical imaging provides human ventricular geometry and can identify ischaemic and infarcted regions, and simulations in human geometries examine the roles of anisotropic and orthotropic architecture in the initiation of arrhythmias.

Complex structure of electrophysiological gradients emerging during long-duration ventricular fibrillation in the canine heart

Long-duration ventricular fibrillation (LDVF) in the globally ischemic heart is a common setting of cardiac arrest. Electrical heterogeneities during LDVF may affect outcomes of defibrillation and resuscitation. Previous studies in large mammalian hearts have investigated the role of Purkinje fibers and electrophysiological gradients between the endocardium (Endo) and epicardium (Epi). Much less is known about gradients between the right ventricle (RV) and left ventricle (LV) and within each chamber during LDVF. We studied the transmural distribution of the VF activation rate (VFR) in the RV and LV and at the junction of RV, LV, and septum (Sep) during LDVF using plunge needle electrodes in opened-chest dogs. We also used optical mapping to analyze the Epi distribution of VFR, action potential duration (APD), and diastolic interval (DI) during LDVF in the RV and LV of isolated hearts. Transmural VFR gradients developed in both the RV and LV, with a faster VFR in Endo. Concurrently, large VFR gradients developed in Epi, with the fastest VFR in the RV-Sep junction, intermediate in the RV, and slowest in the LV. Optical mapping revealed a progressively increasing VFR dispersion within both the LV and RV, with a mosaic presence of fully inexcitable areas after 4-8 min of LDVF. The transmural, interchamber, and intrachamber VFR heterogeneities were of similar magnitude. In both chambers, the inverse of VFR was highly correlated with DI, but not APD, at all time points of LDVF. We conclude that the complex VFR gradients during LDVF in the canine heart cannot be explained solely by the distribution of Purkinje fibers and are related to regional differences in the electrical depression secondary to LDVF.

Models of stretch-activated ventricular arrhythmias

One of the most important components of mechanoelectric coupling is stretch-activated channels, sarcolemmal channels that open upon mechanical stimuli. Uncovering the mechanisms by which stretch-activated channels contribute to ventricular arrhythmogenesis under a variety of pathologic conditions is hampered by the lack of experimental methodologies that can record the 3-dimensional electromechanical activity simultaneously at high spatiotemporal resolution. Computer modeling provides such an opportunity. The goal of this review is to illustrate the utility of sophisticated, physiologically realistic, whole heart computer simulations in determining the role of mechanoelectric coupling in ventricular arrhythmogenesis. We first present the various ways by which stretch-activated channels have been modeled and demonstrate how these channels affect cardiac electrophysiologic properties. Next, we use an electrophysiologic model of the rabbit ventricles to understand how so-called commotio cordis, the mechanical impact to the precordial region of the heart, can initiate ventricular tachycardia via the recruitment of stretch-activated channels. Using the same model, we also provide mechanistic insight into the termination of arrhythmias by precordial thump under normal and globally ischemic conditions. Lastly, we use a novel anatomically realistic dynamic 3-dimensional coupled electromechanical model of the rabbit ventricles to gain insight into the role of electromechanical dysfunction in arrhythmogenesis during acute regional ischemia.

Models of cardiac tissue electrophysiology: progress, challenges and open questions

Models of cardiac tissue electrophysiology are an important component of the Cardiac Physiome Project, which is an international effort to build biophysically based multi-scale mathematical models of the heart. Models of tissue electrophysiology can provide a bridge between electrophysiological cell models at smaller scales, and tissue mechanics, metabolism and blood flow at larger scales. This paper is a critical review of cardiac tissue electrophysiology models, focussing on the micro-structure of cardiac tissue, generic behaviours of action potential propagation, different models of cardiac tissue electrophysiology, the choice of parameter values and tissue geometry, emergent properties in tissue models, numerical techniques and computational issues. We propose a tentative list of information that could be included in published descriptions of tissue electrophysiology models, and used to support interpretation and evaluation of simulation results. We conclude with a discussion of challenges and open questions.

Electromechanical wavebreak in a model of the human left ventricle

In the present report, we introduce an integrative three-dimensional electromechanical model of the left ventricle of the human heart. Electrical activity is represented by the ionic TP06 model for human cardiac cells, and mechanical activity is represented by the Niederer-Hunter-Smith active contractile tension model and the exponential Guccione passive elasticity model. These models were embedded into an anatomic model of the left ventricle that contains a detailed description of cardiac geometry and the fiber orientation field. We demonstrated that fiber shortening and wall thickening during normal excitation were qualitatively similar to experimental recordings. We used this model to study the effect of mechanoelectrical feedback via stretch-activated channels on the stability of reentrant wave excitation. We found that mechanoelectrical feedback can induce the deterioration of an otherwise stable spiral wave into turbulent wave patterns similar to that of ventricular fibrillation. We identified the mechanisms of this transition and studied the three-dimensional organization of this mechanically induced ventricular fibrillation.

In the spotlight: cardiovascular engineering

This paper discusses about the electromechanical effect of rabbit ventricals and abnormalities caused by asynchronous electrical activation in perfusion and pump function. These applications are of paramount importance to therapies that employ pacing of the heart, and particularly cardiac resynchronization therapy (CRT).