Modelling sarcoplasmic reticulum calcium ATPase and its regulation in cardiac myocytes

A detailed mathematical model of the human atrial cardiomyocyte: integration of electrophysiology and cardiomechanics

Mechano-electric regulations (MER) play an important role in the maintenance of cardiac performance. Mechano-calcium and mechano-electric feedback (MCF and MEF) pathways adjust the cardiomyocyte contractile force according to mechanical perturbations and affects electro-mechanical coupling. MER integrates all these regulations in one unit resulting in a complex phenomenon. Computational modelling is a useful tool to accelerate the mechanistic understanding of complex experimental phenomena. We have developed a novel model that integrates the MER loop for human atrial cardiomyocytes with proper consideration of feedforward and feedback pathways. The model couples a modified version of the action potential (AP) Koivumäki model with the contraction model by Quarteroni group. The model simulates iso-sarcometric and isometric twitches and the feedback effects on AP and Ca2+-handling. The model showed a biphasic response of Ca2+ transient (CaT) peak to increasing pacing rates and highlights the possible mechanisms involved. The model has shown a shift of the threshold for AP and CaT alternans from 4.6 to 4 Hz under post-operative atrial fibrillation, induced by depressed SERCA activity. The alternans incidence was dependent on a chain of mechanisms including RyRs availability time, MCF coupling, CaMKII phosphorylation, and the stretch levels. As a result, the model predicted a 10% slowdown of conduction velocity for a 20% stretch, suggesting a role of stretch in creation of substrate formation for atrial fibrillation. Overall, we conclude that the developed model provides a physiological CaT followed by a physiological twitch. This model can open pathways for the future studies of human atrial electromechanics. KEY POINTS: With the availability of human atrial cellular data, interest in atrial-specific model integration has been enhanced. We have developed a detailed mathematical model of human atrial cardiomyocytes including the mechano-electric regulatory loop. The model has gone through calibration and evaluation phases against a wide collection of available human in-vitro data. The usefulness of the model for analysing clinical problems has been preliminaryly tested by simulating the increased incidence of Ca2+ transient and action potential alternans at high rates in post-operative atrial fibrillation condition. The model determines the possible role of mechano-electric feedback in alternans incidence, which can increase vulnerability to atrial arrhythmias by varying stretch levels. We found that our physiologically accurate description of Ca2+ handling can reproduce many experimental phenomena and can help to gain insights into the underlying pathophysiological mechanisms.

Longitudinal diffusion barriers imposed by myofilaments and mitochondria in murine cardiac myocytes

Using optical and electrical methods, we document that diffusion in the cytoplasm of BL6 murine cardiomyocytes becomes restricted >20-fold as molecular weight increases from 30 to 2,000, roughly as expected for pores with porin channel dimensions. Bodipy-FL ATP diffuses >40-fold slower than in free water at 25°C. From several fluorophores analyzed, bound fluorophore fractions range from 0.1 for a 2 kD FITC-labeled polyethylene glycol to 0.93 for sulforhodamine. Unbound fluorophores diffuse at 0.5-8 × 10-7 cm2/s (5-80 μm2/s). Analysis of Na/K pump and veratridine-modified Na channel currents suggests that Na diffusion is nearly unrestricted at 35°C (time constant for equilibration with the pipette tip, ∼20 s). Using multiple strategies, we estimate that at 35°C, ATP diffuses four to eight times slower than in free water. To address whether restrictions are caused more by protein or membrane networks, we verified first that a protein gel, 10 g% gelatin, restricts diffusion with strong dependence on molecular weight. Solute diffusion in membrane-extracted cardiac myofilaments, confined laterally by suction into large-diameter pipette tips, is less restricted than in intact myocytes. Notably, myofilaments extracted similarly from skeletal (diaphragm) myocytes are less restrictive. Solute diffusion in myocytes with sarcolemma permeabilized by β-escin (80 µM) is similar to diffusion in intact myocytes. Restrictions are strain-dependent, being twofold greater in BL6 myocytes than in CD1/J6/129svJ myocytes. Furthermore, longitudinal diffusion is 2.5-fold more restricted in CD1/J6/129svJ myocytes lacking the mitochondrial porin, VDAC1, than in WT CD1/J6/129svJ myocytes. Thus, mitochondria networks restrict long-range diffusion while presumably optimizing nucleotide transfer between myofilaments and mitochondria. We project that diffusion restrictions imposed by both myofilaments and the outer mitochondrial membrane are important determinants of total free cytoplasmic AMP and ADP (∼10 μM). However, the capacity of diffusion to deliver ATP to myofilaments remains ∼100-fold greater than ATP consumption.

Physiology of intracellular calcium buffering

Calcium signaling underlies much of physiology. Almost all the Ca2+ in the cytoplasm is bound to buffers, with typically only ∼1% being freely ionized at resting levels in most cells. Physiological Ca2+ buffers include small molecules and proteins, and experimentally Ca2+ indicators will also buffer calcium. The chemistry of interactions between Ca2+ and buffers determines the extent and speed of Ca2+ binding. The physiological effects of Ca2+ buffers are determined by the kinetics with which they bind Ca2+ and their mobility within the cell. The degree of buffering depends on factors such as the affinity for Ca2+, the Ca2+ concentration, and whether Ca2+ ions bind cooperatively. Buffering affects both the amplitude and time course of cytoplasmic Ca2+ signals as well as changes of Ca2+ concentration in organelles. It can also facilitate Ca2+ diffusion inside the cell. Ca2+ buffering affects synaptic transmission, muscle contraction, Ca2+ transport across epithelia, and the killing of bacteria. Saturation of buffers leads to synaptic facilitation and tetanic contraction in skeletal muscle and may play a role in inotropy in the heart. This review focuses on the link between buffer chemistry and function and how Ca2+ buffering affects normal physiology and the consequences of changes in disease. As well as summarizing what is known, we point out the many areas where further work is required.

Fibrotic Remodeling during Persistent Atrial Fibrillation: In Silico Investigation of the Role of Calcium for Human Atrial Myofibroblast Electrophysiology

During atrial fibrillation, cardiac tissue undergoes different remodeling processes at different scales from the molecular level to the tissue level. One central player that contributes to both electrical and structural remodeling is the myofibroblast. Based on recent experimental evidence on myofibroblasts' ability to contract, we extended a biophysical myofibroblast model with Ca2+ handling components and studied the effect on cellular and tissue electrophysiology. Using genetic algorithms, we fitted the myofibroblast model parameters to the existing in vitro data. In silico experiments showed that Ca2+ currents can explain the experimentally observed variability regarding the myofibroblast resting membrane potential. The presence of an L-type Ca2+ current can trigger automaticity in the myofibroblast with a cycle length of 799.9 ms. Myocyte action potentials were prolonged when coupled to myofibroblasts with Ca2+ handling machinery. Different spatial myofibroblast distribution patterns increased the vulnerable window to induce arrhythmia from 12 ms in non-fibrotic tissue to 22 ± 2.5 ms and altered the reentry dynamics. Our findings suggest that Ca2+ handling can considerably affect myofibroblast electrophysiology and alter the electrical propagation in atrial tissue composed of myocytes coupled with myofibroblasts. These findings can inform experimental validation experiments to further elucidate the role of myofibroblast Ca2+ handling in atrial arrhythmogenesis.

The anti-diabetic drug trelagliptin induces vasodilation via activation of Kv channels and SERCA pumps

AIMS

In this study, we investigated the vasodilatory effects of trelagliptin (a dipeptidyl peptidase-4 inhibitor) and its related mechanisms using rabbit aortic rings.

MAIN METHODS

Arterial tone measurement was performed in rabbit thoracic aortic rings.

KEY FINDINGS

Trelagliptin induced vasodilation in a dose-dependent manner. Pretreatment with the ATP-sensitive K⁺ channel inhibitor glibenclamide, large-conductance Ca²⁺-activated K⁺ channel inhibitor paxilline, and inwardly rectifying K⁺ channel inhibitor Ba²⁺ did not affect the vasodilatory effect of trelagliptin. However, pretreatment with the voltage-dependent K⁺ (Kv) channel inhibitors 4-aminopyridine and tetraethylammonium significantly attenuated the vasodilatory effect of trelagliptin, suggesting that the vasodilatory effect of trelagliptin is associated with Kv channel activation. Although pretreatment with Kv1.5 and Kv2.1 subtype inhibitors did not affect the response to trelagliptin, pretreatment with a Kv7.X subtype inhibitor effectively reduced the vasodilatory effect of trelagliptin. Furthermore, sarco/endoplasmic reticulum Ca²⁺-ATPase (SERCA) pump inhibitors also significantly attenuated the vasodilatory effect of trelagliptin. These effects, however, were not affected by pretreatment with Ca²⁺ channel inhibitors, adenylyl cyclase/PKA inhibitors, guanylyl cyclase/PKG inhibitors, or removal of the endothelium.

SIGNIFICANCE

From these results, we concluded that the vasodilatory effect of trelagliptin was associated with the activation of Kv channels (primary the Kv7.X subtype) and SERCA pump regardless of other K⁺ channels, Ca²⁺ channels, cAMP/PKA-related or cGMP/PKG-related signaling pathways, and the endothelium. Therefore, caution is required when prescribing trelagliptin to the patients with hypotension and diabetes.

Using Machine Learning to Characterize Atrial Fibrotic Substrate From Intracardiac Signals With a Hybrid in silico and in vivo Dataset

In patients with atrial fibrillation, intracardiac electrogram signal amplitude is known to decrease with increased structural tissue remodeling, referred to as fibrosis. In addition to the isolation of the pulmonary veins, fibrotic sites are considered a suitable target for catheter ablation. However, it remains an open challenge to find fibrotic areas and to differentiate their density and transmurality. This study aims to identify the volume fraction and transmurality of fibrosis in the atrial substrate. Simulated cardiac electrograms, combined with a generalized model of clinical noise, reproduce clinically measured signals. Our hybrid dataset approach combines in silico and clinical electrograms to train a decision tree classifier to characterize the fibrotic atrial substrate. This approach captures different in vivo dynamics of the electrical propagation reflected on healthy electrogram morphology and synergistically combines it with synthetic fibrotic electrograms from in silico experiments. The machine learning algorithm was tested on five patients and compared against clinical voltage maps as a proof of concept, distinguishing non-fibrotic from fibrotic tissue and characterizing the patient's fibrotic tissue in terms of density and transmurality. The proposed approach can be used to overcome a single voltage cut-off value to identify fibrotic tissue and guide ablation targeting fibrotic areas.

Altered calcium handling in cardiomyocytes from arginine-glycine amidinotransferase-knockout mice is rescued by creatine

The creatine kinase system facilitates energy transfer between mitochondria and the major ATPases in the heart. Creatine-deficient mice, which lack arginine-glycine amidinotransferase (AGAT) to synthesize creatine and homoarginine, exhibit reduced cardiac contractility. We studied how the absence of a functional CK system influences calcium handling in isolated cardiomyocytes from AGAT-knockouts and wild-type littermates as well as in AGAT-knockout mice receiving lifelong creatine supplementation via the food. Using a combination of whole cell patch clamp and fluorescence microscopy, we demonstrate that the L-type calcium channel (LTCC) current amplitude and voltage range of activation were significantly lower in AGAT-knockout compared with wild-type littermates. Additionally, the inactivation of LTCC and the calcium transient decay were significantly slower. According to our modeling results, these changes can be reproduced by reducing three parameters in knockout mice when compared with wild-type: LTCC conductance, the exchange constant of Ca2+ transfer between subspace and cytosol, and SERCA activity. Because tissue expression of LTCC and SERCA protein were not significantly different between genotypes, this suggests the involvement of posttranslational regulatory mechanisms or structural reorganization. The AGAT-knockout phenotype of calcium handling was fully reversed by dietary creatine supplementation throughout life. Our results indicate reduced calcium cycling in cardiomyocytes from AGAT-knockouts and suggest that the creatine kinase system is important for the development of calcium handling in the heart.NEW & NOTEWORTHY Creatine-deficient mice lacking arginine-glycine amidinotransferase exhibit compromised cardiac function. Here, we show that this is at least partially due to an overall slowing of calcium dynamics. Calcium influx into the cytosol via the L-type calcium current (LTCC) is diminished, and the rate of the sarcoendoplasmic reticulum calcium ATPase (SERCA) pumping calcium back into the sarcoplasmic reticulum is slower. The expression of LTCC and SERCA did not change, suggesting that the changes are regulatory.

Afterdepolarizations and abnormal calcium handling in atrial myocytes with modulated SERCA uptake: a sensitivity analysis of calcium handling channels

Delayed afterdepolarizations (DADs) and spontaneous depolarizations (SDs) are typically triggered by spontaneous diastolic Ca²⁺ release from the sarcoplasmic reticulum (SR) which is caused by an elevated SR Ca²⁺-ATPase (SERCA) uptake and dysfunctional ryanodine receptors. However, recent studies on the T-box transcription factor gene (TBX5) demonstrated that abnormal depolarizations could occur despite a reduced SERCA uptake. Similar findings have also been reported in experimental or clinical studies of diabetes and heart failure. To investigate the sensitivity of SERCA in the genesis of DADs/SDs as well as its dependence on other Ca²⁺ handling channels, we performed systematic analyses using the Maleckar et al. model. Results showed that the modulation of SERCA alone cannot trigger abnormal depolarizations, but can instead affect the interdependency of other Ca²⁺ handling channels in triggering DADs/SDs. Furthermore, we discovered the existence of a threshold value for the intracellular concentration of Ca²⁺ ([Ca²⁺]_i) for abnormal depolarizations, which is modulated by the maximum SERCA uptake and the concentration of Ca²⁺ in the uptake and release compartments in the SR ([Ca²⁺]_up and [Ca²⁺]_rel). For the first time, our modelling study reconciles different mechanisms of abnormal depolarizations in the setting of 'lone' AF, reduced TBX5, diabetes and heart failure, and may lead to more targeted treatment for these patients. This article is part of the theme issue 'Uncertainty quantification in cardiac and cardiovascular modelling and simulation'.

Vildagliptin, an Anti-diabetic Drug of the DPP-4 Inhibitor, Induces Vasodilation via Kv Channel and SERCA Pump Activation in Aortic Smooth Muscle

This study investigated vildagliptin-induced vasodilation and its related mechanisms using phenylephrine induced precontracted rabbit aortic rings. Vildagliptin induced vasodilation in a concentration-dependent manner. Pretreatment with the large-conductance Ca²⁺-activated K⁺ channel blocker paxilline, ATP-sensitive K⁺ channel blocker glibenclamide, and inwardly rectifying K⁺ channel blocker Ba²⁺ did not affect the vasodilatory effects of vildagliptin. However, application of the voltage-dependent K⁺ (Kv) channel inhibitor 4-aminopyridine significantly reduced the vasodilatory effects of vildagliptin. In addition, application of either of two sarcoplasmic/endoplasmic reticulum Ca²⁺-ATPase (SERCA) inhibitors, thapsigargin or cyclopiazonic acid, effectively inhibited the vasodilatory effects of vildagliptin. These vasodilatory effects were not affected by pretreatment with adenylyl cyclase, protein kinase A (PKA), guanylyl cyclase, or protein kinase G (PKG) inhibitors, or by removal of the endothelium. From these results, we concluded that vildagliptin induced vasodilation via activation of Kv channels and the SERCA pump. However, other K⁺ channels, PKA/PKG-related signaling cascades associated with vascular dilation, and the endothelium were not involved in vildagliptin-induced vasodilation.

Fast calcium transients in dendritic spines driven by extreme statistics

Fast calcium transients (<10 ms) remain difficult to analyse in cellular microdomains, yet they can modulate key cellular events such as trafficking, local ATP production by endoplasmic reticulum-mitochondria complex (ER-mitochondria complex), or spontaneous activity in astrocytes. In dendritic spines receiving synaptic inputs, we show here that in the presence of a spine apparatus (SA), which is an extension of the smooth ER, a calcium-induced calcium release (CICR) is triggered at the base of the spine by the fastest calcium ions arriving at a Ryanodyne receptor (RyR). The mechanism relies on the asymmetric distributions of RyRs and sarco/ER calcium-ATPase (SERCA) pumps that we predict using a computational model and further confirm experimentally in culture and slice hippocampal neurons. The present mechanism for which the statistics of the fastest particles arriving at a small target, followed by an amplification, is likely to be generic in molecular transduction across cellular microcompartments, such as thin neuronal processes, astrocytes, endfeets, or protrusions.

Thermodynamics of Cation Binding to the Sarcoendoplasmic Reticulum Calcium ATPase Pump and Impacts on Enzyme Function

Sarcoendoplasmic reticulum Ca²⁺-ATPase (SERCA) is a transmembrane pump that plays an important role in transporting calcium into the sarcoplasmic reticulum (SR). While calcium (Ca²⁺) binds SERCA with micromolar affinity, magnesium (Mg²⁺) and potassium (K⁺) also compete with Ca²⁺ binding. However, the molecular bases for these competing ions' influence on the SERCA function and the selectivity of the pump for Ca²⁺ are not well-established. We therefore used in silico methods to resolve molecular determinants of cation binding in the canonical site I and II Ca²⁺ binding sites via (1) triplicate molecular dynamics (MD) simulations of Mg²⁺, Ca²⁺, and K⁺-bound SERCA, (2) mean spherical approximation (MSA) theory to score the affinity and selectivity of cation binding to the MD-resolved structures, and (3) state models of SERCA turnover informed from MSA-derived affinity data. Our key findings are that (a) coordination at sites I and II is optimized for Ca²⁺ and to a lesser extent for Mg²⁺ and K⁺, as determined by MD-derived cation-amino acid oxygen and bound water configurations, (b) the impaired coordination and high desolvation cost for Mg²⁺ precludes favorable Mg²⁺ binding relative to Ca²⁺, while K⁺ has limited capacity to bind site I, and (c) Mg²⁺ most likely acts as inhibitor and K⁺ as intermediate in SERCA's reaction cycle, based on a best-fit state model of SERCA turnover. These findings provide a quantitative basis for SERCA function that leverages molecular-scale thermodynamic data and rationalizes enzyme activity across broad ranges of K⁺, Ca²⁺, and Mg²⁺ concentrations.

Primary cardiac manifestation of autosomal dominant polycystic kidney disease revealed by patient induced pluripotent stem cell-derived cardiomyocytes

Background

Mutations in PKD1 or PKD2 gene lead to autosomal dominant polycystic kidney disease (ADPKD). The mechanism of ADPKD progression and its link to increased cardiovascular mortality is still elusive.

Methods

We differentiated ADPKD patient induced pluripotent stem cells (iPSCs) to cardiomyocytes (CMs). The electrophysiological properties at the cellular level were analyzed by calcium imaging and whole cell patch clamping.

Findings

The ADPKD patient iPSC-CMs had decreased sarcoplasmic reticulum calcium content compared with Control-CMs. Spontaneous action potential of the PKD2 mutation line-derived CMs demonstrated slower beating rate and longer action potential duration. The PKD1 mutation line-derived CMs showed a comparable dose-dependent shortening of phase II repolarization with the Control-CMs, but a significant increase in beating frequency in response to L-type calcium channel blocker. The PKD1-mutant iPSC-CMs also showed a relatively unstable baseline as a greater percentage of cells exhibited delayed afterdepolarizations (DADs). Both the ADPKD patient iPSC-CMs showed more β-adrenergic agonist-elicited DADs compared with Control-CMs.

Interpretation

Characterization of ADPKD patient iPSC-CMs provides new insights into the increased clinical risk of arrhythmias, and the results enable disease modeling and drug screening for cardiac manifestations of ADPKD.

Fund

Ministry of Science and Technology, National Health Research Institutes, Academia Sinica Program for Technology Supporting Platform Axis Scheme, Thematic Research Program and Summit Research Program, and Kaohsiung Medical University Hospital, Taiwan.

Voltage-Dependent Sarcolemmal Ion Channel Abnormalities in the Dystrophin-Deficient Heart

Mutations in the gene encoding for the intracellular protein dystrophin cause severe forms of muscular dystrophy. These so-called dystrophinopathies are characterized by skeletal muscle weakness and degeneration. Dystrophin deficiency also gives rise to considerable complications in the heart, including cardiomyopathy development and arrhythmias. The current understanding of the pathomechanisms in the dystrophic heart is limited, but there is growing evidence that dysfunctional voltage-dependent ion channels in dystrophin-deficient cardiomyocytes play a significant role. Herein, we summarize the current knowledge about abnormalities in voltage-dependent sarcolemmal ion channel properties in the dystrophic heart, and discuss the potentially underlying mechanisms, as well as their pathophysiological relevance.

ATP triggers a robust intracellular [Ca2+ ]-mediated signalling pathway in human synovial fibroblasts

NEW FINDINGS

What is the central question of this study? What are the main [Ca²⁺ ]_i signalling pathways activated by ATP in human synovial fibroblasts? What is the main finding and its importance? In human synovial fibroblasts ATP acts through a linked G-protein (G_q ) and phospholipase C signalling mechanism to produce IP₃ , which then markedly enhances release of Ca²⁺ from the endoplasmic reticulum. These results provide new information for the detection of early pathophysiology of arthritis.

ABSTRACT

In human articular joints, synovial fibroblasts (HSFs) have essential physiological functions that include synthesis and secretion of components of the extracellular matrix and essential articular joint lubricants, as well as release of paracrine substances such as ATP. Although the molecular and cellular processes that lead to a rheumatoid arthritis (RA) phenotype are not fully understood, HSF cells exhibit significant changes during this disease progression. The effects of ATP on HSFs were studied by monitoring changes in intracellular Ca²⁺ ([Ca²⁺ ]_i ), and measuring electrophysiological properties. ATP application to HSF cell populations that had been enzymatically released from 2-D cell culture revealed that ATP (10-100 μm), or its analogues UTP or ADP, consistently produced a large transient increase in [Ca²⁺ ]_i . These changes (i) were initiated by activation of the P₂ Y purinergic receptor family, (ii) required G_q -mediated signal transduction, (iii) did not involve a transmembrane Ca²⁺ influx, but instead (iv) arose almost entirely from activation of endoplasmic reticulum (ER)-localized inositol 1,4,5-trisphosphate (IP₃ ) receptors that triggered Ca²⁺ release from the ER. Corresponding single cell electrophysiological studies revealed that these ATP effects (i) were insensitive to [Ca²⁺ ]_o removal, (ii) involved an IP₃ -mediated intracellular Ca²⁺ release process, and (iii) strongly turned on Ca²⁺ -activated K⁺ current(s) that significantly hyperpolarized these cells. Application of histamine produced very similar effects in these HSF cells. Since ATP is a known paracrine agonist and histamine is released early in the inflammatory response, these findings may contribute to identification of early steps/defects in the initiation and progression of RA.

Importance of Altered Levels of SERCA, IP3R, and RyR in Vascular Smooth Muscle Cell

Calcium cycling between the sarcoplasmic reticulum (SR) and the cytosol via the sarco-/endoplasmic reticulum Ca-ATPase (SERCA) pump, inositol-1,4,5-triphosphate receptor (IP₃R), and Ryanodine receptor (RyR), plays a major role in agonist-induced intracellular calcium ([Ca²⁺]_cyt) dynamics in vascular smooth muscle cells (VSMC). Levels of these calcium handling proteins in SR get altered under disease conditions. We have developed a mathematical model to understand the significance of altered levels of SERCA, IP₃R, and RyR on the intracellular calcium dynamics of VSMC and to understand how variation in protein levels that arise due to diabetes contribute to different VSMC behavior and thus vascular disease. SR is modeled as a single continuous entity with homogeneous intra-SR calcium. Model results show that agonist-induced intracellular calcium dynamics can be modified by changing the levels of SERCA, IP₃R, and/or RyR. Lowering SERCA level will enable intracellular calcium oscillations at low agonist concentrations whereas lowered levels of IP₃R and RyR need higher agonist concentration for intracellular calcium oscillations. This research suggests that reduced SERCA level is the main factor responsible for the reduced intracellular calcium transients and contractility in VSMCs.

Calcium dynamics in cardiac excitatory and non-excitatory cells and the role of gap junction

Calcium ions aid in the generation of action potential in myocytes and are responsible for the excitation-contraction coupling of heart. The heart muscle has specialized patches of cells, called excitatory cells (EC) such as the Sino-atrial node cells capable of auto-generation of action potential and cells which receive signals from the excitatory cells, called non-excitatory cells (NEC) such as cells of the ventricular and auricular walls. In order to understand cardiac calcium homeostasis, it is, therefore, important to study the calcium dynamics taking into account both types of cardiac cells. Here we have developed a model to capture the calcium dynamics in excitatory and non-excitatory cells taking into consideration the gap junction mediated calcium ion transfer from excitatory cell to non-excitatory cell. Our study revealed that the gap junctional coupling between excitatory and non-excitatory cells plays important role in the calcium dynamics. It is observed that any reduction in the functioning of gap junction may result in abnormal calcium oscillations in NEC, even when the calcium dynamics is normal in EC cell. Sensitivity of gap junction is observed to be independent of the pacing rate and hence a careful monitoring is required to maintain normal cardiomyocyte condition. It also highlights that sarcoplasmic reticulum may not be always able to control the amount of cytoplasmic calcium under the condition of calcium overload.

Analysis of Calcium Transients and Uniaxial Contraction Force in Single Human Embryonic Stem Cell-Derived Cardiomyocytes on Microstructured Elastic Substrate with Spatially Controlled Surface Chemistries

The mechanical activity of cardiomyocytes is the result of a process called excitation-contraction coupling (ECC). A membrane depolarization wave induces a transient cytosolic calcium concentration increase that triggers activation of calcium-sensitive contractile proteins, leading to cell contraction and force generation. An experimental setup capable of acquiring simultaneously all ECC features would have an enormous impact on cardiac drug development and disease study. In this work, we develop a microengineered elastomeric substrate with tailor-made surface chemistry to measure simultaneously the uniaxial contraction force and the calcium transients generated by single human cardiomyocytes in vitro. Microreplication followed by photocuring is used to generate an array consisting of elastomeric micropillars. A second photochemical process is employed to spatially control the surface chemistry of the elastomeric pillar. As result, human embryonic stem cell-derived cardiomyocytes (hESC-CMs) can be confined in rectangular cell-adhesive areas, which induce cell elongation and promote suspended cell anchoring between two adjacent micropillars. In this end-to-end conformation, confocal fluorescence microscopy allows simultaneous detection of calcium transients and micropillar deflection induced by a single-cell uniaxial contraction force. Computational finite elements modeling (FEM) and 3D reconstruction of the cell-pillar interface allow force quantification. The platform is used to follow calcium dynamics and contraction force evolution in hESC-CMs cultures over the course of several weeks. Our results show how a biomaterial-based platform can be a versatile tool for in vitro assaying of cardiac functional properties of single-cell human cardiomyocytes, with applications in both in vitro developmental studies and drug screening on cardiac cultures.

In silico screening of the key cellular remodeling targets in chronic atrial fibrillation

Chronic atrial fibrillation (AF) is a complex disease with underlying changes in electrophysiology, calcium signaling and the structure of atrial myocytes. How these individual remodeling targets and their emergent interactions contribute to cell physiology in chronic AF is not well understood. To approach this problem, we performed in silico experiments in a computational model of the human atrial myocyte. The remodeled function of cellular components was based on a broad literature review of in vitro findings in chronic AF, and these were integrated into the model to define a cohort of virtual cells. Simulation results indicate that while the altered function of calcium and potassium ion channels alone causes a pronounced decrease in action potential duration, remodeling of intracellular calcium handling also has a substantial impact on the chronic AF phenotype. We additionally found that the reduction in amplitude of the calcium transient in chronic AF as compared to normal sinus rhythm is primarily due to the remodeling of calcium channel function, calcium handling and cellular geometry. Finally, we found that decreased electrical resistance of the membrane together with remodeled calcium handling synergistically decreased cellular excitability and the subsequent inducibility of repolarization abnormalities in the human atrial myocyte in chronic AF. We conclude that the presented results highlight the complexity of both intrinsic cellular interactions and emergent properties of human atrial myocytes in chronic AF. Therefore, reversing remodeling for a single remodeled component does little to restore the normal sinus rhythm phenotype. These findings may have important implications for developing novel therapeutic approaches for chronic AF.

Translating intracellular calcium signaling into models

The rich experimental data on intracellular calcium has put theoreticians in an ideal position to derive models of intracellular calcium signaling. Over the last 25 years, a large number of modeling frameworks have been suggested. Here, I will review some of the milestones of intracellular calcium modeling with a special emphasis on calcium-induced calcium release (CICR) through inositol-1,4,5-trisphosphate and ryanodine receptors. I will highlight key features of CICR and how they are represented in models as well as the challenges that theoreticians face when translating our current understanding of calcium signals into equations. The selected examples demonstrate that a successful model provides mechanistic insights into the molecular machinery of the Ca²⁺ signaling toolbox and determines the contribution of local Ca²⁺ release to global Ca²⁺ patterns, which at the moment cannot be resolved experimentally.

Enhanced currents through L-type calcium channels in cardiomyocytes disturb the electrophysiology of the dystrophic heart

Duchenne muscular dystrophy (DMD), induced by mutations in the gene encoding for the cytoskeletal protein dystrophin, is an inherited disease characterized by progressive muscle weakness. Besides the relatively well characterized skeletal muscle degenerative processes, DMD is also associated with cardiac complications. These include cardiomyopathy development and cardiac arrhythmias. The current understanding of the pathomechanisms in the heart is very limited, but recent research indicates that dysfunctional ion channels in dystrophic cardiomyocytes play a role. The aim of the present study was to characterize abnormalities in L-type calcium channel function in adult dystrophic ventricular cardiomyocytes. By using the whole cell patch-clamp technique, the properties of currents through calcium channels in ventricular cardiomyocytes isolated from the hearts of normal and dystrophic adult mice were compared. Besides the commonly used dystrophin-deficient mdx mouse model for human DMD, we also used mdx-utr mice, which are both dystrophin- and utrophin-deficient. We found that calcium channel currents were significantly increased, and channel inactivation was reduced in dystrophic cardiomyocytes. Both effects enhance the calcium influx during an action potential (AP). Whereas the AP in dystrophic mouse cardiomyocytes was nearly normal, implementation of the enhanced dystrophic calcium conductance in a computer model of a human ventricular cardiomyocyte considerably prolonged the AP. Finally, the described dystrophic calcium channel abnormalities entailed alterations in the electrocardiograms of dystrophic mice. We conclude that gain of function in cardiac L-type calcium channels may disturb the electrophysiology of the dystrophic heart and thereby cause arrhythmias.

Targeting intracellular calcium cycling in catecholaminergic polymorphic ventricular tachycardia: a theoretical investigation

Catecholaminergic polymorphic ventricular tachycardia (CPVT) is a malignant arrhythmogenic disorder linked to mutations in the cardiac ryanodine receptor (RyR2) and calsequestrin, predisposing the young to syncope and cardiac arrest. To define the role of β-adrenergic stimulation (BAS) and to identify potential therapeutic targeted sites relating to intracellular calcium cycling, we used a Luo-Rudy dynamic ventricular myocyte model incorporated with interacting Markov models of the L-type Ca2+ channel ( ICa,L) and RyR2 to simulate the heterozygous state of mouse RyR2 R4496C mutation (RyR2R4496C+/−) comparable with CPVT patients with RyR2 R4497C mutation. Characteristically, in simulated cells, pacing at 4 Hz or faster or pacing at 2 Hz under BAS with effects equivalent to those of isoproterenol at ≥0.1 μM could readily induce delayed afterdepolarizations (DADs) and DAD-mediated triggered activity (TA) in RyR2R4496C+/− but not in the wild-type via enhancing both ICa,L and sarcoplasmic reticulum (SR) Ca2+ ATPase ( IUP). Moreover, with the use of steady state values of isolated endocardial (Endo), mid-myocardial (M), and epicardial (Epi) cells as initial data for conducting single cell and one-dimensional strand studies, the M cell was more vulnerable for developing DADs and DAD-mediated TA than Endo and Epi cells, and the gap junction coupling represented by diffusion coefficient ( D) of ≤0.000766*98 cm2/ms was required for generating DAD-mediated TA in RyR2R4496C+/−. Whereas individual reduction of Ca2+ release channel of SR and Na-Ca exchanger up to 50% was ineffective, 30% or more reduction of either ICa,L or IUP could totally suppress the inducibility of arrhythmia under BAS. Of note, 15% reduction of both ICa,L and IUP exerted a synergistic antiarrhythmic efficacy. Findings of this model study confirm that BAS facilitates induction of ventricular tachyarrhythmias via its action on intracellular Ca2+ cycling and a pharmacological regimen capable of reducing ICa,L could be an effective adjunctive to β-adrenergic blockers for suppressing ventricular tachyarrhythmias during CPVT.

Impact of sarcoplasmic reticulum calcium release on calcium dynamics and action potential morphology in human atrial myocytes: a computational study

Electrophysiological studies of the human heart face the fundamental challenge that experimental data can be acquired only from patients with underlying heart disease. Regarding human atria, there exist sizable gaps in the understanding of the functional role of cellular Ca²+ dynamics, which differ crucially from that of ventricular cells, in the modulation of excitation-contraction coupling. Accordingly, the objective of this study was to develop a mathematical model of the human atrial myocyte that, in addition to the sarcolemmal (SL) ion currents, accounts for the heterogeneity of intracellular Ca²+ dynamics emerging from a structurally detailed sarcoplasmic reticulum (SR). Based on the simulation results, our model convincingly reproduces the principal characteristics of Ca²+ dynamics: 1) the biphasic increment during the upstroke of the Ca²+ transient resulting from the delay between the peripheral and central SR Ca²+ release, and 2) the relative contribution of SL Ca²+ current and SR Ca²+ release to the Ca²+ transient. In line with experimental findings, the model also replicates the strong impact of intracellular Ca²+ dynamics on the shape of the action potential. The simulation results suggest that the peripheral SR Ca²+ release sites define the interface between Ca²+ and AP, whereas the central release sites are important for the fire-diffuse-fire propagation of Ca²+ diffusion. Furthermore, our analysis predicts that the modulation of the action potential duration due to increasing heart rate is largely mediated by changes in the intracellular Na+ concentration. Finally, the results indicate that the SR Ca²+ release is a strong modulator of AP duration and, consequently, myocyte refractoriness/excitability. We conclude that the developed model is robust and reproduces many fundamental aspects of the tight coupling between SL ion currents and intracellular Ca²+ signaling. Thus, the model provides a useful framework for future studies of excitation-contraction coupling in human atrial myocytes.

Modeling CICR in rat ventricular myocytes: voltage clamp studies

Background

The past thirty-five years have seen an intense search for the molecular mechanisms underlying calcium-induced calcium-release (CICR) in cardiac myocytes, with voltage clamp (VC) studies being the leading tool employed. Several VC protocols including lowering of extracellular calcium to affect Ca²⁺ loading of the sarcoplasmic reticulum (SR), and administration of blockers caffeine and thapsigargin have been utilized to probe the phenomena surrounding SR Ca²⁺ release. Here, we develop a deterministic mathematical model of a rat ventricular myocyte under VC conditions, to better understand mechanisms underlying the response of an isolated cell to calcium perturbation. Motivation for the study was to pinpoint key control variables influencing CICR and examine the role of CICR in the context of a physiological control system regulating cytosolic Ca²⁺ concentration ([Ca²⁺]_myo).

Methods

The cell model consists of an electrical-equivalent model for the cell membrane and a fluid-compartment model describing the flux of ionic species between the extracellular and several intracellular compartments (cell cytosol, SR and the dyadic coupling unit (DCU), in which resides the mechanistic basis of CICR). The DCU is described as a controller-actuator mechanism, internally stabilized by negative feedback control of the unit's two diametrically-opposed Ca²⁺ channels (trigger-channel and release-channel). It releases Ca²⁺ flux into the cyto-plasm and is in turn enclosed within a negative feedback loop involving the SERCA pump, regulating[Ca²⁺]_myo.

Results

Our model reproduces measured VC data published by several laboratories, and generates graded Ca²⁺ release at high Ca²⁺ gain in a homeostatically-controlled environment where [Ca²⁺]_myo is precisely regulated. We elucidate the importance of the DCU elements in this process, particularly the role of the ryanodine receptor in controlling SR Ca²⁺ release, its activation by trigger Ca²⁺, and its refractory characteristics mediated by the luminal SR Ca²⁺ sensor. Proper functioning of the DCU, sodium-calcium exchangers and SERCA pump are important in achieving negative feedback control and hence Ca²⁺ homeostasis.

Conclusions

We examine the role of the above Ca²⁺ regulating mechanisms in handling various types of induced disturbances in Ca²⁺ levels by quantifying cellular Ca²⁺ balance. Our model provides biophysically-based explanations of phenomena associated with CICR generating useful and testable hypotheses.

Calcium channel blockade in embryonic cardiac progenitor cells disrupts normal cardiac cell differentiation

We suggest that characterization of processes involved in differentiation of the pluripotential cardiac precursor cells in their embryonic environment will permit identifying pathways important for induction of diverse stem cells toward the cardiac phenotype. Phenotypic characteristics of cardiac cells are their contractile and electrical properties. The objective of the present study was to define whether calcium (Ca(++)) has a regulatory role in the pluripotential precursor cell population during commitment into cardiomyocytes. We used the chick embryo model because of ease of staging the embryos and visibility of heart development. Using the Ca(++) indicator Fluo-3/acetoxymethyl and confocal microscopy, we demonstrated the existence of higher free Ca(++) levels in the cardiogenic precursor cells than in neighboring cell populations outside of the heart fields. Subsequently, gastrulation stage 4/5 chick embryos were set up in modified New cultures in the medium containing either the L-type Ca channel blocker, diltiazem, or the N-type Ca channel inhibitor, ω-conotoxin. The embryos were incubated for 22-24 h during which time the control embryos developed, beating looping hearts. At the end of incubation, exposure to the L-type channel blockade with diltiazem resulted in an inhibition of cardiomyogenesis in the most posterior, uncommitted, part of the heart fields. N-type channel blockade with ω-conotoxin was less intense. Cells in the most anterior cardiogenic regions that were already committed at time of exposure continued to differentiate. Thus, regulation and maintenance of normal cytosolic Ca levels are necessary for the early steps of cardiomyocyte specification and commitment leading to differentiation.