Multi-Scale Computational Modeling of Spatial Calcium Handling From Nanodomain to Whole-Heart: Overview and Perspectives

Computational modeling of cardiac electrophysiology and arrhythmogenesis: toward clinical translation

The complexity of cardiac electrophysiology, involving dynamic changes in numerous components across multiple spatial (from ion channel to organ) and temporal (from milliseconds to days) scales, makes an intuitive or empirical analysis of cardiac arrhythmogenesis challenging. Multiscale mechanistic computational models of cardiac electrophysiology provide precise control over individual parameters, and their reproducibility enables a thorough assessment of arrhythmia mechanisms. This review provides a comprehensive analysis of models of cardiac electrophysiology and arrhythmias, from the single cell to the organ level, and how they can be leveraged to better understand rhythm disorders in cardiac disease and to improve heart patient care. Key issues related to model development based on experimental data are discussed, and major families of human cardiomyocyte models and their applications are highlighted. An overview of organ-level computational modeling of cardiac electrophysiology and its clinical applications in personalized arrhythmia risk assessment and patient-specific therapy of atrial and ventricular arrhythmias is provided. The advancements presented here highlight how patient-specific computational models of the heart reconstructed from patient data have achieved success in predicting risk of sudden cardiac death and guiding optimal treatments of heart rhythm disorders. Finally, an outlook toward potential future advances, including the combination of mechanistic modeling and machine learning/artificial intelligence, is provided. As the field of cardiology is embarking on a journey toward precision medicine, personalized modeling of the heart is expected to become a key technology to guide pharmaceutical therapy, deployment of devices, and surgical interventions.

Shedding Light on Cardiac Excitation: In Vitro and In Silico Analysis of Native Ca2+ Channel Activation in Guinea Pig Cardiomyocytes Using Organic Photovoltaic Devices

OBJECTIVE

This study aims to explore the potential of organic electrolytic photocapacitors (OEPCs), an innovative photovoltaic device, in mediating the activation of native voltage-gated Cav1.2 channels (ICa,L) in Guinea pig ventricular cardiomyocytes.

METHODS

Whole-cell patch-clamp recordings were employed to examine light-triggered OEPC mediated ICa,L activation, integrating the channel's kinetic properties into a multicompartment cell model to take intracellular ion concentrations into account. A multidomain model was additionally incorporated to evaluate effects of OEPC-mediated stimulation. The final model combines external stimulation, multicompartmental cell simulation, and a patch-clamp amplifier equivalent circuit to assess the impact on achievable intracellular voltage changes.

RESULTS

Light pulses activated ICa,L, with amplitudes similar to voltage-clamp activation and high sensitivity to the L-type Ca2+ channel blocker, nifedipine. Light-triggered ICa,L inactivation exhibited kinetic parameters comparable to voltage-induced inactivation.

CONCLUSION

OEPC-mediated activation of ICa,L demonstrates their potential for nongenetic optical modulation of cellular physiology potentially paving the way for the development of innovative therapies in cardiovascular health. The integrated model proves the light-mediated activation of ICa,L and advances the understanding of the interplay between the patch-clamp amplifier and external stimulation devices.

SIGNIFICANCE

Treating cardiac conduction disorders by minimal-invasive means without genetic modifications could advance therapeutic approaches increasing patients' quality of life compared with conventional methods employing electronic devices.

Digital twins for cardiac electrophysiology: state of the art and future challenges

Cardiac arrhythmias remain a major cause of death and disability. Current antiarrhythmic therapies are effective to only a limited extent, likely in large part due to their mechanism-independent approach. Precision cardiology aims to deliver targeted therapy for an individual patient to maximize efficacy and minimize adverse effects. In-silico digital twins have emerged as a promising strategy to realize the vision of precision cardiology. While there is no uniform definition of a digital twin, it typically employs digital tools, including simulations of mechanistic computer models, based on patient-specific clinical data to understand arrhythmia mechanisms and/or make clinically relevant predictions. Digital twins have become part of routine clinical practice in the setting of interventional cardiology, where commercially available services use digital twins to non-invasively determine the severity of stenosis (computed tomography-based fractional flow reserve). Although routine clinical application has not been achieved for cardiac arrhythmia management, significant progress towards digital twins for cardiac electrophysiology has been made in recent years. At the same time, significant technical and clinical challenges remain. This article provides a short overview of the history of digital twins for cardiac electrophysiology, including recent applications for the prediction of sudden cardiac death risk and the tailoring of rhythm control in atrial fibrillation. The authors highlight the current challenges for routine clinical application and discuss how overcoming these challenges may allow digital twins to enable a significant precision medicine-based advancement in cardiac arrhythmia management.

Enhanced Ca2+-Driven Arrhythmias in Female Patients with Atrial Fibrillation: Insights from Computational Modeling

BACKGROUND AND AIMS

Substantial sex-based differences have been reported in atrial fibrillation (AF), with female patients experiencing worse symptoms, increased complications from drug side effects or ablation, and elevated risk of AF-related stroke and mortality. Recent studies revealed sex-specific alterations in AF-associated Ca2+ dysregulation, whereby female cardiomyocytes more frequently exhibit potentially proarrhythmic Ca2+-driven instabilities compared to male cardiomyocytes. In this study, we aim to gain a mechanistic understanding of the Ca2+-handling disturbances and Ca2+-driven arrhythmogenic events in males vs females and establish their responses to Ca2+-targeted interventions.

METHODS AND RESULTS

We incorporated known sex differences and AF-associated changes in the expression and phosphorylation of key Ca2+-handling proteins and in ultrastructural properties and dimensions of atrial cardiomyocytes into our recently developed 3D atrial cardiomyocyte model that couples electrophysiology with spatially detailed Ca2+-handling processes. Our simulations of quiescent cardiomyocytes show increased incidence of Ca2+ sparks in female vs male myocytes in AF, in agreement with previous experimental reports. Additionally, our female model exhibited elevated propensity to develop pacing-induced spontaneous Ca2+ releases (SCRs) and augmented beat-to-beat variability in action potential (AP)-elicited Ca2+ transients compared with the male model. Parameter sensitivity analysis uncovered precise arrhythmogenic contributions of each component that was implicated in sex and/or AF alterations. Specifically, increased ryanodine receptor phosphorylation in female AF cardiomyocytes emerged as the major SCR contributor, while reduced L-type Ca2+ current was protective against SCRs for male AF cardiomyocytes. Furthermore, simulations of tentative Ca2+-targeted interventions identified potential strategies to attenuate Ca2+-driven arrhythmogenic events in female atria (e.g., t-tubule restoration, and inhibition of ryanodine receptor and sarcoplasmic/endoplasmic reticulum Ca2+-ATPase), and revealed enhanced efficacy when applied in combination.

CONCLUSIONS

Our sex-specific computational models of human atrial cardiomyocytes uncover increased propensity to Ca2+-driven arrhythmogenic events in female compared to male atrial cardiomyocytes in AF, and point to combined Ca2+-targeted interventions as promising approaches to treat AF in female patients. Our study establishes that AF treatment may benefit from sex-dependent strategies informed by sex-specific mechanisms.

Interactions between calcium-induced arrhythmia triggers and the electrophysiological-anatomical substrate underlying the induction of atrial fibrillation

Atrial fibrillation (AF) is the most common cardiac arrhythmia and is sustained by spontaneous focal excitations and re-entry. Spontaneous electrical firing in the pulmonary vein (PV) sleeves is implicated in AF generation. The aim of this simulation study was to identify the mechanisms determining the localisation of AF triggers in the PVs and their contribution to the genesis of AF. A novel biophysical model of the canine atria was used that integrates stochastic, spontaneous subcellular Ca2+ release events (SCRE) with regional electrophysiological heterogeneity in ionic properties and a detailed three-dimensional model of atrial anatomy, microarchitecture and patchy fibrosis. Simulations highlighted the importance of the smaller inward rectifier potassium current (IK1 ) in PV cells compared to the surrounding atria, which enabled SCRE more readily to result in delayed-afterdepolarisations that induced triggered activity. There was a leftward shift in the dependence of the probability of triggered activity on sarcoplasmic reticulum Ca2+ load. This feature was accentuated in 3D tissue compared to single cells (Δ half-maximal [Ca2+ ]SR  = 58 μM vs. 22 μM). In 3D atria incorporating electrical heterogeneity, excitations preferentially emerged from the PV region. These triggered focal excitations resulted in transient re-entry in the left atrium. Addition of fibrotic patches promoted localised emergence of focal excitations and wavebreaks that had a more substantial impact on generating AF-like patterns than the PVs. Thus, a reduced IK1 , less negative resting membrane potential, and fibrosis-induced changes of the electrotonic load all contribute to the emergence of complex excitation patterns from spontaneous focal triggers. KEY POINTS: Focal excitations in the atria are most commonly associated with the pulmonary veins, but the mechanisms for this localisation are yet to be elucidated. We applied a multi-scale computational modelling approach to elucidate the mechanisms underlying such localisations. Myocytes in the pulmonary vein region of the atria have a less negative resting membrane potential and reduced time-independent potassium current; we demonstrate that both of these factors promote triggered activity in single cells and tissues. The less negative resting membrane potential also contributes to heterogeneous inactivation of the fast sodium current, which can enable re-entrant-like excitation patterns to emerge without traditional conduction block.

Dissecting the roles of calcium cycling and its coupling with voltage in the genesis of early afterdepolarizations in cardiac myocyte models

Early afterdepolarizations (EADs) are abnormal depolarizations during the plateau phase of the action potential, which are known to be associated with lethal arrhythmias in the heart. There are two major hypotheses for EAD genesis based on experimental observations, i.e., the voltage (Vm)-driven and intracellular calcium (Ca)-driven mechanisms. In ventricular myocytes, Ca and Vm are bidirectionally coupled, which can affect each other's dynamics and result in new dynamics, however, the roles of Ca cycling and its coupling with Vm in the genesis of EADs have not been well understood. In this study, we use an action potential model that is capable of independent Vm and Ca oscillations to investigate the roles of Vm and Ca coupling in EAD genesis. Four different mechanisms of EADs are identified, which are either driven by Vm oscillations or Ca oscillations alone, or oscillations caused by their interactions. We also use 5 other ventricular action potential models to assess these EAD mechanisms and show that EADs in these models are mainly Vm-driven. These mechanistic insights from our simulations provide a theoretical base for understanding experimentally observed EADs and EAD-related arrhythmogenesis.

G Protein-Coupled Receptor Signaling: New Insights Define Cellular Nanodomains

G protein-coupled receptors are the largest and pharmacologically most important receptor family and are involved in the regulation of most cell functions. Most of them reside exclusively at the cell surface, from where they signal via heterotrimeric G proteins to control the production of second messengers such as cAMP and IP3 as well as the activity of several ion channels. However, they may also internalize upon agonist stimulation or constitutively reside in various intracellular locations. Recent evidence indicates that their function differs depending on their precise cellular localization. This is because the signals they produce, notably cAMP and Ca2+, are mostly bound to cell proteins that significantly reduce their mobility, allowing the generation of steep concentration gradients. As a result, signals generated by the receptors remain confined to nanometer-sized domains. We propose that such nanometer-sized domains represent the basic signaling units in a cell and a new type of target for drug development.

Diversity of cells and signals in the cardiovascular system

This white paper is the outcome of the seventh UC Davis Cardiovascular Research Symposium on Systems Approach to Understanding Cardiovascular Disease and Arrhythmia. This biannual meeting aims to bring together leading experts in subfields of cardiovascular biomedicine to focus on topics of importance to the field. The theme of the 2022 Symposium was 'Cell Diversity in the Cardiovascular System, cell-autonomous and cell-cell signalling'. Experts in the field contributed their experimental and mathematical modelling perspectives and discussed emerging questions, controversies, and challenges in examining cell and signal diversity, co-ordination and interrelationships involved in cardiovascular function. This paper originates from the topics of formal presentations and informal discussions from the Symposium, which aimed to develop a holistic view of how the multiple cell types in the cardiovascular system integrate to influence cardiovascular function, disease progression and therapeutic strategies. The first section describes the major cell types (e.g. cardiomyocytes, vascular smooth muscle and endothelial cells, fibroblasts, neurons, immune cells, etc.) and the signals involved in cardiovascular function. The second section emphasizes the complexity at the subcellular, cellular and system levels in the context of cardiovascular development, ageing and disease. Finally, the third section surveys the technological innovations that allow the interrogation of this diversity and advancing our understanding of the integrated cardiovascular function and dysfunction.

Modeling Calcium Cycling in the Heart: Progress, Pitfalls, and Challenges

Intracellular calcium (Ca) cycling in the heart plays key roles in excitation-contraction coupling and arrhythmogenesis. In cardiac myocytes, the Ca release channels, i.e., the ryanodine receptors (RyRs), are clustered in the sarcoplasmic reticulum membrane, forming Ca release units (CRUs). The RyRs in a CRU act collectively to give rise to discrete Ca release events, called Ca sparks. A cell contains hundreds to thousands of CRUs, diffusively coupled via Ca to form a CRU network. A rich spectrum of spatiotemporal Ca dynamics is observed in cardiac myocytes, including Ca sparks, spark clusters, mini-waves, persistent whole-cell waves, and oscillations. Models of different temporal and spatial scales have been developed to investigate these dynamics. Due to the complexities of the CRU network and the spatiotemporal Ca dynamics, it is challenging to model the Ca cycling dynamics in the cardiac system, particularly at the tissue sales. In this article, we review the progress of modeling of Ca cycling in cardiac systems from single RyRs to the tissue scale, the pros and cons of the current models and different modeling approaches, and the challenges to be tackled in the future.

Synchronization of spatially discordant voltage and calcium alternans in cardiac tissue

The heart is an excitable medium which is excited by membrane potential depolarization and propagation. Membrane potential depolarization brings in calcium (Ca) through the Ca channels to trigger intracellular Ca release for contraction of the heart. Ca also affects voltage via Ca-dependent ionic currents, and thus, voltage and Ca are bidirectionally coupled. It has been shown that the voltage subsystem or the Ca subsystem can generate its own dynamical instabilities which are affected by their bidirectional couplings, leading to complex dynamics of action potential and Ca cycling. Moreover, the dynamics become spatiotemporal in tissue in which cells are diffusively coupled through voltage. A widely investigated spatiotemporal dynamics is spatially discordant alternans (SDA) in which action potential duration (APD) or Ca amplitude exhibits temporally period-2 and spatially out-of-phase patterns, i.e., APD-SDA and Ca-SDA patterns, respectively. However, the mechanisms of formation, stability, and synchronization of APD-SDA and Ca-SDA patterns remain incompletely understood. In this paper, we use cardiac tissue models described by an amplitude equation, coupled iterated maps, and reaction-diffusion equations with detailed physiology (the ionic model) to perform analytical and computational investigations. We show that, when the Ca subsystem is stable, the Ca-SDA pattern always follows the APD-SDA pattern, and thus, they are always synchronized. When the Ca subsystem is unstable, synchronization of APD-SDA and Ca-SDA patterns depends on the stabilities of both subsystems, their coupling strengths, and the spatial scales of the initial Ca-SDA patterns. Spontaneous (initial condition-independent) synchronization is promoted by enhancing APD instability and reducing Ca instability as well as stronger Ca-to-APD and APD-to-Ca coupling, a pattern formation caused by dynamical instabilities. When Ca is more unstable and APD is less unstable or APD-to-Ca coupling is weak, synchronization of APD-SDA and Ca-SDA patterns is promoted by larger initially synchronized Ca-SDA clusters, i.e., initial condition-dependent synchronization. The synchronized APD-SDA and Ca-SDA patterns can be locked in-phase, antiphase, or quasiperiodic depending on the coupling relationship between APD and Ca. These theoretical and simulation results provide mechanistic insights into the APD-SDA and Ca-SDA dynamics observed in experimental studies.