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

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.

Cell-to-cell heterogeneity in ion channel conductance impacts substrate vulnerability to arrhythmia

Determining whether an ectopic depolarization will lead to a self-perpetuating arrhythmia is of critical importance in determining arrhythmia risk, so it is necessary to understand what factors impact substrate vulnerability. This study sought to explore the impact of cell-to-cell heterogeneity in ion channel conductance on substrate vulnerability to arrhythmia by measuring the duration of the vulnerable window in computational models of one-dimensional cables of ventricular cardiomyocytes. We began by using a population of uniform cable models to determine the mechanisms underlying the vulnerable window phenomenon. We found that in addition to the known importance of GNa, the conductances GCa,L and GKr also play a minor role in determining the vulnerable window duration. We also found that a steeper slope of the repolarizing action potential during the vulnerable window correlated with a shorter vulnerable window duration in uniform cables. We applied our understanding from these initial simulations to an investigation of the vulnerable window in heterogeneous cable models. The heterogeneous cables displayed a great deal of intra-cable variation in vulnerable window duration, highly sensitive to the cardiomyocytes in the local environment of the ectopic stimulus. Coupling strength modulated not only the magnitude of the vulnerable window duration but also the extent of intra-tissue variability in vulnerable window duration.NEW & NOTEWORTHY We investigate the impact of cell-to-cell heterogeneity in ion channel conductance on substrate vulnerability to arrhythmia by measuring the vulnerable window duration in computational cardiomyocyte cable models. We demonstrate a wide range of intra-cable variability in vulnerable window duration (VWD) and show how this is changed by ion channel block and coupling strength perturbations.

In Silico TRials guide optimal stratification of ATrIal FIbrillation patients to Catheter Ablation and pharmacological medicaTION: the i-STRATIFICATION study

AIMS

Patients with persistent atrial fibrillation (AF) experience 50% recurrence despite pulmonary vein isolation (PVI), and no consensus is established for secondary treatments. The aim of our i-STRATIFICATION study is to provide evidence for stratifying patients with AF recurrence after PVI to optimal pharmacological and ablation therapies, through in silico trials.

METHODS AND RESULTS

A cohort of 800 virtual patients, with variability in atrial anatomy, electrophysiology, and tissue structure (low-voltage areas, LVAs), was developed and validated against clinical data from ionic currents to electrocardiogram. Virtual patients presenting AF post-PVI underwent 12 secondary treatments. Sustained AF developed in 522 virtual patients after PVI. Second ablation procedures involving left atrial ablation alone showed 55% efficacy, only succeeding in the small right atria (<60 mL). When additional cavo-tricuspid isthmus ablation was considered, Marshall-PLAN sufficed (66% efficacy) for the small left atria (<90 mL). For the bigger left atria, a more aggressive ablation approach was required, such as anterior mitral line (75% efficacy) or posterior wall isolation plus mitral isthmus ablation (77% efficacy). Virtual patients with LVAs greatly benefited from LVA ablation in the left and right atria (100% efficacy). Conversely, in the absence of LVAs, synergistic ablation and pharmacotherapy could terminate AF. In the absence of ablation, the patient's ionic current substrate modulated the response to antiarrhythmic drugs, being the inward currents critical for optimal stratification to amiodarone or vernakalant.

CONCLUSION

In silico trials identify optimal strategies for AF treatment based on virtual patient characteristics, evidencing the power of human modelling and simulation as a clinical assisting tool.

A Review of Personalised Cardiac Computational Modelling Using Electroanatomical Mapping Data

Computational models of cardiac electrophysiology have gradually matured during the past few decades and are now being personalised to provide patient-specific therapy guidance for improving suboptimal treatment outcomes. The predictive features of these personalised electrophysiology models hold the promise of providing optimal treatment planning, which is currently limited in the clinic owing to reliance on a population-based or average patient approach. The generation of a personalised electrophysiology model entails a sequence of steps for which a range of activation mapping, calibration methods and therapy simulation pipelines have been suggested. However, the optimal methods that can potentially constitute a clinically relevant in silico treatment are still being investigated and face limitations, such as uncertainty of electroanatomical data recordings, generation and calibration of models within clinical timelines and requirements to validate or benchmark the recovered tissue parameters. This paper is aimed at reporting techniques on the personalisation of cardiac computational models, with a focus on calibrating cardiac tissue conductivity based on electroanatomical mapping data.

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.

What determines the optimal pharmacological treatment of atrial fibrillation? Insights from in silico trials in 800 virtual atria

The best pharmacological treatment for each atrial fibrillation (AF) patient is unclear. We aim to exploit AF simulations in 800 virtual atria to identify key patient characteristics that guide the optimal selection of anti-arrhythmic drugs. The virtual cohort considered variability in electrophysiology and low voltage areas (LVA) and was developed and validated against experimental and clinical data from ionic currents to ECG. AF sustained in 494 (62%) atria, with large inward rectifier K+ current (IK1 ) and Na+ /K+ pump (INaK ) densities (IK1 0.11 ± 0.03 vs. 0.07 ± 0.03 S mF-1 ; INaK 0.68 ± 0.15 vs. 0.38 ± 26 S mF-1 ; sustained vs. un-sustained AF). In severely remodelled left atrium, with LVA extensions of more than 40% in the posterior wall, higher IK1 (median density 0.12 ± 0.02 S mF-1 ) was required for AF maintenance, and rotors localized in healthy right atrium. For lower LVA extensions, rotors could also anchor to LVA, in atria presenting short refractoriness (median L-type Ca2+ current, ICaL , density 0.08 ± 0.03 S mF-1 ). This atrial refractoriness, modulated by ICaL and fast Na+ current (INa ), determined pharmacological treatment success for both small and large LVA. Vernakalant was effective in atria presenting long refractoriness (median ICaL density 0.13 ± 0.05 S mF-1 ). For short refractoriness, atria with high INa (median density 8.92 ± 2.59 S mF-1 ) responded more favourably to amiodarone than flecainide, and the opposite was found in atria with low INa (median density 5.33 ± 1.41 S mF-1 ). In silico drug trials in 800 human atria identify inward currents as critical for optimal stratification of AF patient to pharmacological treatment and, together with the left atrial LVA extension, for accurately phenotyping AF dynamics. KEY POINTS: Atrial fibrillation (AF) maintenance is facilitated by small L-type Ca2+ current (ICaL ) and large inward rectifier K+ current (IK1 ) and Na+ /K+ pump. In severely remodelled left atrium, with low voltage areas (LVA) covering more than 40% of the posterior wall, sustained AF requires higher IK1 and rotors localize in healthy right atrium. For lower LVA extensions, rotors can also anchor to LVA, if the atria present short refractoriness (low ICaL ) Vernakalant is effective in atria presenting long refractoriness (high ICaL ). For short refractoriness, atria with fast Na+ current (INa ) up-regulation respond more favourably to amiodarone than flecainide, and the opposite is found in atria with low INa . The inward currents (ICaL and INa ) are critical for optimal stratification of AF patient to pharmacological treatment and, together with the left atrial LVA extension, for accurately phenotyping AF dynamics.