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] [What about the content of this article? (0)] [Affiliation(s)] [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
|
Cumberland MJ, Riebel LL, Roy A, O’Shea C, Holmes AP, Denning C, Kirchhof P, Rodriguez B, Gehmlich K. Basic Research Approaches to Evaluate Cardiac Arrhythmia in Heart Failure and Beyond. Front Physiol 2022; 13:806366. [PMID: 35197863 PMCID: PMC8859441 DOI: 10.3389/fphys.2022.806366] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2021] [Accepted: 01/10/2022] [Indexed: 12/20/2022] Open
Abstract
Patients with heart failure often develop cardiac arrhythmias. The mechanisms and interrelations linking heart failure and arrhythmias are not fully understood. Historically, research into arrhythmias has been performed on affected individuals or in vivo (animal) models. The latter however is constrained by interspecies variation, demands to reduce animal experiments and cost. Recent developments in in vitro induced pluripotent stem cell technology and in silico modelling have expanded the number of models available for the evaluation of heart failure and arrhythmia. An agnostic approach, combining the modalities discussed here, has the potential to improve our understanding for appraising the pathology and interactions between heart failure and arrhythmia and can provide robust and validated outcomes in a variety of research settings. This review discusses the state of the art models, methodologies and techniques used in the evaluation of heart failure and arrhythmia and will highlight the benefits of using them in combination. Special consideration is paid to assessing the pivotal role calcium handling has in the development of heart failure and arrhythmia.
Collapse
Affiliation(s)
- Max J. Cumberland
- Institute of Cardiovascular Sciences, College of Medical and Dental Sciences, University of Birmingham, Birmingham, United Kingdom
| | - Leto L. Riebel
- Department of Computer Science, University of Oxford, Oxford, United Kingdom
| | - Ashwin Roy
- Institute of Cardiovascular Sciences, College of Medical and Dental Sciences, University of Birmingham, Birmingham, United Kingdom
| | - Christopher O’Shea
- Institute of Cardiovascular Sciences, College of Medical and Dental Sciences, University of Birmingham, Birmingham, United Kingdom
| | - Andrew P. Holmes
- Institute of Cardiovascular Sciences, College of Medical and Dental Sciences, University of Birmingham, Birmingham, United Kingdom
- Institute of Clinical Sciences, College of Medical and Dental Sciences, University of Birmingham, Birmingham, United Kingdom
| | - Chris Denning
- Stem Cell Biology Unit, Biodiscovery Institute, British Heart Foundation Centre for Regenerative Medicine, University of Nottingham, Nottingham, United Kingdom
| | - Paulus Kirchhof
- Institute of Cardiovascular Sciences, College of Medical and Dental Sciences, University of Birmingham, Birmingham, United Kingdom
- University Heart and Vascular Center, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
| | - Blanca Rodriguez
- Department of Computer Science, University of Oxford, Oxford, United Kingdom
| | - Katja Gehmlich
- Institute of Cardiovascular Sciences, College of Medical and Dental Sciences, University of Birmingham, Birmingham, United Kingdom
- Cardiovascular Medicine, Radcliffe Department of Medicine, University of Oxford and British Heart Foundation Centre of Research Excellence Oxford, Oxford, United Kingdom
| |
Collapse
|
3
|
Sung E, Etoz S, Zhang Y, Trayanova NA. Whole-heart ventricular arrhythmia modeling moving forward: Mechanistic insights and translational applications. Biophys Rev (Melville) 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] [What about the content of this article? (0)] [Affiliation(s)] [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
|
4
|
Pagani S, Manzoni A. Enabling forward uncertainty quantification and sensitivity analysis in cardiac electrophysiology by reduced order modeling and machine learning. Int J Numer Method Biomed Eng 2021; 37:e3450. [PMID: 33599106 PMCID: PMC8244126 DOI: 10.1002/cnm.3450] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Revised: 02/05/2021] [Accepted: 02/07/2021] [Indexed: 06/12/2023]
Abstract
We present a new, computationally efficient framework to perform forward uncertainty quantification (UQ) in cardiac electrophysiology. We consider the monodomain model to describe the electrical activity in the cardiac tissue, coupled with the Aliev-Panfilov model to characterize the ionic activity through the cell membrane. We address a complete forward UQ pipeline, including both: (i) a variance-based global sensitivity analysis for the selection of the most relevant input parameters, and (ii) a way to perform uncertainty propagation to investigate the impact of intra-subject variability on outputs of interest depending on the cardiac potential. Both tasks exploit stochastic sampling techniques, thus implying overwhelming computational costs because of the huge amount of queries to the high-fidelity, full-order computational model obtained by approximating the coupled monodomain/Aliev-Panfilov system through the finite element method. To mitigate this computational burden, we replace the full-order model with computationally inexpensive projection-based reduced-order models (ROMs) aimed at reducing the state-space dimensionality. Resulting approximation errors on the outputs of interest are finally taken into account through artificial neural network (ANN)-based models, enhancing the accuracy of the whole UQ pipeline. Numerical results show that the proposed physics-based ROMs outperform regression-based emulators relying on ANNs built with the same amount of training data, in terms of both numerical accuracy and overall computational efficiency.
Collapse
Affiliation(s)
- Stefano Pagani
- MOX, Dipartimento di MatematicaPolitecnico di MilanoMilanItaly
| | - Andrea Manzoni
- MOX, Dipartimento di MatematicaPolitecnico di MilanoMilanItaly
| |
Collapse
|
5
|
Ushenin K, Kalinin V, Gitinova S, Sopov O, Solovyova O. Parameter variations in personalized electrophysiological models of human heart ventricles. PLoS One 2021; 16:e0249062. [PMID: 33909606 PMCID: PMC8081243 DOI: 10.1371/journal.pone.0249062] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/18/2020] [Accepted: 03/10/2021] [Indexed: 11/18/2022] Open
Abstract
The objectives of this study were to evaluate the accuracy of personalized numerical simulations of the electrical activity in human ventricles by comparing simulated electrocardiograms (ECGs) with real patients’ ECGs and analyzing the sensitivity of the model output to variations in the model parameters. We used standard 12-lead ECGs and up to 224 unipolar body-surface ECGs to record three patients with cardiac resynchronization therapy devices and three patients with focal ventricular tachycardia. Patient-tailored geometrical models of the ventricles, atria, large vessels, liver, and spine were created using computed tomography data. Ten cases of focal ventricular activation were simulated using the bidomain model and the TNNP 2006 cellular model. The population-based values of electrical conductivities and other model parameters were used for accuracy analysis, and their variations were used for sensitivity analysis. The mean correlation coefficient between the simulated and real ECGs varied significantly (from r = 0.29 to r = 0.86) among the simulated cases. A strong mean correlation (r > 0.7) was found in eight of the ten model cases. The accuracy of the ECG simulation varied widely in the same patient depending on the localization of the excitation origin. The sensitivity analysis revealed that variations in the anisotropy ratio, blood conductivity, and cellular apicobasal heterogeneity had the strongest influence on transmembrane potential, while variation in lung conductivity had the greatest influence on body-surface ECGs. Futhermore, the anisotropy ratio predominantly affected the latest activation time and repolarization time dispersion, while the cellular apicobasal heterogeneity mainly affected the dispersion of action potential duration, and variation in lung conductivity mainly led to changes in the amplitudes of ECGs and cardiac electrograms. We also found that the effects of certain parameter variations had specific regional patterns on the cardiac and body surfaces. These observations are useful for further developing personalized cardiac models.
Collapse
Affiliation(s)
- Konstantin Ushenin
- Institute of Natural Sciences and Mathematics, Ural Federal University, Ekaterinburg, Russia
- Institute of Immunology and Physiology of the Ural Branch of the RAS, Ekaterinburg, Russia
- * E-mail:
| | | | - Sukaynat Gitinova
- Department of Surgical Treatment of Tachyarrhythmias, A.N. Bakulev National Medical Research Center of Cardiovascular Surgery, Moscow, Russia
| | - Oleg Sopov
- Department of Surgical Treatment of Tachyarrhythmias, A.N. Bakulev National Medical Research Center of Cardiovascular Surgery, Moscow, Russia
| | - Olga Solovyova
- Institute of Natural Sciences and Mathematics, Ural Federal University, Ekaterinburg, Russia
- Institute of Immunology and Physiology of the Ural Branch of the RAS, Ekaterinburg, Russia
| |
Collapse
|
6
|
ZHU HONGLEI, JIN LIAN, ZHANG JIAYU, WU XIAOMEI. OPTIMIZATION OF RABBIT VENTRICULAR ELECTROPHYSIOLOGICAL MODEL AND SIMULATION OF SYNTHETIC ELECTROCARDIOGRAM. J MECH MED BIOL 2021. [DOI: 10.1142/s0219519421500019] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
This study aimed to use computer simulation method to study the mechanism of cardiac electrical activities. We optimized an electrophysiological rabbit ventricular model, including myocardial segmentation, heterogeneity and a realistic His-Purkinje network. Simulations of normal state, several types of ventricular premature contractions (VPC), conduction system pacing and right ventricular apical pacing were performed and the detailed cardiac electrical activities were studied from cell level to electrocardiogram (ECG) level. A detailed multiscale optimized ventricular model was obtained. The model effectively simulated various types of electrical activities. The synthetic ECG results were very similar to the real clinical ECG. The duration of QRS of typical VPC is 58[Formula: see text]ms, 71% longer than that of a normal-state synthetic QRS and the amplitude of the QRS is 35% larger, while the QRS duration and amplitude of the real clinical ECG of typical VPC are 69% longer and 36% larger than those of the real normal QRS. The duration of QRS of ventricular fusion beat is 31[Formula: see text]ms, 91% of that of a normal-state synthetic QRS and the amplitude of the QRS is 36% larger, while the QRS duration of the real clinical ECG of a ventricular fusion beat is 92% of the real normal QRS and the amplitude is 37% larger. Therefore, the results indicate that this model is effective and reliable in studying the detailed process of cardiac excitation and pacing.
Collapse
Affiliation(s)
- HONGLEI ZHU
- Department of Electronic Engineering, Fudan University, Shanghai 200433, P. R. China
| | - LIAN JIN
- Department of Electronic Engineering, Fudan University, Shanghai 200433, P. R. China
| | - JIAYU ZHANG
- Department of Electronic Engineering, Fudan University, Shanghai 200433, P. R. China
| | - XIAOMEI WU
- Department of Electronic Engineering, Fudan University, Key Laboratory of Medical Imaging Computing and Computer Assisted Intervention (MICCAI) of Shanghai, Research Center of Assistive Devices, Shanghai, P. R. China
| |
Collapse
|
7
|
Subramaniam DR, Unnikrishnan G, Sundaramurthy A, Rubio JE, Kote VB, Reifman J. The importance of modeling the human cerebral vasculature in blunt trauma. Biomed Eng Online 2021; 20:11. [PMID: 33446217 PMCID: PMC7809851 DOI: 10.1186/s12938-021-00847-x] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2020] [Accepted: 01/04/2021] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Multiple studies describing human head finite element (FE) models have established the importance of including the major cerebral vasculature to improve the accuracy of the model predictions. However, a more detailed network of cerebral vasculature, including the major veins and arteries as well as their branch vessels, can further enhance the model-predicted biomechanical responses and help identify correlates to observed blunt-induced brain injury. METHODS We used an anatomically accurate three-dimensional geometry of a 50th percentile U.S. male head that included the skin, eyes, sinuses, spine, skull, brain, meninges, and a detailed network of cerebral vasculature to develop a high-fidelity model. We performed blunt trauma simulations and determined the intracranial pressure (ICP), the relative displacement (RD), the von Mises stress, and the maximum principal strain. We validated our detailed-vasculature model by comparing the model-predicted ICP and RD values with experimental measurements. To quantify the influence of including a more comprehensive network of brain vessels, we compared the biomechanical responses of our detailed-vasculature model with those of a reduced-vasculature model and a no-vasculature model. RESULTS For an inclined frontal impact, the predicted ICP matched well with the experimental results in the fossa, frontal, parietal, and occipital lobes, with peak-pressure differences ranging from 2.4% to 9.4%. For a normal frontal impact, the predicted ICP matched the experimental results in the frontal lobe and lateral ventricle, with peak-pressure discrepancies equivalent to 1.9% and 22.3%, respectively. For an offset parietal impact, the model-predicted RD matched well with the experimental measurements, with peak RD differences of 27% and 24% in the right and left cerebral hemispheres, respectively. Incorporating the detailed cerebral vasculature did not influence the ICP but redistributed the brain-tissue stresses and strains by as much as 30%. In addition, our detailed-vasculature model predicted strain reductions by as much as 28% when compared to current reduced-vasculature FE models that only include the major cerebral vessels. CONCLUSIONS Our study highlights the importance of including a detailed representation of the cerebral vasculature in FE models to more accurately estimate the biomechanical responses of the human brain to blunt impact.
Collapse
Affiliation(s)
- Dhananjay Radhakrishnan Subramaniam
- Department of Defense Biotechnology High Performance Computing Software Applications Institute, Telemedicine and Advanced Technology Research Center, United States Army Medical Research and Development Command, FCMR-TT, 504 Scott Street, Fort Detrick, MD, 21702-5012, USA.,The Henry M. Jackson Foundation for the Advancement of Military Medicine, Inc, 6720A Rockledge Drive, Bethesda, MD, 20817, USA
| | - Ginu Unnikrishnan
- Department of Defense Biotechnology High Performance Computing Software Applications Institute, Telemedicine and Advanced Technology Research Center, United States Army Medical Research and Development Command, FCMR-TT, 504 Scott Street, Fort Detrick, MD, 21702-5012, USA.,The Henry M. Jackson Foundation for the Advancement of Military Medicine, Inc, 6720A Rockledge Drive, Bethesda, MD, 20817, USA
| | - Aravind Sundaramurthy
- Department of Defense Biotechnology High Performance Computing Software Applications Institute, Telemedicine and Advanced Technology Research Center, United States Army Medical Research and Development Command, FCMR-TT, 504 Scott Street, Fort Detrick, MD, 21702-5012, USA.,The Henry M. Jackson Foundation for the Advancement of Military Medicine, Inc, 6720A Rockledge Drive, Bethesda, MD, 20817, USA
| | - Jose E Rubio
- Department of Defense Biotechnology High Performance Computing Software Applications Institute, Telemedicine and Advanced Technology Research Center, United States Army Medical Research and Development Command, FCMR-TT, 504 Scott Street, Fort Detrick, MD, 21702-5012, USA.,The Henry M. Jackson Foundation for the Advancement of Military Medicine, Inc, 6720A Rockledge Drive, Bethesda, MD, 20817, USA
| | - Vivek Bhaskar Kote
- Department of Defense Biotechnology High Performance Computing Software Applications Institute, Telemedicine and Advanced Technology Research Center, United States Army Medical Research and Development Command, FCMR-TT, 504 Scott Street, Fort Detrick, MD, 21702-5012, USA.,The Henry M. Jackson Foundation for the Advancement of Military Medicine, Inc, 6720A Rockledge Drive, Bethesda, MD, 20817, USA
| | - Jaques Reifman
- Department of Defense Biotechnology High Performance Computing Software Applications Institute, Telemedicine and Advanced Technology Research Center, United States Army Medical Research and Development Command, FCMR-TT, 504 Scott Street, Fort Detrick, MD, 21702-5012, USA.
| |
Collapse
|
8
|
Johnston BM, Johnston PR. Approaches for determining cardiac bidomain conductivity values: progress and challenges. Med Biol Eng Comput 2020; 58:2919-35. [PMID: 33089458 DOI: 10.1007/s11517-020-02272-z] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2020] [Accepted: 09/17/2020] [Indexed: 10/23/2022]
Abstract
Modelling the electrical activity of the heart is an important tool for understanding electrical function in various diseases and conduction disorders. Clearly, for model results to be useful, it is necessary to have accurate inputs for the models, in particular the commonly used bidomain model. However, there are only three sets of four experimentally determined conductivity values for cardiac ventricular tissue and these are inconsistent, were measured around 40 years ago, often produce different results in simulations and do not fully represent the three-dimensional anisotropic nature of cardiac tissue. Despite efforts in the intervening years, difficulties associated with making the measurements and also determining the conductivities from the experimental data have not yet been overcome. In this review, we summarise what is known about the conductivity values, as well as progress to date in meeting the challenges associated with both the mathematical modelling and the experimental techniques. Graphical abstract Epicardial potential distributions, arising from a subendocardial ischaemic region, modelled using conductivity data from the indicated studies.
Collapse
|
9
|
Stella S, Vergara C, Maines M, Catanzariti D, Africa PC, Demattè C, Centonze M, Nobile F, Del Greco M, Quarteroni A. Integration of activation maps of epicardial veins in computational cardiac electrophysiology. Comput Biol Med 2020; 127:104047. [PMID: 33099220 DOI: 10.1016/j.compbiomed.2020.104047] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2020] [Revised: 10/06/2020] [Accepted: 10/06/2020] [Indexed: 12/16/2022]
Abstract
In this work we address the issue of validating the monodomain equation used in combination with the Bueno-Orovio ionic model for the prediction of the activation times in cardiac electro-physiology of the left ventricle. To this aim, we consider four patients who suffered from Left Bundle Branch Block (LBBB). We use activation maps performed at the septum as input data for the model and maps at the epicardial veins for the validation. In particular, a first set (half) of the latter are used to estimate the conductivities of the patient and a second set (the remaining half) to compute the errors of the numerical simulations. We find an excellent agreement between measures and numerical results. Our validated computational tool could be used to accurately predict activation times at the epicardial veins with a short mapping, i.e. by using only a part (the most proximal) of the standard acquisition points, thus reducing the invasive procedure and exposure to radiation.
Collapse
|
10
|
Clayton RH, Aboelkassem Y, Cantwell CD, Corrado C, Delhaas T, Huberts W, Lei CL, Ni H, Panfilov AV, Roney C, dos Santos RW. An audit of uncertainty in multi-scale cardiac electrophysiology models. Philos Trans A Math Phys Eng Sci 2020; 378:20190335. [PMID: 32448070 PMCID: PMC7287340 DOI: 10.1098/rsta.2019.0335] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Accepted: 03/16/2020] [Indexed: 05/21/2023]
Abstract
Models of electrical activation and recovery in cardiac cells and tissue have become valuable research tools, and are beginning to be used in safety-critical applications including guidance for clinical procedures and for drug safety assessment. As a consequence, there is an urgent need for a more detailed and quantitative understanding of the ways that uncertainty and variability influence model predictions. In this paper, we review the sources of uncertainty in these models at different spatial scales, discuss how uncertainties are communicated across scales, and begin to assess their relative importance. We conclude by highlighting important challenges that continue to face the cardiac modelling community, identifying open questions, and making recommendations for future studies. This article is part of the theme issue 'Uncertainty quantification in cardiac and cardiovascular modelling and simulation'.
Collapse
Affiliation(s)
- Richard H. Clayton
- Insigneo institute for in-silico Medicine and Department of Computer Science, University of Sheffield, Sheffield, UK
- e-mail:
| | - Yasser Aboelkassem
- Department of Bioengineering, University of California, San Diego, CA, USA
| | | | - Cesare Corrado
- Division of Imaging Sciences and Biomedical Engineering, King’s College London, London, UK
| | - Tammo Delhaas
- School of Cardiovascular Diseases, Maastricht University, Maastricht, The Netherlands
| | - Wouter Huberts
- School of Cardiovascular Diseases, Maastricht University, Maastricht, The Netherlands
- Department of Biomedical Engineering, Eindhoven University of Technology, Eindhoven, The Netherlands
| | - Chon Lok Lei
- Computational Biology and Health Informatics, Department of Computer Science, University of Oxford, Oxford, UK
| | - Haibo Ni
- Department of Pharmacology, University of California, Davis, CA, USA
| | - Alexander V. Panfilov
- Department of Physics and Astronomy, University of Gent, Gent, Belgium
- Laboratory of Computational Biology and Medicine, Ural Federal University, Ekaterinburg, Russia
| | - Caroline Roney
- Division of Imaging Sciences and Biomedical Engineering, King’s College London, London, UK
| | | |
Collapse
|
11
|
Park JR, Chung SP, Hwang SY, Shin TG, Park JE. Myocardial infarction evaluation from stopping time decision toward interoperable algorithmic states in reinforcement learning. BMC Med Inform Decis Mak 2020; 20:99. [PMID: 32487133 PMCID: PMC7472590 DOI: 10.1186/s12911-020-01133-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2019] [Accepted: 05/17/2020] [Indexed: 12/05/2022] Open
Abstract
BACKGROUND The Elliot wave principle commonly characterizes the impulsive and corrective wave trends for both financial market trends and electrocardiograms. The impulsive wave trends of electrocardiograms can annotate several wave components of heart-beats including pathological heartbeat waveforms. The stopping time inquires which ordinal element satisfies the assumed mathematical condition within a numerical set. The proposed work constitutes several algorithmic states in reinforcement learning from the stopping time decision, which determines the impulsive wave trends. Each proposed algorithmic state is applicable to any relevant algorithmic state in reinforcement learning with fully numerical explanations. Because commercial electrocardiographs still misinterpret myocardial infarctions from extraordinary electrocardiograms, a novel algorithm needs to be developed to evaluate myocardial infarctions. Moreover, differential diagnosis for right ventricle infarction is required to contraindicate a medication such as nitroglycerin. METHODS The proposed work implements the stopping time theory to impulsive wave trend distribution. The searching process of the stopping time theory is equivalent to the actions toward algorithmic states in reinforcement learning. The state value from each algorithmic state represents the numerically deterministic annotated results from the impulsive wave trend distribution. The shape of the impulsive waveform is evaluated from the interoperable algorithmic states via least-first-power approximation and approximate entropy. The annotated electrocardiograms from the impulsive wave trend distribution utilize a structure of neural networks to approximate the isoelectric baseline amplitude value of the electrocardiograms, and detect the conditions of myocardial infarction. The annotated results from the impulsive wave trend distribution consist of another reinforcement learning environment for the evaluation of impulsive waveform direction. RESULTS The accuracy to discern myocardial infarction was found to be 99.2754% for the data from the comma-separated value format files, and 99.3579% for those containing representative beats. The clinical dataset included 276 electrocardiograms from the comma-separated value files and 623 representative beats. CONCLUSIONS Our study aims to support clinical interpretation on 12-channel electrocardiograms. The proposed work is suitable for a differential diagnosis under infarction in the right ventricle to avoid contraindicated medication during emergency. An impulsive waveform that is affected by myocardial infarction or the electrical direction of electrocardiography is represented as an inverse waveform.
Collapse
Affiliation(s)
- Jong-Rul Park
- College of Information and Communication Engineering, Sungkyunkwan University, Suwon, 16419 Republic of Korea
| | - Sung Phil Chung
- Department of Emergency Medicine, Yonsei University Gangnam Severance Hospital, Seoul, 06273 Republic of Korea
| | - Sung Yeon Hwang
- Department of Emergency Medicine, Samsung Medical Center, Sungkyunkwan University School of Medicine, Seoul, 06351 Republic of Korea
| | - Tae Gun Shin
- Department of Emergency Medicine, Samsung Medical Center, Sungkyunkwan University School of Medicine, Seoul, 06351 Republic of Korea
| | - Jong Eun Park
- Department of Emergency Medicine, Samsung Medical Center, Sungkyunkwan University School of Medicine, Seoul, 06351 Republic of Korea
| |
Collapse
|
12
|
Abstract
The treatment of individual patients in cardiology practice increasingly relies on advanced imaging, genetic screening and devices. As the amount of imaging and other diagnostic data increases, paralleled by the greater capacity to personalize treatment, the difficulty of using the full array of measurements of a patient to determine an optimal treatment seems also to be paradoxically increasing. Computational models are progressively addressing this issue by providing a common framework for integrating multiple data sets from individual patients. These models, which are based on physiology and physics rather than on population statistics, enable computational simulations to reveal diagnostic information that would have otherwise remained concealed and to predict treatment outcomes for individual patients. The inherent need for patient-specific models in cardiology is clear and is driving the rapid development of tools and techniques for creating personalized methods to guide pharmaceutical therapy, deployment of devices and surgical interventions.
Collapse
|
13
|
Mincholé A, Zacur E, Ariga R, Grau V, Rodriguez B. MRI-Based Computational Torso/Biventricular Multiscale Models to Investigate the Impact of Anatomical Variability on the ECG QRS Complex. Front Physiol 2019; 10:1103. [PMID: 31507458 PMCID: PMC6718559 DOI: 10.3389/fphys.2019.01103] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2019] [Accepted: 08/08/2019] [Indexed: 01/07/2023] Open
Abstract
AIMS Patient-to-patient anatomical differences are an important source of variability in the electrocardiogram, and they may compromise the identification of pathological electrophysiological abnormalities. This study aims at quantifying the contribution of variability in ventricular and torso anatomies to differences in QRS complexes of the 12-lead ECG using computer simulations. METHODS A computational pipeline is presented that enables computer simulations using human torso/biventricular anatomically based electrophysiological models from clinically standard magnetic resonance imaging (MRI). The ventricular model includes membrane kinetics represented by the biophysically detailed O'Hara Rudy model modified for tissue heterogeneity and includes fiber orientation based on the Streeter rule. A population of 265 torso/biventricular models was generated by combining ventricular and torso anatomies obtained from clinically standard MRIs, augmented with a statistical shape model of the body. 12-lead ECGs were simulated on the 265 human torso/biventricular electrophysiology models, and QRS morphology, duration and amplitude were quantified in each ECG lead for each of the human torso-biventricular models. RESULTS QRS morphologies in limb leads are mainly determined by ventricular anatomy, while in the precordial leads, and especially V1 to V4, they are determined by heart position within the torso. Differences in ventricular orientation within the torso can explain morphological variability from monophasic to biphasic QRS complexes. QRS duration is mainly influenced by myocardial volume, while it is hardly affected by the torso anatomy or position. An average increase of 0.12 ± 0.05 ms in QRS duration is obtained for each cm3 of myocardial volume across all the leads while it hardly changed due to changes in torso volume. CONCLUSION Computer simulations using populations of human torso/biventricular models based on clinical MRI enable quantification of anatomical causes of variability in the QRS complex of the 12-lead ECG. The human models presented also pave the way toward their use as testbeds in silico clinical trials.
Collapse
Affiliation(s)
- Ana Mincholé
- Department of Computer Science, University of Oxford, Oxford, United Kingdom
| | - Ernesto Zacur
- Institute of Biomedical Engineering (IBME), University of Oxford, Oxford, United Kingdom
| | - Rina Ariga
- Division of Cardiovascular Medicine, Radcliffe Department of Medicine, University of Oxford, Oxford, United Kingdom
| | - Vicente Grau
- Institute of Biomedical Engineering (IBME), University of Oxford, Oxford, United Kingdom
| | - Blanca Rodriguez
- Department of Computer Science, University of Oxford, Oxford, United Kingdom
| |
Collapse
|
14
|
Costabal FS, Matsuno K, Yao J, Perdikaris P, Kuhl E. Machine learning in drug development: Characterizing the effect of 30 drugs on the QT interval using Gaussian process regression, sensitivity analysis, and uncertainty quantification. Comput Methods Appl Mech Eng 2019; 348:313-333. [PMID: 32863454 PMCID: PMC7454226 DOI: 10.1016/j.cma.2019.01.033] [Citation(s) in RCA: 46] [Impact Index Per Article: 9.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/05/2023]
Abstract
Prolonged QT intervals are a major risk factor for ventricular arrhythmias and a leading cause of sudden cardiac death. Various drugs are known to trigger QT interval prolongation and increase the proarrhythmic potential. Yet, how precisely the action of drugs on the cellular level translates into QT interval prolongation on the whole organ level remains insufficiently understood. Here we use machine learning techniques to systematically characterize the effect of 30 common drugs on the QT interval. We combine information from high fidelity three-dimensional human heart simulations with low fidelity one-dimensional cable simulations to build a surrogate model for the QT interval using multi-fidelity Gaussian process regression. Once trained and cross-validated, we apply our surrogate model to perform sensitivity analysis and uncertainty quantification. Our sensitivity analysis suggests that compounds that block the rapid delayed rectifier potassium current I Kr have the greatest prolonging effect of the QT interval, and that blocking the L-type calcium current I CaL and late sodium current I NaL shortens the QT interval. Our uncertainty quantification allows us to propagate the experimental variability from individual block-concentration measurements into the QT interval and reveals that QT interval uncertainty is mainly driven by the variability in I Kr block. In a final validation study, we demonstrate an excellent agreement between our predicted QT interval changes and the changes observed in a randomized clinical trial for the drugs dofetilide, quinidine, ranolazine, and verapamil. We anticipate that both the machine learning methods and the results of this study will have great potential in the efficient development of safer drugs.
Collapse
Affiliation(s)
| | - Kristen Matsuno
- Department of Mechanical Engineering, Stanford University, Stanford, CA 94305, USA
| | - Jiang Yao
- Dassault Systèmes Simulia Corporation, Johnston, RI 02919, USA
| | - Paris Perdikaris
- Mechanical Engineering and Applied Mechanics, University of Pennsylvania, Philadelphia, PA 19104, USA
| | - Ellen Kuhl
- Department of Mechanical Engineering, Stanford University, Stanford, CA 94305, USA
| |
Collapse
|
15
|
Johnston BM, Johnston PR. Sensitivity analysis of ST-segment epicardial potentials arising from changes in ischaemic region conductivities in early and late stage ischaemia. Comput Biol Med 2018; 102:288-299. [DOI: 10.1016/j.compbiomed.2018.06.005] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2018] [Revised: 06/07/2018] [Accepted: 06/07/2018] [Indexed: 11/30/2022]
|
16
|
Quaglino A, Pezzuto S, Koutsourelakis PS, Auricchio A, Krause R. Fast uncertainty quantification of activation sequences in patient-specific cardiac electrophysiology meeting clinical time constraints. Int J Numer Method Biomed Eng 2018; 34:e2985. [PMID: 29577657 DOI: 10.1002/cnm.2985] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/16/2017] [Revised: 01/16/2018] [Accepted: 03/15/2018] [Indexed: 06/08/2023]
Abstract
We present a fast, patient-specific methodology for uncertainty quantification in electrophysiology, aimed at meeting the time constraints of clinical practitioners. We focus on computing the statistics of the activation map, given the uncertainties associated with the conductivity tensor modeling the fiber orientation in the heart. We use a fast parallel solution method implemented on a graphics processing unit for the eikonal approximation, in order to compute the activation map and to sample the random fiber field with correlation on the basis of geodesic distances. While this enables to perform uncertainty quantification studies with a manageable computational effort, the required time frame still exceeds clinically suitable time expectations. In order to reduce it further by 2 orders of magnitude, we rely on Bayesian multifidelity methods. In particular, we propose a low-fidelity model that is patient-specific and free from the additional training cost associated with reduced models. This is achieved by a sound physics-based simplification of the full eikonal model. The low-fidelity output is then corrected by the standard multifidelity framework. In practice, the complete procedure only requires approximately 100 new runs of our eikonal graphics processing unit solver for producing the sought estimates and their associated credible intervals, enabling a full online analysis in less than 5 minutes.
Collapse
Affiliation(s)
- A Quaglino
- Center for Computational Medicine in Cardiology, Institute of Computational Science, Università della Svizzera italiana, Lugano, Switzerland
| | - S Pezzuto
- Center for Computational Medicine in Cardiology, Institute of Computational Science, Università della Svizzera italiana, Lugano, Switzerland
| | | | - A Auricchio
- Center for Computational Medicine in Cardiology, Institute of Computational Science, Università della Svizzera italiana, Lugano, Switzerland
- Division of Cardiology, Fondazione Cardiocentro Ticino, Lugano, Switzerland
| | - R Krause
- Center for Computational Medicine in Cardiology, Institute of Computational Science, Università della Svizzera italiana, Lugano, Switzerland
| |
Collapse
|
17
|
Abstract
Cardiac resynchronization therapy (CRT) is an effective treatment for heart failure (HF) patients with an electrical substrate pathology causing ventricular dyssynchrony. However 40-50% of patients do not respond to treatment. Cardiac modeling of the electrophysiology, electromechanics, and hemodynamics of the heart has been used to study mechanisms behind HF pathology and CRT response. Recently, multi-scale dyssynchronous HF models have been used to study optimal device settings and optimal lead locations, investigate the underlying cardiac pathophysiology, as well as investigate emerging technologies proposed to treat cardiac dyssynchrony. However the breadth of patient and experimental data required to create and parameterize these models and the computational resources required currently limits the use of these models to small patient numbers. In the future, once these technical challenges are overcome, biophysically based models of the heart have the potential to become a clinical tool to aid in the diagnosis and treatment of HF.
Collapse
Affiliation(s)
- Angela W C Lee
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK.
| | | | - Marina Strocchi
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
| | | | - Steven A Niederer
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
| |
Collapse
|
18
|
Lyon A, Mincholé A, Martínez JP, Laguna P, Rodriguez B. Computational techniques for ECG analysis and interpretation in light of their contribution to medical advances. J R Soc Interface 2018; 15:20170821. [PMID: 29321268 PMCID: PMC5805987 DOI: 10.1098/rsif.2017.0821] [Citation(s) in RCA: 74] [Impact Index Per Article: 12.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2017] [Accepted: 12/08/2017] [Indexed: 01/09/2023] Open
Abstract
Widely developed for clinical screening, electrocardiogram (ECG) recordings capture the cardiac electrical activity from the body surface. ECG analysis can therefore be a crucial first step to help diagnose, understand and predict cardiovascular disorders responsible for 30% of deaths worldwide. Computational techniques, and more specifically machine learning techniques and computational modelling are powerful tools for classification, clustering and simulation, and they have recently been applied to address the analysis of medical data, especially ECG data. This review describes the computational methods in use for ECG analysis, with a focus on machine learning and 3D computer simulations, as well as their accuracy, clinical implications and contributions to medical advances. The first section focuses on heartbeat classification and the techniques developed to extract and classify abnormal from regular beats. The second section focuses on patient diagnosis from whole recordings, applied to different diseases. The third section presents real-time diagnosis and applications to wearable devices. The fourth section highlights the recent field of personalized ECG computer simulations and their interpretation. Finally, the discussion section outlines the challenges of ECG analysis and provides a critical assessment of the methods presented. The computational methods reported in this review are a strong asset for medical discoveries and their translation to the clinical world may lead to promising advances.
Collapse
Affiliation(s)
- Aurore Lyon
- Department of Computer Science, British Heart Foundation, Oxford, UK
| | - Ana Mincholé
- Department of Computer Science, British Heart Foundation, Oxford, UK
| | - Juan Pablo Martínez
- Biomedical Signal Interpretation and Computational Simulation (BSICoS) Group, University of Zaragoza, CIBER-BBN, Zaragoza, Spain
| | - Pablo Laguna
- Biomedical Signal Interpretation and Computational Simulation (BSICoS) Group, University of Zaragoza, CIBER-BBN, Zaragoza, Spain
| | - Blanca Rodriguez
- Department of Computer Science, British Heart Foundation, Oxford, UK
| |
Collapse
|