Anti-arrhythmic strategies for atrial fibrillation: The role of computational modeling in discovery, development, and optimization

Efficacy of catheter ablation for atrial fibrillation in heart failure: a meta-analysis of randomized controlled trials

The study aims to evaluate whether rhythm control by catheter ablation is superior to medical therapy for the patients with atrial fibrillation (AF) and heart failure (HF). The literatures were searched by using PubMed, Cochrane Library, Embase, and Web of Science databases up to 12 October 2023. The randomized controlled trials (RCTs) comparing rhythm control using catheter ablation vs. medical therapy in AF patients with HF were pooled. The primary outcomes included all-cause mortality, HF re-hospitalization, and stroke, and the secondary outcomes included left ventricular ejection fraction (LVEF), atrial tachyarrythmia recurrence, quality of life (Minnesota Living with Heart Failure Questionnaire score, MLHFQ score), 6 min walking distance (6MWD), the level of N-terminal B-type natriuretic peptide precursor (NT-proBNP), and adverse events. Nine RCTs involving in 2293 patients met the inclusion criteria. Compared with medical therapy, catheter ablation reduced all-cause mortality [10.07% (121/1201) vs. 15.26% (175/1147), risk ratio (RR):0.60, 95% confidence interval (CI): 0.48-0.74, P < 0.00001, I2 = 0%] and the rate of HF re-hospitalization (RR: 0.65, P = 0.02, 95% CI: 0.45 to 0.94, I2 = 74%), but had no obvious difference in incidence of stroke (RR: 0.67, P = 0.27, 95% CI: 0.32 to 1.38, I2 = 0%). Catheter ablation enhanced LVEF [mean difference (MD), 6.26%, P < 0.00001, I2 = 89%], reduced AT recurrence (RR: 0.37, P < 0.00001, 95% CI: 0.26 to 0.52, I2 = 89%), improved the quality of life (MLHFQ score) (MD: -6.83, P = 0.003, I2 = 67%), elevated 6MWD (MD: 15.92, P = 0.006, I2 = 76%), and diminished the level NT-proBNP (MD: -44.19, P < 0.00001, I2 = 75%), but had no significant difference in adverse events [25.81% (310/1201) vs. 30.25% (347/1147), RR: 0.81, 95% CI: 0.65-1.01, P = 0.06, I2 = 55%]. Catheter ablation as rhythm control strategy substantially enhances the survival rate, reduces HF re-hospitalization, increases the rate of sinus rhythm maintenance, improves the left ventricular function and the quality of life for AF patients with HF, and has similar safety, compared with medical therapy. The rhythm control by catheter ablation may be a better strategy for the AF patients with HF.

Multi-faceted potential of sophoridine compound's anti-arrhythmic and antioxidant effects through ROS/CaMKII pathway

Cardiac arrhythmias remain a significant cause of mortality and morbidity, for novel antiarrhythmic therapies. This study states that the first report of sophoridine (SPN), a quinolizidine alkaloid derived from traditional Chinese herbs, shows promise as a potential candidate due to its anti-arrhythmic and antioxidant properties. The study found that cell viability in H9C2 rat cardiomyocytes remained stable even when treated with SPN at a higher dosage of 100 μg/ml. This phenomenon was accompanied by increases in mitochondria-derived reactive oxygen species (ROS) and calcium/calmodulin-dependent protein kinase II (CaMKII) signaling, at 50 and 100 μg/ml. Glucose fluctuations regulate ventricular arrhythmias caused by SPN by activating the ROS/CaMKII pathway. Experimental models using zebrafish provided additional evidence supporting the regulatory effects of SPN on heart rate. In addition, the administration of SPN resulted in substantial deregulation of crucial genes involved in heart development (nppa, nppb, tnnt2a) at the transcriptional level in zebrafish. These findings provide insight into the various pharmacological properties of SPN and this opens up new possibilities for anti-arrhythmic treatment strategies.

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.

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.

Characterization of the Electrophysiological Characteristics of Chronic Atrial Fibrillation for Efficient Simulations

Atrial biophysical simulations have the potential to enhance outcomes by enabling the simulation of pharmacological and ablative strategies. However, the high computational times associated with such simulations render them unsuitable for diagnostic purposes. To address this challenge, discrete models such as cellular automata (CA) have been developed, which consider a finite number of states, thus significantly reducing computational times. Yet, there is a pressing need to determine whether CA can replicate pathological simulations with accuracy. The analysis of simulations under different degrees of electrical remodeling shows an expected increase of Action Potential Duration (APD) with the previous Diastolic Interval (DI) interval, indicating short-term memory of atrial cardiomyocytes: shorter APD0 provoked shorter APD+1, and previous DI has a similar effect on APD+1. Independent prediction using both APD0 and DI was found to provide a far better estimation of APD+1 values, compared to relying on DI alone (p<<0.01). Finally, the CA models were able to replicate reentrant patterns and cycle lengths of different states of atrial remodeling with a high degree of accuracy when compared to biophysical simulations. Overall, the use of atrial CA with short-term memory allows accurate reproduction of arrhythmic behavior in pathological tissue within a clinically relevant timeframe.Clinical Relevance- Discrete electrophysiological models simulate pathological self-sustained arrhythmias in diagnostic times.

Updates on epicardial adipose tissue mechanisms on atrial fibrillation

Obesity is a well-known risk factor for atrial fibrillation (AF). Local epi-myocardial or intra-myocardial adiposity caused by aging, obesity, or cardiovascular disease (CVD) is considered to be a better predictor of the risk of AF than general adiposity. Some of the described mechanisms suggest that epicardial adipose tissue (EAT) participates in structural remodeling owing to its endocrine activity or its infiltration between cardiomyocytes. Epicardial fat also wraps up the ganglionated plexi that reach the myocardium. Although the increment of volume/thickness and activity of EAT might modify autonomic activity, autonomic system dysfunction might also change the endocrine activity of epicardial fat in a feedback response. As a result, new preventive therapeutic strategies are focused on reducing adiposity and weight loss before AF ablation or inhibiting autonomic neurotransmitter secretion on fat pads during open-heart surgery to reduce the recurrence or postoperative risk of AF. In this manuscript, we review some of the novel findings regarding the pathophysiology and associated risk factors of AF, with special emphasis on the role of EAT in the electrical, structural, and molecular mechanisms of AF initiation and maintenance. In addition, we have included a brief note provided on epicardial fat preclinical models that could be useful for identifying new therapeutic targets.

Synergistic antiarrhythmic effect of inward rectifier current inhibition and pulmonary vein isolation in a 3D computer model for atrial fibrillation

AIMS

Recent clinical studies showed that antiarrhythmic drug (AAD) treatment and pulmonary vein isolation (PVI) synergistically reduce atrial fibrillation (AF) recurrences after initially successful ablation. Among newly developed atrial-selective AADs, inhibitors of the G-protein-gated acetylcholine-activated inward rectifier current (IKACh) were shown to effectively suppress AF in an experimental model but have not yet been evaluated clinically. We tested in silico whether inhibition of inward rectifier current or its combination with PVI reduces AF inducibility.

METHODS AND RESULTS

We simulated the effect of inward rectifier current blockade (IK blockade), PVI, and their combination on AF inducibility in a detailed three-dimensional model of the human atria with different degrees of fibrosis. IK blockade was simulated with a 30% reduction of its conductivity. Atrial fibrillation was initiated using incremental pacing applied at 20 different locations, in both atria. IK blockade effectively prevented AF induction in simulations without fibrosis as did PVI in simulations without fibrosis and with moderate fibrosis. Both interventions lost their efficacy in severe fibrosis. The combination of IK blockade and PVI prevented AF in simulations without fibrosis, with moderate fibrosis, and even with severe fibrosis. The combined therapy strongly decreased the number of fibrillation waves, due to a synergistic reduction of wavefront generation rate while the wavefront lifespan remained unchanged.

CONCLUSION

Newly developed blockers of atrial-specific inward rectifier currents, such as IKAch, might prevent AF occurrences and when combined with PVI effectively supress AF recurrences in human.

A Computational Study of the Effects of Tachycardia-Induced Remodeling on Calcium Wave Propagation in Rabbit Atrial Myocytes

In atrial cardiomyocytes without a well-developed T-tubule system, calcium diffuses from the periphery toward the center creating a centripetal wave pattern. During atrial fibrillation, rapid activation of atrial myocytes induces complex remodeling in diffusion properties that result in failure of calcium to propagate in a fully regenerative manner toward the center; a phenomenon termed “calcium silencing.” This has been observed in rabbit atrial myocytes after exposure to prolonged rapid pacing. Although experimental studies have pointed to possible mechanisms underlying calcium silencing, their individual effects and relative importance remain largely unknown. In this study we used computational modeling of the rabbit atrial cardiomyocyte to query the individual and combined effects of the proposed mechanisms leading to calcium silencing and abnormal calcium wave propagation. We employed a population of models obtained from a newly developed model of the rabbit atrial myocyte with spatial representation of intracellular calcium handling. We selected parameters in the model that represent experimentally observed cellular remodeling which have been implicated in calcium silencing, and scaled their values in the population to match experimental observations. In particular, we changed the maximum conductances of ICaL, INCX, and INaK, RyR open probability, RyR density, Serca2a density, and calcium buffering strength. We incorporated remodeling in a population of 16 models by independently varying parameters that reproduce experimentally observed cellular remodeling, and quantified the resulting alterations in calcium dynamics and wave propagation patterns. The results show a strong effect of ICaL in driving calcium silencing, with INCX, INaK, and RyR density also resulting in calcium silencing in some models. Calcium alternans was observed in some models where INCX and Serca2a density had been changed. Simultaneously incorporating changes in all remodeled parameters resulted in calcium silencing in all models, indicating the predominant role of decreasing ICaL in the population phenotype.

In Silico Assessment of Class I Antiarrhythmic Drug Effects on Pitx2-Induced Atrial Fibrillation: Insights from Populations of Electrophysiological Models of Human Atrial Cells and Tissues

Electrical remodelling as a result of homeodomain transcription factor 2 (Pitx2)-dependent gene regulation was linked to atrial fibrillation (AF) and AF patients with single nucleotide polymorphisms at chromosome 4q25 responded favorably to class I antiarrhythmic drugs (AADs). The possible reasons behind this remain elusive. The purpose of this study was to assess the efficacy of the AADs disopyramide, quinidine, and propafenone on human atrial arrhythmias mediated by Pitx2-induced remodelling, from a single cell to the tissue level, using drug binding models with multi-channel pharmacology. Experimentally calibrated populations of human atrial action po-tential (AP) models in both sinus rhythm (SR) and Pitx2-induced AF conditions were constructed by using two distinct models to represent morphological subtypes of AP. Multi-channel pharmaco-logical effects of disopyramide, quinidine, and propafenone on ionic currents were considered. Simulated results showed that Pitx2-induced remodelling increased maximum upstroke velocity (dVdtmax), and decreased AP duration (APD), conduction velocity (CV), and wavelength (WL). At the concentrations tested in this study, these AADs decreased dVdtmax and CV and prolonged APD in the setting of Pitx2-induced AF. Our findings of alterations in WL indicated that disopyramide may be more effective against Pitx2-induced AF than propafenone and quinidine by prolonging WL.

Cardiac transmembrane ion channels and action potentials: cellular physiology and arrhythmogenic behavior

Cardiac arrhythmias are among the leading causes of mortality. They often arise from alterations in the electrophysiological properties of cardiac cells and their underlying ionic mechanisms. It is therefore critical to further unravel the pathophysiology of the ionic basis of human cardiac electrophysiology in health and disease. In the first part of this review, current knowledge on the differences in ion channel expression and properties of the ionic processes that determine the morphology and properties of cardiac action potentials and calcium dynamics from cardiomyocytes in different regions of the heart are described. Then the cellular mechanisms promoting arrhythmias in congenital or acquired conditions of ion channel function (electrical remodeling) are discussed. The focus is on human-relevant findings obtained with clinical, experimental, and computational studies, given that interspecies differences make the extrapolation from animal experiments to human clinical settings difficult. Deepening the understanding of the diverse pathophysiology of human cellular electrophysiology will help in developing novel and effective antiarrhythmic strategies for specific subpopulations and disease conditions.

Populations of in silico myocytes and tissues reveal synergy of multiatrial-predominant K+ -current block in atrial fibrillation

BACKGROUND AND PURPOSE

Pharmacotherapy of atrial fibrillation (AF), the most common cardiac arrhythmia, remains unsatisfactory due to low efficacy and safety concerns. New therapeutic strategies target atrial-predominant ion-channels and involve multichannel block (poly)therapy. As AF is characterized by rapid and irregular atrial activations, compounds displaying potent antiarrhythmic effects at fast and minimal effects at slow rates are desirable. We present a novel systems pharmacology framework to quantitatively evaluate synergistic anti-AF effects of combined block of multiple atrial-predominant K⁺ currents (ultra-rapid delayed rectifier K⁺ current, I_Kur , small conductance Ca²⁺ -activated K⁺ current, I_KCa , K_2P 3.1 2-pore-domain K⁺ current, I_K2P ) in AF.

EXPERIMENTAL APPROACH

We constructed experimentally calibrated populations of virtual atrial myocyte models in normal sinus rhythm and AF-remodelled conditions using two distinct, well-established atrial models. Sensitivity analyses on our atrial populations was used to investigate the rate dependence of action potential duration (APD) changes due to blocking I_Kur , I_K2P or I_KCa and interactions caused by blocking of these currents in modulating APD. Block was simulated in both single myocytes and one-dimensional tissue strands to confirm insights from the sensitivity analyses and examine anti-arrhythmic effects of multi-atrial-predominant K⁺ current block in single cells and coupled tissue.

KEY RESULTS

In both virtual atrial myocytes and tissues, multiple atrial-predominant K⁺ -current block promoted favourable positive rate-dependent APD prolongation and displayed positive rate-dependent synergy, that is, increasing synergistic antiarrhythmic effects at fast pacing versus slow rates.

CONCLUSION AND IMPLICATIONS

Simultaneous block of multiple atrial-predominant K⁺ currents may be a valuable antiarrhythmic pharmacotherapeutic strategy for AF.

Understanding the Beat-to-Beat Variations of P-Waves Morphologies in AF Patients During Sinus Rhythm: A Scoping Review of the Atrial Simulation Studies

The remarkable advances in high-performance computing and the resulting increase of the computational power have the potential to leverage computational cardiology toward improving our understanding of the pathophysiological mechanisms of arrhythmias, such as Atrial Fibrillation (AF). In AF, a complex interaction between various triggers and the atrial substrate is considered to be the leading cause of AF initiation and perpetuation. In electrocardiography (ECG), P-wave is supposed to reflect atrial depolarization. It has been found that even during sinus rhythm (SR), multiple P-wave morphologies are present in AF patients with a history of AF, suggesting a higher dispersion of the conduction route in this population. In this scoping review, we focused on the mechanisms which modify the electrical substrate of the atria in AF patients, while investigating the existence of computational models that simulate the propagation of the electrical signal through different routes. The adopted review methodology is based on a structured analytical framework which includes the extraction of the keywords based on an initial limited bibliographic search, the extensive literature search and finally the identification of relevant articles based on the reference list of the studies. The leading mechanisms identified were classified according to their scale, spanning from mechanisms in the cell, tissue or organ level, and the produced outputs. The computational modeling approaches for each of the factors that influence the initiation and the perpetuation of AF are presented here to provide a clear overview of the existing literature. Several levels of categorization were adopted while the studies which aim to translate their findings to ECG phenotyping are highlighted. The results denote the availability of multiple models, which are appropriate under specific conditions. However, the consideration of complex scenarios taking into account multiple spatiotemporal scales, personalization of electrophysiological and anatomical models and the reproducibility in terms of ECG phenotyping has only partially been tackled so far.

Sensitivity of Ablation Targets Prediction to Electrophysiological Parameter Variability in Image-Based Computational Models of Ventricular Tachycardia in Post-infarction Patients

Ventricular tachycardia (VT), which could lead to sudden cardiac death, occurs frequently in patients with myocardial infarction. Computational modeling has emerged as a powerful platform for the non-invasive investigation of lethal heart rhythm disorders in post-infarction patients and for guiding patient VT ablation. However, it remains unclear how VT dynamics and predicted ablation targets are influenced by inter-patient variability in action potential duration (APD) and conduction velocity (CV). The goal of this study was to systematically assess the effect of changes in the electrophysiological parameters on the induced VTs and predicted ablation targets in personalized models of post-infarction hearts. Simulations were conducted in 5 patient-specific left ventricular models reconstructed from late gadolinium-enhanced magnetic resonance imaging scans. We comprehensively characterized all possible pre-ablation and post-ablation VTs in simulations conducted with either an “average human VT”-based electrophysiological representation (i.e., EP_avg) or with ±10% APD or CV (i.e., EP_var); additional simulations were also executed in some models for an extended range of these parameters. The results showed that: (1) a subset of reentries (76.2–100%, depending on EP parameter set) conducted with ±10% APD/CV was observed in approximately the same locations as reentries observed in EP_avg cases; (2) emergent VTs could be induced sometimes after ablation in EP_avg models, and these emergent VTs often corresponded to the pre-ablation reentries in simulations with EP_var parameter sets. These findings demonstrate that the VT ablation target uncertainty in patient-specific ventricular models with an average representation of VT-remodeled electrophysiology is relatively low and the ablation targets stable, as the localization of the induced VTs was primarily driven by the remodeled structural substrate. Thus, personalized ventricular modeling with an average representation of infarct-remodeled electrophysiology may uncover most targets for VT ablation.

The Association Between Diabetes Mellitus and Atrial Fibrillation: Clinical and Mechanistic Insights

A number of clinical studies have reported that diabetes mellitus (DM) is an independent risk factor for Atrial fibrillation (AF). After adjustment for other known risk factors including age, sex, and cardiovascular risk factors, DM remains a significant if modest risk factor for development of AF. The mechanisms underlying the increased susceptibility to AF in DM are incompletely understood, but are thought to involve electrical, structural, and autonomic remodeling in the atria. Electrical remodeling in DM may involve alterations in gap junction function that affect atrial conduction velocity due to changes in expression or localization of connexins. Electrical remodeling can also occur due to changes in atrial action potential morphology in association with changes in ionic currents, such as sodium or potassium currents, that can affect conduction velocity or susceptibility to triggered activity. Structural remodeling in DM results in atrial fibrosis, which can alter conduction patterns and susceptibility to re-entry in the atria. In addition, increases in atrial adipose tissue, especially in Type II DM, can lead to disruptions in atrial conduction velocity or conduction patterns that may affect arrhythmogenesis. Whether the insulin resistance in type II DM activates unique intracellular signaling pathways independent of obesity requires further investigation. In addition, the relationship between incident AF and glycemic control requires further study.

Computational modeling: What does it tell us about atrial fibrillation therapy?

Atrial fibrillation (AF) is a complex cardiac arrhythmia with diverse etiology that negatively affects morbidity and mortality of millions of patients. Technological and experimental advances have provided a wealth of information on the pathogenesis of AF, highlighting a multitude of mechanisms involved in arrhythmia initiation and maintenance, and disease progression. However, it remains challenging to identify the predominant mechanisms for specific subgroups of AF patients, which, together with an incomplete understanding of the pleiotropic effects of antiarrhythmic therapies, likely contributes to the suboptimal efficacy of current antiarrhythmic approaches. Computer modeling of cardiac electrophysiology has advanced in parallel to experimental research and provides an integrative framework to attempt to overcome some of these challenges. Multi-scale cardiac modeling and simulation integrate structural and functional data from experimental and clinical work with knowledge of atrial electrophysiological mechanisms and dynamics, thereby improving our understanding of AF mechanisms and therapy. In this review, we describe recent advances in our quantitative understanding of AF through mathematical models. We discuss computational modeling of AF mechanisms and therapy using detailed, mechanistic cell/tissue-level models, including approaches to incorporate variability in patient populations. We also highlight efforts using whole-atria models to improve catheter ablation therapies. Finally, we describe recent efforts and suggest future extensions to model clinical concepts of AF using patient-level models.

A Heart for Diversity: Simulating Variability in Cardiac Arrhythmia Research

In cardiac electrophysiology, there exist many sources of inter- and intra-personal variability. These include variability in conditions and environment, and genotypic and molecular diversity, including differences in expression and behavior of ion channels and transporters, which lead to phenotypic diversity (e.g., variable integrated responses at the cell, tissue, and organ levels). These variabilities play an important role in progression of heart disease and arrhythmia syndromes and outcomes of therapeutic interventions. Yet, the traditional in silico framework for investigating cardiac arrhythmias is built upon a parameter/property-averaging approach that typically overlooks the physiological diversity. Inspired by work done in genetics and neuroscience, new modeling frameworks of cardiac electrophysiology have been recently developed that take advantage of modern computational capabilities and approaches, and account for the variance in the biological data they are intended to illuminate. In this review, we outline the recent advances in statistical and computational techniques that take into account physiological variability, and move beyond the traditional cardiac model-building scheme that involves averaging over samples from many individuals in the construction of a highly tuned composite model. We discuss how these advanced methods have harnessed the power of big (simulated) data to study the mechanisms of cardiac arrhythmias, with a special emphasis on atrial fibrillation, and improve the assessment of proarrhythmic risk and drug response. The challenges of using in silico approaches with variability are also addressed and future directions are proposed.

In Silico Assessment of Efficacy and Safety of IKur Inhibitors in Chronic Atrial Fibrillation: Role of Kinetics and State-Dependence of Drug Binding

Current pharmacological therapy against atrial fibrillation (AF), the most common cardiac arrhythmia, is limited by moderate efficacy and adverse side effects including ventricular proarrhythmia and organ toxicity. One way to circumvent the former is to target ion channels that are predominantly expressed in atria vs. ventricles, such as K_V1.5, carrying the ultra-rapid delayed-rectifier K⁺ current (I_Kur). Recently, we used an in silico strategy to define optimal K_V1.5-targeting drug characteristics, including kinetics and state-dependent binding, that maximize AF-selectivity in human atrial cardiomyocytes in normal sinus rhythm (nSR). However, because of evidence for I_Kur being strongly diminished in long-standing persistent (chronic) AF (cAF), the therapeutic potential of drugs targeting I_Kur may be limited in cAF patients. Here, we sought to simulate the efficacy (and safety) of I_Kur inhibitors in cAF conditions. To this end, we utilized sensitivity analysis of our human atrial cardiomyocyte model to assess the importance of I_Kur for atrial cardiomyocyte electrophysiological properties, simulated hundreds of theoretical drugs to reveal those exhibiting anti-AF selectivity, and compared the results obtained in cAF with those in nSR. We found that despite being downregulated, I_Kur contributes more prominently to action potential (AP) and effective refractory period (ERP) duration in cAF vs. nSR, with ideal drugs improving atrial electrophysiology (e.g., ERP prolongation) more in cAF than in nSR. Notably, the trajectory of the AP during cAF is such that more I_Kur is available during the more depolarized plateau potential. Furthermore, I_Kur block in cAF has less cardiotoxic effects (e.g., AP duration not exceeding nSR values) and can increase Ca²⁺ transient amplitude thereby enhancing atrial contractility. We propose that in silico strategies such as that presented here should be combined with in vitro and in vivo assays to validate model predictions and facilitate the ongoing search for novel agents against AF.

Multiple mechanisms mediating carbon monoxide inhibition of the voltage-gated K+ channel Kv1.5

The voltage-gated K+ channel has key roles in the vasculature and in atrial excitability and contributes to apoptosis in various tissues. In this study, we have explored its regulation by carbon monoxide (CO), a product of the cytoprotective heme oxygenase enzymes, and a recognized toxin. CO inhibited recombinant Kv1.5 expressed in HEK293 cells in a concentration-dependent manner that involved multiple signalling pathways. CO inhibition was partially reversed by superoxide dismutase mimetics and by suppression of mitochondrial reactive oxygen species. CO also elevated intracellular nitric oxide (NO) levels. Prevention of NO formation also partially reversed CO inhibition of Kv1.5, as did inhibition of soluble guanylyl cyclase. CO also elevated intracellular peroxynitrite levels, and a peroxynitrite scavenger markedly attenuated the ability of CO to inhibit Kv1.5. CO caused nitrosylation of Kv1.5, an effect that was also observed in C331A and C346A mutant forms of the channel, which had previously been suggested as nitrosylation sites within Kv1.5. Augmentation of Kv1.5 via exposure to hydrogen peroxide was fully reversed by CO. Native Kv1.5 recorded in HL-1 murine atrial cells was also inhibited by CO. Action potentials recorded in HL-1 cells were increased in amplitude and duration by CO, an effect mimicked and occluded by pharmacological inhibition of Kv1.5. Our data indicate that Kv1.5 is a target for modulation by CO via multiple mechanisms. This regulation has important implications for diverse cellular functions, including excitability, contractility and apoptosis.

Mechanistic Systems Modeling to Improve Understanding and Prediction of Cardiotoxicity Caused by Targeted Cancer Therapeutics

Tyrosine kinase inhibitors (TKIs) are highly potent cancer therapeutics that have been linked with serious cardiotoxicity, including left ventricular dysfunction, heart failure, and QT prolongation. TKI-induced cardiotoxicity is thought to result from interference with tyrosine kinase activity in cardiomyocytes, where these signaling pathways help to control critical processes such as survival signaling, energy homeostasis, and excitation–contraction coupling. However, mechanistic understanding is limited at present due to the complexities of tyrosine kinase signaling, and the wide range of targets inhibited by TKIs. Here, we review the use of TKIs in cancer and the cardiotoxicities that have been reported, discuss potential mechanisms underlying cardiotoxicity, and describe recent progress in achieving a more systematic understanding of cardiotoxicity via the use of mechanistic models. In particular, we argue that future advances are likely to be enabled by studies that combine large-scale experimental measurements with Quantitative Systems Pharmacology (QSP) models describing biological mechanisms and dynamics. As such approaches have proven extremely valuable for understanding and predicting other drug toxicities, it is likely that QSP modeling can be successfully applied to cardiotoxicity induced by TKIs. We conclude by discussing a potential strategy for integrating genome-wide expression measurements with models, illustrate initial advances in applying this approach to cardiotoxicity, and describe challenges that must be overcome to truly develop a mechanistic and systematic understanding of cardiotoxicity caused by TKIs.

A computational framework for testing arrhythmia marker sensitivities to model parameters in functionally calibrated populations of atrial cells

Models of cardiac cell electrophysiology are complex non-linear systems which can be used to gain insight into mechanisms of cardiac dynamics in both healthy and pathological conditions. However, the complexity of cardiac models can make mechanistic insight difficult. Moreover, these are typically fitted to averaged experimental data which do not incorporate the variability in observations. Recently, building populations of models to incorporate inter- and intra-subject variability in simulations has been combined with sensitivity analysis (SA) to uncover novel ionic mechanisms and potentially clarify arrhythmogenic behaviors. We used the Koivumäki human atrial cell model to create two populations, representing normal Sinus Rhythm (nSR) and chronic Atrial Fibrillation (cAF), by varying 22 key model parameters. In each population, 14 biomarkers related to the action potential and dynamic restitution were extracted. Populations were calibrated based on distributions of biomarkers to obtain reasonable physiological behavior, and subjected to SA to quantify correlations between model parameters and pro-arrhythmia markers. The two populations showed distinct behaviors under steady state and dynamic pacing. The nSR population revealed greater variability, and more unstable dynamic restitution, as compared to the cAF population, suggesting that simulated cAF remodeling rendered cells more stable to parameter variation and rate adaptation. SA revealed that the biomarkers depended mainly on five ionic currents, with noted differences in sensitivities to these between nSR and cAF. Also, parameters could be selected to produce a model variant with no alternans and unaltered action potential morphology, highlighting that unstable dynamical behavior may be driven by specific cell parameter settings. These results ultimately suggest that arrhythmia maintenance in cAF may not be due to instability in cell membrane excitability, but rather due to tissue-level effects which promote initiation and maintenance of reentrant arrhythmia.

Revealing kinetics and state-dependent binding properties of IKur-targeting drugs that maximize atrial fibrillation selectivity

The K_V1.5 potassium channel, which underlies the ultra-rapid delayed-rectifier current (I_Kur) and is predominantly expressed in atria vs. ventricles, has emerged as a promising target to treat atrial fibrillation (AF). However, while numerous K_V1.5-selective compounds have been screened, characterized, and tested in various animal models of AF, evidence of antiarrhythmic efficacy in humans is still lacking. Moreover, current guidelines for pre-clinical assessment of candidate drugs heavily rely on steady-state concentration-response curves or IC₅₀ values, which can overlook adverse cardiotoxic effects. We sought to investigate the effects of kinetics and state-dependent binding of I_Kur-targeting drugs on atrial electrophysiology in silico and reveal the ideal properties of I_Kur blockers that maximize anti-AF efficacy and minimize pro-arrhythmic risk. To this aim, we developed a new Markov model of I_Kur that describes K_V1.5 gating based on experimental voltage-clamp data in atrial myocytes from patient right-atrial samples in normal sinus rhythm. We extended the I_Kur formulation to account for state-specificity and kinetics of K_V1.5-drug interactions and incorporated it into our human atrial cell model. We simulated 1- and 3-Hz pacing protocols in drug-free conditions and with a [drug] equal to the IC₅₀ value. The effects of binding and unbinding kinetics were determined by examining permutations of the forward (k_on) and reverse (k_off) binding rates to the closed, open, and inactivated states of the K_V1.5 channel. We identified a subset of ideal drugs exhibiting anti-AF electrophysiological parameter changes at fast pacing rates (effective refractory period prolongation), while having little effect on normal sinus rhythm (limited action potential prolongation). Our results highlight that accurately accounting for channel interactions with drugs, including kinetics and state-dependent binding, is critical for developing safer and more effective pharmacological anti-AF options.

Connexin 43 is involved in the sympathetic atrial fibrillation in canine and canine atrial myocytes

Objective:

Atrial fibrillation (AF) is the most common rapid cardiac arrhythmia associated with high morbidity and mortality. Stimulation of the sympathetic nerve is involved in AF occurrence. The gap junction protein connexin 43 (Cx43) plays a key role in electrical conduction velocity in cardiac tissues, and under expression of Cx43 was linked with AF. The aim of this study was to investigate whether Cx43 was involved in sympathetic AF.

Methods:

Fifteen dogs were randomly divided into 3 groups (5 in each group). Sympathetic AF was induced in dogs and isolated canine atrial myocytes by isoproterenol (ISO) perfusion and rapid atrium pacing (RAP). The expression levels of nerve growth factor (NGF) and tyrosine hydroxylase (TH) in the atrial tissues were detected using immunohistochemical staining. The transcription and protein expression of Cx43 in the AF cell model was measured. Subsequently, Cx43 was blocked by short interfering (si) RNA in atrial myocytes and the gap junctional inter-cellular communication was detected using the scrape-loading and dye transfer assay.

Results:

Sympathetic AF was successfully induced by a combination of ISO perfusion and RAP. The expression levels of NGF and TH were increased in the RAP group, and further increased in the RAP + ISO group. Tissue samples from the AF dogs had a lower Cx43 level than those of the control group (p<0.05). The expressions of mRNA and protein of Cx43 in sympathetic AF cell model decreased by 26% and 28%, respectively, when compared with the control group, with p<0.05. Silencing Cx43 in cells by siRNA could also efficiently reduce Cx43 expression. The relative levels of Cx45 mRNA were decreased by 73% compared with unaffected cells. The scrape-loading and dye transfer assay showed that gap junctional intercellular communication was hampered in the sympathetic AF cell model and silencing Cx43 could impede channel conduction.

Conclusion:

The results suggested that low expression of Cx43 was involved in sympathetic AF by influencing intercellular channel conduction. Intervention of Cx43 expression might be an appealing therapy to sympathetic AF.

Atrial fibrillation septal pacing: translation of modelling results

AIMS

Atrial fibrillation (AF) septal pacing consists of rapid pacing from a ring of electrodes around the atrial septum, leading to local capture of both atria during AF. The present model-based study evaluated the impact of the number of stimulation electrodes in the septal ring on AF capture for different types of sustained AF dynamics.

METHODS AND RESULTS

Using a biophysical model of AF based on CT scans from an AF patient, models with different AF substrates (Cholinergic AF and Meandering Wavelets) were created by varying the atrial membrane kinetics. Rapid pacing was applied from the septum area with a ring of 1, 2, 3, 4, 6, 8, or 12 electrodes during 20 seconds at a pacing cycle lengths (PCLs) in the range 60-100% of AF cycle length (AFCL), in 4% steps. Percentage of captured tissue during rapid pacing was determined using 24 sensing electrode pairs evenly distributed on the atrial surface. Results were averaged over 10 AF simulations. For Cholinergic AF, the number of stimulation electrodes on the septal ring had no significant impact on AF capture independently of AF dynamics. For Meandering Wavelets, more electrodes were needed to achieve AF capture in the presence of complex AF.

CONCLUSION

Changes in AF substrate significantly impacted septal pacing outcomes and response to rapid AF pacing may similarly vary patient-to-patient. The number of stimulation electrodes had a lesser impact, suggesting that the design of a ring with 3-4 electrodes around the septum would be sufficient for most AF dynamics.