ABSTRACT
BACKGROUND AND AIMS: The effective refractory period (ERP) is one of the main electrophysiological properties governing arrhythmia, yet ERP personalisation is rarely performed when creating patient-specific computer models of the atria to inform clinical decision-making. This study evaluates the impact of integrating clinical ERP measurements into personalised in silico models on arrhythmia vulnerability. METHODS: Clinical ERP measurements were obtained in seven patients from multiple locations in the atria. Atrial geometries from the electroanatomical mapping system were used to generate personalised anatomical atrial models. The Courtemanche cellular model was adjusted to reproduce patientspecific ERP. Four modelling approaches were compared: homogeneous (A), heterogeneous (B), regional (C), and continuous (D) ERP distributions. Non-personalised approaches (A, B) were based on literature data, while personalised approaches (C, D) were based on patient measurements. Modelling effects were assessed on arrhythmia vulnerability and tachycardia cycle length, with sensitivity analysis on ERP measurement uncertainty. RESULTS: Mean vulnerability was 3.4±4.0%, 7.7±3.4%, 9.0±5.1%, 7.0±3.6% for scenarios A to D, respectively. Mean tachycardia cycle length was 167.1±12.6 ms, 158.4±27.5 ms, 265.2±39.9 ms, and 285.9±77.3 ms for scenarios A to D, respectively. Incorporating perturbations to the measured ERP in the range of 2, 5, 10, 20, and 50ms changed the vulnerability of the model to 5.8±2.7%, 6.1±3.5%, 6.9±3.7%, 5.2±3.5%, 9.7±10.0% respectively. CONCLUSION: Increased ERP dispersion had a greater effect on reentry dynamics than on vulnerability. Inducibility was higher in personalised scenarios compared to scenarios with uniformly reduced ERP; however, this effect was reversed when incorporating fibrosis informed by low voltage areas.ERP measurement uncertainty up to 20 ms slightly influenced vulnerability. Electrophysiological personalisation of atrial in silico models appears essential and requires confirmation in larger cohorts.
ABSTRACT
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.
Subject(s)
Anti-Arrhythmia Agents , Atrial Fibrillation , Catheter Ablation , Pulmonary Veins , Recurrence , Atrial Fibrillation/surgery , Atrial Fibrillation/drug therapy , Atrial Fibrillation/physiopathology , Atrial Fibrillation/diagnosis , Humans , Catheter Ablation/methods , Pulmonary Veins/surgery , Pulmonary Veins/physiopathology , Anti-Arrhythmia Agents/therapeutic use , Treatment Outcome , Models, Cardiovascular , Computer Simulation , Action Potentials , Risk Assessment , Heart Atria/physiopathology , Heart Atria/surgery , Male , Anisoles/therapeutic use , Patient Selection , Female , Patient-Specific Modeling , Middle Aged , Pyrrolidines/therapeutic use , Electrocardiography , Clinical Decision-MakingABSTRACT
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.
Subject(s)
Amiodarone , Atrial Fibrillation , Humans , Atrial Fibrillation/drug therapy , Flecainide/pharmacology , Flecainide/therapeutic use , Heart Atria , Amiodarone/pharmacology , Amiodarone/therapeutic use , Action Potentials/physiologyABSTRACT
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.
ABSTRACT
The high incidence of cardiac arrythmias underlines the need for the assessment of pharmacological therapies. In this field of drug efficacy, as in the field of drug safety highlighted by the Comprehensive in Vitro Proarrhythmia Assay initiative, new pillars for research have become crucial: firstly, the integration of in-silico experiments, and secondly the evaluation of fully integrated biological systems, such as human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs). In this study, we therefore aimed to combine in-vitro experiments and in-silico simulations to evaluate the antiarrhythmic effect of L-type calcium current (ICaL) block in hiPSC-CMs. For this, hiPSC-CM preparations were cultured and an equivalent virtual tissue was modeled. Re-entry patterns of electrical activation were induced and several biomarkers were obtained before and after ICaL block. The virtual hiPSC-CM simulations were also reproduced using a tissue composed of adult ventricular cardiomyocytes (hAdultV-CMs). The analysis of phases, currents and safety factor for propagation showed an increased size of the re-entry core when ICaL was blocked as a result of depressed cellular excitability. The bigger wavefront curvature yielded reductions of 12.2%, 6.9%, and 4.2% in the frequency of the re-entry for hiPSC-CM cultures, virtual hiPSC-CM, and hAdultV-CM tissues, respectively. Furthermore, ICaL block led to a 47.8% shortening of the vulnerable window for re-entry in the virtual hiPSC-CM tissue and to re-entry vanishment in hAdultV-CM tissue. The consistent behavior between in-vitro and in-silico hiPSC-CMs and between in-silico hiPSC-CMs and hAdultV-CMs evidences that virtual hiPSC-CM tissues are suitable for assessing cardiac efficacy, as done in the present study through the analysis of ICaL block.