A model of the ventricular cardiac action potential. Depolarization, repolarization, and their interaction

Observable Atrial and Ventricular Fibrillation Episode Durations Are Conformant With a Power Law Based on System Size and Spatial Synchronization

BACKGROUND

Atrial fibrillation (AF) and ventricular fibrillation (VF) episodes exhibit varying durations, with some spontaneously ending quickly while others persist. A quantitative framework to explain episode durations remains elusive. We hypothesized that observable self-terminating AF and VF episode lengths, whereby durations are known, would conform with a power law based on the ratio of system size and correlation length ([Formula: see text].

METHODS

Using data from computer simulations (2-dimensional sheet and 3-dimensional left-atrial), human ischemic VF recordings (256-electrode sock, n=12 patients), and human AF recordings (64-electrode basket-catheter, n=9 patients; 16-electrode high definition-grid catheter, n=42 patients), conformance with a power law was assessed using the Akaike information criterion, Bayesian information criterion, coefficient of determination (R2, significance=P<0.05) and maximum likelihood estimation. We analyzed fibrillatory episode durations and [Formula: see text], computed by taking the ratio between system size ([Formula: see text], chamber/simulation size) and correlation length (xi, estimated from pairwise correlation coefficients over electrode/node distance).

RESULTS

In all computer models, the relationship between episode durations and [Formula: see text] was conformant with a power law (Aliev-Panfilov R2: 0.90, P<0.001; Courtemanche R2: 0.91, P<0.001; Luo-Rudy R2: 0.61, P<0.001). Observable clinical AF/VF durations were also conformant with a power law relationship (VF R2: 0.86, P<0.001; AF basket R2: 0.91, P<0.001; AF grid R2: 0.92, P<0.001). [Formula: see text] also differentiated between self-terminating and sustained episodes of AF and VF (P<0.001; all systems), as well as paroxysmal versus persistent AF (P<0.001). In comparison, other electrogram metrics showed no statistically significant differences (dominant frequency, Shannon Entropy, mean voltage, peak-peak voltage; P>0.05).

CONCLUSIONS

Observable fibrillation episode durations are conformant with a power law based on system size and correlation length.

Curvature-mediated source and sink effects on the genesis of premature ventricular complexes in long QT syndrome

Premature ventricular complexes (PVCs) are spontaneous excitations occurring in the ventricles of the heart that are associated with ventricular arrhythmias and sudden cardiac death. Under long QT conditions, PVCs can be mediated by repolarization gradient (RG) and early afterdepolarizations (EADs), yet the effects of heterogeneities or geometry of the RG or EAD regions on PVC genesis remain incompletely understood. In this study, we use computer simulation to systematically investigate the effects of the curvature of the RG border region on PVC genesis under long QT conditions. We show that PVCs can be either promoted or suppressed by negative or positive RG border curvature depending on the source and sink conditions. When the origin of oscillation is in the source region and the source is too strong, a positive RG border curvature can promote PVCs by causing the source area to oscillate. When the origin of oscillation is in the sink region, a negative RG border curvature can promote PVCs by causing the sink area to oscillate. Furthermore, EAD-mediated PVCs are also promoted by negative border curvature. We also investigate the effects of wavefront curvature and show that PVCs are promoted by convex but suppressed by concave wavefronts; however, the effect of wavefront curvature is much smaller than that of RG border curvature. In conclusion, besides the increase of RG and occurrence of EADs caused by QT prolongation, the geometry of the RG border plays important roles in PVC genesis, which can greatly increase the risk of arrhythmias in cardiac diseases.NEW & NOTEWORTHY The effects of the curvature or geometry of the repolarization gradient region and wavefront curvature on the genesis of premature ventricular complexes are systematically investigated using computer modeling and simulation. Premature ventricular complexes can be promoted by either positive or negative curvature of the gradient region depending on the source and sink conditions. The underlying mechanisms of the curvature effects are revealed, which provides mechanistic insights into arrhythmogenesis in cardiac diseases.

Applications of model simulation in pharmacological fields and the problems of theoretical reliability

The use of mathematical models has become increasingly prevalent in pharmacological fields, particularly in drug development processes. These models are instrumental in tasks such as designing clinical trials and assessing factors like efficacy, toxicity, and clinical practice. Various types of models have been developed and documented. Nevertheless, emphasizing the reliability of parameter values is crucial, as they play a pivotal role in shaping the behavior of the system. In some instances, parameter values reported previously are treated as fixed values, which can lead to convergence towards values that deviate substantially from those found in actual biological systems. This is especially true when parameter values are determined through fitting to limited observations. To mitigate this risk, the reuse of parameter values from previous reports should be approached with a critical evaluation of their validity. Currently, there is a proposal for a simultaneous search for plausible values for all parameters using comprehensive search algorithms in both pharmacokinetic and pharmacodynamic or systems pharmacological models. Implementing these methodologies can help address issues related to parameter determination. Furthermore, integrating these approaches with methods developed in the field of machine-learning field has the potential to enhance the reliability of parameter values and the resulting model outputs.

Impact of Combined Modulation of Two Potassium Ion Currents on Spiral Waves and Turbulent States in the Heart

In the realm of cardiac research, the control of spiral waves and turbulent states has been a persistent focus for scholars. Among various avenues of investigation, the modulation of ion currents represents a crucial direction. It has been proved that the methods involving combined control of currents are superior to singular approaches. While previous studies have proposed some combination strategies, further reinforcement and supplementation are required, particularly in the context of controlling arrhythmias through the combined regulation of two potassium ion currents. This study employs the Luo-Rudy phase I cardiac model, modulating the maximum conductance of the time-dependent potassium current and the time-independent potassium current, to investigate the effects of this combined modulation on spiral waves and turbulent states. Numerical simulation results indicate that, compared to modulating a single current, combining reductions in the conductance of two potassium ion currents can rapidly control spiral waves and turbulent states in a short duration. This implies that employing blockers for both potassium ion currents concurrently represents a more efficient control strategy. The control outcomes of this study represent a novel and effective combination for antiarrhythmic interventions, offering potential avenues for new antiarrhythmic drug targets.

Cardiac function is regulated by the sodium-dependent inhibition of the sodium-calcium exchanger NCX1

The Na+-Ca2+ exchanger (NCX1) is the dominant Ca2+ extrusion mechanism in cardiac myocytes. NCX1 activity is inhibited by intracellular Na+ via a process known as Na+-dependent inactivation. A central question is whether this inactivation plays a physiological role in heart function. Using CRISPR/Cas9, we inserted the K229Q mutation in the gene (Slc8a1) encoding for NCX1. This mutation removes the Na+-dependent inactivation while preserving transport properties and other allosteric regulations. NCX1 mRNA levels, protein expression, and protein localization are unchanged in K229Q male mice. However, they exhibit reduced left ventricular ejection fraction and fractional shortening, while displaying a prolonged QT interval. K229Q ventricular myocytes show enhanced NCX1 activity, resulting in action potential prolongation, higher incidence of aberrant action potentials, a faster decline of Ca2+ transients, and depressed cell shortening. The results demonstrate that NCX1 Na+-dependent inactivation plays an essential role in heart function by affecting both cardiac excitability and contractility.

Annihilation dynamics during spiral defect chaos revealed by particle models

Pair-annihilation events are ubiquitous in a variety of spatially extended systems and are often studied using computationally expensive simulations. Here, we develop an approach in which we simulate the pair-annihilation of spiral wave tips in cardiac models using a computationally efficient particle model. Spiral wave tips are represented as particles with dynamics governed by diffusive behavior and short-ranged attraction. The parameters for diffusion and attraction are obtained by comparing particle motion to the trajectories of spiral wave tips in cardiac models during spiral defect chaos. The particle model reproduces the annihilation rates of the cardiac models and can determine the statistics of spiral wave dynamics, including its mean termination time. We show that increasing the attraction coefficient sharply decreases the mean termination time, making it a possible target for pharmaceutical intervention.

T-wave inversion through inhomogeneous voltage diffusion within the FK3V cardiac model

The heart beats are due to the synchronized contraction of cardiomyocytes triggered by a periodic sequence of electrical signals called action potentials, which originate in the sinoatrial node and spread through the heart's electrical system. A large body of work is devoted to modeling the propagation of the action potential and to reproducing reliably its shape and duration. Connection of computational modeling of cells to macroscopic phenomenological curves such as the electrocardiogram has been also intense, due to its clinical importance in analyzing cardiovascular diseases. In this work, we simulate the dynamics of action potential propagation using the three-variable Fenton-Karma model that can account for both normal and damaged cells through a the spatially inhomogeneous voltage diffusion coefficient. We monitor the action potential propagation in the cardiac tissue and calculate the pseudo-electrocardiogram that reproduces the R and T waves. The R-wave amplitude varies according to a double exponential law as a function of the (spatially homogeneous, for an isotropic tissue) diffusion coefficient. The addition of spatial inhomogeneity in the diffusion coefficient by means of a defected region representing damaged cardiac cells may result in T-wave inversion in the calculated pseudo-electrocardiogram. The transition from positive to negative polarity of the T-wave is analyzed as a function of the length and the depth of the defected region.

Directed Graph Mapping exceeds Phase Mapping for the detection of simulated 2D meandering rotors in fibrotic tissue with added noise

Cardiac arrhythmias such as atrial fibrillation (AF) are recognised to be associated with re-entry or rotors. A rotor is a wave of excitation in the cardiac tissue that wraps around its refractory tail, causing faster-than-normal periodic excitation. The detection of rotor centres is of crucial importance in guiding ablation strategies for the treatment of arrhythmia. The most popular technique for detecting rotor centres is Phase Mapping (PM), which detects phase singularities derived from the phase of a signal. This method has been proven to be prone to errors, especially in regimes of fibrotic tissue and temporal noise. Recently, a novel technique called Directed Graph Mapping (DGM) was developed to detect rotational activity such as rotors by creating a network of excitation. This research aims to compare the performance of advanced PM techniques versus DGM for the detection of rotors using 64 simulated 2D meandering rotors in the presence of various levels of fibrotic tissue and temporal noise. Four strategies were employed to compare the performances of PM and DGM. These included a visual analysis, a comparison of F2-scores and distance distributions, and calculating p-values using the mid-p McNemar test. Results indicate that in the case of low meandering, fibrosis and noise, PM and DGM yield excellent results and are comparable. However, in the case of high meandering, fibrosis and noise, PM is undeniably prone to errors, mainly in the form of an excess of false positives, resulting in low precision. In contrast, DGM is more robust against these factors as F2-scores remain high, yielding F2≥0.931 as opposed to the best PM F2≥0.635 across all 64 simulations.

Active Micro-Nano-Collaborative Bioelectronic Device for Advanced Electrophysiological Recording

The development of precise and sensitive electrophysiological recording platforms holds the utmost importance for research in the fields of cardiology and neuroscience. In recent years, active micro/nano-bioelectronic devices have undergone significant advancements, thereby facilitating the study of electrophysiology. The distinctive configuration and exceptional functionality of these active micro-nano-collaborative bioelectronic devices offer the potential for the recording of high-fidelity action potential signals on a large scale. In this paper, we review three-dimensional active nano-transistors and planar active micro-transistors in terms of their applications in electro-excitable cells, focusing on the evaluation of the effects of active micro/nano-bioelectronic devices on electrophysiological signals. Looking forward to the possibilities, challenges, and wide prospects of active micro-nano-devices, we expect to advance their progress to satisfy the demands of theoretical investigations and medical implementations within the domains of cardiology and neuroscience research.

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.

Myocardial ischemia simulation based on a multi-scale heart electrophysiology model

BACKGROUND

Myocardial ischemia, caused by insufficient myocardial blood supply, is a leading cause of human death worldwide. Therefore, it is crucial to prioritize the prevention and treatment of this condition. Mathematical modeling is a powerful technique for studying heart diseases.

OBJECTIVE

The aim of this study was to discuss the quantitative relationship between extracellular potassium concentration and the degree of myocardial ischemia directly related to it.

METHODS

A human cardiac electrophysiological multiscale model was developed to calculate action potentials of all cells simultaneously, enhancing efficiency over traditional reaction-diffusion models.

RESULTS

Contrary to the commonly held view that myocardial ischemia is caused by an increase in extracellular potassium concentration, our simulation results indicate that level 1 ischemia is associated with a decrease in extracellular potassium concentration.

CONCLUSION

This unusual finding provides a new perspective on the mechanisms underlying myocardial ischemia and has the potential to lead to the development of new diagnostic and treatment strategies.

Determinants of electrical propagation and propagation block in Arrhythmogenic Cardiomyopathy

Gap junction and ion channel remodeling occur early in Arrhythmogenic Cardiomyopathy (ACM), but their pathogenic consequences have not been elucidated. Here, we identified the arrhythmogenic substrate, consisting of propagation slowing and conduction block, in ACM models expressing two different desmosomal gene variants. Neonatal rat ventricular myocytes were transduced to express variants in genes encoding desmosomal proteins plakoglobin or plakophilin-2. Studies were performed in engineered cells and anisotropic tissues to quantify changes in conduction velocity, formation of unidirectional propagation, cell-cell electrical coupling, and ion currents. Conduction velocity decreased by 71% and 63% in the two ACM models. SB216763, an inhibitor of glycogen synthase kinase-3 beta, restored conduction velocity to near normal levels. Compared to control, both ACM models showed greater propensity for unidirectional conduction block, which increased further at greater stimulation frequencies. Cell-cell electrical conductance measured in cell pairs was reduced by 86% and 87% in the two ACM models. Computer modeling showed close correspondence between simulated and experimentally determined changes in conduction velocity. The simulation identified that reduced cell-cell electrical coupling was the dominant factor leading to slow conduction, while the combination of reduced cell-cell electrical coupling, reduced sodium current and inward rectifier potassium current explained the development of unidirectional block. Expression of two different ACM variants markedly reduced cell-cell electrical coupling and conduction velocity, and greatly increased the likelihood of developing unidirectional block - both key features of arrhythmogenesis. This study provides the first quantitative analysis of cellular electrophysiological changes leading to the substrate of reentrant arrhythmias in early stage ACM.

Mathematical birth of Early Afterdepolarizations in a cardiomyocyte model

Early Afterdepolarizations (EADs) are abnormal behaviors that can lead to cardiac failure and even cardiac death. In this paper we investigate the occurrence and development of these phenomena in a reduced Luo-Rudy cardiac model. Through a comprehensive dynamical analysis, we map out the distinct patterns observed in the parametric plane, differentiating between normal beats without EADs and pathological beats with EADs. By examining the bifurcation structure of the model, we elucidate the dynamical elements associated with these patterns and their transitions. Using a fast-slow analysis, we explore the emergence and evolution of EADs in the model. Notably, our approach combines the two commonly used fast-slow approaches (1-slow-2-fast and 2-slow-1-fast), and we show how both approaches together provide a more complete understanding of this phenomenon.

Bistable spiral wave dynamics in electrically excitable media

We show that a positive feedback loop between sodium current inactivation and wave-front ramp-up speed causes a saddle-node bifurcation to result in bistable planar and spiral waves in electrically excitable media, in which both slow and fast waves are triggered by different stimulation protocols. Moreover, the two types of spiral wave conduction may interact to give rise to more complex spiral wave dynamics. The transitions between different spiral wave behaviors via saddle-node bifurcation can be a candidate mechanism for transitions widely seen in cardiac arrhythmias and neural diseases.

Efficient Patient-Specific Simulations of Ventricular Tachycardia Based on Computed Tomography-Defined Wall Thickness Heterogeneity

BACKGROUND

Electrophysiological mapping of ventricular tachycardia (VT) is tedious and poorly reproducible. Substrate analysis on imaging cannot explicitly display VT circuits.

OBJECTIVES

This study sought to introduce a computed tomography-based model personalization approach, allowing for the simulation of postinfarction VT in a clinically compatible time frame.

METHODS

In 10 patients (age 65 ± 11 years, 9 male) referred for post-VT ablation, computed tomography-derived wall thickness maps were registered to 25 electroanatomical maps (sinus rhythm, paced, and VT). The relationship between wall thickness and electrophysiological characteristics (activation-recovery interval) was analyzed. Wall thickness was then employed to parameterize a fast and tractable organ-scale wave propagation model. Pacing protocols were simulated from multiple sites to test VT induction in silico. In silico VTs were compared to VT circuits mapped clinically.

RESULTS

Clinically, 6 different VTs could be induced with detailed maps in 9 patients. The proposed model allowed for fast simulation (median: 6 min/pacing site). Simulations of steady pacing (600 milliseconds) from 100 different sites/patient never triggered any arrhythmia. Applying S1-S2 or S1-S2-S3 induction schemes allowed for the induction of in silico VTs in the 9 of 10 patients who were clinically inducible. The patient who was not inducible clinically was also noninducible in silico. A total of 42 different VTs were simulated (4.2 ± 2 per patient). Six in silico VTs matched a VT circuit mapped clinically.

CONCLUSIONS

The proposed framework allows for personalized simulations in a matter of hours. In 6 of 9 patients, simulations show re-entrant patterns matching intracardiac recordings.

Numerical study of the drift of scroll waves by optical feedback in cardiac tissue

Nonlinear waves were found in various types of physical, chemical, and biological excitable media, e.g., in heart muscle. They can form three-dimensional (3D) vortices, called scroll waves, that are of particular significance in the heart, as they underlie lethal cardiac arrhythmias. Thus controlling the behavior of scroll waves is interesting and important. Recently, the optical feedback control procedure for two-dimensional vortices, called spiral waves, was developed. It can induce directed linear drift of spiral waves in optogenetically modified cardiac tissue. However, the extension of this methodology to 3D scroll waves is nontrivial, as optogenetic signals only penetrate close to the surface of cardiac tissue. Here we present a study of this extension in a two-variable reaction-diffusion model and in a detailed model of cardiac tissue. We show that the success of the control procedure is determined by the tension of the scroll wave filament. In tissue with positive filament tension the control procedure works in all cases. However, in the case of negative filament tension for a sufficiently large medium, instabilities occur and make drift and control of scroll waves impossible. Because in normal cardiac tissue the filament tension is assumed to be positive, we conclude that the proposed optical feedback scheme can be a robust method in inducing the linear drift of scroll waves that can control their positions in cardiac tissue.

Effect of strain-dependent conduction slowing on the re-entry formation and maintenance in cardiac muscle: 2D computer simulation

The effect of mechano-electrical feedback on re-entry formation and maintenance was studied using a model of myocardial electromechanics that accounts for two components of myocardial conductivity and delayed strain-dependent changes in membrane capacitance that causes a conduction slowing. Two scenarios were simulated in 2D numerical experiments: (i) propagation of an excitation-contraction wave beyond the edge of a nonconductive nonexcitable obstacle; (ii) circulation of a re-entry wave around a nonconductive nonexcitable obstacle. The simulations demonstrated that the delayed strain-dependent deceleration of the conduction waves promotes the detachment of the excitation-contraction waves from the sharp edge of an elongated obstacle and modulates the re-entry waves rotating around a compact obstacle. The data show that the mechano-electrical feedback, together with an increase in the stimulation frequency and an increase in the excitation threshold, is an arrhythmogenic factor that must be taken into account when analyzing the possibility of the re-entry formation.

Low-energy defibrillation using a base-apex epicardial electrode

BACKGROUND

Current implantable cardioverter defibrillators (ICDs) require electric conduction with high voltage and high energy, which can impair cardiac function and induce another malignant arrhythmia. As a result, there has been a demand for an ICD that can effectively operate with lower energy to mitigate the risks of a strong electric shock.

METHODS

A pair of sheet-shaped electrodes covering the heart were analyzed in three configurations (top-bottom, left-right, and front-back) using a heart simulator. We also varied the distance between the two electrodes (clearance) to identify the electrode shape with the lowest defibrillation threshold (DFT). We also investigated the ICD shock waveform, shock direction, and the effect of the backside insulator of the electrode.

RESULTS

The DFT was high when the clearance was too small and the DFT was high even when the clearance was too large, suggesting that an optimal value clearance. The top-bottom electrodes with optimal clearance showed the lowest DFT when the biphasic shocks set the top electrode to a high potential first and then the bottom electrode was set to a high potential. An interval between a first shock waveform and a second shock waveform should be provided for low-energy defibrillation. Because the insulator prevents unnecessary current flow to the backside, the DFT of the electrodes with insulators is less than those without insulators.

CONCLUSION

Painless defibrillation using sheet-shaped electrodes on the epicardium is predicated on the basis of results using a heart simulator.

Electrophysiological effects of stretch-activated ion channels: a systematic computational characterization

Cardiac electrophysiology and mechanics are strongly interconnected. Their interaction is, among others, mediated by mechano-electric feedback through stretch-activated ion channels (SACs). The electrophysiological changes induced by SACs may contribute to arrhythmogenesis, but the precise SAC-induced electrophysiological changes remain incompletely understood. Here, we provide a systematic characterization of stretch effects through three distinguished SACs on cardiac electrophysiology using computational modelling. We implemented potassium-selective, calcium-selective and non-selective SACs in the Tomek-Rodriguez-O'Hara-Rudy model of human ventricular electrophysiology. The model was calibrated to experimental data from isolated cardiomyocytes undergoing stretch, considering inter-species differences, and disease-related remodelling of SACs. SAC-mediated effects on the action potential (AP) were analysed by varying stretch amplitude, application timing and/or duration. Afterdepolarizations of different amplitudes were observed with transient 10-ms stretch stimuli of 15-18% applied during phase 4, while stretch ≥18% during phase 4 elicited triggered APs. Longer stimuli shifted the threshold of AP trigger during phase 4 to lower amplitudes, while shorter stimuli increased it. Continuous stretch provoked electrophysiological remodelling. Furthermore, stretch shortened duration or changed morphology of a subsequent electrically evoked AP, and, if applied during a vulnerable time window with sufficient amplitude, prevented its occurrence because of stretch-induced modulation of sodium and L-type calcium channel gating. These effects were more pronounced with disease-related SAC remodelling due to increased stretch sensitivity of diseased hearts. We showed that SACs may induce afterdepolarizations and triggered activities, and prevent subsequent AP generation or change its morphology. These effects depend on cardiomyocyte stretch characteristics and disease-related SACs remodelling and may contribute to cardiac arrhythmogenesis. KEY POINTS: The interplay between cardiac electrophysiology and mechanics is mediated by mechano-electric feedback through stretch-activated ion channels (SACs). These channels may be pro-arrhythmic, but their precise effect on electrophysiology remains unclear. Here we present a systematic in silico characterization of stretch effects through three SACs by implementing inter-species differences as well as disease-related remodelling of SACs in a novel computational model of human ventricular cardiomyocyte electrophysiology. Our simulations showed that, at the cellular level, SACs may provoke electrophysiological remodelling, afterdepolarizations, triggered activities, change the morphology or shorten subsequent electrically evoked action potentials. The model further suggests that a vulnerable window exists in which stretch prevents the following electrically triggered beat occurrence. The pro-arrhythmic effects of stretch strongly depend on disease-related SAC remodelling as well as on stretch characteristics, such as amplitude, time of application and duration. Our study helps in understanding the role of stretch in cardiac arrhythmogenesis and revealing the underlying cellular mechanisms.

Hypoxia Promotes Atrial Tachyarrhythmias via Opening of ATP-Sensitive Potassium Channels

BACKGROUND

Hypoxia-ischemia predisposes to atrial arrhythmia. Atrial ATP-sensitive potassium channel (KATP) modulation during hypoxia has not been explored. We investigated the effects of hypoxia on atrial electrophysiology in mice with global deletion of KATP pore-forming subunits.

METHODS

Whole heart KATP RNA expression was probed. Whole-cell KATP current and action potentials were recorded in isolated wild-type (WT), Kir6.1 global knockout (6.1-gKO), and Kir6.2 global knockout (6.2-gKO) murine atrial myocytes. Langendorff-perfused hearts were assessed for atrial effective refractory period (ERP), conduction velocity, wavefront path length (WFPL), and arrhymogenicity under normoxia/hypoxia using a microelectrode array and programmed electrical stimulation. Heart histology was assessed.

RESULTS

Expression patterns were essentially identical for all KATP subunit RNA across human heart, whereas in mouse, Kir6.1 and SUR2 (sulphonylurea receptor subunit) were higher in ventricle than atrium, and Kir6.2 and SUR1 were higher in atrium. Compared with WT, 6.2-gKO atrial myocytes had reduced tolbutamide-sensitive current and action potentials were more depolarized with slower upstroke and reduced peak amplitude. Action potential duration was prolonged in 6.1-gKO atrial myocytes, absent of changes in other ion channel gene expression or atrial myocyte hypertrophy. In Langendorff-perfused hearts, baseline atrial ERP was prolonged and conduction velocity reduced in both KATP knockout mice compared with WT, without histological fibrosis. Compared with baseline, hypoxia led to conduction velocity slowing, stable ERP, and WFPL shortening in WT and 6.1-gKO hearts, whereas WFPL was stable in 6.2-gKO hearts due to ERP prolongation with conduction velocity slowing. Tolbutamide reversed hypoxia-induced WFPL shortening in WT and 6.1-gKO hearts through ERP prolongation. Atrial tachyarrhythmias inducible with programmed electrical stimulation during hypoxia in WT and 6.1-gKO mice correlated with WFPL shortening. Spontaneous arrhythmia was not seen.

CONCLUSIONS

KATP block/absence leads to cellular and tissue level atrial electrophysiological modification. Kir6.2 global knockout prevents hypoxia-induced atrial WFPL shortening and atrial arrhythmogenicity to programmed electrical stimulation. This mechanism could be explored translationally to treat ischemically driven atrial arrhythmia.

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.

Sodium channel subpopulations with distinct biophysical properties and subcellular localization enhance cardiac conduction

Sodium (Na+) current is responsible for the rapid depolarization of cardiac myocytes that triggers the cardiac action potential upstroke. Recent studies have illustrated the presence of multiple pools of Na+ channels with distinct biophysical properties and subcellular localization, including clustering of channels at the intercalated disk and along the lateral membrane. Computational studies predict that Na+ channel clusters at the intercalated disk can regulate cardiac conduction via modulation of the narrow intercellular cleft between electrically coupled myocytes. However, these studies have primarily focused on the redistribution of Na+ channels between intercalated disk and lateral membranes and have not considered the distinct biophysical properties of the Na+ channel subpopulations. In this study, we use computational modeling to simulate computational models of single cardiac cells and one-dimensional cardiac tissues and predict the function of distinct Na+ channel subpopulations. Single-cell simulations predict that a subpopulation of Na+ channels with shifted steady-state activation and inactivation voltage dependency promotes an earlier action potential upstroke. In cardiac tissues that account for distinct subcellular spatial localization, simulations predict that shifted Na+ channels contribute to faster and more robust conduction in response to changes in tissue structure (i.e., cleft width), gap junctional coupling, and rapid pacing rates. Simulations predict that the intercalated disk-localized shifted Na+ channels contribute proportionally more to total Na+ charge than lateral membrane-localized Na+ channels. Importantly, our work supports the hypothesis that Na+ channel redistribution may be a critical mechanism by which cells can respond to perturbations to support fast and robust conduction.

Steady-state approximations for Hodgkin-Huxley cell models: Reduction of order for uterine smooth muscle cell model

Multi-scale mathematical bioelectrical models of organs such as the uterus, stomach or heart present challenges both for accuracy and computational tractability. These multi-scale models are typically founded on models of biological cells derived from the classic Hodkgin-Huxley (HH) formalism. Ion channel behaviour is tracked with dynamical variables representing activation or inactivation of currents that relax to steady-state dependencies on cellular membrane voltage. Timescales for relaxation may be orders of magnitude faster than companion ion channel variables or phenomena of physiological interest for the entire cell (such as bursting sequences of action potentials) or the entire organ (such as electromechanical coordination). Exploiting these time scales with steady-state approximations for relatively fast-acting systems is a well-known but often overlooked approach as evidenced by recent published models. We thus investigate feasibility of an extensive reduction of order for an HH-type cell model with steady-state approximations to the full dynamical activation and inactivation ion channel variables. Our effort utilises a published comprehensive uterine smooth muscle cell model that encompasses 19 ordinary differential equations and 105 formulations overall. The numerous ion channel submodels in the published model exhibit relaxation times ranging from order 10-1 to 105 milliseconds. Substitution of the faster dynamic variables with steady-state formulations demonstrates both an accurate reproduction of the full model and substantial improvements in time-to-solve, for test cases performed. Our demonstration here of an effective and relatively straightforward reduction method underlines the particular importance of considering time scales for model simplification before embarking on large-scale computations or parameter sweeps. As a preliminary complement to more intensive reduction of order methods such as parameter sensitivity and bifurcation analysis, this approach can rapidly and accurately improve computational tractability for challenging multi-scale organ modelling efforts.

Graft-host coupling changes can lead to engraftment arrhythmia: a computational study

After myocardial infarction (MI), a significant portion of heart muscle is replaced with scar tissue, progressively leading to heart failure. Human pluripotent stem cell-derived cardiomyocytes (hPSC-CM) offer a promising option for improving cardiac function after MI. However, hPSC-CM transplantation can lead to engraftment arrhythmia (EA). EA is a transient phenomenon arising shortly after transplantation then spontaneously resolving after a few weeks. The underlying mechanism of EA is unknown. We hypothesize that EA may be explained partially by time-varying, spatially heterogeneous, graft-host electrical coupling. Here, we created computational slice models derived from histological images that reflect different configuration of grafts in the infarcted ventricle. We ran simulations with varying degrees of connection imposed upon the graft-host perimeter to assess how heterogeneous electrical coupling affected EA with non-conductive scar, slow-conducting scar and scar replaced by host myocardium. We also quantified the effect of variation in intrinsic graft conductivity. Susceptibility to EA initially increased and subsequently decreased with increasing graft-host coupling, suggesting the waxing and waning of EA is regulated by progressive increases in graft-host coupling. Different spatial distributions of graft, host and scar yielded markedly different susceptibility curves. Computationally replacing non-conductive scar with host myocardium or slow-conducting scar, and increasing intrinsic graft conductivity both demonstrated potential means to blunt EA vulnerability. These data show how graft location, especially relative to scar, along with its dynamic electrical coupling to host, can influence EA burden; moreover, they offer a rational base for further studies aimed to define the optimal delivery of hPSC-CM injection. KEY POINTS: Human pluripotent stem cell-derived cardiomyocytes (hPSC-CM) hold great cardiac regenerative potential but can also cause engraftment arrhythmias (EA). Spatiotemporal evolution in the pattern of electrical coupling between injected hPSC-CMs and surrounding host myocardium may explain the dynamics of EA observed in large animal models. We conducted simulations in histology-derived 2D slice computational models to assess the effects of heterogeneous graft-host electrical coupling on EA propensity, with or without scar tissue. Our findings suggest spatiotemporally heterogeneous graft-host coupling can create an electrophysiological milieu that favours graft-initiated host excitation, a surrogate metric of EA susceptibility. Removing scar from our models reduced but did not abolish the propensity for this phenomenon. Conversely, reduced intra-graft electrical connectedness increased the incidence of graft-initiated host excitation. The computational framework created for this study can be used to generate new hypotheses, targeted delivery of hPSC-CMs.

Bifurcations to transient and oscillatory excitations in inhomogeneous excitable media: Insights into arrhythmogenesis in long QT syndrome

Ventricular arrhythmias are the leading cause of sudden cardiac death. Understanding the mechanisms of arrhythmia initiation is important for developing effective therapeutics for prevention. Arrhythmias can be induced via premature external stimuli or occur spontaneously via dynamical instabilities. Computer simulations have shown that a large repolarization gradient due to regional prolongation of the action potential duration can result in instabilities leading to premature excitations and arrhythmias, but the bifurcation remains to be elucidated. In this study we carry out numerical simulations and linear stability analyses using a one-dimensional heterogeneous cable consisting of the FitzHugh-Nagumo model. We show that a Hopf bifurcation leads to local oscillations, which, once their amplitudes are large enough, lead to spontaneous propagating excitations. Depending on the degree of heterogeneities, these excitations can range from one to many and to be sustained oscillations, manifesting as premature ventricular contractions (PVCs) and sustained arrhythmias. The dynamics depends on the repolarization gradient and the length of the cable. Complex dynamics is also induced by the repolarization gradient. The mechanistic insights from the simple model may help in the understanding of the genesis of PVCs and arrhythmias in long QT syndrome.

Building a search tool for compositely annotated entities using Transformer-based approach: Case study in Biosimulation Model Search Engine (BMSE)

The Transformer-based approaches to solving natural language processing (NLP) tasks such as BERT and GPT are gaining popularity due to their ability to achieve high performance. These approaches benefit from using enormous data sizes to create pre-trained models and the ability to understand the context of words in a sentence. Their use in the information retrieval domain is thought to increase effectiveness and efficiency. This paper demonstrates a BERT-based method (CASBERT) implementation to build a search tool over data annotated compositely using ontologies. The data was a collection of biosimulation models written using the CellML standard in the Physiome Model Repository (PMR). A biosimulation model structurally consists of basic entities of constants and variables that construct higher-level entities such as components, reactions, and the model. Finding these entities specific to their level is beneficial for various purposes regarding variable reuse, experiment setup, and model audit. Initially, we created embeddings representing compositely-annotated entities for constant and variable search (lowest level entity). Then, these low-level entity embeddings were vertically and efficiently combined to create higher-level entity embeddings to search components, models, images, and simulation setups. Our approach was general, so it can be used to create search tools with other data semantically annotated with ontologies - biosimulation models encoded in the SBML format, for example. Our tool is named Biosimulation Model Search Engine (BMSE).

Summation of activation at the branch-stem transition of Mimosa pudica; a comparison with summation in cardiac tissue

In Mimosa pudica plants, local and global responses to environmental stimuli are associated with different types of electrical activity. Non-damaging stimuli (e.g. cooling) generate action potentials (APs), whereas damaging stimuli (e.g. heating) are associated with variation potentials (VPs). Local cooling of Mimosa branches resulted in APs that propagated up to the branch-stem interface and caused drooping of the branch (local response). This electrical activation did not pass the interface. If the branch was triggered by heat, however, a VP was transferred to the stem and caused activation of the entire plant (global response). VPs caused by heat were always preceded by APs and summation of the two types of activation appeared to be necessary for the activation to pass the branch-stem interface. Mechanical cutting of leaves also resulted in VPs preceded by APs, but in those cases a time delay was present between the two activations, which prevented adequate summation and transmission of activation. Simultaneous cold-induced activation of a branch and the stem below the interface occasionally resulted in summation sufficient to activate the stem beyond the interface. To investigate the effect of activation delay on summation, a similar structure of excitable converging pathways, consisting of a star-shaped pattern of neonatal rat heart cells, was used. In this model, summation of activation was not hindered by a small degree of asynchrony. The observations indicate that summation occurs in excitable branching structures and suggest that summation of activation plays a role in the propagation of nocuous stimuli in Mimosa.

Topological charge-density-vector method of identifying filaments of scroll waves

Scroll waves have been found in a variety of three-dimensional excitable media, including physical, chemical, and biological origins. Scroll waves in cardiac tissue are of particular significance as they underlie ventricular fibrillation that can cause sudden death. The behavior of a scroll wave is characterized by a line of phase singularity at its organizing center, known as a filament. A thorough investigation into the filament dynamics is the key to further exploration of the general theory of scroll waves in excitable media and the mechanisms of ventricular fibrillation. In this paper, we propose a method to identify filaments of scroll waves in excitable media. From the definition of the topological charge of filaments, we obtain the discrete expression of the topological charge-density vector, which is useful in calculating the topological charge vectors at each grid in the space directly. The set of starting points of these topological charge vectors represents a set of phase singularities, thereby forming a line of phase singularity, that is, a filament of a scroll wave.

Pacing Dynamics Determines the Arrhythmogenic Mechanism of the CPVT2-Causing CASQ2G112+5X Mutation in a Guinea Pig Ventricular Myocyte Computational Model

Calsequestrin Type 2 (CASQ2) is a high-capacity, low-affinity, Ca2+-binding protein expressed in the sarcoplasmic reticulum (SR) of the cardiac myocyte. Mutations in CASQ2 have been linked to the arrhythmia catecholaminergic polymorphic ventricular tachycardia (CPVT2) that occurs with acute emotional stress or exercise can result in sudden cardiac death (SCD). CASQ2G112+5X is a 16 bp (339-354) deletion CASQ2 mutation that prevents the protein expression due to premature stop codon. Understanding the subcellular mechanisms of CPVT2 is experimentally challenging because the occurrence of arrhythmia is rare. To obtain an insight into the characteristics of this rare disease, simulation studies using a local control stochastic computational model of the Guinea pig ventricular myocyte investigated how the mutant CASQ2s may be responsible for the development of an arrhythmogenic episode under the condition of β-adrenergic stimulation or in the slowing of heart rate afterward once β-adrenergic stimulation ceases. Adjustment of the computational model parameters based upon recent experiments explore the functional changes caused by the CASQ2 mutation. In the simulation studies under rapid pacing (6 Hz), electromechanically concordant cellular alternans appeared under β-adrenergic stimulation in the CPVT mutant but not in the wild-type nor in the non-β-stimulated mutant. Similarly, the simulations of accelerating pacing from slow to rapid and back to the slow pacing did not display alternans but did generate early afterdepolarizations (EADs) during the period of second slow pacing subsequent acceleration of rapid pacing.

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.

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.

An ECG generative model of myocardial infarction

Background and Objective Computer-aided diagnosis (CAD) of Myocardial Infarction (MI) using machine learning depends on a large amount of clinical Electrocardiogram (ECG) data. Existing infarct ECG databases face the problem of class imbalance. Data augmentation using generative simulation models is a new approach to effectively address this problem. Methods A multiscale ECG generative model was established for ECG data augmentation. In the cellular layer, an ischemic Action Potential (AP) model was established to generate APs in cardiomyocytes with different transmural regions of infraction or different ischemic durations. In the tissue layer, a probability-driven cellular automata excitation propagation model was established to simulate the propagation speed and direction of excitation. An infarct tissue model and a coronary artery model were established to describe the spatiotemporal diversity of MI. A ventricle model, a human torso model, and a computational model of surface ECG based on field source theory were established in the heart-torso layer. Results The model generated pathological 12-lead ECGs of MI with different topography and different extent. When simulating different ventricular wall infarction, the lesions appear in the same leads as the clinical 12-lead ECG. The ST-segment decreases and the T-wave amplitude decreases, similar to the clinical ECG features when simulating subendocardial ischemia. The average fidelity of the 12-lead ECG the model generated is 95.6%, according to the designed DTW-GRA distance algorithm. Conclusions The generative model considers the electrophysiological properties of the natural heart, the pathology of myocardial infarction, and the diversity of clinical ECGs. The model can provide many reliable samples for machine learning of MI.

The physics of heart rhythm disorders

The global burden caused by cardiovascular disease is substantial, with heart disease representing the most common cause of death around the world. There remains a need to develop better mechanistic models of cardiac function in order to combat this health concern. Heart rhythm disorders, or arrhythmias, are one particular type of disease which has been amenable to quantitative investigation. Here we review the application of quantitative methodologies to explore dynamical questions pertaining to arrhythmias. We begin by describing single-cell models of cardiac myocytes, from which two and three dimensional models can be constructed. Special focus is placed on results relating to pattern formation across these spatially-distributed systems, especially the formation of spiral waves of activation. Next, we discuss mechanisms which can lead to the initiation of arrhythmias, focusing on the dynamical state of spatially discordant alternans, and outline proposed mechanisms perpetuating arrhythmias such as fibrillation. We then review experimental and clinical results related to the spatio-temporal mapping of heart rhythm disorders. Finally, we describe treatment options for heart rhythm disorders and demonstrate how statistical physics tools can provide insights into the dynamics of heart rhythm disorders.

Noise-induced switching in an oscillator with pulse delayed feedback: A discrete stochastic modeling approach

We study the dynamics of an oscillatory system with pulse delayed feedback and noise of two types: (i) phase noise acting on the oscillator and (ii) stochastic fluctuations of the feedback delay. Using an event-based approach, we reduce the system dynamics to a stochastic discrete map. For weak noise, we find that the oscillator fluctuates around a deterministic state, and we derive an autoregressive model describing the system dynamics. For stronger noise, the oscillator demonstrates noise-induced switching between various deterministic states; our theory provides a good estimate of the switching statistics in the linear limit. We show that the robustness of the system toward this switching is strikingly different depending on the type of noise. We compare the analytical results for linear coupling to numerical simulations of nonlinear coupling and find that the linear model also provides a qualitative explanation for the differences in robustness to both types of noise. Moreover, phase noise drives the system toward higher frequencies, while stochastic delays do not, and we relate this effect to our theoretical results.

Removal of spiral turbulence by virtual electrodes through the use of a circularly polarized electric field

Heart disease is the leading cause of death and is often accompanied by cardiac fibrillation. Defibrillation using the virtual electrode effects is a promising alternative compared to using the high-voltage electric shock in the clinic. Our earlier works [S. L. Murphy, K. D. Kochanek, J. Xu, and E. Arias, NCHS Data Brief 427 (2021); R. A. Gray, A. M. Pertsov, and J. Jalife, Nature 392, 75-78 (1998); F. X. Witkowski, L. J. Leon, P. A. Penkoske, W. R. Giles, M. L. Spano, W. L. Ditto, and A. T. Winfree, Nature 392, 78-82 (1998); M. Santini, C. Pandozi, G. Altamura, G. Gentilucci, M. Villani, M. C. Scianaro, A. Castro, F. Ammirati, and B. Magris, J. Interv. Card. Electrophysiol. 3, 45-51 (1999).] prove that, compared with other external electric fields, a low voltage circularly polarized electric field is more efficient in turning non-excitable defects in cardiac tissue into virtual electrodes. It, therefore, needs lower voltage to stimulate the excitation waves and causes less harm to reset the spiral turbulence of cardiac excitation for defibrillation. In this paper, we investigate the virtual electrode effect of multiple defects by the circularly polarized electric field for the removal of spiral turbulence. Considering different shapes, sizes, and distributions of multiple defects, we reveal the phase locking of stimulated excitations around multiple virtual electrodes. Furthermore, the circularly polarized electric field parameters are optimized to remove the spiral turbulence.

Dynamical mechanism for generation of arrhythmogenic early afterdepolarizations in cardiac myocytes: Insights from in silico electrophysiological models

We analyze the dynamical mechanisms underlying the formation of arrhythmogenic early afterdepolarizations (EADs) in two mathematical models of cardiac cellular electrophysiology: the Sato et al. biophysically detailed model of a rabbit ventricular myocyte of dimension 27 and a reduced version of the Luo-Rudy mammalian myocyte model of dimension 3. Based on a comparison of the two models, with detailed bifurcation analysis using spike-counting techniques and continuation methods in the simple model and numerical explorations in the complex model, we locate the point where the first EAD originates in an unstable branch of periodic orbits. These results serve as a basis to propose a conjectured scheme involving a hysteresis mechanism with the creation of alternans and EADs in the unstable branch. This theoretical scheme fits well with electrophysiological experimental data on EAD generation and hysteresis phenomena. Our findings open the door to the development of novel methods for pro-arrhythmia risk prediction related to EAD generation without actual induction of EADs.

Spiral wave drift under optical feedback in cardiac tissue

Spiral waves occur in various types of excitable media and their dynamics determine the spatial excitation patterns. An important type of spiral wave dynamics is drift, as it can control the position of a spiral wave or eliminate a spiral wave by forcing it to the boundary. In theoretical and experimental studies of the Belousov-Zhabotinsky reaction, it was shown that the most direct way to induce the controlled drift of spiral waves is by application of an external electric field. Mathematically such drift occurs due to the onset of additional gradient terms in the Laplacian operator describing excitable media. However, this approach does not work for cardiac excitable tissue, where an external electric field does not result in gradient terms. In this paper, we propose a method of how to induce a directed linear drift of spiral waves in cardiac tissue, which can be realized as an optical feedback control in tissue where photosensitive ion channels are expressed. We illustrate our method by using the FitzHugh-Nagumo model for cardiac tissue and the generic model of photosensitive ion channels. We show that our method works for continuous and discrete light sources and can effectively move spiral waves in cardiac tissue, or eliminate them by collisions with the boundary or with another spiral wave. We finally implement our method by using a biophysically motivated photosensitive ion channel model included to the Luo-Rudy model for cardiac cells and show that the proposed feedback control also induces directed linear drift of spiral waves in a wide range of light intensities.

Modeling calcium dynamics in neurons with endoplasmic reticulum: existence, uniqueness and an implicit-explicit finite element scheme

Like many other biological processes, calcium dynamics in neurons containing an endoplasmic reticulum is governed by diffusion-reaction equations on interface-separated domains. Interface conditions are typically described by systems of ordinary differential equations that provide fluxes across the interfaces. Using the calcium model as an example of this class of ODE-flux boundary interface problems, we prove the existence, uniqueness and boundedness of the solution by applying comparison theorem, fundamental solution of the parabolic operator and a strategy used in Picard's existence theorem. Then we propose and analyze an efficient implicit-explicit finite element scheme which is implicit for the parabolic operator and explicit for the nonlinear terms. We show that the stability does not depend on the spatial mesh size. Also the optimal convergence rate in H ¹ norm is obtained. Numerical experiments illustrate the theoretical results.

Control of spiral waves in optogenetically modified cardiac tissue by periodic optical stimulation

Resonant drift of nonlinear spiral waves occurs in various types of excitable media under periodic stimulation. Recently a novel methodology of optogenetics has emerged, which allows to affect excitability of cardiac cells and tissues by optical stimuli. In this paper we study if resonant drift of spiral waves in the heart can be induced by a homogeneous weak periodic optical stimulation of cardiac tissue. We use a two-variable and a detailed model of cardiac tissue and add description of depolarizing and hyperpolarizing optogenetic ionic currents. We show that weak periodic optical stimulation at a frequency equal to the natural rotation frequency of the spiral wave induces resonant drift for both depolarizing and hyperpolarizing optogenetic currents. We quantify these effects and study how the speed of the drift and its direction depend on the initial conditions. We also derive analytical formulas based on the response function theory which correctly predict the drift velocity and its trajectory. We conclude that optogenetic methodology can be used for control of spiral waves in cardiac tissue and discuss its possible applications.

Bifurcations and Proarrhythmic Behaviors in Cardiac Electrical Excitations

The heart is a hierarchical dynamic system consisting of molecules, cells, and tissues, and acts as a pump for blood circulation. The pumping function depends critically on the preceding electrical activity, and disturbances in the pattern of excitation propagation lead to cardiac arrhythmia and pump failure. Excitation phenomena in cardiomyocytes have been modeled as a nonlinear dynamical system. Because of the nonlinearity of excitation phenomena, the system dynamics could be complex, and various analyses have been performed to understand the complex dynamics. Understanding the mechanisms underlying proarrhythmic responses in the heart is crucial for developing new ways to prevent and control cardiac arrhythmias and resulting contractile dysfunction. When the heart changes to a pathological state over time, the action potential (AP) in cardiomyocytes may also change to a different state in shape and duration, often undergoing a qualitative change in behavior. Such a dynamic change is called bifurcation. In this review, we first summarize the contribution of ion channels and transporters to AP formation and our knowledge of ion-transport molecules, then briefly describe bifurcation theory for nonlinear dynamical systems, and finally detail its recent progress, focusing on the research that attempts to understand the developing mechanisms of abnormal excitations in cardiomyocytes from the perspective of bifurcation phenomena.

An Integrated Workflow for Building Digital Twins of Cardiac Electromechanics-A Multi-Fidelity Approach for Personalising Active Mechanics

Personalised computer models of cardiac function, referred to as cardiac digital twins, are envisioned to play an important role in clinical precision therapies of cardiovascular diseases. A major obstacle hampering clinical translation involves the significant computational costs involved in the personalisation of biophysically detailed mechanistic models that require the identification of high-dimensional parameter vectors. An important aspect to identify in electromechanics (EM) models are active mechanics parameters that govern cardiac contraction and relaxation. In this study, we present a novel, fully automated, and efficient approach for personalising biophysically detailed active mechanics models using a two-step multi-fidelity solution. In the first step, active mechanical behaviour in a given 3D EM model is represented by a purely phenomenological, low-fidelity model, which is personalised at the organ scale by calibration to clinical cavity pressure data. Then, in the second step, median traces of nodal cellular active stress, intracellular calcium concentration, and fibre stretch are generated and utilised to personalise the desired high-fidelity model at the cellular scale using a 0D model of cardiac EM. Our novel approach was tested on a cohort of seven human left ventricular (LV) EM models, created from patients treated for aortic coarctation (CoA). Goodness of fit, computational cost, and robustness of the algorithm against uncertainty in the clinical data and variations of initial guesses were evaluated. We demonstrate that our multi-fidelity approach facilitates the personalisation of a biophysically detailed active stress model within only a few (2 to 4) expensive 3D organ-scale simulations-a computational effort compatible with clinical model applications.

Automaticity in ventricular myocyte cell pairs with ephaptic and gap junction coupling

Spontaneous electrical activity, or automaticity, in the heart is required for normal physiological function. However, irregular automaticity, in particular, originating from the ventricles, can trigger life-threatening cardiac arrhythmias. Thus, understanding mechanisms of automaticity and synchronization is critical. Recent work has proposed that excitable cells coupled via a shared narrow extracellular cleft can mediate coupling, i.e., ephaptic coupling, that promotes automaticity in cell pairs. However, the dynamics of these coupled cells incorporating both ephaptic and gap junction coupling has not been explored. Here, we show that automaticity and synchronization robustly emerges via a Hopf bifurcation from either (i) increasing the fraction of inward rectifying potassium channels (carrying the I_K1 current) at the junctional membrane or (ii) by decreasing the cleft volume. Furthermore, we explore how heterogeneity in the fraction of potassium channels between coupled cells can produce automaticity of both cells or neither cell, or more rarely in only one cell (i.e., automaticity without synchronization). Interestingly, gap junction coupling generally has minor effects, with only slight changes in regions of parameter space of automaticity. This work provides insight into potentially new mechanisms that promote spontaneous activity and, thus, triggers for arrhythmias in ventricular tissue.

Stochastic termination of spiral wave dynamics in cardiac tissue

Rotating spiral waves are self-organized features in spatially extended excitable media and may play an important role in cardiac arrhythmias including atrial fibrillation (AF). In homogeneous media, spiral wave dynamics are perpetuated through spiral wave breakup, leading to the continuous birth and death of spiral waves, but have a finite probability of termination. In non-homogeneous media, however, heterogeneities can act as anchoring sources that result in sustained spiral wave activity. It is thus unclear how and if AF may terminate following the removal of putative spiral wave sources in patients. Here, we address this question using computer simulations in which a stable spiral wave is trapped by an heterogeneity and is surrounded by spiral wave breakup. We show that, following ablation of spatial heterogeneity to render that region of the medium unexcitable, termination of spiral wave dynamics is stochastic and Poisson-distributed. Furthermore, we show that the dynamics can be accurately described by a master equation using birth and death rates. To validate these predictions in vivo, we mapped spiral wave activity in patients with AF and targeted the locations of spiral wave sources using radiofrequency ablation. Targeted ablation was indeed able to terminate AF, but only after a variable delay of up to several minutes. Furthermore, and consistent with numerical simulations, termination was not accompanied by gradual temporal or spatial organization. Our results suggest that spiral wave sources and tissue heterogeneities play a critical role in the maintenance of AF and that the removal of sources results in spiral wave dynamics with a finite termination time, which could have important clinical implications.

Computational modeling of aberrant electrical activity following remuscularization with intramyocardially injected pluripotent stem cell-derived cardiomyocytes

Acute engraftment arrhythmias (EAs) remain a serious complication of remuscularization therapy. Preliminary evidence suggests that a focal source underlies these EAs stemming from the automaticity of immature pluripotent stem cell-derived cardiomyocytes (PSC-CMs) in nascent myocardial grafts. How these EAs arise though during early engraftment remains unclear. In a series of in silico experiments, we probed the origin of EAs-exploring aspects of altered impulse formation and altered impulse propagation within nascent PSC-CM grafts and at the host-graft interface. To account for poor gap junctional coupling during early PSC-CM engraftment, the voltage dependence of gap junctions and the possibility of ephaptic coupling were incorporated. Inspired by cardiac development, we also studied the contributions of another feature of immature PSC-CMs, circumferential sodium channel (NaCh) distribution in PSC-CMs. Ectopic propagations emerged from nascent grafts of immature PSC-CMs at a rate of <96 bpm. Source-sink effects dictated this rate and contributed to intermittent capture between host and graft. Moreover, ectopic beats emerged from dynamically changing sites along the host-graft interface. The latter arose in part because circumferential NaCh distribution in PSC-CMs contributed to preferential conduction slowing and block of electrical impulses from host to graft myocardium. We conclude that additional mechanisms, in addition to focal ones, contribute to EAs and recognize that their relative contributions are dynamic across the engraftment process.

Finding type and location of the source of cardiac arrhythmias from the averaged flow velocity field using the determinant-trace method

Life threatening cardiac arrhythmias result from abnormal propagation of nonlinear electrical excitation waves in the heart. Finding the locations of the sources of these waves remains a challenging problem. This is mainly due to the low spatial resolution of electrode recordings of these waves. Also, these recordings are subjected to noise. In this paper, we develop a different approach: the AFV-DT method based on an averaged flow velocity (AFV) technique adopted from the analysis of optical flows and the determinant-trace (DT) method used for vector field analysis of dynamical systems. This method can find the location and determine all important types of sources found in excitable media such as focal activity, spiral waves, and waves rotating around obstacles. We test this method on in silico data of various wave excitation patterns obtained using the Luo-Rudy model for cardiac tissue. We show that the method works well for data with low spatial resolutions (up to 8×8) and is stable against noise. Finally, we apply it to two clinical cases and show that it can correctly identify the arrhythmia type and location. We discuss further steps on the development and improvement of this approach.

Cardiac Alternans Occurs through the Synergy of Voltage- and Calcium-Dependent Mechanisms

Cardiac alternans is characterized by alternating weak and strong beats of the heart. This signaling at the cellular level may appear as alternating long and short action potentials (APs) that occur in synchrony with alternating large and small calcium transients, respectively. Previous studies have suggested that alternans manifests itself through either a voltage dependent mechanism based upon action potential restitution or as a calcium dependent mechanism based on refractoriness of calcium release. We use a novel model of cardiac excitation-contraction (EC) coupling in the rat ventricular myocyte that includes 20,000 calcium release units (CRU) each with 49 ryanodine receptors (RyR2s) and 7 L-type calcium channels that are all stochastically gated. The model suggests that at the cellular level in the case of alternans produced by rapid pacing, the mechanism requires a synergy of voltage- and calcium-dependent mechanisms. The rapid pacing reduces AP duration and magnitude reducing the number of L-type calcium channels activating individual CRUs during each AP and thus increases the population of CRUs that can be recruited stochastically. Elevated myoplasmic and sarcoplasmic reticulum (SR) calcium, [Ca2+]myo and [Ca2+]SR respectively, increases ryanodine receptor open probability (Po) according to our model used in this simulation and this increased the probability of activating additional CRUs. A CRU that opens in one beat is less likely to open the subsequent beat due to refractoriness caused by incomplete refilling of the junctional sarcoplasmic reticulum (jSR). Furthermore, the model includes estimates of changes in Na+ fluxes and [Na+]i and thus provides insight into how changes in electrical activity, [Na+]i and sodium-calcium exchanger activity can modulate alternans. The model thus tracks critical elements that can account for rate-dependent changes in [Na+]i and [Ca2+]myo and how they contribute to the generation of Ca2+ signaling alternans in the heart.

Current progress of computational modeling for guiding clinical atrial fibrillation ablation

Atrial fibrillation (AF) is one of the most common arrhythmias, associated with high morbidity, mortality, and healthcare costs, and it places a significant burden on both individuals and society. Anti-arrhythmic drugs are the most commonly used strategy for treating AF. However, drug therapy faces challenges because of its limited efficacy and potential side effects. Catheter ablation is widely used as an alternative treatment for AF. Nevertheless, because the mechanism of AF is not fully understood, the recurrence rate after ablation remains high. In addition, the outcomes of ablation can vary significantly between medical institutions and patients, especially for persistent AF. Therefore, the issue of which ablation strategy is optimal is still far from settled. Computational modeling has the advantages of repeatable operation, low cost, freedom from risk, and complete control, and is a useful tool for not only predicting the results of different ablation strategies on the same model but also finding optimal personalized ablation targets for clinical reference and even guidance. This review summarizes three-dimensional computational modeling simulations of catheter ablation for AF, from the early-stage attempts such as Maze III or circumferential pulmonary vein isolation to the latest advances based on personalized substrate-guided ablation. Finally, we summarize current developments and challenges and provide our perspectives and suggestions for future directions.

Detecting spiral wave tips using deep learning

The chaotic spatio-temporal electrical activity during life-threatening cardiac arrhythmias like ventricular fibrillation is governed by the dynamics of vortex-like spiral or scroll waves. The organizing centers of these waves are called wave tips (2D) or filaments (3D) and they play a key role in understanding and controlling the complex and chaotic electrical dynamics. Therefore, in many experimental and numerical setups it is required to detect the tips of the observed spiral waves. Most of the currently used methods significantly suffer from the influence of noise and are often adjusted to a specific situation (e.g. a specific numerical cardiac cell model). In this study, we use a specific type of deep neural networks (UNet), for detecting spiral wave tips and show that this approach is robust against the influence of intermediate noise levels. Furthermore, we demonstrate that if the UNet is trained with a pool of numerical cell models, spiral wave tips in unknown cell models can also be detected reliably, suggesting that the UNet can in some sense learn the concept of spiral wave tips in a general way, and thus could also be used in experimental situations in the future (ex-vivo, cell-culture or optogenetic experiments).

Localization of Na+ channel clusters in narrowed perinexi of gap junctions enhances cardiac impulse transmission via ephaptic coupling: a model study

Abstract

It has been proposed that when gap junctional coupling is reduced in cardiac tissue, action potential propagation can be supported via ephaptic coupling, a mechanism mediated by negative electric potentials occurring in narrow intercellular clefts of intercalated discs (IDs). Recent studies showed that sodium (Na⁺) channels form clusters near gap junction plaques in nanodomains called perinexi, where the ID cleft is even narrower. To examine the electrophysiological relevance of Na⁺ channel clusters being located in perinexi, we developed a 3D finite element model of two longitudinally abutting cardiomyocytes, with a central Na⁺ channel cluster on the ID membranes. When this cluster was located in the perinexus of a closely positioned gap junction plaque, varying perinexal width greatly modulated impulse transmission from one cell to the other, with narrow perinexi potentiating ephaptic coupling. This modulation occurred via the interplay of Na⁺ currents, extracellular potentials in the cleft and patterns of current flow within the cleft. In contrast, when the Na⁺ channel cluster was located remotely from the gap junction plaque, this modulation by perinexus width largely disappeared. Interestingly, the Na⁺ current in the ID membrane of the pre‐junctional cell switched from inward to outward during excitation, thus contributing ions to the activating channels on the post‐junctional ID membrane. In conclusion, these results indicate that the localization of Na⁺ channel clusters in the perinexi of gap junction plaques is crucial for ephaptic coupling, which is furthermore greatly modulated by perinexal width. These findings are relevant for a comprehensive understanding of cardiac excitation.

Key points

Ephaptic coupling is a cardiac conduction mechanism involving nanoscale‐level interactions between the sodium (Na⁺) current and the extracellular potential in narrow intercalated disc clefts.

When gap junctional coupling is reduced, ephaptic coupling acts in conjunction with the classical cardiac conduction mechanism based on gap junctional current flow.

In intercalated discs, Na⁺ channels form clusters that are preferentially located in the periphery of gap junction plaques, in nanodomains known as perinexi, but the electrophysiological role of these perinexi has never been examined.

In our new 3D finite element model of two cardiac cells abutting each other with their intercalated discs, a Na⁺ channel cluster located inside a narrowed perinexus facilitated impulse transmission via ephaptic coupling.

Our simulations demonstrate the role of narrowed perinexi as privileged sites for ephaptic coupling in pathological situations when gap junctional coupling is decreased.

The openCARP simulation environment for cardiac electrophysiology

BACKGROUND AND OBJECTIVE

Cardiac electrophysiology is a medical specialty with a long and rich tradition of computational modeling. Nevertheless, no community standard for cardiac electrophysiology simulation software has evolved yet. Here, we present the openCARP simulation environment as one solution that could foster the needs of large parts of this community.

METHODS AND RESULTS

openCARP and the Python-based carputils framework allow developing and sharing simulation pipelines which automate in silico experiments including all modeling and simulation steps to increase reproducibility and productivity. The continuously expanding openCARP user community is supported by tailored infrastructure. Documentation and training material facilitate access to this complementary research tool for new users. After a brief historic review, this paper summarizes requirements for a high-usability electrophysiology simulator and describes how openCARP fulfills them. We introduce the openCARP modeling workflow in a multi-scale example of atrial fibrillation simulations on single cell, tissue, organ and body level and finally outline future development potential.

CONCLUSION

As an open simulator, openCARP can advance the computational cardiac electrophysiology field by making state-of-the-art simulations accessible. In combination with the carputils framework, it offers a tailored software solution for the scientific community and contributes towards increasing use, transparency, standardization and reproducibility of in silico experiments.