Noninvasive reconstruction of the three-dimensional ventricular activation sequence during pacing and ventricular tachycardia in the canine heart

Single-beat imaging of myocardial activation promises to aid in both cardiovascular research and clinical medicine. In the present study we validate a three-dimensional (3D) cardiac electrical imaging (3DCEI) technique with the aid of simultaneous 3D intracardiac mapping to assess its capability to localize endocardial and epicardial initiation sites and image global activation sequences during pacing and ventricular tachycardia (VT) in the canine heart. Body surface potentials were measured simultaneously with bipolar electrical recordings in a closed-chest condition in healthy canines. Computed tomography images were obtained after the mapping study to construct realistic geometry models. Data analysis was performed on paced rhythms and VTs induced by norepinephrine (NE). The noninvasively reconstructed activation sequence was in good agreement with the simultaneous measurements from 3D cardiac mapping with a correlation coefficient of 0.74 ± 0.06, a relative error of 0.29 ± 0.05, and a root mean square error of 9 ± 3 ms averaged over 460 paced beats and 96 ectopic beats including premature ventricular complexes, couplets, and nonsustained monomorphic VTs and polymorphic VTs. Endocardial and epicardial origins of paced beats were successfully predicted in 72% and 86% of cases, respectively, during left ventricular pacing. The NE-induced ectopic beats initiated in the subendocardium by a focal mechanism. Sites of initial activation were estimated to be ～7 mm from the measured initiation sites for both the paced beats and ectopic beats. For the polymorphic VTs, beat-to-beat dynamic shifts of initiation site and activation pattern were characterized by the reconstruction. The present results suggest that 3DCEI can noninvasively image the 3D activation sequence and localize the origin of activation of paced beats and NE-induced VTs in the canine heart with good accuracy. This 3DCEI technique offers the potential to aid interventional therapeutic procedures for treating ventricular arrhythmias arising from epicardial or endocardial sites and to noninvasively assess the mechanisms of these arrhythmias.     

Noninvasive electrocardiographic imaging for cardiac electrophysiology and arrhythmia

Over 7 million people worldwide die annually from erratic heart rhythms (cardiac arrhythmias), and many more are disabled. Yet there is no imaging modality to identify patients at risk, provide accurate diagnosis and guide therapy. Standard diagnostic techniques such as the electrocardiogram (ECG) provide only low-resolution projections of cardiac electrical activity on the body surface. Here we demonstrate the successful application in humans of a new imaging modality called electrocardiographic imaging (ECGI), which noninvasively images cardiac electrical activity in the heart. In ECGI, a multielectrode vest records 224 body-surface electrocardiograms; electrical potentials, electrograms and isochrones are then reconstructed on the heart's surface using geometrical information from computed tomography (CT) and a mathematical algorithm. We provide examples of ECGI application during atrial and ventricular activation and ventricular repolarization in (i) normal heart (ii) heart with a conduction disorder (right bundle branch block) (iii) focal activation initiated by right or left ventricular pacing, and (iv) atrial flutter.     

Differences in the pattern of ventricular activation in small rodents determined by morphological organization of the cardiac ventricular conduction system

Ventricular activation of the mouse heart differs significantly compared to activation in larger mammals. Knowledge of structural and functional characteristics of laboratory animals is essential for evaluation of results obtained from experiments. The present study was performed to evaluate whether the different pattern of activation is common to small rodents or unique for mice. Hearts of adult Wistar rats were isolated and Langendorff perfused. After removing the right and left ventricular free wall, extracellular activity of the septum and bundle branches (BB) was determined using a multi-terminal electrode harboring 247 terminals. Immunolabeling on cryosections was performed to assess expression and distribution of the gap junction proteins Connexin40 (Cx40), Cx43, Cx45, contractile (Desmin, alpha-actinin) and intercalated disk-related (N-cadherin, beta-catenin) proteins. Collagen distribution was assessed by Sirius Red staining. Reconstruction of the left and right bundle branch (LBB and RBB) using immuno-labeling revealed that the LBB spreads all over the septal surface. The RBB too is broad, albeit to a lesser extend than LBB. A sheet of connective tissue electrically separates the common bundle and proximal BB from the septal working myocardium. Immunolabeling revealed clear differences between the conduction system and the working myocardium with respect to expression level and distribution of the different proteins analyzed. The morphological organization of the area resulted in an electrical activation pattern of the septum comparable to what is common in larger mammals: earliest activation at the midseptum via the bundle branches. From our data we conclude that the pattern of ventricular activation in the rat heart and the structure of the conduction system fit to data described for larger mammals and differ from the different pattern previously found in mouse heart.     

Biophysical models of the heart electrical activity

Based on the fundamental knowledge of the space-temporal organization of extracellular electrical fields of the myocardium, a system for 3-D computer modeling of the cardiac electrical activity at different structural levels of the object is being developed at the Institute of Theoretical and Experimental Biophysics. The system is based on the earlier proposed and modified biophysical model of the electrocardiosignal genesis represented by a double electrical layer along the surface of the electrically active myocardium. The system combines the model for activation and repolarization of the heart ventricles; the advanced model for the evaluation of parameters of the cardiac electric field, which makes it possible to derive model electrocardiosignals both in the direct regime of calculation of the potentials and in the regime of calculation of electrocardiosignals from preliminarily determined components of the multipole equivalent electrical heart generator; a database for the model parameters and their combinations in the form of cards of simulated "patients", and a database of modeled electrocardiosignals. In the present paper (first from three within the framework of the problem), simulation methods in electrocardiology are briefly described and a biophysical model of the heart electrical activity is presented, which has made up the basis of the system for computer modeling of forward and inverse problems of the cardiac electric field. The parameters of the model are electrophysiological, anatomical, and biophysical characteristics of the heart.     

Protective Role of False Tendon in Subjects with Left Bundle Branch Block: A Virtual Population Study

False tendons (FTs) are fibrous or fibromuscular bands that can be found in both the normal and abnormal human heart in various anatomical forms depending on their attachment points, tissue types, and geometrical properties. While FTs are widely considered to affect the function of the heart, their specific roles remain largely unclear and unexplored. In this paper, we present an in silico study of the ventricular activation time of the human heart in the presence of FTs. This study presents the first computational model of the human heart that includes a FT, Purkinje network, and papillary muscles. Based on this model, we perform simulations to investigate the effect of different types of FTs on hearts with the electrical conduction abnormality of a left bundle branch block (LBBB). We employ a virtual population of 70 human hearts derived from a statistical atlas, and run a total of 560 simulations to assess ventricular activation time with different FT configurations. The obtained results indicate that, in the presence of a LBBB, the FT reduces the total activation time that is abnormally augmented due to a branch block, to such an extent that surgical implant of cardiac resynchronisation devices might not be recommended by international guidelines. Specifically, the simulation results show that FTs reduce the QRS duration at least 10 ms in 80% of hearts, and up to 45 ms for FTs connecting to the ventricular free wall, suggesting a significant reduction of cardiovascular mortality risk. In further simulation studies we show the reduction in the QRS duration is more sensitive to the shape of the heart then the size of the heart or the exact location of the FT. Finally, the model suggests that FTs may contribute to reducing the activation time difference between the left and right ventricles from 12 ms to 4 ms. We conclude that FTs may provide an alternative conduction pathway that compensates for the propagation delay caused by the LBBB. Further investigation is needed to quantify the clinical impact of FTs on cardiovascular mortality risk.     

Electromechanical analysis of infarct border zone in chronic myocardial infarction

To test the hypothesis that alterations in electrical activation sequence contribute to depressed systolic function in the infarct border zone, we examined the anatomic correlation of abnormal electromechanics and infarct geometry in the canine post-myocardial infarction (MI) heart, using a high-resolution MR-based cardiac electromechanical mapping technique. Three to eight weeks after an MI was created in six dogs, a 247-electrode epicardial sock was placed over the ventricular epicardium under thoracotomy. MI location and geometry were evaluated with delayed hyperenhancement MRI. Three-dimensional systolic strains in epicardial and endocardial layers were measured in five short-axis slices with motion-tracking MRI (displacement encoding with stimulated echoes). Epicardial electrical activation was determined from sock recordings immediately before and after the MR scans. The electrodes and MR images were spatially registered to create a total of 160 nodes per heart that contained mechanical, transmural infarct extent, and electrical data. The average depth of the infarct was 55% (SD 11), and the infarct covered 28% (SD 6) of the left ventricular mass. Significantly delayed activation (>mean + 2SD) was observed within the infarct zone. The strain map showed abnormal mechanics, including abnormal stretch and loss of the transmural gradient of radial, circumferential, and longitudinal strains, in the region extending far beyond the infarct zone. We conclude that the border zone is characterized by abnormal mechanics directly coupled with normal electrical depolarization. This indicates that impaired function in the border zone is not contributed by electrical factors but results from mechanical interaction between ischemic and normal myocardium.     

Cooperative interactions between R531 and acidic residues in the voltage sensing module of hERG1 channels.

HERG1 K(+) channels are critical for modulating the duration of the cardiac action potential. The role of hERG1 channels in maintaining electrical stability in the heart derives from their unusual gating properties: slow activation and fast inactivation. HERG1 channel inactivation is intrinsically voltage sensitive and is not coupled to activation in the same way as in the Shaker family of K(+) channels. We recently proposed that the S4 transmembrane domain functions as the primary voltage sensor for hERG1 activation and inactivation and that distinct regions of S4 contribute to each gating process. In this study, we tested the hypothesis that S4 rearrangements underlying activation and inactivation gating may be associated with distinct cooperative interactions between a key residue in the S4 domain (R531) and acidic residues in neighboring regions (S1 - S3 domains) of the voltage sensing module. Using double-mutant cycle analysis, we found that R531 was energetically coupled to all acidic residues in S1-S3 during activation, but was coupled only to acidic residues near the extracellular portion of S2 and S3 (D456, D460 and D509) during inactivation. We propose that hERG1 activation involves a cooperative conformational change involving the entire voltage sensing module, while inactivation may involve a more limited interaction between R531 and D456, D460 and D509.     

Biophysical models of cardiac electrical activity

A system for 3D simulation of heart electrical activity at different structural levels based on fundamental knowledge on the spatiotemporal organization of extracellular electric fields in the myocardium is being developed at the Institute of Theoretical and Experimental Biophysics of the Russian Academy of Sciences. The system is based on a biophysical model of the genesis of electrocardiosignals (ECSs) in the form of a double electric layer on the surface of the electrically active myocardium, which was proposed earlier and then modified. The system combines a model of the activation and repolarization of the heart ventricles, an advanced model for determining the parameters of the heart electric field, which makes it possible to obtain model ECSs both by direct calculation of the potentials and calculation of ECSs from preliminarily determined components of a multipole equivalent heart generator, a database of model parameters and their combinations in the form of cards of simulated “patients,” and a database of simulated ECSs. This paper (the first in a series of three on the subject) briefly describes simulation methods used in electrocardiology and the biophysical model of heart electrical activity that forms the basis of the system for computer simulation of direct and inverse problems concerning the heart electric field. Electrophysiological, anatomical, and biophysical characteristics of the heart are the parameters of the model.     

Three-dimensional analysis of regional cardiac function: a model of rabbit ventricular anatomy

The three-dimensional geometry and anisotropic properties of the heart give rise to nonhomogeneous distributions of stress, strain, electrical activation and repolarization. In this article we review the ventricular geometry and myofiber architecture of the heart, and the experimental and modeling studies of three-dimensional cardiac mechanics and electrophysiology. The development of a three-dimensional finite element model of the rabbit ventricular geometry and fiber architecture is described in detail. Finally, we review the experimental results, from the level of the cell to the intact organ, which motivate the development of coupled three-dimensional models of cardiac electromechanics and mechanoelectric feedback.     

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.     

Fast Simulation of Mechanical Heterogeneity in the Electrically Asynchronous Heart Using the MultiPatch Module

Cardiac electrical asynchrony occurs as a result of cardiac pacing or conduction disorders such as left bundle-branch block (LBBB). Electrically asynchronous activation causes myocardial contraction heterogeneity that can be detrimental for cardiac function. Computational models provide a tool for understanding pathological consequences of dyssynchronous contraction. Simulations of mechanical dyssynchrony within the heart are typically performed using the finite element method, whose computational intensity may present an obstacle to clinical deployment of patient-specific models. We present an alternative based on the CircAdapt lumped-parameter model of the heart and circulatory system, called the MultiPatch module. Cardiac walls are subdivided into an arbitrary number of patches of homogeneous tissue. Tissue properties and activation time can differ between patches. All patches within a wall share a common wall tension and curvature. Consequently, spatial location within the wall is not required to calculate deformation in a patch. We test the hypothesis that activation time is more important than tissue location for determining mechanical deformation in asynchronous hearts. We perform simulations representing an experimental study of myocardial deformation induced by ventricular pacing, and a patient with LBBB and heart failure using endocardial recordings of electrical activation, wall volumes, and end-diastolic volumes. Direct comparison between simulated and experimental strain patterns shows both qualitative and quantitative agreement between model fibre strain and experimental circumferential strain in terms of shortening and rebound stretch during ejection. Local myofibre strain in the patient simulation shows qualitative agreement with circumferential strain patterns observed in the patient using tagged MRI. We conclude that the MultiPatch module produces realistic regional deformation patterns in the asynchronous heart and that activation time is more important than tissue location within a wall for determining myocardial deformation. The CircAdapt model is therefore capable of fast and realistic simulations of dyssynchronous myocardial deformation embedded within the closed-loop cardiovascular system.     

Activation of the right atrioventricular valve in the avian heart

Using intramural needle electrode, in situ determinations have been made of the moment of activation of the right muscular atrioventricular valve during electrical systole of the heart. In the hen, Gallus domesticus, the electrical activation of the valve was observed 8.8 +/- 1.4 ms (mean QRS duration 25 ms) whereas in the pigeon Columba livia--5 +/- 1 ms (QRS duration 17 ms) after the beginning of ventricular depolarization, i.e. during excitation of the main bulk of the myocardium of the free wall of the right ventricle.     

Cardiac baroreflex facilitation evoked by hypothalamus and prefrontal cortex stimulation: role of the nucleus tractus solitarius 5-HT2A receptors

We previously showed that serotonin (5-HT2) receptor activation in the nucleus of the tractus solitarius (NTS) produced hypotension, bradycardia, and facilitation of the baroreflex bradycardia. Activation of the preoptic area (POA) of the hypothalamus, which is involved in shock-evoked passive behaviors, induces similar modifications. In addition, previous studies showed that blockade of the infralimbic (IL) part of the medial prefrontal cortex, which sends projections to POA, produced an inhibitory influence on the baroreflex cardiac response. Thus, to assess the possible implication of NTS 5-HT2 receptors in passive cardiovascular responses, we analyzed in anesthetized rats the effects of NTS inhibition and NTS 5-HT2 receptor blockade on the cardiovascular modifications induced by chemical (0.3 M D,L-homocysteic acid) and electrical (50 Hz, 150-200 microA) stimulation of IL or POA. Intra-NTS microinjections of muscimol, a GABAA receptor agonist, prevented the decreases in blood pressure and heart rate normally evoked by IL or POA activation. In addition, we found that intra-NTS microinjection of R(+)-alpha-(2,3-dimethoxyphenyl)-1-[2-(4-fluorophenylethyl)]-4-piperidine-methanol, a specific 5-HT2A receptor antagonist, did not affect the decreases in cardiovascular baseline parameters induced by IL or POA stimulation but prevented the facilitation of the aortic baroreflex bradycardia normally observed during IL (+65 and +60%) or POA (+70 and +69%) electrical and chemical stimulation, respectively. These results show that NTS 5-HT2A receptors play a key role in the enhancement of the cardiac response of the baroreflex but not in the changes in basal heart rate and blood pressure induced by IL or POA stimulation.     

Synthesis of voltage-sensitive fluorescence signals from three-dimensional myocardial activation patterns

Voltage-sensitive fluorescent dyes are commonly used to measure cardiac electrical activity. Recent studies indicate, however, that optical action potentials (OAPs) recorded from the myocardial surface originate from a widely distributed volume beneath the surface and may contain useful information regarding intramural activation. The first step toward obtaining this information is to predict OAPs from known patterns of three-dimensional (3-D) electrical activity. To achieve this goal, we developed a two-stage model in which the output of a 3-D ionic model of electrical excitation serves as the input to an optical model of light scattering and absorption inside heart tissue. The two-stage model permits unique optical signatures to be obtained for given 3-D patterns of electrical activity for direct comparison with experimental data, thus yielding information about intramural electrical activity. To illustrate applications of the model, we simulated surface fluorescence signals produced by 3-D electrical activity during epicardial and endocardial pacing. We discovered that OAP upstroke morphology was highly sensitive to the transmural component of wave front velocity and could be used to predict wave front orientation with respect to the surface. These findings demonstrate the potential of the model for obtaining useful 3-D information about intramural propagation.     

Physiological heart activation by adrenaline involves parallel activation of ATP usage and supply

During low-to-high work transition in adult mammalian heart in vivo the concentrations of free ADP, ATP, PCr (phosphocreatine), P(i) and NADH are essentially constant, in striking contrast with skeletal muscle. The direct activation by calcium ions of ATP usage and feedback activation of ATP production by ADP (and P(i)) alone cannot explain this perfect homoeostasis. A comparison of the response to adrenaline (increase in rate-pressure product and [PCr]) of the intact beating perfused rat heart with the elasticities of the PCr producer and consumer to PCr concentration demonstrated that both the ATP/PCr-producing block and ATP/PCr-consuming block are directly activated to a similar extent during physiological heart activation. Our finding constitutes a direct evidence for the parallel-activation mechanism of the regulation of oxidative phosphorylation in heart postulated previously in a theoretical way.     

Patient-specific modeling of dyssynchronous heart failure: A case study

The development and clinical use of patient-specific models of the heart is now a feasible goal. Models have the potential to aid in diagnosis and support decision-making in clinical cardiology. Several groups are now working on developing multi-scale models of the heart for understanding therapeutic mechanisms and better predicting clinical outcomes of interventions such as cardiac resynchronization therapy. Here we describe the methodology for generating a patient-specific model of the failing heart with a myocardial infarct and left ventricular bundle branch block. We discuss some of the remaining challenges in developing reliable patient-specific models of cardiac electromechanical activity, and identify some of the main areas for focusing future research efforts. Key challenges include: efficiently generating accurate patient-specific geometric meshes and mapping regional myofiber architecture to them; modeling electrical activation patterns based on cellular alterations in human heart failure, and estimating regional tissue conductivities based on clinically available electrocardiographic recordings; estimating unloaded ventricular reference geometry and material properties for biomechanical simulations; and parameterizing systemic models of circulatory dynamics from available hemodynamic measurements.     

CFD simulation of flow through heart: a perspective review

The heart is an organ which pumps blood around the body by contraction of muscular wall. There is a coupled system in the heart containing the motion of wall and the motion of blood fluid; both motions must be computed simultaneously, which make biological computational fluid dynamics (CFD) difficult. The wall of the heart is not rigid and hence proper boundary conditions are essential for CFD modelling. Fluid-wall interaction is very important for real CFD modelling. There are many assumptions for CFD simulation of the heart that make it far from a real model. A realistic fluid-structure interaction modelling the structure by the finite element method and the fluid flow by CFD use more realistic coupling algorithms. This type of method is very powerful to solve the complex properties of the cardiac structure and the sensitive interaction of fluid and structure. The final goal of heart modelling is to simulate the total heart function by integrating cardiac anatomy, electrical activation, mechanics, metabolism and fluid mechanics together, as in the computational framework.     

Optical mapping of activation patterns in an animal model of congenital heart block

Congenital heart block (CHB) is associated with high mortality and affects children of mothers with autoantibodies (IgG) to ribonucleoproteins SSB/La and SSA/Ro. IgG from mothers of children with CHB (positive IgG) was used to assess activation patterns in both the right atrium (RA) and right ventricle (RV) of Langendorff-perfused young rabbit hearts. Optical action potentials (AP) were obtained by using a 124-site photodiode array with 4-[-[2-(di-n-butylamino)-6-naphthyl]vinyl]pyridinium. Optical APs were recorded to simultaneously image activation patterns from the RA and RV. Perfusion of positive IgG (800--1,200 micro resulted in sinus bradycardia and varying degrees of heart block. Activation maps revealed marked conduction delay at the sinoatrial junction but only minor changes in overall atrial and ventricular activation patterns. No conduction disturbances were seen in the presence of IgG from mothers with healthy children. In conclusion, besides atrioventricular (AV) block, positive IgG induces sinus bradycardia. These results establish that the sequelae of CHB are associated with impaired intrasinus and/or sinoatrial conduction. The findings raise the possibility that sinus bradycardia in the developing heart may indicate the potential for AV conduction disturbances.