1
|
Salvador M, Strocchi M, Regazzoni F, Augustin CM, Dede' L, Niederer SA, Quarteroni A. Whole-heart electromechanical simulations using Latent Neural Ordinary Differential Equations. NPJ Digit Med 2024; 7:90. [PMID: 38605089 PMCID: PMC11009296 DOI: 10.1038/s41746-024-01084-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2023] [Accepted: 03/22/2024] [Indexed: 04/13/2024] Open
Abstract
Cardiac digital twins provide a physics and physiology informed framework to deliver personalized medicine. However, high-fidelity multi-scale cardiac models remain a barrier to adoption due to their extensive computational costs. Artificial Intelligence-based methods can make the creation of fast and accurate whole-heart digital twins feasible. We use Latent Neural Ordinary Differential Equations (LNODEs) to learn the pressure-volume dynamics of a heart failure patient. Our surrogate model is trained from 400 simulations while accounting for 43 parameters describing cell-to-organ cardiac electromechanics and cardiovascular hemodynamics. LNODEs provide a compact representation of the 3D-0D model in a latent space by means of an Artificial Neural Network that retains only 3 hidden layers with 13 neurons per layer and allows for numerical simulations of cardiac function on a single processor. We employ LNODEs to perform global sensitivity analysis and parameter estimation with uncertainty quantification in 3 hours of computations, still on a single processor.
Collapse
Affiliation(s)
- Matteo Salvador
- Institute for Computational and Mathematical Engineering, Stanford University, California, CA, USA.
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK.
- MOX, Department of Mathematics, Politecnico di Milano, Milan, Italy.
| | - Marina Strocchi
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
- National Heart and Lung Institute, Imperial College London, London, UK
| | | | - Christoph M Augustin
- Institute of Biophysics, Medical University of Graz, Graz, Austria
- BioTechMed-Graz, Graz, Austria
| | - Luca Dede'
- MOX, Department of Mathematics, Politecnico di Milano, Milan, Italy
| | - Steven A Niederer
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
- National Heart and Lung Institute, Imperial College London, London, UK
- The Alan Turing Institute, London, UK
| | - Alfio Quarteroni
- MOX, Department of Mathematics, Politecnico di Milano, Milan, Italy
- École Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
| |
Collapse
|
2
|
Ernault AC, Al-Shama RFM, Li J, Devalla HD, de Groot JR, Coronel R, Vigmond E, Boukens BJ. Interpretation of field and LEAP potentials recorded from cardiomyocyte monolayers. Am J Physiol Heart Circ Physiol 2024; 326:H800-H811. [PMID: 38180452 DOI: 10.1152/ajpheart.00463.2023] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 07/27/2023] [Revised: 12/04/2023] [Accepted: 01/02/2024] [Indexed: 01/06/2024]
Abstract
Multielectrode arrays (MEAs) are the method of choice for electrophysiological characterization of cardiomyocyte monolayers. The field potentials recorded using an MEA are like extracellular electrograms recorded from the myocardium using conventional electrodes. Nevertheless, different criteria are used to interpret field potentials and extracellular electrograms, which hamper correct interpretation and translation to the patient. To validate the criteria for interpretation of field potentials, we used neonatal rat cardiomyocytes to generate monolayers. We recorded field potentials using an MEA and simultaneously recorded action potentials using sharp microelectrodes. In parallel, we recreated our experimental setting in silico and performed simulations. We show that the amplitude of the local RS complex of a field potential correlated with conduction velocity in silico but not in vitro. The peak time of the T wave in field potentials exhibited a strong correlation with APD90 while the steepest upslope correlated well with APD50. However, this relationship only holds when the T wave displayed a biphasic pattern. Next, we simulated local extracellular action potentials (LEAPs). The shape of the LEAP differed markedly from the shape of the local action potential, but the final duration of the LEAP coincided with APD90. Criteria for interpretation of extracellular electrograms should be applied to field potentials. This will provide a strong basis for the analysis of heterogeneity in conduction velocity and repolarization in cultured monolayers of cardiomyocytes. Finally, a LEAP is not a recording of the local action potential but is generated by intracellular current provided by neighboring cardiomyocytes and is superior to field potential duration in estimating APD90.NEW & NOTEWORTHY We present a physiological basis for the interpretation of multielectrode array-derived, extracellular, electrical signals.
Collapse
Affiliation(s)
- Auriane C Ernault
- Heart Center, Department of Clinical and Experimental Cardiology, Amsterdam Cardiovascular Sciences, Amsterdam University Medical Centers, Amsterdam, The Netherlands
| | - Rushd F M Al-Shama
- Heart Center, Department of Clinical and Experimental Cardiology, Amsterdam Cardiovascular Sciences, Amsterdam University Medical Centers, Amsterdam, The Netherlands
| | - Jiuru Li
- Department of Medical Biology, Amsterdam Cardiovascular Sciences, Amsterdam University Medical Centers, Amsterdam, The Netherlands
| | - Harsha D Devalla
- Department of Medical Biology, Amsterdam Cardiovascular Sciences, Amsterdam University Medical Centers, Amsterdam, The Netherlands
| | - Joris R de Groot
- Heart Center, Department of Clinical and Experimental Cardiology, Amsterdam Cardiovascular Sciences, Amsterdam University Medical Centers, Amsterdam, The Netherlands
| | - Ruben Coronel
- Heart Center, Department of Clinical and Experimental Cardiology, Amsterdam Cardiovascular Sciences, Amsterdam University Medical Centers, Amsterdam, The Netherlands
| | - Edward Vigmond
- IHU Liryc, Fondation Bordeaux Université, Bordeaux, France
- University of Bordeaux, Talence, France
| | - Bastiaan J Boukens
- Department of Medical Biology, Amsterdam Cardiovascular Sciences, Amsterdam University Medical Centers, Amsterdam, The Netherlands
- Department of Physiology, Cardiovascular Research Institute Maastricht, Maastricht University Medical Center, Maastricht, The Netherlands
| |
Collapse
|
3
|
Qian S, Ugurlu D, Fairweather E, Strocchi M, Toso LD, Deng Y, Plank G, Vigmond E, Razavi R, Young A, Lamata P, Bishop M, Niederer S. Developing Cardiac Digital Twins at Scale: Insights from Personalised Myocardial Conduction Velocity. MEDRXIV : THE PREPRINT SERVER FOR HEALTH SCIENCES 2024:2023.12.05.23299435. [PMID: 38106072 PMCID: PMC10723499 DOI: 10.1101/2023.12.05.23299435] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/19/2023]
Abstract
Large-cohort studies using cardiovascular imaging and diagnostic datasets have assessed cardiac anatomy, function, and outcomes, but typically do not reveal underlying biological mechanisms. Cardiac digital twins (CDTs) provide personalized physics- and physiology-constrained in-silico representations, enabling inference of multi-scale properties tied to these mechanisms. We constructed 3464 anatomically-accurate CDTs using cardiac magnetic resonance images from UK biobank and personalised their myocardial conduction velocities (CVs) from electrocardiograms (ECG), through an automated framework. We found well-known sex-specific differences in QRS duration were fully explained by myocardial anatomy, as CV remained consistent across sexes. Conversely, significant associations of CV with ageing and increased BMI suggest myocardial tissue remodelling. Novel associations were observed with left ventricular ejection fraction and mental-health phenotypes, through a phenome-wide association study, and CV was also linked with adverse clinical outcomes. Our study highlights the utility of population-based CDTs in assessing intersubject variability and uncovering strong links with mental health.
Collapse
|
4
|
Ochs AR, Boyle PM. Optogenetic Modulation of Arrhythmia Triggers: Proof-of-Concept from Computational Modeling. Cell Mol Bioeng 2023; 16:243-259. [PMID: 37810996 PMCID: PMC10550900 DOI: 10.1007/s12195-023-00781-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2023] [Accepted: 08/14/2023] [Indexed: 10/10/2023] Open
Abstract
Introduction Early afterdepolarizations (EADs) are secondary voltage depolarizations associated with reduced repolarization reserve (RRR) that can trigger lethal arrhythmias. Relating EADs to triggered activity is difficult to study, so the ability to suppress or provoke EADs would be experimentally useful. Here, we use computational simulations to assess the feasibility of subthreshold optogenetic stimulation modulating the propensity for EADs (cell-scale) and EAD-associated ectopic beats (organ-scale). Methods We modified a ventricular ionic model by reducing rapid delayed rectifier potassium (0.25-0.1 × baseline) and increasing L-type calcium (1.0-3.5 × baseline) currents to create RRR conditions with varying severity. We ran simulations in models of single cardiomyocytes and left ventricles from post-myocardial infarction patient MRI scans. Optogenetic stimulation was simulated using either ChR2 (depolarizing) or GtACR1 (repolarizing) opsins. Results In cell-scale simulations without illumination, EADs were seen for 164 of 416 RRR conditions. Subthreshold stimulation of GtACR1 reduced EAD incidence by up to 84.8% (25/416 RRR conditions; 0.1 μW/mm2); in contrast, subthreshold ChR2 excitation increased EAD incidence by up to 136.6% (388/416 RRR conditions; 50 μW/mm2). At the organ scale, we assumed simultaneous, uniform illumination of the epicardial and endocardial surfaces. GtACR1-mediated suppression (10-50 μW/mm2) and ChR2-mediated unmasking (50-100 μW/mm2) of EAD-associated ectopic beats were feasible in three distinct ventricular models. Conclusions Our findings suggest that optogenetics could be used to silence or provoke both EADs and EAD-associated ectopic beats. Validation in animal models could lead to exciting new experimental regimes and potentially to novel anti-arrhythmia treatments. Supplementary Information The online version contains supplementary material available at 10.1007/s12195-023-00781-z.
Collapse
Affiliation(s)
- Alexander R. Ochs
- Department of Bioengineering, UW Bioengineering, University of Washington, 3720 15th Ave NE N107, UW Mailbox 355061, Seattle, WA 98195 USA
| | - Patrick M. Boyle
- Department of Bioengineering, UW Bioengineering, University of Washington, 3720 15th Ave NE N107, UW Mailbox 355061, Seattle, WA 98195 USA
- Institute for Stem Cell and Regenerative Medicine, University of Washington, Seattle, WA USA
- Center for Cardiovascular Biology, University of Washington, Seattle, WA USA
| |
Collapse
|
5
|
Strocchi M, Longobardi S, Augustin CM, Gsell MAF, Petras A, Rinaldi CA, Vigmond EJ, Plank G, Oates CJ, Wilkinson RD, Niederer SA. Cell to whole organ global sensitivity analysis on a four-chamber heart electromechanics model using Gaussian processes emulators. PLoS Comput Biol 2023; 19:e1011257. [PMID: 37363928 DOI: 10.1371/journal.pcbi.1011257] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2022] [Accepted: 06/09/2023] [Indexed: 06/28/2023] Open
Abstract
Cardiac pump function arises from a series of highly orchestrated events across multiple scales. Computational electromechanics can encode these events in physics-constrained models. However, the large number of parameters in these models has made the systematic study of the link between cellular, tissue, and organ scale parameters to whole heart physiology challenging. A patient-specific anatomical heart model, or digital twin, was created. Cellular ionic dynamics and contraction were simulated with the Courtemanche-Land and the ToR-ORd-Land models for the atria and the ventricles, respectively. Whole heart contraction was coupled with the circulatory system, simulated with CircAdapt, while accounting for the effect of the pericardium on cardiac motion. The four-chamber electromechanics framework resulted in 117 parameters of interest. The model was broken into five hierarchical sub-models: tissue electrophysiology, ToR-ORd-Land model, Courtemanche-Land model, passive mechanics and CircAdapt. For each sub-model, we trained Gaussian processes emulators (GPEs) that were then used to perform a global sensitivity analysis (GSA) to retain parameters explaining 90% of the total sensitivity for subsequent analysis. We identified 45 out of 117 parameters that were important for whole heart function. We performed a GSA over these 45 parameters and identified the systemic and pulmonary peripheral resistance as being critical parameters for a wide range of volumetric and hemodynamic cardiac indexes across all four chambers. We have shown that GPEs provide a robust method for mapping between cellular properties and clinical measurements. This could be applied to identify parameters that can be calibrated in patient-specific models or digital twins, and to link cellular function to clinical indexes.
Collapse
Affiliation(s)
- Marina Strocchi
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
- National Heart and Lung Institute, Imperial College London, London, United Kingdom
| | - Stefano Longobardi
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
| | | | | | - Argyrios Petras
- Johann Radon Institute for Computational and Applied Mathematics (RICAM), Linz, Austria
| | - Christopher A Rinaldi
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
- Guy's and St Thomas' NHS Foundation Trust, London, United Kingdom
| | - Edward J Vigmond
- University of Bordeaux, CNRS, Bordeaux, Talence, France
- IHU Liryc, Bordeaux, Talence, France
| | - Gernot Plank
- Medical University of Graz, Graz, Austria
- BioTechMed-Graz, Graz, Austria
| | - Chris J Oates
- Newcastle University, Newcastle upon Tyne, United Kingdom
| | | | - Steven A Niederer
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
- National Heart and Lung Institute, Imperial College London, London, United Kingdom
- Alan Turing Institute, London, United Kingdom
| |
Collapse
|
6
|
Qian S, Monaci S, Mendonca-Costa C, Campos F, Gemmell P, Zaidi HA, Rajani R, Whitaker J, Rinaldi CA, Bishop MJ. Additional coils mitigate elevated defibrillation threshold in right-sided implantable cardioverter defibrillator generator placement: a simulation study. Europace 2023; 25:euad146. [PMID: 37314196 PMCID: PMC10265967 DOI: 10.1093/europace/euad146] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/20/2023] [Accepted: 05/13/2023] [Indexed: 06/15/2023] Open
Abstract
AIMS The standard implantable cardioverter defibrillator (ICD) generator (can) is placed in the left pectoral area; however, in certain circumstances, right-sided cans may be required which may increase defibrillation threshold (DFT) due to suboptimal shock vectors. We aim to quantitatively assess whether the potential increase in DFT of right-sided can configurations may be mitigated by alternate positioning of the right ventricular (RV) shocking coil or adding coils in the superior vena cava (SVC) and coronary sinus (CS). METHODS AND RESULTS A cohort of CT-derived torso models was used to assess DFT of ICD configurations with right-sided cans and alternate positioning of RV shock coils. Efficacy changes with additional coils in the SVC and CS were evaluated. A right-sided can with an apical RV shock coil significantly increased DFT compared to a left-sided can [19.5 (16.4, 27.1) J vs. 13.3 (11.7, 19.9) J, P < 0.001]. Septal positioning of the RV coil led to a further DFT increase when using a right-sided can [26.7 (18.1, 36.1) J vs. 19.5 (16.4, 27.1) J, P < 0.001], but not a left-sided can [12.1 (8.1, 17.6) J vs. 13.3 (11.7, 19.9) J, P = 0.099). Defibrillation threshold of a right-sided can with apical or septal coil was reduced the most by adding both SVC and CS coils [19.5 (16.4, 27.1) J vs. 6.6 (3.9, 9.9) J, P < 0.001, and 26.7 (18.1, 36.1) J vs. 12.1 (5.7, 13.5) J, P < 0.001]. CONCLUSION Right-sided, compared to left-sided, can positioning results in a 50% increase in DFT. For right-sided cans, apical shock coil positioning produces a lower DFT than septal positions. Elevated right-sided can DFTs may be mitigated by utilizing additional coils in SVC and CS.
Collapse
Affiliation(s)
- Shuang Qian
- Department of Biomedical Engineering, School of Imaging Sciences and Biomedical Engineering, Kings College London, 4th North Wing, St Thomas’ Hospital, London SE1 7EH, UK
| | - Sofia Monaci
- Department of Biomedical Engineering, School of Imaging Sciences and Biomedical Engineering, Kings College London, 4th North Wing, St Thomas’ Hospital, London SE1 7EH, UK
| | - Caroline Mendonca-Costa
- Department of Biomedical Engineering, School of Imaging Sciences and Biomedical Engineering, Kings College London, 4th North Wing, St Thomas’ Hospital, London SE1 7EH, UK
| | - Fernando Campos
- Department of Biomedical Engineering, School of Imaging Sciences and Biomedical Engineering, Kings College London, 4th North Wing, St Thomas’ Hospital, London SE1 7EH, UK
| | - Philip Gemmell
- Department of Biomedical Engineering, School of Imaging Sciences and Biomedical Engineering, Kings College London, 4th North Wing, St Thomas’ Hospital, London SE1 7EH, UK
| | - Hassan A Zaidi
- Department of Biomedical Engineering, School of Imaging Sciences and Biomedical Engineering, Kings College London, 4th North Wing, St Thomas’ Hospital, London SE1 7EH, UK
| | - Ronak Rajani
- Department of Cardiology, Guy’s and St Thomas’ Hospital, Westminster Bridge Rd, London SE1 7EH, UK
| | - John Whitaker
- Department of Biomedical Engineering, School of Imaging Sciences and Biomedical Engineering, Kings College London, 4th North Wing, St Thomas’ Hospital, London SE1 7EH, UK
- Department of Cardiology, Guy’s and St Thomas’ Hospital, Westminster Bridge Rd, London SE1 7EH, UK
| | - Christopher A Rinaldi
- Department of Biomedical Engineering, School of Imaging Sciences and Biomedical Engineering, Kings College London, 4th North Wing, St Thomas’ Hospital, London SE1 7EH, UK
- Department of Cardiology, Guy’s and St Thomas’ Hospital, Westminster Bridge Rd, London SE1 7EH, UK
| | - Martin J Bishop
- Department of Biomedical Engineering, School of Imaging Sciences and Biomedical Engineering, Kings College London, 4th North Wing, St Thomas’ Hospital, London SE1 7EH, UK
| |
Collapse
|
7
|
Alexander C, Bishop MJ, Gilchrist RJ, Burton FL, Smith GL, Myles RC. Initiation of ventricular arrhythmia in the acquired long QT syndrome. Cardiovasc Res 2023; 119:465-476. [PMID: 35727943 PMCID: PMC10064840 DOI: 10.1093/cvr/cvac103] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 01/10/2022] [Revised: 04/25/2022] [Accepted: 06/02/2022] [Indexed: 11/15/2022] Open
Abstract
AIMS Long QT syndrome (LQTS) carries a risk of life-threatening polymorphic ventricular tachycardia (Torsades de Pointes, TdP) and is a major cause of premature sudden cardiac death. TdP is induced by R-on-T premature ventricular complexes (PVCs), thought to be generated by cellular early-afterdepolarisations (EADs). However, EADs in tissue require cellular synchronisation, and their role in TdP induction remains unclear. We aimed to determine the mechanism of TdP induction in rabbit hearts with acquired LQTS (aLQTS). METHODS AND RESULTS Optical mapping of action potentials (APs) and intracellular Ca2+ was performed in Langendorff-perfused rabbit hearts (n = 17). TdP induced by R-on-T PVCs was observed during aLQTS (50% K+/Mg++ & E4031) conditions in all hearts (P < 0.0001 vs. control). Islands of AP prolongation bounded by steep voltage gradients (VGs) were consistently observed before arrhythmia and peak VGs were more closely related to the PVC upstroke than EADs, both temporally (7 ± 5 ms vs. 44 ± 27 ms, P < 0.0001) and spatially (1.0 ± 0.7 vs. 3.6 ± 0.9 mm, P < 0.0001). PVCs were initiated at estimated voltages of ∼ -40 mV and had upstroke dF/dtmax and Vm-Ca2+ dynamics compatible with ICaL activation. Computational simulations demonstrated that PVCs could arise directly from VGs, through electrotonic triggering of ICaL. In experiments and the model, sub-maximal L-type Ca2+ channel (LTCC) block (200 nM nifedipine and 90% gCaL, respectively) abolished both PVCs and TdP in the continued presence of aLQTS. CONCLUSION These data demonstrate that ICaL activation at sites displaying steep VGs generates the PVCs which induce TdP, providing a mechanism and rationale for LTCC blockers as a novel therapeutic approach in LQTS.
Collapse
Affiliation(s)
- Cherry Alexander
- Institute of Cardiovascular and Medical Sciences, University of Glasgow, Glasgow, G12 8QQ, UK
| | - Martin J Bishop
- School of Biomedical Engineering and Imaging Sciences, King’s College London, St Thomas' Hospital, Westminster Bridge Road, London SE1 7EH, UK
| | - Rebecca J Gilchrist
- Institute of Cardiovascular and Medical Sciences, University of Glasgow, Glasgow, G12 8QQ, UK
| | - Francis L Burton
- Institute of Cardiovascular and Medical Sciences, University of Glasgow, Glasgow, G12 8QQ, UK
| | - Godfrey L Smith
- Institute of Cardiovascular and Medical Sciences, University of Glasgow, Glasgow, G12 8QQ, UK
| | - Rachel C Myles
- Institute of Cardiovascular and Medical Sciences, University of Glasgow, Glasgow, G12 8QQ, UK
| |
Collapse
|
8
|
Qian S, Connolly A, Mendonca-Costa C, Campos F, Rodero C, Whitaker J, Rinaldi CA, Bishop MJ. Optimization of anti-tachycardia pacing efficacy through scar-specific delivery and minimization of re-initiation: a virtual study on a cohort of infarcted porcine hearts. Europace 2023; 25:716-725. [PMID: 36197749 PMCID: PMC9935023 DOI: 10.1093/europace/euac165] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2022] [Accepted: 08/31/2022] [Indexed: 11/15/2022] Open
Abstract
AIMS Anti-tachycardia pacing (ATP) is a reliable electrotherapy to painlessly terminate ventricular tachycardia (VT). However, ATP is often ineffective, particularly for fast VTs. The efficacy may be enhanced by optimized delivery closer to the re-entrant circuit driving the VT. This study aims to compare ATP efficacy for different delivery locations with respect to the re-entrant circuit, and further optimize ATP by minimizing failure through re-initiation. METHODS AND RESULTS Seventy-three sustained VTs were induced in a cohort of seven infarcted porcine ventricular computational models, largely dominated by a single re-entrant pathway. The efficacy of burst ATP delivered from three locations proximal to the re-entrant circuit (septum) and three distal locations (lateral/posterior left ventricle) was compared. Re-initiation episodes were used to develop an algorithm utilizing correlations between successive sensed electrogram morphologies to automatically truncate ATP pulse delivery. Anti-tachycardia pacing was more efficacious at terminating slow compared with fast VTs (65 vs. 46%, P = 0.000039). A separate analysis of slow VTs showed that the efficacy was significantly higher when delivered from distal compared with proximal locations (distal 72%, proximal 59%), being reversed for fast VTs (distal 41%, proximal 51%). Application of our early termination detection algorithm (ETDA) accurately detected VT termination in 79% of re-initiated cases, improving the overall efficacy for proximal delivery with delivery inside the critical isthmus (CI) itself being overall most effective. CONCLUSION Anti-tachycardia pacing delivery proximal to the re-entrant circuit is more effective at terminating fast VTs, but less so slow VTs, due to frequent re-initiation. Attenuating re-initiation, through ETDA, increases the efficacy of delivery within the CI for all VTs.
Collapse
Affiliation(s)
- Shuang Qian
- Department of Biomedical Engineering, School of Imaging Sciences and Biomedical Engineering, Kings College London, London, UK
| | | | - Caroline Mendonca-Costa
- Department of Biomedical Engineering, School of Imaging Sciences and Biomedical Engineering, Kings College London, London, UK
| | - Fernando Campos
- Department of Biomedical Engineering, School of Imaging Sciences and Biomedical Engineering, Kings College London, London, UK
| | - Cristobal Rodero
- Department of Biomedical Engineering, School of Imaging Sciences and Biomedical Engineering, Kings College London, London, UK
| | - John Whitaker
- Department of Biomedical Engineering, School of Imaging Sciences and Biomedical Engineering, Kings College London, London, UK
- Department of Cardiology, Guy’s and St Thomas’ Hospital, London, UK
| | - Christopher A Rinaldi
- Department of Biomedical Engineering, School of Imaging Sciences and Biomedical Engineering, Kings College London, London, UK
- Department of Cardiology, Guy’s and St Thomas’ Hospital, London, UK
| | - Martin J Bishop
- Department of Biomedical Engineering, School of Imaging Sciences and Biomedical Engineering, Kings College London, London, UK
| |
Collapse
|
9
|
Monaci S, Qian S, Gillette K, Puyol-Antón E, Mukherjee R, Elliott MK, Whitaker J, Rajani R, O’Neill M, Rinaldi CA, Plank G, King AP, Bishop MJ. Non-invasive localization of post-infarct ventricular tachycardia exit sites to guide ablation planning: a computational deep learning platform utilizing the 12-lead electrocardiogram and intracardiac electrograms from implanted devices. Europace 2023; 25:469-477. [PMID: 36369980 PMCID: PMC9935046 DOI: 10.1093/europace/euac178] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/21/2022] [Accepted: 09/08/2022] [Indexed: 11/13/2022] Open
Abstract
AIMS Existing strategies that identify post-infarct ventricular tachycardia (VT) ablation target either employ invasive electrophysiological (EP) mapping or non-invasive modalities utilizing the electrocardiogram (ECG). Their success relies on localizing sites critical to the maintenance of the clinical arrhythmia, not always recorded on the 12-lead ECG. Targeting the clinical VT by utilizing electrograms (EGM) recordings stored in implanted devices may aid ablation planning, enhancing safety and speed and potentially reducing the need of VT induction. In this context, we aim to develop a non-invasive computational-deep learning (DL) platform to localize VT exit sites from surface ECGs and implanted device intracardiac EGMs. METHODS AND RESULTS A library of ECGs and EGMs from simulated paced beats and representative post-infarct VTs was generated across five torso models. Traces were used to train DL algorithms to localize VT sites of earliest systolic activation; first tested on simulated data and then on a clinically induced VT to show applicability of our platform in clinical settings. Localization performance was estimated via localization errors (LEs) against known VT exit sites from simulations or clinical ablation targets. Surface ECGs successfully localized post-infarct VTs from simulated data with mean LE = 9.61 ± 2.61 mm across torsos. VT localization was successfully achieved from implanted device intracardiac EGMs with mean LE = 13.10 ± 2.36 mm. Finally, the clinically induced VT localization was in agreement with the clinical ablation volume. CONCLUSION The proposed framework may be utilized for direct localization of post-infarct VTs from surface ECGs and/or implanted device EGMs, or in conjunction with efficient, patient-specific modelling, enhancing safety and speed of ablation planning.
Collapse
Affiliation(s)
- Sofia Monaci
- Biomedical Engineering & Imaging Sciences, King’s College London, London SE1 7EH, United Kingdom
| | - Shuang Qian
- Biomedical Engineering & Imaging Sciences, King’s College London, London SE1 7EH, United Kingdom
| | | | - Esther Puyol-Antón
- Biomedical Engineering & Imaging Sciences, King’s College London, London SE1 7EH, United Kingdom
| | - Rahul Mukherjee
- Biomedical Engineering & Imaging Sciences, King’s College London, London SE1 7EH, United Kingdom
- Guy’s and St Thomas’ Hospital, London SE1 7EH, United Kingdom
| | - Mark K Elliott
- Biomedical Engineering & Imaging Sciences, King’s College London, London SE1 7EH, United Kingdom
- Guy’s and St Thomas’ Hospital, London SE1 7EH, United Kingdom
| | - John Whitaker
- Biomedical Engineering & Imaging Sciences, King’s College London, London SE1 7EH, United Kingdom
- Guy’s and St Thomas’ Hospital, London SE1 7EH, United Kingdom
| | - Ronak Rajani
- Biomedical Engineering & Imaging Sciences, King’s College London, London SE1 7EH, United Kingdom
- Guy’s and St Thomas’ Hospital, London SE1 7EH, United Kingdom
| | - Mark O’Neill
- Biomedical Engineering & Imaging Sciences, King’s College London, London SE1 7EH, United Kingdom
| | - Christopher A Rinaldi
- Biomedical Engineering & Imaging Sciences, King’s College London, London SE1 7EH, United Kingdom
- Guy’s and St Thomas’ Hospital, London SE1 7EH, United Kingdom
| | | | - Andrew P King
- Biomedical Engineering & Imaging Sciences, King’s College London, London SE1 7EH, United Kingdom
| | - Martin J Bishop
- Biomedical Engineering & Imaging Sciences, King’s College London, London SE1 7EH, United Kingdom
| |
Collapse
|
10
|
Kamali R, Gillete K, Tate J, Abhyankar DA, Dosdall DJ, Plank G, Bunch TJ, Macleod RS, Ranjan R. Treatment Planning for Atrial Fibrillation Using Patient-Specific Models Showing the Importance of Fibrillatory-Areas. Ann Biomed Eng 2023; 51:329-342. [PMID: 35930093 PMCID: PMC10440744 DOI: 10.1007/s10439-022-03029-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2022] [Accepted: 07/18/2022] [Indexed: 01/25/2023]
Abstract
Computational models have made it possible to study the effect of fibrosis and scar on atrial fibrillation (AF) and plan future personalized treatments. Here, we study the effect of area available for fibrillatory waves to sustain AF. Then we use it to plan for AF ablation to improve procedural outcomes. CARPentry was used to create patient-specific models to determine the association between the size of residual contiguous areas available for AF wavefronts to propagate and sustain AF [fibrillatory area (FA)] after ablation with procedural outcomes. The FA was quantified in a novel manner accounting for gaps in ablation lines. We selected 30 persistent AF patients with known ablation outcomes. We divided the atrial surface into five areas based on ablation scar pattern and anatomical landmarks and calculated the FAs. We validated the models based on clinical outcomes and suggested future ablation lines that minimize the FAs and terminate rotor activities in simulations. We also simulated the effects of three common antiarrhythmic drugs. In the patient-specific models, the predicted arrhythmias matched the clinical outcomes in 25 of 30 patients (accuracy 83.33%). The average largest FA (FAmax) in the recurrence group was 8517 ± 1444 vs. 6772 ± 1531 mm2 in the no recurrence group (p < 0.004). The final FAs after adding the suggested ablation lines in the AF recurrence group reduced the average FAmax from 8517 ± 1444 to 6168 ± 1358 mm2 (p < 0.001) and stopped the sustained rotor activity. Simulations also correctly anticipated the effect of antiarrhythmic drugs in 5 out of 6 patients who used drug therapy post unsuccessful ablation (accuracy 83.33%). Sizes of FAs available for AF wavefronts to propagate are important determinants for ablation outcomes. FA size in combination with computational simulations can be used to direct ablation in persistent AF to minimize the critical mass required to sustain recurrent AF.
Collapse
Affiliation(s)
- Roya Kamali
- Department of Biomedical Engineering, University of Utah, Salt Lake City, UT, USA
- Division of Cardiovascular Medicine, Department of Internal Medicine, University of Utah, Salt Lake City, UT, USA
- Nora Eccles Harrison Cardiovascular Research and Training Institute, University of Utah, Salt Lake City, UT, USA
- Department of Biomedical Engineering, Georgia Institute of Technology, Atlanta, GA, USA
| | - Karli Gillete
- Institute of Biophysics, Medical University of Graz, Graz, Austria
| | - Jess Tate
- Scientific Computing and Imaging Institute, University of Utah, Salt Lake City, UT, USA
| | | | - Derek J Dosdall
- Department of Biomedical Engineering, University of Utah, Salt Lake City, UT, USA
- Division of Cardiovascular Medicine, Department of Internal Medicine, University of Utah, Salt Lake City, UT, USA
- Nora Eccles Harrison Cardiovascular Research and Training Institute, University of Utah, Salt Lake City, UT, USA
| | - Gernot Plank
- Institute of Biophysics, Medical University of Graz, Graz, Austria
| | - T Jared Bunch
- Division of Cardiovascular Medicine, Department of Internal Medicine, University of Utah, Salt Lake City, UT, USA
| | - Rob S Macleod
- Department of Biomedical Engineering, University of Utah, Salt Lake City, UT, USA
- Nora Eccles Harrison Cardiovascular Research and Training Institute, University of Utah, Salt Lake City, UT, USA
| | - Ravi Ranjan
- Department of Biomedical Engineering, University of Utah, Salt Lake City, UT, USA.
- Division of Cardiovascular Medicine, Department of Internal Medicine, University of Utah, Salt Lake City, UT, USA.
- Nora Eccles Harrison Cardiovascular Research and Training Institute, University of Utah, Salt Lake City, UT, USA.
| |
Collapse
|
11
|
Ramlugun GS, Kulkarni K, Pallares-Lupon N, Boukens BJ, Efimov IR, Vigmond EJ, Bernus O, Walton RD. A comprehensive framework for evaluation of high pacing frequency and arrhythmic optical mapping signals. Front Physiol 2023; 14:734356. [PMID: 36755791 PMCID: PMC9901579 DOI: 10.3389/fphys.2023.734356] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2021] [Accepted: 01/09/2023] [Indexed: 01/24/2023] Open
Abstract
Introduction: High pacing frequency or irregular activity due to arrhythmia produces complex optical mapping signals and challenges for processing. The objective is to establish an automated activation time-based analytical framework applicable to optical mapping images of complex electrical behavior. Methods: Optical mapping signals with varying complexity from sheep (N = 7) ventricular preparations were examined. Windows of activation centered on each action potential upstroke were derived using Hilbert transform phase. Upstroke morphology was evaluated for potential multiple activation components and peaks of upstroke signal derivatives defined activation time. Spatially and temporally clustered activation time points were grouped in to wave fronts for individual processing. Each activation time point was evaluated for corresponding repolarization times. Each wave front was subsequently classified based on repetitive or non-repetitive events. Wave fronts were evaluated for activation time minima defining sites of wave front origin. A visualization tool was further developed to probe dynamically the ensemble activation sequence. Results: Our framework facilitated activation time mapping during complex dynamic events including transitions to rotor-like reentry and ventricular fibrillation. We showed that using fixed AT windows to extract AT maps can impair interpretation of the activation sequence. However, the phase windowing of action potential upstrokes enabled accurate recapitulation of repetitive behavior, providing spatially coherent activation patterns. We further demonstrate that grouping the spatio-temporal distribution of AT points in to coherent wave fronts, facilitated interpretation of isolated conduction events, such as conduction slowing, and to derive dynamic changes in repolarization properties. Focal origins precisely detected sites of stimulation origin and breakthrough for individual wave fronts. Furthermore, a visualization tool to dynamically probe activation time windows during reentry revealed a critical single static line of conduction slowing associated with the rotation core. Conclusion: This comprehensive analytical framework enables detailed quantitative assessment and visualization of complex electrical behavior.
Collapse
Affiliation(s)
- Girish S. Ramlugun
- IHU-Liryc, Fondation Bordeaux Université, Pessac-Bordeaux, France,Univ. Bordeaux, Inserm, Centre de Recherche Cardio-Thoracique, Bordeaux, France
| | - Kanchan Kulkarni
- IHU-Liryc, Fondation Bordeaux Université, Pessac-Bordeaux, France,Univ. Bordeaux, Inserm, Centre de Recherche Cardio-Thoracique, Bordeaux, France
| | - Nestor Pallares-Lupon
- IHU-Liryc, Fondation Bordeaux Université, Pessac-Bordeaux, France,Univ. Bordeaux, Inserm, Centre de Recherche Cardio-Thoracique, Bordeaux, France
| | - Bastiaan J. Boukens
- Department of Physiology, Cardiovascular Research Institute Maastricht, University Maastricht, Maastricht, Netherlands,Department of Medical Biology, Amsterdam Cardiovascular Sciences, Amsterdam University Medical Center, University of Amsterdam, Amsterdam, Netherlands
| | - Igor R. Efimov
- Department of Biomedical Engineering, The George Washington University, Washington, DC, United States,Department of Biomedical Engineering, Northwestern University, Chicago, IL, United States,Department of Medicine, Northwestern University, Chicago, IL, United States
| | - Edward J. Vigmond
- IHU-Liryc, Fondation Bordeaux Université, Pessac-Bordeaux, France,Univ. Bordeaux, Centre National de la Recherche Scientifique (CNRS), Institut de Mathématiques de Bordeaux, UMR5251, Bordeaux, France
| | - Olivier Bernus
- IHU-Liryc, Fondation Bordeaux Université, Pessac-Bordeaux, France,Univ. Bordeaux, Inserm, Centre de Recherche Cardio-Thoracique, Bordeaux, France
| | - Richard D. Walton
- IHU-Liryc, Fondation Bordeaux Université, Pessac-Bordeaux, France,Univ. Bordeaux, Inserm, Centre de Recherche Cardio-Thoracique, Bordeaux, France,*Correspondence: Richard D. Walton,
| |
Collapse
|
12
|
Calibration of Cohorts of Virtual Patient Heart Models Using Bayesian History Matching. Ann Biomed Eng 2023; 51:241-252. [PMID: 36271218 PMCID: PMC9832095 DOI: 10.1007/s10439-022-03095-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2022] [Accepted: 09/29/2022] [Indexed: 01/28/2023]
Abstract
Previous patient-specific model calibration techniques have treated each patient independently, making the methods expensive for large-scale clinical adoption. In this work, we show how we can reuse simulations to accelerate the patient-specific model calibration pipeline. To represent anatomy, we used a Statistical Shape Model and to represent function, we ran electrophysiological simulations. We study the use of 14 biomarkers to calibrate the model, training one Gaussian Process Emulator (GPE) per biomarker. To fit the models, we followed a Bayesian History Matching (BHM) strategy, wherein each iteration a region of the parameter space is ruled out if the emulation with that set of parameter values produces is "implausible". We found that without running any extra simulations we can find 87.41% of the non-implausible parameter combinations. Moreover, we showed how reducing the uncertainty of the measurements from 10 to 5% can reduce the final parameter space by 6 orders of magnitude. This innovation allows for a model fitting technique, therefore reducing the computational load of future biomedical studies.
Collapse
|
13
|
Li G, Brumback BD, Huang L, Zhang DM, Yin T, Lipovsky CE, Hicks SC, Jimenez J, Boyle PM, Rentschler SL. Acute Glycogen Synthase Kinase-3 Inhibition Modulates Human Cardiac Conduction. JACC Basic Transl Sci 2022; 7:1001-1017. [PMID: 36337924 PMCID: PMC9626903 DOI: 10.1016/j.jacbts.2022.04.007] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 07/20/2021] [Revised: 04/11/2022] [Accepted: 04/13/2022] [Indexed: 01/14/2023]
Abstract
Glycogen synthase kinase 3 (GSK-3) inhibition has emerged as a potential therapeutic target for several diseases, including cancer. However, the role for GSK-3 regulation of human cardiac electrophysiology remains ill-defined. We demonstrate that SB216763, a GSK-3 inhibitor, can acutely reduce conduction velocity in human cardiac slices. Combined computational modeling and experimental approaches provided mechanistic insight into GSK-3 inhibition-mediated changes, revealing that decreased sodium-channel conductance and tissue conductivity may underlie the observed phenotypes. Our study demonstrates that GSK-3 inhibition in human myocardium alters electrophysiology and may predispose to an arrhythmogenic substrate; therefore, monitoring for adverse arrhythmogenic events could be considered.
Collapse
Key Words
- ABC, active β-catenin
- APD, action potential duration
- BDM, 2,3-butanedione monoxime
- CV, conduction velocity
- Cx43, connexin 43
- GNa, sodium-channel conductance
- GOF, gain of function
- GSK-3 inhibitor
- GSK-3, glycogen synthase kinase 3
- INa, sodium current
- LV, left ventricle
- NaV1.5, pore-forming α-subunit protein of the voltage-gated cardiac sodium channel
- PCR, polymerase chain reaction
- RMP, resting membrane potential
- RT-qPCR, reverse transcription-quantitative polymerase chain reaction
- SB2, SB216763
- SB216763
- cDNA, complementary DNA
- dVm/dtmax, maximum upstroke velocity
- electrophysiology
- human cardiac slices
Collapse
Affiliation(s)
- Gang Li
- Department of Medicine, Cardiovascular Division, Washington University School of Medicine in St. Louis, Missouri, USA
- Department of Biomedical Engineering, Washington University McKelvey School of Engineering in St. Louis, Missouri, USA
| | - Brittany D. Brumback
- Department of Medicine, Cardiovascular Division, Washington University School of Medicine in St. Louis, Missouri, USA
- Department of Biomedical Engineering, Washington University McKelvey School of Engineering in St. Louis, Missouri, USA
| | - Lei Huang
- Department of Medicine, Cardiovascular Division, Washington University School of Medicine in St. Louis, Missouri, USA
| | - David M. Zhang
- Department of Medicine, Cardiovascular Division, Washington University School of Medicine in St. Louis, Missouri, USA
| | - Tiankai Yin
- Department of Medicine, Cardiovascular Division, Washington University School of Medicine in St. Louis, Missouri, USA
| | - Catherine E. Lipovsky
- Department of Medicine, Cardiovascular Division, Washington University School of Medicine in St. Louis, Missouri, USA
- Department of Developmental Biology, Washington University School of Medicine in St. Louis, Missouri, USA
| | - Stephanie C. Hicks
- Department of Medicine, Cardiovascular Division, Washington University School of Medicine in St. Louis, Missouri, USA
| | - Jesus Jimenez
- Department of Medicine, Cardiovascular Division, Washington University School of Medicine in St. Louis, Missouri, USA
| | - Patrick M. Boyle
- Department of Bioengineering, Center for Cardiovascular Biology, and Institute for Stem Cell and Regenerative Medicine, University of Washington, Seattle, Washington, USA
| | - Stacey L. Rentschler
- Department of Medicine, Cardiovascular Division, Washington University School of Medicine in St. Louis, Missouri, USA
- Department of Biomedical Engineering, Washington University McKelvey School of Engineering in St. Louis, Missouri, USA
- Department of Developmental Biology, Washington University School of Medicine in St. Louis, Missouri, USA
| |
Collapse
|
14
|
Nazari-Vanani R, Mohammadpour R, Asadian E, Rafii-Tabar H, Sasanpour P. A computational modelling study of excitation of neuronal cells with triboelectric nanogenerators. Sci Rep 2022; 12:13411. [PMID: 35927441 PMCID: PMC9352766 DOI: 10.1038/s41598-022-17050-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2022] [Accepted: 07/20/2022] [Indexed: 11/09/2022] Open
Abstract
Neurological disorders and nerve injuries, such as spinal cord injury, stroke, and multiple sclerosis can result in the loss of muscle function. Electrical stimulation of the neuronal cells is the currently available clinical treatment in this regard. As an effective energy harvester, the triboelectric nanogenerators (TENG) can be used for self-powered neural/muscle stimulations because the output of the TENG provides stimulation pulses for nerves. In the present study, using a computational modelling approach, the effect of surface micropatterns on the electric field distribution, induced voltage and capacitance of the TENG structures have been investigated. By incorporating the effect of the TENG inside the mathematical model of neuron’s electrical behavior (cable equation with Hodgkin-Huxley model), its impact on the electrical behavior of the neurons has been studied. The results show that the TENG operates differently with various surface modifications. The performance of the TENG in excitation of neurons depends on the contact and release speed of its electrodes accordingly.
Collapse
Affiliation(s)
- Razieh Nazari-Vanani
- Department of Medical Physics and Biomedical Engineering, School of Medicine, Shahid Beheshti University of Medical Sciences, Tehran, Iran
| | - Raheleh Mohammadpour
- Institute for Nanoscience and Nanotechnology (INST), Sharif University of Technology, Tehran, Iran.
| | - Elham Asadian
- Department of Medical Physics and Biomedical Engineering, School of Medicine, Shahid Beheshti University of Medical Sciences, Tehran, Iran
| | - Hashem Rafii-Tabar
- Department of Medical Physics and Biomedical Engineering, School of Medicine, Shahid Beheshti University of Medical Sciences, Tehran, Iran.,The Physics Branch of Iran Academy of Sciences, Tehran, Iran
| | - Pezhman Sasanpour
- Department of Medical Physics and Biomedical Engineering, School of Medicine, Shahid Beheshti University of Medical Sciences, Tehran, Iran. .,School of Nanoscience, Institute for Research in Fundamental Sciences (IPM), P. O. Box 19395-5531, Tehran, Iran.
| |
Collapse
|
15
|
Campos FO, Neic A, Mendonca Costa C, Whitaker J, O'Neill M, Razavi R, Rinaldi CA, DanielScherr, Niederer SA, Plank G, Bishop MJ. An automated near-real time computational method for induction and treatment of scar-related ventricular tachycardias. Med Image Anal 2022; 80:102483. [PMID: 35667328 PMCID: PMC10114098 DOI: 10.1016/j.media.2022.102483] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2021] [Revised: 04/22/2022] [Accepted: 05/20/2022] [Indexed: 02/05/2023]
Abstract
Catheter ablation is currently the only curative treatment for scar-related ventricular tachycardias (VTs). However, not only are ablation procedures long, with relatively high risk, but success rates are punitively low, with frequent VT recurrence. Personalized in-silico approaches have the opportunity to address these limitations. However, state-of-the-art reaction diffusion (R-D) simulations of VT induction and subsequent circuits used for in-silico ablation target identification require long execution times, along with vast computational resources, which are incompatible with the clinical workflow. Here, we present the Virtual Induction and Treatment of Arrhythmias (VITA), a novel, rapid and fully automated computational approach that uses reaction-Eikonal methodology to induce VT and identify subsequent ablation targets. The rationale for VITA is based on finding isosurfaces associated with an activation wavefront that splits in the ventricles due to the presence of an isolated isthmus of conduction within the scar; once identified, each isthmus may be assessed for their vulnerability to sustain a reentrant circuit, and the corresponding exit site automatically identified for potential ablation targeting. VITA was tested on a virtual cohort of 7 post-infarcted porcine hearts and the results compared to R-D simulations. Using only a standard desktop machine, VITA could detect all scar-related VTs, simulating activation time maps and ECGs (for clinical comparison) as well as computing ablation targets in 48 minutes. The comparable VTs probed by the R-D simulations took 68.5 hours on 256 cores of high-performance computing infrastructure. The set of lesions computed by VITA was shown to render the ventricular model VT-free. VITA could be used in near real-time as a complementary modality aiding in clinical decision-making in the treatment of post-infarction VTs.
Collapse
Affiliation(s)
- Fernando O Campos
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom.
| | | | - Caroline Mendonca Costa
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
| | - John Whitaker
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom; Guy's and St. Thomas' NHS Foundation Trust, Cardiovascular Directorate
| | - Mark O'Neill
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom; Guy's and St. Thomas' NHS Foundation Trust, Cardiovascular Directorate
| | - Reza Razavi
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
| | - Christopher A Rinaldi
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom; Guy's and St. Thomas' NHS Foundation Trust, Cardiovascular Directorate
| | - DanielScherr
- Division of Cardiology, Department of Internal Medicine, Medical University of Graz, Austria
| | - Steven A Niederer
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
| | - Gernot Plank
- Gottfried Schatz Research Center Biophysics, Medical University of Graz, Graz, Austria; BioTechMed-Graz, Graz, Austria
| | - Martin J Bishop
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
| |
Collapse
|
16
|
Sánchez J, Loewe A. A Review of Healthy and Fibrotic Myocardium Microstructure Modeling and Corresponding Intracardiac Electrograms. Front Physiol 2022; 13:908069. [PMID: 35620600 PMCID: PMC9127661 DOI: 10.3389/fphys.2022.908069] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2022] [Accepted: 04/21/2022] [Indexed: 11/13/2022] Open
Abstract
Computational simulations of cardiac electrophysiology provide detailed information on the depolarization phenomena at different spatial and temporal scales. With the development of new hardware and software, in silico experiments have gained more importance in cardiac electrophysiology research. For plane waves in healthy tissue, in vivo and in silico electrograms at the surface of the tissue demonstrate symmetric morphology and high peak-to-peak amplitude. Simulations provided insight into the factors that alter the morphology and amplitude of the electrograms. The situation is more complex in remodeled tissue with fibrotic infiltrations. Clinically, different changes including fractionation of the signal, extended duration and reduced amplitude have been described. In silico, numerous approaches have been proposed to represent the pathological changes on different spatial and functional scales. Different modeling approaches can reproduce distinct subsets of the clinically observed electrogram phenomena. This review provides an overview of how different modeling approaches to incorporate fibrotic and structural remodeling affect the electrogram and highlights open challenges to be addressed in future research.
Collapse
|
17
|
An Automata-Based Cardiac Electrophysiology Simulator to Assess Arrhythmia Inducibility. MATHEMATICS 2022. [DOI: 10.3390/math10081293] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
Abstract
Personalized cardiac electrophysiology simulations have demonstrated great potential to study cardiac arrhythmias and help in therapy planning of radio-frequency ablation. Its application to analyze vulnerability to ventricular tachycardia and sudden cardiac death in infarcted patients has been recently explored. However, the detailed multi-scale biophysical simulations used in these studies are very demanding in terms of memory and computational resources, which prevents their clinical translation. In this work, we present a fast phenomenological system based on cellular automata (CA) to simulate personalized cardiac electrophysiology. The system is trained on biophysical simulations to reproduce cellular and tissue dynamics in healthy and pathological conditions, including action potential restitution, conduction velocity restitution and cell safety factor. We show that a full ventricular simulation can be performed in the order of seconds, emulate the results of a biophysical simulation and reproduce a patient’s ventricular tachycardia in a model that includes a heterogeneous scar region. The system could be used to study the risk of arrhythmia in infarcted patients for a large number of scenarios.
Collapse
|
18
|
Feng Y, Roney CH, Bayer JD, Niederer SA, Hocini M, Vigmond EJ. Detection of focal source and arrhythmogenic substrate from body surface potentials to guide atrial fibrillation ablation. PLoS Comput Biol 2022; 18:e1009893. [PMID: 35312675 PMCID: PMC8970486 DOI: 10.1371/journal.pcbi.1009893] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2021] [Revised: 03/31/2022] [Accepted: 02/02/2022] [Indexed: 11/18/2022] Open
Abstract
Focal sources (FS) are believed to be important triggers and a perpetuation mechanism for paroxysmal atrial fibrillation (AF). Detecting FS and determining AF sustainability in atrial tissue can help guide ablation targeting. We hypothesized that sustained rotors during FS-driven episodes indicate an arrhythmogenic substrate for sustained AF, and that non-invasive electrical recordings, like electrocardiograms (ECGs) or body surface potential maps (BSPMs), could be used to detect FS and AF sustainability. Computer simulations were performed on five bi-atrial geometries. FS were induced by pacing at cycle lengths of 120-270 ms from 32 atrial sites and four pulmonary veins. Self-sustained reentrant activities were also initiated around the same 32 atrial sites with inexcitable cores of radii of 0, 0.5 and 1 cm. FS fired for two seconds and then AF inducibility was tested by whether activation was sustained for another second. ECGs and BSPMs were simulated. Equivalent atrial sources were extracted using second-order blind source separation, and their cycle length, periodicity and contribution, were used as features for random forest classifiers. Longer rotor duration during FS-driven episodes indicates higher AF inducibility (area under ROC curve = 0.83). Our method had accuracy of 90.6±1.0% and 90.6±0.6% in detecting FS presence, and 93.1±0.6% and 94.2±1.2% in identifying AF sustainability, and 80.0±6.6% and 61.0±5.2% in determining the atrium of the focal site, from BSPMs and ECGs of five atria. The detection of FS presence and AF sustainability were insensitive to vest placement (±9.6%). On pre-operative BSPMs of 52 paroxysmal AF patients, patients classified with initiator-type FS on a single atrium resulted in improved two-to-three-year AF-free likelihoods (p-value < 0.01, logrank tests). Detection of FS and arrhythmogenic substrate can be performed from ECGs and BSPMs, enabling non-invasive mapping towards mechanism-targeted AF treatment, and malignant ectopic beat detection with likely AF progression.
Collapse
Affiliation(s)
- Yingjing Feng
- IHU Liryc, Electrophysiology and Heart Modeling Institute, fondation Bordeaux Université, Pessac-Bordeaux, France
- Univ. Bordeaux, IMB, UMR 5251, Talence, France
| | - Caroline H. Roney
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom
| | - Jason D. Bayer
- IHU Liryc, Electrophysiology and Heart Modeling Institute, fondation Bordeaux Université, Pessac-Bordeaux, France
- Univ. Bordeaux, IMB, UMR 5251, Talence, France
| | - Steven A. Niederer
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom
| | - Mélèze Hocini
- IHU Liryc, Electrophysiology and Heart Modeling Institute, fondation Bordeaux Université, Pessac-Bordeaux, France
- Bordeaux University Hospital (CHU), Electrophysiology and Ablation Unit, Pessac, France
| | - Edward J. Vigmond
- IHU Liryc, Electrophysiology and Heart Modeling Institute, fondation Bordeaux Université, Pessac-Bordeaux, France
- Univ. Bordeaux, IMB, UMR 5251, Talence, France
| |
Collapse
|
19
|
Electro-anatomical computational cardiology in humans and experimental animal models. TRANSLATIONAL RESEARCH IN ANATOMY 2022. [DOI: 10.1016/j.tria.2022.100162] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022] Open
|
20
|
Mendonca Costa C, Gemmell P, Elliott MK, Whitaker J, Campos FO, Strocchi M, Neic A, Gillette K, Vigmond E, Plank G, Razavi R, O'Neill M, Rinaldi CA, Bishop MJ. Determining anatomical and electrophysiological detail requirements for computational ventricular models of porcine myocardial infarction. Comput Biol Med 2022; 141:105061. [PMID: 34915331 PMCID: PMC8819160 DOI: 10.1016/j.compbiomed.2021.105061] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2021] [Revised: 11/04/2021] [Accepted: 11/20/2021] [Indexed: 12/01/2022]
Abstract
BACKGROUND Computational models of the heart built from cardiac MRI and electrophysiology (EP) data have shown promise for predicting the risk of and ablation targets for myocardial infarction (MI) related ventricular tachycardia (VT), as well as to predict paced activation sequences in heart failure patients. However, most recent studies have relied on low resolution imaging data and little or no EP personalisation, which may affect the accuracy of model-based predictions. OBJECTIVE To investigate the impact of model anatomy, MI scar morphology, and EP personalisation strategies on paced activation sequences and VT inducibility to determine the level of detail required to make accurate model-based predictions. METHODS Imaging and EP data were acquired from a cohort of six pigs with experimentally induced MI. Computational models of ventricular anatomy, incorporating MI scar, were constructed including bi-ventricular or left ventricular (LV) only anatomy, and MI scar morphology with varying detail. Tissue conductivities and action potential duration (APD) were fitted to 12-lead ECG data using the QRS duration and the QT interval, respectively, in addition to corresponding literature parameters. Paced activation sequences and VT induction were simulated. Simulated paced activation and VT inducibility were compared between models and against experimental data. RESULTS Simulations predict that the level of model anatomical detail has little effect on simulated paced activation, with all model predictions comparing closely with invasive EP measurements. However, detailed scar morphology from high-resolution images, bi-ventricular anatomy, and personalized tissue conductivities are required to predict experimental VT outcome. CONCLUSION This study provides clear guidance for model generation based on clinical data. While a representing high level of anatomical and scar detail will require high-resolution image acquisition, EP personalisation based on 12-lead ECG can be readily incorporated into modelling pipelines, as such data is widely available.
Collapse
Affiliation(s)
- Caroline Mendonca Costa
- Department of Biomedical Engineering, School of Biomedical Engineering & Imaging Sciences, King's College London, UK.
| | - Philip Gemmell
- Department of Biomedical Engineering, School of Biomedical Engineering & Imaging Sciences, King's College London, UK
| | - Mark K Elliott
- Department of Biomedical Engineering, School of Biomedical Engineering & Imaging Sciences, King's College London, UK
| | - John Whitaker
- Department of Biomedical Engineering, School of Biomedical Engineering & Imaging Sciences, King's College London, UK
| | - Fernando O Campos
- Department of Biomedical Engineering, School of Biomedical Engineering & Imaging Sciences, King's College London, UK
| | - Marina Strocchi
- Department of Biomedical Engineering, School of Biomedical Engineering & Imaging Sciences, King's College London, UK
| | | | - Karli Gillette
- Gottfried Schatz Research Center, Biophysics, Medical University of Graz, Austria; Medical University of Graz, Austria and BioTechMed, Graz, Austria
| | - Edward Vigmond
- Institut de Rythmologie et de modélisation cardiaque (LIRYC), University of Bordeaux, France
| | - Gernot Plank
- Medical University of Graz, Austria and BioTechMed, Graz, Austria
| | - Reza Razavi
- Department of Biomedical Engineering, School of Biomedical Engineering & Imaging Sciences, King's College London, UK
| | - Mark O'Neill
- Department of Cardiology, Guy's and St Thomas' Hospital, London, UK
| | - Christopher A Rinaldi
- Department of Biomedical Engineering, School of Biomedical Engineering & Imaging Sciences, King's College London, UK; Department of Cardiology, Guy's and St Thomas' Hospital, London, UK
| | - Martin J Bishop
- Department of Biomedical Engineering, School of Biomedical Engineering & Imaging Sciences, King's College London, UK
| |
Collapse
|
21
|
Multiscale Modeling of the Mitochondrial Origin of Cardiac Reentrant and Fibrillatory Arrhythmias. Methods Mol Biol 2022; 2399:247-259. [PMID: 35604560 PMCID: PMC10186263 DOI: 10.1007/978-1-0716-1831-8_11] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Abstract
While mitochondrial dysfunction has been implicated in the pathogenesis of cardiac arrhythmias, how the abnormality occurring at the organelle level escalates to influence the rhythm of the heart remains incompletely understood. This is due, in part, to the complexity of the interactions formed by cardiac electrical, mechanical, and metabolic subsystems at various spatiotemporal scales that is difficult to fully comprehend solely with experiments. Computational models have emerged as a powerful tool to explore complicated and highly dynamic biological systems such as the heart, alone or in combination with experimental measurements. Here, we describe a strategy of integrating computer simulations with optical mapping of cardiomyocyte monolayers to examine how regional mitochondrial dysfunction elicits abnormal electrical activity, such as rebound and spiral waves, leading to reentry and fibrillation in cardiac tissue. We anticipate that this advanced modeling technology will enable new insights into the mechanisms by which changes in subcellular organelles can impact organ function.
Collapse
|
22
|
Qian S, Connolly A, Mendonca-Costa C, Campos F, Williams SE, Whitaker J, Rinaldi CA, Bishop MJ. An in-silico assessment of efficacy of two novel intra-cardiac electrode configurations versus traditional anti-tachycardia pacing therapy for terminating sustained ventricular tachycardia. Comput Biol Med 2021; 139:104987. [PMID: 34741904 PMCID: PMC8669079 DOI: 10.1016/j.compbiomed.2021.104987] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2021] [Revised: 10/24/2021] [Accepted: 10/24/2021] [Indexed: 11/06/2022]
Abstract
The implanted cardioverter defibrillator (ICD) is an effective direct therapy for the treatment of cardiac arrhythmias, including ventricular tachycardia (VT). Anti-tachycardia pacing (ATP) is often applied by the ICD as the first mode of therapy, but is often found to be ineffective, particularly for fast VTs. In such cases, strong, painful and damaging backup defibrillation shocks are applied by the device. Here, we propose two novel electrode configurations: "bipolar" and "transmural" which both combine the concept of targeted shock delivery with the advantage of reduced energy required for VT termination. We perform an in silico study to evaluate the efficacy of VT termination by applying one single (low-energy) monophasic shock from each novel configuration, comparing with conventional ATP therapy. Both bipolar and transmural configurations are able to achieve a higher efficacy (93% and 85%) than ATP (45%), with energy delivered similar to and two orders of magnitudes smaller than conventional ICD defibrillation shocks, respectively. Specifically, the transmural configuration (which applies the shock vector directly across the scar substrate sustaining the VT) is most efficient, requiring typically less than 1 J shock energy to achieve a high efficacy. The efficacy of both bipolar and transmural configurations are higher when applied to slow VTs (100% and 97%) compared to fast VTs (57% and 29%). Both novel electrode configurations introduced are able to improve electrotherapy efficacy while reducing the overall number of required therapies and need for strong backup shocks.
Collapse
Affiliation(s)
- Shuang Qian
- School of Biomedical Engineering and Imaging Sciences, Rayne Institute, King's College London, 4th Floor, Lambeth Wing, St. Thomas' Hospital, Westminster Bridge Road, London, SE1 7EH, United Kingdom.
| | - Adam Connolly
- Invicro, Burlington Danes Building, Du Cane Rd, London, W12 0N, United Kingdom
| | - Caroline Mendonca-Costa
- School of Biomedical Engineering and Imaging Sciences, Rayne Institute, King's College London, 4th Floor, Lambeth Wing, St. Thomas' Hospital, Westminster Bridge Road, London, SE1 7EH, United Kingdom
| | - Fernando Campos
- School of Biomedical Engineering and Imaging Sciences, Rayne Institute, King's College London, 4th Floor, Lambeth Wing, St. Thomas' Hospital, Westminster Bridge Road, London, SE1 7EH, United Kingdom
| | - Steven E Williams
- School of Biomedical Engineering and Imaging Sciences, Rayne Institute, King's College London, 4th Floor, Lambeth Wing, St. Thomas' Hospital, Westminster Bridge Road, London, SE1 7EH, United Kingdom
| | - John Whitaker
- School of Biomedical Engineering and Imaging Sciences, Rayne Institute, King's College London, 4th Floor, Lambeth Wing, St. Thomas' Hospital, Westminster Bridge Road, London, SE1 7EH, United Kingdom; Department of Cardiology, Guy's and St Thomas' Hospital, London, SE1 7EH, United Kingdom
| | - Christopher A Rinaldi
- School of Biomedical Engineering and Imaging Sciences, Rayne Institute, King's College London, 4th Floor, Lambeth Wing, St. Thomas' Hospital, Westminster Bridge Road, London, SE1 7EH, United Kingdom; Department of Cardiology, Guy's and St Thomas' Hospital, London, SE1 7EH, United Kingdom
| | - Martin J Bishop
- School of Biomedical Engineering and Imaging Sciences, Rayne Institute, King's College London, 4th Floor, Lambeth Wing, St. Thomas' Hospital, Westminster Bridge Road, London, SE1 7EH, United Kingdom
| |
Collapse
|
23
|
Roney CH, Child N, Porter B, Sim I, Whitaker J, Clayton RH, Laughner JI, Shuros A, Neuzil P, Williams SE, Razavi RS, O'Neill M, Rinaldi CA, Taggart P, Wright M, Gill JS, Niederer SA. Time-Averaged Wavefront Analysis Demonstrates Preferential Pathways of Atrial Fibrillation, Predicting Pulmonary Vein Isolation Acute Response. Front Physiol 2021; 12:707189. [PMID: 34646149 PMCID: PMC8503618 DOI: 10.3389/fphys.2021.707189] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2021] [Accepted: 08/24/2021] [Indexed: 11/13/2022] Open
Abstract
Electrical activation during atrial fibrillation (AF) appears chaotic and disorganised, which impedes characterisation of the underlying substrate and treatment planning. While globally chaotic, there may be local preferential activation pathways that represent potential ablation targets. This study aimed to identify preferential activation pathways during AF and predict the acute ablation response when these are targeted by pulmonary vein isolation (PVI). In patients with persistent AF (n = 14), simultaneous biatrial contact mapping with basket catheters was performed pre-ablation and following each ablation strategy (PVI, roof, and mitral lines). Unipolar wavefront activation directions were averaged over 10 s to identify preferential activation pathways. Clinical cases were classified as responders or non-responders to PVI during the procedure. Clinical data were augmented with a virtual cohort of 100 models. In AF pre-ablation, pathways originated from the pulmonary vein (PV) antra in PVI responders (7/7) but not in PVI non-responders (6/6). We proposed a novel index that measured activation waves from the PV antra into the atrial body. This index was significantly higher in PVI responders than non-responders (clinical: 16.3 vs. 3.7%, p = 0.04; simulated: 21.1 vs. 14.1%, p = 0.02). Overall, this novel technique and proof of concept study demonstrated that preferential activation pathways exist during AF. Targeting patient-specific activation pathways that flowed from the PV antra to the left atrial body using PVI resulted in AF termination during the procedure. These PV activation flow pathways may correspond to the presence of drivers in the PV regions.
Collapse
Affiliation(s)
- Caroline H. Roney
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
| | - Nicholas Child
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
| | - Bradley Porter
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
| | - Iain Sim
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
| | - John Whitaker
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
| | - Richard H. Clayton
- INSIGNEO Institute for In Silico Medicine and Department of Computer Science, University of Sheffield, Sheffield, United Kingdom
| | | | - Allan Shuros
- Boston Scientific Corp, St. Paul, MN, United States
| | - Petr Neuzil
- Department of Cardiology, Na Holmolce Hospital, Prague, Czechia
| | - Steven E. Williams
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
| | - Reza S. Razavi
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
| | - Mark O'Neill
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
| | | | - Peter Taggart
- Institute of Cardiovascular Science, University College London, London, United Kingdom
| | - Matt Wright
- Department of Cardiology, Guy's and St Thomas' Hospital, London, United Kingdom
| | - Jaswinder S. Gill
- Department of Cardiology, Guy's and St Thomas' Hospital, London, United Kingdom
| | - Steven A. Niederer
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
| |
Collapse
|
24
|
Ochs AR, Karathanos TV, Trayanova NA, Boyle PM. Optogenetic Stimulation Using Anion Channelrhodopsin (GtACR1) Facilitates Termination of Reentrant Arrhythmias With Low Light Energy Requirements: A Computational Study. Front Physiol 2021; 12:718622. [PMID: 34526912 PMCID: PMC8435849 DOI: 10.3389/fphys.2021.718622] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2021] [Accepted: 07/23/2021] [Indexed: 12/24/2022] Open
Abstract
Optogenetic defibrillation of hearts expressing light-sensitive cation channels (e.g., ChR2) has been proposed as an alternative to conventional electrotherapy. Past modeling work has shown that ChR2 stimulation can depolarize enough myocardium to interrupt arrhythmia, but its efficacy is limited by light attenuation and high energy needs. These shortcomings may be mitigated by using new optogenetic proteins like Guillardia theta Anion Channelrhodopsin (GtACR1), which produces a repolarizing outward current upon illumination. Accordingly, we designed a study to assess the feasibility of GtACR1-based optogenetic arrhythmia termination in human hearts. We conducted electrophysiological simulations in MRI-based atrial or ventricular models (n = 3 each), with pathological remodeling from atrial fibrillation or ischemic cardiomyopathy, respectively. We simulated light sensitization via viral gene delivery of three different opsins (ChR2, red-shifted ChR2, GtACR1) and uniform endocardial illumination at the appropriate wavelengths (blue, red, or green light, respectively). To analyze consistency of arrhythmia termination, we varied pulse timing (three evenly spaced intervals spanning the reentrant cycle) and intensity (atrial: 0.001–1 mW/mm2; ventricular: 0.001–10 mW/mm2). In atrial models, GtACR1 stimulation with 0.005 mW/mm2 green light consistently terminated reentry; this was 10–100x weaker than the threshold levels for ChR2-mediated defibrillation. In ventricular models, defibrillation was observed in 2/3 models for GtACR1 stimulation at 0.005 mW/mm2 (100–200x weaker than ChR2 cases). In the third ventricular model, defibrillation failed in nearly all cases, suggesting that attenuation issues and patient-specific organ/scar geometry may thwart termination in some cases. Across all models, the mechanism of GtACR1-mediated defibrillation was voltage forcing of illuminated tissue toward the modeled channel reversal potential of −40 mV, which made propagation through affected regions impossible. Thus, our findings suggest GtACR1-based optogenetic defibrillation of the human heart may be feasible with ≈2–3 orders of magnitude less energy than ChR2.
Collapse
Affiliation(s)
- Alexander R Ochs
- Department of Bioengineering, University of Washington, Seattle, WA, United States
| | - Thomas V Karathanos
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, United States
| | - Natalia A Trayanova
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, United States.,Alliance for Cardiovascular Diagnostic and Treatment Innovation, Johns Hopkins University, Baltimore, MD, United States
| | - Patrick M Boyle
- Department of Bioengineering, University of Washington, Seattle, WA, United States.,Institute for Stem Cell and Regenerative Medicine, University of Washington, Seattle, WA, United States.,Center for Cardiovascular Biology, University of Washington, Seattle, WA, United States
| |
Collapse
|
25
|
Tate JD, Elhabian S, Zemzemi N, Good WW, van Dam P, Brooks DH, MacLeod RS. A Cardiac Shape Model for Segmentation Uncertainty Quantification. COMPUTING IN CARDIOLOGY 2021; 48:10.23919/cinc53138.2021.9662917. [PMID: 35479610 PMCID: PMC9039803 DOI: 10.23919/cinc53138.2021.9662917] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Segmentation of cardiac images is a variable component of many patient specific computational pipelines, yet its impact on simulated results are still not fully understood. A hurdle to to exploring the impact of the segmentation variability is the technical challenge of building a statistical shape model of the ventricles. In this study, we improved open our previous shape analysis by creating a unified shape model including both the epicardium and endocardium. We tested four techniques within ShapeWorks to generate a ventricular shape model: standard, multidomain, hybrid multidomain, and geodesic distance. The multidomain and hybrid multidomain generated a shape model using all eleven segmentations, and the geodesic distance method generated a shape model using a subset of four segmentations. Each of the shape models captured spatially dependent characteristics of the segmentation variability, including wall thickness, annular diameter, and basal truncation. While each of the three methods have benefits, the hybrid multidomain approach provided the most accurate shape model with fewest points and may be most useful in a majority of applications.
Collapse
|
26
|
Plank G, Loewe A, Neic A, Augustin C, Huang YL, Gsell MAF, Karabelas E, Nothstein M, Prassl AJ, Sánchez J, Seemann G, Vigmond EJ. The openCARP simulation environment for cardiac electrophysiology. COMPUTER METHODS AND PROGRAMS IN BIOMEDICINE 2021; 208:106223. [PMID: 34171774 DOI: 10.1016/j.cmpb.2021.106223] [Citation(s) in RCA: 61] [Impact Index Per Article: 20.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/15/2021] [Accepted: 05/28/2021] [Indexed: 06/13/2023]
Abstract
BACKGROUND AND OBJECTIVE Cardiac electrophysiology is a medical specialty with a long and rich tradition of computational modeling. Nevertheless, no community standard for cardiac electrophysiology simulation software has evolved yet. Here, we present the openCARP simulation environment as one solution that could foster the needs of large parts of this community. METHODS AND RESULTS openCARP and the Python-based carputils framework allow developing and sharing simulation pipelines which automate in silico experiments including all modeling and simulation steps to increase reproducibility and productivity. The continuously expanding openCARP user community is supported by tailored infrastructure. Documentation and training material facilitate access to this complementary research tool for new users. After a brief historic review, this paper summarizes requirements for a high-usability electrophysiology simulator and describes how openCARP fulfills them. We introduce the openCARP modeling workflow in a multi-scale example of atrial fibrillation simulations on single cell, tissue, organ and body level and finally outline future development potential. CONCLUSION As an open simulator, openCARP can advance the computational cardiac electrophysiology field by making state-of-the-art simulations accessible. In combination with the carputils framework, it offers a tailored software solution for the scientific community and contributes towards increasing use, transparency, standardization and reproducibility of in silico experiments.
Collapse
Affiliation(s)
- Gernot Plank
- Gottfried Schatz Research Center, Division of Biophysics, Medical University of Graz, Graz, Austria.
| | - Axel Loewe
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | | | - Christoph Augustin
- Gottfried Schatz Research Center, Division of Biophysics, Medical University of Graz, Graz, Austria
| | - Yung-Lin Huang
- Institute for Experimental Cardiovascular Medicine, University Heart Center Freiburg. Bad Krozingen, Medical Center - University of Freiburg, Freiburg, Germany; Faculty of Medicine, University of Freiburg, Freiburg, Germany
| | - Matthias A F Gsell
- Gottfried Schatz Research Center, Division of Biophysics, Medical University of Graz, Graz, Austria
| | - Elias Karabelas
- Institute of Mathematics and Scientific Computing, University of Graz, Graz, Austria
| | - Mark Nothstein
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | - Anton J Prassl
- Gottfried Schatz Research Center, Division of Biophysics, Medical University of Graz, Graz, Austria
| | - Jorge Sánchez
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | - Gunnar Seemann
- Institute for Experimental Cardiovascular Medicine, University Heart Center Freiburg. Bad Krozingen, Medical Center - University of Freiburg, Freiburg, Germany; Faculty of Medicine, University of Freiburg, Freiburg, Germany
| | - Edward J Vigmond
- IHU Liryc, Electrophysiology and Heart Modeling Institute, Fondation Bordeaux Université, F-33600 Pessac-Bordeaux, France; Université Bordeaux, IMB, UMR 5251, F-33400 Talence, France
| |
Collapse
|
27
|
Kaboudian A, Cherry EM, Fenton FH. Real-Time Interactive Simulations of Complex Ionic Cardiac Cell Models in 2D and 3D Heart Structures with GPUs on Personal Computers. COMPUTING IN CARDIOLOGY 2021; 48:10.23919/cinc53138.2021.9662759. [PMID: 35754517 PMCID: PMC9228612 DOI: 10.23919/cinc53138.2021.9662759] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
AIMS Cardiac modeling in heart structures for the study of arrhythmia mechanisms requires the use of software that runs on supercomputers. Therefore, computational studies are limited to groups with access to computer clusters and personnel with high-performance computing experience. We present how to use and implement WebGL programs via a custom-written library to run and visualize simulations of the most complex ionic models in 2D and 3D, in real time, interactively using the multi-core GPU of a single computer. METHODS We use Abubu.js, a library we developed for solving partial differential equations such as those describing crystal growth and fluid flow, along with a newly implemented visualization algorithm, to simulate complex ionic cell models. By combining this library with JavaScript, we allow direct real-time interactions with simulations. We implemented: 1) modification of any model parameters and equations at any time, with direct access to the code while it runs, 2) electrode stimulation anywhere in the 2D/3D tissue with a mouse click, 3) saving the solution of the system at any time to re-initiate the dynamics from saved initial conditions, and 4) rotation/visualization of 3D structures at any angle. RESULTS As examples of this modeling platform, we implemented a phenomenological cell model and the human ventricular OVVR model (41 variables). In 2D we illustrate the dynamics in an annulus, disk, and square tissue; in 2D and 3D porcine ventricles, we show the initiation of functional/anatomical reentry, spiral wave dynamics in different regimes, initiation of early afterdepolarizations (EADs), and the effects of model parameter variations in real time. CONCLUSIONS We present the first simulations of complex models in anatomical structures with enhanced visualization and extended interactivity that run on a single PC, without software downloads, and as fast as in real-time even for 3D full ventricles.
Collapse
Affiliation(s)
- Abouzar Kaboudian
- School of Physics, Georgia Institute of Technology, Atlanta, GA, USA
| | - Elizabeth M Cherry
- School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, GA, USA
| | - Flavio H Fenton
- School of Physics, Georgia Institute of Technology, Atlanta, GA, USA
| |
Collapse
|
28
|
Vigmond EJ, Bouyssier J, Bayer J, Haïssaguerre M, Ashikaga H. On the nature of delays allowing anatomical re-entry involving the Purkinje network: a simulation study. Europace 2021; 23:i71-i79. [PMID: 33463686 DOI: 10.1093/europace/euaa395] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2020] [Accepted: 12/03/2020] [Indexed: 12/31/2022] Open
Abstract
AIMS Clinical observations suggest that the Purkinje network can be part of anatomical re-entry circuits in monomorphic or polymorphic ventricular arrhythmias. However, significant conduction delay is needed to support anatomical re-entry given the high conduction velocity within the Purkinje network. METHODS AND RESULTS We investigated, in computer models, whether damage rendering the Purkinje network as either an active lesion with slow conduction or a passive lesion with no excitable ionic channel, could explain clinical observations. Active lesions had compromised sodium current and a severe reduction in gap junction coupling, while passive lesions remained coupled by gap junctions, but modelled the membrane as a fixed resistance. Both types of tissue could provide significant delays of over 100 ms. Electrograms consistent with those obtained clinically were reproduced. However, passive tissue could not support re-entry as electrotonic coupling across the delay effectively increased the proximal refractory period to an extremely long interval. Active tissue, conversely, could robustly maintain re-entry. CONCLUSION Formation of anatomical re-entry using the Purkinje network is possible through highly reduced gap junctional coupling leading to slowed conduction.
Collapse
Affiliation(s)
- Edward J Vigmond
- IHU Liryc, Electrophysiology and Heart Modeling Institute, Fondation Bordeaux Université, F-33600 Pessac-Bordeaux, France.,Univ. Bordeaux, IMB, UMR 5251, F-33400 Talence, France
| | - Julien Bouyssier
- IHU Liryc, Electrophysiology and Heart Modeling Institute, Fondation Bordeaux Université, F-33600 Pessac-Bordeaux, France.,Univ. Bordeaux, IMB, UMR 5251, F-33400 Talence, France
| | - Jason Bayer
- IHU Liryc, Electrophysiology and Heart Modeling Institute, Fondation Bordeaux Université, F-33600 Pessac-Bordeaux, France.,Univ. Bordeaux, IMB, UMR 5251, F-33400 Talence, France
| | - Michel Haïssaguerre
- IHU Liryc, Electrophysiology and Heart Modeling Institute, Fondation Bordeaux Université, F-33600 Pessac-Bordeaux, France.,Univ. Bordeaux, Centre de recherche Cardio-Thoracique de Bordeaux, U1045, F-33000 Bordeaux, France.,Bordeaux University Hospital (CHU), Electrophysiology and Ablation Unit, F-33600 Pessac, France
| | - Hiroshi Ashikaga
- IHU Liryc, Electrophysiology and Heart Modeling Institute, Fondation Bordeaux Université, F-33600 Pessac-Bordeaux, France.,Cardiac Arrhythmia Service, Johns Hopkins University School of Medicine, Baltimore, MD, USA
| |
Collapse
|
29
|
Salvador M, Fedele M, Africa PC, Sung E, Dede' L, Prakosa A, Chrispin J, Trayanova N, Quarteroni A. Electromechanical modeling of human ventricles with ischemic cardiomyopathy: numerical simulations in sinus rhythm and under arrhythmia. Comput Biol Med 2021; 136:104674. [PMID: 34340126 DOI: 10.1016/j.compbiomed.2021.104674] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/04/2021] [Revised: 07/15/2021] [Accepted: 07/19/2021] [Indexed: 12/16/2022]
Abstract
We developed a novel patient-specific computational model for the numerical simulation of ventricular electromechanics in patients with ischemic cardiomyopathy (ICM). This model reproduces the activity both in sinus rhythm (SR) and in ventricular tachycardia (VT). The presence of scars, grey zones and non-remodeled regions of the myocardium is accounted for by the introduction of a spatially heterogeneous coefficient in the 3D electromechanics model. This 3D electromechanics model is firstly coupled with a 2-element Windkessel afterload model to fit the pressure-volume (PV) loop of a patient-specific left ventricle (LV) with ICM in SR. Then, we employ the coupling with a 0D closed-loop circulation model to analyze a VT circuit over multiple heartbeats on the same LV. We highlight similarities and differences on the solutions obtained by the electrophysiology model and those of the electromechanics model, while considering different scenarios for the circulatory system. We observe that very different parametrizations of the circulation model induce the same hemodynamical considerations for the patient at hand. Specifically, we classify this VT as unstable. We conclude by stressing the importance of combining electrophysiological, mechanical and hemodynamical models to provide relevant clinical indicators in how arrhythmias evolve and can potentially lead to sudden cardiac death.
Collapse
Affiliation(s)
- Matteo Salvador
- MOX-Dipartimento di Matematica, Politecnico di Milano, Milan, Italy.
| | - Marco Fedele
- MOX-Dipartimento di Matematica, Politecnico di Milano, Milan, Italy
| | | | - Eric Sung
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, USA
| | - Luca Dede'
- MOX-Dipartimento di Matematica, Politecnico di Milano, Milan, Italy
| | - Adityo Prakosa
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, USA
| | | | - Natalia Trayanova
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, USA
| | - Alfio Quarteroni
- MOX-Dipartimento di Matematica, Politecnico di Milano, Milan, Italy; École Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
| |
Collapse
|
30
|
Sánchez J, Luongo G, Nothstein M, Unger LA, Saiz J, Trenor B, Luik A, Dössel O, Loewe A. Using Machine Learning to Characterize Atrial Fibrotic Substrate From Intracardiac Signals With a Hybrid in silico and in vivo Dataset. Front Physiol 2021; 12:699291. [PMID: 34290623 PMCID: PMC8287829 DOI: 10.3389/fphys.2021.699291] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2021] [Accepted: 06/08/2021] [Indexed: 11/15/2022] Open
Abstract
In patients with atrial fibrillation, intracardiac electrogram signal amplitude is known to decrease with increased structural tissue remodeling, referred to as fibrosis. In addition to the isolation of the pulmonary veins, fibrotic sites are considered a suitable target for catheter ablation. However, it remains an open challenge to find fibrotic areas and to differentiate their density and transmurality. This study aims to identify the volume fraction and transmurality of fibrosis in the atrial substrate. Simulated cardiac electrograms, combined with a generalized model of clinical noise, reproduce clinically measured signals. Our hybrid dataset approach combines in silico and clinical electrograms to train a decision tree classifier to characterize the fibrotic atrial substrate. This approach captures different in vivo dynamics of the electrical propagation reflected on healthy electrogram morphology and synergistically combines it with synthetic fibrotic electrograms from in silico experiments. The machine learning algorithm was tested on five patients and compared against clinical voltage maps as a proof of concept, distinguishing non-fibrotic from fibrotic tissue and characterizing the patient's fibrotic tissue in terms of density and transmurality. The proposed approach can be used to overcome a single voltage cut-off value to identify fibrotic tissue and guide ablation targeting fibrotic areas.
Collapse
Affiliation(s)
- Jorge Sánchez
- Institute of Biomedical Engineering, Karlsruhe Institute for Technology, Karlsruhe, Germany
- Centro de Investigación e Innovación en Bioingeniería (Ci2B), Universitàt Politècnica de València, Valencia, Spain
| | - Giorgio Luongo
- Institute of Biomedical Engineering, Karlsruhe Institute for Technology, Karlsruhe, Germany
| | - Mark Nothstein
- Institute of Biomedical Engineering, Karlsruhe Institute for Technology, Karlsruhe, Germany
| | - Laura A. Unger
- Institute of Biomedical Engineering, Karlsruhe Institute for Technology, Karlsruhe, Germany
| | - Javier Saiz
- Centro de Investigación e Innovación en Bioingeniería (Ci2B), Universitàt Politècnica de València, Valencia, Spain
| | - Beatriz Trenor
- Centro de Investigación e Innovación en Bioingeniería (Ci2B), Universitàt Politècnica de València, Valencia, Spain
| | - Armin Luik
- Medizinische Klinik IV, Städtisches Klinikum Karlsruhe, Karlsruhe, Germany
| | - Olaf Dössel
- Institute of Biomedical Engineering, Karlsruhe Institute for Technology, Karlsruhe, Germany
| | - Axel Loewe
- Institute of Biomedical Engineering, Karlsruhe Institute for Technology, Karlsruhe, Germany
| |
Collapse
|
31
|
Monaci S, Gillette K, Puyol-Antón E, Rajani R, Plank G, King A, Bishop M. Automated Localization of Focal Ventricular Tachycardia From Simulated Implanted Device Electrograms: A Combined Physics-AI Approach. Front Physiol 2021; 12:682446. [PMID: 34276403 PMCID: PMC8281305 DOI: 10.3389/fphys.2021.682446] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2021] [Accepted: 05/31/2021] [Indexed: 11/13/2022] Open
Abstract
Background: Focal ventricular tachycardia (VT) is a life-threating arrhythmia, responsible for high morbidity rates and sudden cardiac death (SCD). Radiofrequency ablation is the only curative therapy against incessant VT; however, its success is dependent on accurate localization of its source, which is highly invasive and time-consuming. Objective: The goal of our study is, as a proof of concept, to demonstrate the possibility of utilizing electrogram (EGM) recordings from cardiac implantable electronic devices (CIEDs). To achieve this, we utilize fast and accurate whole torso electrophysiological (EP) simulations in conjunction with convolutional neural networks (CNNs) to automate the localization of focal VTs using simulated EGMs. Materials and Methods: A highly detailed 3D torso model was used to simulate ∼4000 focal VTs, evenly distributed across the left ventricle (LV), utilizing a rapid reaction-eikonal environment. Solutions were subsequently combined with lead field computations on the torso to derive accurate electrocardiograms (ECGs) and EGM traces, which were used as inputs to CNNs to localize focal sources. We compared the localization performance of a previously developed CNN architecture (Cartesian probability-based) with our novel CNN algorithm utilizing universal ventricular coordinates (UVCs). Results: Implanted device EGMs successfully localized VT sources with localization error (8.74 mm) comparable to ECG-based localization (6.69 mm). Our novel UVC CNN architecture outperformed the existing Cartesian probability-based algorithm (errors = 4.06 mm and 8.07 mm for ECGs and EGMs, respectively). Overall, localization was relatively insensitive to noise and changes in body compositions; however, displacements in ECG electrodes and CIED leads caused performance to decrease (errors 16-25 mm). Conclusion: EGM recordings from implanted devices may be used to successfully, and robustly, localize focal VT sources, and aid ablation planning.
Collapse
Affiliation(s)
| | - Karli Gillette
- Division of Biophysics, Medical University of Graz, Graz, Austria
| | | | | | - Gernot Plank
- Division of Biophysics, Medical University of Graz, Graz, Austria
| | - Andrew King
- King’s College London, London, United Kingdom
| | | |
Collapse
|
32
|
Gaur N, Qi XY, Benoist D, Bernus O, Coronel R, Nattel S, Vigmond EJ. A computational model of pig ventricular cardiomyocyte electrophysiology and calcium handling: Translation from pig to human electrophysiology. PLoS Comput Biol 2021; 17:e1009137. [PMID: 34191797 PMCID: PMC8277015 DOI: 10.1371/journal.pcbi.1009137] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2021] [Revised: 07/13/2021] [Accepted: 06/01/2021] [Indexed: 12/01/2022] Open
Abstract
The pig is commonly used as an experimental model of human heart disease, including for the study of mechanisms of arrhythmia. However, there exist differences between human and porcine cellular electrophysiology: The pig action potential (AP) has a deeper phase-1 notch, a longer duration at 50% repolarization, and higher plateau potentials than human. Ionic differences underlying the AP include larger rapid delayed-rectifier and smaller inward-rectifier K+-currents (IKr and IK1 respectively) in humans. AP steady-state rate-dependence and restitution is steeper in pigs. Porcine Ca2+ transients can have two components, unlike human. Although a reliable computational model for human ventricular myocytes exists, one for pigs is lacking. This hampers translation from results obtained in pigs to human myocardium. Here, we developed a computational model of the pig ventricular cardiomyocyte AP using experimental datasets of the relevant ionic currents, Ca2+-handling, AP shape, AP duration restitution, and inducibility of triggered activity and alternans. To properly capture porcine Ca2+ transients, we introduced a two-step process with a faster release in the t-tubular region, followed by a slower diffusion-induced release from a non t-tubular subcellular region. The pig model behavior was compared with that of a human ventricular cardiomyocyte (O’Hara-Rudy) model. The pig, but not the human model, developed early afterdepolarizations (EADs) under block of IK1, while IKr block led to EADs in the human but not in the pig model. At fast rates (pacing cycle length = 400 ms), the human cell model was more susceptible to spontaneous Ca2+ release-mediated delayed afterdepolarizations (DADs) and triggered activity than pig. Fast pacing led to alternans in human but not pig. Developing species-specific models incorporating electrophysiology and Ca2+-handling provides a tool to aid translating antiarrhythmic and arrhythmogenic assessment from the bench to the clinic. The pig is an animal commonly used experimentally to study diseases of the heart, as well as investigate therapies to treat them, such as drugs. However, although similar, pigs differ from humans in certain aspects which may mean experimental results do not always directly translate between species. We propose a mathematical model of porcine electrophysiology which can serve as a tool to understand differences between the species and translate responses. Using new measurements along with values from literature, we built a computer model of porcine cardiac myocyte which replicated voltage and calcium behaviour over a range of pacing frequencies. The pig cell had a two-stage calcium release, unlike humans with a single stage. We predict that pigs and humans differ in the type of potassium current block that makes them most susceptible to cardiac arrhythmia. The model we developed can elucidate important differences between human and pig arrhythmia response.
Collapse
Affiliation(s)
- Namit Gaur
- IHU Liryc, Electrophysiology and Heart Modeling Institute, Fondation Bordeaux Université, Pessac- Bordeaux, France
- Univ. Bordeaux, IMB, UMR 5251, Talence, France
| | - Xiao-Yan Qi
- Montreal Heart Institute and Université de Montréal, Montreal, Canada
| | - David Benoist
- IHU Liryc, Electrophysiology and Heart Modeling Institute, Fondation Bordeaux Université, Pessac- Bordeaux, France
- Univ. Bordeaux, Inserm, CRCTB, U1045, Pessac, France
| | - Olivier Bernus
- IHU Liryc, Electrophysiology and Heart Modeling Institute, Fondation Bordeaux Université, Pessac- Bordeaux, France
- Univ. Bordeaux, Inserm, CRCTB, U1045, Pessac, France
| | - Ruben Coronel
- IHU Liryc, Electrophysiology and Heart Modeling Institute, Fondation Bordeaux Université, Pessac- Bordeaux, France
- Department of Experimental Cardiology, Academic Medical Center, Amsterdam, The Netherlands
| | - Stanley Nattel
- IHU Liryc, Electrophysiology and Heart Modeling Institute, Fondation Bordeaux Université, Pessac- Bordeaux, France
- Montreal Heart Institute and Université de Montréal, Montreal, Canada
- Institute of Pharmacology, West German Heart and Vascular Center, Faculty of Medicine, University Duisburg-Essen, Essen, Germany
| | - Edward J. Vigmond
- IHU Liryc, Electrophysiology and Heart Modeling Institute, Fondation Bordeaux Université, Pessac- Bordeaux, France
- Univ. Bordeaux, IMB, UMR 5251, Talence, France
- * E-mail:
| |
Collapse
|
33
|
Nothstein M, Luik A, Jadidi A, Sánchez J, Unger LA, Wülfers EM, Dössel O, Seemann G, Schmitt C, Loewe A. CVAR-Seg: An Automated Signal Segmentation Pipeline for Conduction Velocity and Amplitude Restitution. Front Physiol 2021; 12:673047. [PMID: 34108887 PMCID: PMC8181407 DOI: 10.3389/fphys.2021.673047] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2021] [Accepted: 04/30/2021] [Indexed: 11/13/2022] Open
Abstract
BACKGROUND Rate-varying S1S2 stimulation protocols can be used for restitution studies to characterize atrial substrate, ionic remodeling, and atrial fibrillation risk. Clinical restitution studies with numerous patients create large amounts of these data. Thus, an automated pipeline to evaluate clinically acquired S1S2 stimulation protocol data necessitates consistent, robust, reproducible, and precise evaluation of local activation times, electrogram amplitude, and conduction velocity. Here, we present the CVAR-Seg pipeline, developed focusing on three challenges: (i) No previous knowledge of the stimulation parameters is available, thus, arbitrary protocols are supported. (ii) The pipeline remains robust under different noise conditions. (iii) The pipeline supports segmentation of atrial activities in close temporal proximity to the stimulation artifact, which is challenging due to larger amplitude and slope of the stimulus compared to the atrial activity. METHODS AND RESULTS The S1 basic cycle length was estimated by time interval detection. Stimulation time windows were segmented by detecting synchronous peaks in different channels surpassing an amplitude threshold and identifying time intervals between detected stimuli. Elimination of the stimulation artifact by a matched filter allowed detection of local activation times in temporal proximity. A non-linear signal energy operator was used to segment periods of atrial activity. Geodesic and Euclidean inter electrode distances allowed approximation of conduction velocity. The automatic segmentation performance of the CVAR-Seg pipeline was evaluated on 37 synthetic datasets with decreasing signal-to-noise ratios. Noise was modeled by reconstructing the frequency spectrum of clinical noise. The pipeline retained a median local activation time error below a single sample (1 ms) for signal-to-noise ratios as low as 0 dB representing a high clinical noise level. As a proof of concept, the pipeline was tested on a CARTO case of a paroxysmal atrial fibrillation patient and yielded plausible restitution curves for conduction speed and amplitude. CONCLUSION The proposed openly available CVAR-Seg pipeline promises fast, fully automated, robust, and accurate evaluations of atrial signals even with low signal-to-noise ratios. This is achieved by solving the proximity problem of stimulation and atrial activity to enable standardized evaluation without introducing human bias for large data sets.
Collapse
Affiliation(s)
- Mark Nothstein
- Institute of Biomedical Engineering (IBT), Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | - Armin Luik
- Medizinische Klinik IV, Städtisches Klinikum Karlsruhe, Karlsruhe, Germany
| | - Amir Jadidi
- Klinik für Kardiologie und Angiologie II, University Heart Center Freiburg-Bad Krozingen, Bad Krozingen, Germany
- Faculty of Medicine, University of Freiburg, Freiburg, Germany
| | - Jorge Sánchez
- Institute of Biomedical Engineering (IBT), Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | - Laura A. Unger
- Institute of Biomedical Engineering (IBT), Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | - Eike M. Wülfers
- Faculty of Medicine, University of Freiburg, Freiburg, Germany
- Institute for Experimental Cardiovascular Medicine, University Heart Center Freiburg-Bad Krozingen, Freiburg, Germany
| | - Olaf Dössel
- Institute of Biomedical Engineering (IBT), Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | - Gunnar Seemann
- Faculty of Medicine, University of Freiburg, Freiburg, Germany
- Institute for Experimental Cardiovascular Medicine, University Heart Center Freiburg-Bad Krozingen, Freiburg, Germany
| | - Claus Schmitt
- Medizinische Klinik IV, Städtisches Klinikum Karlsruhe, Karlsruhe, Germany
| | - Axel Loewe
- Institute of Biomedical Engineering (IBT), Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| |
Collapse
|
34
|
Bifulco SF, Scott GD, Sarairah S, Birjandian Z, Roney CH, Niederer SA, Mahnkopf C, Kuhnlein P, Mitlacher M, Tirschwell D, Longstreth WT, Akoum N, Boyle PM. Computational modeling identifies embolic stroke of undetermined source patients with potential arrhythmic substrate. eLife 2021; 10:e64213. [PMID: 33942719 PMCID: PMC8143793 DOI: 10.7554/elife.64213] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/21/2020] [Accepted: 04/16/2021] [Indexed: 12/25/2022] Open
Abstract
Cardiac magnetic resonance imaging (MRI) has revealed fibrosis in embolic stroke of undetermined source (ESUS) patients comparable to levels seen in atrial fibrillation (AFib). We used computational modeling to understand the absence of arrhythmia in ESUS despite the presence of putatively pro-arrhythmic fibrosis. MRI-based atrial models were reconstructed for 45 ESUS and 45 AFib patients. The fibrotic substrate's arrhythmogenic capacity in each patient was assessed computationally. Reentrant drivers were induced in 24/45 (53%) ESUS and 22/45 (49%) AFib models. Inducible models had more fibrosis (16.7 ± 5.45%) than non-inducible models (11.07 ± 3.61%; p<0.0001); however, inducible subsets of ESUS and AFib models had similar fibrosis levels (p=0.90), meaning that the intrinsic pro-arrhythmic substrate properties of fibrosis in ESUS and AFib are indistinguishable. This suggests that some ESUS patients have latent pre-clinical fibrotic substrate that could be a future source of arrhythmogenicity. Thus, our work prompts the hypothesis that ESUS patients with fibrotic atria are spared from AFib due to an absence of arrhythmia triggers.
Collapse
Affiliation(s)
- Savannah F Bifulco
- Department of Bioengineering, University of WashingtonSeattleUnited States
| | - Griffin D Scott
- Department of Bioengineering, University of WashingtonSeattleUnited States
| | - Sakher Sarairah
- Division of Cardiology, University of WashingtonSeattleUnited States
| | - Zeinab Birjandian
- Division of Cardiology, University of WashingtonSeattleUnited States
- Department of Neurology, University of WashingtonSeattleUnited States
| | - Caroline H Roney
- School of Biomedical Engineering and Imaging Sciences, King’s College LondonLondonUnited Kingdom
| | - Steven A Niederer
- School of Biomedical Engineering and Imaging Sciences, King’s College LondonLondonUnited Kingdom
| | | | | | | | - David Tirschwell
- Department of Neurology, University of WashingtonSeattleUnited States
| | - WT Longstreth
- Department of Neurology, University of WashingtonSeattleUnited States
- Department of Epidemiology, University of WashingtonSeattleUnited States
| | - Nazem Akoum
- Division of Cardiology, University of WashingtonSeattleUnited States
| | - Patrick M Boyle
- Department of Bioengineering, University of WashingtonSeattleUnited States
- Center for Cardiovascular Biology, University of WashingtonSeattleUnited States
- Institute for Stem Cell and Regenerative Medicine, University of WashingtonSeattleUnited States
| |
Collapse
|
35
|
Boyle PM, Yu J, Klimas A, Williams JC, Trayanova NA, Entcheva E. OptoGap is an optogenetics-enabled assay for quantification of cell-cell coupling in multicellular cardiac tissue. Sci Rep 2021; 11:9310. [PMID: 33927252 PMCID: PMC8085001 DOI: 10.1038/s41598-021-88573-1] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2020] [Accepted: 03/31/2021] [Indexed: 12/23/2022] Open
Abstract
Intercellular electrical coupling is an essential means of communication between cells. It is important to obtain quantitative knowledge of such coupling between cardiomyocytes and non-excitable cells when, for example, pathological electrical coupling between myofibroblasts and cardiomyocytes yields increased arrhythmia risk or during the integration of donor (e.g., cardiac progenitor) cells with native cardiomyocytes in cell-therapy approaches. Currently, there is no direct method for assessing heterocellular coupling within multicellular tissue. Here we demonstrate experimentally and computationally a new contactless assay for electrical coupling, OptoGap, based on selective illumination of inexcitable cells that express optogenetic actuators and optical sensing of the response of coupled excitable cells (e.g., cardiomyocytes) that are light-insensitive. Cell-cell coupling is quantified by the energy required to elicit an action potential via junctional current from the light-stimulated cell(s). The proposed technique is experimentally validated against the standard indirect approach, GapFRAP, using light-sensitive cardiac fibroblasts and non-transformed cardiomyocytes in a two-dimensional setting. Its potential applicability to the complex three-dimensional setting of the native heart is corroborated by computational modelling and proper calibration. Lastly, the sensitivity of OptoGap to intrinsic cell-scale excitability is robustly characterized via computational analysis.
Collapse
Affiliation(s)
- Patrick M Boyle
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, USA
- Institute for Computational Medicine, Johns Hopkins University, Baltimore, MD, USA
- Department of Bioengineering, University of Washington, Seattle, WA, USA
- Institute for Stem Cell and Regenerative Medicine, University of Washington, Seattle, WA, USA
- Center for Cardiovascular Biology, University of Washington, Seattle, WA, USA
| | - Jinzhu Yu
- Department of Biomedical Engineering, Stony Brook University, Stony Brook, NY, USA
| | - Aleksandra Klimas
- Department of Biomedical Engineering, Stony Brook University, Stony Brook, NY, USA
- Department of Biomedical Engineering, George Washington University, 800 22nd Street NW, Suite 5000, Washington, DC, 20052, USA
| | - John C Williams
- Department of Biomedical Engineering, Stony Brook University, Stony Brook, NY, USA
| | - Natalia A Trayanova
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, USA
- Institute for Computational Medicine, Johns Hopkins University, Baltimore, MD, USA
- Alliance for Cardiovascular Diagnostic and Treatment Innovation, Johns Hopkins University, Baltimore, MD, USA
| | - Emilia Entcheva
- Department of Biomedical Engineering, Stony Brook University, Stony Brook, NY, USA.
- Department of Biomedical Engineering, George Washington University, 800 22nd Street NW, Suite 5000, Washington, DC, 20052, USA.
| |
Collapse
|
36
|
Yu JK, Liang JA, Franceschi WH, Huang Q, Pashakhanloo F, Sung E, Boyle PM, Trayanova NA. Assessment of arrhythmia mechanism and burden of the infarcted ventricles following remuscularization with pluripotent stem cell-derived cardiomyocyte patches using patient-derived models. Cardiovasc Res 2021; 118:1247-1261. [PMID: 33881518 DOI: 10.1093/cvr/cvab140] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 07/06/2020] [Revised: 01/14/2021] [Accepted: 04/19/2021] [Indexed: 12/24/2022] Open
Abstract
AIMS Direct remuscularization with pluripotent stem cell-derived cardiomyocytes (PSC-CMs) seeks to address the onset of heart failure post-myocardial infarction (MI) by treating the persistent muscle deficiency that underlies it. However, direct remuscularization with PSC-CMs could potentially be arrhythmogenic. We investigated two possible mechanisms of arrhythmogenesis-focal vs reentrant-arising from direct remuscularization with PSC-CM patches in two personalized, human ventricular computer models of post-MI. Moreover, we developed a principled approach for evaluating arrhythmogenicity of direct remuscularization that factors in the VT propensity of the patient-specific post-MI fibrotic substrate and use it to investigate different conditions of patch remuscularization. METHODS & RESULTS Two personalized, human ventricular models of post-MI (P1 & P2) were constructed from late gadolinium enhanced (LGE)-magnetic resonance images (MRI). In each model, remuscularization with PSC-CM patches were simulated under different treatment conditions that included patch engraftment, patch myofibril orientation, remuscularization site, patch size (thickness and diameter), and patch maturation. To determine arrhythmogenicity of treatment conditions, VT burden of heart models was quantified prior to and after simulated remuscularization and compared. VT burden was quantified based on inducibility (i.e., weighted sum of pacing sites that induced) and severity (i.e., the number of distinct VT morphologies induced). Prior to remuscularization, VT burden was significant in P1 (0.275) and not in P2 (0.0, not VT inducible). We highlight that reentrant VT mechanisms would dominate over focal mechanisms; spontaneous beats emerging from PSC-CM grafts were always a fraction of resting sinus rate. Moreover, incomplete patch engraftment can be particularly arrhythmogenic, giving rise to particularly aberrant electrical activation and conduction slowing across the PSC-CM patches along with elevated VT burden when compared to complete engraftment. Under conditions of complete patch engraftment, remuscularization was almost always arrhythmogenic in P2 but certain treatment conditions could be anti-arrhythmogenic in P1. Moreover, the remuscularization site was the most important factor affecting VT burden in both P1 and P2. Complete maturation of PSC-CM patches, both ionically and electrotonically, at the appropriate site could completely alleviate VT burden. CONCLUSION We identified that reentrant VT would be the primary VT mechanism in patch remuscularization. To evaluate the arrhythmogenicity of remuscularization, we developed a principled approach that factors in the propensity of the patient-specific fibrotic substrate for VT. We showed that arrhythmogenicity is sensitive to the patient-specific fibrotic substrate and remuscularization site. We demonstrate that targeted remuscularization can be safe in the appropriate individual and holds the potential to nondestructively eliminate VT post-MI in addition to addressing muscle deficiency underlying heart failure progression. TRANSLATIONAL PERSPECTIVE If safety from ventricular arrhythmias can be addressed, direct remuscularization with PSC-CMs-achieved either through engineered myocardial patches or intramyocardial injections-holds the potential to halt heart failure progression post-MI. Using personalized 3 D models of the post-MI ventricles derived from LGE-MRI, we provide evidence that arrhythmogenesis following remuscularization with PSC-CM patches is driven by a reentrant as opposed to focal VT mechanism. Moreover, the existing patient-specific fibrotic substrate together with the remuscularization site were primary determinants of arrhythmogenesis. These results suggest that the clinical safety of remuscularization can be achieved through patient-specific optimization guided in-part by computational modeling.
Collapse
Affiliation(s)
- Joseph K Yu
- Institute for Computational Medicine, Johns Hopkins University, 3400 N Charles Street, 208 Hackerman, Baltimore, MD, 21218, USA.,Department of Biomedical Engineering, Johns Hopkins University, 3400 N Charles Street, 208 Hackerman, Baltimore, MD, 21218, USA.,Department of Medicine, Johns Hopkins University School of Medicine, Baltimore, MD, USA.,Alliance for Cardiovascular Diagnostic and Treatment Innovation (ADVANCE), Johns Hopkins University, 3400 N Charles Street, 216 Hackerman, Baltimore, MD, USA
| | - Jialiu A Liang
- Institute for Computational Medicine, Johns Hopkins University, 3400 N Charles Street, 208 Hackerman, Baltimore, MD, 21218, USA.,Department of Biomedical Engineering, Johns Hopkins University, 3400 N Charles Street, 208 Hackerman, Baltimore, MD, 21218, USA
| | - William H Franceschi
- Institute for Computational Medicine, Johns Hopkins University, 3400 N Charles Street, 208 Hackerman, Baltimore, MD, 21218, USA.,Department of Biomedical Engineering, Johns Hopkins University, 3400 N Charles Street, 208 Hackerman, Baltimore, MD, 21218, USA
| | - Qinwen Huang
- Institute for Computational Medicine, Johns Hopkins University, 3400 N Charles Street, 208 Hackerman, Baltimore, MD, 21218, USA.,Department of Biomedical Engineering, Johns Hopkins University, 3400 N Charles Street, 208 Hackerman, Baltimore, MD, 21218, USA
| | - Farhad Pashakhanloo
- Department of Biomedical Engineering, Johns Hopkins University, 3400 N Charles Street, 208 Hackerman, Baltimore, MD, 21218, USA
| | - Eric Sung
- Department of Biomedical Engineering, Johns Hopkins University, 3400 N Charles Street, 208 Hackerman, Baltimore, MD, 21218, USA.,Department of Medicine, Johns Hopkins University School of Medicine, Baltimore, MD, USA.,Alliance for Cardiovascular Diagnostic and Treatment Innovation (ADVANCE), Johns Hopkins University, 3400 N Charles Street, 216 Hackerman, Baltimore, MD, USA
| | - Patrick M Boyle
- Department of Biomedical Engineering, Johns Hopkins University, 3400 N Charles Street, 208 Hackerman, Baltimore, MD, 21218, USA.,Department of Bioengineering, University of Washington, Seattle, WA, USA.,Institute for Stem Cell and Regenerative Medicine, University of Washington, Seattle, WA, USA.,Center for Cardiovascular Biology, University of Washington, Seattle, WA, USA
| | - Natalia A Trayanova
- Department of Biomedical Engineering, Johns Hopkins University, 3400 N Charles Street, 208 Hackerman, Baltimore, MD, 21218, USA.,Department of Medicine, Johns Hopkins University School of Medicine, Baltimore, MD, USA.,Alliance for Cardiovascular Diagnostic and Treatment Innovation (ADVANCE), Johns Hopkins University, 3400 N Charles Street, 216 Hackerman, Baltimore, MD, USA
| |
Collapse
|
37
|
Vigmond EJ, Neic A, Blauer J, Swenson D, Plank G. How Electrode Position Affects Selective His Bundle Capture: A Modelling Study. IEEE Trans Biomed Eng 2021; 68:3410-3416. [PMID: 33835914 DOI: 10.1109/tbme.2021.3072334] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
In certain cardiac conduction system pathologies, like bundle branch block, block may be proximal, allowing for electrical stimulation of the more distal His bundle to most effectively restore activation. While selective stimulation of the His bundle is sought, surrounding myocardium may also be excited, resulting in nonselective pacing. The myocardium and His bundle have distinct capture thresholds, but the factors affecting whether His bundle pacing is selective or nonselective remain unelucidated. OBJECTIVE We investigated the properties which affect the capture thresholds in order to improve selective pacing. METHODS We performed biophysically detailed, computer simulations of a His fibre running through a septal wedge preparation to compute capture thresholds under various configurations of electrode polarity and orientation. RESULTS The myocardial capture threshold was close to that of the His bundle. The His fibre needed to intersect with the electrode tip to favor its activation. Inserting the electrode fully within the septum increased the myocardial capture threshold. Reversing polarity, to rely on anode break excitation, also increased the ease of selective pacing. CONCLUSION Model results were consistent with clinical observations. For selective pacing, the tip needs to be in contact with the His fibre and anodal stimulation is preferable. SIGNIFICANCE This study provides insight into helping establish electrode and stimulation parameters for selective His bundle pacing in patients.
Collapse
|
38
|
Azzolin L, Schuler S, Dössel O, Loewe A. A Reproducible Protocol to Assess Arrhythmia Vulnerability in silico: Pacing at the End of the Effective Refractory Period. Front Physiol 2021; 12:656411. [PMID: 33868025 PMCID: PMC8047415 DOI: 10.3389/fphys.2021.656411] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/20/2021] [Accepted: 03/10/2021] [Indexed: 01/31/2023] Open
Abstract
In both clinical and computational studies, different pacing protocols are used to induce arrhythmia and non-inducibility is often considered as the endpoint of treatment. The need for a standardized methodology is urgent since the choice of the protocol used to induce arrhythmia could lead to contrasting results, e.g., in assessing atrial fibrillation (AF) vulnerabilty. Therefore, we propose a novel method-pacing at the end of the effective refractory period (PEERP)-and compare it to state-of-the-art protocols, such as phase singularity distribution (PSD) and rapid pacing (RP) in a computational study. All methods were tested by pacing from evenly distributed endocardial points at 1 cm inter-point distance in two bi-atrial geometries. Seven different atrial models were implemented: five cases without specific AF-induced remodeling but with decreasing global conduction velocity and two persistent AF cases with an increasing amount of fibrosis resembling different substrate remodeling stages. Compared with PSD and RP, PEERP induced a larger variety of arrhythmia complexity requiring, on average, only 2.7 extra-stimuli and 3 s of simulation time to initiate reentry. Moreover, PEERP and PSD were the protocols which unveiled a larger number of areas vulnerable to sustain stable long living reentries compared to RP. Finally, PEERP can foster standardization and reproducibility, since, in contrast to the other protocols, it is a parameter-free method. Furthermore, we discuss its clinical applicability. We conclude that the choice of the inducing protocol has an influence on both initiation and maintenance of AF and we propose and provide PEERP as a reproducible method to assess arrhythmia vulnerability.
Collapse
Affiliation(s)
- Luca Azzolin
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | | | | | | |
Collapse
|
39
|
Luongo G, Azzolin L, Schuler S, Rivolta MW, Almeida TP, Martínez JP, Soriano DC, Luik A, Müller-Edenborn B, Jadidi A, Dössel O, Sassi R, Laguna P, Loewe A. Machine learning enables noninvasive prediction of atrial fibrillation driver location and acute pulmonary vein ablation success using the 12-lead ECG. CARDIOVASCULAR DIGITAL HEALTH JOURNAL 2021; 2:126-136. [PMID: 33899043 PMCID: PMC8053175 DOI: 10.1016/j.cvdhj.2021.03.002] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022] Open
Abstract
BACKGROUND Atrial fibrillation (AF) is the most common supraventricular arrhythmia, characterized by disorganized atrial electrical activity, maintained by localized arrhythmogenic atrial drivers. Pulmonary vein isolation (PVI) allows to exclude PV-related drivers. However, PVI is less effective in patients with additional extra-PV arrhythmogenic drivers. OBJECTIVES To discriminate whether AF drivers are located near the PVs vs extra-PV regions using the noninvasive 12-lead electrocardiogram (ECG) in a computational and clinical framework, and to computationally predict the acute success of PVI in these cohorts of data. METHODS AF drivers were induced in 2 computerized atrial models and combined with 8 torso models, resulting in 1128 12-lead ECGs (80 ECGs with AF drivers located in the PVs and 1048 in extra-PV areas). A total of 103 features were extracted from the signals. Binary decision tree classifier was trained on the simulated data and evaluated using hold-out cross-validation. The PVs were subsequently isolated in the models to assess PVI success. Finally, the classifier was tested on a clinical dataset (46 patients: 23 PV-dependent AF and 23 with additional extra-PV sources). RESULTS The classifier yielded 82.6% specificity and 73.9% sensitivity for detecting PV drivers on the clinical data. Consistency analysis on the 46 patients resulted in 93.5% results match. Applying PVI on the simulated AF cases terminated AF in 100% of the cases in the PV class. CONCLUSION Machine learning-based classification of 12-lead-ECG allows discrimination between patients with PV drivers vs those with extra-PV drivers of AF. The novel algorithm may aid to identify patients with high acute success rates to PVI.
Collapse
Affiliation(s)
- Giorgio Luongo
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology, Karlsruhe, Germany
- Address reprint requests and correspondence: Mr Giorgio Luongo, Fritz-Haber-Weg 1, 76131 Karlsruhe, Germany.
| | - Luca Azzolin
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology, Karlsruhe, Germany
| | - Steffen Schuler
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology, Karlsruhe, Germany
| | - Massimo W. Rivolta
- Dipartimento di Informatica, Università degli Studi di Milano, Milan, Italy
| | - Tiago P. Almeida
- Department of Cardiovascular Sciences, University of Leicester, Leicester, United Kingdom
| | | | - Diogo C. Soriano
- Engineering, Modelling and Applied Social Sciences Centre, ABC Federal University, São Bernardo do Campo, Brazil
| | - Armin Luik
- Medizinische Klinik IV, Städtisches Klinikum Karlsruhe, Karlsruhe, Germany
| | - Björn Müller-Edenborn
- Department of Electrophysiology, University-Heart-Center Freiburg-Bad Krozingen, Bad Krozingen Campus, Bad Krozingen, Germany
| | - Amir Jadidi
- Department of Electrophysiology, University-Heart-Center Freiburg-Bad Krozingen, Bad Krozingen Campus, Bad Krozingen, Germany
| | - Olaf Dössel
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology, Karlsruhe, Germany
| | - Roberto Sassi
- Dipartimento di Informatica, Università degli Studi di Milano, Milan, Italy
| | - Pablo Laguna
- I3A, Universidad de Zaragoza, and CIBER-BNN, Zaragoza, Spain
| | - Axel Loewe
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology, Karlsruhe, Germany
| |
Collapse
|
40
|
Rodero C, Strocchi M, Marciniak M, Longobardi S, Whitaker J, O’Neill MD, Gillette K, Augustin C, Plank G, Vigmond EJ, Lamata P, Niederer SA. Linking statistical shape models and simulated function in the healthy adult human heart. PLoS Comput Biol 2021; 17:e1008851. [PMID: 33857152 PMCID: PMC8049237 DOI: 10.1371/journal.pcbi.1008851] [Citation(s) in RCA: 24] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2020] [Accepted: 03/03/2021] [Indexed: 01/09/2023] Open
Abstract
Cardiac anatomy plays a crucial role in determining cardiac function. However, there is a poor understanding of how specific and localised anatomical changes affect different cardiac functional outputs. In this work, we test the hypothesis that in a statistical shape model (SSM), the modes that are most relevant for describing anatomy are also most important for determining the output of cardiac electromechanics simulations. We made patient-specific four-chamber heart meshes (n = 20) from cardiac CT images in asymptomatic subjects and created a SSM from 19 cases. Nine modes captured 90% of the anatomical variation in the SSM. Functional simulation outputs correlated best with modes 2, 3 and 9 on average (R = 0.49 ± 0.17, 0.37 ± 0.23 and 0.34 ± 0.17 respectively). We performed a global sensitivity analysis to identify the different modes responsible for different simulated electrical and mechanical measures of cardiac function. Modes 2 and 9 were the most important for determining simulated left ventricular mechanics and pressure-derived phenotypes. Mode 2 explained 28.56 ± 16.48% and 25.5 ± 20.85, and mode 9 explained 12.1 ± 8.74% and 13.54 ± 16.91% of the variances of mechanics and pressure-derived phenotypes, respectively. Electrophysiological biomarkers were explained by the interaction of 3 ± 1 modes. In the healthy adult human heart, shape modes that explain large portions of anatomical variance do not explain equivalent levels of electromechanical functional variation. As a result, in cardiac models, representing patient anatomy using a limited number of modes of anatomical variation can cause a loss in accuracy of simulated electromechanical function.
Collapse
Affiliation(s)
- Cristobal Rodero
- Cardiac Electromechanics Research Group, Biomedical Engineering Department, King´s College London, London, United Kingdom
- Cardiac Modelling and Imaging Biomarkers, Biomedical Engineering Department, King´s College London, London, United Kingdom
- * E-mail:
| | - Marina Strocchi
- Cardiac Electromechanics Research Group, Biomedical Engineering Department, King´s College London, London, United Kingdom
| | - Maciej Marciniak
- Cardiac Modelling and Imaging Biomarkers, Biomedical Engineering Department, King´s College London, London, United Kingdom
| | - Stefano Longobardi
- Cardiac Electromechanics Research Group, Biomedical Engineering Department, King´s College London, London, United Kingdom
| | - John Whitaker
- Cardiovascular Imaging Department, King’s College London, London, United Kingdom
| | - Mark D. O’Neill
- Department of Cardiology, St Thomas’ Hospital, London, United Kingdom
| | - Karli Gillette
- Institute of Biophysics, Medical University of Graz, Graz, Austria
| | | | - Gernot Plank
- Institute of Biophysics, Medical University of Graz, Graz, Austria
| | - Edward J. Vigmond
- Institute of Electrophysiology and Heart Modeling, Foundation Bordeaux University, Bordeaux, France
- Bordeaux Institute of Mathematics, University of Bordeaux, Bordeaux, France
| | - Pablo Lamata
- Cardiac Modelling and Imaging Biomarkers, Biomedical Engineering Department, King´s College London, London, United Kingdom
| | - Steven A. Niederer
- Cardiac Electromechanics Research Group, Biomedical Engineering Department, King´s College London, London, United Kingdom
| |
Collapse
|
41
|
Balaban G, Halliday BP, Porter B, Bai W, Nygåard S, Owen R, Hatipoglu S, Ferreira ND, Izgi C, Tayal U, Corden B, Ware J, Pennell DJ, Rueckert D, Plank G, Rinaldi CA, Prasad SK, Bishop MJ. Late-Gadolinium Enhancement Interface Area and Electrophysiological Simulations Predict Arrhythmic Events in Patients With Nonischemic Dilated Cardiomyopathy. JACC Clin Electrophysiol 2021; 7:238-249. [PMID: 33602406 PMCID: PMC7900608 DOI: 10.1016/j.jacep.2020.08.036] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2020] [Revised: 08/18/2020] [Accepted: 08/19/2020] [Indexed: 11/16/2022]
Abstract
OBJECTIVES This study sought to investigate whether shape-based late gadolinium enhancement (LGE) metrics and simulations of re-entrant electrical activity are associated with arrhythmic events in patients with nonischemic dilated cardiomyopathy (NIDCM). BACKGROUND The presence of LGE predicts life-threatening ventricular arrhythmias in NIDCM; however, risk stratification remains imprecise. LGE shape and simulations of electrical activity may be able to provide additional prognostic information. METHODS Cardiac magnetic resonance (CMR)-LGE shape metrics were computed for a cohort of 156 patients with NIDCM and visible LGE and tested retrospectively for an association with an arrhythmic composite endpoint of sudden cardiac death and ventricular tachycardia. Computational models were created from images and used in conjunction with simulated stimulation protocols to assess the potential for re-entry induction in each patient's scar morphology. A mechanistic analysis of the simulations was carried out to explain the associations. RESULTS During a median follow-up of 1,611 (interquartile range: 881 to 2,341) days, 16 patients (10.3%) met the primary endpoint. In an inverse probability weighted Cox regression, the LGE-myocardial interface area (hazard ratio [HR]: 1.75; 95% confidence interval [CI]: 1.24 to 2.47; p = 0.001), number of simulated re-entries (HR: 1.40; 95% CI: 1.23 to 1.59; p < 0.01) and LGE volume (HR: 1.44; 95% CI: 1.07 to 1.94; p = 0.02) were associated with arrhythmic events. Computational modeling revealed repolarization heterogeneity and rate-dependent block of electrical wavefronts at the LGE-myocardial interface as putative arrhythmogenic mechanisms directly related to the LGE interface area. CONCLUSIONS The area of interface between scar and surviving myocardium, as well as simulated re-entrant activity, are associated with an elevated risk of major arrhythmic events in patients with NIDCM and LGE and represent novel risk predictors.
Collapse
Affiliation(s)
- Gabriel Balaban
- Department of Biomedical Engineering, School of Biomedical & Imaging Sciences, King's College London, United Kingdom; Department of Informatics, University of Oslo, Oslo, Norway
| | - Brian P Halliday
- Cardiovascular Magnetic Resonance Unit, Royal Brompton Hospital, London, United Kingdom; National Heart and Lung Institute, Imperial College London, London, United Kingdom
| | - Bradley Porter
- Department of Biomedical Engineering, School of Biomedical & Imaging Sciences, King's College London, United Kingdom; Department of Cardiology, St Thomas' Hospital, London, United Kingdom
| | - Wenjia Bai
- Department of Computer Science, Imperial College London, United Kingdom
| | - Ståle Nygåard
- Department of Informatics, University of Oslo, Oslo, Norway
| | - Ruth Owen
- London School of Hygiene and Tropical Medicine, London, United Kingdom
| | - Suzan Hatipoglu
- Cardiovascular Magnetic Resonance Unit, Royal Brompton Hospital, London, United Kingdom
| | - Nuno Dias Ferreira
- Cardiovascular Magnetic Resonance Unit, Royal Brompton Hospital, London, United Kingdom
| | - Cemil Izgi
- Cardiovascular Magnetic Resonance Unit, Royal Brompton Hospital, London, United Kingdom
| | - Upasana Tayal
- Cardiovascular Magnetic Resonance Unit, Royal Brompton Hospital, London, United Kingdom
| | - Ben Corden
- National Heart and Lung Institute, Imperial College London, London, United Kingdom
| | - James Ware
- National Heart and Lung Institute, Imperial College London, London, United Kingdom
| | - Dudley J Pennell
- Cardiovascular Magnetic Resonance Unit, Royal Brompton Hospital, London, United Kingdom; National Heart and Lung Institute, Imperial College London, London, United Kingdom
| | - Daniel Rueckert
- Department of Computer Science, Imperial College London, United Kingdom
| | - Gernot Plank
- Institute of Biophysics, Medical University of Graz, Graz, Austria
| | - Christopher A Rinaldi
- Department of Biomedical Engineering, School of Biomedical & Imaging Sciences, King's College London, United Kingdom; Department of Cardiology, St Thomas' Hospital, London, United Kingdom
| | - Sanjay K Prasad
- Cardiovascular Magnetic Resonance Unit, Royal Brompton Hospital, London, United Kingdom; National Heart and Lung Institute, Imperial College London, London, United Kingdom
| | - Martin J Bishop
- Department of Biomedical Engineering, School of Biomedical & Imaging Sciences, King's College London, United Kingdom.
| |
Collapse
|
42
|
Building Models of Patient-Specific Anatomy and Scar Morphology from Clinical MRI Data. SYSTEMS MEDICINE 2021. [DOI: 10.1016/b978-0-12-801238-3.11663-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022] Open
|
43
|
Albatat M, Arevalo H, Bergsland J, Strøm V, Balasingham I, Odland HH. Optimal pacing sites in cardiac resynchronization by left ventricular activation front analysis. Comput Biol Med 2020; 128:104159. [PMID: 33301952 DOI: 10.1016/j.compbiomed.2020.104159] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2020] [Revised: 11/14/2020] [Accepted: 11/29/2020] [Indexed: 10/22/2022]
Abstract
Cardiac resynchronization therapy (CRT) can substantially improve dyssynchronous heart failure and reduce mortality. However, about one-third of patients who are implanted, derive no measurable benefit from CRT. Non-response may partly be due to suboptimal activation of the left ventricle (LV) caused by electrophysiological heterogeneities. The goal of this study is to investigate the performance of a newly developed method used to analyze electrical wavefront propagation in a heart model including myocardial scar and compare this to clinical benchmark studies. We used computational models to measure the maximum activation front (MAF) in the LV during different pacing scenarios. Different heart geometries and scars were created based on cardiac MR images of three patients. The right ventricle (RV) was paced from the apex and the LV was paced from 12 different sites, single site, dual-site and triple site. Our results showed that for single LV site pacing, the pacing site with the largest MAF corresponded with the latest activated regions of the LV demonstrated during RV pacing, which also agrees with previous markers used for predicting optimal single-site pacing location. We then demonstrated the utility of MAF in predicting optimal electrode placements in more complex scenarios including scar and multi-site LV pacing. This study demonstrates the potential value of computational simulations in understanding and planning CRT.
Collapse
Affiliation(s)
- Mohammad Albatat
- Intervention Centre, Oslo University Hospital, Oslo, Norway; Institute of Clinical Medicine, University of Oslo, Oslo, Norway.
| | - Hermenegild Arevalo
- Department of Computational Physiology, Simula Research Laboratory, Fornebu, Norway
| | | | - Vilde Strøm
- Department of Computational Physiology, Simula Research Laboratory, Fornebu, Norway
| | - Ilangko Balasingham
- Intervention Centre, Oslo University Hospital, Oslo, Norway; Department of Electronic Systems, Norwegian University of Science and Technology, Trondheim, Norway
| | | |
Collapse
|
44
|
Balaban G, Costa CM, Porter B, Halliday B, Rinaldi CA, Prasad S, Plank G, Ismail TF, Bishop MJ. 3D Electrophysiological Modeling of Interstitial Fibrosis Networks and Their Role in Ventricular Arrhythmias in Non-Ischemic Cardiomyopathy. IEEE Trans Biomed Eng 2020; 67:3125-3133. [PMID: 32275581 PMCID: PMC7116885 DOI: 10.1109/tbme.2020.2976924] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
Abstract
OBJECTIVE Interstitial fibrosis is a pathological expansion of the heart's inter-cellular collagen matrix. It is a potential complication of nonischemic cardiomyopathy (NICM), a class of diseases involving electrical and or mechanical dysfunction of cardiac tissue not caused by atherosclerosis. Patients with NICM and interstitial fibrosis often suffer from life threatening arrhythmias, which we aim to simulate in this study. METHODS Our methodology builds on an efficient discrete finite element (DFE) method which allows for the representation of fibrosis as infinitesimal splits in a mesh. We update the DFE method with a local connectivity analysis which creates a consistent topology in the fibrosis network. This is particularly important in nonischemic disease due to the potential presence of large and contiguous fibrotic regions and therefore potentially complex fibrosis networks. RESULTS In experiments with an image-based model, we demonstrate that our methodology is able to simulate reentrant electrical events associated with cardiac arrhythmias. These reentries depended crucially upon sufficient fibrosis density, which was marked by conduction slowing at high pacing rates. We also created a 2D test-case which demonstrated that fibrosis topologies can modulate transient conduction block, and thereby reentrant activations. CONCLUSION Ventricular arrhythmias due to interstitial fibrosis in NICM can be efficiently simulated using our methods in medical image based geometries. Furthermore, fibrosis topology modulates transient conduction block, and should be accounted for in electrophysiological simulations with interstitial fibrosis. SIGNIFICANCE Our study provides methodology which has the potential to predict arrhythmias and to optimize treatments non-invasively for nonischemic cardiomyopathies.
Collapse
|
45
|
Del Corso G, Verzicco R, Viola F. Sensitivity analysis of an electrophysiology model for the left ventricle. J R Soc Interface 2020; 17:20200532. [PMID: 33109017 DOI: 10.1098/rsif.2020.0532] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
Modelling the cardiac electrophysiology entails dealing with the uncertainties related to the input parameters such as the heart geometry and the electrical conductivities of the tissues, thus calling for an uncertainty quantification (UQ) of the results. Since the chambers of the heart have different shapes and tissues, in order to make the problem affordable, here we focus on the left ventricle with the aim of identifying which of the uncertain inputs mostly affect its electrophysiology. In a first phase, the uncertainty of the input parameters is evaluated using data available from the literature and the output quantities of interest (QoIs) of the problem are defined. According to the polynomial chaos expansion, a training dataset is then created by sampling the parameter space using a quasi-Monte Carlo method whereas a smaller independent dataset is used for the validation of the resulting metamodel. The latter is exploited to run a global sensitivity analysis with nonlinear variance-based indices and thus reduce the input parameter space accordingly. Thereafter, the uncertainty probability distribution of the QoIs are evaluated using a direct UQ strategy on a larger dataset and the results discussed in the light of the medical knowledge.
Collapse
Affiliation(s)
| | - Roberto Verzicco
- Gran Sasso Science Institute (GSSI), L'Aquila, Italy.,University of Rome Tor Vergata, Rome, Italy.,POF Group, University of Twente, Enschede, The Netherlands
| | | |
Collapse
|
46
|
Luongo G, Azzolin L, Rivolta MW, Sassi R, Martinez JP, Laguna P, Dossel O, Loewe A. Non-Invasive Identification of Atrial Fibrillation Driver Location Using the 12-lead ECG: Pulmonary Vein Rotors vs. other Locations. ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL INTERNATIONAL CONFERENCE 2020; 2020:410-413. [PMID: 33018015 DOI: 10.1109/embc44109.2020.9176135] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
Atrial fibrillation (AF) is an irregular heart rhythm due to disorganized atrial electrical activity, often sustained by rotational drivers called rotors. In the present work, we sought to characterize and discriminate whether simulated single stable rotors are located in the pulmonary veins (PVs) or not, only by using non-invasive signals (i.e., the 12-lead ECG). Several features have been extracted from the signals, such as Hjort descriptors, recurrence quantification analysis (RQA), and principal component analysis. All the extracted features have shown significant discriminatory power, with particular emphasis to the RQA parameters. A decision tree classifier achieved 98.48% accuracy, 83.33% sensitivity, and 100% specificity on simulated data.Clinical Relevance-This study might guide ablation procedures, suggesting doctors to proceed directly in some patients with a pulmonary veins isolation, and avoiding the prior use of an invasive atrial mapping system.
Collapse
|
47
|
Boyle PM, Del Álamo JC, Akoum N. Fibrosis, atrial fibrillation and stroke: clinical updates and emerging mechanistic models. Heart 2020; 107:99-105. [PMID: 33097562 DOI: 10.1136/heartjnl-2020-317455] [Citation(s) in RCA: 28] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 05/27/2020] [Revised: 09/21/2020] [Accepted: 09/25/2020] [Indexed: 02/06/2023] Open
Abstract
The current paradigm of stroke risk assessment and mitigation in patients with atrial fibrillation (AF) is centred around clinical risk factors which, in the presence of AF, lead to thrombus formation. The mechanisms by which these clinical risk factors lead to thromboembolism, including any role played by atrial fibrosis, are not understood. In patients who had embolic stroke of undetermined source (ESUS), the problem is compounded by the absence of AF in a majority of patients despite long-term monitoring. Atrial fibrosis has emerged as a unifying mechanism that independently provides a substrate for arrhythmia and thrombus formation. Fibrosis-based computational models of AF initiation and maintenance promise to identify therapeutic targets in catheter ablation. In ESUS, fibrosis is also increasingly recognised as a major risk factor, but the underlying mechanism of this correlation is unclear. Simulations have uncovered potential vulnerability to arrhythmia induction in patients who had ESUS. Likewise, computational models of fluid dynamics representing blood flow in the left atrium and left atrium appendage have improved our understanding of thrombus formation, in particular left atrium appendage shapes and blood flow changes influenced by atrial remodelling. Multiscale modelling of blood flow dynamics based on structural fibrotic and morphological changes with associated cellular and tissue electrical remodelling leading to electromechanical abnormalities holds tremendous promise in providing a mechanistic understanding of the clinical problem of thromboembolisation. We present a review of clinical knowledge alongside computational modelling frameworks and conclude with a vision of a future paradigm integrating simulations in formulating personalised treatment plans for each patient.
Collapse
Affiliation(s)
- Patrick M Boyle
- Bioengineering, University of Washington, Seattle, Washington, USA
| | - Juan Carlos Del Álamo
- Mechanical Engineering, University of Washington College of Engineering, Seattle, Washington, USA
| | - Nazem Akoum
- Cardiology, University of Washington School of Medicine, Seattle, Washington, USA
| |
Collapse
|
48
|
Monaci S, Strocchi M, Rodero C, Gillette K, Whitaker J, Rajani R, Rinaldi CA, O'Neill M, Plank G, King A, Bishop MJ. In-silico pace-mapping using a detailed whole torso model and implanted electronic device electrograms for more efficient ablation planning. Comput Biol Med 2020; 125:104005. [PMID: 32971325 DOI: 10.1016/j.compbiomed.2020.104005] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2020] [Revised: 09/07/2020] [Accepted: 09/07/2020] [Indexed: 10/23/2022]
Abstract
BACKGROUND Pace-mapping is a commonly used electrophysiological (EP) procedure which aims to identify exit sites of ventricular tachycardia (VT) by matching ventricular activation patterns (assessed by QRS morphology) at specific pacing locations with activation during VT. However, long procedure durations and the need for VT induction render this technique non-optimal. To demonstrate the potential of in-silico pace-mapping, using stored electrogram (EGM) recordings of clinical VT from implanted devices to guide pre-procedural ablation planning. METHOD Six scar-related VT episodes were simulated in a 3D torso model reconstructed from computed tomography (CT) imaging data, including three different infarct anatomies mapped from infarcted porcine imaging data. In-silico pace-mapping was performed to localise VT exit sites and isthmuses by using 12-lead electrocardiogram (ECG) signals and different combinations of EGM sensing vectors from implanted devices, through the creation of conventional correlation maps and reference-less maps. RESULTS Our in-silico platform was successful in identifying VT exit sites for a variety of different VT morphologies from both ECG correlation maps and corresponding EGM maps, with the latter dependent upon the number of sensing vectors used. We also showed the added utility of both ECG and EGM reference-less pace-mapping for the identification of slow-conducting isthmuses, uncovering the optimal algorithm parameters. Finally, EGM-based pace-mapping was shown to be more dependent upon the mapped surface (epicardial/endocardial), relative to the VT origin. CONCLUSIONS In-silico pace-mapping can be used along with EGMs from implanted devices to localise VT ablation targets in pre-procedural planning.
Collapse
Affiliation(s)
| | | | | | | | | | - Ronak Rajani
- King's College London, London, United Kingdom; Guy's and St Thomas' Hospital, London, United Kingdom
| | - Christopher A Rinaldi
- King's College London, London, United Kingdom; Guy's and St Thomas' Hospital, London, United Kingdom
| | | | | | - Andrew King
- King's College London, London, United Kingdom
| | | |
Collapse
|
49
|
Roney CH, Beach ML, Mehta AM, Sim I, Corrado C, Bendikas R, Solis-Lemus JA, Razeghi O, Whitaker J, O’Neill L, Plank G, Vigmond E, Williams SE, O’Neill MD, Niederer SA. In silico Comparison of Left Atrial Ablation Techniques That Target the Anatomical, Structural, and Electrical Substrates of Atrial Fibrillation. Front Physiol 2020; 11:1145. [PMID: 33041850 PMCID: PMC7526475 DOI: 10.3389/fphys.2020.572874] [Citation(s) in RCA: 28] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2020] [Accepted: 08/18/2020] [Indexed: 12/17/2022] Open
Abstract
Catheter ablation therapy for persistent atrial fibrillation (AF) typically includes pulmonary vein isolation (PVI) and may include additional ablation lesions that target patient-specific anatomical, electrical, or structural features. Clinical centers employ different ablation strategies, which use imaging data together with electroanatomic mapping data, depending on data availability. The aim of this study was to compare ablation techniques across a virtual cohort of AF patients. We constructed 20 paroxysmal and 30 persistent AF patient-specific left atrial (LA) bilayer models incorporating fibrotic remodeling from late-gadolinium enhancement (LGE) MRI scans. AF was simulated and post-processed using phase mapping to determine electrical driver locations over 15 s. Six different ablation approaches were tested: (i) PVI alone, modeled as wide-area encirclement of the pulmonary veins; PVI together with: (ii) roof and inferior lines to model posterior wall box isolation; (iii) isolating the largest fibrotic area (identified by LGE-MRI); (iv) isolating all fibrotic areas; (v) isolating the largest driver hotspot region [identified as high simulated phase singularity (PS) density]; and (vi) isolating all driver hotspot regions. Ablation efficacy was assessed to predict optimal ablation therapies for individual patients. We subsequently trained a random forest classifier to predict ablation response using (a) imaging metrics alone, (b) imaging and electrical metrics, or (c) imaging, electrical, and ablation lesion metrics. The optimal ablation approach resulting in termination, or if not possible atrial tachycardia (AT), varied among the virtual patient cohort: (i) 20% PVI alone, (ii) 6% box ablation, (iii) 2% largest fibrosis area, (iv) 4% all fibrosis areas, (v) 2% largest driver hotspot, and (vi) 46% all driver hotspots. Around 20% of cases remained in AF for all ablation strategies. The addition of patient-specific and ablation pattern specific lesion metrics to the trained random forest classifier improved predictive capability from an accuracy of 0.73 to 0.83. The trained classifier results demonstrate that the surface areas of pre-ablation driver regions and of fibrotic tissue not isolated by the proposed ablation strategy are both important for predicting ablation outcome. Overall, our study demonstrates the need to select the optimal ablation strategy for each patient. It suggests that both patient-specific fibrosis properties and driver locations are important for planning ablation approaches, and the distribution of lesions is important for predicting an acute response.
Collapse
Affiliation(s)
- Caroline H. Roney
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom
| | - Marianne L. Beach
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom
| | - Arihant M. Mehta
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom
| | - Iain Sim
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom
| | - Cesare Corrado
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom
| | - Rokas Bendikas
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom
| | - Jose A. Solis-Lemus
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom
| | - Orod Razeghi
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom
| | - John Whitaker
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom
| | - Louisa O’Neill
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom
| | - Gernot Plank
- Department of Biophysics, Medical University of Graz, Graz, Austria
| | - Edward Vigmond
- IHU Liryc, Electrophysiology and Heart Modeling Institute, Fondation Bordeaux Université, Bordeaux, France
| | - Steven E. Williams
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom
| | - Mark D. O’Neill
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom
| | - Steven A. Niederer
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom
| |
Collapse
|
50
|
Gemmell PM, Gillette K, Balaban G, Rajani R, Vigmond EJ, Plank G, Bishop MJ. A computational investigation into rate-dependant vectorcardiogram changes due to specific fibrosis patterns in non-ischæmic dilated cardiomyopathy. Comput Biol Med 2020; 123:103895. [PMID: 32741753 PMCID: PMC7429989 DOI: 10.1016/j.compbiomed.2020.103895] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2020] [Revised: 06/12/2020] [Accepted: 06/27/2020] [Indexed: 01/13/2023]
Abstract
Patients with scar-associated fibrotic tissue remodelling are at greater risk of ventricular arrhythmic events, but current methods to detect the presence of such remodelling require invasive procedures. We present here a potential method to detect the presence, location and dimensions of scar using pacing-dependent changes in the vectorcardiogram (VCG). Using a clinically-derived whole-torso computational model, simulations were conducted at both slow and rapid pacing for a variety of scar patterns within the myocardium, with various VCG-derived metrics being calculated, with changes in these metrics being assessed for their ability to discern the presence and size of scar. Our results indicate that differences in the dipole angle at the end of the QRS complex and differences in the QRS area and duration may be used to predict scar properties. Using machine learning techniques, we were also able to predict the location of the scar to high accuracy, using only these VCG-derived rate-dependent changes as input. Such a non-invasive predictive tool for the presence of scar represents a potentially useful clinical tool for identifying patients at arrhythmic risk.
Collapse
Affiliation(s)
- Philip M Gemmell
- King's College London, St. Thomas' Hospital North Wing, London, SE1 7EH, UK.
| | - Karli Gillette
- Medical University of Graz, Division of Biophysics, Neue Stiftingtalstraße 6(MC1.D.)/IV, 8010 Graz, Austria
| | - Gabriel Balaban
- University of Oslo, Research Group for Biomedical Infomatics, Gaustadalléen 23B 0373 Oslo, Norway
| | - Ronak Rajani
- King's College London, St. Thomas' Hospital North Wing, London, SE1 7EH, UK
| | - Edward J Vigmond
- University of Bordeaux, IHU Liryc, Site Hopital Xavier Arnozan, Avenue de Haut-Leveque, 33604 Pessac, France
| | - Gernot Plank
- Medical University of Graz, Division of Biophysics, Neue Stiftingtalstraße 6(MC1.D.)/IV, 8010 Graz, Austria
| | - Martin J Bishop
- King's College London, St. Thomas' Hospital North Wing, London, SE1 7EH, UK
| |
Collapse
|