Mechanisms of mechanically induced spontaneous arrhythmias in acute regional ischemia

Ectopic foci do not co-locate with ventricular epicardial stretch during early acute regional ischemia in isolated pig hearts

Ectopic activation during early acute regional ischemia may initiate fatal reentrant arrhythmias. However, the origin of this ectopy remains poorly understood. Studies suggest that systolic stretch arising from dyskinesia in ischemic tissue may cause ectopic depolarization due to cardiac mechanosensitivity. The aim of this study was to investigate the link between mechanical stretch and ectopic electrical activation during early acute regional ischemia. We used a recently developed optical mapping technique capable of simultaneous imaging of mechanical deformation and electrical activation in isolated hearts. Eight domestic swine hearts were prepared in left ventricular working mode (LVW), in which the left ventricle was loaded and contracting. In an additional eight non-working (NW) hearts, contraction was pharmacologically suppressed with blebbistatin and the left ventricle was not loaded. In both groups, the left anterior descending coronary artery was tied below the first diagonal branch. Positive mechanical stretch (bulging) during systole was observed in the ischemic zones of LVW, but not NW, hearts. During ischemia phase 1a (0-15 min post-occlusion), LVW hearts had more ectopic beats than NW hearts (median: 19, interquartile range: 10-28 vs. median: 2, interquartile range: 1-6; p = 0.02); but the difference during phase 1b (15-60 min post-occlusion) was not significant (median: 27, interquartile range: 22-42 vs. median: 16, interquartile range: 12-31; p = 0.37). Ectopic beats arose preferentially from the ischemic border zone in both groups (p < 0.01). In LVW hearts, local mechanical stretch was only occasionally co-located with ectopic foci (9 of 69 ectopic beats). Despite the higher rate of ectopy observed in LVW hearts during ischemia phase 1a, the ectopic beats generally did not arise by the hypothesized mechanism in which ectopic foci are generated by co-local epicardial mechanical stretch.

Computational modelling of mechano-electric feedback and its arrhythmogenic effects in human ventricular models

The current study aims to investigate the role of mechano-electric feedback (MEF) in healthy cardiac cycles and in cardiac arrhythmia using human ventricular models. The numerical formulation of stretch-activated channels (SACs) in terms of the fibre stretch of the myocardium is incorporated into the modified Hill model that describes the myocardium as an electro-visco-active material. Additionally, we propose models of SACs formulated in terms of the rate of stretch along fibre direction and the stretch along sheet direction. We analyze the effect of the three different models for SACs and different material properties on the regular cycles by using electrocardiogram and volume-time curves, and show that the each model of SACs has regionally different influences on the heart model. Moreover, we simulate 'commotio cordis' and 'precordial thump' and demonstrate that MEF plays a major role in the occurrence of fibrillation and defibrillation in the absence of the structural cardiac damage. Furthermore, we study the role of MEF in premature ventricular contraction when the blood pressure is disturbed.

Cardiac electro-mechanical activity in a deforming human cardiac tissue: modeling, existence-uniqueness, finite element computation and application to multiple ischemic disease

In this study, the cardiac electro-mechanical model in a deforming domain is taken with the addition of mechanical feedback and stretch-activated channel current coupled with the ten Tusscher human ventricular cell level model that results in a coupled PDE-ODE system. The existence and uniqueness of such a coupled system in a deforming domain is proved. At first, the existence of a solution is proved in the deformed domain. The local existence of the solution is proved using the regularization and the Faedo-Galerkin technique. Then, the global existence is proved using the energy estimates in appropriate Banach spaces, Gronwall lemma, and the compactness procedure. The existence of the solution in an undeformed domain is proved using the lower semi-continuity of the norms. Uniqueness is proved using Young's inequality, Gronwall lemma, and the Cauchy-Schwartz inequality. For the application purpose, this model is applied to understand the electro-mechanical activity in ischemic cardiac tissue. It also takes care of the development of active tension, conductive, convective, and ionic feedback. The Second Piola-Kirchoff stress tensor arising in Lagrangian mapping between reference and moving frames is taken as a combination of active, passive, and volumetric components. We investigated the effect of varying strength of hyperkalemia and hypoxia, in the ischemic subregions of human cardiac tissue with local multiple ischemic subregions, on the electro-mechanical activity of healthy and ischemic zones. This system is solved numerically using the [Formula: see text] finite element method in space and the implicit-explicit Euler method in time. Discontinuities arising with the modeled multiple ischemic regions are treated to the desired order of accuracy by a simple regularization technique using the interpolating polynomials. We examined the cardiac electro-mechanical activity for several cases in multiple hyperkalemic and hypoxic human cardiac tissue. We concluded that local multiple ischemic subregions severely affect the cardiac electro-mechanical activity more, in terms of action potential (v) and mechanical parameters, intracellular calcium ion concentration [Formula: see text], active tension ([Formula: see text]), stretch ([Formula: see text]) and stretch rate ([Formula: see text]), of a healthy cell in its vicinity, compared to a single Hyperkalemic or Hypoxic subregion. The four moderate hypoxically generated ischemic subregions affect the waveform of the stretch along the fiber and the stretch rate more than a single severe ischemic subregion.

Roles of stretch-activated channels and NADPH oxidase 2 in the induction of twitch contraction by muscle stretching in rat ventricular muscle

Mechano-electric feedback means that muscle stretching causes depolarization of membrane potential. We investigated whether muscle stretching induces action potential and twitch contraction with a threshold of sarcomere length (SL) and what roles stretch-activated channels (SACs) and stretch-activated NADPH oxidase (X-ROS signaling) play in the induction. Trabeculae were obtained from the right ventricles of rat hearts. Force, SL, and [Ca²⁺]_i were measured. Various degrees of stretching from the SL of 2.0 μm were applied 0.5 s after the last stimulus of the electrical train with 0.4-s intervals for 7.5 s. The SL_twitch was defined as the minimal SL at which twitch contraction was induced by the stretching. Muscle stretching induced twitch contraction with a threshold of SL at 0.4-s stimulus intervals ([Ca²⁺]_o = 0.7 mmol/L). The SL_twitch was not changed by increasing the stimulus intervals and [Ca²⁺]_o and by adding 1 μmol/L isoproterenol. The SL_twitch was not changed by adding 10 μmol/L Gd³⁺, 100 μmol/L or 200 μmol/L streptomycin, and 5 μmol/L GsMTx4. The SL_twitch was not changed by adding 1 μmol/L ryanodine and 3 μmol/L diphenyleneiodonium chloride. In contrast, the SL_twitch was increased by elevating extracellular K⁺ from 5 to 10 mmol/L and by adding the stretching during the refractory period of membrane potential. The addition of the stretching-induced twitch contraction more frequently induced arrhythmias. These results suggest that muscle stretching can induce twitch contraction with a threshold of SL and concern the occurrence of arrhythmias and that SACs and X-ROS signaling play no roles in the induction.

Heart-on-Chip for Combined Cellular Dynamics Measurements and Computational Modeling Towards Clinical Applications

Organ-on-chip or micro-engineered three-dimensional cellular or tissue models are increasingly implemented in the study of cardiovascular pathophysiology as alternatives to traditional in vitro cell culture. Drug induced cardiotoxicity is a key issue in drug development pipelines, but the current in vitro and in vivo studies suffer from inter-species differences, high costs, and lack of reliability and accuracy in predicting cardiotoxicity. Microfluidic heart-on-chip devices can impose a paradigm shift to the current tools. They can not only recapitulate cardiac tissue level functionality and the communication between cells and extracellular matrices but also allow higher throughput studies conducive to drug screening especially with their added functionalities or sensors that extract disease-specific phenotypic, genotypic, and electrophysiological information in real-time. Such electrical and mechanical components can tailor the electrophysiology and mechanobiology of the experiment to better mimic the in vivo condition as well. Recent advancements and challenges are reviewed in the fabrication, functionalization and sensor assisted mechanical and electrophysiological measurements, numerical and computational modeling of cardiomyocytes' behavior, and the clinical applications in drug screening and disease modeling. This review concludes with the current challenges and perspectives on the future of such organ-on-chip platforms.

A clinical-in silico study on the effectiveness of multipoint bicathodic and cathodic-anodal pacing in cardiac resynchronization therapy

Up to one-third of patients undergoing cardiac resynchronization therapy (CRT) are nonresponders. Multipoint bicathodic and cathodic-anodal left ventricle (LV) stimulations could overcome this clinical challenge, but their effectiveness remains controversial. Here we evaluate the performance of such stimulations through both in vivo and in silico experiments, the latter based on computer electromechanical modeling. Seven patients, all candidates for CRT, received a quadripolar LV lead. Four stimulations were tested: right ventricular (RVS); conventional single point biventricular (S-BS); multipoint biventricular bicathodic (CC-BS) and multipoint biventricular cathodic-anodal (CA-BS). The following parameters were processed: QRS duration; maximal time derivative of arterial pressure (dPdt_max); systolic arterial pressure (P_sys); and stroke volume (SV). Echocardiographic data of each patient were then obtained to create an LV geometric model. Numerical simulations were based on a strongly coupled Bidomain electromechanical coupling model. Considering the in vivo parameters, when comparing S-BS to RVS, there was no significant decrease in SV (from 45 ± 11 to 44 ± 20 ml) and 6% and 4% increases of dPdt_max and P_sys, respectively. Focusing on in silico parameters, with respect to RVS, S-BS exhibited a significant increase of SV, dPdt_max and P_sys. Neither the in vivo nor in silico results showed any significant hemodynamic and electrical difference among S-BS, CC-BS and CA-BS configurations. These results show that CC-BS and CA-BS yield a comparable CRT performance, but they do not always yield improvement in terms of hemodynamic parameters with respect to S-BS. The computational results confirmed the in vivo observations, thus providing theoretical support to the clinical experiments.

Cardiac Mechano-Electric Coupling: Acute Effects of Mechanical Stimulation on Heart Rate and Rhythm

The heart is vital for biological function in almost all chordates, including humans. It beats continually throughout our life, supplying the body with oxygen and nutrients while removing waste products. If it stops, so does life. The heartbeat involves precise coordination of the activity of billions of individual cells, as well as their swift and well-coordinated adaption to changes in physiological demand. Much of the vital control of cardiac function occurs at the level of individual cardiac muscle cells, including acute beat-by-beat feedback from the local mechanical environment to electrical activity (as opposed to longer term changes in gene expression and functional or structural remodeling). This process is known as mechano-electric coupling (MEC). In the current review, we present evidence for, and implications of, MEC in health and disease in human; summarize our understanding of MEC effects gained from whole animal, organ, tissue, and cell studies; identify potential molecular mediators of MEC responses; and demonstrate the power of computational modeling in developing a more comprehensive understanding of ‟what makes the heart tick.ˮ.

Mechano-electric and mechano-chemo-transduction in cardiomyocytes

Cardiac excitation-contraction (E-C) coupling is influenced by (at least) three dynamic systems that couple and feedback to one another (see Abstract Figure). Here we review the mechanical effects on cardiomyocytes that include mechano-electro-transduction (commonly referred to as mechano-electric coupling, MEC) and mechano-chemo-transduction (MCT) mechanisms at cell and molecular levels which couple to Ca2+ -electro and E-C coupling reviewed elsewhere. These feedback loops from muscle contraction and mechano-transduction to the Ca2+ homeodynamics and to the electrical excitation are essential for understanding the E-C coupling dynamic system and arrhythmogenesis in mechanically loaded hearts. This white paper comprises two parts, each reflecting key aspects from the 2018 UC Davis symposium: MEC (how mechanical load influences electrical dynamics) and MCT (how mechanical load alters cell signalling and Ca2+ dynamics). Of course, such separation is artificial since Ca2+ dynamics profoundly affect ion channels and electrogenic transporters and vice versa. In time, these dynamic systems and their interactions must become fully integrated, and that should be a goal for a comprehensive understanding of how mechanical load influences cell signalling, Ca2+ homeodynamics and electrical dynamics. In this white paper we emphasize current understanding, consensus, controversies and the pressing issues for future investigations. Space constraints make it impossible to cover all relevant articles in the field, so we will focus on the topics discussed at the symposium.

Computational prediction of the effect of D172N KCNJ2 mutation on ventricular pumping during sinus rhythm and reentry

The understanding of cardiac arrhythmia under genetic mutations has grown in interest among researchers. Previous studies focused on the effect of the D172N mutation on electrophysiological behavior. In this study, we analyzed not only the electrophysiological activity but also the mechanical responses during normal sinus rhythm and reentry conditions by using computational modeling. We simulated four different ventricular conditions including normal case of ten Tusscher model 2006 (TTM), wild-type (WT), heterozygous (WT/D172N), and homozygous D172N mutation. The 2D simulation result (in wire-shaped mesh) showed the WT/D172N and D172N mutation shortened the action potential duration by 14%, and by 23%, respectively. The 3D electrophysiological simulation results showed that the electrical wavelength between TTM and WT conditions were identical. Under sinus rhythm condition, the WT/D172N and D172N reduced the pumping efficacy with a lower left ventricle (LV) and aortic pressures, stroke volume, ejection fraction, and cardiac output. Under the reentry conditions, the WT condition has a small probability of reentry. However, in the event of reentry, WT has shown the most severe condition. Furthermore, we found that the position of the rotor or the scroll wave substantially influenced the ventricular pumping efficacy during arrhythmia. If the rotor stays in the LV, it will cause very poor pumping performance. Graphical Abstract A model of a ventricular electromechanical system. This whole model was established to observe the effect of D172N KCNJ2 mutation on ventricular pumping behavior during sinus rhythm and reentry conditions. The model consists of two components; electrical component and mechanical component. The electrophysiological model based on ten Tusscher et al. with the IK1 D172N KCNJ2 mutation, and the myofilament dynamic (cross-bridge) model based on Rice et al. study. The 3D electrical component is a ventricular geometry based on MRI which composed of nodes representing single-cell with electrophysiological activation. The 3D ventricular mechanic is a finite element mesh composed of single-cells myofilament dynamic model. Both components were coupled with Ca2+ concentration. We used Gaussian points for the calcium interpolation from the electrical mesh to the mechanical mesh.

The role of desmin alterations in mechanical electrical feedback in heart failure

AIM

Mechanoelectric feedback (MEF) was related to malignant arrhythmias in heart failure (HF). Desmin is a cytoskeleton protein and could be involved in MEF as a mechanoelectrical transducer. In this study, we will discuss the role of desmin alterations in mechanical electrical feedback in heart failure and its mechanisms.

METHODS

We used both an in vivo rat model and an in vitro cardiomyocyte model to address this issue. For the in vivo experiments, we establish a sham group, an HF group, streptomycin (SM) group, and an MDL-28170 group. The occurrence of ventricular arrhythmias (VA) was recorded in each group. For the in vitro cardiomyocyte model, we established an NC group, a si-desmin group, and a si-desmin + NBD IKK group. The expression of desmin, IKKβ, p-IKKβ, IKBα, p-NF-κB, and SERCA2 were detected in both in vivo and in vitro experiments. The content of Ca²⁺ in cytoplasm and sarcoplasmic were detected by confocal imaging in vitro experiments.

RESULTS

An increased number of VAs were found in the HF group. SM and MDL-28170 can reduce desmin breakdown and the number of VAs in heart failure. The knockdown of desmin in the cardiomyocyte can activate the NF-κB pathway, decrease the level of SERCA2, and result in abnormal distribution of Ca²⁺. While treatment with NF-κB inhibitor can elevate the level of SERCA2 and alleviate the abnormal distribution of Ca²⁺.

SIGNIFICANCE

Overall, desmin may participate in MEF through the NF-κB pathway. This study provides a potential therapeutic target for VA in HF.

Interaction of the Mechano-Electrical Feedback With Passive Mechanical Models on a 3D Rat Left Ventricle: A Computational Study

In this paper, we are investigating the interaction between different passive material models and the mechano-electrical feedback (MEF) in cardiac modeling. Various types of passive mechanical laws (nearly incompressible/compressible, polynomial/exponential-type, transversally isotropic/orthotropic material models) are integrated in a fully coupled electromechanical model in order to study their specific influence on the overall MEF behavior. Our computational model is based on a three-dimensional (3D) geometry of a healthy rat left ventricle reconstructed from magnetic resonance imaging (MRI). The electromechanically coupled problem is solved using a fully implicit finite element-based approach. The effects of different passive material models on the MEF are studied with the help of numerical examples. It turns out that there is a significant difference between the behavior of the MEF for compressible and incompressible material models. Numerical results for the incompressible models exhibit that a change in the electrophysiology can be observed such that the transmembrane potential (TP) is unable to reach the resting state in the repolarization phase, and this leads to non-zero relaxation deformations. The most significant and strongest effects of the MEF on the rat cardiac muscle response are observed for the exponential passive material law.

4D cardiac electromechanical activation imaging

Cardiac abnormalities, a major cause of morbidity and mortality, affect millions of people worldwide. Despite the urgent clinical need for early diagnosis, there is currently no noninvasive technique that can infer to the electrical function of the whole heart in 3D and thereby localize abnormalities at the point of care. Here we present a new method for noninvasive 4D mapping of the cardiac electromechanical activity in a single heartbeat for heart disease characterization such as arrhythmia and infarction. Our novel technique captures the 3D activation wave of the heart in vivo using high volume-rate (500 volumes per second) ultrasound with a 32 × 32 matrix array. Electromechanical activation maps are first presented in a normal and infarcted cardiac model in silico and in canine heart during pacing and re-entrant ventricular tachycardia in vivo. Noninvasive 4D electromechanical activation mapping in a healthy volunteer and a heart failure patient are also determined. The technique described herein allows for direct, simultaneous and noninvasive visualization of electromechanical activation in 3D, which provides complementary information on myocardial viability and/or abnormality to clinical imaging.

ECG Imaging to Detect the Site of Ventricular Ischemia Using Torso Electrodes: A Computational Study

Electrocardiography provides some information useful for ischemic diagnosis. However, more recently there has been substantial growth in the area of ECG imaging, which by solving the inverse problem of electrocardiography aims to produce high-resolution mapping of the electrical and magnetic dynamics of the heart. Most inverse studies use the full resolution of the body surface potential (BSP) to reconstruct the epicardial potentials, however using a limited number of torso electrodes to interpolate the BSP is more clinically relevant and has an important effect on the reconstruction which must be quantified. A circular ischemic lesion on the right ventricle lateral wall 27 mm in radius is reconstructed using three Tikhonov methods along with 6 different electrode configurations ranging from 32 leads to 1,024 leads. The 2nd order Tikhonov solution performed the most accurately (~80% lesion identified) followed by the 1st (~50% lesion identified) and then the 0 order Tikhonov solution performed the worst with a maximum of ~30% lesion identified regardless of how many leads were used. With an increasing number of leads the solution produces less error, and the error becomes more localised around the lesion for all three regularisation methods. In noisy conditions, the relative performance gap of the 1st and 2nd order Tikhonov solutions was reduced, and determining an accurate regularisation parameter became relatively more difficult. Lesions located on the left ventricle walls were also able to be identified but comparatively to the right ventricle lateral wall performed marginally worse with lesions located on the interventricular septum being able to be indicated by the reconstructions but not successfully identified against the error. The quality of reconstruction was found to decrease as the lesion radius decreased, with a lesion radius of <20 mm becoming difficult to correctly identify against the error even when using >512 torso electrodes.

Fully coupled fluid-electro-mechanical model of the human heart for supercomputers

In this work, we present a fully coupled fluid-electro-mechanical model of a 50th percentile human heart. The model is implemented on Alya, the BSC multi-physics parallel code, capable of running efficiently in supercomputers. Blood in the cardiac cavities is modeled by the incompressible Navier-Stokes equations and an arbitrary Lagrangian-Eulerian (ALE) scheme. Electrophysiology is modeled with a monodomain scheme and the O'Hara-Rudy cell model. Solid mechanics is modeled with a total Lagrangian formulation for discrete strains using the Holzapfel-Ogden cardiac tissue material model. The three problems are simultaneously and bidirectionally coupled through an electromechanical feedback and a fluid-structure interaction scheme. In this paper, we present the scheme in detail and propose it as a computational cardiac workbench.

Computational Prediction of the Combined Effect of CRT and LVAD on Cardiac Electromechanical Delay in LBBB and RBBB

Two case reports showed that the combination of CRT and LVAD benefits the end-stage heart failure patients with prolonged QRS interval significantly. In one of the reports, the patient had the LVAD removed due to the recovery of the heart function. However, the quantification of the combined devices has yet to be conducted. This study aimed at computationally predicting the effects of CRT-only or combined with LVAD on electromechanical behaviour in the failing ventricle with left bundle branch blocked (LBBB) and right bundle branch blocked (RBBB) conditions. The subjects are normal sinus rhythm, LBBB, RBBB, LBBB with CRT-only, RBBB with CRT-only, LBBB with CRT + LVAD, and RBBB with CRT + LVAD. The results showed that the CRT-only shortened the total electrical activation time (EAT) in the LBBB and RBBB conditions by 20.2% and 17.1%, respectively. The CRT-only reduced the total mechanical activation time (MAT) and electromechanical delay (EMD) of the ventricle under LBBB by 21.3% and 10.1%, respectively. Furthermore, the CRT-only reduced the contractile adenosine triphosphate (ATP) consumption by 5%, increased left ventricular (LV) pressure by 6%, and enhanced cardiac output (CO) by 0.2 L/min under LBBB condition. However, CRT-only barely affects the ventricle under RBBB condition. Under the LBBB condition, CRT + LVAD increased LV pressure and CO by 10.5% and by 0.9 L/min, respectively. CRT + LVAD reduced ATP consumption by 15%, shortened the MAT by 23.4%, and shortened the EMD by 15.2%. In conclusion, we computationally predicted and quantified that the CRT + LVAD implementation is superior to CRT-only implementation particularly in HF with LBBB condition.

Gaussian Process Regressions for Inverse Problems and Parameter Searches in Models of Ventricular Mechanics

Patient specific models of ventricular mechanics require the optimization of their many parameters under the uncertainties associated with imaging of cardiac function. We present a strategy to reduce the complexity of parametric searches for 3-D FE models of left ventricular contraction. The study employs automatic image segmentation and analysis of an image database to gain geometric features for several classes of patients. Statistical distributions of geometric parameters are then used to design parametric studies investigating the effects of: (1) passive material properties during ventricular filling, and (2) infarct geometry on ventricular contraction in patients after a heart attack. Gaussian Process regression is used in both cases to build statistical models trained on the results of biophysical FEM simulations. The first statistical model estimates unloaded configurations based on either the intraventricular pressure or the end-diastolic fiber strain. The technique provides an alternative to the standard fixed-point iteration algorithm, which is more computationally expensive when used to unload more than 10 ventricles. The second statistical model captures the effects of varying infarct geometries on cardiac output. For training, we designed high resolution models of non-transmural infarcts including refinements of the border zone around the lesion. This study is a first effort in developing a platform combining HPC models and machine learning to investigate cardiac function in heart failure patients with the goal of assisting clinical diagnostics.

Effects of gadolinium on cardiac mechanosensitivity in whole isolated swine hearts

Mechanical stimulation can elicit electrical activation of the heart. This mechanosensitivity can start life-threatening arrhythmias (commotio cordis) or terminate them (precordial thump). Mechanosensitivity may also be involved in arrhythmogenesis in other settings. Stretch-activated ion channels (SACs) are thought to be important in mechanosensitivity and a number of agents that block them have been identified. Such agents could potentially be used as tools in experimental investigation of mechanosensitivity. However, studies using them in intact-heart preparations have yielded inconsistent results. In the present study, we used isolated, perfused hearts from 25-35 kg pigs and a computer-controlled device that repeatably delivered focal mechanical stimuli. The concentration-dependent ability of the SAC blocker gadolinium to suppress mechanical activation was assessed by the success rate of mechanical stimulation and by the delay between successful mechanical stimulation and electrical activation. In six hearts, perfusate was recirculated. In an additional six hearts, perfusate was not recirculated to prevent gadolinium from forming complexes with metabolic waste and possibly precipitating. Gadolinium did not suppress mechanically-induced activation. Although gadolinium has been shown to be an effective SAC blocker in isolated cells, using it to probe the role of mechanical stimulation in whole heart preparations should be done with great caution.

Mechano-electric feedback effects in a three-dimensional (3D) model of the contracting cardiac ventricle

Mechano-electric feedback affects the electrophysiological and mechanical function of the heart and the cellular, tissue, and organ properties. To determine the main factors that contribute to this effect, this study investigated the changes in the action potential characteristics of the ventricle during contraction. A model of stretch-activated channels was incorporated into a three-dimensional multiscale model of the contracting ventricle to assess the effect of different preload lengths on the electrophysiological behavior. The model describes the initiation and propagation of the electrical impulse, as well as the passive (stretch) and active (contraction) changes in the cardiac mechanics. Simulations were performed to quantify the relationship between the cellular activation and recovery patterns as well as the action potential durations at different preload lengths in normal and heart failure pathological conditions. The simulation results showed that heart failure significantly affected the excitation propagation parameters compared to normal condition. The results showed that the mechano-electrical feedback effects appear to be most important in failing hearts with low ejection fraction.

Effects of mechanical feedback on the stability of cardiac scroll waves: A bidomain electro-mechanical simulation study

In this work, we investigate the influence of cardiac tissue deformation on re-entrant wave dynamics. We have developed a 3D strongly coupled electro-mechanical Bidomain model posed on an ideal monoventricular geometry, including fiber direction anisotropy and stretch-activated currents (SACs). The cardiac mechanical deformation influences the bioelectrical activity with two main mechanical feedback: (a) the geometric feedback (GEF) due to the presence of the deformation gradient in the diffusion coefficients and in a convective term depending on the deformation rate and (b) the mechano-electric feedback (MEF) due to SACs. Here, we investigate the relative contribution of these two factors with respect to scroll wave stability. We extend the previous works [Keldermann et al., Am. J. Physiol. Heart Circ. Physiol. 299, H134-H143 (2010) and Hu et al., PLoS One 8(4), e60287 (2013)] that were based on the Monodomain model and a simple non-selective linear SAC, while here we consider the full Bidomain model and both selective and non-selective components of SACs. Our simulation results show that the stability of cardiac scroll waves is influenced by MEF, which in case of low reversal potential of non-selective SACs might be responsible for the onset of ventricular fibrillation; GEF increases the scroll wave meandering but does not determine the scroll wave stability.

The importance of mechano-electrical feedback and inertia in cardiac electromechanics

In the past years, a number cardiac electromechanics models have been developed to better understand the excitation-contraction behavior of the heart. However, there is no agreement on whether inertial forces play a role in this system. In this study, we assess the influence of mass in electromechanical simulations, using a fully coupled finite element model. We include the effect of mechano-electrical feedback via stretch activated currents. We compare five different models: electrophysiology, electromechanics, electromechanics with mechano-electrical feedback, electromechanics with mass, and electromechanics with mass and mechano-electrical feedback. We simulate normal conduction to study conduction velocity and spiral waves to study fibrillation. During normal conduction, mass in conjunction with mechano-electrical feedback increased the conduction velocity by 8.12% in comparison to the plain electrophysiology case. During the generation of a spiral wave, mass and mechano-electrical feedback generated secondary wavefronts, which were not present in any other model. These secondary wavefronts were initiated in tensile stretch regions that induced electrical currents. We expect that this study will help the research community to better understand the importance of mechanoelectrical feedback and inertia in cardiac electromechanics.

Early regional wall distension is strongly associated with vulnerability to ventricular fibrillation but not arrhythmia triggers following coronary occlusion in vivo

Wall stress may favor ischemic ventricular arrhythmias, yet its association with ventricular fibrillation (VF) or ventricular ectopy has been inconsistent among studies and its potential arrhythmogenicity across the cardiac cycle is unclear. In 91 open-chest pigs undergoing 40-50 min left anterior descending artery occlusion, we assessed the association between diastolic or systolic distension of the ischemic area and the incidence of ventricular premature beats (VPBs) and VF. End-diastolic segment length (EDL) and systolic bulging ([maximum systolic length-EDL] × 100/EDL) were measured by ultrasonic crystals. Fifteen minutes after occlusion, EDL increased to 112.7 ± 5.6% of baseline (P < 0.001) and systolic bulging averaged 3.4 ± 2.2%. Median VPB number was 52 (IQR, 16-110), 2 (0-7) in phase Ia and 49 (13-94) in phase Ib. VF occurred in 26 animals (28.6%), the first episode appearing 24 ± 6 min after occlusion. EDL increase was associated with subsequent VF (115.9 ± 5.7 and 111.4 ± 5.1% in animals with and without VF, P < 0.001) and with the number of VF episodes (P = 0.001) but not with VPB number, overall (r = 0.028, P = 0.801) or in phases Ia or Ib. Systolic bulging was related neither to VF occurrence (3.2 ± 2.2 and 3.5 ± 2.2%, respectively, P = 0.561) nor to VBP number (r = 0.095, P = 0.397). EDL increase predicted VF after adjusting for ischemic area size and K⁺ levels (odds ratio for 1% increase: 1.17, 95%CI 1.06-1.29, P = 0.001). Thus, diastolic regional ventricular distension predicts VF occurrence after coronary occlusion whereas neither diastolic nor systolic distension is associated with ventricular ectopy, which suggests that distension favors VF by acting on the arrhythmic substrate but not on arrhythmia triggers.

Mechanisms Underlying the Emergence of Post-acidosis Arrhythmia at the Tissue Level: A Theoretical Study

Acidosis has complex electrophysiological effects, which are associated with a high recurrence of ventricular arrhythmias. Through multi-scale cardiac computer modeling, this study investigated the mechanisms underlying the emergence of post-acidosis arrhythmia at the tissue level. In simulations, ten Tusscher-Panfilov ventricular model was modified to incorporate various data on acidosis-induced alterations of cellular electrophysiology and intercellular electrical coupling. The single cell models were incorporated into multicellular one-dimensional (1D) fiber and 2D sheet tissue models. Electrophysiological effects were quantified as changes of action potential profile, sink-source interactions of fiber tissue, and the vulnerability of tissue to the genesis of unidirectional conduction that led to initiation of re-entry. It was shown that acidosis-induced sarcoplasmic reticulum (SR) calcium load contributed to delayed afterdepolarizations (DADs) in single cells. These DADs may be synchronized to overcome the source-sink mismatch arising from intercellular electrotonic coupling, and produce a premature ventricular complex (PVC) at the tissue level. The PVC conduction can be unidirectionally blocked in the transmural ventricular wall with altered electrical heterogeneity, resulting in the genesis of re-entry. In conclusion, altered source-sink interactions and electrical heterogeneity due to acidosis-induced cellular electrophysiological alterations may increase susceptibility to post-acidosis ventricular arrhythmias.

The role of infarct transmural extent in infarct extension: A computational study

Infarct extension, a process involving progressive extension of the infarct zone (IZ) into the normally perfused border zone (BZ), leads to continuous degradation of the myocardial function and adverse remodelling. Despite carrying a high risk of mortality, detailed understanding of the mechanisms leading to BZ hypoxia and infarct extension remains unexplored. In the present study, we developed a 3D truncated ellipsoidal left ventricular model incorporating realistic electromechanical properties and fibre orientation to examine the mechanical interaction among the remote, infarct and BZs in the presence of varying infarct transmural extent (TME). Localized highly abnormal systolic fibre stress was observed at the BZ, owing to the simultaneous presence of moderately increased stiffness and fibre strain at this region, caused by the mechanical tethering effect imposed by the overstretched IZ. Our simulations also demonstrated the greatest tethering effect and stress in BZ regions with fibre direction tangential to the BZ-remote zone boundary. This can be explained by the lower stiffness in the cross-fibre direction, which gave rise to a greater stretching of the IZ in this direction. The average fibre strain of the IZ, as well as the maximum stress in the sub-endocardial layer, increased steeply from 10% to 50% infarct TME, and slower thereafter. Based on our stress-strain loop analysis, we found impairment in the myocardial energy efficiency and elevated energy expenditure with increasing infarct TME, which we believe to place the BZ at further risk of hypoxia. Copyright © 2016 John Wiley & Sons, Ltd.

From Incidental, Mechanically-Induced Arrhythmias to Reflex-Defined Arrhythmogenicity: On The Track of The Ternary Reflex System Resemblance to The "Infancy" of New Era or Rediscovery

The underlying pathophysiology of supraventricular and ventricular arrhythmias remains a matter of intense investigation. Though evolving, the contemporary explanations do not encompass all aspects of arrhythmogenicity. An improved understanding of arrhythmia substrate is needed to augment therapeutic capabilities. Our observation and literature sources demonstrate relatively high incidence of transitory arrhythmias which are non-intentionally generated by the endocardial lead/catheter manipulation. These findings are interesting and potentially may crystallize the reflex-dependent proarrhythmic cardiac activity. Herein we suggest the "reflexogenic arrhythmogenicity" concept extending an overall spectrum of known hypotheses. Cardiovascular reflex action can be categorized into three-tiered levels - intra-cellular, inter-cellular and inter-organic. The first two levels of the triplicate system reside within the cardiac anatomical landmarks (in fact intramurally, intra-organically), however the third one implicates central (cerebral) activity which boomerangs back via centripetal and centrifugal connections. These levels likely compose synoptic ternary reflex set system which may be validated in future studies. To hypothesize, coordinated mutual reciprocity of reflex activity results in stabilization of heart rhythm in robust heart. Any stressful cardiac event may lead to the shift of the rhythm toward unfavorable clinical entity probably via the loss of the influence of dominant reflex. Overall, an interaction and likely intrinsic inter-tiered competition along with possible interplay between physiological and pathological reflexes may be treated as contributing factors for the inception and maintaining of arrhythmias and cardiac performance as well. These assumptions await further documentation. If such a tenet were recognized, the changes in the clinical approach to arrhythmia management might be anticipated, preferably by selective reflex suppression or activation strategy.

Comparison of Electric- and Magnetic-Cardiograms Produced by Myocardial Ischemia in Models of the Human Ventricle and Torso

Myocardial ventricular ischemia arises from a lack of blood supply to the heart, which may cause abnormal repolarization and excitation wave conduction patterns in the tissue, leading to cardiac arrhythmias and even sudden death. Current diagnosis of cardiac ischemia by the 12-lead electrocardiogram (ECG) has limitations as they are insensitive in many cases and may show unnoticeable differences to normal patterns. As the magnetic field provides extra information on cardiac excitation and is more sensitive to tangential currents to the surface of the chest, whereas the electric field is more sensitive to flux currents, it has been hypothesized that the magnetocardiogram (MCG) may provide a complementary method to the ECG in ischemic diagnosis. However, it is unclear yet about the differences in sensitivity regions of body surface ECG and MCG signals to ischemic conditions. The aim of this study was to investigate such differences by using 12-, 36- ECG and 36-MCG computed from multi-scale biophysically detailed computational models of the human ventricles and torso in both control and ischemic conditions. It was shown that ischemia produced changes in the ECG and MCG signals in the QRS complex, T-wave and ST-segment, with greater relative differences seen in the 36-lead ECG and MCG as compared to the 12-leads ECG (34% and 37% vs 26%, respectively). The 36-lead ECG showed more averaged sensitivity than the MCG in the change of T-wave due to ischemia (37% vs 32%, respectively), whereas the MCG showed greater sensitivity than the ECG in the change of the ST-segment (50% vs 40%, respectively). In addition, both MCG and ECG showed regional-dependent changes to ischemia, but with MCG showing a stronger correlation between ischemic region in the heart. In conclusion, MCG shows more sensitivity than ECG in response to ischemia, which may provide an alternative method for the diagnosis of ischemia.

Association Between Left Atrial Compression And Atrial Fibrillation: A Case Presentation And A Short Review Of Literature

This case report describes a patient who developed palpitations and chest pain and was found to be in atrial fibrillation, which was likely due to the presence of an extra-cardiac mass. This was compressing the left atrium. The mass was related to small cell carcinoma, which decreased significantly in size after chemotherapy. Resolution of the atrial fibrillation correlated temporally with reduction in the size of the mass and alleviation of the left atrial compression.

GsMTx4-D is a cardioprotectant against myocardial infarction during ischemia and reperfusion

GsMTx4 is a selective inhibitor of cationic mechanosensitive ion channels (MSCs) and has helped establish the role of MSCs in cardiac physiology. Inhomogeneous local mechanical stresses due to hypercontracture and swelling during ischemic reperfusion injury (IRI) likely induce elevated MSC activity that can contribute to cation imbalance. The aim of this study was to determine if the D enantiomer of GsMTx4 can act as a cardioprotectant in a mouse IRI model. Ischemia and reperfusion involved ligating a coronary artery followed by release of the ligature. GsMTx4-D was tested by either acute intravenous injection during the ischemic event or by two day pretreatment by intraperitoneal injection, both methods achieving similar results. Based on pharmacokinetic studies, GsMTx4-D dosage was set to achieve expected plasma concentrations between 50 and 5000nM and heart tissue concentrations between 1 and 200nM by intravenous injection. Relative to vehicle injected animals, GsMTx4-D reduced infarct area by ~40% for acute and pretreated animals for both 20 and 45min ischemic challenges. Many indicators of cardiac output were indistinguishable from sham-treated control hearts after GsMTx4-D treatment showing improvement at both 4 and 48h post ischemia, and premature ventricular beats immediately following reperfusion were also significantly reduced. To determine if GsMTx4-D cardioprotection could act directly at the level of cardiomyocytes, we tested its effects in vitro on indicators of IRI damage like cation influx and activation of inflammatory kinases in isolated myocytes cultured under hypoxic conditions. Hypoxia challenged cardiomyocytes treated with 10μM GsMTx4-D showed improved contractility and near normal contraction-related Ca(2+) influx. GsMTx4-D inhibited indicators of ischemic damage such as the apoptotic signaling system JNK/c-Jun, but also inhibited the energy response signaling system Akt kinase. We conclude that GsMTx4-D is a potent cardioprotectant in vivo that may act directly on cardiomyocytes and potentially be useful in multidrug strategies to treat IRI.

Rabbit models of cardiac mechano-electric and mechano-mechanical coupling

Cardiac auto-regulation involves integrated regulatory loops linking electrics and mechanics in the heart. Whereas mechanical activity is usually seen as 'the endpoint' of cardiac auto-regulation, it is important to appreciate that the heart would not function without feed-back from the mechanical environment to cardiac electrical (mechano-electric coupling, MEC) and mechanical (mechano-mechanical coupling, MMC) activity. MEC and MMC contribute to beat-by-beat adaption of cardiac output to physiological demand, and they are involved in various pathological settings, potentially aggravating cardiac dysfunction. Experimental and computational studies using rabbit as a model species have been integral to the development of our current understanding of MEC and MMC. In this paper we review this work, focusing on physiological and pathological implications for cardiac function.

Cardiac Mechano-Gated Ion Channels and Arrhythmias

Mechanical forces will have been omnipresent since the origin of life, and living organisms have evolved mechanisms to sense, interpret, and respond to mechanical stimuli. The cardiovascular system in general, and the heart in particular, is exposed to constantly changing mechanical signals, including stretch, compression, bending, and shear. The heart adjusts its performance to the mechanical environment, modifying electrical, mechanical, metabolic, and structural properties over a range of time scales. Many of the underlying regulatory processes are encoded intracardially and are, thus, maintained even in heart transplant recipients. Although mechanosensitivity of heart rhythm has been described in the medical literature for over a century, its molecular mechanisms are incompletely understood. Thanks to modern biophysical and molecular technologies, the roles of mechanical forces in cardiac biology are being explored in more detail, and detailed mechanisms of mechanotransduction have started to emerge. Mechano-gated ion channels are cardiac mechanoreceptors. They give rise to mechano-electric feedback, thought to contribute to normal function, disease development, and, potentially, therapeutic interventions. In this review, we focus on acute mechanical effects on cardiac electrophysiology, explore molecular candidates underlying observed responses, and discuss their pharmaceutical regulation. From this, we identify open research questions and highlight emerging technologies that may help in addressing them.

Computational rabbit models to investigate the initiation, perpetuation, and termination of ventricular arrhythmia

Current understanding of cardiac electrophysiology has been greatly aided by computational work performed using rabbit ventricular models. This article reviews the contributions of multiscale models of rabbit ventricles in understanding cardiac arrhythmia mechanisms. This review will provide an overview of multiscale modeling of the rabbit ventricles. It will then highlight works that provide insights into the role of the conduction system, complex geometric structures, and heterogeneous cellular electrophysiology in diseased and healthy rabbit hearts to the initiation and maintenance of ventricular arrhythmia. Finally, it will provide an overview on the contributions of rabbit ventricular modeling on understanding the mechanisms underlying shock-induced defibrillation.

Rabbit-specific computational modelling of ventricular cell electrophysiology: Using populations of models to explore variability in the response to ischemia

Computational modelling, combined with experimental investigations, is a powerful method for investigating complex cardiac electrophysiological behaviour. The use of rabbit-specific models, due to the similarities of cardiac electrophysiology in this species with human, is especially prevalent. In this paper, we first briefly review rabbit-specific computational modelling of ventricular cell electrophysiology, multi-cellular simulations including cellular heterogeneity, and acute ischemia. This mini-review is followed by an original computational investigation of variability in the electrophysiological response of two experimentally-calibrated populations of rabbit-specific ventricular myocyte action potential models to acute ischemia. We performed a systematic exploration of the response of the model populations to varying degrees of ischemia and individual ischemic parameters, to investigate their individual and combined effects on action potential duration and refractoriness. This revealed complex interactions between model population variability and ischemic factors, which combined to enhance variability during ischemia. This represents an important step towards an improved understanding of the role that physiological variability may play in electrophysiological alterations during acute ischemia.

Anatomically accurate high resolution modeling of human whole heart electromechanics: A strongly scalable algebraic multigrid solver method for nonlinear deformation

Electromechanical (EM) models of the heart have been used successfully to study fundamental mechanisms underlying a heart beat in health and disease. However, in all modeling studies reported so far numerous simplifications were made in terms of representing biophysical details of cellular function and its heterogeneity, gross anatomy and tissue microstructure, as well as the bidirectional coupling between electrophysiology (EP) and tissue distension. One limiting factor is the employed spatial discretization methods which are not sufficiently flexible to accommodate complex geometries or resolve heterogeneities, but, even more importantly, the limited efficiency of the prevailing solver techniques which are not sufficiently scalable to deal with the incurring increase in degrees of freedom (DOF) when modeling cardiac electromechanics at high spatio-temporal resolution. This study reports on the development of a novel methodology for solving the nonlinear equation of finite elasticity using human whole organ models of cardiac electromechanics, discretized at a high para-cellular resolution. Three patient-specific, anatomically accurate, whole heart EM models were reconstructed from magnetic resonance (MR) scans at resolutions of 220 μm, 440 μm and 880 μm, yielding meshes of approximately 184.6, 24.4 and 3.7 million tetrahedral elements and 95.9, 13.2 and 2.1 million displacement DOF, respectively. The same mesh was used for discretizing the governing equations of both electrophysiology (EP) and nonlinear elasticity. A novel algebraic multigrid (AMG) preconditioner for an iterative Krylov solver was developed to deal with the resulting computational load. The AMG preconditioner was designed under the primary objective of achieving favorable strong scaling characteristics for both setup and solution runtimes, as this is key for exploiting current high performance computing hardware. Benchmark results using the 220 μm, 440 μm and 880 μm meshes demonstrate efficient scaling up to 1024, 4096 and 8192 compute cores which allowed the simulation of a single heart beat in 44.3, 87.8 and 235.3 minutes, respectively. The efficiency of the method allows fast simulation cycles without compromising anatomical or biophysical detail.

Mechano-electric feedback in one-dimensional model of myocardium

We utilized our earlier developed 1D mathematical model of the heart muscle strand to study contribution of the bilateral interactions between excitation and contraction on the cellular and tissue levels to the local and global myocardium function. Numerical experiments on the model showed that an initially uniform strand, formed on the inherently identical cells, became functionally heterogeneous due to the asynchronous excitation via the electrical wave spread. Mechanical interactions between the cells and the mechano-electric feedback beat-to-beat affect the functional characteristics of coupled cardiomyocytes further, adjusting their electrical and mechanical heterogeneity to the activation timing. Model simulations showed that functional heterogeneity increases with an enlarged spatial extension of the myocardial strand (in terms of the longer slack length not a higher stretch of the strand), demonstrating a special role of the heart size in its function. Model analysis suggests that cooperative mechanisms of myofilament calcium activation contribute essentially to the generation of cellular functional heterogeneity in contracting cardiac tissue.

Effects of the Selective Stretch-Activated Channel Blocker GsMtx4 on Stretch-Induced Changes in Refractoriness in Isolated Rat Hearts and on Ventricular Premature Beats and Arrhythmias after Coronary Occlusion in Swine

Mechanical factors may contribute to ischemic ventricular arrhythmias. GsMtx4 peptide, a selective stretch-activated channel blocker, inhibits stretch-induced atrial arrhythmias. We aimed to assess whether GsMtx4 protects against ventricular ectopy and arrhythmias following coronary occlusion in swine. First, the effects of 170-nM GsMtx4 on the changes in the effective refractory period (ERP) induced by left ventricular (LV) dilatation were assessed in 8 isolated rat hearts. Then, 44 anesthetized, open-chest pigs subjected to 50-min left anterior descending artery occlusion and 2-h reperfusion were blindly allocated to GsMtx4 (57 μg/kg iv. bolus and 3.8 μg/kg/min infusion, calculated to attain the above concentration in plasma) or saline, starting 5-min before occlusion and continuing until after reflow. In rat hearts, LV distension induced progressive reductions in ERP (35±2, 32±2, and 29±2 ms at 0, 20, and 40 mmHg of LV end-diastolic pressure, respectively, P<0.001) that were prevented by GsMTx4 (33±2, 33±2, and 32±2 ms, respectively, P=0.002 for the interaction with LV end-diastolic pressure). Pigs receiving GsMtx4 had similar number of ventricular premature beats during the ischemic period as control pigs (110±28 vs. 103±21, respectively, P=0.842). There were not significant differences among treated and untreated animals in the incidence of ventricular fibrillation (13.6 vs. 22.7%, respectively, P=0.696) or tachycardia (36.4 vs. 50.0%, P=0.361) or in the number of ventricular tachycardia episodes during the occlusion period (1.8±0.7 vs. 5.5±2.6, P=0.323). Thus, GsMtx4 administered under these conditions does not suppress ventricular ectopy following coronary occlusion in swine. Whether it might protect against malignant arrhythmias should be tested in studies powered for these outcomes.

Variations in local calcium signaling in adjacent cardiac myocytes of the intact mouse heart detected with two-dimensional confocal microscopy

Dyssynchronous local Ca release within individual cardiac myocytes has been linked to cellular contractile dysfunction. Differences in Ca kinetics in adjacent cells may also provide a substrate for inefficient contraction and arrhythmias. In a new approach we quantify variation in local Ca transients between adjacent myocytes in the whole heart. Langendorff-perfused mouse hearts were loaded with Fluo-8 AM to detect Ca and Di-4-ANEPPS to visualize cell membranes. A spinning disc confocal microscope with a fast camera allowed us to record Ca signals within an area of 465 μm by 315 μm with an acquisition speed of 55 fps. Images from multiple transients recorded at steady state were registered to their time point in the cardiac cycle to restore averaged local Ca transients with a higher temporal resolution. Local Ca transients within and between adjacent myocytes were compared with regard to amplitude, time to peak and decay at steady state stimulation (250 ms cycle length). Image registration from multiple sequential Ca transients allowed reconstruction of high temporal resolution (2.4 ± 1.3 ms) local CaT in 2D image sets (N = 4 hearts, n = 8 regions). During steady state stimulation, spatial Ca gradients were homogeneous within cells in both directions and independent of distance between measured points. Variation in CaT amplitudes was similar across the short and the long side of neighboring cells. Variations in TAU and TTP were similar in both directions. Isoproterenol enhanced the CaT but not the overall pattern of spatial heterogeneities. Here we detected and analyzed local Ca signals in intact mouse hearts with high temporal and spatial resolution, taking into account 2D arrangement of the cells. We observed significant differences in the variation of CaT amplitude along the long and short axis of cardiac myocytes. Variations of Ca signals between neighboring cells may contribute to the substrate of cardiac remodeling.

Exploring susceptibility to atrial and ventricular arrhythmias resulting from remodeling of the passive electrical properties in the heart: a simulation approach

Under diseased conditions, remodeling of the cardiac tissue properties (“passive properties”) takes place; these are aspects of electrophysiological behavior that are not associated with active ion transport across cell membranes. Remodeling of the passive electrophysiological properties most often results from structural remodeling, such as gap junction down-regulation and lateralization, fibrotic growth infiltrating the myocardium, or the development of an infarct scar. Such structural remodeling renders atrial or ventricular tissue as a major substrate for arrhythmias. The current review focuses on these aspects of cardiac arrhythmogenesis. Due to the inherent complexity of cardiac arrhythmias, computer simulations have provided means to elucidate interactions pertinent to this spatial scale. Here we review the current state-of-the-art in modeling atrial and ventricular arrhythmogenesis as arising from the disease-induced changes in the passive tissue properties, as well as the contributions these modeling studies have made to our understanding of the mechanisms of arrhythmias in the heart. Because of the rapid advance of structural imaging methodologies in cardiac electrophysiology, we chose to present studies that have used such imaging methodologies to construct geometrically realistic models of cardiac tissue, or the organ itself, where the regional remodeling properties of the myocardium can be represented in a realistic way. We emphasize how the acquired knowledge can be used to pave the way for clinical applications of cardiac organ modeling under the conditions of structural remodeling.

Multiscale computational analysis of the bioelectric consequences of myocardial ischaemia and infarction

Ischaemic heart disease is considered as the single most frequent cause of death, provoking more than 7 000 000 deaths every year worldwide. A high percentage of patients experience sudden cardiac death, caused in most cases by tachyarrhythmic mechanisms associated to myocardial ischaemia and infarction. These diseases are difficult to study using solely experimental means due to their complex dynamics and unstable nature. In the past decades, integrative computational simulation techniques have become a powerful tool to complement experimental and clinical research when trying to elucidate the intimate mechanisms of ischaemic electrophysiological processes and to aid the clinician in the improvement and optimization of therapeutic procedures. The purpose of this paper is to briefly review some of the multiscale computational models of myocardial ischaemia and infarction developed in the past 20 years, ranging from the cellular level to whole-heart simulations.

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.

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.

Advances in modeling ventricular arrhythmias: from mechanisms to the clinic

Modern cardiovascular research has increasingly recognized that heart models and simulation can help interpret an array of experimental data and dissect important mechanisms and interrelationships, with developments rooted in the iterative interaction between modeling and experimentation. This article reviews the progress made in simulating cardiac electrical behavior at the level of the organ and, specifically, in the development of models of ventricular arrhythmias and fibrillation, as well as their termination (defibrillation). The ability to construct multiscale models of ventricular arrhythmias, representing integrative behavior from the molecule to the entire organ, has enabled mechanistic inquiry into the dynamics of ventricular arrhythmias in the diseased myocardium, in understanding drug-induced proarrhythmia, and in the development of new modalities for defibrillation, to name a few. In this article, we also review the initial use of ventricular models of arrhythmia in personalized diagnosis, treatment planning, and prevention of sudden cardiac death. Implementing individualized cardiac simulations at the patient bedside is poised to become one of the most thrilling examples of computational science and engineering approaches in translational medicine.

Efficient preloading of the ventricles by a properly timed atrial contraction underlies stroke work improvement in the acute response to cardiac resynchronization therapy

BACKGROUND

The acute response to cardiac resynchronization therapy (CRT) has been shown to be due to 3 mechanisms: resynchronization of ventricular contraction, efficient preloading of the ventricles by a properly timed atrial contraction, and mitral regurgitation reduction. However, the contribution of each of the 3 mechanisms to the acute response to CRT, specifically stroke work improvement, has not been quantified.

OBJECTIVE

To use a magnetic resonance image-based anatomically accurate 3-dimensional model of failing canine ventricular electromechanics to quantify the contribution of each of the 3 mechanisms to stroke work improvement and identify the predominant mechanisms.

METHODS

An MRI-based electromechanical model of the failing canine ventricles assembled previously by our group was further developed and modified. Three different protocols were used to dissect the contribution of each of the 3 mechanisms to stroke work improvement.

RESULTS

Resynchronization of ventricular contraction did not lead to a significant stroke work improvement. Efficient preloading of the ventricles by a properly timed atrial contraction was the predominant mechanism underlying stroke work improvement. Stroke work improvement peaked at an intermediate atrioventricular delay, as it allowed ventricular filling by atrial contraction to occur at a low diastolic left ventricular pressure but also provided adequate time for ventricular filling before ventricular contraction. Reduction of mitral regurgitation by CRT led to stroke work worsening instead of improvement.

CONCLUSION

Efficient preloading of the ventricles by a properly timed atrial contraction is responsible for a significant stroke work improvement in the acute CRT response.

Mechanistic insight into prolonged electromechanical delay in dyssynchronous heart failure: a computational study

In addition to the left bundle branch block type of electrical activation, there are further remodeling aspects associated with dyssynchronous heart failure (HF) that affect the electromechanical behavior of the heart. Among the most important are altered ventricular structure (both geometry and fiber/sheet orientation), abnormal Ca(2+) handling, slowed conduction, and reduced wall stiffness. In dyssynchronous HF, the electromechanical delay (EMD), the time interval between local myocyte depolarization and myofiber shortening onset, is prolonged. However, the contributions of the four major HF remodeling aspects in extending EMD in the dyssynchronous failing heart remain unknown. The goal of this study was to determine the individual and combined contributions of HF-induced remodeling aspects to EMD prolongation. We used MRI-based models of dyssynchronous nonfailing and HF canine electromechanics and constructed additional models in which varying combinations of the four remodeling aspects were represented. A left bundle branch block electrical activation sequence was simulated in all models. The simulation results revealed that deranged Ca(2+) handling is the primary culprit in extending EMD in dyssynchronous HF, with the other aspects of remodeling contributing insignificantly. Mechanistically, we found that abnormal Ca(2+) handling in dyssynchronous HF slows myofiber shortening velocity at the early-activated septum and depresses both myofiber shortening and stretch rate at the late-activated lateral wall. These changes in myofiber dynamics delay the onset of myofiber shortening, thus giving rise to prolonged EMD in dyssynchronous HF.

Effects of mechano-electric feedback on scroll wave stability in human ventricular fibrillation

Recruitment of stretch-activated channels, one of the mechanisms of mechano-electric feedback, has been shown to influence the stability of scroll waves, the waves that underlie reentrant arrhythmias. However, a comprehensive study to examine the effects of recruitment of stretch-activated channels with different reversal potentials and conductances on scroll wave stability has not been undertaken; the mechanisms by which stretch-activated channel opening alters scroll wave stability are also not well understood. The goals of this study were to test the hypothesis that recruitment of stretch-activated channels affects scroll wave stability differently depending on stretch-activated channel reversal potential and channel conductance, and to uncover the relevant mechanisms underlying the observed behaviors. We developed a strongly-coupled model of human ventricular electromechanics that incorporated human ventricular geometry and fiber and sheet orientation reconstructed from MR and diffusion tensor MR images. Since a wide variety of reversal potentials and channel conductances have been reported for stretch-activated channels, two reversal potentials, −60 mV and −10 mV, and a range of channel conductances (0 to 0.07 mS/µF) were implemented. Opening of stretch-activated channels with a reversal potential of −60 mV diminished scroll wave breakup for all values of conductances by flattening heterogeneously the action potential duration restitution curve. Opening of stretch-activated channels with a reversal potential of −10 mV inhibited partially scroll wave breakup at low conductance values (from 0.02 to 0.04 mS/µF) by flattening heterogeneously the conduction velocity restitution relation. For large conductance values (>0.05 mS/µF), recruitment of stretch-activated channels with a reversal potential of −10 mV did not reduce the likelihood of scroll wave breakup because Na channel inactivation in regions of large stretch led to conduction block, which counteracted the increased scroll wave stability due to an overall flatter conduction velocity restitution.

The effects of κ-opioid receptor on stretch-induced electrophysiological changes in infarcted rat hearts

INTRODUCTION

Kappa-opioid receptors (κ-OR) and mechanoelectric feedback seem to have common pathways that influence electrophysiological changes resulting from acute myocardial infarction (MI). This study aims to determine the effects of the κ-OR on stretch-induced electrophysiological changes after acute MI.

METHODS

Male Sprague-Dawley rats were randomly divided into 4 groups: sham operated, MI, U-50488H (a selective κ-OR agonist) -treated MI (MI+U-50488H) and nor-BNI (a selective κ-OR antagonist) -treated MI (MI+nor-BNI). After Langendorff perfusion to maintain stabilization, a transient stretch (5 seconds) was delivered early in diastole. Electrophysiological changes were recorded for 1 minute before and after stretch. Similarly, the 20%, 50% and 90% monophasic action potential duration (MAPD20, MAPD50 and MAPD90, respectively) and stretch-induced arrhythmias were recorded.

RESULTS

MAPD90 significantly increased in all 4 groups. MAPD90 in the MI and MI+nor-BNI groups increased significantly before stretch (P < 0.05) and after stretch (P < 0.01) but was reversed in the MI+U-50488H group (P > 0.05). MAPD90 in the MI group was increased compared with that of the MI+U-50488H group but decreased compared with that of the MI+ nor-BNI group after stretch (P < 0.01). The arrhythmia score in the MI and MI+nor-BNI groups was higher than that of the sham-operated group (P < 0.01), and the arrhythmia score in the MI+nor-BNI group was higher than that in MI group after stretch (P < 0.01). The arrhythmia score of the MI+U-50488H group was lower than that of MI group after stretch (P < 0.01).

CONCLUSIONS

The κ-OR could influence the stretch-induced electrophysiological changes and play an antiarrhythmic role in stretch-induced arrhythmias after acute MI.

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.