ATP-sensitive K+ channels and cellular K+ loss in hypoxic and ischaemic mammalian ventricle

Electro-metabolic signaling

Precise matching of energy substrate delivery to local metabolic needs is essential for the health and function of all tissues. Here, we outline a mechanistic framework for understanding this critical process, which we refer to as electro-metabolic signaling (EMS). All tissues exhibit changes in metabolism over varying spatiotemporal scales and have widely varying energetic needs and reserves. We propose that across tissues, common signatures of elevated metabolism or increases in energy substrate usage that exceed key local thresholds rapidly engage mechanisms that generate hyperpolarizing electrical signals in capillaries that then relax contractile elements throughout the vasculature to quickly adjust blood flow to meet changing needs. The attendant increase in energy substrate delivery serves to meet local metabolic requirements and thus avoids a mismatch in supply and demand and prevents metabolic stress. We discuss in detail key examples of EMS that our laboratories have discovered in the brain and the heart, and we outline potential further EMS mechanisms operating in tissues such as skeletal muscle, pancreas, and kidney. We suggest that the energy imbalance evoked by EMS uncoupling may be central to cellular dysfunction from which the hallmarks of aging and metabolic diseases emerge and may lead to generalized organ failure states-such as diverse flavors of heart failure and dementia. Understanding and manipulating EMS may be key to preventing or reversing these dysfunctions.

The mechanisms of potassium loss in acute myocardial ischemia: New insights from computational simulations

Acute myocardial ischemia induces hyperkalemia (accumulation of extracellular potassium), a major perpetrator of lethal reentrant ventricular arrhythmias. Despite considerable experimental efforts to explain this pathology in the last decades, the intimate mechanisms behind hyperkalemia remain partially unknown. In order to investigate these mechanisms, we developed a novel computational model of acute myocardial ischemia which couples a) an electrophysiologically detailed human cardiomyocyte model that incorporates modifications to account for ischemia-induced changes in transmembrane currents, with b) a model of cardiac tissue and extracellular K ⁺ transport. The resulting model is able to reproduce and explain the triphasic time course of extracellular K ⁺ concentration within the ischemic zone, with values of [ K + ] o close to 14 mmol/L in the central ischemic zone after 30 min. In addition, the formation of a [ K + ] o border zone of approximately 1.2 cm 15 min after the onset of ischemia is predicted by the model. Our results indicate that the primary rising phase of [ K + ] o is mainly due to the imbalance between K ⁺ efflux, that increases slightly, and K ⁺ influx, that follows a reduction of the NaK pump activity by more than 50%. The onset of the plateau phase is caused by the appearance of electrical alternans (a novel mechanism identified by the model), which cause an abrupt reduction in the K ⁺ efflux. After the plateau, the secondary rising phase of [ K + ] o is caused by a subsequent imbalance between the K ⁺ influx, which continues to decrease slowly, and the K ⁺ efflux, which remains almost constant. Further, the study shows that the modulation of these mechanisms by the electrotonic coupling is the main responsible for the formation of the ischemic border zone in tissue, with K ⁺ transport playing only a minor role. Finally, the results of the model indicate that the injury current established between the healthy and the altered tissue is not sufficient to depolarize non-ischemic cells within the healthy tissue.

Mechanisms of ventricular arrhythmias elicited by coexistence of multiple electrophysiological remodeling in ischemia: A simulation study

Myocardial ischemia, injury and infarction (MI) are the three stages of acute coronary syndrome (ACS). In the past two decades, a great number of studies focused on myocardial ischemia and MI individually, and showed that the occurrence of reentrant arrhythmias is often associated with myocardial ischemia or MI. However, arrhythmogenic mechanisms in the tissue with various degrees of remodeling in the ischemic heart have not been fully understood. In this study, biophysical detailed single-cell models of ischemia 1a, 1b, and MI were developed to mimic the electrophysiological remodeling at different stages of ACS. 2D tissue models with different distributions of ischemia and MI areas were constructed to investigate the mechanisms of the initiation of reentrant waves during the progression of ischemia. Simulation results in 2D tissues showed that the vulnerable windows (VWs) in simultaneous presence of multiple ischemic conditions were associated with the dynamics of wave propagation in the tissues with each single pathological condition. In the tissue with multiple pathological conditions, reentrant waves were mainly induced by two different mechanisms: one is the heterogeneity along the excitation wavefront, especially the abrupt variation in conduction velocity (CV) across the border of ischemia 1b and MI, and the other is the decreased safe factor (SF) for conduction at the edge of the tissue in MI region which is attributed to the increased excitation threshold of MI region. Finally, the reentrant wave was observed in a 3D model with a scar reconstructed from MRI images of a MI patient. These comprehensive findings provide novel insights for understanding the arrhythmic risk during the progression of myocardial ischemia and highlight the importance of the multiple pathological stages in designing medical therapies for arrhythmias in ischemia.

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.

Analysis of vulnerability to reentry in acute myocardial ischemia using a realistic human heart model

Electrophysiological alterations of the myocardium caused by acute ischemia constitute a pro-arrhythmic substrate for the generation of potentially lethal arrhythmias. Experimental evidence has shown that the main components of acute ischemia that induce these electrophysiological alterations are hyperkalemia, hypoxia (or anoxia in complete artery occlusion), and acidosis. However, the influence of each ischemic component on the likelihood of reentry is not completely established. Moreover, the role of the His-Purkinje system (HPS) in the initiation and maintenance of arrhythmias is not completely understood. In the present work, we investigate how the three components of ischemia affect the vulnerable window (VW) for reentry using computational simulations. In addition, we analyze the role of the HPS on arrhythmogenesis. A 3D biventricular/torso human model that includes a realistic geometry of the central and border ischemic zones with one of the most electrophysiologically detailed model of ischemia to date, as well as a realistic cardiac conduction system, were used to assess the VW for reentry. Four scenarios of ischemic severity corresponding to different minutes after coronary artery occlusion were simulated. Our results suggest that ischemic severity plays an important role in the generation of reentries. Indeed, this is the first 3D simulation study to show that ventricular arrhythmias could be generated under moderate ischemic conditions, but not in mild and severe ischemia. Moreover, our results show that anoxia is the ischemic component with the most significant effect on the width of the VW. Thus, a change in the level of anoxia from moderate to severe leads to a greater increment in the VW (40 ms), in comparison with the increment of 20 ms and 35 ms produced by the individual change in the level of hyperkalemia and acidosis, respectively. Finally, the HPS was a necessary element for the generation of approximately 17% of reentries obtained. The retrograde conduction from the myocardium to HPS in the ischemic region, conduction blocks in discrete sections of the HPS, and the degree of ischemia affecting Purkinje cells, are suggested as mechanisms that favor the generation of ventricular arrhythmias.

Protonic conductor: better understanding neural resting and action potential

With the employment of the transmembrane electrostatic proton localization theory with a new membrane potential equation, neural resting and action potential is now much better understood as the voltage contributed by the localized protons/cations at a neural liquid- membrane interface. Accordingly, the neural resting/action potential is essentially a protonic/cationic membrane capacitor behavior. It is now understood with a newly formulated action potential equation: when action potential is <0 (negative number), the localized protons/cations charge density at the liquid-membrane interface along the periplasmic side is >0 (positive number); when the action potential is >0, the concentration of the localized protons and localized nonproton cations is <0, indicating a "depolarization" state. The nonlinear curve of the localized protons/cations charge density in the real-time domain of an action potential spike appears as an inverse mirror image to the action potential. The newly formulated action potential equation provides biophysical insights for neuron electrophysiology, which may represent a complementary development to the classic Goldman-Hodgkin-Katz equation. With the use of the action potential equation, the biological significance of axon myelination is now also elucidated as to provide protonic insulation and prevent any ions both inside and outside of the neuron from interfering with the action potential signal, so that the action potential can quickly propagate along the axon with minimal (e.g., 40 times less) energy requirement.NEW & NOTEWORTHY The newly formulated action potential equation provides biophysical insights for neuron electrophysiology, which may represent a complementary development to the classic Goldman-Hodgkin-Katz equation. The nonlinear curve of the localized protons/cations charge density in the real-time domain of an action potential spike appears as an inverse mirror image to the action potential. The biological significance of axon myelination is now elucidated as to provide protonic insulation and prevent any ions from interfering with action potential signal.

Effects of β-adrenergic stimulation on fetal heart rate, heart rate variability, and T-wave elevation during brief umbilical cord occlusions in fetal sheep

Circulating catecholamines are critical for fetal adaptation to hypoxia by regulating fetal heart rate (FHR) and promoting myocardial contractility and peripheral vasoconstriction. They have been hypothesized to contribute to changes in FHR variability (FHRV) and T-wave morphology, clinical indexes of fetal well-being during labor. β-Adrenergic blockade with propranolol does not affect FHRV during labor-like hypoxemia and only attenuated the increase in T-wave height between the episodes of hypoxemia. To further investigate the potential role of catecholamines, we investigated whether pharmacological β-adrenergic stimulation could increase FHRV and T-wave elevation during intermittent labor-like hypoxemia. Nineteen chronically instrumented fetal sheep at 0.85 of gestation received isoprenaline hydrochloride (n = 7) or saline (control, n = 12), followed by three 1-min complete umbilical cord occlusions (UCOs) separated by 4-min reperfusion periods. Before the UCOs, infusion of isoprenaline increased FHR (P < 0.001), absolute-T/QRS ratio (P < 0.001), and one measure of FHRV [root-mean-square of successive RR interval differences (RMSSD), P < 0.05]. UCOs triggered deep FHR decelerations. During UCOs, isoprenaline was associated with increased FHR (P < 0.001) and absolute-T/QRS ratio (P < 0.05), but no effect on T/QRS ratio was observed when normalized to baseline before UCOs (normalized-T/QRS ratio). Between UCOs, isoprenaline increased FHR (P < 0.001) and absolute-T/QRS ratio (P < 0.05) but did not affect normalized-T/QRS ratio or any measures of FHRV. Arterial pressure was not affected by isoprenaline at any point. Our findings indicate that circulating catecholamines regulate FHR but not FHRV during labor-like hypoxemia and promote T-wave elevation between but not during intermittent fetal hypoxemia.

Circulating catecholamines partially regulate T-wave morphology but not heart rate variability during repeated umbilical cord occlusions in fetal sheep

Fetal heart rate (FHR) variability (FHRV) and ST segment morphology are potential clinical indices of fetal well-being during labor. β-Adrenergic stimulation by circulating catecholamines has been hypothesized to contribute to both FHRV and ST segment morphology during labor, but this has not been tested during brief repeated fetal hypoxemia that is characteristic of labor. Near-term fetal sheep (0.85 gestation) received propranolol (β-adrenergic blockade; n = 10) or saline (n = 7) 30 min before being exposed to three 2-min complete umbilical cord occlusions (UCOs) separated by 3-min reperfusions. T/QRS ratio was calculated throughout UCOs and reperfusion periods, and measures of FHRV (RMSSD, SDNN, and STV) were calculated between UCOs. During the baseline period, before the start of UCOs, propranolol was associated with reduced FHR, SDNN, and STV but did not affect RMSSD or T/QRS ratio. UCOs were associated with rapid FHR decelerations and increased T/QRS ratio; propranolol significantly reduced FHR during UCOs and was associated with a slower rise in T/QRS ratio during the first UCOs, without affecting the maximal rise or T/QRS ratio during the second and third UCO. Between UCOs propranolol reduced FHR and T/QRS ratio but did not affect any measure of FHRV. These data demonstrate that circulating catecholamines do not contribute to FHRV during labor-like hypoxemia. Furthermore, circulating catecholamines did not contribute to the major rise in T/QRS ratio during labor-like hypoxemia but may regulate T/QRS ratio between brief hypoxemia.

Pro-arrhythmic Effects of Hydrogen Sulfide in Healthy and Ischemic Cardiac Tissues: Insight From a Simulation Study

Hydrogen sulfide (H₂S), an ambient air pollutant, has been reported to increase cardiac events in patients with cardiovascular diseases, but the underlying mechanisms remain not elucidated. This study investigated the pro-arrhythmic effects of H₂S in healthy and ischemic conditions. Experimental data of H₂S effects on ionic channels (including the L-type Ca²⁺ channel and ATP-sensitive K⁺ channel) were incorporated into a virtual heart model to evaluate their integral action on cardiac arrhythmogenesis. It was shown that H₂S depressed cellular excitability, abbreviated action potential duration, and augmented tissue’s transmural dispersion of repolarization, resulting in an increase in tissue susceptibility to initiation and maintenance of reentry. The observed effects of H₂S on cardiac excitation are more remarkable in the ischemic condition than in the healthy condition. This study provides mechanistic insights into the pro-arrhythmic effects of air pollution (H₂S), especially in the case with extant ischemic conditions.

Modeling and simulation of cardiac electric activity in a human cardiac tissue with multiple ischemic zones

In this work, a human ventricular model (ten Tusscher and Panfilov model) coupled with the tissue level monodomain model is used to analyze the influence of multiple myocardial ischemia on the human cardiac tissue. The existence and uniqueness of the ischemic model comprising the monodomain model with a discontinuous ionic model for the human cardiac tissue is discussed. The coupled system of partial differential equation and ordinary differential equations are solved numerically using [Formula: see text] finite elements in space and Backward Euler finite difference scheme in time. The apriori finite element error estimate for the numerical scheme has been shown to be of [Formula: see text]. Essentially, we evaluate the impact of the increasing size of the ischemic region and the presence of the multiple ischemic regions having equal or different intensities on the neighboring healthy part of the cardiac tissue. We examine both the individual and the combined influence of two types of ischemia, Hyperkalemia (with the variation of the extracellular potassium ion concentration, [Formula: see text]) and Hypoxia (with the variation of intracellular Adenosine triphosphate (ATP) concentration via parameter [Formula: see text]) on the cardiac electrical activity of cardiac tissue. We observe that with the increase in the ischemic region size by a factor five times, there is an additional almost 10% drop in the action potential duration (APD) in the neighboring healthy regions. The combined effect of Hyperkalemia and Hypoxia brings an additional 12% drop in APD in the ischemic subregions and an additional 5% drop in APD in the neighboring healthy part of the cardic tissue in comparison to the only Hyperkalemic ischemia. When the Hyperkalemic and/or Hypoxic degeneracy of a ischemic zone is non-uniform then innercore degeneracy has greater influence on resting potential and APD of outercore of variable intensity ischemic zone than the other way. Also, increasing the number of ischemic subregions from 2 to 4 leads to a 4% drop in APD.

β-Adrenergic Inhibition Prevents Action Potential and Calcium Handling Changes during Regional Myocardial Ischemia

β-adrenergic receptor (β-AR) blockers may be administered during acute myocardial infarction (MI), as they reduce energy demand through negative chronotropic and inotropic effects and prevent ischemia-induced arrhythmogenesis. However, the direct effects of β-AR blockers on ventricular electrophysiology and intracellular Ca²⁺ handling during ischemia remain unknown. Using optical mapping of transmembrane potential (with RH237) and sarcoplasmic reticulum (SR) Ca²⁺ (with the low-affinity indicator Fluo-5N AM), the effects of 15 min of regional ischemia were assessed in isolated rabbit hearts (n = 19). The impact of β-AR inhibition on isolated hearts was assessed by pre-treatment with 100 nM propranolol (Prop) prior to ischemia (n = 7). To control for chronotropy and inotropy, hearts were continuously paced at 3.3 Hz and contraction was inhibited with 20 μM blebbistatin. Untreated ischemic hearts displayed prototypical shortening of action potential duration (APD₈₀) in the ischemic zone (IZ) compared to the non-ischemic zone (NI) at 10 and 15 min ischemia, whereas APD shortening was prevented with Prop. Untreated ischemic hearts also displayed significant changes in SR Ca²⁺ handling in the IZ, including prolongation of SR Ca²⁺ reuptake and SR Ca²⁺ alternans, which were prevented with Prop pre-treatment. At 5 min ischemia, Prop pre-treated hearts also showed larger SR Ca²⁺ release amplitude in the IZ compared to untreated hearts. These results suggest that even when controlling for chronotropic and inotropic effects, β-AR inhibition has a favorable effect during acute regional ischemia via direct effects on APD and Ca²⁺ handling.

KATP channel inhibition blunts electromechanical decline during hypoxia in left ventricular working rabbit hearts

KEY POINTS

Heart function is critically dependent upon the balance of energy production and utilization. Sarcolemmal ATP-sensitive potassium channels (K_ATP channels) in cardiac myocytes adjust contractile function to compensate for the level of available energy. Understanding the activation of K_ATP channels in working myocardium during high-stress situations is crucial to the treatment of cardiovascular disease, especially ischaemic heart disease. Using a new optical mapping approach, we measured action potentials from the surface of excised contracting rabbit hearts to assess when sarcolemmal K_ATP channels were activated during physiologically relevant workloads and during gradual reductions in myocardial oxygenation. We demonstrate that left ventricular pressure is closely linked to K_ATP channel activation and that K_ATP channel inhibition with a low concentration of tolbutamide prevents electromechanical decline when oxygen availability is reduced. As a result, K_ATP channel inhibition probably exacerbates a mismatch between energy demand and energy production when myocardial oxygenation is low.

ABSTRACT

Sarcolemmal ATP-sensitive potassium channel (K_ATP channel) activation in isolated cells is generally understood, although the relationship between myocardial oxygenation and K_ATP activation in excised working rabbit hearts remains unknown. We optically mapped action potentials (APs) in excised rabbit hearts to test the hypothesis that hypoxic changes would be more severe in left ventricular (LV) working hearts (LWHs) than Langendorff (LANG) perfused hearts. We further hypothesized that K_ATP inhibition would prevent those changes. Optical APs were mapped when measuring LV developed pressure (LVDP), coronary flow rate and oxygen consumption in LANG and LWHs. Hearts were paced to increase workload and perfusate was deoxygenated to study the effects of myocardial hypoxia. A subset of hearts was perfused with 1 μm tolbutamide (TOLB) to identify the level of AP duration (APD) shortening attributed to K_ATP channel activation. During sinus rhythm, APD was shorter in LWHs compared to LANG hearts. APD in both LWHs and LANG hearts dropped steadily during deoxygenation. With TOLB, APDs in LWHs were longer at all workloads and APD reductions during deoxygenation were blunted in both LWHs and LANG hearts. At 50% perfusate oxygenation, APD and LVDP were significantly higher in LWHs perfused with TOLB (199 ± 16 ms; 92 ± 5.3 mmHg) than in LWHs without TOLB (109 ± 14 ms, P = 0.005; 65 ± 6.5 mmHg, P = 0.01). Our results indicate that K_ATP channels are activated to a greater extent in perfused hearts when the LV performs pressure-volume work. The results of the present study demonstrate the critical role of K_ATP channels in modulating myocardial function over a wide range of physiological conditions.

Electrophysiological properties of computational human ventricular cell action potential models under acute ischemic conditions

Acute myocardial ischemia is one of the main causes of sudden cardiac death. The mechanisms have been investigated primarily in experimental and computational studies using different animal species, but human studies remain scarce. In this study, we assess the ability of four human ventricular action potential models (ten Tusscher and Panfilov, 2006; Grandi et al., 2010; Carro et al., 2011; O'Hara et al., 2011) to simulate key electrophysiological consequences of acute myocardial ischemia in single cell and tissue simulations. We specifically focus on evaluating the effect of extracellular potassium concentration and activation of the ATP-sensitive inward-rectifying potassium current on action potential duration, post-repolarization refractoriness, and conduction velocity, as the most critical factors in determining reentry vulnerability during ischemia. Our results show that the Grandi and O'Hara models required modifications to reproduce expected ischemic changes, specifically modifying the intracellular potassium concentration in the Grandi model and the sodium current in the O'Hara model. With these modifications, the four human ventricular cell AP models analyzed in this study reproduce the electrophysiological alterations in repolarization, refractoriness, and conduction velocity caused by acute myocardial ischemia. However, quantitative differences are observed between the models and overall, the ten Tusscher and modified O'Hara models show closest agreement to experimental data.

Depth Attenuation Degree Based Visualization for Cardiac Ischemic Electrophysiological Feature Exploration

Although heart researches and acquirement of clinical and experimental data are progressively open to public use, cardiac biophysical functions are still not well understood. Due to the complex and fine structures of the heart, cardiac electrophysiological features of interest may be occluded when there is a necessity to demonstrate cardiac electrophysiological behaviors. To investigate cardiac abnormal electrophysiological features under the pathological condition, in this paper, we implement a human cardiac ischemic model and acquire the electrophysiological data of excitation propagation. A visualization framework is then proposed which integrates a novel depth weighted optic attenuation model into the pathological electrophysiological model. The hidden feature of interest in pathological tissue can be revealed from sophisticated overlapping biophysical information. Experiment results verify the effectiveness of the proposed method for intuitively exploring and inspecting cardiac electrophysiological activities, which is fundamental in analyzing and explaining biophysical mechanisms of cardiac functions for doctors and medical staff.

Non-cell autonomous cues for enhanced functionality of human embryonic stem cell-derived cardiomyocytes via maturation of sarcolemmal and mitochondrial KATP channels

Human embryonic stem cells (hESCs) is a potential unlimited ex vivo source of ventricular (V) cardiomyocytes (CMs), but hESC-VCMs and their engineered tissues display immature traits. In adult VCMs, sarcolemmal (sarc) and mitochondrial (mito) ATP-sensitive potassium (K_ATP) channels play crucial roles in excitability and cardioprotection. In this study, we aim to investigate the biological roles and use of sarcK_ATP and mitoK_ATP in hESC-VCM. We showed that SarcI_{K, ATP} in single hESC-VCMs was dormant under baseline conditions, but became markedly activated by cyanide (CN) or the known opener P1075 with a current density that was ~8-fold smaller than adult; These effects were reversible upon washout or the addition of GLI or HMR1098. Interestingly, sarcI_{K, ATP} displayed a ~3-fold increase after treatment with hypoxia (5% O₂). MitoI_{K, ATP} was absent in hESC-VCMs. However, the thyroid hormone T3 up-regulated mitoI_{K, ATP,} conferring diazoxide protective effect on T3-treated hESC-VCMs. When assessed using a multi-cellular engineered 3D ventricular cardiac micro-tissue (hvCMT) system, T3 substantially enhanced the developed tension by 3-folds. Diazoxide also attenuated the decrease in contractility induced by simulated ischemia (1% O₂). We conclude that hypoxia and T3 enhance the functionality of hESC-VCMs and their engineered tissues by selectively acting on sarc and mitoI_{K, ATP}.

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.

Restrictions in ATP diffusion within sarcomeres can provoke ATP-depleted zones impairing exercise capacity in chronic obstructive pulmonary disease

BACKGROUND

Chronic obstructive pulmonary disease (COPD) is characterized by the inability of patients to sustain a high level of ventilation resulting in perceived exertional discomfort and limited exercise capacity of leg muscles at average intracellular ATP levels sufficient to support contractility.

METHODS

Myosin ATPase activity in biopsy samples from healthy and COPD individuals was implemented as a local nucleotide sensor to determine ATP diffusion coefficients within myofibrils. Ergometric parameters clinically measured during maximal exercise tests in both groups were used to define the rates of myosin ATPase reaction and aerobic ATP re-synthesis. The obtained parameters in combination with AK- and CK-catalyzed reactions were implemented to compute the kinetic and steady-state spatial ATP distributions within control and COPD sarcomeres.

RESULTS

The developed reaction-diffusion model of two-dimensional sarcomeric space identified similar, yet extremely low nucleotide diffusion in normal and COPD myofibrils. The corresponding spatio-temporal ATP distributions, constructed during imposed exercise, predicted in COPD sarcomeres a depletion of ATP in the zones of overlap between actin and myosin filaments along the center axis at average cytosolic ATP levels similar to healthy muscles.

CONCLUSIONS

ATP-depleted zones can induce rigor tension foci impairing muscle contraction and increase a risk for sarcomere damages. Thus, intra-sarcomeric diffusion restrictions at limited aerobic ATP re-synthesis can be an additional risk factor contributing to the muscle contractile deficiency experienced by COPD patients.

GENERAL SIGNIFICANCE

This study demonstrates how restricted substrate mobility within a cellular organelle can provoke an energy imbalance state paradoxically occurring at abounding average metabolic resources.

Conduction abnormalities and ventricular arrhythmogenesis: The roles of sodium channels and gap junctions

Ventricular arrhythmias arise from disruptions in the normal orderly sequence of electrical activation and recovery of the heart. They can be categorized into disorders affecting predominantly cellular depolarization or repolarization, or those involving action potential (AP) conduction. This article briefly discusses the factors causing conduction abnormalities in the form of unidirectional conduction block and reduced conduction velocity (CV). It then examines the roles that sodium channels and gap junctions play in AP conduction. Finally, it synthesizes experimental results to illustrate molecular mechanisms of how abnormalities in these proteins contribute to such conduction abnormalities and hence ventricular arrhythmogenesis, in acquired pathologies such as acute ischaemia and heart failure, as well as inherited arrhythmic syndromes.

Histone deacetylase inhibitors modulate KATP subunit transcription in HL-1 cardiomyocytes through effects on cholesterol homeostasis

Histone deacetylase inhibitors (HDIs) are under investigation for the treatment of a number of human health problems. HDIs have proven therapeutic value in refractory cases of cutaneous T-cell lymphoma. Electrocardiographic ST segment morphological changes associated with HDIs were observed during development. Because ST segment morphology is typically linked to changes in ATP sensitive potassium (K_ATP) channel activity, we tested the hypothesis that HDIs affect cardiac K_ATP channel subunit expression. Two different HDIs, romidepsin and trichostatin A, caused ~20-fold increase in SUR2 (Abcc9) subunit mRNA expression in HL-1 cardiomyocytes. The effect was specific for the SUR2 subunit as neither compound causes a marked change in SUR1 (Abcc8) expression. Moreover, the effect was cell specific as neither HDI markedly altered K_ATP subunit expression in MIN6 pancreatic β-cells. We observe significant enrichment of the H3K9Ac histone mark specifically at the SUR2 promoter consistent with the conclusion that chromatin remodeling at this locus plays a role in increasing SUR2 gene expression. Unexpectedly, however, we also discovered that HDI-dependent depletion of cellular cholesterol is required for the observed effects on SUR2 expression. Taken together, the data in the present study demonstrate that K_ATP subunit expression can be epigenetically regulated in cardiomyocytes, defines a role for cholesterol homeostasis in mediating epigenetic regulation and suggests a potential molecular basis for the cardiac effects of the HDIs.

EFFECTS OF ACUTE GLOBAL ISCHEMIA ON RE-ENTRANT ARRHYTHMOGENESIS: A SIMULATION STUDY

Sudden cardiac death is mainly caused by arrhythmogenesis. For a functional abnormal heart, such as an ischemic heart, the probability of arrhythmia occurring is greatly increased. During myocardial ischemia, re-entry is prone to degenerate into ventricular fibrillation (VF). Therefore it has important meaning to investigate the intricate mechanisms underlying VF under an ischemic condition in order to better facilitate therapeutic interventions. In this paper, to analyze the functional influence of acute global ischemia on cardiac electrical activity and subsequently on re-entrant arrhythmogenesis, we take into account three main pathophysiological consequences of ischemia: hyperkalaemia, acidosis, and anoxia, and develop a 3D human ventricular ischemic model that combines a detailed biophysical description of the excitation kinetics of human ventricular cells with an integrated geometry of human ventricular tissue which incorporates fiber direction anisotropy and the stimulation activation sequence. The results show that under acute global ischemia, the tissue excitability and the slope of ventricular cellular action potential duration restitution (APDR) are greatly decreased. As a result, the complexity of VF activation patterns is reduced. For the three components of ischemia, hyperkalaemia is the dominant contributor to the stability of re-entry under acute global ischemia. Increasing [K+]o acts to prolong the cell refractory period, reduce the tissue excitability and slow the conduction velocity. Our results also show that VF can be eliminated by decreasing cellular excitability, primarily by elevating the concentration value of extracellular K+.

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.

Kir6.2 limits Ca(2+) overload and mitochondrial oscillations of ventricular myocytes in response to metabolic stress

ATP-sensitive K(+) (KATP) channels are abundant membrane proteins in cardiac myocytes that are directly gated by intracellular ATP and form a signaling complex with metabolic enzymes, such as creatine kinase. KATP channels are known to be essential for adaption to cardiac stress, such as ischemia; however, how all the molecular components of the stress response interact is not fully understood. We examined the effects of decreasing the KATP current density on Ca(2+) and mitochondrial homeostasis and ischemic preconditioning. Acute knockdown of the pore-forming subunit, Kir6.2, was achieved using adenoviral delivery of short hairpin RNA targeted to Kir6.2. The acute nature of the knockdown of Kir6.2 accurately shows the effects of Kir6.2 depletion without any compensatory effects that may arise in transgenic studies. We also investigated the effect of reducing the KATP current while maintaining KATP channel protein in the sarcolemmal membrane using a nonconducting Kir6.2 construct. Only 50% KATP current remained after Kir6.2 knockdown, yet there were profound effects on myocyte responses to metabolic stress. Kir6.2 was essential for cardiac myocyte Ca(2+) homeostasis under both baseline conditions before any metabolic stress and after metabolic stress. Expression of nonconducting Kir6.2 also resulted in increased Ca(2+) overload, showing the importance of K(+) conductance in the protective response. Both ischemic preconditioning and protection during ischemia were lost when Kir6.2 was knocked down. KATP current density was also important for the mitochondrial membrane potential at rest and prevented mitochondrial membrane potential oscillations during oxidative stress. KATP channel density is important for adaption to metabolic stress.

Heterogeneity and function of K(ATP) channels in canine hearts

BACKGROUND

The concept that pore-forming Kir6.2 and regulatory SUR2A subunits form cardiac ATP-sensitive potassium (K(ATP)) channels is challenged by recent reports that SUR1 is predominant in mouse atrial K(ATP) channels.

OBJECTIVE

To assess SUR subunit composition of K(ATP) channels and consequence of K(ATP) activation for action potential duration (APD) in dog hearts.

METHODS

Patch-clamp techniques were used on isolated dog cardiomyocytes to investigate K(ATP) channel properties. Dynamic current clamp, by injection of a linear K(+) conductance to simulate activation of the native current, was used to study the consequences of K(ATP) activation on APD.

RESULTS

Metabolic inhibitor (MI)-activated current was not significantly different from pinacidil (SUR2A-specific)-activated current, and both currents were larger than diazoxide (SUR1-specific)-activated current in both the atrium and the ventricle. Mean K(ATP) conductance (activated by MI) did not differ significantly between chambers, although, within the ventricle, both MI-induced and pinacidil-induced currents tended to decrease from the epicardium to the endocardium. Dynamic current-clamp results indicate that myocytes with longer baseline APDs are more susceptible to injected K(ATP) current, a result reproduced in silico by using a canine action potential model (Hund-Rudy) to simulate epicardial and endocardial myocytes.

CONCLUSIONS

Even a small fraction of K(ATP) activation significantly shortens APD in a manner that depends on existing heterogeneity in K(ATP) current and APD.

Mechanical effects on KATP channel gating in rat ventricular myocytes

Cardiac K_ATP channels link metabolism with electrical activity. They are implicated in arrhythmias, secretion of atrial natriuretic peptide and protection of the heart from hypertrophy and failure. These processes may involve mechanosensitivity. K_ATP channels can be activated by mechanical stimulation and disrupting the cortical actin increases the activity. We propose that K_ATP channels are modulated by local bilayer tension and this tension is affected by cortical F-actin. Here we measured K_ATP background activity and stretch sensitivity with inside-out patches of rat ventricular myocytes before and after disrupting F-actin. Disrupting F-actin potentiated background activity but did not influence the slope sensitivity in the semilog relationship of NP_o vs. suction that is a measure of the change in dimensions between closed and open states. Thus actin alters prestress on the channel probably by parallel elastic sharing of mean cortical tension with the bilayer.

The ischemic heart: what causes ectopic beating?

The mechanisms by which spontaneous electrical activity originates in the ischemic heart and leads to arrhythmia remain unknown, however mechanical stretch of the diseased region has been hypothesized to play a role. The goal of this study is to investigate the conditions that favor the initiation of stretch-induced premature beats in the ischemic heart. We employ a mathematical model of the ischemic cell subjected to stretch. The study found that upon stretch, spontaneous beats occur in the ischemic cell, which are due to the stretch-induced re-activation of the L-type calcium current.

Toward an integrative computational model of the Guinea pig cardiac myocyte

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

Compartmentation of membrane processes and nucleotide dynamics in diffusion-restricted cardiac cell microenvironment

Orchestrated excitation-contraction coupling in heart muscle requires adequate spatial arrangement of systems responsible for ion movement and metabolite turnover. Co-localization of regulatory and transporting proteins into macromolecular complexes within an environment of microanatomical cell components raises intracellular diffusion barriers that hamper the mobility of metabolites and signaling molecules. Compared to substrate diffusion in the cytosol, diffusional restrictions underneath the sarcolemma are much larger and could impede ion and nucleotide movement by a factor of 10(3)-10(5). Diffusion barriers thus seclude metabolites within the submembrane space enabling rapid and vectorial effector targeting, yet hinder energy supply from the bulk cytosolic space implicating the necessity for a shunting transfer mechanism. Here, we address principles of membrane protein compartmentation, phosphotransfer enzyme-facilitated interdomain energy transfer, and nucleotide signal dynamics at the subsarcolemma-cytosol interface. This article is part of a Special Issue entitled "Local Signaling in Myocytes".

Mitochondrial uncoupling agents trigger ventricular fibrillation in isolated rat hearts

Sudden cardiac death resulting from ventricular fibrillation (VF) remains a major cause of mortality. The purpose of this study was to investigate the roles of loss of oxidative phosphorylation and activation of the mitochondrial ATP-sensitive K+ channel and permeability transition pore in VF development during myocardial ischemia by using mitochondrial uncoupling agents (carbonyl cyanide m-chlorophenylhydrazone and 2,4-dinitrophenol) and channel blockers (5-hydroxydecanoate and cyclosporine A) at concentrations that have been demonstrated to block the intended targets selectively. Isolated rat hearts (n = 8 per group) were perfused with 0.3 μM carbonyl cyanide m-chlorophenylhydrazone, 100 μM 2,4-dinitrophenol, 0.2 μM cyclosporine A, 100 μM 5-hydroxydecanoate, or vehicle solution and regional ischemia induced after 10 minutes. Carbonyl cyanide m-chlorophenylhydrazone and 2,4 dinitrophenol caused profound QT shortening and triggered VF in 100% of hearts before ischemia. During ischemia, neither cyclosporine A (88%) nor 5-hydroxydecanoate (100%) reduced VF incidence compared with control (100% VF). In separate hearts, carbonyl cyanide m-chlorophenylhydrazone decreased tissue ATP content, and glibenclamide or glimepiride delayed the QT shortening and onset of VF triggered by carbonyl cyanide m-chlorophenylhydrazone. In conclusion, mitochondrial uncoupling agents trigger VF, likely as a result of ATP depletion with subsequent activation of sarcolemmal ATP-sensitive K+ currents. The mechanism of VF in ischemia does not involve activation of the mitochondrial ATP-sensitive K+ channel or permeability transition pore.

Experiment-model interaction for analysis of epicardial activation during human ventricular fibrillation with global myocardial ischaemia

We describe a combined experiment-modelling framework to investigate the effects of ischaemia on the organisation of ventricular fibrillation in the human heart. In a series of experimental studies epicardial activity was recorded from 10 patients undergoing routine cardiac surgery. Ventricular fibrillation was induced by burst pacing, and recording continued during 2.5 min of global cardiac ischaemia followed by 30 s of coronary reflow. Modelling used a 2D description of human ventricular tissue. Global cardiac ischaemia was simulated by (i) decreased intracellular ATP concentration and subsequent activation of an ATP sensitive K⁺ current, (ii) elevated extracellular K⁺ concentration, and (iii) acidosis resulting in reduced magnitude of the L-type Ca²⁺ current I(Ca,L). Simulated ischaemia acted to shorten action potential duration, reduce conduction velocity, increase effective refractory period, and flatten restitution. In the model, these effects resulted in slower re-entrant activity that was qualitatively consistent with our observations in the human heart. However, the flattening of restitution also resulted in the collapse of many re-entrant waves to several stable re-entrant waves, which was different to the overall trend we observed in the experimental data. These findings highlight a potential role for other factors, such as structural or functional heterogeneity in sustaining wavebreak during human ventricular fibrillation with global myocardial ischaemia.

K(ATP) channels process nucleotide signals in muscle thermogenic response

Uniquely gated by intracellular adenine nucleotides, sarcolemmal ATP-sensitive K(+) (K(ATP)) channels have been typically assigned to protective cellular responses under severe energy insults. More recently, K(ATP) channels have been instituted in the continuous control of muscle energy expenditure under non-stressed, physiological states. These advances raised the question of how K(ATP) channels can process trends in cellular energetics within a milieu where each metabolic system is set to buffer nucleotide pools. Unveiling the mechanistic basis of the K(ATP) channel-driven thermogenic response in muscles thus invites the concepts of intracellular compartmentalization of energy and proteins, along with nucleotide signaling over diffusion barriers. Furthermore, it requires gaining insight into the properties of reversibility of intrinsic ATPase activity associated with K(ATP) channel complexes. Notwithstanding the operational paradigm, the homeostatic role of sarcolemmal K(ATP) channels can be now broadened to a wider range of environmental cues affecting metabolic well-being. In this way, under conditions of energy deficit such as ischemic insult or adrenergic stress, the operation of K(ATP) channel complexes would result in protective energy saving, safeguarding muscle performance and integrity. Under energy surplus, downregulation of K(ATP) channel function may find potential implications in conditions of energy imbalance linked to obesity, cold intolerance and associated metabolic disorders.

Muscle KATP channels: recent insights to energy sensing and myoprotection

ATP-sensitive potassium (K(ATP)) channels are present in the surface and internal membranes of cardiac, skeletal, and smooth muscle cells and provide a unique feedback between muscle cell metabolism and electrical activity. In so doing, they can play an important role in the control of contractility, particularly when cellular energetics are compromised, protecting the tissue against calcium overload and fiber damage, but the cost of this protection may be enhanced arrhythmic activity. Generated as complexes of Kir6.1 or Kir6.2 pore-forming subunits with regulatory sulfonylurea receptor subunits, SUR1 or SUR2, the differential assembly of K(ATP) channels in different tissues gives rise to tissue-specific physiological and pharmacological regulation, and hence to the tissue-specific pharmacological control of contractility. The last 10 years have provided insights into the regulation and role of muscle K(ATP) channels, in large part driven by studies of mice in which the protein determinants of channel activity have been deleted or modified. As yet, few human diseases have been correlated with altered muscle K(ATP) activity, but genetically modified animals give important insights to likely pathological roles of aberrant channel activity in different muscle types.

Modeling cardiac action potential shortening driven by oxidative stress-induced mitochondrial oscillations in guinea pig cardiomyocytes

Ischemia-induced shortening of the cardiac action potential and its heterogeneous recovery upon reperfusion are thought to set the stage for reentrant arrhythmias and sudden cardiac death. We have recently reported that the collapse of mitochondrial membrane potential (DeltaPsi(m)) through a mechanism triggered by reactive oxygen species (ROS), coupled to the opening of sarcolemmal ATP-sensitive potassium (K(ATP)) channels, contributes to electrical dysfunction during ischemia-reperfusion. Here we present a computational model of excitation-contraction coupling linked to mitochondrial bioenergetics that incorporates mitochondrial ROS-induced ROS release with coupling between the mitochondrial energy state and electrical excitability mediated by the sarcolemmal K(ATP) current (I(K,ATP)). Whole-cell model simulations demonstrate that increasing the fraction of oxygen diverted from the respiratory chain to ROS production triggers limit-cycle oscillations of DeltaPsi(m), redox potential, and mitochondrial respiration through the activation of a ROS-sensitive inner membrane anion channel. The periods of transient mitochondrial uncoupling decrease the cytosolic ATP/ADP ratio and activate I(K,ATP), consequently shortening the cellular action potential duration and ultimately suppressing electrical excitability. The model simulates emergent behavior observed in cardiomyocytes subjected to metabolic stress and provides a new tool for examining how alterations in mitochondrial oxidative phosphorylation will impact the electrophysiological, contractile, and Ca(2+) handling properties of the cardiac cell. Moreover, the model is an important step toward building multiscale models that will permit investigation of the role of spatiotemporal heterogeneity of mitochondrial metabolism in the mechanisms of arrhythmogenesis and contractile dysfunction in cardiac muscle.

Cardiac sarcolemmal K(ATP) channels: Latest twists in a questing tale!

Reconstitution of K(ATP) channel activity from coexpression of members of the pore-forming inward rectifier gene family (Kir6.1, KCNJ8, and Kir6.2 KCNJ11) with sulfonylurea receptors (SUR1, ABCC8, and SUR2, ABCC9) of the ABCC protein sub-family, has led to the elucidation of many details of channel gating and pore properties, as well as the essential roles of Kir6.2 and SUR2 subunits in generating cardiac ventricular K(ATP). However, despite this extensive body of knowledge, there remain significant holes in our understanding of the physiological role of the cardiac K(ATP) channel, and surprising new findings keep emerging. Recent findings from genetically modified animals include the apparent insensitivity of cardiac sarcolemmal channels to nucleotide levels, and unenvisioned complexities of the subunit make-up of the cardiac channels. This topical review focuses on these new findings and considers their implications.

The relative role of refractoriness and source-sink relationship in reentry generation during simulated acute ischemia

During acute myocardial ischemia, reentrant episodes may lead to ventricular fibrillation (VF), giving rise to potentially mortal arrhythmias. VF has been traditionally related to dispersion of refractoriness and more recently to the source-sink relationship. Our goal is to theoretically investigate the relative role of dispersion of refractoriness and source-sink mismatch in vulnerability to reentry in the specific situation of regional myocardial acute ischemia. The electrical activity of a regionally ischemic tissue was simulated using a modified version of the Luo-Rudy dynamic model. Ischemic conditions were varied to simulate the time-course of acute ischemia. Our results showed that dispersion of refractoriness increased with the severity of ischemia. However, no correlation between dispersion of refractoriness and the width of the vulnerable window was found. Additionally, in approximately 50% of the reentries, unidirectional block (UDB) took place in cells completely recovered from refractoriness. We examined patterns of activation after premature stimulation and they were intimately related to the source-sink relationship, quantified by the safety factor (SF). Moreover, the isoline where the SF dropped below unity matched the area where propagation failed. It was concluded that the mismatch of the source-sink relationship, rather than solely refractoriness, was the ultimate cause of the UDB leading to reentry. The SF represents a very powerful tool to study the mechanisms responsible for reentry.

Ventricular fibrillation with prominent early repolarization associated with a rare variant of KCNJ8/KATP channel

BACKGROUND

Early repolarization in the inferolateral leads has been recently recognized as a frequent syndrome associated with idiopathic ventricular fibrillation (VF). We report the case of a patient presenting dramatic changes in the ECG in association with recurrent VF in whom a novel genetic variant has been identified.

CASE REPORT

This young female (14 years) was resuscitated in 2001 following an episode of sudden death due to VF. All examinations including coronary angiogram with ergonovine injection, MRI, and flecainide or isoproterenol infusion were normal. The patient had multiple (>100) recurrences of VF unresponsive to beta-blockers, lidocaine/mexiletine, verapamil, and amiodarone. Recurrences of VF were associated with massive accentuation of the early repolarization pattern at times mimicking acute myocardial ischemia. Coronary angiography during an episode with 1.2 mV J/ST elevation was normal. Isoproterenol infusion acutely suppressed electrical storms, while quinidine eliminated all recurrences of VF and restored a normal ECG over a follow-up of 65 months. Genomic DNA sequencing of K(ATP) channel genes showed missense variant in exon 3 (NC_000012) of the KCNJ8 gene, a subunit of the K(ATP) channel, conferring predisposition to dramatic repolarization changes and ventricular vulnerability.

Mechanistic investigation into the arrhythmogenic role of transmural heterogeneities in regional ischaemia phase 1A

AIMS

Studies of arrhythmogenesis during ischemia have focused primarily on reentrant mechanisms manifested on the epicardial surface. The goal of this study was to use a physiologically-accurate model of acute regional ischemia phase 1A to determine the contribution of ischaemia-induced transmural electrophysiological heterogeneities to arrhythmogenesis following left anterior descending artery occlusion.

METHODS AND RESULTS

A slice through a geometrical model of the rabbit ventricles was extracted and a model of regional ischaemia developed. The model included a central ischaemic zone incorporating transmural gradients of I(K(ATP)) activation and [K+]o, surrounded by ischaemic border zones (BZs), with the degree of ischaemic effects varied to represent progression of ischaemia 2-10 min post-occlusion. Premature stimulation was applied over a range of coupling intervals to induce re-entry. The presence of ischaemic BZs and a transmural gradient in I(K(ATP)) activation provided the substrate for re-entrant arrhythmias. Increased dispersion of refractoriness and conduction velocity in the BZs with time post-occlusion led to a progressive increase in arrhythmogenesis. In the absence of a transmural gradient of I(K(ATP)) activation, re-entry was rarely sustained.

CONCLUSION

Knowledge of the mechanism by which specific electrophysiological heterogeneities underlie arrhythmogenesis during acute ischaemia could be useful in developing preventative treatments for patients at risk of coronary vascular disease.

ATP-sensitive potassium channels: metabolic sensing and cardioprotection

The cardiovascular system operates under a wide scale of demands, ranging from conditions of rest to extreme stress. How the heart muscle matches rates of ATP production with utilization is an area of active investigation. ATP-sensitive potassium (K(ATP)) channels serve a critical role in the orchestration of myocardial energetic well-being. K(ATP) channel heteromultimers consist of inwardly-rectifying K(+) channel 6.2 and ATP-binding cassette sulfonylurea receptor 2A that translates local ATP/ADP levels, set by ATPases and phosphotransfer reactions, to the channel pore function. In cells in which the mobility of metabolites between intracellular microdomains is limited, coupling of phosphotransfer pathways with K(ATP) channels permits a high-fidelity transduction of nucleotide fluxes into changes in membrane excitability, matching energy demands with metabolic resources. This K(ATP) channel-dependent optimization of cardiac action potential duration preserves cellular energy balance at varying workloads. Mutations of K(ATP) channels result in disruption of the nucleotide signaling network and generate a stress-vulnerable phenotype with excessive susceptibility to injury, development of cardiomyopathy, and arrhythmia. Solving the mechanisms underlying the integration of K(ATP) channels into the cellular energy network will advance the understanding of endogenous cardioprotection and the development of strategies for the management of cardiovascular injury and disease progression.

Vulnerability to reentry in a regionally ischemic tissue: a simulation study

Sudden cardiac death is mainly provoked by arrhythmogenic processes. During myocardial ischemia many malignant arrhythmias, such as reentry, take place and can degenerate into ventricular fibrillation. It is thus of great interest to unravel the intricate mechanisms underlying the initiation and maintenance of a reentry. In this computational study, we analyze the probability of reentry during different stages of the acute phase of ischemia. We also aimed at the understanding of the role of its main components: hypoxia, hyperkalemia, and acidosis analyzing the intricate ionic mechanisms responsible for reentry generation. We simulated the electrical activity of a ventricular tissue affected by regional ischemia based on a modified version of the Luo-Rudy model (LRd00). The ischemic conditions were varied to simulate different stages of this pathology. After premature stimulation, we evaluated the vulnerability to reentry. We obtained an unimodal behavior for the vulnerable window as ischemia progressed, peaking at the eighth minute after the onset of ischemia where the vulnerable window yielded 58 ms. Under more severe conditions the vulnerable window decreased and became zero for minute 8.75. The present work provides insight into the mechanisms of reentry generation during ischemia, highlighting the role of acidosis and hypoxia when hyperkalemia is present.

Effect of Ca2+ on cardiac mitochondrial energy production is modulated by Na+ and H+ dynamics

The energy production of mitochondria in heart increases during exercise. Several works have suggested that calcium acts at multiple control points to activate net ATP production in what is termed "parallel activation". To study this, a computational model of mitochondrial energy metabolism in the heart has been developed that integrates the Dudycha-Jafri model for the tricarboxylic acid cycle with the Magnus-Keizer model for mitochondrial energy metabolism and calcium dynamics. The model improves upon the previous formulation by including an updated formulation for calcium dynamics, and new descriptions of sodium, hydrogen, phosphate, and ATP balance. To this end, it incorporates new formulations for the calcium uniporter, sodium-calcium exchange, sodium-hydrogen exchange, the F(1)F(0)-ATPase, and potassium-hydrogen exchange. The model simulates a wide range of experimental data, including steady-state and simulated pacing protocols. The model suggests that calcium is a potent activator of net ATP production and that as pacing increases energy production due to calcium goes up almost linearly. Furthermore, it suggests that during an extramitochondrial calcium transient, calcium entry and extrusion cause a transient depolarization that serve to increase NADH production by the tricarboxylic acid cycle and NADH consumption by the respiration driven proton pumps. The model suggests that activation of the F(1)F(0)-ATPase by calcium is essential to increase ATP production. In mitochondria very close to the release sites, the depolarization is more severe causing a temporary loss of ATP production. However, due to the short duration of the depolarization the net ATP production is also increased.

Modeling of spatial metabolite distributions in the cardiac sarcomere

Although a high ATP diffusion rate implies homogeneous distribution of the principal energetic currency in the cytosol, local diffusion barriers represented by macromolecular structures can render ATP concentrations to be inhomogeneous. A method is presented here that provides apparent diffusion coefficient values in local intracellular regions and allows the estimation of spatial metabolite distribution. The apparent local diffusion coefficient for ATP in cardiac myofibrils was determined from the analysis of diffusion-dependent rightward shift of the substrate dependence for actomyosin ATPase activity using the reaction-diffusion model, which accounted for the properties of phosphotransfer reactions. This functional analysis, which took into account the local diffusional ATP delivery to the active sites, provided an apparent value that was three orders of magnitude lower than that defined by direct methods for the cytosol. The low value of the diffusion coefficient was shown to define unusual properties of the intracellular space in working heart, where small reductions in ATP levels in the surrounding cytosol result in a large drop in [ATP] inside myofibrils. This drop is critical for vital cellular functions, and the analysis presented here defines its physical basis. The diffusion barriers thus defined explain the coexistence of pathological energy deficit with almost normal average ATP levels.

Modeling cardiac ischemia

Myocardial ischemia is one of the main causes of sudden cardiac death, with 80% of victims suffering from coronary heart disease. In acute myocardial ischemia, the obstruction of coronary flow leads to the interruption of oxygen flow, glucose, and washout in the affected tissue. Cellular metabolism is impaired and severe electrophysiological changes in ionic currents and concentrations ensue, which favor the development of lethal cardiac arrhythmias such as ventricular fibrillation. Due to the burden imposed by ischemia in our societies, a large body of research has attempted to unravel the mechanisms of initiation, sustenance, and termination of cardiac arrhythmias in acute ischemia, but the rapidity and complexity of ischemia-induced changes as well as the limitations in current experimental techniques have hampered evaluation of ischemia-induced alterations in cardiac electrical activity and understanding of the underlying mechanisms. Over the last decade, computer simulations have demonstrated the ability to provide insight, with high spatiotemporal resolution, into ischemic abnormalities in cardiac electrophysiological behavior from the ionic channel to the whole organ. This article aims to review and summarize the results of these studies and to emphasize the role of computer simulations in improving the understanding of ischemia-related arrhythmias and how to efficiently terminate them.

Ketanserin, a 5-HT2 antagonist, directly inhibits the ATP-sensitive potassium channel in mouse ventricular myocytes

The effects of ketanserin, a 5-HT2 antagonist, on the ATP-sensitive K+ (K(ATP)) channels were studied in mouse ventricular myocytes using patch clamp technique. Under the whole-cell voltage clamp conditions, ketanserin (1-100 microM) reversibly inhibited pinacidil-induced K(ATP) current in a concentration-dependent fashion with a Ki value of 9.36 microM and the Hill coefficient was 0.67. This inhibition was developed even with the presence of 5-hydroxytryptamine (100 microM) in the bath. Prazosin, a selective alpha1-antagonist, also failed to mimic the effect of ketanserin. Ketanserin did not affect the channel activity in inside-out configuration under the ATP-free internal solution. Furthermore, ketanserin applied to the external solution did not affect the pinacidil-induced channel activity in the cell-attached patches, but did inhibit it when applied into the pipette. These results suggest that the inhibitory action of ketanserin observed in this study was probably due to a direct action on the K(ATP) channel rather than to an action through the 5-HT2 receptor or alpha1-adrenoceptor blockade, and that the antiarrhythmic activity of ketanserin against cardiac arrhythmias induced in the ischemic/reperfused heart is at least in part attributable to its inhibition of the K(ATP) channel.

Effects of pinacidil on reentrant arrhythmias generated during acute regional ischemia: a simulation study

Many experimental studies have pointed out the controversy involving the arrhythmogenic effects of potassium channel openers (KCOs) in ischemia. KCOs activate the ATP-sensitive potassium current [IK(ATP)], resulting in action potential duration (APD) shortening, especially under pathological conditions such as ischemia. Acute myocardial ischemia leads to electrophysiological inhomogeneities in APD, conduction velocity, and refractoriness, which provide the substrate for reentry initiation and maintenance and may lead to malignant arrhythmias. The aim of this work is to analyze the effect of the KCO pinacidil on vulnerability to reentry during acute regional ischemia using computer simulations. We use a two-dimensional virtual heart tissue with implementation of acute regional ischemia conditions. Membrane kinetics are represented by a modified version of Luo-Rudy (phase II) action potential model that incorporates the effect of pinacidil on IK(ATP). The vulnerable window (VW) for reentry is quantified for different doses of pinacidil. Our results show that for doses below 3 micromol/l the VW widens with increasing pinacidil concentration, whereas for higher doses of pinacidil the VW decreases, becoming zero for concentrations above 10 micromol/l. The ionic mechanisms involved in this behavior are explored. This study demonstrates that the effect of pinacidil on arrhythmogenesis is strongly dose-dependent, and that high doses of pinacidil exert a strong antiarrhythmic effect.

Modeling regulation of cardiac KATP and L-type Ca2+ currents by ATP, ADP, and Mg2+

Changes in cytosolic free Mg(2+) and adenosine nucleotide phosphates affect cardiac excitability and contractility. To investigate how modulation by Mg(2+), ATP, and ADP of K(ATP) and L-type Ca(2+) channels influences excitation-contraction coupling, we incorporated equations for intracellular ATP and MgADP regulation of the K(ATP) current and MgATP regulation of the L-type Ca(2+) current in an ionic-metabolic model of the canine ventricular myocyte. The new model: 1), quantitatively reproduces a dose-response relationship for the effects of changes in ATP on K(ATP) current, 2), simulates effects of ADP in modulating ATP sensitivity of K(ATP) channel, 3), predicts activation of Ca(2+) current during rapid increase in MgATP, and 4), demonstrates that decreased ATP/ADP ratio with normal total Mg(2+) or increased free Mg(2+) with normal ATP and ADP activate K(ATP) current, shorten action potential, and alter ionic currents and intracellular Ca(2+) signals. The model predictions are in agreement with experimental data measured under normal and a variety of pathological conditions.