1
|
Pagani S, Manzoni A. Enabling forward uncertainty quantification and sensitivity analysis in cardiac electrophysiology by reduced order modeling and machine learning. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2021; 37:e3450. [PMID: 33599106 PMCID: PMC8244126 DOI: 10.1002/cnm.3450] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Revised: 02/05/2021] [Accepted: 02/07/2021] [Indexed: 06/12/2023]
Abstract
We present a new, computationally efficient framework to perform forward uncertainty quantification (UQ) in cardiac electrophysiology. We consider the monodomain model to describe the electrical activity in the cardiac tissue, coupled with the Aliev-Panfilov model to characterize the ionic activity through the cell membrane. We address a complete forward UQ pipeline, including both: (i) a variance-based global sensitivity analysis for the selection of the most relevant input parameters, and (ii) a way to perform uncertainty propagation to investigate the impact of intra-subject variability on outputs of interest depending on the cardiac potential. Both tasks exploit stochastic sampling techniques, thus implying overwhelming computational costs because of the huge amount of queries to the high-fidelity, full-order computational model obtained by approximating the coupled monodomain/Aliev-Panfilov system through the finite element method. To mitigate this computational burden, we replace the full-order model with computationally inexpensive projection-based reduced-order models (ROMs) aimed at reducing the state-space dimensionality. Resulting approximation errors on the outputs of interest are finally taken into account through artificial neural network (ANN)-based models, enhancing the accuracy of the whole UQ pipeline. Numerical results show that the proposed physics-based ROMs outperform regression-based emulators relying on ANNs built with the same amount of training data, in terms of both numerical accuracy and overall computational efficiency.
Collapse
Affiliation(s)
- Stefano Pagani
- MOX, Dipartimento di MatematicaPolitecnico di MilanoMilanItaly
| | - Andrea Manzoni
- MOX, Dipartimento di MatematicaPolitecnico di MilanoMilanItaly
| |
Collapse
|
2
|
Campos FO, Orini M, Arnold R, Whitaker J, O'Neill M, Razavi R, Plank G, Hanson B, Porter B, Rinaldi CA, Gill J, Lambiase PD, Taggart P, Bishop MJ. Assessing the ability of substrate mapping techniques to guide ventricular tachycardia ablation using computational modelling. Comput Biol Med 2021; 130:104214. [PMID: 33476992 DOI: 10.1016/j.compbiomed.2021.104214] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2020] [Revised: 01/05/2021] [Accepted: 01/05/2021] [Indexed: 01/19/2023]
Abstract
BACKGROUND Identification of targets for ablation of post-infarction ventricular tachycardias (VTs) remains challenging, often requiring arrhythmia induction to delineate the reentrant circuit. This carries a risk for the patient and may not be feasible. Substrate mapping has emerged as a safer strategy to uncover arrhythmogenic regions. However, VT recurrence remains common. GOAL To use computer simulations to assess the ability of different substrate mapping approaches to identify VT exit sites. METHODS A 3D computational model of the porcine post-infarction heart was constructed to simulate VT and paced rhythm. Electroanatomical maps were constructed based on endocardial electrogram features and the reentry vulnerability index (RVI - a metric combining activation (AT) and repolarization timings to identify tissue susceptibility to reentry). Since scar transmurality in our model was not homogeneous, parameters derived from all signals (including dense scar regions) were used in the analysis. Potential ablation targets obtained from each electroanatomical map during pacing were compared to the exit site detected during VT mapping. RESULTS Simulation data showed that voltage cut-offs applied to bipolar electrograms could delineate the scar, but not the VT circuit. Electrogram fractionation had the highest correlation with scar transmurality. The RVI identified regions closest to VT exit site but was outperformed by AT gradients combined with voltage cut-offs. The performance of all metrics was affected by pacing location. CONCLUSIONS Substrate mapping could provide information about the infarct, but the directional dependency on activation should be considered. Activation-repolarization metrics have utility in safely identifying VT targets, even with non-transmural scars.
Collapse
Affiliation(s)
- Fernando O Campos
- School of Biomedical Engineering and Imaging Sciences, Rayne Institute, 4th Floor, Lambeth Wing, St. Thomas' Hospital, Westminster Bridge Road, London, SE1 7EH, United Kingdom.
| | - Michele Orini
- Institute of Cardiovascular Science, University College London, London, United Kingdom; Electrophysiology Department, Barts Heart Centre, St Bartholomew's Hospital, London, United Kingdom
| | - Robert Arnold
- Gottfried Schatz Research Center (for Cell Signaling, Metabolism and Aging), Division of Biophysics, Graz, Austria
| | - John Whitaker
- School of Biomedical Engineering and Imaging Sciences, Rayne Institute, 4th Floor, Lambeth Wing, St. Thomas' Hospital, Westminster Bridge Road, London, SE1 7EH, United Kingdom
| | - Mark O'Neill
- School of Biomedical Engineering and Imaging Sciences, Rayne Institute, 4th Floor, Lambeth Wing, St. Thomas' Hospital, Westminster Bridge Road, London, SE1 7EH, United Kingdom
| | - Reza Razavi
- School of Biomedical Engineering and Imaging Sciences, Rayne Institute, 4th Floor, Lambeth Wing, St. Thomas' Hospital, Westminster Bridge Road, London, SE1 7EH, United Kingdom
| | - Gernot Plank
- Gottfried Schatz Research Center (for Cell Signaling, Metabolism and Aging), Division of Biophysics, Graz, Austria
| | - Ben Hanson
- Department of Mechanical Engineering, University College London, London, United Kingdom
| | - Bradley Porter
- School of Biomedical Engineering and Imaging Sciences, Rayne Institute, 4th Floor, Lambeth Wing, St. Thomas' Hospital, Westminster Bridge Road, London, SE1 7EH, United Kingdom; Department of Cardiology, Guys and St Thomas' NHS Trust, London, United Kingdom
| | | | - Jaswinder Gill
- School of Biomedical Engineering and Imaging Sciences, Rayne Institute, 4th Floor, Lambeth Wing, St. Thomas' Hospital, Westminster Bridge Road, London, SE1 7EH, United Kingdom; Department of Cardiology, Guys and St Thomas' NHS Trust, London, United Kingdom
| | - Pier D Lambiase
- Institute of Cardiovascular Science, University College London, London, United Kingdom; Electrophysiology Department, Barts Heart Centre, St Bartholomew's Hospital, London, United Kingdom
| | - Peter Taggart
- Institute of Cardiovascular Science, University College London, London, United Kingdom; Electrophysiology Department, Barts Heart Centre, St Bartholomew's Hospital, London, United Kingdom
| | - Martin J Bishop
- School of Biomedical Engineering and Imaging Sciences, Rayne Institute, 4th Floor, Lambeth Wing, St. Thomas' Hospital, Westminster Bridge Road, London, SE1 7EH, United Kingdom
| |
Collapse
|
3
|
Corral-Acero J, Margara F, Marciniak M, Rodero C, Loncaric F, Feng Y, Gilbert A, Fernandes JF, Bukhari HA, Wajdan A, Martinez MV, Santos MS, Shamohammdi M, Luo H, Westphal P, Leeson P, DiAchille P, Gurev V, Mayr M, Geris L, Pathmanathan P, Morrison T, Cornelussen R, Prinzen F, Delhaas T, Doltra A, Sitges M, Vigmond EJ, Zacur E, Grau V, Rodriguez B, Remme EW, Niederer S, Mortier P, McLeod K, Potse M, Pueyo E, Bueno-Orovio A, Lamata P. The 'Digital Twin' to enable the vision of precision cardiology. Eur Heart J 2020; 41:4556-4564. [PMID: 32128588 PMCID: PMC7774470 DOI: 10.1093/eurheartj/ehaa159] [Citation(s) in RCA: 206] [Impact Index Per Article: 51.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 09/06/2019] [Revised: 11/29/2019] [Accepted: 02/24/2020] [Indexed: 12/26/2022] Open
Abstract
Providing therapies tailored to each patient is the vision of precision medicine, enabled by the increasing ability to capture extensive data about individual patients. In this position paper, we argue that the second enabling pillar towards this vision is the increasing power of computers and algorithms to learn, reason, and build the 'digital twin' of a patient. Computational models are boosting the capacity to draw diagnosis and prognosis, and future treatments will be tailored not only to current health status and data, but also to an accurate projection of the pathways to restore health by model predictions. The early steps of the digital twin in the area of cardiovascular medicine are reviewed in this article, together with a discussion of the challenges and opportunities ahead. We emphasize the synergies between mechanistic and statistical models in accelerating cardiovascular research and enabling the vision of precision medicine.
Collapse
Affiliation(s)
| | - Francesca Margara
- Department of Computer Science, British Heart Foundation Centre of Research Excellence, University of Oxford, Oxford, UK
| | - Maciej Marciniak
- Department of Biomedical Engineering, Division of Imaging Sciences and Biomedical Engineering, King’s College London, London, UK
| | - Cristobal Rodero
- Department of Biomedical Engineering, Division of Imaging Sciences and Biomedical Engineering, King’s College London, London, UK
| | - Filip Loncaric
- Institut Clínic Cardiovascular, Hospital Clínic, Universitat de Barcelona, Institut d’Investigacions Biomèdiques August Pi i Sunyer (IDIBAPS), Barcelona, Spain
| | - Yingjing Feng
- IHU Liryc, Electrophysiology and Heart Modeling Institute, fondation Bordeaux Université, Pessac-Bordeaux F-33600, France
- IMB, UMR 5251, University of Bordeaux, Talence F-33400, France
| | | | - Joao F Fernandes
- Department of Biomedical Engineering, Division of Imaging Sciences and Biomedical Engineering, King’s College London, London, UK
| | - Hassaan A Bukhari
- IMB, UMR 5251, University of Bordeaux, Talence F-33400, France
- Aragón Institute of Engineering Research, Universidad de Zaragoza, IIS Aragón, Zaragoza, Spain
| | - Ali Wajdan
- The Intervention Centre, Oslo University Hospital, Rikshospitalet, Oslo, Norway
| | | | | | - Mehrdad Shamohammdi
- CARIM School for Cardiovascular Diseases, Maastricht University, Maastricht, The Netherlands
| | - Hongxing Luo
- CARIM School for Cardiovascular Diseases, Maastricht University, Maastricht, The Netherlands
| | - Philip Westphal
- Medtronic PLC, Bakken Research Center, Maastricht, the Netherlands
| | - Paul Leeson
- Radcliffe Department of Medicine, Division of Cardiovascular Medicine, Oxford Cardiovascular Clinical Research Facility, John Radcliffe Hospital, University of Oxford, Oxford, UK
| | - Paolo DiAchille
- Healthcare and Life Sciences Research, IBM T.J. Watson Research Center, Yorktown Heights, NY, USA
| | - Viatcheslav Gurev
- Healthcare and Life Sciences Research, IBM T.J. Watson Research Center, Yorktown Heights, NY, USA
| | - Manuel Mayr
- King’s British Heart Foundation Centre, King’s College London, London, UK
| | - Liesbet Geris
- Virtual Physiological Human Institute, Leuven, Belgium
| | - Pras Pathmanathan
- Center for Devices and Radiological Health, U.S. Food and Drug Administration, Silver Spring, MD, USA
| | - Tina Morrison
- Center for Devices and Radiological Health, U.S. Food and Drug Administration, Silver Spring, MD, USA
| | | | - Frits Prinzen
- CARIM School for Cardiovascular Diseases, Maastricht University, Maastricht, The Netherlands
| | - Tammo Delhaas
- CARIM School for Cardiovascular Diseases, Maastricht University, Maastricht, The Netherlands
| | - Ada Doltra
- Institut Clínic Cardiovascular, Hospital Clínic, Universitat de Barcelona, Institut d’Investigacions Biomèdiques August Pi i Sunyer (IDIBAPS), Barcelona, Spain
| | - Marta Sitges
- Institut Clínic Cardiovascular, Hospital Clínic, Universitat de Barcelona, Institut d’Investigacions Biomèdiques August Pi i Sunyer (IDIBAPS), Barcelona, Spain
- CIBERCV, Instituto de Salud Carlos III, (CB16/11/00354), CERCA Programme/Generalitat de, Catalunya, Spain
| | - Edward J Vigmond
- IHU Liryc, Electrophysiology and Heart Modeling Institute, fondation Bordeaux Université, Pessac-Bordeaux F-33600, France
- IMB, UMR 5251, University of Bordeaux, Talence F-33400, France
| | - Ernesto Zacur
- Department of Engineering Science, University of Oxford, Oxford, UK
| | - Vicente Grau
- Department of Engineering Science, University of Oxford, Oxford, UK
| | - Blanca Rodriguez
- Department of Computer Science, British Heart Foundation Centre of Research Excellence, University of Oxford, Oxford, UK
| | - Espen W Remme
- The Intervention Centre, Oslo University Hospital, Rikshospitalet, Oslo, Norway
| | - Steven Niederer
- Department of Biomedical Engineering, Division of Imaging Sciences and Biomedical Engineering, King’s College London, London, UK
| | | | | | - Mark Potse
- IHU Liryc, Electrophysiology and Heart Modeling Institute, fondation Bordeaux Université, Pessac-Bordeaux F-33600, France
- IMB, UMR 5251, University of Bordeaux, Talence F-33400, France
- Inria Bordeaux Sud-Ouest, CARMEN team, Talence F-33400, France
| | - Esther Pueyo
- Aragón Institute of Engineering Research, Universidad de Zaragoza, IIS Aragón, Zaragoza, Spain
- CIBER in Bioengineering, Biomaterials and Nanomedicine (CIBER‐BBN), Madrid, Spain
| | - Alfonso Bueno-Orovio
- Department of Computer Science, British Heart Foundation Centre of Research Excellence, University of Oxford, Oxford, UK
| | - Pablo Lamata
- Department of Biomedical Engineering, Division of Imaging Sciences and Biomedical Engineering, King’s College London, London, UK
| |
Collapse
|
4
|
Ramírez WA, Gizzi A, Sack KL, Guccione JM, Hurtado DE. In-silico study of the cardiac arrhythmogenic potential of biomaterial injection therapy. Sci Rep 2020; 10:12990. [PMID: 32737400 PMCID: PMC7395773 DOI: 10.1038/s41598-020-69900-4] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2019] [Accepted: 06/19/2020] [Indexed: 02/06/2023] Open
Abstract
Biomaterial injection is a novel therapy to treat ischemic heart failure (HF) that has shown to reduce remodeling and restore cardiac function in recent preclinical studies. While the effect of biomaterial injection in reducing mechanical wall stress has been recently demonstrated, the influence of biomaterials on the electrical behavior of treated hearts has not been elucidated. In this work, we developed computational models of swine hearts to study the electrophysiological vulnerability associated with biomaterial injection therapy. The propagation of action potentials on realistic biventricular geometries was simulated by numerically solving the monodomain electrophysiology equations on anatomically-detailed models of normal, HF untreated, and HF treated hearts. Heart geometries were constructed from high-resolution magnetic resonance images (MRI) where the healthy, peri-infarcted, infarcted and gel regions were identified, and the orientation of cardiac fibers was informed from diffusion-tensor MRI. Regional restitution properties in each case were evaluated by constructing a probability density function of the action potential duration (APD) at different cycle lengths. A comparative analysis of the ventricular fibrillation (VF) dynamics for every heart was carried out by measuring the number of filaments formed after wave braking. Our results suggest that biomaterial injection therapy does not affect the regional dispersion of repolarization when comparing untreated and treated failing hearts. Further, we found that the treated failing heart is more prone to sustain VF than the normal heart, and is at least as susceptible to sustained VF as the untreated failing heart. Moreover, we show that the main features of VF dynamics in a treated failing heart are not affected by the level of electrical conductivity of the biogel injectates. This work represents a novel proof-of-concept study demonstrating the feasibility of computer simulations of the heart in understanding the arrhythmic behavior in novel therapies for HF.
Collapse
Affiliation(s)
- William A Ramírez
- Department of Structural and Geotechnical Engineering, School of Engineering, Pontificia Universidad Católica de Chile, Santiago, Chile
| | - Alessio Gizzi
- Nonlinear Physics and Mathematical Modeling Lab, Department of Engineering, Campus Bio-Medico University of Rome, Rome, Italy
| | - Kevin L Sack
- Department of Surgery, University of California at San Francisco, San Francisco, CA, USA
- Division of Biomedical Engineering, Department of Human Biology, University of Cape Town, Cape Town, South Africa
| | - Julius M Guccione
- Department of Surgery, University of California at San Francisco, San Francisco, CA, USA
| | - Daniel E Hurtado
- Department of Structural and Geotechnical Engineering, School of Engineering, Pontificia Universidad Católica de Chile, Santiago, Chile.
- Institute for Biological and Medical Engineering, Schools of Engineering, Medicine and Biological Sciences, Pontificia Universidad Católica de Chile, Santiago, Chile.
- Millennium Nucleus for Cardiovascular Magnetic Resonance, Santiago, Chile.
| |
Collapse
|
5
|
Evaluation of the reentry vulnerability index to predict ventricular tachycardia circuits using high-density contact mapping. Heart Rhythm 2019; 17:576-583. [PMID: 31751771 PMCID: PMC7105818 DOI: 10.1016/j.hrthm.2019.11.013] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 08/11/2019] [Indexed: 11/24/2022]
Abstract
Background Identifying arrhythmogenic sites to improve ventricular tachycardia (VT) ablation outcomes remains unresolved. The reentry vulnerability index (RVI) combines activation and repolarization timings to identify sites critical for reentrant arrhythmia initiation without inducing VT. Objective The purpose of this study was to provide the first assessment of RVI’s capability to identify VT sites of origin using high-density contact mapping and comparison with other activation-repolarization markers of functional substrate. Methods Eighteen VT ablation patients (16 male; 72% ischemic) were studied. Unipolar electrograms were recorded during ventricular pacing and analyzed offline. Activation time (AT), activation–recovery interval (ARI), and repolarization time (RT) were measured. Vulnerability to reentry was mapped based on RVI and spatial distribution of AT, ARI, and RT. The distance from sites identified as vulnerable to reentry to the VT site of origin was measured, with distances <10 mm and >20 mm indicating accurate and inaccurate localization, respectively. Results The origins of 18 VTs (6 entrainment, 12 pace-mapping) were identified. RVI maps included 1012 (408–2098) (median, 1st–3rd quartiles) points per patient. RVI accurately localized 72.2% VT sites of origin, with median distance of 5.1 (3.2–10.1) mm. Inaccurate localization was significantly less frequent for RVI than AT (5.6% vs 33.3%; odds ratio 0.12; P = .035). Compared to RVI, distance to VT sites of origin was significantly larger for sites showing prolonged RT and ARI and were nonsignificantly larger for sites showing highest AT and ARI gradients. Conclusion RVI identifies vulnerable regions closest to VT sites of origin. Activation-repolarization metrics may improve VT substrate delineation and inform novel ablation strategies.
Collapse
|
6
|
Factors Promoting Conduction Slowing as Substrates for Block and Reentry in Infarcted Hearts. Biophys J 2019; 117:2361-2374. [PMID: 31521328 PMCID: PMC6990374 DOI: 10.1016/j.bpj.2019.08.008] [Citation(s) in RCA: 30] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2019] [Revised: 07/03/2019] [Accepted: 08/05/2019] [Indexed: 01/11/2023] Open
Abstract
The development of effective and safe therapies for scar-related ventricular tachycardias requires a detailed understanding of the mechanisms underlying the conduction block that initiates electrical re-entries associated with these arrhythmias. Conduction block has been often associated with electrophysiological changes that prolong action potential duration (APD) within the border zone (BZ) of chronically infarcted hearts. However, experimental evidence suggests that remodeling processes promoting conduction slowing as opposed to APD prolongation mark the chronic phase. In this context, the substrate for the initial block at the mouth of an isthmus/diastolic channel leading to ventricular tachycardia is unclear. The goal of this study was to determine whether electrophysiological parameters associated with conduction slowing can cause block and re-entry in the BZ. In silico experiments were conducted on two-dimensional idealized infarct tissue as well as on a cohort of postinfarction porcine left ventricular models constructed from ex vivo magnetic resonance imaging scans. Functional conduction slowing in the BZ was modeled by reducing sodium current density, whereas structural conduction slowing was represented by decreasing tissue conductivity and including fibrosis. The arrhythmogenic potential of APD prolongation was also tested as a basis for comparison. Within all models, the combination of reduced sodium current with structural remodeling more often degenerated into re-entry and, if so, was more likely to be sustained for more cycles. Although re-entries were also detected in experiments with prolonged APD, they were often not sustained because of the subsequent block caused by long-lasting repolarization. Functional and structural conditions associated with slow conduction rather than APD prolongation form a potent substrate for arrhythmogenesis at the isthmus/BZ of chronically infarcted hearts. Reduced excitability led to block while slow conduction shortened the wavelength of propagation, facilitating the sustenance of re-entries. These findings provide important insights for models of patient-specific risk stratification and therapy planning.
Collapse
|
7
|
Lopez-Perez A, Sebastian R, Izquierdo M, Ruiz R, Bishop M, Ferrero JM. Personalized Cardiac Computational Models: From Clinical Data to Simulation of Infarct-Related Ventricular Tachycardia. Front Physiol 2019; 10:580. [PMID: 31156460 PMCID: PMC6531915 DOI: 10.3389/fphys.2019.00580] [Citation(s) in RCA: 41] [Impact Index Per Article: 8.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2018] [Accepted: 04/25/2019] [Indexed: 12/20/2022] Open
Abstract
In the chronic stage of myocardial infarction, a significant number of patients develop life-threatening ventricular tachycardias (VT) due to the arrhythmogenic nature of the remodeled myocardium. Radiofrequency ablation (RFA) is a common procedure to isolate reentry pathways across the infarct scar that are responsible for VT. Unfortunately, this strategy show relatively low success rates; up to 50% of patients experience recurrent VT after the procedure. In the last decade, intensive research in the field of computational cardiac electrophysiology (EP) has demonstrated the ability of three-dimensional (3D) cardiac computational models to perform in-silico EP studies. However, the personalization and modeling of certain key components remain challenging, particularly in the case of the infarct border zone (BZ). In this study, we used a clinical dataset from a patient with a history of infarct-related VT to build an image-based 3D ventricular model aimed at computational simulation of cardiac EP, including detailed patient-specific cardiac anatomy and infarct scar geometry. We modeled the BZ in eight different ways by combining the presence or absence of electrical remodeling with four different levels of image-based patchy fibrosis (0, 10, 20, and 30%). A 3D torso model was also constructed to compute the ECG. Patient-specific sinus activation patterns were simulated and validated against the patient's ECG. Subsequently, the pacing protocol used to induce reentrant VTs in the EP laboratory was reproduced in-silico. The clinical VT was induced with different versions of the model and from different pacing points, thus identifying the slow conducting channel responsible for such VT. Finally, the real patient's ECG recorded during VT episodes was used to validate our simulation results and to assess different strategies to model the BZ. Our study showed that reduced conduction velocities and heterogeneity in action potential duration in the BZ are the main factors in promoting reentrant activity. Either electrical remodeling or fibrosis in a degree of at least 30% in the BZ were required to initiate VT. Moreover, this proof-of-concept study confirms the feasibility of developing 3D computational models for cardiac EP able to reproduce cardiac activation in sinus rhythm and during VT, using exclusively non-invasive clinical data.
Collapse
Affiliation(s)
- Alejandro Lopez-Perez
- Center for Research and Innovation in Bioengineering (Ci2B), Universitat Politècnica de València, Valencia, Spain
| | - Rafael Sebastian
- Computational Multiscale Simulation Lab (CoMMLab), Universitat de València, Valencia, Spain
| | - M Izquierdo
- INCLIVA Health Research Institute, Valencia, Spain.,Arrhythmia Unit, Cardiology Department, Hospital Clínico Universitario de Valencia, Valencia, Spain
| | - Ricardo Ruiz
- INCLIVA Health Research Institute, Valencia, Spain.,Arrhythmia Unit, Cardiology Department, Hospital Clínico Universitario de Valencia, Valencia, Spain
| | - Martin Bishop
- Division of Imaging Sciences & Biomedical Engineering, Department of Biomedical Engineering, King's College London, London, United Kingdom
| | - Jose M Ferrero
- Center for Research and Innovation in Bioengineering (Ci2B), Universitat Politècnica de València, Valencia, Spain
| |
Collapse
|
8
|
Campos FO, Orini M, Taggart P, Hanson B, Lambiase PD, Porter B, Rinaldi CA, Gill J, Bishop MJ. Characterizing the clinical implementation of a novel activation-repolarization metric to identify targets for catheter ablation of ventricular tachycardias using computational models. Comput Biol Med 2019; 108:263-275. [PMID: 31009930 PMCID: PMC6538827 DOI: 10.1016/j.compbiomed.2019.03.018] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2018] [Revised: 03/08/2019] [Accepted: 03/19/2019] [Indexed: 11/24/2022]
Abstract
Identification of targets for catheter ablation of ventricular tachycardias (VTs) remains a significant challenge. VTs are often driven by re-entrant circuits resulting from a complex interaction between the front (activation) and tail (repolarization) of the electrical wavefront. Most mapping techniques do not take into account the tissue repolarization which may hinder the detection of ablation targets. The re-entry vulnerability index (RVI), a recently proposed mapping procedure, incorporates both activation and repolarization times to uncover VT circuits. The method showed potential in a series of experiments, but it still requires further development to enable its incorporation into a clinical protocol. Here, in-silico experiments were conducted to thoroughly assess RVI maps constructed under clinically-relevant mapping conditions. Within idealized as well as anatomically realistic infarct models, we show that parameters of the algorithm such as the search radius can significantly alter the specificity and sensitivity of the RVI maps. When constructed on sparse grids obtained following various placements of clinical recording catheters, RVI maps can identify vulnerable regions as long as two electrodes were placed on both sides of the line of block. Moreover, maps computed during pacing without inducing VT can reveal areas of abnormal repolarization and slow conduction but not directly vulnerability. In conclusion, the RVI algorithm can detect re-entrant circuits during VT from low resolution mapping grids resembling the clinical setting. Furthermore, RVI maps may provide information about the underlying tissue electrophysiology to guide catheter ablation without the need of inducing potentially harmful VT during the clinical procedure. Finally, the ability of the RVI maps to identify vulnerable regions with specificity in high resolution computer models could potentially improve the prediction of optimal ablation targets of simulation-based strategies. Safe and accurate detection of targets for catheter ablation remains a challenge. We conducted a thorough assessment of the Re-entry Vulnerability Index (RVI). Parameters of the algorithm can alter the specificity and sensitivity of RVI maps. When constructed on sparse grids RVI maps could still detect arrhythmogenic sites. In absence of arrhythmia, RVI maps revealed abnormal sites, but not vulnerability.
Collapse
Affiliation(s)
- Fernando O Campos
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
| | - Michele Orini
- The Heart Hospital, University College London, London, United Kingdom; Institute of Cardiovascular Science, University College London, London, United Kingdom
| | - Peter Taggart
- The Heart Hospital, University College London, London, United Kingdom; Institute of Cardiovascular Science, University College London, London, United Kingdom
| | - Ben Hanson
- Department of Mechanical Engineering, University College London, London, United Kingdom
| | - Pier D Lambiase
- Institute of Cardiovascular Science, University College London, London, United Kingdom; Electrophysiology Department, Barts Heart Centre, St Bartholomew's Hospital, London, United Kingdom
| | - Bradley Porter
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom; Department of Cardiology, Guys and St Thomas' NHS Trust, London, United Kingdom
| | | | - Jaswinder Gill
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom; Department of Cardiology, Guys and St Thomas' NHS Trust, London, United Kingdom
| | - Martin J Bishop
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom. https://kclpure.kcl.ac.uk/portal/martin.bishop.html
| |
Collapse
|
9
|
Loewe A, Poremba E, Oesterlein T, Luik A, Schmitt C, Seemann G, Dössel O. Patient-Specific Identification of Atrial Flutter Vulnerability-A Computational Approach to Reveal Latent Reentry Pathways. Front Physiol 2019; 9:1910. [PMID: 30692934 PMCID: PMC6339942 DOI: 10.3389/fphys.2018.01910] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2018] [Accepted: 12/18/2018] [Indexed: 11/23/2022] Open
Abstract
Atypical atrial flutter (AFlut) is a reentrant arrhythmia which patients frequently develop after ablation for atrial fibrillation (AF). Indeed, substrate modifications during AF ablation can increase the likelihood to develop AFlut and it is clinically not feasible to reliably and sensitively test if a patient is vulnerable to AFlut. Here, we present a novel method based on personalized computational models to identify pathways along which AFlut can be sustained in an individual patient. We build a personalized model of atrial excitation propagation considering the anatomy as well as the spatial distribution of anisotropic conduction velocity and repolarization characteristics based on a combination of a priori knowledge on the population level and information derived from measurements performed in the individual patient. The fast marching scheme is employed to compute activation times for stimuli from all parts of the atria. Potential flutter pathways are then identified by tracing loops from wave front collision sites and constricting them using a geometric snake approach under consideration of the heterogeneous wavelength condition. In this way, all pathways along which AFlut can be sustained are identified. Flutter pathways can be instantiated by using an eikonal-diffusion phase extrapolation approach and a dynamic multifront fast marching simulation. In these dynamic simulations, the initial pattern eventually turns into the one driven by the dominant pathway, which is the only pathway that can be observed clinically. We assessed the sensitivity of the flutter pathway maps with respect to conduction velocity and its anisotropy. Moreover, we demonstrate the application of tailored models considering disease-specific repolarization properties (healthy, AF-remodeled, potassium channel mutations) as well as applicabiltiy on a clinical dataset. Finally, we tested how AFlut vulnerability of these substrates is modulated by exemplary antiarrhythmic drugs (amiodarone, dronedarone). Our novel method allows to assess the vulnerability of an individual patient to develop AFlut based on the personal anatomical, electrophysiological, and pharmacological characteristics. In contrast to clinical electrophysiological studies, our computational approach provides the means to identify all possible AFlut pathways and not just the currently dominant one. This allows to consider all relevant AFlut pathways when tailoring clinical ablation therapy in order to reduce the development and recurrence of AFlut.
Collapse
Affiliation(s)
- Axel Loewe
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | - Emanuel Poremba
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | - Tobias Oesterlein
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | - Armin Luik
- Medizinische Klinik IV, Städtisches Klinikum Karlsruhe, Karlsruhe, Germany
| | - Claus Schmitt
- Medizinische Klinik IV, Städtisches Klinikum Karlsruhe, Karlsruhe, Germany
| | - Gunnar Seemann
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
- Institute for Experimental Cardiovascular Medicine, University Heart Center Freiburg Bad Krozingen, Freiburg, Germany
- Faculty of Medicine, Albert-Ludwigs University, Freiburg, Germany
| | - Olaf Dössel
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| |
Collapse
|
10
|
Martin CA, Orini M, Srinivasan NT, Bhar-Amato J, Honarbakhsh S, Chow AW, Lowe MD, Ben-Simon R, Elliott PM, Taggart P, Lambiase PD. Assessment of a conduction-repolarisation metric to predict Arrhythmogenesis in right ventricular disorders. Int J Cardiol 2018; 271:75-80. [PMID: 29871808 PMCID: PMC6152588 DOI: 10.1016/j.ijcard.2018.05.063] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 03/29/2018] [Revised: 05/16/2018] [Accepted: 05/18/2018] [Indexed: 11/22/2022]
Abstract
BACKGROUND The re-entry vulnerability index (RVI) is a recently proposed activation-repolarization metric designed to quantify tissue susceptibility to re-entry. This study aimed to test feasibility of an RVI-based algorithm to predict the earliest endocardial activation site of ventricular tachycardia (VT) during electrophysiological studies and occurrence of haemodynamically significant ventricular arrhythmias in follow-up. METHODS Patients with Arrhythmogenic Right Ventricular Cardiomyopathy (ARVC) (n = 11), Brugada Syndrome (BrS) (n = 13) and focal RV outflow tract VT (n = 9) underwent programmed stimulation with unipolar electrograms recorded from a non-contact array in the RV. RESULTS Lowest values of RVI co-localised with VT earliest activation site in ARVC/BrS but not in focal VT. The distance between region of lowest RVI and site of VT earliest site (Dmin) was lower in ARVC/BrS than in focal VT (6.8 ± 6.7 mm vs 26.9 ± 13.3 mm, p = 0.005). ARVC/BrS patients with inducible VT had lower Global-RVI (RVIG) than those who were non-inducible (-54.9 ± 13.0 ms vs -35.9 ± 8.6 ms, p = 0.005) or those with focal VT (-30.6 ± 11.5 ms, p = 0.001). Patients were followed up for 112 ± 19 months. Those with clinical VT events had lower Global-RVI than both ARVC and BrS patients without VT (-54.5 ± 13.5 ms vs -36.2 ± 8.8 ms, p = 0.007) and focal VT patients (-30.6 ± 11.5 ms, p = 0.002). CONCLUSIONS RVI reliably identifies the earliest RV endocardial activation site of VT in BrS and ARVC but not focal ventricular arrhythmias and predicts the incidence of haemodynamically significant arrhythmias. Therefore, RVI may be of value in predicting VT exit sites and hence targeting of re-entrant arrhythmias.
Collapse
Affiliation(s)
- C A Martin
- Barts Heart Centre, Barts Health NHS Trust, West Smithfield, London EC1A 7BE, UK; Institute of Cardiovascular Science, University College London, Gower Street, London WC1E 6BT, UK
| | - M Orini
- Barts Heart Centre, Barts Health NHS Trust, West Smithfield, London EC1A 7BE, UK; Institute of Cardiovascular Science, University College London, Gower Street, London WC1E 6BT, UK
| | - N T Srinivasan
- Barts Heart Centre, Barts Health NHS Trust, West Smithfield, London EC1A 7BE, UK; Institute of Cardiovascular Science, University College London, Gower Street, London WC1E 6BT, UK
| | - J Bhar-Amato
- Institute of Cardiovascular Science, University College London, Gower Street, London WC1E 6BT, UK
| | - S Honarbakhsh
- Barts Heart Centre, Barts Health NHS Trust, West Smithfield, London EC1A 7BE, UK
| | - A W Chow
- Barts Heart Centre, Barts Health NHS Trust, West Smithfield, London EC1A 7BE, UK
| | - M D Lowe
- Barts Heart Centre, Barts Health NHS Trust, West Smithfield, London EC1A 7BE, UK
| | - R Ben-Simon
- Barts Heart Centre, Barts Health NHS Trust, West Smithfield, London EC1A 7BE, UK
| | - P M Elliott
- Barts Heart Centre, Barts Health NHS Trust, West Smithfield, London EC1A 7BE, UK
| | - P Taggart
- Institute of Cardiovascular Science, University College London, Gower Street, London WC1E 6BT, UK
| | - P D Lambiase
- Barts Heart Centre, Barts Health NHS Trust, West Smithfield, London EC1A 7BE, UK; Institute of Cardiovascular Science, University College London, Gower Street, London WC1E 6BT, UK.
| |
Collapse
|
11
|
Deng D, Nikolov P, Arevalo HJ, Trayanova NA. Optimal contrast-enhanced MRI image thresholding for accurate prediction of ventricular tachycardia using ex-vivo high resolution models. Comput Biol Med 2018; 102:426-432. [PMID: 30301573 PMCID: PMC6218273 DOI: 10.1016/j.compbiomed.2018.09.031] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2018] [Revised: 09/12/2018] [Accepted: 09/30/2018] [Indexed: 11/23/2022]
Abstract
Patient specific models created from contrast-enhanced (i.e. late-gadolinium, LGE) MRI images can be used for prediction of reentry location and clinical ablation planning. However, there is still a need for direct and systematic comparison between characteristics of ventricular tachycardia (VT) morphologies predicted in computational models and those acquired in clinical or experimental protocols. In this study, we aimed to: 1) assess the differences in VT morphologies predicted by modeling and recorded in experiments in terms of patterns and location of reentries, earliest and latest activation sites, and cycle lengths; and 2) define the optimal range of infarct tissue threshold values which provide best match between simulation and experimental results. To achieve these goals, we utilized LGE-MRI images from 4 swine hearts with inducible monomorphic VT. The images were segmented to identify non-infarcted myocardium, semi viable gray zone (GZ), and core scar based on pixel intensity. Several models were reconstructed from each LGE-MRI scan, with voxels of intensity between that of non-infarcted myocardium and 20-50% of the maximum intensity (in 10% increments) in the infarct region classified as GZ. VT induction was simulated in each model. Our simulation results showed that using GZ intensity thresholds of 20% or 30% resulted in the best match of simulated propagation patterns and reentry locations with those from the experiment. Overall, we matched 70% (7/10) morphologies for all the hearts. Our simulation shows that MRI-based computational models of hearts with myocardial infarction can accurately reproduce the majority of experimentally recorded post-infarction VTs.
Collapse
Affiliation(s)
- Dongdong Deng
- Institute for Computational Medicine, Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, USA
| | - Plamen Nikolov
- Institute for Computational Medicine, Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, USA
| | - Hermenegild J Arevalo
- Institute for Computational Medicine, Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, USA; Cardiac Modelling Department, Simula Research Laboratory, Fornebu, Norway
| | - Natalia A Trayanova
- Institute for Computational Medicine, Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, USA.
| |
Collapse
|
12
|
|
13
|
Connolly AJ, Bishop MJ. Computational Representations of Myocardial Infarct Scars and Implications for Arrhythmogenesis. CLINICAL MEDICINE INSIGHTS-CARDIOLOGY 2016; 10:27-40. [PMID: 27486348 PMCID: PMC4962962 DOI: 10.4137/cmc.s39708] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 04/19/2016] [Revised: 05/17/2016] [Accepted: 05/27/2016] [Indexed: 11/30/2022]
Abstract
Image-based computational modeling is becoming an increasingly used clinical tool to provide insight into the mechanisms of reentrant arrhythmias. In the context of ischemic heart disease, faithful representation of the electrophysiological properties of the infarct region within models is essential, due to the scars known for arrhythmic properties. Here, we review the different computational representations of the infarcted region, summarizing the experimental measurements upon which they are based. We then focus on the two most common representations of the scar core (complete insulator or electrically passive tissue) and perform simulations of electrical propagation around idealized infarct geometries. Our simulations highlight significant differences in action potential duration and focal effective refractory period (ERP) around the scar, driven by differences in electrotonic loading, depending on the choice of scar representation. Finally, a novel mechanism for arrhythmia induction, following a focal ectopic beat, is demonstrated, which relies on localized gradients in ERP directly caused by the electrotonic sink effects of the neighboring passive scar.
Collapse
Affiliation(s)
- Adam J Connolly
- Department of Imaging Sciences and Bioengineering, King's College London, St Thomas' Hospital, London, UK
| | - Martin J Bishop
- Department of Imaging Sciences and Bioengineering, King's College London, St Thomas' Hospital, London, UK
| |
Collapse
|
14
|
Arevalo HJ, Boyle PM, Trayanova NA. Computational rabbit models to investigate the initiation, perpetuation, and termination of ventricular arrhythmia. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2016; 121:185-94. [PMID: 27334789 DOI: 10.1016/j.pbiomolbio.2016.06.004] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/25/2016] [Accepted: 06/13/2016] [Indexed: 12/29/2022]
Abstract
Current understanding of cardiac electrophysiology has been greatly aided by computational work performed using rabbit ventricular models. This article reviews the contributions of multiscale models of rabbit ventricles in understanding cardiac arrhythmia mechanisms. This review will provide an overview of multiscale modeling of the rabbit ventricles. It will then highlight works that provide insights into the role of the conduction system, complex geometric structures, and heterogeneous cellular electrophysiology in diseased and healthy rabbit hearts to the initiation and maintenance of ventricular arrhythmia. Finally, it will provide an overview on the contributions of rabbit ventricular modeling on understanding the mechanisms underlying shock-induced defibrillation.
Collapse
Affiliation(s)
- Hermenegild J Arevalo
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, USA; Simula Research Laboratory, Oslo, Norway
| | - Patrick M Boyle
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, USA
| | - Natalia A Trayanova
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, USA.
| |
Collapse
|