Jussi T. Koivumäki

Cardiac cell modelling: observations from the heart of the cardiac physiome project.

In this manuscript we review the state of cardiac cell modelling in the context of international initiatives such as the IUPS Physiome and Virtual Physiological Human Projects, which aim to integrate computational models across scales and physics. In particular we focus on the relationship between experimental data and model parameterisation across a range of model types and cellular physiological systems. Finally, in the context of parameter identification and model reuse within the Cardiac Physiome, we suggest some future priority areas for this field.

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 Ca2+ 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 Ca2+ dynamics emerging from a structurally detailed sarcoplasmic reticulum (SR). Based on the simulation results, our model convincingly reproduces the principal characteristics of Ca2+ dynamics: 1) the biphasic increment during the upstroke of the Ca2+ transient resulting from the delay between the peripheral and central SR Ca2+ release, and 2) the relative contribution of SL Ca2+ current and SR Ca2+ release to the Ca2+ transient. In line with experimental findings, the model also replicates the strong impact of intracellular Ca2+ dynamics on the shape of the action potential. The simulation results suggest that the peripheral SR Ca2+ release sites define the interface between Ca2+ and AP, whereas the central release sites are important for the fire-diffuse-fire propagation of Ca2+ 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 Ca2+ 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 Ca2+ signaling. Thus, the model provides a useful framework for future studies of excitation-contraction coupling in human atrial myocytes.

Ca2+–calmodulin‐dependent protein kinase II represses cardiac transcription of the L‐type calcium channel α1C‐subunit gene (Cacna1c) by DREAM translocation

Non‐technical summary  In heart muscle cells, fluctuations of intracellular calcium (Ca2+) concentration ([Ca2+]i) at the frequency defined by the heart rate induce contractions of the cells. Over a longer timescale the same fluctuations define the properties of the cells by regulating expressions of specific genes through the activation of a variety of cellular enzymes. In this study, we have characterized a specific cell signalling pathway, explaining how [Ca2+]i regulates the expression of the L‐type calcium channel (LTCC). We show that [Ca2+]i‐activated calmodulin kinase II (CaMKII) activates downstream regulatory element (DRE) binding transcription factor DREAM, which consequently suppresses the expression of LTCCs. By experiments and mathematical modelling we demonstrate that the LTCC downregulation through the Ca2+–CaMKII–DREAM cascade constitutes a physiological feedback mechanism enabling cardiomyocytes to adjust the calcium intrusion through LTCCs to the amount of intracellular calcium detected by CaMKII.

Variable t-tubule organization and Ca2+ homeostasis across the atria

Although t-tubules have traditionally been thought to be absent in atrial cardiomyocytes, recent studies have suggested that t-tubules exist in the atria of large mammals. However, it is unclear whether regional differences in t-tubule organization exist that define cardiomyocyte function across the atria. We sought to investigate regional t-tubule density in pig and rat atria and the consequences for cardiomyocyte Ca(2+) homeostasis. We observed t-tubules in approximately one-third of rat atrial cardiomyocytes, in both tissue cryosections and isolated cardiomyocytes. In a minority (≈10%) of atrial cardiomyocytes, the t-tubular network was well organized, with a transverse structure resembling that of ventricular cardiomyocytes. In both rat and pig atrial tissue, we observed higher t-tubule density in the epicardium than in the endocardium. Consistent with high variability in the distribution of t-tubules and Ca(2+) channels among cells, L-type Ca(2+) current amplitude was also highly variable and steeply dependent on capacitance and t-tubule density. Accordingly, Ca(2+) transients showed great variability in Ca(2+) release synchrony. Simultaneous imaging of the cell membrane and Ca(2+) transients confirmed t-tubule functionality. Results from mathematical modeling indicated that a transmural gradient in t-tubule organization and Ca(2+) release kinetics supports synchronization of contraction across the atrial wall and may underlie transmural differences in the refractory period. In conclusion, our results indicate that t-tubule density is highly variable across the atria. We propose that higher t-tubule density in cells localized in the epicardium may promote synchronization of contraction across the atrial wall.

Calcium signalling in developing cardiomyocytes: implications for model systems and disease

Adult cardiomyocytes exhibit complex Ca2+ homeostasis, enabling tight control of contraction and relaxation. This intricate regulatory system develops gradually, with progressive maturation of specialized structures and increasing capacity of Ca2+ sources and sinks. In this review, we outline current understanding of these developmental processes, and draw parallels to pathophysiological conditions where cardiomyocytes exhibit a striking regression to an immature state of Ca2+ homeostasis. We further highlight the importance of understanding developmental physiology when employing immature cardiomyocyte models such as cultured neonatal cells and stem cells.

Modelling sarcoplasmic reticulum calcium ATPase and its regulation in cardiac myocytes.

When developing large-scale mathematical models of physiology, some reduction in complexity is necessarily required to maintain computational efficiency. A prime example of such an intricate cell is the cardiac myocyte. For the predictive power of the cardiomyocyte models, it is vital to accurately describe the calcium transport mechanisms, since they essentially link the electrical activation to contractility. The removal of calcium from the cytoplasm takes place mainly by the Na+/Ca2+ exchanger, and the sarcoplasmic reticulum Ca2+ ATPase (SERCA). In the present study, we review the properties of SERCA, its frequency-dependent and β-adrenergic regulation, and the approaches of mathematical modelling that have been used to investigate its function. Furthermore, we present novel theoretical considerations that might prove useful for the elucidation of the role of SERCA in cardiac function, achieving a reduction in model complexity, but at the same time retaining the central aspects of its function. Our results indicate that to faithfully predict the physiological properties of SERCA, we should take into account the calcium-buffering effect and reversible function of the pump. This ‘uncomplicated’ modelling approach could be useful to other similar transport mechanisms as well.

Regulation of excitation-contraction coupling in mouse cardiac myocytes: integrative analysis with mathematical modelling

BackgroundThe cardiomyocyte is a prime example of inherently complex biological system with inter- and cross-connected feedback loops in signalling, forming the basic properties of intracellular homeostasis. Functional properties of cells and tissues have been studied e.g. with powerful tools of genetic engineering, combined with extensive experimentation. While this approach provides accurate information about the physiology at the endpoint, complementary methods, such as mathematical modelling, can provide more detailed information about the processes that have lead to the endpoint phenotype.ResultsIn order to gain novel mechanistic information of the excitation-contraction coupling in normal myocytes and to analyze sophisticated genetically engineered heart models, we have built a mathematical model of a mouse ventricular myocyte. In addition to the fundamental components of membrane excitation, calcium signalling and contraction, our integrated model includes the calcium-calmodulin-dependent enzyme cascade and the regulation it imposes on the proteins involved in excitation-contraction coupling. With the model, we investigate the effects of three genetic modifications that interfere with calcium signalling: 1) ablation of phospholamban, 2) disruption of the regulation of L-type calcium channels by calcium-calmodulin-dependent kinase II (CaMK) and 3) overexpression of CaMK. We show that the key features of the experimental phenotypes involve physiological compensatory and autoregulatory mechanisms that bring the system to a state closer to the original wild-type phenotype in all transgenic models. A drastic phenotype was found when the genetic modification disrupts the regulatory signalling system itself, i.e. the CaMK overexpression model.ConclusionThe novel features of the presented cardiomyocyte model enable accurate description of excitation-contraction coupling. The model is thus an applicable tool for further studies of both normal and defective cellular physiology. We propose that integrative modelling as in the present work is a valuable complement to experiments in understanding the causality within complex biological systems such as cardiac myocytes.

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.

Local Ca2+ releases enable rapid heart rates in developing cardiomyocytes

The ability to generate homogeneous intracellular Ca2+ oscillations at high frequency is the basis of the rhythmic contractions of mammalian cardiac myocytes. While the specific mechanisms and structures enabling homogeneous high‐frequency Ca2+ signals in adult cardiomyocytes are well characterized, it is not known how these kind of Ca2+ signals are produced in developing cardiomyocytes. Here we investigated the mechanisms reducing spatial and temporal heterogeneity of cytosolic Ca2+ signals in mouse embryonic ventricular cardiomyocytes. We show that in developing cardiomyocytes the propagating Ca2+ signals are amplified in cytosol by local Ca2+ releases. Local releases are based on regular 3‐D sarcoplasmic reticulum (SR) structures containing SR Ca2+ uptake ATPases (SERCA) and Ca2+ release channels (ryanodine receptors, RyRs) at regular intervals throughout the cytosol. By evoking [Ca2+]i‐induced Ca2+ sparks, the local release sites promote a 3‐fold increase in the cytosolic Ca2+ propagation speed. We further demonstrate by mathematical modelling that without these local release sites the developing cardiomyocytes lose their ability to generate homogeneous global Ca2+ signals at a sufficiently high frequency. The mechanism described here is robust and indispensable for normal mammalian cardiomyocyte function from the first heartbeats during the early embryonic phase till terminal differentiation after birth. These results suggest that local cytosolic Ca2+ releases are indispensable for normal cardiomyocyte development and function of developing heart.

Investigations of the Navβ1b sodium channel subunit in human ventricle; functional characterization of the H162P Brugada syndrome mutant

Brugada syndrome (BrS) is a rare inherited disease that can give rise to ventricular arrhythmia and ultimately sudden cardiac death. Numerous loss-of-function mutations in the cardiac sodium channel Nav1.5 have been associated with BrS. However, few mutations in the auxiliary Navβ1-4 subunits have been linked to this disease. Here we investigated differences in expression and function between Navβ1 and Navβ1b and whether the H162P/Navβ1b mutation found in a BrS patient is likely to be the underlying cause of disease. The impact of Navβ subunits was investigated by patch-clamp electrophysiology, and the obtained in vitro values were used for subsequent in silico modeling. We found that Navβ1b transcripts were expressed at higher levels than Navβ1 transcripts in the human heart. Navβ1 and Navβ1b coexpressed with Nav1.5 induced a negative shift on steady state of activation and inactivation compared with Nav1.5 alone. Furthermore, Navβ1b was found to increase the current level when coexpressed with Nav1.5, Navβ1b/H162P mutated subunit peak current density was reduced by 48% (-645 ± 151 vs. -334 ± 71 pA/pF), V1/2 steady-state inactivation shifted by -6.7 mV (-70.3 ± 1.5 vs. -77.0 ± 2.8 mV), and time-dependent recovery from inactivation slowed by >50% compared with coexpression with Navβ1b wild type. Computer simulations revealed that these electrophysiological changes resulted in a reduction in both action potential amplitude and maximum upstroke velocity. The experimental data thereby indicate that Navβ1b/H162P results in reduced sodium channel activity functionally affecting the ventricular action potential. This result is an important replication to support the notion that BrS can be linked to the function of Navβ1b and is associated with loss-of-function of the cardiac sodium channel.