Human atrial action potential and Ca2+ model: sinus rhythm and chronic atrial fibrillation

Characterization of hERG K+ channel inhibition by the new class III antiarrhythmic drug cavutilide

Cavutilide (niferidil, refralon) is a new class III antiarrhythmic drug which effectively terminates persistent atrial fibrillation (AF; 84.6% of patients, mean AF duration 3 months) and demonstrates low risk of torsade de pointes (1.7%). ERG channels of rapid delayed rectifier current(IKr) are the primary target of cavutilide, but the particular reasons of higher effectiveness and lower proarrhythmic risk in comparison with other class III IKr blockers are unclear. The inhibition of hERG channels expressed in CHO-K1 cells by cavutilide was studied using whole-cell patch-clamp. The present study demonstrates high sensitivity of IhERG expressed in CHO-K1 cells to cavutilide (IC50 = 12.8 nM). Similarly to methanesulfonanilide class III agents, but unlike amiodarone and related drugs, cavutilide does not bind to hERG channels in their resting state. However, in contrast to dofetilide, cavutilide binds not only to opened, but also to inactivated channels. Moreover, at positive constantly set membrane potential (+ 60 mV) inhibition of IhERG by 100 nM cavutilide develops faster than at 0 mV and, especially, - 30 mV (τ of inhibition was 78.8, 103, and 153 ms, respectively). Thereby, cavutilide produces IhERG inhibition only when the cell is depolarized. During the same period of time, cavutilide produces greater block of IhERG when the cell is depolarized with 2 Hz frequency, if compared to 0.2 Hz. We suggest that, during the limited time after injection, cavutilide produces stronger inhibition of IKr in fibrillating atrium than in non-fibrillating ventricle. This leads to beneficial combination of antiarrhythmic effectiveness and low proarrhythmicity of cavutilide.

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.

Influence of genetic mutations to atria vulnerability to atrial fibrillation: An in-silico 3D human atria study

BACKGROUND AND OBJECTIVE

Personalized 3D computer models of atria have been extensively implemented in the last yearsas a tool to facilitate the understanding of the mechanisms underlying different forms of arrhythmia, such as atrial fibrillation (AF). Meanwhile, genetic mutations acting on potassium channel dynamics were demonstrated to induce fibrillatory episodes in asymptomatic patients. This research study aims at assessing the effects and the atrial susceptibility to AF of three gain-of-function mutations - namely, KCNH2 T895M, KCNH2 T436M, and KCNE3-V17M - associated with AF outbreaks, using highly detailed 3D atrial models with realistic wall thickness and heterogenous histological properties.

METHODS

The 3D atrial model was generated by reconstructing segmented anatomical structures from CT scans of an AF patient. Modified versions of the Courtemanche human atrial myocyte model were used to reproduce the electrophysiological activity of the WT and of the three mutant cells. Ectopic foci (EF) were simulated in sixteen locations across the atrial mesh using an S1-S2 protocol with two S2 basic cycle lengths (BCL) and eleven coupling intervals in order to induce arrhythmias.

RESULTS

The three genetic mutations at 3D level reduced the APD90. The KCNE3-V17M mutation provoked the highest shortening (55 % in RA and LA with respect to WT), followed by KCNH2 T895M (14 % in RA and 18 % LA with respect to WT)and KCNH2 T436M (7 % in RA and 9 % LA with respect to WT). The KCNE3-V17M mutation led to arrhythmia in 67 % of the cases simulated and in 94 % of ectopic foci considered, at S2 BCL equal to 100 ms. The KCNH2 T436M and KCNH2 T895M mutations increased the vulnerability to AF in a similar way, leading to arrhythmic episodes in 7 % of the simulated conditions, at S2 BCL set to 160 ms. Overall, 60 % of the arrhythmic events generated arise in the left atrium. Spiral waves, multiple rotors and disordered electrical pattern were elicited in the presence of the KCNE3-V17M mutation, exhibiting an instantaneous mean frequency of 7.6 Hz with a mean standard deviation of 1.12 Hz. The scroll waves induced in the presence of the KCNH2 T436M and KCNH2 T895M mutations showed steadiness and regularity with an instantaneous mean frequencies in the range of 4.9 - 5.1 Hz and a mean standard deviation within 0.19 - 0.53 Hz.

CONCLUSIONS

The pro-arrhythmogenicity of the KCNE3-V17M, KCNH2 T895M and KCNH2 T436M mutations was studied and proved on personalized 3D cardiac models. The three genetic mutations were demonstrated to increase the predisposition of atrial tissue to the formation of AF-susceptible substrate in different ways based on their effects on electrophysiological properties of the atria.

Connecting Transcriptomics with Computational Modeling to Reveal Developmental Adaptations in the Human Pediatric Myocardium

Background

Nearly 1% or 1.3 million babies are born with congenital heart disease (CHD) globally each year - many of whom will require palliative or corrective heart surgery within the first few years of life. A detailed understanding of cardiac maturation can help to expand our knowledge on cardiac diseases that develop during gestation, identify age-appropriate cardiovascular drug therapies, and inform clinical care decisions related to surgical repair, myocardial preservation, or postoperative management. Yet, to date, our knowledge of the temporal changes that cardiomyocytes undergo during postnatal development is largely limited to animal models.

Methods

Right atrial tissue samples were collected from n=117 neonatal, infant, and pediatric patients undergoing correct surgery due to (acyanotic) CHD. Patients were stratified into five age groups: neonate (0-30 days), infant (31-364 days), toddler to preschool (1-5 years), school age (6-11 years), and adolescent to young adults (12-32 years). We measured age-dependent adaptations in cardiac gene expression, and used computational modeling to simulate action potential and calcium transients.

Results

Enrichment of differentially expressed genes (DEG) was explored, revealing age-dependent changes in several key biological processes (cell cycle, cell division, mitosis), cardiac ion channels, and calcium handling genes. Gene-associated changes in ionic currents exhibited both linear trends and sudden shifts across developmental stages, with changes in calcium handling ( I NCX ) and repolarization ( I K1 ) most strongly associated with an age-dependent decrease in the action potential plateau potential and increase in triangulation, respectively. We also note a shift in repolarization reserve, with lower I Kr expression in younger patients, a finding likely tied to the increased amplitude of I Ks triggered by elevated sympathetic activation in pediatric patients.

Conclusion

This study provides valuable insights into age-dependent changes in human cardiac gene expression and electrophysiology among patients with CHD, shedding light on molecular mechanisms underlying cardiac development and function across different developmental stages.

The spiral wave frequency effect in atrial fibrillation

A spiral wavefront (WF), generated by a cardiac rotor that drifts between surface electrodes during atrial fibrillation, exhibits frequency changes inconsistent with classical Doppler effect (CDE) phenomena. Recent clinical studies reveal three repeatedly observed events--1) side-dependent frequency changes across the path of the rotor, 2) one additional WF strike on the higher frequency side, and 3) a reversal of WF strike sequence--which constitute a diametrical property of spinning WF sources. A linear ray model is first used to reveal and develop the diametrical phenomena. Mathematical models of an Archimedean spiral and a spiral generated by the diffusion equation are developed and compared. Each formulation predicts the diametrical property that CDE does not capture and illuminates the occurrence of a strong side and weak side with respect to the rotor path. Whereas CDE exhibits higher and lower frequencies from approaching and receding sources of WFs, respectively, spiral rotors generate higher and lower frequencies on opposite sides of the migration path. This motivates the reconsideration of mapping and ablation strategies that have traditionally been based on identifying sites of the dominant frequency. While this research aims to characterize the path of a spiral rotor during atrial fibrillation accurately, the results are applicable in other fields of science and engineering in which rotating spiral waves occur.

Advances in cardiac tissue engineering and heart-on-a-chip

Recent advances in both cardiac tissue engineering and hearts-on-a-chip are grounded in new biomaterial development as well as the employment of innovative fabrication techniques that enable precise control of the mechanical, electrical, and structural properties of the cardiac tissues being modelled. The elongated structure of cardiomyocytes requires tuning of substrate properties and application of biophysical stimuli to drive its mature phenotype. Landmark advances have already been achieved with induced pluripotent stem cell-derived cardiac patches that advanced to human testing. Heart-on-a-chip platforms are now commonly used by a number of pharmaceutical and biotechnology companies. Here, we provide an overview of cardiac physiology in order to better define the requirements for functional tissue recapitulation. We then discuss the biomaterials most commonly used in both cardiac tissue engineering and heart-on-a-chip, followed by the discussion of recent representative studies in both fields. We outline significant challenges common to both fields, specifically: scalable tissue fabrication and platform standardization, improving cellular fidelity through effective tissue vascularization, achieving adult tissue maturation, and ultimately developing cryopreservation protocols so that the tissues are available off the shelf.

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.

Caveolin-3 and Caveolin-1 Interaction Decreases Channel Dysfunction Due to Caveolin-3 Mutations

Caveolae constitute membrane microdomains where receptors and ion channels functionally interact. Caveolin-3 (cav-3) is the key structural component of muscular caveolae. Mutations in CAV3 lead to caveolinopathies, which result in both muscular dystrophies and cardiac diseases. In cardiomyocytes, cav-1 participates with cav-3 to form caveolae; skeletal myotubes and adult skeletal fibers do not express cav-1. In the heart, the absence of cardiac alterations in the majority of cases may depend on a conserved organization of caveolae thanks to the expression of cav-1. We decided to focus on three specific cav-3 mutations (Δ62-64YTT; T78K and W101C) found in heterozygosis in patients suffering from skeletal muscle disorders. We overexpressed both the WT and mutated cav-3 together with ion channels interacting with and modulated by cav-3. Patch-clamp analysis conducted in caveolin-free cells (MEF-KO), revealed that the T78K mutant is dominant negative, causing its intracellular retention together with cav-3 WT, and inducing a significant reduction in current densities of all three ion channels tested. The other cav-3 mutations did not cause significant alterations. Mathematical modelling of the effects of cav-3 T78K would impair repolarization to levels incompatible with life. For this reason, we decided to compare the effects of this mutation in other cell lines that endogenously express cav-1 (MEF-STO and CHO cells) and to modulate cav-1 expression with an shRNA approach. In these systems, the membrane localization of cav-3 T78K was rescued in the presence of cav-1, and the current densities of hHCN4, hKv1.5 and hKir2.1 were also rescued. These results constitute the first evidence of a compensatory role of cav-1 in the heart, justifying the reduced susceptibility of this organ to caveolinopathies.

Computational Modeling and Simulation of the Fibrotic Human Atria

Patient-specific modeling of atrial electrical activity enables the execution of simulations that can provide mechanistic insights and provide novel solutions to vexing clinical problems. The geometry and fibrotic remodeling of the heart can be reconstructed from clinical-grade medical scans and used to inform personalized models with detail incorporated at the cell- and tissue-scale to represent changes in image-identified diseased regions. Here, we provide a rubric for the reconstruction of realistic atrial models from pre-segmented 3D renderings of the left atrium with fibrotic tissue regions delineated, which are the output from clinical-grade systems for quantifying fibrosis. We then provide a roadmap for using those models to carry out patient-specific characterization of the fibrotic substrate in terms of its potential to harbor reentrant drivers via cardiac electrophysiology simulations.

Computational Modeling of Cardiac Electrophysiology

Mathematical modeling and simulation are well-established and powerful tools to integrate experimental data of individual components of cardiac electrophysiology, excitation-contraction coupling, and regulatory signaling pathways, to gain quantitative and mechanistic insight into pathophysiological processes and guide therapeutic strategies. Here, we briefly describe the processes governing cardiac myocyte electrophysiology and Ca2+ handling and their regulation, as well as action potential propagation in tissue. We discuss the models and methods used to describe these phenomena, including procedures for model parameterization and validation, in addition to protocols for model interrogation and analysis and techniques that account for phenotypic variability and parameter uncertainty. Our objective is to provide a summary of basic concepts and approaches as a resource for scientists training in this discipline and for all researchers aiming to gain an understanding of cardiac modeling studies.

Mechanisms of pulmonary arterial hypertension-induced atrial fibrillation: insights from multi-scale models of the human atria

This study aimed to use multi-scale atrial models to investigate pulmonary arterial hypertension (PAH)-induced atrial fibrillation mechanisms. The results of our computer simulations revealed that, at the single-cell level, PAH-induced remodelling led to a prolonged action potential (AP) (ΔAPD: 49.6 ms in the right atria (RA) versus 41.6 ms in the left atria (LA)) and an increased calcium transient (CaT) (ΔCaT: 7.5 × 10-2 µM in the RA versus 0.9 × 10-3 µM in the LA). Moreover, heterogeneous remodelling increased susceptibility to afterdepolarizations, particularly in the RA. At the tissue level, we observed a significant reduction in conduction velocity (CV) (ΔCV: -0.5 m s-1 in the RA versus -0.05 m s-1 in the LA), leading to a shortened wavelength in the RA, but not in the LA. Additionally, afterdepolarizations in the RA contributed to enhanced repolarization dispersion and facilitated unidirectional conduction block. Furthermore, the increased fibrosis in the RA amplified the likelihood of excitation wave breakdown and the occurrence of sustained re-entries. Our results indicated that the RA is characterized by increased susceptibility to afterdepolarizations, slow conduction, reduced wavelength and upregulated fibrosis. These findings shed light on the underlying factors that may promote atrial fibrillation in patients with PAH.

Mechanistic insight into the functional role of human sinoatrial node conduction pathways and pacemaker compartments heterogeneity: A computer model analysis

The sinoatrial node (SAN), the primary pacemaker of the heart, is responsible for the initiation and robust regulation of sinus rhythm. 3D mapping studies of the ex-vivo human heart suggested that the robust regulation of sinus rhythm relies on specialized fibrotically-insulated pacemaker compartments (head, center and tail) with heterogeneous expressions of key ion channels and receptors. They also revealed up to five sinoatrial conduction pathways (SACPs), which electrically connect the SAN with neighboring right atrium (RA). To elucidate the role of these structural-molecular factors in the functional robustness of human SAN, we developed comprehensive biophysical computer models of the SAN based on 3D structural, functional and molecular mapping of ex-vivo human hearts. Our key finding is that the electrical insulation of the SAN except SACPs, the heterogeneous expression of If, INa currents and adenosine A1 receptors (A1R) across SAN pacemaker-conduction compartments are required to experimentally reproduce observed SAN activation patterns and important phenomena such as shifts of the leading pacemaker and preferential SACP. In particular, we found that the insulating border between the SAN and RA, is required for robust SAN function and protection from SAN arrest during adenosine challenge. The heterogeneity in the expression of A1R within the human SAN compartments underlies the direction of pacemaker shift and preferential SACPs in the presence of adenosine. Alterations of INa current and fibrotic remodelling in SACPs can significantly modulate SAN conduction and shift the preferential SACP/exit from SAN. Finally, we show that disease-induced fibrotic remodeling, INa suppression or increased adenosine make the human SAN vulnerable to pacing-induced exit blocks and reentrant arrhythmia. In summary, our computer model recapitulates the structural and functional features of the human SAN and can be a valuable tool for investigating mechanisms of SAN automaticity and conduction as well as SAN arrhythmia mechanisms under different pathophysiological conditions.

Transition mechanisms from atrial flutter to atrial fibrillation during anti-tachycardia pacing therapy

BACKGROUND

Atrial anti-tachycardia pacing (aATP) has been shown to be effective for the termination of atrial tachyarrhythmias, but its success rate is still not high enough.

OBJECTIVE

The main objective of this study was to investigate the mechanisms of atrial flutter (AFL) termination by aATP and the transition from AFL to atrial fibrillation (AF) during aATP.

METHODS

We developed a multi-scale model of the human atrium based on magnetic resonance images and examined the atrial electrophysiology of AFL during aATP with a ramp protocol.

RESULTS

In successful cases of aATP, paced excitation entered the excitable gap and collided with the leading edge of the reentrant wave front. Furthermore, the excitation propagating in the opposite direction collided with the trailing edge of the reentrant wave to terminate AFL. The second collision was made possible by the distribution of the wave propagation velocity in the atria. The detailed analysis revealed that the slowing of propagation velocity occurred at the exit of the sub-Eustachian isthmus, probably due to source-sink mismatch. During the transition from AFL to AF, the excitation collided with the refractory zone of the preceding wave and broke into multiple wave fronts to induce AF. A similar observation was made for the transition from AF to sinus rhythm. In both cases, the complex anatomy of the atria played an essential role.

CONCLUSION

The complex anatomy of atria plays an essential role in the maintenance of stable AFL and its termination by aATP, which were revealed by the realistic three-dimensional simulation model.

Integrative human atrial modelling unravels interactive protein kinase A and Ca2+/calmodulin-dependent protein kinase II signalling as key determinants of atrial arrhythmogenesis

AIMS

Atrial fibrillation (AF), the most prevalent clinical arrhythmia, is associated with atrial remodelling manifesting as acute and chronic alterations in expression, function, and regulation of atrial electrophysiological and Ca2+-handling processes. These AF-induced modifications crosstalk and propagate across spatial scales creating a complex pathophysiological network, which renders AF resistant to existing pharmacotherapies that predominantly target transmembrane ion channels. Developing innovative therapeutic strategies requires a systems approach to disentangle quantitatively the pro-arrhythmic contributions of individual AF-induced alterations.

METHODS AND RESULTS

Here, we built a novel computational framework for simulating electrophysiology and Ca2+-handling in human atrial cardiomyocytes and tissues, and their regulation by key upstream signalling pathways [i.e. protein kinase A (PKA), and Ca2+/calmodulin-dependent protein kinase II (CaMKII)] involved in AF-pathogenesis. Populations of atrial cardiomyocyte models were constructed to determine the influence of subcellular ionic processes, signalling components, and regulatory networks on atrial arrhythmogenesis. Our results reveal a novel synergistic crosstalk between PKA and CaMKII that promotes atrial cardiomyocyte electrical instability and arrhythmogenic triggered activity. Simulations of heterogeneous tissue demonstrate that this cellular triggered activity is further amplified by CaMKII- and PKA-dependent alterations of tissue properties, further exacerbating atrial arrhythmogenesis.

CONCLUSIONS

Our analysis reveals potential mechanisms by which the stress-associated adaptive changes turn into maladaptive pro-arrhythmic triggers at the cellular and tissue levels and identifies potential anti-AF targets. Collectively, our integrative approach is powerful and instrumental to assemble and reconcile existing knowledge into a systems network for identifying novel anti-AF targets and innovative approaches moving beyond the traditional ion channel-based strategy.

Effects of hydroxychloroquine on atrial electrophysiology in in silico wild-type and PITX2+/- atrial cardiomyocytes

BACKGROUND

Hydroxychloroquine (HCQ) is commonly used in the treatment of autoimmune diseases and increases the risk of QT interval prolongation. However, it is unclear how HCQ affects atrial electrophysiology and the risk of atrial fibrillation (AF).

METHODS

We quantitatively examined the potential atrial arrhythmogenic effects of HCQ on AF using a computational model of human atrial cardiomyocytes. We measured atrial electrophysiological markers after systematically varying HCQ concentrations.

RESULTS

The HCQ concentrations were positively correlated with the action potential duration (APD), resting membrane potential, refractory period, APD alternans threshold, and calcium transient alternans threshold (p < 0.05). By contrast, HCQ concentrations were inversely correlated with the maximum upstroke velocity and calcium transient amplitude (p < 0.05). When the therapeutic concentration (Cmax) of HCQ was applied, HCQ increased APD90 by 1.4% in normal sinus rhythm, 1.8% in wild-type AF, and 2.6% in paired-like homeodomain transcription factor 2 (PITX2)+/- AF, but did not affect the alternans thresholds. The overall in silico results suggest no significant atrial arrhythmogenic effects of HCQ at Cmax, instead implying a potential antiarrhythmic role of low-dose HCQ in AF. However, at an HCQ concentration of fourfold Cmax, a rapid pacing rate of 4 Hz induced prominent APD alternans, particularly in the PITX2+/- AF model.

CONCLUSION

Our in silico analysis suggests a potential antiarrhythmic role of low-dose HCQ in AF. Concomitant PITX2 mutations and high-dose HCQ treatments may increase the risk of AF, and this potential genotype/dose-dependent arrhythmogenic effect of HCQ should be investigated further.

Dual effects of the small-conductance Ca2+-activated K+ current on human atrial electrophysiology and Ca2+-driven arrhythmogenesis: an in silico study

By sensing changes in intracellular Ca2+, small-conductance Ca2+-activated K+ (SK) channels dynamically regulate the dynamics of the cardiac action potential (AP) on a beat-to-beat basis. Given their predominance in atria versus ventricles, SK channels are considered a promising atrial-selective pharmacological target against atrial fibrillation (AF), the most common cardiac arrhythmia. However, the precise contribution of SK current (ISK) to atrial arrhythmogenesis is poorly understood, and may potentially involve different mechanisms that depend on species, heart rates, and degree of AF-induced atrial remodeling. Both reduced and enhanced ISK have been linked to AF. Similarly, both SK channel up- and downregulation have been reported in chronic AF (cAF) versus normal sinus rhythm (nSR) patient samples. Here, we use our multiscale modeling framework to obtain mechanistic insights into the contribution of ISK in human atrial cardiomyocyte electrophysiology. We simulate several protocols to quantify how ISK modulation affects the regulation of AP duration (APD), Ca2+ transient, refractoriness, and occurrence of alternans and delayed afterdepolarizations (DADs). Our simulations show that ISK activation shortens the APD and atrial effective refractory period, limits Ca2+ cycling, and slightly increases the propensity for alternans in both nSR and cAF conditions. We also show that increasing ISK counteracts DAD development by enhancing the repolarization force that opposes the Ca2+-dependent depolarization. Taken together, our results suggest that increasing ISK in human atrial cardiomyocytes could promote reentry while protecting against triggered activity. Depending on the leading arrhythmogenic mechanism, ISK inhibition may thus be a beneficial or detrimental anti-AF strategy.NEW & NOTEWORTHY Using our established framework for human atrial myocyte simulations, we investigated the role of the small-conductance Ca2+-activated K+ current (ISK) in the regulation of cell function and the development of Ca2+-driven arrhythmias. We found that ISK inhibition, a promising atrial-selective pharmacological strategy against atrial fibrillation, counteracts the reentry-promoting abbreviation of atrial refractoriness, but renders human atrial myocytes more vulnerable to delayed afterdepolarizations, thus potentially increasing the propensity for ectopic (triggered) activity.

Qi-Po-Sheng-Mai granule ameliorates Ach-CaCl2 -induced atrial fibrillation by regulating calcium homeostasis in cardiomyocytes

BACKGROUND

Atrial fibrillation (AF) is one of the most common arrhythmias encountered in clinical settings. Currently, the pathophysiology of AF remains unclear, which severely limits the effectiveness and safety of medical therapies. The Chinese herbal formula Qi-Po-Sheng-Mai Granule (QPSM) has been widely used in China to treat AF. However, its pharmacological and molecular mechanisms remain unknown.

PURPOSE

The purpose of this study was to investigate the molecular mechanisms and potential targets of QPSM for AF.

STUDY DESIGN AND METHODS

The AF model was induced by Ach (66 μg/ml) and CaCl2 (10 mg/kg), and the dose of 0.1 ml/100 g was injected into the tail vein for 5 weeks. QPSM was administered daily at doses of 4.42 and 8.84 g/kg, and amiodarone (0.18 g/kg) was used as the positive control. The effect of QPSM on AF was assessed by electrocardiogram, echocardiography, and histopathological analysis. Then, we employed network pharmacology with single nucleus RNA sequencing (snRNA-Seq) to investigate the molecular mechanisms and potential targets of QPSM for AF. Furthermore, high performance liquid chromatography (HPLC) method was used for component analysis of QPSM, and molecular docking was used to verify the potential targets. Using the IonOptix single cell contraction and ion synchronization test equipment, single myocyte length and calcium ion variations were observed in real time. The expression levels of calcium Transporter-related proteins were detected by western blot and immunohistochemistry.

RESULTS

Based on an Ach-CaCl2-induced AF model, we found that QPSM treatment significantly reduced atrial electrical remodeling-related markers, such as AF inducibility and duration, and attenuated atrial dilation and fibrosis. Network pharmacology identified 52 active ingredients and 119 potential targets for QPSM in the treatment of AF, and 45 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were enriched, among which calcium pathway had the greatest impact. Using single nucleus sequencing (snRNA-seq), we identified cardiomyocytes as the most differentially expressed in response to drug treatment, with nine differentially expressed genes enriched in calcium signaling pathways. High performance liquid chromatography and molecular docking confirmed that the core components of QPSM strongly bind to the key factors in the calcium signaling pathway. Additional experiments have shown that QPSM increases calcium transients (CaT) and contractility in the individual cardiomyocyte. This was accomplished by increasing the expression of CACNA1C and SERCA2a and decreasing the expression of CAMK2B and NCX1.

CONCLUSION

The present study has systematically elucidated the role of QPSM in maintaining calcium homeostasis in cardiomyocytes through the regulation of calcium transporters, which could lead to new drug development ideas for AF.

Physiological and pathophysiological roles of the KCNK3 potassium channel in the pulmonary circulation and the heart

Potassium channel subfamily K member 3 (KCNK3), encoded by the KCNK3 gene, is part of the two-pore domain potassium channel family, constitutively active at resting membrane potentials in excitable cells, including smooth muscle and cardiac cells. Several physiological and pharmacological mediators, such as intracellular signalling pathways, extracellular pH, hypoxia and anaesthetics, regulate KCNK3 channel function. Recent studies show that modulation of KCNK3 channel expression and function strongly influences pulmonary vascular cell and cardiomyocyte function. The altered activity of KCNK3 in pathological situations such as atrial fibrillation, pulmonary arterial hypertension and right ventricular dysfunction demonstrates the crucial role of KCNK3 in cardiovascular homeostasis. Furthermore, loss of function variants of KCNK3 have been identified in patients suffering from pulmonary arterial hypertension and atrial fibrillation. This review focuses on current knowledge of the role of the KCNK3 channel in pulmonary circulation and the heart, in healthy and pathological conditions.

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

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

Pro-arrhythmic role of adrenergic spatial densities in the human atria: An in-silico study

Chronic stress among young patients (≤ 45 years old) could result in autonomic dysfunction. Autonomic dysfunction could be exhibited via sympathetic hyperactivity, sympathetic nerve sprouting, and diffuse adrenergic stimulation in the atria. Adrenergic spatial densities could alter atrial electrophysiology and increase arrhythmic susceptibility. Therefore, we examined the role of adrenergic spatial densities in creating arrhythmogenic substrates in silico. We simulated three 25 cm2 atrial sheets with varying adrenergic spatial densities (ASD), activation rates, and external transmembrane currents. We measured their effects on spatial and temporal heterogeneity of action potential durations (APD) at 50% and 20%. Increasing ASD shortens overall APD, and maximum spatial heterogeneity (31%) is achieved at 15% ASD. The addition of a few (5% to 10%) adrenergic elements decreases the excitation threshold, below 18 μA/cm2, while ASDs greater than 10% increase their excitation threshold up to 22 μA/cm2. Increase in ASD during rapid activation increases APD50 and APD20 by 21% and 41%, respectively. Activation times of captured beats during rapid activation could change by as much as 120 ms from the baseline cycle length. Rapidly activated atrial sheets with high ASDs significantly increase temporal heterogeneity of APD50 and APD20. Rapidly activated atrial sheets with 10% ASD have a high likelihood (0.7 ± 0.06) of fragmenting otherwise uniform wavefronts due to the transient inexcitability of adrenergically stimulated elements, producing an effective functional block. The likelihood of wave fragmentation due to ASD highly correlates with the spatial variations of APD20 (ρ = 0.90, p = 0.04). Our simulations provide a novel insight into the contributions of ASD to spatial and temporal heterogeneities of APDs, changes in excitation thresholds, and a potential explanation for wave fragmentation in the human atria due to sympathetic hyperactivity. Our work may aid in elucidating an electrophysiological link to arrhythmia initiation due to chronic stress among young patients.

Multiplatform modeling of atrial fibrillation identifies phospholamban as a central regulator of cardiac rhythm

Atrial fibrillation (AF) is a common and genetically inheritable form of cardiac arrhythmia; however, it is currently not known how these genetic predispositions contribute to the initiation and/or maintenance of AF-associated phenotypes. One major barrier to progress is the lack of experimental systems to investigate the effects of gene function on rhythm parameters in models with human atrial and whole-organ relevance. Here, we assembled a multi-model platform enabling high-throughput characterization of the effects of gene function on action potential duration and rhythm parameters using human induced pluripotent stem cell-derived atrial-like cardiomyocytes and a Drosophila heart model, and validation of the findings using computational models of human adult atrial myocytes and tissue. As proof of concept, we screened 20 AF-associated genes and identified phospholamban loss of function as a top conserved hit that shortens action potential duration and increases the incidence of arrhythmia phenotypes upon stress. Mechanistically, our study reveals that phospholamban regulates rhythm homeostasis by functionally interacting with L-type Ca2+ channels and NCX. In summary, our study illustrates how a multi-model system approach paves the way for the discovery and molecular delineation of gene regulatory networks controlling atrial rhythm with application to AF.

Mechanisms of spontaneous Ca2+ release-mediated arrhythmia in a novel 3D human atrial myocyte model: I. Transverse-axial tubule variation

Intracellular calcium (Ca2+ ) cycling is tightly regulated in the healthy heart ensuring effective contraction. This is achieved by transverse (t)-tubule membrane invaginations that facilitate close coupling of key Ca2+ -handling proteins such as the L-type Ca2+ channel and Na+ -Ca2+ exchanger (NCX) on the cell surface with ryanodine receptors (RyRs) on the intracellular Ca2+ store. Although less abundant and regular than in the ventricle, t-tubules also exist in atrial myocytes as a network of transverse invaginations with axial extensions known as the transverse-axial tubule system (TATS). In heart failure and atrial fibrillation, there is TATS remodelling that is associated with aberrant Ca2+ -handling and Ca2+ -induced arrhythmic activity; however, the mechanism underlying this is not fully understood. To address this, we developed a novel 3D human atrial myocyte model that couples electrophysiology and Ca2+ -handling with variable TATS organization and density. We extensively parameterized and validated our model against experimental data to build a robust tool examining TATS regulation of subcellular Ca2+ release. We found that varying TATS density and thus the localization of key Ca2+ -handling proteins has profound effects on Ca2+ handling. Following TATS loss, there is reduced NCX that results in increased cleft Ca2+ concentration through decreased Ca2+ extrusion. This elevated Ca2+ increases RyR open probability causing spontaneous Ca2+ releases and the promotion of arrhythmogenic waves (especially in the cell interior) leading to voltage instabilities through delayed afterdepolarizations. In summary, the present study demonstrates a mechanistic link between TATS remodelling and Ca2+ -driven proarrhythmic behaviour that probably reflects the arrhythmogenic state observed in disease. KEY POINTS: Transverse-axial tubule systems (TATS) modulate Ca2+ handling and excitation-contraction coupling in atrial myocytes, with TATS remodelling in heart failure and atrial fibrillation being associated with altered Ca2+ cycling and subsequent arrhythmogenesis. To investigate the poorly understood mechanisms linking TATS variation and spontaneous Ca2+ release, we built, parameterized and validated a 3D human atrial myocyte model coupling electrophysiology and spatially-detailed subcellular Ca2+ handling governed by the TATS. Simulated TATS loss causes diastolic Ca2+ and voltage instabilities through reduced Na+ -Ca2+ exchanger-mediated Ca2+ removal, cleft Ca2+ accumulation and increased ryanodine receptor open probability, resulting in spontaneous Ca2+ release and promotion of arrhythmogenic waves and delayed afterdepolarizations. At fast electrical rates typical of atrial tachycardia/fibrillation, spontaneous Ca2+ releases are larger and more frequent in the cell interior than at the periphery. Our work provides mechanistic insight into how atrial TATS remodelling can lead to Ca2+ -driven instabilities that may ultimately contribute to the arrhythmogenic state in disease.

Human atrial fibrillation and genetic defects in transient outward currents: mechanistic insights from multi-scale computational models

Previous studies have linked dysfunctional I_to arising from mutations to KCND3-encoded Kv4.3 and KCND2-encoded Kv4.2 to atrial fibrillation. Using computational models, this study aimed to investigate the mechanisms underlying pro-arrhythmic effects of the gain-of-function Kv4.3 (T361S, A545P) and Kv4.2 (S447R) mutations. Wild-type and mutant I_to formulations were developed from and validated against experimental data and incorporated into the Colman et al. model of human atrial cells. Single-cell models were incorporated into one- (1D) and two-dimensional (2D) models of atrial tissue, and a three-dimensional (3D) realistic model of the human atria. The three gain-of-function mutations had similar, albeit quantitatively different, effects: shortening of the action potential duration; lowering the plateau membrane potential, abbreviating the effective refractory period (ERP) and the wavelength (WL) of atrial excitation at the tissue level. Restitution curves for the WL, the ERP and the conduction velocity were leftward shifted, facilitating the conduction of atrial excitation waves at high excitation rates. The mutations also increased lifespan and stationarity of re-entry in both 2D and 3D simulations, which further highlighted a mutation-induced increase in spatial dispersion of repolarization. Collectively, these changes account for pro-arrhythmic effects of these Kv4.3 and Kv4.2 mutations in facilitating AF. This article is part of the theme issue 'The heartbeat: its molecular basis and physiological mechanisms'.

Timing mechanisms to control heart rhythm and initiate arrhythmias: roles for intracellular organelles, signalling pathways and subsarcolemmal Ca2

Rhythms of electrical activity in all regions of the heart can be influenced by a variety of intracellular membrane bound organelles. This is true both for normal pacemaker activity and for abnormal rhythms including those caused by early and delayed afterdepolarizations under pathological conditions. The influence of the sarcoplasmic reticulum (SR) on cardiac electrical activity is widely recognized, but other intracellular organelles including lysosomes and mitochondria also contribute. Intracellular organelles can provide a timing mechanism (such as an SR clock driven by cyclic uptake and release of Ca²⁺, with an important influence of intraluminal Ca²⁺), and/or can act as a Ca²⁺ store involved in signalling mechanisms. Ca²⁺ plays many diverse roles including carrying electric current, driving electrogenic sodium-calcium exchange (NCX) particularly when Ca²⁺ is extruded across the surface membrane causing depolarization, and activation of enzymes which target organelles and surface membrane proteins. Heart function is also influenced by Ca²⁺ mobilizing agents (cADP-ribose, nicotinic acid adenine dinucleotide phosphate and inositol trisphosphate) acting on intracellular organelles. Lysosomal Ca²⁺ release exerts its effects via calcium/calmodulin-dependent protein kinase II to promote SR Ca²⁺ uptake, and contributes to arrhythmias resulting from excessive beta-adrenoceptor stimulation. A separate arrhythmogenic mechanism involves lysosomes, mitochondria and SR. Interacting intracellular organelles, therefore, have profound effects on heart rhythms and NCX plays a central role. This article is part of the theme issue 'The heartbeat: its molecular basis and physiological mechanisms'.

Generation of cardiomyocytes from human-induced pluripotent stem cells resembling atrial cells with ability to respond to adrenoceptor agonists

Atrial fibrillation (AF) is the most common chronic arrhythmia presenting a heavy disease burden. We report a new approach for generating cardiomyocytes (CMs) resembling atrial cells from human-induced pluripotent stem cells (hiPSCs) using a combination of Gremlin 2 and retinoic acid treatment. More than 40% of myocytes showed rod-shaped morphology, expression of CM proteins (including ryanodine receptor 2, α-actinin-2 and F-actin) and striated appearance, all of which were broadly similar to the characteristics of adult atrial myocytes (AMs). Isolated myocytes were electrically quiescent until stimulated to fire action potentials with an AM profile and an amplitude of approximately 100 mV, arising from a resting potential of approximately -70 mV. Single-cell RNA sequence analysis showed a high level of expression of several atrial-specific transcripts including NPPA, MYL7, HOXA3, SLN, KCNJ4, KCNJ5 and KCNA5. Amplitudes of calcium transients recorded from spontaneously beating cultures were increased by the stimulation of α-adrenoceptors (activated by phenylephrine and blocked by prazosin) or β-adrenoceptors (activated by isoproterenol and blocked by CGP20712A). Our new approach provides human AMs with mature characteristics from hiPSCs which will facilitate drug discovery by enabling the study of human atrial cell signalling pathways and AF. This article is part of the theme issue 'The heartbeat: its molecular basis and physiological mechanisms'.

Enhanced Ca2+-Dependent SK-Channel Gating and Membrane Trafficking in Human Atrial Fibrillation

BACKGROUND

Small-conductance Ca2+-activated K+ (SK)-channel inhibitors have antiarrhythmic effects in animal models of atrial fibrillation (AF), presenting a potential novel antiarrhythmic option. However, the regulation of SK-channels in human atrial cardiomyocytes and its modification in patients with AF are poorly understood and were the object of this study.

METHODS

Apamin-sensitive SK-channel current (ISK) and action potentials were recorded in human right-atrial cardiomyocytes from sinus rhythm control (Ctl) patients or patients with (long-standing persistent) chronic AF (cAF).

RESULTS

ISK was significantly higher, and apamin caused larger action potential prolongation in cAF- versus Ctl-cardiomyocytes. Sensitivity analyses in an in silico human atrial cardiomyocyte model identified IK1 and ISK as major regulators of repolarization. Increased ISK in cAF was not associated with increases in mRNA/protein levels of SK-channel subunits in either right- or left-atrial tissue homogenates or right-atrial cardiomyocytes, but the abundance of SK2 at the sarcolemma was larger in cAF versus Ctl in both tissue-slices and cardiomyocytes. Latrunculin-A and primaquine (anterograde and retrograde protein-trafficking inhibitors) eliminated the differences in SK2 membrane levels and ISK between Ctl- and cAF-cardiomyocytes. In addition, the phosphatase-inhibitor okadaic acid reduced ISK amplitude and abolished the difference between Ctl- and cAF-cardiomyocytes, indicating that reduced calmodulin-Thr80 phosphorylation due to increased protein phosphatase-2A levels in the SK-channel complex likely contribute to the greater ISK in cAF-cardiomyocytes. Finally, rapid electrical activation (5 Hz, 10 minutes) of Ctl-cardiomyocytes promoted SK2 membrane-localization, increased ISK and reduced action potential duration, effects greatly attenuated by apamin. Latrunculin-A or primaquine prevented the 5-Hz-induced ISK-upregulation.

CONCLUSIONS

ISK is upregulated in patients with cAF due to enhanced channel function, mediated by phosphatase-2A-dependent calmodulin-Thr80 dephosphorylation and tachycardia-dependent enhanced trafficking and targeting of SK-channel subunits to the sarcolemma. The observed AF-associated increases in ISK, which promote reentry-stabilizing action potential duration shortening, suggest an important role for SK-channels in AF auto-promotion and provide a rationale for pursuing the antiarrhythmic effects of SK-channel inhibition in humans.

Steady-state and transient effects of SK channel block and adrenergic stimulation to counteract acetylcholine-induced arrhythmogenic effects in the human atria: A computational study

Hyperactivity of the parasympathetic nervous system has been linked to the development of paroxysmal atrial fibrillation (AF). The parasympathetic neurotransmitter acetylcholine (ACh) causes a reduction in action potential (AP) duration (APD) and an increase in resting membrane potential (RMP), both of which contribute to enhance the risk for reentry. Research suggests that small-conductance calcium activated potassium (SK) channels may be an effective target for treating AF. Therapies targeting the autonomic nervous system, either alone or in combination with other drugs, have been explored and have been shown to decrease the incidence of atrial arrhythmias. This study uses computational modeling and simulation to examine the impact of SK channel block (SKb) and β-adrenergic stimulation through Isoproterenol (Iso) on countering the negative effects of cholinergic activity in human atrial cell and 2D tissue models. The steady-state effects of Iso and/or SKb on AP shape, APD at 90% repolarization (APD₉₀) and RMP were evaluated. The ability to terminate stable rotational activity in cholinergically-stimulated 2D tissue models of AF was also investigated. A range of SKb and Iso application kinetics, which reflect varying drug binding rates, were taken into consideration. The results showed that SKb alone prolonged APD₉₀ and was able to stop sustained rotors in the presence of ACh concentrations up to 0.01 μM. Iso terminated rotors under all tested ACh concentrations, but resulted in highly-variable steady-state outcomes depending on baseline AP morphology. Importantly, the combination of SKb and Iso resulted in greater APD₉₀ prolongation and showed promising anti-arrhythmic potential by stopping stable rotors and preventing re-inducibility.

Models of the cardiac L-type calcium current: A quantitative review

The L-type calcium current (I CaL ) plays a critical role in cardiac electrophysiology, and models ofI CaL are vital tools to predict arrhythmogenicity of drugs and mutations. Five decades of measuring and modelingI CaL have resulted in several competing theories (encoded in mathematical equations). However, the introduction of new models has not typically been accompanied by a data-driven critical comparison with previous work, so that it is unclear which model is best suited for any particular application. In this review, we describe and compare 73 published mammalianI CaL models and use simulated experiments to show that there is a large variability in their predictions, which is not substantially diminished when grouping by species or other categories. We provide model code for 60 models, list major data sources, and discuss experimental and modeling work that will be required to reduce this huge list of competing theories and ultimately develop a community consensus model ofI CaL . This article is categorized under: Cardiovascular Diseases > Computational Models Cardiovascular Diseases > Molecular and Cellular Physiology.

Pacemaker activity and ion channels in the sinoatrial node cells: MicroRNAs and arrhythmia

The primary pacemaking activity of the heart is determined by a spontaneous action potential (AP) within sinoatrial node (SAN) cells. This unique AP generation relies on two mechanisms: membrane clocks and calcium clocks. Nonhomologous arrhythmias are caused by several functional and structural changes in the myocardium. MicroRNAs (miRNAs) are essential regulators of gene expression in cardiomyocytes. These miRNAs play a vital role in regulating the stability of cardiac conduction and in the remodeling process that leads to arrhythmias. Although it remains unclear how miRNAs regulate the expression and function of ion channels in the heart, these regulatory mechanisms may support the development of emerging therapies. This study discusses the spread and generation of AP in the SAN as well as the regulation of miRNAs and individual ion channels. Arrhythmogenicity studies on ion channels will provide a research basis for miRNA modulation as a new therapeutic target.

Sensitivity Analysis of Cardiac Alternans and Tachyarrhythmia to Ion Channel Conductance Using Population Modeling

Action potential duration (APD) alternans, an alternating phenomenon between action potentials in cardiomyocytes, causes heart arrhythmia when the heart rate is high. However, some of the APD alternans observed in clinical trials occurs under slow heart rate conditions of 100 to 120 bpm, increasing the likelihood of heart arrhythmias such as atrial fibrillation. Advanced studies have identified the occurrence of this type of APD alternans in terms of electrophysiological ion channel currents in cells. However, they only identified physiological phenomena, such as action potential due to random changes in a particular ion channel’s conductivity through ion models specializing in specific ion channel currents. In this study, we performed parameter sensitivity analysis via population modeling using a validated human ventricular physiology model to check the sensitivity of APD alternans to ion channel conductances. Through population modeling, we expressed the changes in alternans onset cycle length (AOCL) and mean APD in AOCL (AO meanAPD) according to the variations in ion channel conductance. Finally, we identified the ion channel that maximally affected the occurrence of APD alternans. AOCL and AO meanAPD were sensitive to changes in the plateau Ca²⁺ current. Accordingly, it was expected that APD alternans would be vulnerable to changes in intracellular calcium concentration.

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.

In-silico drug trials for precision medicine in atrial fibrillation: From ionic mechanisms to electrocardiogram-based predictions in structurally-healthy human atria

Atrial fibrillation (AF) inducibility, sustainability and response to pharmacological treatment of individual patients are expected to be determined by their ionic current properties, especially in structurally-healthy atria. Mechanisms underlying AF and optimal cardioversion are however still unclear. In this study, in-silico drug trials were conducted using a population of human structurally-healthy atria models to 1) identify key ionic current properties determining AF inducibility, maintenance and pharmacological cardioversion, and 2) compare the prognostic value for predicting individual AF cardioversion of ionic current properties and electrocardiogram (ECG) metrics. In the population of structurally-healthy atria, 477 AF episodes were induced in ionic current profiles with both steep action potential duration (APD) restitution (eliciting APD alternans), and high excitability (enabling propagation at fast rates that transformed alternans into discordant). High excitability also favored 211 sustained AF episodes, so its decrease, through prolonged refractoriness, explained pharmacological cardioversion. In-silico trials over 200 AF episodes, 100 ionic profiles and 10 antiarrhythmic compounds were consistent with previous clinical trials, and identified optimal treatments for individual electrophysiological properties of the atria. Algorithms trained on 211 simulated AF episodes exhibited >70% accuracy in predictions of cardioversion for individual treatments using either ionic current profiles or ECG metrics. In structurally-healthy atria, AF inducibility and sustainability are enabled by discordant alternans, under high excitability and steep restitution conditions. Successful pharmacological cardioversion is predicted with 70% accuracy from either ionic or ECG properties, and it is optimal for treatments maximizing refractoriness (thus reducing excitability) for the given ionic current profile of the atria.

Interpretable machine learning of action potential duration restitution kinetics in single-cell models of atrial cardiomyocytes

Action potential duration (APD) restitution curve and its maximal slope (Smax) reflect single cell-level dynamic instability for inducing chaotic heart rhythms. However, conventional parameter sensitivity analysis often fails to describe nonlinear relationships between ion channel parameters and electrophysiological phenotypes, such as Smax. We explored the parameter-phenotype mapping in a population of 5000 single-cell atrial cell models through interpretable machine learning (ML) approaches. Parameter sensitivity analyses could explain the linear relationships between parameters and electrophysiological phenotypes, including APD₉₀, resting membrane potential, Vmax, refractory period, and APD/calcium alternans threshold, but not for Smax. However, neural network models had better prediction performance for Smax. To interpret the ML model, we evaluated the parameter importance at the global and local levels by computing the permutation feature importance and the local interpretable model-agnostic explanations (LIME) values, respectively. Increases in I_CaL, I_NCX, and I_Kr, and decreases in I_K1, I_b,Cl, I_Kur, I_SERCA, and I_to are correlated with higher Smax values. The LIME algorithm determined that I_NaK plays a significant role in determining Smax as well as I_to and I_Kur. The atrial cardiomyocyte population was hierarchically clustered into three distinct groups based on the LIME values and the single-cell simulation confirmed that perturbations in I_NaK resulted in different behaviors of APD restitution curves in three clusters. Our combined top-down interpretable ML and bottom-up mechanistic simulation approaches uncovered the role of I_NaK in heterogeneous behaviors of Smax in the atrial cardiomyocyte population.

Ablation and antiarrhythmic drug effects on PITX2+/− deficient atrial fibrillation: A computational modeling study

Introduction

Atrial fibrillation (AF) is a heritable disease, and the paired-like homeodomain transcription factor 2 (PITX2) gene is highly associated with AF. We explored the differences in the circumferential pulmonary vein isolation (CPVI), which is the cornerstone procedure for AF catheter ablation, additional high dominant frequency (DF) site ablation, and antiarrhythmic drug (AAD) effects according to the patient genotype (wild-type and PITX2+/− deficient) using computational modeling.

Methods

We included 25 patients with AF (68% men, 59.8 ± 9.8 years of age, 32% paroxysmal AF) who underwent AF catheter ablation to develop a realistic computational AF model. The ion currents for baseline AF and the amiodarone, dronedarone, and flecainide AADs according to the patient genotype (wild type and PITX2+/− deficient) were defined by relevant publications. We tested the virtual CPVI (V-CPVI) with and without DF ablation (±DFA) and three virtual AADs (V-AADs, amiodarone, dronedarone, and flecainide) and evaluated the AF defragmentation rates (AF termination or changes to regular atrial tachycardia (AT), DF, and maximal slope of the action potential duration restitution curves (Smax), which indicates the vulnerability of wave-breaks.

Results

At the baseline AF, mean DF (p = 0.003), and Smax (p &lt; 0.001) were significantly lower in PITX2+/− deficient patients than wild-type patients. In the overall AF episodes, V-CPVI (±DFA) resulted in a higher AF defragmentation relative to V-AADs (65 vs. 42%, p &lt; 0.001) without changing the DF or Smax. Although a PITX2+/− deficiency did not affect the AF defragmentation rate after the V-CPVI (±DFA), V-AADs had a higher AF defragmentation rate (p = 0.014), lower DF (p &lt; 0.001), and lower Smax (p = 0.001) in PITX2+/− deficient AF than in wild-type patients. In the clinical setting, the PITX2+/− genetic risk score did not affect the AF ablation rhythm outcome (Log-rank p = 0.273).

Conclusion

Consistent with previous clinical studies, the V-CPVI had effective anti-AF effects regardless of the PITX2 genotype, whereas V-AADs exhibited more significant defragmentation or wave-dynamic change in the PITX2+/− deficient patients.

SNTA1 gene rescues ion channel function and is antiarrhythmic in cardiomyocytes derived from induced pluripotent stem cells from muscular dystrophy patients

Background

Patients with cardiomyopathy of Duchenne Muscular Dystrophy (DMD) are at risk of developing life-threatening arrhythmias, but the mechanisms are unknown. We aimed to determine the role of ion channels controlling cardiac excitability in the mechanisms of arrhythmias in DMD patients.

Methods

To test whether dystrophin mutations lead to defective cardiac NaV1.5-Kir2.1 channelosomes and arrhythmias, we generated iPSC-CMs from two hemizygous DMD males, a heterozygous female, and two unrelated control males. We conducted studies including confocal microscopy, protein expression analysis, patch-clamping, non-viral piggy-bac gene expression, optical mapping and contractility assays.

Results

Two patients had abnormal ECGs with frequent runs of ventricular tachycardia. iPSC-CMs from all DMD patients showed abnormal action potential profiles, slowed conduction velocities, and reduced sodium (INa) and inward rectifier potassium (IK1) currents. Membrane NaV1.5 and Kir2.1 protein levels were reduced in hemizygous DMD iPSC-CMs but not in heterozygous iPSC-CMs. Remarkably, transfecting just one component of the dystrophin protein complex (α1-syntrophin) in hemizygous iPSC-CMs from one patient restored channelosome function, INa and IK1 densities, and action potential profile in single cells. In addition, α1-syntrophin expression restored impulse conduction and contractility and prevented reentrant arrhythmias in hiPSC-CM monolayers.

Conclusions

We provide the first demonstration that iPSC-CMs reprogrammed from skin fibroblasts of DMD patients with cardiomyopathy have a dysfunction of the NaV1.5-Kir2.1 channelosome, with consequent reduction of cardiac excitability and conduction. Altogether, iPSC-CMs from patients with DMD cardiomyopathy have a NaV1.5-Kir2.1 channelosome dysfunction, which can be rescued by the scaffolding protein α1-syntrophin to restore excitability and prevent arrhythmias.

Funding

Supported by National Institutes of Health R01 HL122352 grant; 'la Caixa' Banking Foundation (HR18-00304); Fundación La Marató TV3: Ayudas a la investigación en enfermedades raras 2020 (LA MARATO-2020); Instituto de Salud Carlos III/FEDER/FSE; Horizon 2020 - Research and Innovation Framework Programme GA-965286 to JJ; the CNIC is supported by the Instituto de Salud Carlos III (ISCIII), the Ministerio de Ciencia e Innovación (MCIN) and the Pro CNIC Foundation), and is a Severo Ochoa Center of Excellence (grant CEX2020-001041-S funded by MICIN/AEI/10.13039/501100011033). American Heart Association postdoctoral fellowship 19POST34380706s to JVEN. Israel Science Foundation to OB and MA [824/19]. Rappaport grant [01012020RI]; and Niedersachsen Foundation [ZN3452] to OB; US-Israel Binational Science Foundation (BSF) to OB and TH [2019039]; Dr. Bernard Lublin Donation to OB; and The Duchenne Parent Project Netherlands (DPPNL 2029771) to OB. National Institutes of Health R01 AR068428 to DM and US-Israel Binational Science Foundation Grant [2013032] to DM and OB.

Arrhythmogenic influence of mutations in a myocyte-based computational model of the pulmonary vein sleeve

In the heart, electrophysiological dysregulation arises from defects at many biological levels (from point mutations in ion channel proteins to gross structural abnormalities). These defects disrupt the normal pattern of electrical activation, producing ectopic activity and reentrant arrhythmia. To interrogate mechanisms that link these primary biological defects to macroscopic electrophysiologic dysregulation most prior computational studies have utilized either (i) detailed models of myocyte ion channel dynamics at limited spatial scales, or (ii) homogenized models of action potential conduction that reproduce arrhythmic activity at tissue and organ levels. Here we apply our recent model (EMI), which integrates electrical activation and propagation across these scales, to study human atrial arrhythmias originating in the pulmonary vein (PV) sleeves. These small structures initiate most supraventricular arrhythmias and include pronounced myocyte-to-myocyte heterogeneities in ion channel expression and intercellular coupling. To test EMI's cell-based architecture in this physiological context we asked whether ion channel mutations known to underlie atrial fibrillation are capable of initiating arrhythmogenic behavior via increased excitability or reentry in a schematic PV sleeve geometry. Our results illustrate that EMI's improved spatial resolution can directly interrogate how electrophysiological changes at the individual myocyte level manifest in tissue and as arrhythmia in the PV sleeve.

Anti-atrial Fibrillation Effects of Pulmonary Vein Isolation With or Without Ablation Gaps: A Computational Modeling Study

Background

Although pulmonary vein isolation (PVI) gaps contribute to recurrence after atrial fibrillation (AF) catheter ablation, the mechanism is unclear. We used realistic computational human AF modeling to explore the AF wave-dynamic changes of PVI with gaps (PVI-gaps).

Methods

We included 40 patients (80% male, 61.0 ± 9.8 years old, 92.5% persistent AF) who underwent AF catheter ablation to develop our realistic computational AF model. We compared the effects of a complete PVI (CPVI) and PVI-gap (2-mm × 4) on the AF wave-dynamics by evaluating the dominant frequency (DF), spatial change of DF, maximal slope of the action potential duration restitution curve (Smax), and AF defragmentation rate (termination or change to atrial tachycardia), and tested the effects of additional virtual interventions and flecainide on ongoing AF with PVI-gaps.

Results

Compared with the baseline AF, CPVIs significantly reduced extra-PV DFs (p < 0.001), but PVI-gaps did not. COV-DFs were greater after CPVIs than PVI-gaps (p < 0.001). Neither CPVIs nor PVI-gaps changed the mean Smax. CPVIs resulted in higher AF defragmentation rates (80%) than PVI-gaps (12.5%, p < 0.001). In ongoing AF after PVI-gaps, the AF defragmentation rates after a wave-breaking gap ablation, extra-PV DF ablation, or flecainide were 60.0, 34.3, and 25.7%, respectively (p = 0.010).

Conclusion

CPVIs effectively reduced the DF, increased its spatial heterogeneity in extra-PV areas, and offered better anti-AF effects than extra-PV DF ablation or additional flecainide in PVI-gap conditions.

A Review on Atrial Fibrillation (Computer Simulation and Clinical Perspectives)

Atrial fibrillation (AF), a heart condition, has been a well-researched topic for the past few decades. This multidisciplinary field of study deals with signal processing, finite element analysis, mathematical modeling, optimization, and clinical procedure. This article is focused on a comprehensive review of journal articles published in the field of AF. Topics from the age-old fundamental concepts to specialized modern techniques involved in today’s AF research are discussed. It was found that a lot of research articles have already been published in modeling and simulation of AF. In comparison to that, the diagnosis and post-operative procedures for AF patients have not yet been totally understood or explored by the researchers. The simulation and modeling of AF have been investigated by many researchers in this field. Cellular model, tissue model, and geometric model among others have been used to simulate AF. Due to a very complex nature, the causes of AF have not been fully perceived to date, but the simulated results are validated with real-life patient data. Many algorithms have been proposed to detect the source of AF in human atria. There are many ablation strategies for AF patients, but the search for more efficient ablation strategies is still going on. AF management for patients with different stages of AF has been discussed in the literature as well but is somehow limited mostly to the patients with persistent AF. The authors hope that this study helps to find existing research gaps in the analysis and the diagnosis of AF.

Basic Research Approaches to Evaluate Cardiac Arrhythmia in Heart Failure and Beyond

Patients with heart failure often develop cardiac arrhythmias. The mechanisms and interrelations linking heart failure and arrhythmias are not fully understood. Historically, research into arrhythmias has been performed on affected individuals or in vivo (animal) models. The latter however is constrained by interspecies variation, demands to reduce animal experiments and cost. Recent developments in in vitro induced pluripotent stem cell technology and in silico modelling have expanded the number of models available for the evaluation of heart failure and arrhythmia. An agnostic approach, combining the modalities discussed here, has the potential to improve our understanding for appraising the pathology and interactions between heart failure and arrhythmia and can provide robust and validated outcomes in a variety of research settings. This review discusses the state of the art models, methodologies and techniques used in the evaluation of heart failure and arrhythmia and will highlight the benefits of using them in combination. Special consideration is paid to assessing the pivotal role calcium handling has in the development of heart failure and arrhythmia.

Voltage-mediated mechanism for calcium wave synchronization and arrhythmogenesis in atrial tissue

A wide range of atrial arrythmias are caused by molecular defects in proteins that regulate calcium (Ca) cycling. In many cases, these defects promote the propagation of subcellular Ca waves in the cell, which can perturb the voltage time course and induce dangerous perturbations of the action potential (AP). However, subcellular Ca waves occur randomly in cells and, therefore, electrical coupling between cells substantially decreases their effect on the AP. In this study, we present evidence that Ca waves in atrial tissue can synchronize in-phase owing to an order-disorder phase transition. In particular, we show that, below a critical pacing rate, Ca waves are desynchronized and therefore do not induce substantial AP fluctuations in tissue. However, above this critical pacing rate, Ca waves gradually synchronize over millions of cells, which leads to a dramatic amplification of AP fluctuations. We exploit an underlying Ising symmetry of paced cardiac tissue to show that this transition exhibits universal properties common to a wide range of physical systems in nature. Finally, we show that in the heart, phase synchronization induces spatially out-of-phase AP duration alternans which drives wave break and reentry. These results suggest that cardiac tissue exhibits a phase transition that is required for subcellular Ca cycling defects to induce a life-threatening arrhythmia.

Heart-on-Chip for Combined Cellular Dynamics Measurements and Computational Modeling Towards Clinical Applications

Organ-on-chip or micro-engineered three-dimensional cellular or tissue models are increasingly implemented in the study of cardiovascular pathophysiology as alternatives to traditional in vitro cell culture. Drug induced cardiotoxicity is a key issue in drug development pipelines, but the current in vitro and in vivo studies suffer from inter-species differences, high costs, and lack of reliability and accuracy in predicting cardiotoxicity. Microfluidic heart-on-chip devices can impose a paradigm shift to the current tools. They can not only recapitulate cardiac tissue level functionality and the communication between cells and extracellular matrices but also allow higher throughput studies conducive to drug screening especially with their added functionalities or sensors that extract disease-specific phenotypic, genotypic, and electrophysiological information in real-time. Such electrical and mechanical components can tailor the electrophysiology and mechanobiology of the experiment to better mimic the in vivo condition as well. Recent advancements and challenges are reviewed in the fabrication, functionalization and sensor assisted mechanical and electrophysiological measurements, numerical and computational modeling of cardiomyocytes' behavior, and the clinical applications in drug screening and disease modeling. This review concludes with the current challenges and perspectives on the future of such organ-on-chip platforms.

Individual Contributions of Cardiac Ion Channels on Atrial Repolarization and Reentrant Waves: A Multiscale In-Silico Study

The excitation, contraction, and relaxation of an atrial cardiomyocyte are maintained by the activation and inactivation of numerous cardiac ion channels. Their collaborative efforts cause time-dependent changes of membrane potential, generating an action potential (AP), which is a surrogate marker of atrial arrhythmias. Recently, computational models of atrial electrophysiology emerged as a modality to investigate arrhythmia mechanisms and to predict the outcome of antiarrhythmic therapies. However, the individual contribution of atrial ion channels on atrial action potential and reentrant arrhythmia is not yet fully understood. Thus, in this multiscale in-silico study, perturbations of individual atrial ionic currents (I_Na, I_to, I_CaL, I_Kur, I_Kr, I_Ks, I_K1, I_NCX and I_NaK) in two in-silico models of human atrial cardiomyocyte (i.e., Courtemanche-1998 and Grandi-2011) were performed at both cellular and tissue levels. The results show that the inhibition of I_CaL and I_NCX resulted in AP shortening, while the inhibition of I_Kur, I_Kr, I_Ks, I_K1 and I_NaK prolonged AP duration (APD). Particularly, in-silico perturbations (inhibition and upregulation) of I_Kr and I_Ks only minorly affected atrial repolarization in the Grandi model. In contrast, in the Courtemanche model, the inhibition of I_Kr and I_Ks significantly prolonged APD and vice versa. Additionally, a 50% reduction of I_to density abbreviated APD in the Courtemanche model, while the same perturbation prolonged APD in the Grandi model. Similarly, a strong model dependence was also observed at tissue scale, with an observable I_K1-mediated reentry stabilizing effect in the Courtemanche model but not in the Grandi atrial model. Moreover, the Grandi model was highly sensitive to a change on intracellular Ca²⁺ concentration, promoting a repolarization failure in I_CaL upregulation above 150% and facilitating reentrant spiral waves stabilization by I_CaL inhibition. Finally, by incorporating the previously published atrial fibrillation (AF)-associated ionic remodeling in the Courtemanche atrial model, in-silico modeling revealed the antiarrhythmic effect of I_Kr inhibition in both acute and chronic settings. Overall, our multiscale computational study highlights the strong model-dependent effects of ionic perturbations which could affect the model’s accuracy, interpretability, and prediction. This observation also suggests the need for a careful selection of in-silico models of atrial electrophysiology to achieve specific research aims.

Clinical Spectrum of SCN5A Channelopathy in Children with Primary Electrical Disease and Structurally Normal Hearts

Sodium voltage-gated channel α subunit 5 (SCN5A)-mutations may cause an array of arrhythmogenic syndromes most frequently as an autosomal dominant trait, with incomplete penetrance, variable expressivity and male predominance. In the present study, we retrospectively describe a group of Mexican patients with SCN5A-disease causing variants in whom the onset of symptoms occurred in the pediatric age range. The study included 17 patients with clinical diagnosis of primary electrical disease, at least one SCN5A pathogenic or likely pathogenic mutation and age of onset <18 years, and all available first- and second-degree relatives. Fifteen patients (88.2%) were male, and sixteen independent variants were found (twelve missense, three truncating and one complex inframe deletion/insertion). The frequency of compound heterozygosity was remarkably high (3/17, 17.6%), with early childhood onset and severe disease. Overall, 70.6% of pediatric patients presented with overlap syndrome, 11.8% with isolated sick sinus syndrome, 11.8% with isolated Brugada syndrome (BrS) and 5.9% with isolated type 3 long QT syndrome (LQTS). A total of 24/45 SCN5A mutation carriers were affected (overall penetrance 53.3%), and penetrance was higher in males (63.3%, 19 affected/30 mutation carriers) than in females (33.3%, 5 affected/15 carriers). In conclusion, pediatric patients with SCNA-disease causing variants presented mainly as overlap syndrome, with predominant loss-of-function phenotypes of sick sinus syndrome (SSS), progressive cardiac conduction disease (PCCD) and ventricular arrhythmias.

Human Atrial Cardiac Microtissues for Chamber-Specific Arrhythmic Risk Assessment

Introduction

Although atrial fibrillation is the most prevalent disorder of electrical conduction, the mechanisms behind atrial arrhythmias remain elusive. To address this challenge, we developed a robust in vitro model of 3D atrial microtissue from human induced pluripotent stem cell (hiPSC)-derived cardiomyocytes and evaluated chamber-specific chemical responses experimentally and computationally.

Methods

We differentiated atrial and ventricular cardiomyocytes (aCMs/vCMs) from GCaMP6f-expressing hiPSCs and assessed spontaneous AP activity using fluorescence imaging. Self-assembling 3D microtissues were formed with lactate purified CMs and 5% human cardiac fibroblasts and electrically stimulated for one week before high resolution action potential (AP) optical mapping. AP responses to the atrial-specific potassium repolarizing current I _Kur-blocker 4-Aminopyridine (4-AP) and funny current I _f-blocker Ivabradine were characterized within their therapeutic window. Finally, we expanded upon a published hiPSC-CM computational model by incorporating the atrial-specific I _Kur current, modifying ion channel conductances to match the AP waveforms of our microtissues, and employing the updated model to reinforce our experimental findings.

Results

High purity CMs (> 75% cTnT⁺) demonstrated subtype specification by MLC2v expression. Spontaneous beating rates significantly decreased following 3D microtissue formation, with atrial microtissues characterized by their faster spontaneous beating rate, slower AP rise time, and shorter AP duration (APD) compared to ventricular microtissues. We measured atrial-specific responses, including dose-dependent APD prolongation with 4-AP treatment and dose-dependent reduction in spontaneous activity post-Ivabradine treatment.

Conclusion

The presented in vitro platform for screening atrial-specific responses is both robust and sensitive, with high throughput, enabling studies focused at elucidating the mechanisms underlying atrial arrhythmias.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12195-021-00703-x.

Spatial Changes in the Atrial Fibrillation Wave-Dynamics After Using Antiarrhythmic Drugs: A Computational Modeling Study

Background: We previously reported that a computational modeling-guided antiarrhythmic drug (AAD) test was feasible for evaluating multiple AADs in patients with atrial fibrillation (AF). We explored the anti-AF mechanisms of AADs and spatial change in the AF wave-dynamics by a realistic computational model. Methods: We used realistic computational modeling of 25 AF patients (68% male, 59.8 ± 9.8 years old, 32.0% paroxysmal AF) reflecting the anatomy, histology, and electrophysiology of the left atrium (LA) to characterize the effects of five AADs (amiodarone, sotalol, dronedarone, flecainide, and propafenone). We evaluated the spatial change in the AF wave-dynamics by measuring the mean dominant frequency (DF) and its coefficient of variation [dominant frequency-coefficient of variation (DF-COV)] in 10 segments of the LA. The mean DF and DF-COV were compared according to the pulmonary vein (PV) vs. extra-PV, maximal slope of the restitution curves (Smax), and defragmentation of AF. Results: The mean DF decreased after the administration of AADs in the dose dependent manner (p < 0.001). Under AADs, the DF was significantly lower (p < 0.001) and COV-DF higher (p = 0.003) in the PV than extra-PV region. The mean DF was significantly lower at a high Smax (≥1.4) than a lower Smax condition under AADs. During the episodes of AF defragmentation, the mean DF was lower (p < 0.001), but the COV-DF was higher (p < 0.001) than that in those without defragmentation. Conclusions: The DF reduction with AADs is predominant in the PVs and during a high Smax condition and causes AF termination or defragmentation during a lower DF and spatially unstable (higher DF-COV) condition.

Recent Advances in Maturation of Pluripotent Stem Cell-Derived Cardiomyocytes Promoted by Mechanical Stretch

Stem cells have significant potential use in tissue regeneration, especially for treating cardiac diseases because of their multi-directional differentiation capability. By mimicking the in vivo physiological environment of native cardiomyocytes during their development and maturation, researchers have been able to induce pluripotent stem cell-derived cardiomyocytes (PSC-CMs) at high purity. However, the phenotype of these PSC-CMs is immature compared with that of adult cardiomyocytes. Various strategies have been explored to improve the maturity of PSC-CMs, such as long-term culturing, mechanical stimuli, chemical stimuli, and combinations of these strategies. Among these strategies, mechanical stretch as a key mechanical stimulus plays an important role in PSC-CM maturation. In this review, the optimal parameters of mechanical stretch, the effects of mechanical stretch on maturation of PSC-CMs, underlying molecular mechanisms as well as existing problems are discussed. Mechanical stretch is a powerful approach to promote the maturation of SC-CMs in terms of morphology, structure, and functionality. Nonetheless, further research efforts are needed to reach a satisfactory standard for clinical applications of PSC-CMs in treating cardiac diseases.

Neural Network Differential Equations For Ion Channel Modelling

Mathematical models of cardiac ion channels have been widely used to study and predict the behaviour of ion currents. Typically models are built using biophysically-based mechanistic principles such as Hodgkin-Huxley or Markov state transitions. These models provide an abstract description of the underlying conformational changes of the ion channels. However, due to the abstracted conformation states and assumptions for the rates of transition between them, there are differences between the models and reality-termed model discrepancy or misspecification. In this paper, we demonstrate the feasibility of using a mechanistically-inspired neural network differential equation model, a hybrid non-parametric model, to model ion channel kinetics. We apply it to the hERG potassium ion channel as an example, with the aim of providing an alternative modelling approach that could alleviate certain limitations of the traditional approach. We compare and discuss multiple ways of using a neural network to approximate extra hidden states or alternative transition rates. In particular we assess their ability to learn the missing dynamics, and ask whether we can use these models to handle model discrepancy. Finally, we discuss the practicality and limitations of using neural networks and their potential applications.

Understanding PITX2-Dependent Atrial Fibrillation Mechanisms through Computational Models

Atrial fibrillation (AF) is a common arrhythmia. Better prevention and treatment of AF are needed to reduce AF-associated morbidity and mortality. Several major mechanisms cause AF in patients, including genetic predispositions to AF development. Genome-wide association studies have identified a number of genetic variants in association with AF populations, with the strongest hits clustering on chromosome 4q25, close to the gene for the homeobox transcription PITX2. Because of the inherent complexity of the human heart, experimental and basic research is insufficient for understanding the functional impacts of PITX2 variants on AF. Linking PITX2 properties to ion channels, cells, tissues, atriums and the whole heart, computational models provide a supplementary tool for achieving a quantitative understanding of the functional role of PITX2 in remodelling atrial structure and function to predispose to AF. It is hoped that computational approaches incorporating all we know about PITX2-related structural and electrical remodelling would provide better understanding into its proarrhythmic effects leading to development of improved anti-AF therapies. In the present review, we discuss advances in atrial modelling and focus on the mechanistic links between PITX2 and AF. Challenges in applying models for improving patient health are described, as well as a summary of future perspectives.

A Comparison of Regional Classification Strategies Implemented for the Population Based Approach to Modelling Atrial Fibrillation

(1) Background: in silico models are increasingly relied upon to study the mechanisms of atrial fibrillation. Due to the complexity associated with atrial models, cellular variability is often ignored. Recent studies have shown that cellular variability may have a larger impact on electrophysiological behaviour than previously expected. This paper compares two methods for AF remodelling using regional populations. (2) Methods: using 200,000 action potentials, experimental data was used to calibrate healthy atrial regional populations with two cellular models. AF remodelling was applied by directly adjusting maximum channel conductances. AF remodelling was also applied through adjusting biomarkers. The methods were compared upon replication of experimental data. (3) Results: compared to the percentage method, the biomarker approach resulted in smaller changes. RMP, APD20, APD50, and APD90 were changed in the percentage method by up to 11%, 500%, 50%, and 60%, respectively. In the biomarker approach, RMP, APD20, APD50, and APD90 were changed by up to 4.5%, 132%, 50%, and 35%, respectively. (4) Conclusion: applying AF remodelling through biomarker-based clustering resulted in channel conductance changes that were consistent with experimental data, while maintaining the highly non-linear relationships between channel conductances and biomarkers. Directly changing conductances in the healthy regional populations impacted the non-linear relationships and resulted in non-physiological APD20 and APD50 values.