1
|
Trayanova NA, Lyon A, Shade J, Heijman J. Computational modeling of cardiac electrophysiology and arrhythmogenesis: toward clinical translation. Physiol Rev 2024; 104:1265-1333. [PMID: 38153307 DOI: 10.1152/physrev.00017.2023] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2023] [Revised: 12/19/2023] [Accepted: 12/21/2023] [Indexed: 12/29/2023] Open
Abstract
The complexity of cardiac electrophysiology, involving dynamic changes in numerous components across multiple spatial (from ion channel to organ) and temporal (from milliseconds to days) scales, makes an intuitive or empirical analysis of cardiac arrhythmogenesis challenging. Multiscale mechanistic computational models of cardiac electrophysiology provide precise control over individual parameters, and their reproducibility enables a thorough assessment of arrhythmia mechanisms. This review provides a comprehensive analysis of models of cardiac electrophysiology and arrhythmias, from the single cell to the organ level, and how they can be leveraged to better understand rhythm disorders in cardiac disease and to improve heart patient care. Key issues related to model development based on experimental data are discussed, and major families of human cardiomyocyte models and their applications are highlighted. An overview of organ-level computational modeling of cardiac electrophysiology and its clinical applications in personalized arrhythmia risk assessment and patient-specific therapy of atrial and ventricular arrhythmias is provided. The advancements presented here highlight how patient-specific computational models of the heart reconstructed from patient data have achieved success in predicting risk of sudden cardiac death and guiding optimal treatments of heart rhythm disorders. Finally, an outlook toward potential future advances, including the combination of mechanistic modeling and machine learning/artificial intelligence, is provided. As the field of cardiology is embarking on a journey toward precision medicine, personalized modeling of the heart is expected to become a key technology to guide pharmaceutical therapy, deployment of devices, and surgical interventions.
Collapse
Affiliation(s)
- Natalia A Trayanova
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, Maryland, United States
- Alliance for Cardiovascular Diagnostic and Treatment Innovation, Johns Hopkins University, Baltimore, Maryland, United States
| | - Aurore Lyon
- Department of Biomedical Engineering, CARIM School for Cardiovascular Diseases, Maastricht University, Maastricht, The Netherlands
- Division of Heart and Lungs, Department of Medical Physiology, University Medical Centre Utrecht, Utrecht, The Netherlands
| | - Julie Shade
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, Maryland, United States
- Alliance for Cardiovascular Diagnostic and Treatment Innovation, Johns Hopkins University, Baltimore, Maryland, United States
| | - Jordi Heijman
- Department of Cardiology, CARIM School for Cardiovascular Diseases, Maastricht University, Maastricht, The Netherlands
| |
Collapse
|
2
|
Jaffery OA, Melki L, Slabaugh G, Good WW, Roney CH. A Review of Personalised Cardiac Computational Modelling Using Electroanatomical Mapping Data. Arrhythm Electrophysiol Rev 2024; 13:e08. [PMID: 38807744 PMCID: PMC11131150 DOI: 10.15420/aer.2023.25] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 10/11/2023] [Accepted: 12/27/2023] [Indexed: 05/30/2024] Open
Abstract
Computational models of cardiac electrophysiology have gradually matured during the past few decades and are now being personalised to provide patient-specific therapy guidance for improving suboptimal treatment outcomes. The predictive features of these personalised electrophysiology models hold the promise of providing optimal treatment planning, which is currently limited in the clinic owing to reliance on a population-based or average patient approach. The generation of a personalised electrophysiology model entails a sequence of steps for which a range of activation mapping, calibration methods and therapy simulation pipelines have been suggested. However, the optimal methods that can potentially constitute a clinically relevant in silico treatment are still being investigated and face limitations, such as uncertainty of electroanatomical data recordings, generation and calibration of models within clinical timelines and requirements to validate or benchmark the recovered tissue parameters. This paper is aimed at reporting techniques on the personalisation of cardiac computational models, with a focus on calibrating cardiac tissue conductivity based on electroanatomical mapping data.
Collapse
Affiliation(s)
- Ovais A Jaffery
- School of Engineering and Materials Science, Queen Mary University of London London, UK
| | - Lea Melki
- R&D Algorithms, Acutus Medical Carlsbad, CA, US
| | - Gregory Slabaugh
- Digital Environment Research Institute, Queen Mary University of London London, UK
| | | | - Caroline H Roney
- School of Engineering and Materials Science, Queen Mary University of London London, UK
| |
Collapse
|
3
|
Abramochkin D, Li B, Zhang H, Kravchuk E, Nesterova T, Glukhov G, Shestak A, Zaklyazminskaya E, Sokolova OS. Novel Gain-of-Function Mutation in the Kv11.1 Channel Found in the Patient with Brugada Syndrome and Mild QTc Shortening. BIOCHEMISTRY. BIOKHIMIIA 2024; 89:543-552. [PMID: 38648771 DOI: 10.1134/s000629792403012x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/06/2023] [Revised: 02/20/2024] [Accepted: 02/27/2024] [Indexed: 04/25/2024]
Abstract
Brugada syndrome (BrS) is an inherited disease characterized by right precordial ST-segment elevation in the right precordial leads on electrocardiograms (ECG), and high risk of life-threatening ventricular arrhythmia and sudden cardiac death (SCD). Mutations in the responsible genes have not been fully characterized in the BrS patients, except for the SCN5A gene. We identified a new genetic variant, c.1189C>T (p.R397C), in the KCNH2 gene in the asymptomatic male proband diagnosed with BrS and mild QTc shortening. We hypothesize that this variant could alter IKr-current and may be causative for the rare non-SCN5A-related form of BrS. To assess its pathogenicity, we performed patch-clamp analysis on IKr reconstituted with this KCNH2 mutation in the Chinese hamster ovary cells and compared the phenotype with the wild type. It appeared that the R397C mutation does not affect the IKr density, but facilitates activation, hampers inactivation of the hERG channels, and increases magnitude of the window current suggesting that the p.R397C is a gain-of-function mutation. In silico modeling demonstrated that this missense mutation potentially leads to the shortening of action potential in the heart.
Collapse
Affiliation(s)
- Denis Abramochkin
- Shenzhen MSU-BIT University, Shenzhen, China.
- Lomonosov Moscow State University, 119234, Moscow, Russia
| | - Bowen Li
- Shenzhen MSU-BIT University, Shenzhen, China.
| | - Han Zhang
- Shenzhen MSU-BIT University, Shenzhen, China.
| | | | - Tatiana Nesterova
- Institute of Immunology and Physiology, Ural Branch of Russian Academy of Sciences, Ekaterinburg, 620049, Russia.
- Institute of Natural Sciences and Mathematics, Ural Federal University, Ekaterinburg, 620075, Russia
| | - Grigory Glukhov
- Shenzhen MSU-BIT University, Shenzhen, China.
- Lomonosov Moscow State University, 119234, Moscow, Russia
| | - Anna Shestak
- Petrovsky National Research Center of Surgery, Moscow, 119991, Russia.
| | | | - Olga S Sokolova
- Shenzhen MSU-BIT University, Shenzhen, China.
- Lomonosov Moscow State University, 119234, Moscow, Russia
| |
Collapse
|
4
|
Llopis-Lorente J, Baroudi S, Koloskoff K, Mora MT, Basset M, Romero L, Benito S, Dayan F, Saiz J, Trenor B. Combining pharmacokinetic and electrophysiological models for early prediction of drug-induced arrhythmogenicity. COMPUTER METHODS AND PROGRAMS IN BIOMEDICINE 2023; 242:107860. [PMID: 37844488 DOI: 10.1016/j.cmpb.2023.107860] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/30/2023] [Revised: 09/28/2023] [Accepted: 10/10/2023] [Indexed: 10/18/2023]
Abstract
BACKGROUND AND OBJECTIVE In silico methods are gaining attention for predicting drug-induced Torsade de Pointes (TdP) in different stages of drug development. However, many computational models tended not to account for inter-individual response variability due to demographic covariates, such as sex, or physiologic covariates, such as renal function, which may be crucial when predicting TdP. This study aims to compare the effects of drugs in male and female populations with normal and impaired renal function using in silico methods. METHODS Pharmacokinetic models considering sex and renal function as covariates were implemented from data published in pharmacokinetic studies. Drug effects were simulated using an electrophysiologically calibrated population of cellular models of 300 males and 300 females. The population of models was built by modifying the endocardial action potential model published by O'Hara et al. (2011) according to the experimentally measured gene expression levels of 12 ion channels. RESULTS Fifteen pharmacokinetic models for CiPA drugs were implemented and validated in this study. Eight pharmacokinetic models included the effect of renal function and four the effect of sex. The mean difference in action potential duration (APD) between male and female populations was 24.9 ms (p<0.05). Our simulations indicated that women with impaired renal function were particularly susceptible to drug-induced arrhythmias, whereas healthy men were less prone to TdP. Differences between patient groups were more pronounced for high TdP-risk drugs. The proposed in silico tool also revealed that individuals with impaired renal function, electrophysiologically simulated with hyperkalemia (extracellular potassium concentration [K+]o = 7 mM) exhibited less pronounced APD prolongation than individuals with normal potassium levels. The pharmacokinetic/electrophysiological framework was used to determine the maximum safe dose of dofetilide in different patient groups. As a proof of concept, 3D simulations were also run for dofetilide obtaining QT prolongation in accordance with previously reported clinical values. CONCLUSIONS This study presents a novel methodology that combines pharmacokinetic and electrophysiological models to incorporate the effects of sex and renal function into in silico drug simulations and highlights their impact on TdP-risk assessment. Furthermore, it may also help inform maximum dose regimens that ensure TdP-related safety in a specific sub-population of patients.
Collapse
Affiliation(s)
- Jordi Llopis-Lorente
- Centro de Investigación e Innovación en Bioingeniería (Ci(2)B), Universitat Politècnica de València, camino de Vera, s/n, 46022, Valencia, Spain
| | | | | | - Maria Teresa Mora
- Centro de Investigación e Innovación en Bioingeniería (Ci(2)B), Universitat Politècnica de València, camino de Vera, s/n, 46022, Valencia, Spain
| | | | - Lucía Romero
- Centro de Investigación e Innovación en Bioingeniería (Ci(2)B), Universitat Politècnica de València, camino de Vera, s/n, 46022, Valencia, Spain
| | | | | | - Javier Saiz
- Centro de Investigación e Innovación en Bioingeniería (Ci(2)B), Universitat Politècnica de València, camino de Vera, s/n, 46022, Valencia, Spain
| | - Beatriz Trenor
- Centro de Investigación e Innovación en Bioingeniería (Ci(2)B), Universitat Politècnica de València, camino de Vera, s/n, 46022, Valencia, Spain.
| |
Collapse
|
5
|
Filgueiras-Rama D, Ramos-Prada A, Cluitmans MJM. Arrhythmogenic vulnerability of reentrant pathways in post-infarct ventricular tachycardia assessed by advanced computational modelling. Europace 2023; 25:euad258. [PMID: 37647101 PMCID: PMC10481246 DOI: 10.1093/europace/euad258] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2023] [Accepted: 08/24/2023] [Indexed: 09/01/2023] Open
Affiliation(s)
- David Filgueiras-Rama
- Centro Nacional de Investigaciones Cardiovasculares (CNIC), Novel Arrhythmogenic Mechanisms Program, Melchor Fernández Almagro, 3, 28029 Madrid, Spain
- Instituto de Investigación Sanitaria del Hospital Clínico San Carlos (IdISSC), Cardiovascular Institute, Profesor Martín Lagos s/n, 28040Madrid, Spain
- Centro de Investigación Biomédica en Red de Enfermedades Cardiovasculares (CIBERCV), Monforte de Lemos 3-5, 28029 Madrid, Spain
| | - Alba Ramos-Prada
- Centro Nacional de Investigaciones Cardiovasculares (CNIC), Novel Arrhythmogenic Mechanisms Program, Melchor Fernández Almagro, 3, 28029 Madrid, Spain
- Fundación Interhospitalaria para la Investigación Cardiovascular (FIC), Madrid, Spain
| | - Matthijs J M Cluitmans
- Cardiovascular Research Institute Maastricht (CARIM), Maastricht University Medical Center, Department of Cardiology, Maastricht, The Netherlands
- Philips Research, Eindhoven, The Netherlands
| |
Collapse
|
6
|
Berg LA, Rocha BM, Oliveira RS, Sebastian R, Rodriguez B, de Queiroz RAB, Cherry EM, Dos Santos RW. Enhanced optimization-based method for the generation of patient-specific models of Purkinje networks. Sci Rep 2023; 13:11788. [PMID: 37479707 PMCID: PMC10362015 DOI: 10.1038/s41598-023-38653-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2023] [Accepted: 07/12/2023] [Indexed: 07/23/2023] Open
Abstract
Cardiac Purkinje networks are a fundamental part of the conduction system and are known to initiate a variety of cardiac arrhythmias. However, patient-specific modeling of Purkinje networks remains a challenge due to their high morphological complexity. This work presents a novel method based on optimization principles for the generation of Purkinje networks that combines geometric and activation accuracy in branch size, bifurcation angles, and Purkinje-ventricular-junction activation times. Three biventricular meshes with increasing levels of complexity are used to evaluate the performance of our approach. Purkinje-tissue coupled monodomain simulations are executed to evaluate the generated networks in a realistic scenario using the most recent Purkinje/ventricular human cellular models and physiological values for the Purkinje-ventricular-junction characteristic delay. The results demonstrate that the new method can generate patient-specific Purkinje networks with controlled morphological metrics and specified local activation times at the Purkinje-ventricular junctions.
Collapse
Affiliation(s)
- Lucas Arantes Berg
- Graduate Program in Computational Modeling, Federal University of Juiz de Fora, Juiz de Fora, Brazil.
- Department of Computer Science, University of Oxford, Oxford, UK.
| | - Bernardo Martins Rocha
- Graduate Program in Computational Modeling, Federal University of Juiz de Fora, Juiz de Fora, Brazil
| | - Rafael Sachetto Oliveira
- Department of Computer Science, Federal University of São João del-Rei, São João del-Rei, Brazil
| | - Rafael Sebastian
- Department of Computer Science, Universitat de Valencia, Valencia, Spain
| | - Blanca Rodriguez
- Department of Computer Science, University of Oxford, Oxford, UK
| | - Rafael Alves Bonfim de Queiroz
- Graduate Program in Computational Modeling, Federal University of Juiz de Fora, Juiz de Fora, Brazil
- Department of Computer Science, Federal University of Ouro Preto, Ouro Preto, Brazil
| | - Elizabeth M Cherry
- School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, GA, USA
| | - Rodrigo Weber Dos Santos
- Graduate Program in Computational Modeling, Federal University of Juiz de Fora, Juiz de Fora, Brazil
| |
Collapse
|
7
|
Dokuchaev A, Chumarnaya T, Bazhutina A, Khamzin S, Lebedeva V, Lyubimtseva T, Zubarev S, Lebedev D, Solovyova O. Combination of personalized computational modeling and machine learning for optimization of left ventricular pacing site in cardiac resynchronization therapy. Front Physiol 2023; 14:1162520. [PMID: 37497440 PMCID: PMC10367108 DOI: 10.3389/fphys.2023.1162520] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2023] [Accepted: 06/26/2023] [Indexed: 07/28/2023] Open
Abstract
Introduction: The 30-50% non-response rate to cardiac resynchronization therapy (CRT) calls for improved patient selection and optimized pacing lead placement. The study aimed to develop a novel technique using patient-specific cardiac models and machine learning (ML) to predict an optimal left ventricular (LV) pacing site (ML-PS) that maximizes the likelihood of LV ejection fraction (LVEF) improvement in a given CRT candidate. To validate the approach, we evaluated whether the distance DPS between the clinical LV pacing site (ref-PS) and ML-PS is associated with improved response rate and magnitude. Materials and methods: We reviewed retrospective data for 57 CRT recipients. A positive response was defined as a more than 10% LVEF improvement. Personalized models of ventricular activation and ECG were created from MRI and CT images. The characteristics of ventricular activation during intrinsic rhythm and biventricular (BiV) pacing with ref-PS were derived from the models and used in combination with clinical data to train supervised ML classifiers. The best logistic regression model classified CRT responders with a high accuracy of 0.77 (ROC AUC = 0.84). The LR classifier, model simulations and Bayesian optimization with Gaussian process regression were combined to identify an optimal ML-PS that maximizes the ML-score of CRT response over the LV surface in each patient. Results: The optimal ML-PS improved the ML-score by 17 ± 14% over the ref-PS. Twenty percent of the non-responders were reclassified as positive at ML-PS. Selection of positive patients with a max ML-score >0.5 demonstrated an improved clinical response rate. The distance DPS was shorter in the responders. The max ML-score and DPS were found to be strong predictors of CRT response (ROC AUC = 0.85). In the group with max ML-score > 0.5 and DPS< 30 mm, the response rate was 83% compared to 14% in the rest of the cohort. LVEF improvement in this group was higher than in the other patients (16 ± 8% vs. 7 ± 8%). Conclusion: A new technique combining clinical data, personalized heart modelling and supervised ML demonstrates the potential for use in clinical practice to assist in optimizing patient selection and predicting optimal LV pacing lead position in HF candidates for CRT.
Collapse
Affiliation(s)
- Arsenii Dokuchaev
- Institute of Immunology and Physiology, Ural Branch of the Russian Academy of Sciences, Ekaterinburg, Russia
| | - Tatiana Chumarnaya
- Institute of Immunology and Physiology, Ural Branch of the Russian Academy of Sciences, Ekaterinburg, Russia
- Laboratory of Mathematical Modeling in Physiology and Medicine Based on Supercomputers, Ural Federal University, Ekaterinburg, Russia
| | - Anastasia Bazhutina
- Institute of Immunology and Physiology, Ural Branch of the Russian Academy of Sciences, Ekaterinburg, Russia
- Laboratory of Mathematical Modeling in Physiology and Medicine Based on Supercomputers, Ural Federal University, Ekaterinburg, Russia
| | - Svyatoslav Khamzin
- Institute of Immunology and Physiology, Ural Branch of the Russian Academy of Sciences, Ekaterinburg, Russia
| | | | - Tamara Lyubimtseva
- Institute of Immunology and Physiology, Ural Branch of the Russian Academy of Sciences, Ekaterinburg, Russia
- Almazov National Medical Research Centre, Saint Petersburg, Russia
| | - Stepan Zubarev
- Institute of Immunology and Physiology, Ural Branch of the Russian Academy of Sciences, Ekaterinburg, Russia
- Almazov National Medical Research Centre, Saint Petersburg, Russia
| | - Dmitry Lebedev
- Institute of Immunology and Physiology, Ural Branch of the Russian Academy of Sciences, Ekaterinburg, Russia
- Almazov National Medical Research Centre, Saint Petersburg, Russia
| | - Olga Solovyova
- Institute of Immunology and Physiology, Ural Branch of the Russian Academy of Sciences, Ekaterinburg, Russia
- Laboratory of Mathematical Modeling in Physiology and Medicine Based on Supercomputers, Ural Federal University, Ekaterinburg, Russia
| |
Collapse
|
8
|
Biasi N, Seghetti P, Mercati M, Tognetti A. A smoothed boundary bidomain model for cardiac simulations in anatomically detailed geometries. PLoS One 2023; 18:e0286577. [PMID: 37294777 PMCID: PMC10256234 DOI: 10.1371/journal.pone.0286577] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2022] [Accepted: 05/18/2023] [Indexed: 06/11/2023] Open
Abstract
This manuscript presents a novel finite difference method to solve cardiac bidomain equations in anatomical models of the heart. The proposed method employs a smoothed boundary approach that represents the boundaries between the heart and the surrounding medium as a spatially diffuse interface of finite thickness. The bidomain boundary conditions are implicitly implemented in the smoothed boundary bidomain equations presented in the manuscript without the need of a structured mesh that explicitly tracks the heart-torso boundaries. We reported some significant examples assessing the method's accuracy using nontrivial test geometries and demonstrating the applicability of the method to complex anatomically detailed human cardiac geometries. In particular, we showed that our approach could be employed to simulate cardiac defibrillation in a human left ventricle comprising fiber architecture. The main advantage of the proposed method is the possibility of implementing bidomain boundary conditions directly on voxel structures, which makes it attractive for three dimensional, patient specific simulations based on medical images. Moreover, given the ease of implementation, we believe that the proposed method could provide an interesting and feasible alternative to finite element methods, and could find application in future cardiac research guiding electrotherapy with computational models.
Collapse
Affiliation(s)
- Niccolò Biasi
- Information Engineering Department, University of Pisa, Pisa, Italy
| | - Paolo Seghetti
- Health Science Interdisciplinary Center, Scuola Superiore Sant’Anna, Pisa, Italy
- National Research Council, Institute of Clinical Physiology, Pisa, Italy
| | - Matteo Mercati
- Information Engineering Department, University of Pisa, Pisa, Italy
| | - Alessandro Tognetti
- Information Engineering Department, University of Pisa, Pisa, Italy
- Research Centre “E. Piaggio”, University of Pisa, Pisa, Italy
| |
Collapse
|
9
|
Que W, Han C, Zhao X, Shi L. An ECG generative model of myocardial infarction. COMPUTER METHODS AND PROGRAMS IN BIOMEDICINE 2022; 225:107062. [PMID: 35994870 DOI: 10.1016/j.cmpb.2022.107062] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/21/2022] [Revised: 08/02/2022] [Accepted: 08/06/2022] [Indexed: 06/15/2023]
Abstract
Background and Objective Computer-aided diagnosis (CAD) of Myocardial Infarction (MI) using machine learning depends on a large amount of clinical Electrocardiogram (ECG) data. Existing infarct ECG databases face the problem of class imbalance. Data augmentation using generative simulation models is a new approach to effectively address this problem. Methods A multiscale ECG generative model was established for ECG data augmentation. In the cellular layer, an ischemic Action Potential (AP) model was established to generate APs in cardiomyocytes with different transmural regions of infraction or different ischemic durations. In the tissue layer, a probability-driven cellular automata excitation propagation model was established to simulate the propagation speed and direction of excitation. An infarct tissue model and a coronary artery model were established to describe the spatiotemporal diversity of MI. A ventricle model, a human torso model, and a computational model of surface ECG based on field source theory were established in the heart-torso layer. Results The model generated pathological 12-lead ECGs of MI with different topography and different extent. When simulating different ventricular wall infarction, the lesions appear in the same leads as the clinical 12-lead ECG. The ST-segment decreases and the T-wave amplitude decreases, similar to the clinical ECG features when simulating subendocardial ischemia. The average fidelity of the 12-lead ECG the model generated is 95.6%, according to the designed DTW-GRA distance algorithm. Conclusions The generative model considers the electrophysiological properties of the natural heart, the pathology of myocardial infarction, and the diversity of clinical ECGs. The model can provide many reliable samples for machine learning of MI.
Collapse
Affiliation(s)
- Wenge Que
- Department of Automation, Tsinghua University, Beijing 100084, China.
| | - Chuang Han
- School of Computer and Communication Engineering, Zhengzhou University of Light Industry, Zhengzhou 450000, China
| | - Xiliang Zhao
- Center for Coronary Artery Disease, Division of Cardiology, Beijing Anzhen Hospital, Capital Medical University, Beijing 100029, China.
| | - Li Shi
- Department of Automation, Tsinghua University, Beijing 100084, China; Beijing National Research Center for Information Science and Technology, Beijing 100084, China.
| |
Collapse
|
10
|
Gillette K, Gsell MAF, Strocchi M, Grandits T, Neic A, Manninger M, Scherr D, Roney CH, Prassl AJ, Augustin CM, Vigmond EJ, Plank G. A personalized real-time virtual model of whole heart electrophysiology. Front Physiol 2022; 13:907190. [PMID: 36213235 PMCID: PMC9539798 DOI: 10.3389/fphys.2022.907190] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2022] [Accepted: 08/15/2022] [Indexed: 01/12/2023] Open
Abstract
Computer models capable of representing the intrinsic personal electrophysiology (EP) of the heart in silico are termed virtual heart technologies. When anatomy and EP are tailored to individual patients within the model, such technologies are promising clinical and industrial tools. Regardless of their vast potential, few virtual technologies simulating the entire organ-scale EP of all four-chambers of the heart have been reported and widespread clinical use is limited due to high computational costs and difficulty in validation. We thus report on the development of a novel virtual technology representing the electrophysiology of all four-chambers of the heart aiming to overcome these limitations. In our previous work, a model of ventricular EP embedded in a torso was constructed from clinical magnetic resonance image (MRI) data and personalized according to the measured 12 lead electrocardiogram (ECG) of a single subject under normal sinus rhythm. This model is then expanded upon to include whole heart EP and a detailed representation of the His-Purkinje system (HPS). To test the capacities of the personalized virtual heart technology to replicate standard clinical morphological ECG features under such conditions, bundle branch blocks within both the right and the left ventricles under two different conduction velocity settings are modeled alongside sinus rhythm. To ensure clinical viability, model generation was completely automated and simulations were performed using an efficient real-time cardiac EP simulator. Close correspondence between the measured and simulated 12 lead ECG was observed under normal sinus conditions and all simulated bundle branch blocks manifested relevant clinical morphological features.
Collapse
Affiliation(s)
- Karli Gillette
- Gottfried Schatz Research Center—Biophysics, Medical University of Graz, Graz, Austria
- BioTechMed-Graz, Graz, Austria
| | - Matthias A. F. Gsell
- Gottfried Schatz Research Center—Biophysics, Medical University of Graz, Graz, Austria
- NAWI Graz, Institute of Mathematics and Scientific Computing, University of Graz, Graz, Austria
| | | | - Thomas Grandits
- Gottfried Schatz Research Center—Biophysics, Medical University of Graz, Graz, Austria
- NAWI Graz, Institute of Mathematics and Scientific Computing, University of Graz, Graz, Austria
| | | | - Martin Manninger
- Division of Cardiology, Department of Internal Medicine, Medical University of Graz, Graz, Austria
| | - Daniel Scherr
- Division of Cardiology, Department of Internal Medicine, Medical University of Graz, Graz, Austria
| | | | - Anton J. Prassl
- Gottfried Schatz Research Center—Biophysics, Medical University of Graz, Graz, Austria
- BioTechMed-Graz, Graz, Austria
| | - Christoph M. Augustin
- Gottfried Schatz Research Center—Biophysics, Medical University of Graz, Graz, Austria
- BioTechMed-Graz, Graz, Austria
| | | | - Gernot Plank
- Gottfried Schatz Research Center—Biophysics, Medical University of Graz, Graz, Austria
- BioTechMed-Graz, Graz, Austria
- *Correspondence: Gernot Plank,
| |
Collapse
|
11
|
Abstract
The global burden caused by cardiovascular disease is substantial, with heart disease representing the most common cause of death around the world. There remains a need to develop better mechanistic models of cardiac function in order to combat this health concern. Heart rhythm disorders, or arrhythmias, are one particular type of disease which has been amenable to quantitative investigation. Here we review the application of quantitative methodologies to explore dynamical questions pertaining to arrhythmias. We begin by describing single-cell models of cardiac myocytes, from which two and three dimensional models can be constructed. Special focus is placed on results relating to pattern formation across these spatially-distributed systems, especially the formation of spiral waves of activation. Next, we discuss mechanisms which can lead to the initiation of arrhythmias, focusing on the dynamical state of spatially discordant alternans, and outline proposed mechanisms perpetuating arrhythmias such as fibrillation. We then review experimental and clinical results related to the spatio-temporal mapping of heart rhythm disorders. Finally, we describe treatment options for heart rhythm disorders and demonstrate how statistical physics tools can provide insights into the dynamics of heart rhythm disorders.
Collapse
Affiliation(s)
- Wouter-Jan Rappel
- Department of Physics, University of California San Diego, La Jolla, CA 92037
| |
Collapse
|
12
|
de Lepper AGW, Buck CMA, van 't Veer M, Huberts W, van de Vosse FN, Dekker LRC. From evidence-based medicine to digital twin technology for predicting ventricular tachycardia in ischaemic cardiomyopathy. JOURNAL OF THE ROYAL SOCIETY, INTERFACE 2022; 19:20220317. [PMID: 36128708 DOI: 10.1098/rsif.2022.0317] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
Survivors of myocardial infarction are at risk of life-threatening ventricular tachycardias (VTs) later in their lives. Current guidelines for implantable cardioverter defibrillators (ICDs) implantation to prevent VT-related sudden cardiac death is solely based on symptoms and left ventricular ejection fraction. Catheter ablation of scar-related VTs is performed following ICD therapy, reducing VTs, painful shocks, anxiety, depression and worsening heart failure. We postulate that better prediction of the occurrence and circuit of VT, will improve identification of patients at risk for VT and boost preventive ablation, reducing mortality and morbidity. For this purpose, multiple time-evolving aspects of the underlying pathophysiology, including the anatomical substrate, triggers and modulators, should be part of VT prediction models. We envision digital twins as a solution combining clinical expertise with three prediction approaches: evidence-based medicine (clinical practice), data-driven models (data science) and mechanistic models (biomedical engineering). This paper aims to create a mutual understanding between experts in the different fields by providing a comprehensive description of the clinical problem and the three approaches in an understandable manner, leveraging future collaborations and technological innovations for clinical decision support. Moreover, it defines open challenges and gains for digital twin solutions and discusses the potential of hybrid modelling.
Collapse
Affiliation(s)
| | - Carlijn M A Buck
- Department of Biomedical Engineering, Eindhoven University of Technology, Eindhoven, The Netherlands
| | - Marcel van 't Veer
- Department of Cardiology, Catharina Hospital, Eindhoven, The Netherlands.,Department of Biomedical Engineering, Eindhoven University of Technology, Eindhoven, The Netherlands
| | - Wouter Huberts
- Department of Biomedical Engineering, Eindhoven University of Technology, Eindhoven, The Netherlands.,Department of Biomedical Engineering, CARIM School for Cardiovascular Diseases, Maastricht University, Maastricht, The Netherlands
| | - Frans N van de Vosse
- Department of Biomedical Engineering, Eindhoven University of Technology, Eindhoven, The Netherlands
| | - Lukas R C Dekker
- Department of Cardiology, Catharina Hospital, Eindhoven, The Netherlands.,Department of Biomedical Engineering, Eindhoven University of Technology, Eindhoven, The Netherlands
| |
Collapse
|
13
|
Doste R, Lozano M, Jimenez-Perez G, Mont L, Berruezo A, Penela D, Camara O, Sebastian R. Training machine learning models with synthetic data improves the prediction of ventricular origin in outflow tract ventricular arrhythmias. Front Physiol 2022; 13:909372. [PMID: 36035489 PMCID: PMC9412034 DOI: 10.3389/fphys.2022.909372] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2022] [Accepted: 07/07/2022] [Indexed: 11/13/2022] Open
Abstract
In order to determine the site of origin (SOO) in outflow tract ventricular arrhythmias (OTVAs) before an ablation procedure, several algorithms based on manual identification of electrocardiogram (ECG) features, have been developed. However, the reported accuracy decreases when tested with different datasets. Machine learning algorithms can automatize the process and improve generalization, but their performance is hampered by the lack of large enough OTVA databases. We propose the use of detailed electrophysiological simulations of OTVAs to train a machine learning classification model to predict the ventricular origin of the SOO of ectopic beats. We generated a synthetic database of 12-lead ECGs (2,496 signals) by running multiple simulations from the most typical OTVA SOO in 16 patient-specific geometries. Two types of input data were considered in the classification, raw and feature ECG signals. From the simulated raw 12-lead ECG, we analyzed the contribution of each lead in the predictions, keeping the best ones for the training process. For feature-based analysis, we used entropy-based methods to rank the obtained features. A cross-validation process was included to evaluate the machine learning model. Following, two clinical OTVA databases from different hospitals, including ECGs from 365 patients, were used as test-sets to assess the generalization of the proposed approach. The results show that V2 was the best lead for classification. Prediction of the SOO in OTVA, using both raw signals or features for classification, presented high accuracy values (>0.96). Generalization of the network trained on simulated data was good for both patient datasets (accuracy of 0.86 and 0.84, respectively) and presented better values than using exclusively real ECGs for classification (accuracy of 0.84 and 0.76 for each dataset). The use of simulated ECG data for training machine learning-based classification algorithms is critical to obtain good SOO predictions in OTVA compared to real data alone. The fast implementation and generalization of the proposed methodology may contribute towards its application to a clinical routine.
Collapse
Affiliation(s)
- Ruben Doste
- Department of Computer Science, University of Oxford, Oxford, United Kingdom
- *Correspondence: Ruben Doste,
| | - Miguel Lozano
- Computational Multiscale Simulation Lab (CoMMLab), Department of Computer Science, Universitat de Valencia, Valencia, Spain
| | - Guillermo Jimenez-Perez
- Physense, BCN Medtech, Department of Information and Communication Technologies, Universitat Pompeu Fabra, Barcelona, Spain
| | - Lluis Mont
- Arrhythmia Section, Cardiology Department, Cardiovascular Clinical Institute, Hospital Clínic, Universitat de Barcelona - IDIBAPS, Barcelona, Spain
| | - Antonio Berruezo
- Cardiology Department, Heart Institute, Teknon Medical Center, Barcelona, Spain
| | - Diego Penela
- Cardiology Department, Heart Institute, Teknon Medical Center, Barcelona, Spain
| | - Oscar Camara
- Physense, BCN Medtech, Department of Information and Communication Technologies, Universitat Pompeu Fabra, Barcelona, Spain
| | - Rafael Sebastian
- Computational Multiscale Simulation Lab (CoMMLab), Department of Computer Science, Universitat de Valencia, Valencia, Spain
| |
Collapse
|
14
|
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
|
15
|
Meshless Electrophysiological Modeling of Cardiac Resynchronization Therapy—Benchmark Analysis with Finite-Element Methods in Experimental Data. APPLIED SCIENCES-BASEL 2022. [DOI: 10.3390/app12136438] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/10/2022]
Abstract
Computational models of cardiac electrophysiology are promising tools for reducing the rates of non-response patients suitable for cardiac resynchronization therapy (CRT) by optimizing electrode placement. The majority of computational models in the literature are mesh-based, primarily using the finite element method (FEM). The generation of patient-specific cardiac meshes has traditionally been a tedious task requiring manual intervention and hindering the modeling of a large number of cases. Meshless models can be a valid alternative due to their mesh quality independence. The organization of challenges such as the CRT-EPiggy19, providing unique experimental data as open access, enables benchmarking analysis of different cardiac computational modeling solutions with quantitative metrics. We present a benchmark analysis of a meshless-based method with finite-element methods for the prediction of cardiac electrical patterns in CRT, based on a subset of the CRT-EPiggy19 dataset. A data assimilation strategy was designed to personalize the most relevant parameters of the electrophysiological simulations and identify the optimal CRT lead configuration. The simulation results obtained with the meshless model were equivalent to FEM, with the most relevant aspect for accurate CRT predictions being the parameter personalization strategy (e.g., regional conduction velocity distribution, including the Purkinje system and CRT lead distribution).
Collapse
|
16
|
Molinari L, Zaltieri M, Massaroni C, Filippi S, Gizzi A, Schena E. Multiscale and Multiphysics Modeling of Anisotropic Cardiac RFCA: Experimental-Based Model Calibration via Multi-Point Temperature Measurements. Front Physiol 2022; 13:845896. [PMID: 35514332 PMCID: PMC9062295 DOI: 10.3389/fphys.2022.845896] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/30/2021] [Accepted: 03/22/2022] [Indexed: 11/13/2022] Open
Abstract
Radiofrequency catheter ablation (RFCA) is the mainstream treatment for drug-refractory cardiac fibrillation. Multiple studies demonstrated that incorrect dosage of radiofrequency energy to the myocardium could lead to uncontrolled tissue damage or treatment failure, with the consequent need for unplanned reoperations. Monitoring tissue temperature during thermal therapy and predicting the extent of lesions may improve treatment efficacy. Cardiac computational modeling represents a viable tool for identifying optimal RFCA settings, though predictability issues still limit a widespread usage of such a technology in clinical scenarios. We aim to fill this gap by assessing the influence of the intrinsic myocardial microstructure on the thermo-electric behavior at the tissue level. By performing multi-point temperature measurements on ex-vivo swine cardiac tissue samples, the experimental characterization of myocardial thermal anisotropy allowed us to assemble a fine-tuned thermo-electric material model of the cardiac tissue. We implemented a multiphysics and multiscale computational framework, encompassing thermo-electric anisotropic conduction, phase-lagging for heat transfer, and a three-state dynamical system for cellular death and lesion estimation. Our analysis resulted in a remarkable agreement between ex-vivo measurements and numerical results. Accordingly, we identified myocardium anisotropy as the driving effect on the outcomes of hyperthermic treatments. Furthermore, we characterized the complex nonlinear couplings regulating tissue behavior during RFCA, discussing model calibration, limitations, and perspectives.
Collapse
Affiliation(s)
- Leonardo Molinari
- Department of Mathematics and Computer Science, Emory University, Atlanta, GA, United States
| | - Martina Zaltieri
- Laboratory of Measurement and Biomedical Instrumentation, Department of Engineering, University of Rome Campus Bio-Medico, Rome, Italy
| | - Carlo Massaroni
- Laboratory of Measurement and Biomedical Instrumentation, Department of Engineering, University of Rome Campus Bio-Medico, Rome, Italy
| | - Simonetta Filippi
- Nonlinear Physics and Mathematical Modeling Lab, Department of Engineering, University of Rome Campus Bio-Medico, Rome, Italy
| | - Alessio Gizzi
- Nonlinear Physics and Mathematical Modeling Lab, Department of Engineering, University of Rome Campus Bio-Medico, Rome, Italy
| | - Emiliano Schena
- Laboratory of Measurement and Biomedical Instrumentation, Department of Engineering, University of Rome Campus Bio-Medico, Rome, Italy
| |
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
|
Fu Z, Zhang J, Luo R, Sun Y, Deng D, Xia L. TF-Unet:An automatic cardiac MRI image segmentation method. MATHEMATICAL BIOSCIENCES AND ENGINEERING : MBE 2022; 19:5207-5222. [PMID: 35430861 DOI: 10.3934/mbe.2022244] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Personalized heart models are widely used to study the mechanisms of cardiac arrhythmias and have been used to guide clinical ablation of different types of arrhythmias in recent years. MRI images are now mostly used for model building. In cardiac modeling studies, the degree of segmentation of the heart image determines the success of subsequent 3D reconstructions. Therefore, a fully automated segmentation is needed. In this paper, we combine U-Net and Transformer as an alternative approach to perform powerful and fully automated segmentation of medical images. On the one hand, we use convolutional neural networks for feature extraction and spatial encoding of inputs to fully exploit the advantages of convolution in detail grasping; on the other hand, we use Transformer to add remote dependencies to high-level features and model features at different scales to fully exploit the advantages of Transformer. The results show that, the average dice coefficients for ACDC and Synapse datasets are 91.72 and 85.46%, respectively, and compared with Swin-Unet, the segmentation accuracy are improved by 1.72% for ACDC dataset and 6.33% for Synapse dataset.
Collapse
Affiliation(s)
- Zhenyin Fu
- Key Laboratory for Biomedical Engineering of Ministry of Education, Institute of Biomedical Engineering, Zhejiang University, Hangzhou 310027, China
| | - Jin Zhang
- Key Laboratory for Biomedical Engineering of Ministry of Education, Institute of Biomedical Engineering, Zhejiang University, Hangzhou 310027, China
| | - Ruyi Luo
- Hangzhou Science and Technology Information Institute, Hangzhou 310026, China
| | - Yutong Sun
- School of Biomedical Engineering, Dalian University of Technology, Dalian 116024, China
| | - Dongdong Deng
- School of Biomedical Engineering, Dalian University of Technology, Dalian 116024, China
| | - Ling Xia
- Key Laboratory for Biomedical Engineering of Ministry of Education, Institute of Biomedical Engineering, Zhejiang University, Hangzhou 310027, China
- Research Center for Healthcare Data Science, Zhejiang Lab, Hangzhou 310026, China
| |
Collapse
|
19
|
Khamzin S, Dokuchaev A, Bazhutina A, Chumarnaya T, Zubarev S, Lyubimtseva T, Lebedeva V, Lebedev D, Gurev V, Solovyova O. Machine Learning Prediction of Cardiac Resynchronisation Therapy Response From Combination of Clinical and Model-Driven Data. Front Physiol 2022; 12:753282. [PMID: 34970154 PMCID: PMC8712879 DOI: 10.3389/fphys.2021.753282] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2021] [Accepted: 11/22/2021] [Indexed: 11/17/2022] Open
Abstract
Background: Up to 30–50% of chronic heart failure patients who underwent cardiac resynchronization therapy (CRT) do not respond to the treatment. Therefore, patient stratification for CRT and optimization of CRT device settings remain a challenge. Objective: The main goal of our study is to develop a predictive model of CRT outcome using a combination of clinical data recorded in patients before CRT and simulations of the response to biventricular (BiV) pacing in personalized computational models of the cardiac electrophysiology. Materials and Methods: Retrospective data from 57 patients who underwent CRT device implantation was utilized. Positive response to CRT was defined by a 10% increase in the left ventricular ejection fraction in a year after implantation. For each patient, an anatomical model of the heart and torso was reconstructed from MRI and CT images and tailored to ECG recorded in the participant. The models were used to compute ventricular activation time, ECG duration and electrical dyssynchrony indices during intrinsic rhythm and BiV pacing from the sites of implanted leads. For building a predictive model of CRT response, we used clinical data recorded before CRT device implantation together with model-derived biomarkers of ventricular excitation in the left bundle branch block mode of activation and under BiV stimulation. Several Machine Learning (ML) classifiers and feature selection algorithms were tested on the hybrid dataset, and the quality of predictors was assessed using the area under receiver operating curve (ROC AUC). The classifiers on the hybrid data were compared with ML models built on clinical data only. Results: The best ML classifier utilizing a hybrid set of clinical and model-driven data demonstrated ROC AUC of 0.82, an accuracy of 0.82, sensitivity of 0.85, and specificity of 0.78, improving quality over that of ML predictors built on clinical data from much larger datasets by more than 0.1. Distance from the LV pacing site to the post-infarction zone and ventricular activation characteristics under BiV pacing were shown as the most relevant model-driven features for CRT response classification. Conclusion: Our results suggest that combination of clinical and model-driven data increases the accuracy of classification models for CRT outcomes.
Collapse
Affiliation(s)
- Svyatoslav Khamzin
- Institute of Immunology and Physiology Ural Branch of the Russian Academy of Sciences, Yekaterinburg, Russia
| | - Arsenii Dokuchaev
- Institute of Immunology and Physiology Ural Branch of the Russian Academy of Sciences, Yekaterinburg, Russia
| | - Anastasia Bazhutina
- Institute of Immunology and Physiology Ural Branch of the Russian Academy of Sciences, Yekaterinburg, Russia.,Ural Federal University, Yekaterinburg, Russia
| | - Tatiana Chumarnaya
- Institute of Immunology and Physiology Ural Branch of the Russian Academy of Sciences, Yekaterinburg, Russia
| | - Stepan Zubarev
- Almazov National Medical Research Centre, Saint Petersburg, Russia
| | | | | | - Dmitry Lebedev
- Almazov National Medical Research Centre, Saint Petersburg, Russia
| | | | - Olga Solovyova
- Institute of Immunology and Physiology Ural Branch of the Russian Academy of Sciences, Yekaterinburg, Russia.,Ural Federal University, Yekaterinburg, Russia
| |
Collapse
|
20
|
Tong L, Zhao C, Fu Z, Dong R, Wu Z, Wang Z, Zhang N, Wang X, Cao B, Sun Y, Zheng D, Xia L, Deng D. Preliminary Study: Learning the Impact of Simulation Time on Reentry Location and Morphology Induced by Personalized Cardiac Modeling. Front Physiol 2021; 12:733500. [PMID: 35002750 PMCID: PMC8739986 DOI: 10.3389/fphys.2021.733500] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2021] [Accepted: 11/30/2021] [Indexed: 11/13/2022] Open
Abstract
Personalized cardiac modeling is widely used for studying the mechanisms of cardiac arrythmias. Due to the high demanding of computational resource of modeling, the arrhythmias induced in the models are usually simulated for just a few seconds. In clinic, it is common that arrhythmias last for more than several minutes and the morphologies of reentries are not always stable, so it is not clear that whether the simulation of arrythmias for just a few seconds is long enough to match the arrhythmias detected in patients. This study aimed to observe how long simulation of the induced arrhythmias in the personalized cardiac models is sufficient to match the arrhythmias detected in patients. A total of 5 contrast enhanced MRI datasets of patient hearts with myocardial infarction were used in this study. Then, a classification method based on Gaussian mixture model was used to detect the infarct tissue. For each reentry, 3 s and 10 s were simulated. The characteristics of each reentry simulated for different duration were studied. Reentries were induced in all 5 ventricular models and sustained reentries were induced at 39 stimulation sites in the model. By analyzing the simulation results, we found that 41% of the sustained reentries in the 3 s simulation group terminated in the longer simulation groups (10 s). The second finding in our simulation was that only 23.1% of the sustained reentries in the 3 s simulation did not change location and morphology in the extended 10 s simulation. The third finding was that 35.9% reentries were stable in the 3 s simulation and should be extended for the simulation time. The fourth finding was that the simulation results in 10 s simulation matched better with the clinical measurements than the 3 s simulation. It was shown that 10 s simulation was sufficient to make simulation results stable. The findings of this study not only improve the simulation accuracy, but also reduce the unnecessary simulation time to achieve the optimal use of computer resources to improve the simulation efficiency and shorten the simulation time to meet the time node requirements of clinical operation on patients.
Collapse
Affiliation(s)
- Lv Tong
- School of Biomedical Engineering, Dalian University of Technology, Dalian, China
| | - Caiming Zhao
- Department of Cardiology, The First Affiliated Hospital of Soochow University, Suzhou, China
| | - Zhenyin Fu
- College of Biomedical Engineering and Instrument Science, Zhejiang University, Hangzhou, China
| | - Ruiqing Dong
- Department of Cardiology, Dushu Lake Hospital Affiliated to Soochow University, Suzhou, China
| | - Zhenghong Wu
- College of Biomedical Engineering and Instrument Science, Zhejiang University, Hangzhou, China
| | - Zefeng Wang
- Department of Cardiology, Beijing Anzhen Hospital Affiliated to Capital Medical University, Beijing, China
| | - Nan Zhang
- Department of Radiology, Beijing Anzhen Hospital Affiliated to Capital Medical University, Beijing, China
| | - Xinlu Wang
- Department of Cardiology, Beijing Anzhen Hospital Affiliated to Capital Medical University, Beijing, China
| | - Boyang Cao
- School of Biomedical Engineering, Dalian University of Technology, Dalian, China
| | - Yutong Sun
- School of Biomedical Engineering, Dalian University of Technology, Dalian, China
| | - Dingchang Zheng
- Research Centre for Intelligent Healthcare, Faculty of Health and Life Science, Coventry University, Coventry, United Kingdom
| | - Ling Xia
- College of Biomedical Engineering and Instrument Science, Zhejiang University, Hangzhou, China
| | - Dongdong Deng
- School of Biomedical Engineering, Dalian University of Technology, Dalian, China
| |
Collapse
|
21
|
Analysis of vulnerability to reentry in acute myocardial ischemia using a realistic human heart model. Comput Biol Med 2021; 141:105038. [PMID: 34836624 DOI: 10.1016/j.compbiomed.2021.105038] [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] [Received: 07/31/2021] [Revised: 10/25/2021] [Accepted: 11/12/2021] [Indexed: 11/21/2022]
Abstract
Electrophysiological alterations of the myocardium caused by acute ischemia constitute a pro-arrhythmic substrate for the generation of potentially lethal arrhythmias. Experimental evidence has shown that the main components of acute ischemia that induce these electrophysiological alterations are hyperkalemia, hypoxia (or anoxia in complete artery occlusion), and acidosis. However, the influence of each ischemic component on the likelihood of reentry is not completely established. Moreover, the role of the His-Purkinje system (HPS) in the initiation and maintenance of arrhythmias is not completely understood. In the present work, we investigate how the three components of ischemia affect the vulnerable window (VW) for reentry using computational simulations. In addition, we analyze the role of the HPS on arrhythmogenesis. A 3D biventricular/torso human model that includes a realistic geometry of the central and border ischemic zones with one of the most electrophysiologically detailed model of ischemia to date, as well as a realistic cardiac conduction system, were used to assess the VW for reentry. Four scenarios of ischemic severity corresponding to different minutes after coronary artery occlusion were simulated. Our results suggest that ischemic severity plays an important role in the generation of reentries. Indeed, this is the first 3D simulation study to show that ventricular arrhythmias could be generated under moderate ischemic conditions, but not in mild and severe ischemia. Moreover, our results show that anoxia is the ischemic component with the most significant effect on the width of the VW. Thus, a change in the level of anoxia from moderate to severe leads to a greater increment in the VW (40 ms), in comparison with the increment of 20 ms and 35 ms produced by the individual change in the level of hyperkalemia and acidosis, respectively. Finally, the HPS was a necessary element for the generation of approximately 17% of reentries obtained. The retrograde conduction from the myocardium to HPS in the ischemic region, conduction blocks in discrete sections of the HPS, and the degree of ischemia affecting Purkinje cells, are suggested as mechanisms that favor the generation of ventricular arrhythmias.
Collapse
|
22
|
Sung E, Etoz S, Zhang Y, Trayanova NA. Whole-heart ventricular arrhythmia modeling moving forward: Mechanistic insights and translational applications. BIOPHYSICS REVIEWS 2021; 2:031304. [PMID: 36281224 PMCID: PMC9588428 DOI: 10.1063/5.0058050] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/19/2023]
Abstract
Ventricular arrhythmias are the primary cause of sudden cardiac death and one of the leading causes of mortality worldwide. Whole-heart computational modeling offers a unique approach for studying ventricular arrhythmias, offering vast potential for developing both a mechanistic understanding of ventricular arrhythmias and clinical applications for treatment. In this review, the fundamentals of whole-heart ventricular modeling and current methods of personalizing models using clinical data are presented. From this foundation, the authors summarize recent advances in whole-heart ventricular arrhythmia modeling. Efforts in gaining mechanistic insights into ventricular arrhythmias are discussed, in addition to other applications of models such as the assessment of novel therapeutics. The review emphasizes the unique benefits of computational modeling that allow for insights that are not obtainable by contemporary experimental or clinical means. Additionally, the clinical impact of modeling is explored, demonstrating how patient care is influenced by the information gained from ventricular arrhythmia models. The authors conclude with future perspectives about the direction of whole-heart ventricular arrhythmia modeling, outlining how advances in neural network methodologies hold the potential to reduce computational expense and permit for efficient whole-heart modeling.
Collapse
Affiliation(s)
- Eric Sung
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, Maryland 21218, USA
- Alliance for Cardiovascular Diagnostic and Treatment Innovation, Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Sevde Etoz
- Alliance for Cardiovascular Diagnostic and Treatment Innovation, Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Yingnan Zhang
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, Maryland 21218, USA
- Alliance for Cardiovascular Diagnostic and Treatment Innovation, Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Natalia A. Trayanova
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, Maryland 21218, USA
- Alliance for Cardiovascular Diagnostic and Treatment Innovation, Johns Hopkins University, Baltimore, Maryland 21218, USA
- Author to whom correspondence should be addressed: . Tel.: 410-516-4375
| |
Collapse
|
23
|
Romero P, Lozano M, Martínez-Gil F, Serra D, Sebastián R, Lamata P, García-Fernández I. Clinically-Driven Virtual Patient Cohorts Generation: An Application to Aorta. Front Physiol 2021; 12:713118. [PMID: 34539438 PMCID: PMC8440937 DOI: 10.3389/fphys.2021.713118] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2021] [Accepted: 08/03/2021] [Indexed: 12/20/2022] Open
Abstract
The combination of machine learning methods together with computational modeling and simulation of the cardiovascular system brings the possibility of obtaining very valuable information about new therapies or clinical devices through in-silico experiments. However, the application of machine learning methods demands access to large cohorts of patients. As an alternative to medical data acquisition and processing, which often requires some degree of manual intervention, the generation of virtual cohorts made of synthetic patients can be automated. However, the generation of a synthetic sample can still be computationally demanding to guarantee that it is clinically meaningful and that it reflects enough inter-patient variability. This paper addresses the problem of generating virtual patient cohorts of thoracic aorta geometries that can be used for in-silico trials. In particular, we focus on the problem of generating a cohort of patients that meet a particular clinical criterion, regardless the access to a reference sample of that phenotype. We formalize the problem of clinically-driven sampling and assess several sampling strategies with two goals, sampling efficiency, i.e., that the generated individuals actually belong to the target population, and that the statistical properties of the cohort can be controlled. Our results show that generative adversarial networks can produce reliable, clinically-driven cohorts of thoracic aortas with good efficiency. Moreover, non-linear predictors can serve as an efficient alternative to the sometimes expensive evaluation of anatomical or functional parameters of the organ of interest.
Collapse
Affiliation(s)
- Pau Romero
- Computational Multiscale Simulation Lab, Department of Computer Science, Universitat de Valencia, Valencia, Spain
| | - Miguel Lozano
- Computational Multiscale Simulation Lab, Department of Computer Science, Universitat de Valencia, Valencia, Spain
| | - Francisco Martínez-Gil
- Computational Multiscale Simulation Lab, Department of Computer Science, Universitat de Valencia, Valencia, Spain
| | - Dolors Serra
- Computational Multiscale Simulation Lab, Department of Computer Science, Universitat de Valencia, Valencia, Spain
| | - Rafael Sebastián
- Computational Multiscale Simulation Lab, Department of Computer Science, Universitat de Valencia, Valencia, Spain
| | - Pablo Lamata
- Department of Biomedical Engineering, School of Biomedical Engineering and Imaging Sciences, Kings College London, London, United Kingdom
| | - Ignacio García-Fernández
- Computational Multiscale Simulation Lab, Department of Computer Science, Universitat de Valencia, Valencia, Spain
| |
Collapse
|
24
|
Barber F, Langfield P, Lozano M, Garcia-Fernandez I, Duchateau J, Hocini M, Haissaguerre M, Vigmond E, Sebastian R. Estimation of Personalized Minimal Purkinje Systems From Human Electro-Anatomical Maps. IEEE TRANSACTIONS ON MEDICAL IMAGING 2021; 40:2182-2194. [PMID: 33856987 DOI: 10.1109/tmi.2021.3073499] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
The Purkinje system is a heart structure responsible for transmitting electrical impulses through the ventricles in a fast and coordinated way to trigger mechanical contraction. Estimating a patient-specific compatible Purkinje Network from an electro-anatomical map is a challenging task, that could help to improve models for electrophysiology simulations or provide aid in therapy planning, such as radiofrequency ablation. In this study, we present a methodology to inversely estimate a Purkinje network from a patient's electro-anatomical map. First, we carry out a simulation study to assess the accuracy of the method for different synthetic Purkinje network morphologies and myocardial junction densities. Second, we estimate the Purkinje network from a set of 28 electro-anatomical maps from patients, obtaining an optimal conduction velocity in the Purkinje network of 1.95 ± 0.25 m/s, together with the location of their Purkinje-myocardial junctions, and Purkinje network structure. Our results showed an average local activation time error of 6.8±2.2 ms in the endocardium. Finally, using the personalized Purkinje network, we obtained correlations higher than 0.85 between simulated and clinical 12-lead ECGs.
Collapse
|
25
|
Abstract
Computer modeling of the electrophysiology of the heart has undergone significant progress. A healthy heart can be modeled starting from the ion channels via the spread of a depolarization wave on a realistic geometry of the human heart up to the potentials on the body surface and the ECG. Research is advancing regarding modeling diseases of the heart. This article reviews progress in calculating and analyzing the corresponding electrocardiogram (ECG) from simulated depolarization and repolarization waves. First, we describe modeling of the P-wave, the QRS complex and the T-wave of a healthy heart. Then, both the modeling and the corresponding ECGs of several important diseases and arrhythmias are delineated: ischemia and infarction, ectopic beats and extrasystoles, ventricular tachycardia, bundle branch blocks, atrial tachycardia, flutter and fibrillation, genetic diseases and channelopathies, imbalance of electrolytes and drug-induced changes. Finally, we outline the potential impact of computer modeling on ECG interpretation. Computer modeling can contribute to a better comprehension of the relation between features in the ECG and the underlying cardiac condition and disease. It can pave the way for a quantitative analysis of the ECG and can support the cardiologist in identifying events or non-invasively localizing diseased areas. Finally, it can deliver very large databases of reliably labeled ECGs as training data for machine learning.
Collapse
|
26
|
Fedele M, Quarteroni A. Polygonal surface processing and mesh generation tools for the numerical simulation of the cardiac function. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2021; 37:e3435. [PMID: 33415829 PMCID: PMC8244076 DOI: 10.1002/cnm.3435] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/02/2020] [Revised: 01/01/2021] [Accepted: 01/02/2021] [Indexed: 06/05/2023]
Abstract
In order to simulate the cardiac function for a patient-specific geometry, the generation of the computational mesh is crucially important. In practice, the input is typically a set of unprocessed polygonal surfaces coming either from a template geometry or from medical images. These surfaces need ad-hoc processing to be suitable for a volumetric mesh generation. In this work we propose a set of new algorithms and tools aiming to facilitate the mesh generation process. In particular, we focus on different aspects of a cardiac mesh generation pipeline: (1) specific polygonal surface processing for cardiac geometries, like connection of different heart chambers or segmentation outputs; (2) generation of accurate boundary tags; (3) definition of mesh-size functions dependent on relevant geometric quantities; (4) processing and connecting together several volumetric meshes. The new algorithms-implemented in the open-source software vmtk-can be combined with each other allowing the creation of personalized pipelines, that can be optimized for each cardiac geometry or for each aspect of the cardiac function to be modeled. Thanks to these features, the proposed tools can significantly speed-up the mesh generation process for a large range of cardiac applications, from single-chamber single-physics simulations to multi-chambers multi-physics simulations. We detail all the proposed algorithms motivating them in the cardiac context and we highlight their flexibility by showing different examples of cardiac mesh generation pipelines.
Collapse
Affiliation(s)
- Marco Fedele
- MOX, Department of MathematicsPolitecnico di MilanoMilanItaly
| | - Alfio Quarteroni
- MOX, Department of MathematicsPolitecnico di MilanoMilanItaly
- Institute of MathematicsÉcole Polytechnique Fédérale de LausanneLausanneSwitzerland
| |
Collapse
|
27
|
Jeong DU, Lim KM. Prediction of Cardiac Mechanical Performance From Electrical Features During Ventricular Tachyarrhythmia Simulation Using Machine Learning Algorithms. Front Physiol 2020; 11:591681. [PMID: 33329041 PMCID: PMC7732497 DOI: 10.3389/fphys.2020.591681] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2020] [Accepted: 10/28/2020] [Indexed: 11/13/2022] Open
Abstract
In ventricular tachyarrhythmia, electrical instability features including action potential duration, dominant frequency, phase singularity, and filaments are associated with mechanical contractility. However, there are insufficient studies on estimated mechanical contractility based on electrical features during ventricular tachyarrhythmia using a stochastic model. In this study, we predicted cardiac mechanical performance from features of electrical instability during ventricular tachyarrhythmia simulation using machine learning algorithms, including support vector regression (SVR) and artificial neural network (ANN) models. We performed an electromechanical tachyarrhythmia simulation and extracted 12 electrical instability features and two mechanical properties, including stroke volume and the amplitude of myocardial tension (ampTens). We compared predictive performance according to kernel types of the SVR model and the number of hidden layers of the ANN model. In the SVR model, the prediction accuracies of stroke volume and ampTens were the highest when using the polynomial kernel and linear kernel, respectively. The predictive performance of the ANN model was better than that of the SVR model. The prediction accuracies were the highest when the ANN model consisted of three hidden layers. Accordingly, we propose the ANN model with three hidden layers as an optimal model for predicting cardiac mechanical contractility in ventricular tachyarrhythmia. The results of this study are expected to be used to indirectly estimate the hemodynamic response from the electrical cardiac map measured by the optical mapping system during cardiac surgery, as well as cardiac contractility under normal sinus rhythm conditions.
Collapse
Affiliation(s)
- Da Un Jeong
- Computational Medicine Lab, Department of IT Convergence Engineering, Kumoh National Institute of Technology, Gumi, South Korea
| | - Ki Moo Lim
- Computational Medicine Lab, Department of IT Convergence Engineering, Kumoh National Institute of Technology, Gumi, South Korea.,Computational Medicine Lab, Department of Medical IT Convergence Engineering, Kumoh National Institute of Technology, Gumi, South Korea
| |
Collapse
|
28
|
Wan Ab Naim WN, Mokhtarudin MJM, Lim E, Chan BT, Ahmad Bakir A, Nik Mohamed NA. The study of border zone formation in ischemic heart using electro-chemical coupled computational model. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2020; 36:e3398. [PMID: 32857480 DOI: 10.1002/cnm.3398] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/10/2020] [Revised: 08/22/2020] [Accepted: 08/24/2020] [Indexed: 06/11/2023]
Abstract
Myocardial infarction (MI) is the most common cause of a heart failure, which occurs due to myocardial ischemia leading to left ventricular (LV) remodeling. LV remodeling particularly occurs at the ischemic area and the region surrounds it, known as the border zone. The role of the border zone in initiating LV remodeling process urges the investigation on the correlation between early border zone changes and remodeling outcome. Thus, this study aims to simulate a preliminary conceptual work of the border zone formation and evolution during onset of MI and its effect towards early LV remodeling processes by incorporating the oxygen concentration effect on the electrophysiology of an idealized three-dimensional LV through electro-chemical coupled mathematical model. The simulation result shows that the region of border zone, represented by the distribution of electrical conductivities, keeps expanding over time. Based on this result, the border zone is also proposed to consist of three sub-regions, namely mildly, moderately, and seriously impaired conductivity regions, which each region categorized depending on its electrical conductivities. This division could be used as a biomarker for classification of reversible and irreversible myocardial injury and will help to identify the different risks for the survival of patient. Larger ischemic size and complete occlusion of the coronary artery can be associated with an increased risk of developing irreversible injury, in particular if the reperfusion treatment is delayed. Increased irreversible injury area can be related with cardiovascular events and will further deteriorate the LV function over time.
Collapse
Affiliation(s)
- Wan N Wan Ab Naim
- Faculty of Mechanical and Automotive Engineering Technology, University Malaysia Pahang, Pekan, Malaysia
| | - Mohd J Mohamed Mokhtarudin
- Department of Mechanical Engineering, College of Engineering, University Malaysia Pahang, Kuantan, Malaysia
| | - Einly Lim
- Department of Biomedical Engineering, Faculty of Engineering, University of Malaya, Kuala Lumpur, Malaysia
| | - Bee T Chan
- Department of Mechanical, Materials and Manufacturing Engineering, Faculty of Science and Engineering, University of Nottingham, Semenyih, Malaysia
| | - Azam Ahmad Bakir
- University of Southampton Malaysia Campus, Iskandar Puteri, Malaysia
| | - Nik A Nik Mohamed
- Faculty of Mechanical and Automotive Engineering Technology, University Malaysia Pahang, Pekan, Malaysia
| |
Collapse
|
29
|
Nguyen TD, Kadri OE, Voronov RS. An Introductory Overview of Image-Based Computational Modeling in Personalized Cardiovascular Medicine. Front Bioeng Biotechnol 2020; 8:529365. [PMID: 33102452 PMCID: PMC7546862 DOI: 10.3389/fbioe.2020.529365] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2020] [Accepted: 08/31/2020] [Indexed: 02/05/2023] Open
Abstract
Cardiovascular diseases account for the number one cause of deaths in the world. Part of the reason for such grim statistics is our limited understanding of the underlying mechanisms causing these devastating pathologies, which is made difficult by the invasiveness of the procedures associated with their diagnosis (e.g., inserting catheters into the coronal artery to measure blood flow to the heart). Likewise, it is also difficult to design and test assistive devices without implanting them in vivo. However, with the recent advancements made in biomedical scanning technologies and computer simulations, image-based modeling (IBM) has arisen as the next logical step in the evolution of non-invasive patient-specific cardiovascular medicine. Yet, due to its novelty, it is still relatively unknown outside of the niche field. Therefore, the goal of this manuscript is to review the current state-of-the-art and the limitations of the methods used in this area of research, as well as their applications to personalized cardiovascular investigations and treatments. Specifically, the modeling of three different physics – electrophysiology, biomechanics and hemodynamics – used in the cardiovascular IBM is discussed in the context of the physiology that each one of them describes and the mechanisms of the underlying cardiac diseases that they can provide insight into. Only the “bare-bones” of the modeling approaches are discussed in order to make this introductory material more accessible to an outside observer. Additionally, the imaging methods, the aspects of the unique cardiac anatomy derived from them, and their relation to the modeling algorithms are reviewed. Finally, conclusions are drawn about the future evolution of these methods and their potential toward revolutionizing the non-invasive diagnosis, virtual design of treatments/assistive devices, and increasing our understanding of these lethal cardiovascular diseases.
Collapse
Affiliation(s)
- Thanh Danh Nguyen
- Otto H. York Department of Chemical and Materials Engineering, Newark College of Engineering, New Jersey Institute of Technology, Newark, NJ, United States
| | - Olufemi E Kadri
- Otto H. York Department of Chemical and Materials Engineering, Newark College of Engineering, New Jersey Institute of Technology, Newark, NJ, United States.,UC-P&G Simulation Center, University of Cincinnati, Cincinnati, OH, United States
| | - Roman S Voronov
- Otto H. York Department of Chemical and Materials Engineering, Newark College of Engineering, New Jersey Institute of Technology, Newark, NJ, United States.,Department of Biomedical Engineering, Newark College of Engineering, New Jersey Institute of Technology, Newark, NJ, United States
| |
Collapse
|
30
|
Wan Ab Naim WN, Mohamed Mokhtarudin MJ, Chan BT, Lim E, Ahmad Bakir A, Nik Mohamed NA. The study of myocardial ischemia-reperfusion treatment through computational modelling. J Theor Biol 2020; 509:110527. [PMID: 33096094 DOI: 10.1016/j.jtbi.2020.110527] [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: 03/31/2020] [Revised: 09/25/2020] [Accepted: 10/15/2020] [Indexed: 10/23/2022]
Abstract
Reperfusion of the blood flow to ischemic myocardium is the standard treatment for patients suffering myocardial infarction. However, the reperfusion itself can also induce myocardial injury, in which the actual mechanism and its risk factors remain unclear. This work aims to study the mechanism of ischemia-reperfusion treatment using a three-dimensional (3D) oxygen diffusion model. An electrical model is then coupled to an oxygen model to identify the possible region of myocardial damage. Our findings show that the value of oxygen exceeds its optimum (>1.0) at the ischemic area during early reperfusion period. This complication was exacerbated in a longer ischemic period. While a longer reperfusion time causes a continuous excessive oxygen supply to the ischemic area throughout the reperfusion time. This work also suggests the use of less than 0.8 of initial oxygen concentration in the reperfusion treatment to prevent undesired upsurge at the early reperfusion period and further myocardial injury. We also found the region at risk for myocardial injury is confined in the ischemic vicinity revealed by its electrical conductivity impairment. Although there is a risk that reperfusion leads to myocardial injury for excessive oxygen accumulation, the reperfusion treatment is helpful in reducing the infarct size.
Collapse
Affiliation(s)
- Wan Naimah Wan Ab Naim
- Faculty of Mechanical and Automotive Engineering Technology, University Malaysia Pahang, 26600 Pekan, Pahang, Malaysia
| | - Mohd Jamil Mohamed Mokhtarudin
- Department of Mechanical Engineering, College of Engineering, University Malaysia Pahang, Lebuhraya Tun Razak, 26300 Gambang, Kuantan, Pahang, Malaysia.
| | - Bee Ting Chan
- Department of Mechanical, Materials and Manufacturing Engineering, Faculty of Science and Engineering, University of Nottingham, 43500 Selangor, Malaysia
| | - Einly Lim
- Department of Biomedical Engineering, Faculty of Engineering, University of Malaya, 50603 Kuala Lumpur, Malaysia
| | - Azam Ahmad Bakir
- University of Southampton Malaysia Campus, No 3, Persiaran Canselor 1, Kota Ilmu Educity, 79200 Iskandar Puteri, Johor, Malaysia
| | - Nik Abdullah Nik Mohamed
- Faculty of Mechanical and Automotive Engineering Technology, University Malaysia Pahang, 26600 Pekan, Pahang, Malaysia
| |
Collapse
|
31
|
Corrado C, Avezzù A, Lee AWC, Mendoca Costa C, Roney CH, Strocchi M, Bishop M, Niederer SA. Using cardiac ionic cell models to interpret clinical data. WIREs Mech Dis 2020; 13:e1508. [PMID: 33027553 DOI: 10.1002/wsbm.1508] [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: 05/31/2020] [Revised: 08/27/2020] [Accepted: 09/04/2020] [Indexed: 01/24/2023]
Abstract
For over 100 years cardiac electrophysiology has been measured in the clinic. The electrical signals that can be measured span from noninvasive ECG and body surface potentials measurements through to detailed invasive measurements of local tissue electrophysiology. These electrophysiological measurements form a crucial component of patient diagnosis and monitoring; however, it remains challenging to quantitatively link changes in clinical electrophysiology measurements to biophysical cellular function. Multi-scale biophysical computational models represent one solution to this problem. These models provide a formal framework for linking cellular function through to emergent whole organ function and routine clinical diagnostic signals. In this review, we describe recent work on the use of computational models to interpret clinical electrophysiology signals. We review the simulation of human cardiac myocyte electrophysiology in the atria and the ventricles and how these models are being used to link organ scale function to patient disease mechanisms and therapy response in patients receiving implanted defibrillators, \cardiac resynchronisation therapy or suffering from atrial fibrillation and ventricular tachycardia. There is a growing use of multi-scale biophysical models to interpret clinical data. This allows cardiologists to link clinical observations with cellular mechanisms to better understand cardiopathophysiology and identify novel treatment strategies. This article is categorized under: Cardiovascular Diseases > Computational Models Cardiovascular Diseases > Biomedical Engineering Cardiovascular Diseases > Molecular and Cellular Physiology.
Collapse
|
32
|
Sutanto H, Cluitmans MJM, Dobrev D, Volders PGA, Bébarová M, Heijman J. Acute effects of alcohol on cardiac electrophysiology and arrhythmogenesis: Insights from multiscale in silico analyses. J Mol Cell Cardiol 2020; 146:69-83. [PMID: 32710981 DOI: 10.1016/j.yjmcc.2020.07.007] [Citation(s) in RCA: 27] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 05/01/2020] [Revised: 06/27/2020] [Accepted: 07/16/2020] [Indexed: 12/11/2022]
Abstract
Acute excessive ethyl alcohol (ethanol) consumption alters cardiac electrophysiology and can evoke cardiac arrhythmias, e.g., in 'holiday heart syndrome'. Ethanol acutely modulates numerous targets in cardiomyocytes, including ion channels, Ca2+-handling proteins and gap junctions. However, the mechanisms underlying ethanol-induced arrhythmogenesis remain incompletely understood and difficult to study experimentally due to the multiple electrophysiological targets involved and their potential interactions with preexisting electrophysiological or structural substrates. Here, we employed cellular- and tissue-level in-silico analyses to characterize the acute effects of ethanol on cardiac electrophysiology and arrhythmogenesis. Acute electrophysiological effects of ethanol were incorporated into human atrial and ventricular cardiomyocyte computer models: reduced INa, ICa,L, Ito, IKr and IKur, dual effects on IK1 and IK,ACh (inhibition at low and augmentation at high concentrations), and increased INCX and SR Ca2+ leak. Multiscale simulations in the absence or presence of preexistent atrial fibrillation or heart-failure-related remodeling demonstrated that low ethanol concentrations prolonged atrial action-potential duration (APD) without effects on ventricular APD. Conversely, high ethanol concentrations abbreviated atrial APD and prolonged ventricular APD. High ethanol concentrations promoted reentry in tissue simulations, but the extent of reentry promotion was dependent on the presence of altered intercellular coupling, and the degree, type, and pattern of fibrosis. Taken together, these data provide novel mechanistic insight into the potential proarrhythmic interactions between a preexisting substrate and acute changes in cardiac electrophysiology. In particular, acute ethanol exposure has concentration-dependent electrophysiological effects that differ between atria and ventricles, and between healthy and diseased hearts. Low concentrations of ethanol can have anti-fibrillatory effects in atria, whereas high concentrations promote the inducibility and maintenance of reentrant atrial and ventricular arrhythmias, supporting a role for limiting alcohol intake as part of cardiac arrhythmia management.
Collapse
Affiliation(s)
- Henry Sutanto
- Department of Cardiology, CARIM School for Cardiovascular Diseases, Maastricht University, the Netherlands
| | - Matthijs J M Cluitmans
- Department of Cardiology, CARIM School for Cardiovascular Diseases, Maastricht University, the Netherlands
| | - Dobromir Dobrev
- Institute of Pharmacology, West German Heart and Vascular Center, University Duisburg-Essen, Essen, Germany
| | - Paul G A Volders
- Department of Cardiology, CARIM School for Cardiovascular Diseases, Maastricht University, the Netherlands
| | - Markéta Bébarová
- Department of Physiology, Faculty of Medicine, Masaryk University, Brno, Czech Republic
| | - Jordi Heijman
- Department of Cardiology, CARIM School for Cardiovascular Diseases, Maastricht University, the Netherlands.
| |
Collapse
|
33
|
Human Purkinje in silico model enables mechanistic investigations into automaticity and pro-arrhythmic abnormalities. J Mol Cell Cardiol 2020; 142:24-38. [PMID: 32251669 PMCID: PMC7294239 DOI: 10.1016/j.yjmcc.2020.04.001] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 10/24/2019] [Revised: 03/30/2020] [Accepted: 04/01/2020] [Indexed: 02/06/2023]
Abstract
Cardiac Purkinje cells (PCs) are implicated in lethal arrhythmias caused by cardiac diseases, mutations, and drug action. However, the pro-arrhythmic mechanisms in PCs are not entirely understood, particularly in humans, as most investigations are conducted in animals. The aims of this study are to present a novel human PCs electrophysiology biophysically-detailed computational model, and to disentangle ionic mechanisms of human Purkinje-related electrophysiology, pacemaker activity and arrhythmogenicity. The new Trovato2020 model incorporates detailed Purkinje-specific ionic currents and Ca2+ handling, and was developed, calibrated and validated using human experimental data acquired at multiple frequencies, both in control conditions and following drug application. Multiscale investigations were performed in a Purkinje cell, in fibre and using an experimentally-calibrated population of PCs to evaluate biological variability. Simulations demonstrate the human Purkinje Trovato2020 model is the first one to yield: (i) all key AP features consistent with human Purkinje recordings; (ii) Automaticity with funny current up-regulation (iii) EADs at slow pacing and with 85% hERG block; (iv) DADs following fast pacing; (v) conduction velocity of 160 cm/s in a Purkinje fibre, as reported in human. The human in silico PCs population highlights that: (1) EADs are caused by ICaL reactivation in PCs with large inward currents; (2) DADs and triggered APs occur in PCs experiencing Ca2+ accumulation, at fast pacing, caused by large L-type calcium current and small Na+/Ca2+ exchanger. The novel human Purkinje model unlocks further investigations into the role of cardiac Purkinje in ventricular arrhythmias through computer modeling and multiscale simulations. A human in silico AP model was developed to investigate arrhythmia in cardiac Purkinje. The new Purkinje model enables multiscale investigations from single cell to tissue. Populations of human Purkinje models reproduce and explain experimental variability. Ca2+-current reactivation triggers EADs in virtual Purkinje cells with weak repolarisation reserve. Ca2+ accumulation caused by increased Ca2+ and NCX currents triggers DADs.
Collapse
|
34
|
Morissette P, Polak S, Chain A, Zhai J, Imredy JP, Wildey MJ, Travis J, Fitzgerald K, Fanelli P, Passini E, Rodriguez B, Sannajust F, Regan C. Combining an in silico proarrhythmic risk assay with a tPKPD model to predict QTc interval prolongation in the anesthetized guinea pig assay. Toxicol Appl Pharmacol 2020; 390:114883. [PMID: 31981640 PMCID: PMC7322544 DOI: 10.1016/j.taap.2020.114883] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2019] [Revised: 01/03/2020] [Accepted: 01/14/2020] [Indexed: 12/12/2022]
Abstract
Human-based in silico models are emerging as important tools to study the effects of integrating inward and outward ion channel currents to predict clinical proarrhythmic risk. The aims of this study were 2-fold: 1) Evaluate the capacity of an in silico model to predict QTc interval prolongation in the in vivo anesthetized cardiovascular guinea pig (CVGP) assay for new chemical entities (NCEs) and; 2) Determine if a translational pharmacokinetic/pharmacodynamic (tPKPD) model can improve the predictive capacity. In silico simulations for NCEs were performed using a population of human ventricular action potential (AP) models. PatchXpress® (PX) or high throughput screening (HTS) ion channel data from respectively n = 73 and n = 51 NCEs were used as inputs for the in silico population. These NCEs were also tested in the CVGP (n = 73). An M5 pruned decision tree-based regression tPKPD model was used to evaluate the concentration at which an NCE is liable to prolong the QTc interval in the CVGP. In silico results successfully predicted the QTc interval prolongation outcome observed in the CVGP with an accuracy/specificity of 85%/73% and 75%/77%, when using PX and HTS ion channel data, respectively. Considering the tPKPD predicted concentration resulting in QTc prolongation (EC5%) increased accuracy/specificity to 97%/95% using PX and 88%/97% when using HTS. Our results support that human-based in silico simulations in combination with tPKPD modeling can provide correlative results with a commonly used early in vivo safety assay, suggesting a path toward more rapid NCE assessment with reduced resources, cycle time, and animal use. Cardiac electrophysiological in silico model predicts QTc interval prolongation in the guinea pig. PKPD model predicts relevant QTc interval prolongation concentration in guinea pig. Combining the models improves the accuracy of predicting guinea pig QTc effects. Combining models accelerates assessment of QTc with lower resources and animal use.
Collapse
Affiliation(s)
- Pierre Morissette
- Safety Assessment & Laboratory Animal Resources (SALAR), Merck & Co., Inc., West Point, PA, USA.
| | - Sebastian Polak
- Certara UK Limited, Simcyp Division, Sheffield, UK; Jagiellonian University Medical College, Faculty of Pharmacy, Krakow, Poland
| | - Anne Chain
- Pharmacokinetics, Pharmacodynamics and Drug Metabolism (PPDM), Merck & Co., Inc., Rahway, NJ, USA
| | - Jin Zhai
- Safety Assessment & Laboratory Animal Resources (SALAR), Merck & Co., Inc., West Point, PA, USA
| | - John P Imredy
- Safety Assessment & Laboratory Animal Resources (SALAR), Merck & Co., Inc., West Point, PA, USA
| | - Mary Jo Wildey
- Pharmacology, Screening and Informatics, Merck & Co., Kenilworth, NJ, USA
| | - Jeffrey Travis
- Safety Assessment & Laboratory Animal Resources (SALAR), Merck & Co., Inc., West Point, PA, USA
| | - Kevin Fitzgerald
- Safety Assessment & Laboratory Animal Resources (SALAR), Merck & Co., Inc., West Point, PA, USA
| | - Patrick Fanelli
- Safety Assessment & Laboratory Animal Resources (SALAR), Merck & Co., Inc., West Point, PA, USA
| | - Elisa Passini
- Computational Cardiovascular Science Group, Department of Computer Science, BHF Centre of Research Excellence, University of Oxford, Oxford, UK
| | - Blanca Rodriguez
- Computational Cardiovascular Science Group, Department of Computer Science, BHF Centre of Research Excellence, University of Oxford, Oxford, UK
| | - Frederick Sannajust
- Safety Assessment & Laboratory Animal Resources (SALAR), Merck & Co., Inc., West Point, PA, USA
| | - Christopher Regan
- Safety Assessment & Laboratory Animal Resources (SALAR), Merck & Co., Inc., West Point, PA, USA
| |
Collapse
|
35
|
Trayanova NA, Doshi AN, Prakosa A. How personalized heart modeling can help treatment of lethal arrhythmias: A focus on ventricular tachycardia ablation strategies in post-infarction patients. WILEY INTERDISCIPLINARY REVIEWS-SYSTEMS BIOLOGY AND MEDICINE 2020; 12:e1477. [PMID: 31917524 DOI: 10.1002/wsbm.1477] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/12/2019] [Revised: 12/16/2019] [Accepted: 12/17/2019] [Indexed: 12/18/2022]
Abstract
Precision Cardiology is a targeted strategy for cardiovascular disease prevention and treatment that accounts for individual variability. Computational heart modeling is one of the novel approaches that have been developed under the umbrella of Precision Cardiology. Personalized computational modeling of patient hearts has made strides in the development of models that incorporate the individual geometry and structure of the heart as well as other patient-specific information. Of these developments, one of the potentially most impactful is the research aimed at noninvasively predicting the targets of ablation of lethal arrhythmia, ventricular tachycardia (VT), using patient-specific models. The approach has been successfully applied to patients with ischemic cardiomyopathy in proof-of-concept studies. The goal of this paper is to review the strategies for computational VT ablation guidance in ischemic cardiomyopathy patients, from model developments to the intricacies of the actual clinical application. To provide context in describing the road these computational modeling applications have undertaken, we first review the state of the art in VT ablation in the clinic, emphasizing the benefits that personalized computational prediction of ablation targets could bring to the clinical electrophysiology practice. This article is characterized under: Analytical and Computational Methods > Computational Methods Models of Systems Properties and Processes > Organ, Tissue, and Physiological Models Translational, Genomic, and Systems Medicine > Translational Medicine.
Collapse
Affiliation(s)
- Natalia A Trayanova
- Alliance for Cardiovascular Diagnostic and Treatment Innovation, Johns Hopkins University, Baltimore, Maryland.,Department of Biomedical Engineering, Johns Hopkins University, Baltimore, Maryland
| | - Ashish N Doshi
- Alliance for Cardiovascular Diagnostic and Treatment Innovation, Johns Hopkins University, Baltimore, Maryland
| | - Adityo Prakosa
- Alliance for Cardiovascular Diagnostic and Treatment Innovation, Johns Hopkins University, Baltimore, Maryland
| |
Collapse
|