Ionic Mechanisms of Propagated Repolarization in a One-Dimensional Strand of Human Ventricular Myocyte Model.

Although repolarization has been suggested to propagate in cardiac tissue both theoretically and experimentally, it has been challenging to estimate how and to what extent the propagation of repolarization contributes to relaxation because repolarization only occurs in the course of membrane excitation in normal hearts. We established a mathematical model of a 1D strand of 600 myocytes stabilized at an equilibrium potential near the plateau potential level by introducing a sustained component of the late sodium current (INaL). By applying a hyperpolarizing stimulus to a small part of the strand, we succeeded in inducing repolarization which propagated along the strand at a velocity of 1~2 cm/s. The ionic mechanisms responsible for repolarization at the myocyte level, i.e., the deactivation of both the INaL and the L-type calcium current (ICaL), and the activation of the rapid component of delayed rectifier potassium current (IKr) and the inward rectifier potassium channel (IK1), were found to be important for the propagation of repolarization in the myocyte strand. Using an analogy with progressive activation of the sodium current (INa) in the propagation of excitation, regenerative activation of the predominant magnitude of IK1 makes the myocytes at the wave front start repolarization in succession through the electrical coupling via gap junction channels.

Local Control Model of a Human Ventricular Myocyte: An Exploration of Frequency-Dependent Changes and Calcium Sparks.

Calcium (Ca2+) sparks are the elementary events of excitation-contraction coupling, yet they are not explicitly represented in human ventricular myocyte models. A stochastic ventricular cardiomyocyte human model that adapts to intracellular Ca2+ ([Ca2+]i) dynamics, spark regulation, and frequency-dependent changes in the form of locally controlled Ca2+ release was developed. The 20,000 CRUs in this model are composed of 9 individual LCCs and 49 RyRs that function as couplons. The simulated action potential duration at 1 Hz steady-state pacing is ~0.280 s similar to human ventricular cell recordings. Rate-dependence experiments reveal that APD shortening mechanisms are largely contributed by the L-type calcium channel inactivation, RyR open fraction, and [Ca2+]myo concentrations. The dynamic slow-rapid-slow pacing protocol shows that RyR open probability during high pacing frequency (2.5 Hz) switches to an adapted "nonconducting" form of Ca2+-dependent transition state. The predicted force was also observed to be increased in high pacing, but the SR Ca2+ fractional release was lower due to the smaller difference between diastolic and systolic [Ca2+]SR. Restitution analysis through the S1S2 protocol and increased LCC Ca2+-dependent activation rate show that the duration of LCC opening helps modulate its effects on the APD restitution at different diastolic intervals. Ultimately, a longer duration of calcium sparks was observed in relation to the SR Ca2+ loading at high pacing rates. Overall, this study demonstrates the spontaneous Ca2+ release events and ion channel responses throughout various stimuli.

Activation of small conductance Ca2+ -activated K+ channels suppresses Ca2+ transient and action potential alternans in ventricular myocytes.

At the cellular level, cardiac alternans is observed as beat-to-beat alternations in contraction strength, action potential (AP) morphology and Ca2+ transient (CaT) amplitude, and is a risk factor for cardiac arrhythmia. The (patho)physiological roles of small conductance Ca2+ -activated K+ (SK) channels in ventricles are poorly understood. We tested the hypothesis that in single rabbit ventricular myocytes pharmacological modulation of SK channels plays a causative role for the development of pacing-induced CaT and AP duration (APD) alternans. SK channel blockers (apamin, UCL1684) had only a minor effect on AP repolarization. However, SK channel activation by NS309 resulted in significant APD shortening, demonstrating that functional SK channels are well expressed in ventricular myocytes. The effects of NS309 were prevented or reversed by apamin and UCL1684, indicating that NS309 acted on SK channels. SK channel activation abolished or reduced the degree of pacing-induced CaT and APD alternans. Inhibition of KV 7.1 (with HMR1556) and KV 11.1 (with E4031) channels was used to mimic conditions of long QT syndromes type-1 and type-2, respectively. Both HMR1556 and E4031 enhanced CaT alternans that was prevented by SK channel activation. In AP voltage-clamped cells the SK channel activator had no effect on CaT alternans, confirming that suppression of CaT alternans was caused by APD shortening. APD shortening contributed to protection from alternans by lowering sarcoplasmic reticulum Ca2+ content and curtailing Ca2+ release. The data suggest that SK activation could be a potential intervention to avert development of alternans with important ramifications for arrhythmia prevention and therapy for patients with long QT syndrome. KEY POINTS: At the cellular level, cardiac alternans is observed as beat-to-beat alternations in contraction strength, action potential (AP) morphology and intracellular Ca2+ release amplitude, and is a risk factor for cardiac arrhythmia. The (patho)physiological roles of small conductance Ca2+ -activated K+ (SK) channels in ventricles are poorly understood. We investigated whether pharmacological modulation of SK channels affects the development of cardiac alternans in normal ventricular cells and in cells with drug-induced long QT syndrome (LQTS). While SK channel blockers have only a minor effect on AP morphology, their activation leads to AP shortening and abolishes or reduces the degree of pacing-induced Ca2+ and AP alternans. AP shortening contributed to protection against alternans by lowering sarcoplasmic reticulum Ca2+ content and curtailing Ca2+ release. The data suggest SK activation as a potential intervention to avert the development of alternans with important ramifications for arrhythmia prevention for patients with LQTS.

The Role of Ca2+ Sparks in Force Frequency Relationships in Guinea Pig Ventricular Myocytes.

Calcium sparks are the elementary Ca2+ release events in excitation-contraction coupling that underlie the Ca2+ transient. The frequency-dependent contractile force generated by cardiac myocytes depends upon the characteristics of the Ca2+ transients. A stochastic computational local control model of a guinea pig ventricular cardiomyocyte was developed, to gain insight into mechanisms of force-frequency relationship (FFR). This required the creation of a new three-state RyR2 model that reproduced the adaptive behavior of RyR2, in which the RyR2 channels transition into a different state when exposed to prolonged elevated subspace [Ca2+]. The model simulations agree with previous experimental and modeling studies on interval-force relations. Unlike previous common pool models, this local control model displayed stable action potential trains at 7 Hz. The duration and the amplitude of the [Ca2+]myo transients increase in pacing rates consistent with the experiments. The [Ca2+]myo transient reaches its peak value at 4 Hz and decreases afterward, consistent with experimental force-frequency curves. The model predicts, in agreement with previous modeling studies of Jafri and co-workers, diastolic sarcoplasmic reticulum, [Ca2+]sr, and RyR2 adaptation increase with the increased stimulation frequency, producing rising, rather than falling, amplitude of the myoplasmic [Ca2+] transients. However, the local control model also suggests that the reduction of the L-type Ca2+ current, with an increase in pacing frequency due to Ca2+-dependent inactivation, also plays a role in the negative slope of the FFR. In the simulations, the peak Ca2+ transient in the FFR correlated with the highest numbers of SR Ca2+ sparks: the larger average amplitudes of those sparks, and the longer duration of the Ca2+ sparks.

Cardiac TRPV4 channels.

Transient Receptor Potential Vanilloid 4 (TRPV4) is expressed in numerous cell types within the heart, yet the expression levels, subcellular localization, and functional relevance of TRPV4 in cardiac myocytes is under-appreciated. Recent data indicate a critical role of TRPV4 in both atrial and ventricular myocyte biology, with expression levels and channel function increasing following pathological scenarios including ischemia, myocardial infarction, mechanical stress, and inflammation. Excessive activation of TRPV4 at the cellular level contributes to enhanced Ca2+ entry which predisposes the cardiac myocyte to pro-arrhythmic Ca2+ overload and electrophysiological abnormalities. At the organ level, excessive TRPV4 activity associates with cardiac hypercontractility, cardiac damage, ventricular arrhythmia, and atrial fibrillation. This manuscript chapter describes the emerging literature on TRPV4 in cardiac myocytes in physiology and disease.

Quantitative mapping of force-pCa curves to whole-heart contraction and relaxation.

The force-pCa (F-pCa) curve is used to characterize steady-state contractile properties of cardiac muscle cells in different physiological, pathological and pharmacological conditions. This provides a reduced preparation in which to isolate sarcomere mechanisms. However, it is unclear how changes in the F-pCa curve impact emergent whole-heart mechanics quantitatively. We study the link between sarcomere and whole-heart function using a multiscale mathematical model of rat biventricular mechanics that describes sarcomere, tissue, anatomy, preload and afterload properties quantitatively. We first map individual cell-level changes in sarcomere-regulating parameters to organ-level changes in the left ventricular function described by pressure-volume loop characteristics (e.g. end-diastolic and end-systolic volumes, ejection fraction and isovolumetric relaxation time). We next map changes in the sarcomere-regulating parameters to changes in the F-pCa curve. We demonstrate that a change in the F-pCa curve can be caused by multiple different changes in sarcomere properties. We demonstrate that changes in sarcomere properties cause non-linear and, importantly, non-monotonic changes in left ventricular function. As a result, a change in sarcomere properties yielding changes in the F-pCa curve that improve contractility does not guarantee an improvement in whole-heart function. Likewise, a desired change in whole-heart function (i.e. ejection fraction or relaxation time) is not caused by a unique shift in the F-pCa curve. Changes in the F-pCa curve alone cannot be used to predict the impact of a compound on whole-heart function. KEY POINTS: The force-pCa (F-pCa) curve is used to assess myofilament calcium sensitivity after pharmacological modulation and to infer pharmacological effects on whole-heart function. We demonstrate that there is a non-unique mapping from changes in F-pCa curves to changes in left ventricular (LV) function. The effect of changes in F-pCa on LV function depend on the state of the heart and could be different for different pathological conditions. Screening of compounds to impact whole-heart function by F-pCa should be combined with active tension and calcium transient measurements to predict better how changes in muscle function will impact whole-heart physiology.

Critical Requirements for the Initiation of a Cardiac Arrhythmia in Rat Ventricle: How Many Myocytes?

Cardiovascular disease is the leading cause of death worldwide due in a large part to arrhythmia. In order to understand how calcium dynamics play a role in arrhythmogenesis, normal and dysfunctional Ca2+ signaling in a subcellular, cellular, and tissued level is examined using cardiac ventricular myocytes at a high temporal and spatial resolution using multiscale computational modeling. Ca2+ sparks underlie normal excitation-contraction coupling. However, under pathological conditions, Ca2+ sparks can combine to form Ca2+ waves. These propagating elevations of (Ca2+)i can activate an inward Na+-Ca2+ exchanger current (INCX) that contributes to early after-depolarization (EADs) and delayed after-depolarizations (DADs). However, how cellular currents lead to full depolarization of the myocardium and how they initiate extra systoles is still not fully understood. This study explores how many myocytes must be entrained to initiate arrhythmogenic depolarizations in biophysically detailed computational models. The model presented here suggests that only a small number of myocytes must activate in order to trigger an arrhythmogenic propagating action potential. These conditions were examined in 1-D, 2-D, and 3-D considering heart geometry. The depolarization of only a few hundred ventricular myocytes is required to trigger an ectopic depolarization. The number decreases under disease conditions such as heart failure. Furthermore, in geometrically restricted parts of the heart such as the thin muscle strands found in the trabeculae and papillary muscle, the number of cells needed to trigger a propagating depolarization falls even further to less than ten myocytes.

Multiscale Modeling of the Mitochondrial Origin of Cardiac Reentrant and Fibrillatory Arrhythmias.

While mitochondrial dysfunction has been implicated in the pathogenesis of cardiac arrhythmias, how the abnormality occurring at the organelle level escalates to influence the rhythm of the heart remains incompletely understood. This is due, in part, to the complexity of the interactions formed by cardiac electrical, mechanical, and metabolic subsystems at various spatiotemporal scales that is difficult to fully comprehend solely with experiments. Computational models have emerged as a powerful tool to explore complicated and highly dynamic biological systems such as the heart, alone or in combination with experimental measurements. Here, we describe a strategy of integrating computer simulations with optical mapping of cardiomyocyte monolayers to examine how regional mitochondrial dysfunction elicits abnormal electrical activity, such as rebound and spiral waves, leading to reentry and fibrillation in cardiac tissue. We anticipate that this advanced modeling technology will enable new insights into the mechanisms by which changes in subcellular organelles can impact organ function.

Acute Detubulation of Ventricular Myocytes Amplifies the Inhibitory Effect of Cholinergic Agonist on Intracellular Ca2+ Transients.

Muscarinic receptors expressed in cardiac myocytes play a critical role in the regulation of heart function by the parasympathetic nervous system. How the structural organization of cardiac myocytes affects the regulation of Ca2+ handling by muscarinic receptors is not well-defined. Using confocal Ca2+ imaging, patch-clamp techniques, and immunocytochemistry, the relationship between t-tubule density and cholinergic regulation of intracellular Ca2+ in normal murine ventricular myocytes and myocytes with acute disruption of the t-tubule system caused by formamide treatment was studied. The inhibitory effect of muscarinic receptor agonist carbachol (CCh, 10 µM) on the amplitude of Ca2+ transients, evoked by field-stimulation in the presence of 100 nM isoproterenol (Iso), a ß-adrenergic agonist, was directly proportional to the level of myocyte detubulation. The timing of the maximal rate of fluorescence increase of fluo-4, a Ca2+-sensitive dye, was used to classify image pixels into the regions functionally coupled or uncoupled to the sarcolemmal Ca2+ influx (ICa). CCh decreased the fraction of coupled regions and suppressed Ca2+ propagation from sarcolemma inside the cell. Formamide treatment reduced ICa density and decreased sarcoplasmic reticulum (SR) Ca2+ content. CCh did not change SR Ca2+ content in Iso-stimulated control and formamide-treated myocytes. CCh inhibited peak ICa recorded in the presence of Iso by ∼20% in both the control and detubulated myocytes. Reducing ICa amplitude up to 40% by changing the voltage step levels from 0 to -25 mV decreased Ca2+ transients in formamide-treated but not in control myocytes in the presence of Iso. CCh inhibited CaMKII activity, whereas CaMKII inhibition with KN93 mimicked the effect of CCh on Ca2+ transients in formamide-treated myocytes. It was concluded that the downregulation of t-tubules coupled with the diminished efficiency of excitation-contraction coupling, increases the sensitivity of Ca2+ release and propagation to muscarinic receptor-mediated inhibition of both ICa and CaMKII activity.

Inhibitory effects of aloperine on voltage-gated Na+ channels in rat ventricular myocytes.

Aloperine (ALO), a quinolizidine alkaloid extracted from Sophora alopecuroides L., modulates hypertension, ventricular remodeling, and myocardial ischemia. However, few studies have evaluated the effects of ALO on other cardiovascular parameters. Accordingly, in this study, we used a rat model of aconitine-induced ventricular arrhythmia to assess the effects of ALO. Notably, ALO pretreatment delayed the onset of ventricular premature and ventricular tachycardia and reduced the incidence of fatal ventricular fibrillation. Moreover, whole-cell patch-clamp assays in rats' ventricular myocyte showed that ALO (3, 10, and 30 µM) significantly reduced the peak sodium current density of voltage-gated Na+ channel currents (INa) in a concentration-dependent manner. The gating kinetics characteristics showed that the steady-state activation and recovery curve were shifted in positive direction along the voltage axis, respectively, and the steady-state inactivation curve was shifted in negative direction along the voltage axis, i.e., which was similar to the inhibitory effects of amiodarone. These results indicated that ALO had anti-arrhythmic effects, partly attributed to INa inhibition. ALO may act as a class I sodium channel anti-arrhythmia agent.

Compartmentation of ß2 -adrenoceptor stimulated cAMP responses by phosphodiesterase types 2 and 3 in cardiac ventricular myocytes.

BACKGROUND AND PURPOSE: In cardiac myocytes, cyclic AMP (cAMP) produced by both ß1 - and ß2 -adrenoceptors increases L-type Ca2+ channel activity and myocyte contraction. However, only cAMP produced by ß1 -adrenoceptors enhances myocyte relaxation through phospholamban-dependent regulation of the sarco/endoplasmic reticulum Ca2+ ATPase 2 (SERCA2). Here we have tested the hypothesis that stimulation of ß2 -adrenoceptors produces a cAMP signal that is unable to reach SERCA2 and determine what role, if any, phosphodiesterase (PDE) activity plays in this compartmentation. EXPERIMENTAL APPROACH: The cAMP responses produced by ß1 -and ß2 -adrenoceptor stimulation were studied in adult rat ventricular myocytes using two different fluorescence resonance energy transfer (FRET)-based biosensors, the Epac2-camps, which is expressed uniformly throughout the cytoplasm of the entire cell and the Epac2-αKAP, which is targeted to the SERCA2 signalling complex. KEY RESULTS: Selective activation of ß1 - or ß2 -adrenoceptors produced cAMP responses detected by Epac2-camps. However, only stimulation of ß1 -adrenoceptors produced a cAMP response detected by Epac2-αKAP. Yet, stimulation of ß2 -adrenoceptors was able to produce a cAMP signal detected by Epac2-αKAP in the presence of selective inhibitors of PDE2 or PDE3, but not PDE4. CONCLUSION AND IMPLICATIONS: These results support the conclusion that cAMP produced by ß2 -adrenoceptor stimulation was not able to reach subcellular locations where the SERCA2 pump is located. Furthermore, this compartmentalized response is due at least in part to PDE2 and PDE3 activity. This discovery could lead to novel PDE-based therapeutic treatments aimed at correcting cardiac relaxation defects associated with certain forms of heart failure.

Investigation of the Effects of the Short QT Syndrome D172N Kir2.1 Mutation on Ventricular Action Potential Profile Using Dynamic Clamp.

The congenital short QT syndrome (SQTS) is a cardiac condition that leads to abbreviated ventricular repolarization and an increased susceptibility to arrhythmia and sudden death. The SQT3 form of the syndrome is due to mutations to the KCNJ2 gene that encodes Kir2.1, a critical component of channels underlying cardiac inwardly rectifying K+ current, IK1. The first reported SQT3 KCNJ2 mutation gives rise to the D172N Kir2.1 mutation, the consequences of which have been studied on recombinant channels in vitro and in ventricular cell and tissue simulations. The aim of this study was to establish the effects of the D172N mutation on ventricular repolarization through real-time replacement of IK1 using the dynamic clamp technique. Whole-cell patch-clamp recordings were made from adult guinea-pig left ventricular myocytes at physiological temperature. Action potentials (APs) were elicited at 1 Hz. Intrinsic IK1 was inhibited with a low concentration (50 µM) of Ba2+ ions, which led to AP prolongation and triangulation, accompanied by a ∼6 mV depolarization of resting membrane potential. Application of synthetic IK1 through dynamic clamp restored AP duration, shape and resting potential. Replacement of wild-type (WT) IK1 with heterozygotic (WT-D172N) or homozygotic (D172N) mutant formulations under dynamic clamp significantly abbreviated AP duration (APD90) and accelerated maximal AP repolarization velocity, with no significant hyperpolarization of resting potential. Across stimulation frequencies from 0.5 to 3 Hz, the relationship between APD90 and cycle length was downward shifted, reflecting AP abbreviation at all stimulation frequencies tested. In further AP measurements at 1 Hz from hiPSC cardiomyocytes, the D172N mutation produced similar effects on APD and repolarization velocity; however, resting potential was moderately hyperpolarized by application of mutant IK1 to these cells. Overall, the results of this study support the major changes in ventricular cell AP repolarization with the D172N predicted from prior AP modelling and highlight the potential utility of using adult ventricular cardiomyocytes for dynamic clamp exploration of functional consequences of Kir2.1 mutations.

Reduced nNOS activity is responsible for impaired fatty acid-dependent mitochondrial oxygen consumption in atrial myocardium from hypertensive rat.

Fatty acid (FA)-dependent mitochondrial activities of atrial myocardium in hypertension (HTN) and its regulation by nitric oxide (NO) remain unidentified. Here, we have studied palmitic acid (PA) regulation of cardiac mitochondrial oxygen consumption rate (OCR) in left atrial (LA) myocardium of sham and angiotensin II-induced HTN rats and their regulations by endothelial NO synthase (eNOS) and neuronal NO synthase (nNOS). The effects were compared with those of left ventricular (LV) myocytes. Our results showed that OCR was greater in HTN-LA compared with that in sham-LA. PA increased OCR in sham-LA, sham-LV, and HTN-LV but reduced it in HTN-LA. Inhibition of nNOS (S-methyl-L-thiocitrulline, SMTC) or eNOS/nNOS (Nω-nitro-L-arginine methyl ester hydrochloride, L-NAME) reduced PA increment of OCR in sham-LA but exerted no effect on OCR in HTN-LA. SMTC reduced OCR in HTN-LV and L-NAME reduced OCR in sham-LV. nNOS was the predominant source of NO in LA and LV. nNOS-derived NO was increased in HTN-LA and HTN-LV. PA reduced eNOSSer1177, nNOSSer1417, and NO level in HTN-LA but exerted no effect in sham-LA. In contrast, PA increased NO in HTN-LV and enhanced nNOSSer1417 but reduced NO level in sham-LV without affecting eNOSSer1177, eNOSThr495, or nNOSSer1417. 2-Bromopalmitate (2BP), which blocks the S-palmitoylation of target proteins, prevented PA-dependent decrease of nNOSSer1417 and OCR in HTN-LA. In HTN-LV, 2BP prevented PA-induced OCR without affecting nNOSSer1417. Our results reveal that FA-induced mitochondrial activity in atrial myocardium is impaired in HTN which is mediated by reduced nNOS activity and NO bioavailability. Metabolic dysregulation may underlie diastolic dysfunction of atrial myocardium in HTN.

Quantification of model and data uncertainty in a network analysis of cardiac myocyte mechanosignalling.

Cardiac myocytes transduce changes in mechanical loading into cellular responses via interacting cell signalling pathways. We previously reported a logic-based ordinary differential equation model of the myocyte mechanosignalling network that correctly predicts 78% of independent experimental results not used to formulate the original model. Here, we use Monte Carlo and polynomial chaos expansion simulations to examine the effects of uncertainty in parameter values, model logic and experimental validation data on the assessed accuracy of that model. The prediction accuracy of the model was robust to parameter changes over a wide range being least sensitive to uncertainty in time constants and most affected by uncertainty in reaction weights. Quantifying epistemic uncertainty in the reaction logic of the model showed that while replacing 'OR' with 'AND' reactions greatly reduced model accuracy, replacing 'AND' with 'OR' reactions was more likely to maintain or even improve accuracy. Finally, data uncertainty had a modest effect on assessment of model accuracy. This article is part of the theme issue 'Uncertainty quantification in cardiac and cardiovascular modelling and simulation'.

Alterations of Ca2+ signaling and Ca2+ release sites in cultured ventricular myocytes with intact internal Ca2+ storage.

Although cultured adult cardiac myocytes in combination with cell-level genetic modifications have been adopted for the study of protein function, the cellular alterations caused by the culture conditions themselves need to be clarified before we can interpret the effects of genetically altered proteins. We systematically compared the cellular morphology, global Ca2+ signaling, elementary Ca2+ release (sparks), and arrangement of ryanodine receptor (RyR) clusters in short-term (2 days)-cultured adult rat ventricular myocytes with those of freshly isolated myocytes. The transverse (t)-tubules were remarkably decreased (to ∼25%) by culture, and whole-cell capacitance was reduced by ∼35%. The magnitude of depolarization-induced Ca2+ transients decreased to ∼50%, and Ca2+ transient decay was slowed by culture. The culture did not affect sarcoplasmic reticulum (SR) Ca2+ loading. Therefore, fractional Ca2+ release was attenuated by culture. In the cultured cells, the L-type Ca2+ current (ICa) was smaller (∼50% of controls) and its inactivation was slower. In cultured myocytes, there were significantly fewer (∼50% of control) Ca2+ sparks, the local Ca2+ releases through RyR clusters, compared with in freshly isolated cells. Amplitude and kinetics (duration and time-to-peak) of individual sparks were similar, but they showed greater width in cultured cells. Immunolocalization analysis revealed that the cross-striation of RyRs distribution became weaker and less organized, and that the density of RyR clusters decreased in cultured myocytes. Our data suggest that the loss of t-tubules and generation of compromised Ca2+ transients and ICa in short-term adult ventricular cell culture are independent of SR Ca2+ loading status. In addition, the deteriorated arrangement of the RyR-clusters and their decreased density after short-term culture may be partly responsible for fewer Ca2+ sparks and a decrease in global Ca2+ release.

Location and function of transient receptor potential canonical channel 1 in ventricular myocytes.

Transient receptor potential canonical 1 (TRPC1) protein is abundantly expressed in cardiomyocytes. While TRPC1 is supposed to be critically involved in cardiac hypertrophy, its physiological role in cardiomyocytes is poorly understood. We investigated the subcellular location of TRPC1 and its contribution to Ca2+ signaling in mammalian ventricular myocytes. Immunolabeling, three-dimensional scanning confocal microscopy and quantitative colocalization analysis revealed an abundant intracellular location of TRPC1 in neonatal rat ventricular myocytes (NRVMs) and adult rabbit ventricular myocytes. TRPC1 was colocalized with intracellular proteins including sarco/endoplasmic reticulum Ca2+ ATPase 2 in the sarcoplasmic reticulum (SR). Colocalization with wheat germ agglutinin, which labels the glycocalyx and thus marks the sarcolemma including the transverse tubular system, was low. Super-resolution and immunoelectron microscopy supported the intracellular location of TRPC1. We investigated Ca2+ signaling in NRVMs after adenoviral TRPC1 overexpression or silencing. In NRVMs bathed in Na+ and Ca2+ free solution, TRPC1 overexpression and silencing was associated with a decreased and increased SR Ca2+ content, respectively. In isolated rabbit cardiomyocytes bathed in Na+ and Ca2+ free solution, we found an increased decay of the cytosolic Ca2+ concentration [Ca2+]i and increased SR Ca2+ content in the presence of the TRPC channel blocker SKF-96365. In a computational model of rabbit ventricular myocytes at physiological pacing rates, Ca2+ leak through SR TRPC channels increased the systolic and diastolic [Ca2+]i with only minor effects on the action potential and SR Ca2+ content. Our studies suggest that TRPC1 channels are localized in the SR, and not present in the sarcolemma of ventricular myocytes. The studies provide evidence for a role of TRPC1 as a contributor to SR Ca2+ leak in cardiomyocytes, which was previously explained by ryanodine receptors only. We propose that the findings will guide us to an understanding of TRPC1 channels as modulators of [Ca2+]i and contractility in cardiomyocytes.

Development, calibration, and validation of a novel human ventricular myocyte model in health, disease, and drug block.

Human-based modelling and simulations are becoming ubiquitous in biomedical science due to their ability to augment experimental and clinical investigations. Cardiac electrophysiology is one of the most advanced areas, with cardiac modelling and simulation being considered for virtual testing of pharmacological therapies and medical devices. Current models present inconsistencies with experimental data, which limit further progress. In this study, we present the design, development, calibration and independent validation of a human-based ventricular model (ToR-ORd) for simulations of electrophysiology and excitation-contraction coupling, from ionic to whole-organ dynamics, including the electrocardiogram. Validation based on substantial multiscale simulations supports the credibility of the ToR-ORd model under healthy and key disease conditions, as well as drug blockade. In addition, the process uncovers new theoretical insights into the biophysical properties of the L-type calcium current, which are critical for sodium and calcium dynamics. These insights enable the reformulation of L-type calcium current, as well as replacement of the hERG current model.

Decades of intensive experimental and clinical research have revealed much about how the human heart works. Though incomplete, this knowledge has been used to construct computer models that represent the activity of this organ as a whole, and of its individual chambers (the atria and ventricles), tissues and cells. Such models have been used to better understand life-threatening irregular heartbeats; they are also beginning to be used to guide decisions about the treatment of patients and the development of new drugs by the pharmaceutical industry. Yet existing computer models of the electrical activity of the human heart are sometimes inconsistent with experimental data. This problem led Tomek et al. to try to create a new model that was consistent with established biophysical knowledge and experimental data for a wide range of conditions including disease and drug action. Tomek et al. designed a strategy that explicitly separated the construction and validation of a model that could recreate the electrical activity of the ventricles in a human heart. This model was able to integrate and explain a wide range of properties of both healthy and diseased hearts, including their response to different drugs. The development of the model also uncovered and resolved theoretical inconsistencies that have been present in almost all models of the heart from the last 25 years. Tomek et al. hope that their new human heart model will enable more basic, translational and clinical research into a range of heart diseases and accelerate the development of new therapies.

BK channels regulate calcium oscillations in ventricular myocytes on different substrate stiffness.

Substrate stiffness is essential for cell functions, but the mechanisms by which cell sense mechanical cues are still unclear. Here we show that the frequency and the amplitude of spontaneous Ca2+ oscillations were greater in chick cardiomyocytes cultured on the stiff substrates than that on the soft substrates. The spontaneous Ca2+ oscillations were increased on stiff substrates. However, an eliminated dependence of the Ca2+ oscillations on substrate stiffness was observed after applying blocker of the large-conductance Ca2+-activated K+ (BK) channels. In addition, the activity of BK channels in cardiomyocytes cultured on the stiff substrates was decreased. These results provide compelling evidences to show that BK channels are crucial in substrate stiffness-dependent regulation of the Ca2+ oscillation in cardiomyocytes.

Prediction of the mechanical response of cardiac alternans by using an electromechanical model of human ventricular myocytes.

PURPOSE: Although the quantitative analysis of electromechanical alternans is important, previous studies have focused on electrical alternans, and there is a lack quantitative analysis of mechanical alternans at the subcellular level according to various basic cycle lengths (BCLs). Therefore, we used the excitation-contraction (E-C) coupling model of human ventricular cells to quantitatively analyze the mechanical alternans of ventricular cells according to various BCLs. METHODS: To implement E-C coupling, we used calcium transient data, which is the output data of electrical simulation using the electrophysiological model of human ventricular myocytes, as the input data of mechanical simulation using the contractile myofilament dynamics model. Moreover, we applied various loads on ventricular cells for implementation of isotonic and isometric contraction. RESULTS: As the BCL was reduced from 1000 to 200 ms at 30 ms increments, mechanical alternans, as well as electrical alternans, were observed. At this time, the myocardial diastolic tension increased, and the contractile ATP consumption rate remained greater than zero even in the resting state. Furthermore, the time of peak tension, equivalent cell length, and contractile ATP consumption rate were all reduced. There are two tendencies that endocardial, mid-myocardial, and epicardial cells have the maximum amplitude of tension and the peak systolic tension begins to appear at a high rate under the isometric condition at a particular BCL. CONCLUSIONS: We observed mechanical alternans of ventricular myocytes as well as electrical alternans, and identified unstable conditions associated with mechanical alternans. We also determined the amount of BCL given to each ventricular cell to generate stable and high tension state in the case of isometric contraction.

Role of spatial dispersion of repolarization in reentry around a functional core versus reentry around a fixed anatomical core.

INTRODUCTION: Successful initiation of spiral wave reentry in the neonatal rat ventricular myocyte (NRVM) monolayer implicitly assumes the presence of spatial dispersion of repolarization (DR), which is difficult to quantify. We recently introduced a NRVM monolayer that utilizes anthopleurin-A to impart a prolonged plateau to the NRVM action potential. This was associated with a significant degree of spatial DR that lends itself to accurate quantification. METHODS AND RESULTS: We utilized the monolayer and fluorescence optical mapping of intracellular calcium transients (FCai ) to systematically study and compare the contribution of spatial dispersion of the duration of FCai (as a surrogate of DR) to induction of spiral wave reentry around a functional core versus reentry around a fixed anatomical obstacle. We show that functional reentry could be initiated by a premature stimulus acting on a substrate of spatial DR resulting in a functional line of propagation block. Subsequent wave fronts circulated around a central core of functional obstacle created by sustained depolarization from the circulating wave front. Both initiation and termination of spiral wave reentry around an anatomical obstacle consistently required participation of a region of functional propagation block. This region was similarly based on spatial DR. Spontaneous termination of spiral wave reentry also resulted from block in the functional component of the circuit obstacle, usually preceded by beat-to-beat slowing of propagation. CONCLUSIONS: The study demonstrates the critical contribution of DR to spiral wave reentry around a purely functional core as well as reentry around a fixed anatomical core.