Viatcheslav Gurev

Models of cardiac electromechanics based on individual hearts imaging data Image-based electromechanical models of the heart

Current multi-scale computational models of ventricular electromechanics describe the full process of cardiac contraction on both the micro- and macro- scales including: the depolarization of cardiac cells, the release of calcium from intracellular stores, tension generation by cardiac myofilaments, and mechanical contraction of the whole heart. Such models are used to reveal basic mechanisms of cardiac contraction as well as the mechanisms of cardiac dysfunction in disease conditions. In this paper, we present a methodology to construct finite element electromechanical models of ventricular contraction with anatomically accurate ventricular geometry based on magnetic resonance and diffusion tensor magnetic resonance imaging of the heart. The electromechanical model couples detailed representations of the cardiac cell membrane, cardiac myofilament dynamics, electrical impulse propagation, ventricular contraction, and circulation to simulate the electrical and mechanical activity of the ventricles. The utility of the model is demonstrated in an example simulation of contraction during sinus rhythm using a model of the normal canine ventricles.

Towards predictive modelling of the electrophysiology of the heart

The simulation of cardiac electrical function is an example of a successful integrative multiscale modelling approach that is directly relevant to human disease. Today we stand at the threshold of a new era, in which anatomically detailed, tomographically reconstructed models are being developed that integrate from the ion channel to the electromechanical interactions in the intact heart. Such models hold high promise for interpretation of clinical and physiological measurements, for improving the basic understanding of the mechanisms of dysfunction in disease, such as arrhythmias, myocardial ischaemia and heart failure, and for the development and performance optimization of medical devices. The goal of this article is to present an overview of current state‐of‐art advances towards predictive computational modelling of the heart as developed recently by the authors of this article. We first outline the methodology for constructing electrophysiological models of the heart. We then provide three examples that demonstrate the use of these models, focusing specifically on the mechanisms for arrhythmogenesis and defibrillation in the heart. These include: (1) uncovering the role of ventricular structure in defibrillation; (2) examining the contribution of Purkinje fibres to the failure of the shock; and (3) using magnetic resonance imaging reconstructed heart models to investigate the re‐entrant circuits formed in the presence of an infarct scar.

Mechanisms of Mechanically Induced Spontaneous Arrhythmias in Acute Regional Ischemia

Rationale: Although ventricular premature beats (VPBs) during acute regional ischemia have been linked to mechanical stretch of ischemic tissue, whether and how ischemia-induced mechanical dysfunction can induce VPBs and facilitate their degradation into reentrant arrhythmias has not been yet addressed. Objective: This study used a novel multiscale electromechanical model of the rabbit ventricles to investigate the origin of and the substrate for spontaneous arrhythmias arising from ischemia-induced electrophysiological and mechanical changes. Methods and Results: Two stages of ischemia were simulated. Dynamic mechanoelectrical feedback was modeled as spatially and temporally nonuniform membrane currents through mechanosensitive channels, the conductances of which depended on local strain rate. Our results reveal that both strains and strain rates were significantly larger in the central ischemic zone than in the border zone. However, in both ischemia stages, a VPB originated from the ischemic border in the left ventricular apical endocardium because of mechanically induced suprathreshold depolarizations. It then traveled fully intramurally until emerging from the ischemic border on the anterior epicardium. Reentry was formed only in the advanced ischemia stage as the result of a widened temporal excitable gap. Mechanically induced delayed afterdepolarization-like events contributed to the formation of reentry by further decreasing the already reduced-by-hyperkalemia local excitability, causing extended conduction block lines and slowed conduction in the ischemic region. Conclusions: Mechanically induced membrane depolarizations in the ischemic region are the mechanism by which mechanical activity contributes to both the origin of and substrate for spontaneous arrhythmias under the conditions of acute regional ischemia.

Distribution of Electromechanical Delay in the Heart: Insights from a Three-Dimensional Electromechanical Model

In the intact heart, the distribution of electromechanical delay (EMD), the time interval between local depolarization and myocyte shortening onset, depends on the loading conditions. The distribution of EMD throughout the heart remains, however, unknown because current experimental techniques are unable to evaluate three-dimensional cardiac electromechanical behavior. The goal of this study was to determine the three-dimensional EMD distributions in the intact ventricles for sinus rhythm (SR) and epicardial pacing (EP) by using a new, to our knowledge, electromechanical model of the rabbit ventricles that incorporates a biophysical representation of myofilament dynamics. Furthermore, we aimed to ascertain the mechanisms that underlie the specific three-dimensional EMD distributions. The results revealed that under both conditions, the three-dimensional EMD distribution is nonuniform. During SR, EMD is longer at the epicardium than at the endocardium, and is greater near the base than at the apex. After EP, the three-dimensional EMD distribution is markedly different; it also changes with the pacing rate. For both SR and EP, late-depolarized regions were characterized with significant myofiber prestretch caused by the contraction of the early-depolarized regions. This prestretch delays myofiber-shortening onset, and results in a longer EMD, giving rise to heterogeneous three-dimensional EMD distributions.

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.

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.

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.

Verification of cardiac mechanics software: benchmark problems and solutions for testing active and passive material behaviour

Models of cardiac mechanics are increasingly used to investigate cardiac physiology. These models are characterized by a high level of complexity, including the particular anisotropic material properties of biological tissue and the actively contracting material. A large number of independent simulation codes have been developed, but a consistent way of verifying the accuracy and replicability of simulations is lacking. To aid in the verification of current and future cardiac mechanics solvers, this study provides three benchmark problems for cardiac mechanics. These benchmark problems test the ability to accurately simulate pressure-type forces that depend on the deformed objects geometry, anisotropic and spatially varying material properties similar to those seen in the left ventricle and active contractile forces. The benchmark was solved by 11 different groups to generate consensus solutions, with typical differences in higher-resolution solutions at approximately 0.5%, and consistent results between linear, quadratic and cubic finite elements as well as different approaches to simulating incompressible materials. Online tools and solutions are made available to allow these tests to be effectively used in verification of future cardiac mechanics software.

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.

Towards real-time simulation of cardiac electrophysiology in a human heart at high resolution

We have developed the capability to rapidly simulate cardiac electrophysiological phenomena in a human heart discretised at a resolution comparable with the length of a cardiac myocyte. Previous scientific investigation has generally invoked simplified geometries or coarse-resolution hearts, with simulation duration limited to 10s of heartbeats. Using state-of-the-art high-performance computing techniques coupled with one of the most powerful computers available (the 20 PFlop/s IBM BlueGene/Q at Lawrence Livermore National Laboratory), high-resolution simulation of the human heart can now be carried out over 1200 times faster compared with published results in the field. We demonstrate the utility of this capability by simulating, for the first time, the formation of transmural re-entrant waves in a 3D human heart. Such wave patterns are thought to underlie Torsades de Pointes, an arrhythmia that indicates a high risk of sudden cardiac death. Our new simulation capability has the potential to impact a multitude of applications in medicine, pharmaceuticals and implantable devices.