A computational model of pig ventricular cardiomyocyte electrophysiology and calcium handling: Translation from pig to human electrophysiology.

The pig is commonly used as an experimental model of human heart disease, including for the study of mechanisms of arrhythmia. However, there exist differences between human and porcine cellular electrophysiology: The pig action potential (AP) has a deeper phase-1 notch, a longer duration at 50% repolarization, and higher plateau potentials than human. Ionic differences underlying the AP include larger rapid delayed-rectifier and smaller inward-rectifier K+-currents (IKr and IK1 respectively) in humans. AP steady-state rate-dependence and restitution is steeper in pigs. Porcine Ca2+ transients can have two components, unlike human. Although a reliable computational model for human ventricular myocytes exists, one for pigs is lacking. This hampers translation from results obtained in pigs to human myocardium. Here, we developed a computational model of the pig ventricular cardiomyocyte AP using experimental datasets of the relevant ionic currents, Ca2+-handling, AP shape, AP duration restitution, and inducibility of triggered activity and alternans. To properly capture porcine Ca2+ transients, we introduced a two-step process with a faster release in the t-tubular region, followed by a slower diffusion-induced release from a non t-tubular subcellular region. The pig model behavior was compared with that of a human ventricular cardiomyocyte (O'Hara-Rudy) model. The pig, but not the human model, developed early afterdepolarizations (EADs) under block of IK1, while IKr block led to EADs in the human but not in the pig model. At fast rates (pacing cycle length = 400 ms), the human cell model was more susceptible to spontaneous Ca2+ release-mediated delayed afterdepolarizations (DADs) and triggered activity than pig. Fast pacing led to alternans in human but not pig. Developing species-specific models incorporating electrophysiology and Ca2+-handling provides a tool to aid translating antiarrhythmic and arrhythmogenic assessment from the bench to the clinic.

Validation of quantitative measure of repolarization reserve as a novel marker of drug induced proarrhythmia.

Repolarization reserve, the robustness of a cell to repolarize even when one of the repolarization mechanisms is failing, has been described qualitatively in terms of ionic currents, but has not been quantified by a generic metric that is applicable to drug screening. Prolonged repolarization leading to repolarization failure is highly arrhythmogenic. It may lead to ventricular tachycardia caused by triggered activity from early afterdepolarizations (EADs), or it may promote the occurrence of unidirectional conduction block and reentry. Both types of arrhythmia may deteriorate into ventricular fibrillation (VF) and death. We define the Repolarization Reserve Current (RRC) as the minimum constant current necessary to prevent normal repolarization of a cell. After developing and testing RRC for nine computational ionic models of various species, we applied it experimentally to atrial and ventricular human induced pluripotent stem cell-derived cardiomyocyte (hiPSC-CM), and isolated guinea-pig ventricular cardiomyocytes. In simulations, repolarization was all-or-none with a precise, model-dependent critical RRC, resulting in a discrete shift in the Action Potential Duration (APD) - RRC relation, in the occurrence of EADs and repolarization failure. These data were faithfully reproduced in cellular experiments. RRC allows simple, fast, unambiguous quantification of the arrhythmogenic propensity in cardiac cells of various origins and species without the need of prior knowledge of underlying currents and is suitable for high throughput applications, and personalized medicine applications.

Propagation Failure by TRPM4 Overexpression.

Transient receptor potential melastatin member 4 (TRPM4) channels are nonselective monovalent cationic channels found in human atria and conduction system. Overexpression of TRPM4 channels has been found in families suffering from inherited cardiac arrhythmias, notably heart block. In this study, we integrate a mathematical formulation of the TRPM4 channel into a Purkinje cell model (Pan-Rudy model). Instead of simply adding the channel to the model, a combination of existing currents equivalent to the TRPM4 current was constructed, based on TRPM4 current dynamics. The equivalent current was then replaced by the TRPM4 current to preserve the model action potential. Single-cell behavior showed early afterdepolarizations for increases in TRPM4 channel expression above twofold. In a homogeneous strand of tissue, propagation conducted faithfully for lower expression levels but failed completely for more than a doubling of TRPM4 channel expression. Only with a heterogeneous distribution of channel expression was intermittent heart block seen. This study suggests that in Purkinje fibers, TRPM4 channels may account for sodium background current (INab), and that a heterogeneous expression of TRPM4 channels in the His/Purkinje system is required for type II heart block, as seen clinically.

Disrupted junctional membrane complexes and hyperactive ryanodine receptors after acute junctophilin knockdown in mice.

BACKGROUND: Excitation-contraction coupling in striated muscle requires proper communication of plasmalemmal voltage-activated Ca2+ channels and Ca2+ release channels on sarcoplasmic reticulum within junctional membrane complexes. Although previous studies revealed a loss of junctional membrane complexes and embryonic lethality in germ-line junctophilin-2 (JPH2) knockout mice, it has remained unclear whether JPH2 plays an essential role in junctional membrane complex formation and the Ca(2+)-induced Ca(2+) release process in the heart. Our recent work demonstrated loss-of-function mutations in JPH2 in patients with hypertrophic cardiomyopathy. METHODS AND RESULTS: To elucidate the role of JPH2 in the heart, we developed a novel approach to conditionally reduce JPH2 protein levels using RNA interference. Cardiac-specific JPH2 knockdown resulted in impaired cardiac contractility, which caused heart failure and increased mortality. JPH2 deficiency resulted in loss of excitation-contraction coupling gain, precipitated by a reduction in the number of junctional membrane complexes and increased variability in the plasmalemma-sarcoplasmic reticulum distance. CONCLUSIONS: Loss of JPH2 had profound effects on Ca2+ release channel inactivation, suggesting a novel functional role for JPH2 in regulating intracellular Ca2+ release channels in cardiac myocytes. Thus, our novel approach of cardiac-specific short hairpin RNA-mediated knockdown of junctophilin-2 has uncovered a critical role for junctophilin in intracellular Ca2+ release in the heart.

Multiscale modeling of calcium cycling in cardiac ventricular myocyte: macroscopic consequences of microscopic dyadic function.

In cardiac ventricular myocytes, calcium (Ca) release occurs at distinct structures (dyads) along t-tubules, where L-type Ca channels (LCCs) appose sarcoplasmic reticulum (SR) Ca release channels (RyR2s). We developed a model of the cardiac ventricular myocyte that simulates local stochastic Ca release processes. At the local Ca release level, the model reproduces Ca spark properties. At the whole-cell level, the model reproduces the action potential, Ca currents, and Ca transients. Changes in microscopic dyadic properties (e.g., during detubulation in heart failure) affect whole-cell behavior in complex ways, which we investigated by simulating changes in the dyadic volume and number of LCCs/RyR2s in the dyad, and effects of calsequestrin (CSQN) as a Ca buffer (CSQN buffer) or a luminal Ca sensor (CSQN regulator). We obtained the following results: 1), Increased dyadic volume and reduced LCCs/RyR2s decrease excitation-contraction coupling gain and cause asynchrony of SR Ca release, and interdyad coupling partially compensates for the reduced synchrony. 2), Impaired CSQN buffer depresses Ca transients without affecting the synchrony of SR Ca release. 3), When CSQN regulator function is impaired, interdyad coupling augments diastolic Ca release activity to form Ca waves and long-lasting Ca release events.

Contributions of ion channel currents to ventricular action potential changes and induction of early afterdepolarizations during acute hypoxia.

RATIONALE: Variability in delivery of oxygen can lead to electric instability in the myocardium and the generation of arrhythmias. In addition ischemic heart disease and angina are associated with an increase in circulating catecholamines that further increases the risk of developing ventricular tachyarrhythmias. OBJECTIVE: We investigated the net effects of acute hypoxia and catecholamines on the cardiac action potential. METHODS AND RESULTS: We incorporated all published data on the effects of hypoxia on the late Na(+) current (I(Na-L)), the fast Na(+) current (I(Na)), the basal L-type Ca(2+) channel current (I(Ca-L)), and the slow (I(Ks)) and rapid components of the delayed rectifier K(+)-current (I(Kr)) in the absence and presence of beta-adrenergic receptor (beta-AR) stimulation into the Luo-Rudy model of the action potential. Hypoxia alone had little effect on the action potential configuration or action potential duration. However in the presence of beta-AR stimulation, hypoxia caused a prolongation of the action potential and early afterdepolarizations (EADs) and spontaneous tachycardia were induced. Experiments performed in guinea pig ventricular myocytes confirmed the modeling results. CONCLUSIONS: EADs occur predominantly because of the increased sensitivity of I(Ca-L) to beta-AR stimulation during hypoxia. beta-AR stimulation is necessary to induce EADs as EADs are never observed during hypoxia in the absence of beta-AR stimulation.