Nonlinear physics of electrical wave propagation in the heart: a review

The beating of the heart is a synchronized contraction of muscle cells (myocytes) that are triggered by a periodic sequence of electrical waves (action potentials) originating in the sino-atrial node and propagating over the atria and the ventricles. Cardiac arrhythmias like atrial and ventricular fibrillation (AF,VF) or ventricular tachycardia (VT) are caused by disruptions and instabilities of these electrical excitations, that lead to the emergence of rotating waves (VT) and turbulent wave patterns (AF,VF). Numerous simulation and experimental studies during the last 20 years have addressed these topics. In this review we focus on the nonlinear dynamics of wave propagation in the heart with an emphasis on the theory of pulses, spirals and scroll waves and their instabilities in excitable media and their application to cardiac modeling. After an introduction into electrophysiological models for action potential propagation, the modeling and analysis of spatiotemporal alternans, spiral and scroll meandering, spiral breakup and scroll wave instabilities like negative line tension and sproing are reviewed in depth and discussed with emphasis on their impact in cardiac arrhythmias.Peer ReviewedPreprin

Three-dimensional cardiac computational modelling: methods, features and applications

[EN] The combination of computational models and biophysical simulations can help to interpret an array of experimental data and contribute to the understanding, diagnosis and treatment of complex diseases such as cardiac arrhythmias. For this reason, three-dimensional (3D) cardiac computational modelling is currently a rising field of research. The advance of medical imaging technology over the last decades has allowed the evolution from generic to patient-specific 3D cardiac models that faithfully represent the anatomy and different cardiac features of a given alive subject. Here we analyse sixty representative 3D cardiac computational models developed and published during the last fifty years, describing their information sources, features, development methods and online availability. This paper also reviews the necessary components to build a 3D computational model of the heart aimed at biophysical simulation, paying especial attention to cardiac electrophysiology (EP), and the existing approaches to incorporate those components. We assess the challenges associated to the different steps of the building process, from the processing of raw clinical or biological data to the final application, including image segmentation, inclusion of substructures and meshing among others. We briefly outline the personalisation approaches that are currently available in 3D cardiac computational modelling. Finally, we present examples of several specific applications, mainly related to cardiac EP simulation and model-based image analysis, showing the potential usefulness of 3D cardiac computational modelling into clinical environments as a tool to aid in the prevention, diagnosis and treatment of cardiac diseases.This work was partially supported by the "VI Plan Nacional de Investigacion Cientifica, Desarrollo e Innovacion Tecnologica" from the Ministerio de Economia y Competitividad of Spain (TIN2012-37546-C03-01 and TIN2011-28067) and the European Commission (European Regional Development Funds - ERDF - FEDER) and by "eTorso project" (GVA/2013-001404) from the Generalitat Valenciana (Spain). ALP is financially supported by the program "Ayudas para contratos predoctorales para la formacion de doctores" from the Ministerio de Economia y Competitividad of Spain (BES-2013-064089).López Pérez, AD.; Sebastián Aguilar, R.; Ferrero De Loma-Osorio, JM. (2015). Three-dimensional cardiac computational modelling: methods, features and applications. BioMedical Engineering OnLine. 14(35):1-31. https://doi.org/10.1186/s12938-015-0033-5S1311435Koushanpour E, Collings W: Validation and dynamic applications of an ellipsoid model of the left ventricle. J Appl Physiol 1966, 21: 1655–61.Ghista D, Sandler H: An analytic elastic-viscoelastic model for the shape and the forces in the left ventricle. J Biomech 1969, 2: 35–47.Janz RF, Grimm AF: Finite-Element Model for the Mechanical Behavior of the Left Ventricle: prediction of deformation in the potassium-arrested rat heart. Circ Res 1972, 30: 244–52.Van den Broek JHJM, Van den Broek MHLM: Application of an ellipsoidal heart model in studying left ventricular contractions. J Biomech 1980, 13: 493–503.Colli Franzone P, Guerri L, Pennacchio M, Taccardi B: Spread of excitation in 3-D models of the anisotropic cardiac tissue. II. Effects of fiber architecture and ventricular geometry. Math Biosci 1998, 147: 131–71.Kerckhoffs RCP, Bovendeerd PHM, Kotte JCS, Prinzen FW, Smits K, Arts T: Homogeneity of cardiac contraction despite physiological asynchrony of depolarization: a model study. Ann Biomed Eng 2003, 31: 536–47.Sermesant M, Moireau P, Camara O, Sainte-Marie J, Andriantsimiavona R, Cimrman R, et al.: Cardiac function estimation from MRI using a heart model and data assimilation: advances and difficulties. Med Image Anal 2006, 10: 642–56.Okajima M, Fujino T, Kobayashi T, Yamada K: Computer simulation of the propagation process in excitation of the ventricles. Circ Res 1968, 23: 203–11.Horan LG, Hand RC, Johnson JC, Sridharan MR, Rankin TB, Flowers NC: A theoretical examination of ventricular repolarization and the secondary T wave. Circ Res 1978, 42: 750–7.Miller WT, Geselowitz DB: Simulation studies of the electrocardiogram. I. The normal heart. Circ Res 1978, 43: 301–15.Vetter FJ, McCulloch AD: Three-dimensional analysis of regional cardiac function: a model of rabbit ventricular anatomy. Prog Biophys Mol Biol 1998, 69: 157–83.Nielsen PMF, LeGrice IJ, Smaill BH, Hunter PJ: Mathematical model of geometry and fibrous structure of the heart. Am J Physiol Heart Circ Physiol 1991, 260: H1365–78.Stevens C, Remme E, LeGrice I, Hunter P: Ventricular mechanics in diastole: material parameter sensitivity. J Biomech 2003, 36: 737–48.Aoki M, Okamoto Y, Musha T, Harumi KI: Three-dimensional simulation of the ventricular depolarization and repolarization processes and body surface potentials: normal heart and bundle branch block. IEEE Trans Biomed Eng 1987, 34: 454–62.Thakor NV, Eisenman LN: Three-dimensional computer model of the heart: fibrillation induced by extrastimulation. Comput Biomed Res 1989, 22: 532–45.Freudenberg J, Schiemann T, Tiede U, Höhne KH: Simulation of cardiac excitation patterns in a three-dimensional anatomical heart atlas. Comput Biol Med 2000, 30: 191–205.Trunk P, Mocnik J, Trobec R, Gersak B: 3D heart model for computer simulations in cardiac surgery. Comput Biol Med 2007, 37: 1398–403.Siregar P, Sinteff JP, Julen N, Le Beux P: An interactive 3D anisotropic cellular automata model of the heart. Comput Biomed Res 1998, 31: 323–47.Harrild DM, Henriquez CS: A computer model of normal conduction in the human atria. Circ Res 2000, 87: e25–36.Bodin ON, Kuz’min AV: Synthesis of a realistic model of the surface of the heart. Biomed Eng (NY) 2006, 40: 280–3.Ruiz-Villa CA, Tobón C, Rodríguez JF, Ferrero JM, Hornero F, Saíz J: Influence of atrial dilatation in the generation of re-entries caused by ectopic activity in the left atrium. Comput Cardiol 2009, 36: 457–60.Blanc O, Virag N, Vesin JM, Kappenberger L: A computer model of human atria with reasonable computation load and realistic anatomical properties. IEEE Trans Biomed Eng 2001, 48: 1229–37.Zemlin CW, Herzel H, Ho SY, Panfilov AV: A realistic and efficient model of excitation propagation in the human atria. In Comput Simul Exp Assess Card Electrophysiol. Edited by: Virag N, Kappenberger L, Blanc O. Futura Publishing Company, Inc, Arkmonk, New York; 2001:29–34.Seemann G, Höper C, Sachse FB, Dössel O, Holden AV, Zhang H: Heterogeneous three-dimensional anatomical and electrophysiological model of human atria. Philos Trans R Soc A Math Phys Eng Sci 2006, 364: 1465–81.Zhao J, Butters TD, Zhang H, LeGrice IJ, Sands GB, Smaill BH: Image-based model of atrial anatomy and electrical activation: a computational platform for investigating atrial arrhythmia. IEEE Trans Med Imaging 2013, 32: 18–27.Creswell LL, Wyers SG, Pirolo JS, Perman WH, Vannier MW, Pasque MK: Mathematical modeling of the heart using magnetic resonance imaging. IEEE Trans Med Imaging 1992, 11: 581–9.Lorange M, Gulrajani RM: A computer heart model incorporating anisotropic propagation: I. Model construction and simulation of normal activation. J Electrocardiol 1993, 26: 245–61.Winslow RL, Scollan DF, Holmes A, Yung CK, Zhang J, Jafri MS: Electrophysiological modeling of cardiac ventricular function: from cell to organ. Annu Rev Biomed Eng 2000, 2: 119–55.Virag N, Jacquemet V, Henriquez CS, Zozor S, Blanc O, Vesin JM, et al.: Study of atrial arrhythmias in a computer model based on magnetic resonance images of human atria. Chaos 2002, 12: 754–63.Helm PA, Tseng HJ, Younes L, McVeigh ER, Winslow RL: Ex vivo 3D diffusion tensor imaging and quantification of cardiac laminar structure. Magn Reson Med 2005, 54: 850–9.Arevalo HJ, Helm PA, Trayanova NA: Development of a model of the infarcted canine heart that predicts arrhythmia generation from specific cardiac geometry and scar distribution. Comput Cardiol 2008, 35: 497–500.Plotkowiak M, Rodriguez B, Plank G, Schneider JE, Gavaghan D, Kohl P, et al.: High performance computer simulations of cardiac electrical function based on high resolution MRI datasets. In Int Conf Comput Sci 2008, LNCS 5101. Springer–Verlag, Berlin Heidelberg; 2008:571–80.Heidenreich EA, Ferrero JM, Doblaré M, Rodríguez JF: Adaptive macro finite elements for the numerical solution of monodomain equations in cardiac electrophysiology. Ann Biomed Eng 2010, 38: 2331–45.Gurev V, Lee T, Constantino J, Arevalo H, Trayanova NA: Models of cardiac electromechanics based on individual hearts imaging data: Image-based electromechanical models of the heart. Biomech Model Mechanobiol 2011, 10: 295–306.Deng D, Jiao P, Ye X, Xia L: An image-based model of the whole human heart with detailed anatomical structure and fiber orientation. Comput Math Methods Med 2012, 2012: 16.Aslanidi OV, Nikolaidou T, Zhao J, Smaill BH, Gilbert SH, Holden AV, et al.: Application of micro-computed tomography with iodine staining to cardiac imaging, segmentation, and computational model development. IEEE Trans Med Imaging 2013, 32: 8–17.Haddad R, Clarysse P, Orkisz M, Croisille P, Revel D, Magnin IE: A realistic anthropomorphic numerical model of the beating heart. In Funct Imaging Model Heart 2005, LNCS 3504. Springer–Verlag, Berlin Heidelberg; 2005:384–93.Appleton B, Wei Q, Liu N, Xia L, Crozier S, Liu F, et al.: An electrical heart model incorporating real geometry and motion. In 27th Annu Int Conf Eng Med Biol Soc (IEEE-EMBS 2005). IEEE, Shanghai, China; 2006:345–8.Niederer S, Rhode K, Razavi R, Smith N: The importance of model parameters and boundary conditions in whole organ models of cardiac contraction. In Funct Imaging Model Heart 2009, LNCS 5528. Springer–Verlag, Berlin Heidelberg; 2009:348–56.Yang G, Toumoulin C, Coatrieux JL, Shu H, Luo L, Boulmier D: A 3D static heart model from a MSCT data set. In 27th Annu Int Conf IEEE Eng Med Biol Soc (IEEE-EMBS 2005). IEEE, Shangai, China; 2006:5499–502.Romero D, Sebastian R, Bijnens BH, Zimmerman V, Boyle PM, Vigmond EJ, et al.: Effects of the purkinje system and cardiac geometry on biventricular pacing: a model study. Ann Biomed Eng 2010, 38: 1388–98.Lorenzo-Valdés M, Sanchez-Ortiz GI, Mohiaddin R, Rueckert D: Atlas-based segmentation and tracking of 3D cardiac MR images using non-rigid registration. In Med Image Comput Comput Assist Interv 2002, LNCS 2488. Springer–Verlag, Berlin Heidelberg; 2002:642–50.Ordas S, Oubel E, Sebastian R, Frangi AF: Computational anatomy atlas of the heart. In 5th Int Symp Image Signal Process Anal (ISPA 2007). IEEE, Istanbul, Turkey; 2007:338–42.Burton RAB, Plank G, Schneider JE, Grau V, Ahammer H, Keeling SL, et al.: Three-dimensional models of individual cardiac histoanatomy: tools and challenges. Ann N Y Acad Sci 2006, 1080: 301–19.Plank G, Burton RAB, Hales P, Bishop M, Mansoori T, Bernabeu MO, et al.: Generation of histo-anatomically representative models of the individual heart: tools and application. Philos Trans R Soc A Math Phys Eng Sci 2009, 367: 2257–92.Bishop MJ, Plank G, Burton RAB, Schneider JE, Gavaghan DJ, Grau V, et al.: Development of an anatomically detailed MRI-derived rabbit ventricular model and assessment of its impact on simulations of electrophysiological function. Am J Physiol - Heart Circ Physiol 2010, 298: H699–718.Ecabert O, Peters J, Schramm H, Lorenz C, von Berg J, Walker MJ, et al.: Automatic model-based segmentation of the heart in CT images. IEEE Trans Med Imaging 2008, 27: 1189–201.Ecabert O, Peters J, Walker MJ, Ivanc T, Lorenz C, von Berg J, et al.: Segmentation of the heart and great vessels in CT images using a model-based adaptation framework. Med Image Anal 2011, 15: 863–76.Schulte RF, Sands GB, Sachse FB, Dössel O, Pullan AJ: Creation of a human heart model and its customisation using ultrasound images. Biomed Tech Eng 2001, 46: 26–8.Wenk JF, Zhang Z, Cheng G, Malhotra D, Acevedo-Bolton G, Burger M, et al.: First finite element model of the left ventricle with mitral valve: insights into ischemic mitral regurgitation. Ann Thorac Surg 2010, 89: 1546–53.Frangi AF, Rueckert D, Schnabel JA, Niessen WJ: Automatic construction of multiple-object three-dimensional statistical shape models: application to cardiac modeling. IEEE Trans Med Imaging 2002, 21: 1151–66.Hoogendoorn C, Duchateau N, Sánchez-Quintana D, Whitmarsh T, Sukno FM, De Craene M, et al.: A high-resolution atlas and statistical model of the human heart from multislice CT. IEEE Trans Med Imaging 2013, 32: 28–44.Vadakkumpadan F, Rantner LJ, Tice B, Boyle P, Prassl AJ, Vigmond E, et al.: Image-based models of cardiac structure with applications in arrhythmia and defibrillation studies. J Electrocardiol 2009, 42: 157.Perperidis D, Mohiaddin R, Rueckert D: Construction of a 4D statistical atlas of the cardiac anatomy and its use in classification. In Med Image Comput Comput Interv 2005, LNCS 3750. Springer–Verlag, Berlin Heidelberg; 2005:402–10.Lötjönen J, Kivistö S, Koikkalainen J, Smutek D, Lauerma K: Statistical shape model of atria, ventricles and epicardium from short- and long-axis MR images. Med Image Anal 2004, 8: 371–86.Lorenz C, von Berg J: A comprehensive shape model of the heart. Med Image Anal 2006, 10: 657–70.Mansoori T, Plank G, Burton R, Schneider J, Khol P, Gavaghan D, et al.: An iterative method for registration of high-resolution cardiac histoanatomical and MRI images. In 4th IEEE Int Symp Biomed Imaging: From Nano to Macro (ISBI 2007). IEEE, Arlington, VA (USA); 2007:572–5.Gibb M, Burton RAB, Bollensdorff C, Afonso C, Mansoori T, Schotten U, et al.: Resolving the three-dimensional histology of the heart. In Comput Methods Syst Biol - Lect Notes Comput Sci 7605. Springer, Berlin Heidelberg; 2012:2–16.Burton RAB, Lee P, Casero R, Garny A, Siedlecka U, Schneider JE, et al.: Three-dimensional histology: tools and application to quantitative assessment of cell-type distribution in rabbit heart. Europace 2014,16(Suppl 4):iv86–95.Niederer SA, Shetty AK, Plank G, Bostock J, Razavi R, Smith NP, et al.: Biophysical modeling to simulate the response to multisite left ventricular stimulation using a quadripolar pacing lead. Pacing Clin Electrophysiol 2012, 35: 204–14.Weese J, Groth A, Nickisch H, Barschdorf H, Weber FM, Velut J, et al.: Generating anatomical models of the heart and the aorta from medical images for personalized physiological simulations. Med Biol Eng Comput 2013, 51: 1209–19.Gibb M, Bishop M, Burton R, Kohl P, Grau V, Plank G, et al.: The role of blood vessels in rabbit propagation dynamics and cardiac arrhythmias. In Funct Imaging Model Heart - FIMH 2009, LNCS 5528. Springer, Berlin Heidelberg; 2009:268–76.Prassl AJ, Kickinger F, Ahammer H, Grau V, Schneider JE, Hofer E, et al.: Automatically generated, anatomically accurate meshes for cardiac electrophysiology problems. IEEE Trans Biomed Eng 2009, 56: 1318–30.Dux-Santoy L, Sebastian R, Felix-Rodriguez J, Ferrero JM, Saiz J: Interaction of specialized cardiac conduction system with antiarrhythmic drugs: a simulation study. IEEE Trans Biomed Eng 2011, 58: 3475–8.Lamata P, Niederer S, Nordsletten D, Barber DC, Roy I, Hose DR, et al.: An accurate, fast and robust method to generate patient-specific cubic Hermite meshes. Med Image Anal 2011, 15: 801–13.Pathmanathan P, Cooper J, Fletcher A, Mirams G, Murray P, Osborne J, et al.: A computational study of discrete mechanical tissue models. Phys Biol 2009, 6: 036001.Niederer SA, Kerfoot E, Benson AP, Bernabeu MO, Bernus O, Bradley C, et al.: Verification of cardiac tissue electrophysiology simulators using an N-version benchmark. Philos Trans R Soc A Math Phys Eng Sci 2011, 369: 4331–51.Ten Tusscher KHWJ, Panfilov AV: Cell model for efficient simulation of wave propagation in human ventricular tissue under normal and pathological conditions. Phys Med Biol 2006, 51: 6141–56.LeGrice I, Smaill B, Chai L, Edgar S, Gavin J, Hunter P: Laminar structure of the heart: ventricular myocyte arrangement and connective tissue architecture in the dog. Am J Physiol Heart Circ Physiol 1995, 269: H571–82.Anderson RH, Smerup M, Sanchez-Quintana D, Loukas M, Lunkenheimer PP: The three-dimensional arrangement of the myocytes in the ventricular walls. Clin Anat 2009, 22: 64–76.Clerc L: Directional differences of impulse spread in trabecular muscle from mammalian heart. J Physiol 1976, 255: 335–46.Streeter DD Jr, Spotnitz HM, Patel DP, Ross J Jr, Sonnenblick EH: Fiber orientation in the canine left ventricle during diastole and systole. Circ Res 1969, 24: 339–47.Scollan D, Holmes A, Winslow R, Forder J: Histological validation of myocardial microstructure obtained from diffusion tensor magnetic resonance imaging. Am J Physiol Heart Circ Physiol 1998, 275: H2308–18.Hsu EW, Muzikant AL, Matulevicius SA, Penland RC, Henriquez CS: Magnetic resonance myocardial fiber-orientation mapping with direct histological correlation. Am J Physiol Heart Circ Physiol 1998, 274: H1627–34.Holmes AA, Scollan DF, Winslow RL: Direct histological validation of diffusion tensor MRI in formaldehyde-fixed myocardium. Magn Reson Med 2000, 44: 157–61.Sermesant M, Forest C, Pennec X, Delingette H, Ayache N: Deformable biomechanical models: application to 4D cardiac image analysis. Med Image Anal 2003, 7: 475–88.Peyrat JM, Sermesant M, Pennec X, Delingette H, Xu C, McVeigh ER, et al.: A computational framework for the statistical analysis of cardiac diffusion tensors: application to a small database of canine hearts. IEEE Trans Med Imaging 2007, 26: 1500–14.Toussaint N, Sermesant M, Stoeck CT, Kozerke S, Batchelor PG: In vivo human 3D cardiac fibre architecture: reconstruction using curvilinear interpolation of diffusion tensor images. Med Image Comput Comput Assist Interv 2010,13(Pt 1):418–25.Toussaint N, Stoeck CT, Schaeffter T, Kozerke S, Sermesant M, Batchelor PG: In vivo human cardiac fibre architecture estimation using shape-based diffusion tensor processing. Med Image Anal 2013, 17: 1243–55.Bishop MJ, Hales P, Plank G, Gavaghan DJ, Scheider J, Grau V: Comparison of rule-based and DTMRI-derived fibre architecture in a whole rat ventricular computational model. In Funct Imaging Model Heart 2009, LNCS 5528. Springer–Verlag, Berlin Heidelberg; 2009:87–96.Bayer JD, Blake RC, Plank G, Trayanova NA: A novel rule-based algorithm for assigning myocardial fiber orientation to computational heart models. Ann Biomed Eng 2012, 40: 2243–54.Dobrzynski H, Anderson RH, Atkinson A, Borbas Z, D’Souza A, Fraser JF, et al.: Structure, function and clinical relevance of the cardiac conduction system, including the atrioventricular ring and outflow tract tissues. Pharmacol Ther 2013, 139: 260–88.Tranum-Jensen J, Wilde AA, Vermeulen JT, Janse MJ: Morphology of electrophysiologically identified junctions between Purkinje fibers and ventricular muscle in rabbit and pig hearts. Circ Res 1991, 69: 429–37.Boyle PM, Deo M, Plank G, Vigmond EJ: Purkinje-mediated effects in the response of quiescent ventricles to defibrillation shocks. Ann Biomed Eng 2010, 38: 456–68.Behradfar E, Nygren A, Vigmond EJ: The role of Purkinje-myocardial coupling during ventricular arrhythmia: a modeling study. PLoS One 2014., 9: Article ID e88000DiFrancesco D, Noble D: A model of cardiac electrical activity incorporating ionic pumps and concentration changes. Philos Trans R Soc B Biol Sci 1985, 307: 353–98.Stewart P, Aslanidi OV, Noble D, Noble PJ, Boyett MR, Zhang H: Mathematical models of the electrical action potential of Purkinje fibre cells. Philos Trans R Soc A Math Phys Eng Sci 2009, 367: 2225–55.Li P, Rudy Y: A model of canine purkinje cell electrophysiology and Ca(2+) cycling: rate dependence, triggered activity, and comparison to ventricular myocytes. Circ Res 2011, 109: 71–9.Chinchapatnam P, Rhode KS, Ginks M, Mansi T, Peyrat JM, Lambiase P, et al.: Estimation of volumetric myocardial apparent conductivity from endocardial electro-anatomical mapping. In 31st Annu Int Conf IEEE Eng Med Biol Soc (EMBC 2009). IEEE, Minneapolis, MN (USA); 2009:2907–10.Durrer D, Van Dam RT, Freud GE, Janse MJ, Meijler FL, Arzbaecher RC: Total excitation of the isolated human heart. Circulation 1970, 41: 899–912.Pollard AE, Barr RC: Computer simulations of activation in an anatomically based model of the human ventricular conduction system. IEEE Trans Biomed Eng 1991, 38: 982–96.Abboud S, Berenfeld O, Sadeh D: Simulation of high-resolution QRS complex using a ventricular model with a fractal conduction system. Effects of ischemia on high-frequency QRS potentials. Circ Res 1991, 68: 1751–60.Sebastian R, Zimmerman V, Romero D, Sanchez-Quintana D, Frangi AF: Characterization and modeling of the peripheral cardiac conduction system. IEEE Trans Med Imaging 2013, 32: 45–55.Bordas R, Gillow K, Lou Q, Efimov IR, Gavaghan D, Kohl P, et al.: Rabbit-specific ventricular model of cardiac electrophysiological function including specialized conduction system. Prog Biophys Mol Biol 2011, 107: 90–100.Stephenson RS, Boyett MR, Hart G, Nikolaidou T, Cai X, Corno AF, et al.: Contrast enhanced micro-computed tomography resolves the 3-dimensional morphology of the cardiac conduction system in mammalian hearts. PLoS One 2012., 7: Article ID e35299Berenfeld O, Jalife J: Purkinje-Muscle reentry as a mechanism of polymorphic ventricular arrhythmias in a 3-dimensional model of the ventricles. Circ Res 1998, 82: 1063–77.Azzouzi A, Coudière Y, Turpault R, Zemzemi N: A mathematical model of the Purkinje-muscle junctions. Math Biosci Eng MBE 2011, 8: 915–30.Dux-Santoy L, Sebastian R, Rodriguez JF, Ferrero JM: Modeling the different sections of the cardiac conduction system to obtain realistic electrocardiograms. In 35th Annu Int Conf IEEE Eng Med Biol Soc (EMBC 2013). IEEE, Osaka, Japan; 2013:6846–9.Cardenes R, Sebastian R, Berruezo A, Camara O: Inverse

Fractional diffusion models of cardiac electrical propagation: role of structural heterogeneity in dispersion of repolarization

Structural heterogeneity constitutes one of the main substrates influencing impulse propagation in living tissues. In cardiac muscle, improved understanding on its role is key to advancing our interpretation of cell-to-cell coupling, and how tissue structure modulates electrical propagation and arrhythmogenesis in the intact and diseased heart. We propose fractional diffusion models as a novel mathematical description of structurally heterogeneous excitable media, as a mean of representing the modulation of the total electric field by the secondary electrical sources associated with tissue inhomogeneities. Our results, validated against in-vivo human recordings and experimental data of different animal species, indicate that structural heterogeneity underlies many relevant characteristics of cardiac propagation, including the shortening of action potential duration along the activation pathway, and the progressive modulation by premature beats of spatial patterns of dispersion of repolarization. The proposed approach may also have important implications in other research fields involving excitable complex media

Novel Cardiac Mapping Approaches and Multimodal Techniques to Unravel Multidomain Dynamics of Complex Arrhythmias Towards a Framework for Translational Mechanistic-Based Therapeutic Strategies

[ES] Las arritmias cardíacas son un problema importante para los sistemas de salud en el mundo desarrollado debido a su alta incidencia y prevalencia a medida que la población envejece. La fibrilación auricular (FA) y la fibrilación ventricular (FV) se encuentran entre las arritmias más complejas observadas en la práctica clínica. Las consecuencias clínicas de tales alteraciones arrítmicas incluyen el desarrollo de eventos cardioembólicos complejos en la FA, y repercusiones dramáticas debido a procesos fibrilatorios sostenidos que amenazan la vida infringiendo daño neurológico tras paro cardíaco por FV, y que pueden provocar la muerte súbita cardíaca (MSC). Sin embargo, a pesar de los avances tecnológicos de las últimas décadas, sus mecanismos intrínsecos se comprenden de forma incompleta y, hasta la fecha, las estrategias terapéuticas carecen de una base mecanicista suficiente y poseen bajas tasas de éxito. Entre los mecanismos implicados en la inducción y perpetuación de arritmias cardíacas, como la FA, se cree que las dinámicas de las fuentes focales y reentrantes de alta frecuencia, en sus diferentes modalidades, son las fuentes primarias que mantienen la arritmia. Sin embargo, se sabe poco sobre los atractores, así como, de la dinámica espacio-temporal de tales fuentes fibrilatorias primarias, específicamente, las fuentes focales o rotacionales dominantes que mantienen la arritmia. Por ello, se ha desarrollado una plataforma computacional, para comprender los factores (activos, pasivos y estructurales) determinantes, y moduladores de dicha dinámica. Esto ha permitido establecer un marco para comprender la compleja dinámica de los rotores con énfasis en sus propiedades deterministas para desarrollar herramientas basadas en los mecanismos para ayuda diagnóstica y terapéutica. Comprender los procesos fibrilatorios es clave para desarrollar marcadores y herramientas fisiológica- y clínicamente relevantes para la ayuda de diagnóstico temprano. Específicamente, las propiedades espectrales y de tiempo-frecuencia de los procesos fibrilatorios han demostrado resaltar el comportamiento determinista principal de los mecanismos intrínsecos subyacentes a las arritmias y el impacto de tales eventos arrítmicos. Esto es especialmente relevante para determinar el pronóstico temprano de los supervivientes comatosos después de un paro cardíaco debido a fibrilación ventricular (FV). Las técnicas de mapeo electrofisiológico, el mapeo eléctrico y óptico cardíaco, han demostrado ser recursos muy valiosos para dar forma a nuevas hipótesis y desarrollar nuevos enfoques mecanicistas y estrategias terapéuticas mejoradas. Esta tecnología permite además el trabajo multidisciplinar entre clínicos y bioingenieros, para el desarrollo y validación de dispositivos y metodologías para identificar biomarcadores multi-dominio que permitan rastrear con precisión la dinámica de las arritmias identificando fuentes dominantes y atractores con alta precisión para ser dianas de estrategias terapeúticas innovadoras. Es por ello que uno de los objetivos fundamentales ha sido la implantación y validación de nuevos sistemas de mapeo en distintas configuraciones que sirvan de plataforma de desarrollo de nuevas estrategias terapeúticas. Aunque el mapeo panorámico es el método principal y más completo para rastrear simultáneamente biomarcadores electrofisiológicos, su adopción por la comunidad científica es limitada principalmente debido al coste elevado de la tecnología. Aprovechando los avances tecnológicos recientes, nos hemos enfocado en desarrollar, y validar, sistemas de mapeo óptico de alta resolución para registro panorámico cardíaco, utilizando modelos clínicamente relevantes para la investigación básica y la bioingeniería.[CA] Les arítmies cardíaques són un problema important per als sistemes de salut del món desenvolupat a causa de la seva alta incidència i prevalença a mesura que la població envelleix. La fibril·lació auricular (FA) i la fibril·lació ventricular (FV), es troben entre les arítmies més complexes observades a la pràctica clínica. Les conseqüències clíniques d'aquests trastorns arítmics inclouen el desenvolupament d'esdeveniments cardioembòlics complexos en FA i repercussions dramàtiques a causa de processos fibril·latoris sostinguts que posen en perill la vida amb danys neurològics posteriors a la FV, que condueixen a una aturada cardíaca i a la mort cardíaca sobtada (SCD). Tanmateix, malgrat els avanços tecnològics de les darreres dècades, els seus mecanismes intrínsecs s'entenen de forma incompleta i, fins a la data, les estratègies terapèutiques no tenen una base mecanicista suficient i tenen baixes taxes d'èxit. La majoria dels avenços en el desenvolupament de biomarcadors òptims i noves estratègies terapèutiques en aquest camp provenen de tècniques valuoses en la investigació de mecanismes d'arítmia. Entre els mecanismes implicats en la inducció i perpetuació de les arítmies cardíaques, es creu que les fonts primàries subjacents a l'arítmia són les fonts focals reingressants d'alta freqüència dinàmica i AF, en les seves diferents modalitats. Tot i això, se sap poc sobre els atractors i la dinàmica espaciotemporal d'aquestes fonts primàries fibril·ladores, específicament les fonts rotacionals o focals dominants que mantenen l'arítmia. Per tant, s'ha desenvolupat una plataforma computacional per entendre determinants actius, passius, estructurals i moduladors d'aquestes dinàmiques. Això va permetre establir un marc per entendre la complexa dinàmica multidomini dels rotors amb ènfasi en les seves propietats deterministes per desenvolupar enfocaments mecanicistes per a l'ajuda i la teràpia diagnòstiques. La comprensió dels processos fibril·latoris és clau per desenvolupar puntuacions i eines rellevants fisiològicament i clínicament per ajudar al diagnòstic precoç. Concretament, les propietats espectrals i de temps-freqüència dels processos fibril·latoris han demostrat destacar un comportament determinista important dels mecanismes intrínsecs subjacents a les arítmies i l'impacte d'aquests esdeveniments arítmics. Mitjançant coneixements previs, processament de senyals, tècniques d'aprenentatge automàtic i anàlisi de dades, es va desenvolupar una puntuació de risc mecanicista a la aturada cardíaca per FV. Les tècniques de cartografia òptica cardíaca i electrofisiològica han demostrat ser recursos inestimables per donar forma a noves hipòtesis i desenvolupar nous enfocaments mecanicistes i estratègies terapèutiques. Aquesta tecnologia ha permès durant molts anys provar noves estratègies terapèutiques farmacològiques o ablatives i desenvolupar mètodes multidominis per fer un seguiment precís de la dinàmica d'arrímies que identifica fonts i atractors dominants. Tot i que el mapatge panoràmic és el mètode principal per al seguiment simultani de paràmetres electrofisiològics, la seva adopció per part de la comunitat multidisciplinària d'investigació cardiovascular està limitada principalment pel cost de la tecnologia. Aprofitant els avenços tecnològics recents, ens centrem en el desenvolupament i la validació de sistemes de mapes òptics de baix cost per a imatges panoràmiques mitjançant models clínicament rellevants per a la investigació bàsica i la bioenginyeria.[EN] Cardiac arrhythmias are a major problem for health systems in the developed world due to their high incidence and prevalence as the population ages. Atrial fibrillation (AF) and ventricular fibrillation (VF), are amongst the most complex arrhythmias seen in the clinical practice. Clinical consequences of such arrhythmic disturbances include developing complex cardio-embolic events in AF, and dramatic repercussions due to sustained life-threatening fibrillatory processes with subsequent neurological damage under VF, leading to cardiac arrest and sudden cardiac death (SCD). However, despite the technological advances in the last decades, their intrinsic mechanisms are incompletely understood, and, to date, therapeutic strategies lack of sufficient mechanistic basis and have low success rates. Most of the progress for developing optimal biomarkers and novel therapeutic strategies in this field has come from valuable techniques in the research of arrhythmia mechanisms. Amongst the mechanisms involved in the induction and perpetuation of cardiac arrhythmias such AF, dynamic high-frequency re-entrant and focal sources, in its different modalities, are thought to be the primary sources underlying the arrhythmia. However, little is known about the attractors and spatiotemporal dynamics of such fibrillatory primary sources, specifically dominant rotational or focal sources maintaining the arrhythmia. Therefore, a computational platform for understanding active, passive and structural determinants, and modulators of such dynamics was developed. This allowed stablishing a framework for understanding the complex multidomain dynamics of rotors with enphasis in their deterministic properties to develop mechanistic approaches for diagnostic aid and therapy. Understanding fibrillatory processes is key to develop physiologically and clinically relevant scores and tools for early diagnostic aid. Specifically, spectral and time-frequency properties of fibrillatory processes have shown to highlight major deterministic behaviour of intrinsic mechanisms underlying the arrhythmias and the impact of such arrhythmic events. Using prior knowledge, signal processing, machine learning techniques and data analytics, we aimed at developing a reliable mechanistic risk-score for comatose survivors of cardiac arrest due to VF. Cardiac optical mapping and electrophysiological mapping techniques have shown to be unvaluable resources to shape new hypotheses and develop novel mechanistic approaches and therapeutic strategies. This technology has allowed for many years testing new pharmacological or ablative therapeutic strategies, and developing multidomain methods to accurately track arrhymia dynamics identigying dominant sources and attractors. Even though, panoramic mapping is the primary method for simultaneously tracking electrophysiological parameters, its adoption by the multidisciplinary cardiovascular research community is limited mainly due to the cost of the technology. Taking advantage of recent technological advances, we focus on developing and validating low-cost optical mapping systems for panoramic imaging using clinically relevant models for basic research and bioengineering.Calvo Saiz, CJ. (2022). Novel Cardiac Mapping Approaches and Multimodal Techniques to Unravel Multidomain Dynamics of Complex Arrhythmias Towards a Framework for Translational Mechanistic-Based Therapeutic Strategies [Tesis doctoral]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/182329TESI

Arrhythmia mechanisms in acute ischaemia and chronic infarction in rabbit heart

In this thesis, a method for studying the electrophysiological consequences of acute regional ischaemia in rabbit heart was established using a combination of a novel snare technique and optical mapping. The purpose of this approach was to discover the mechanistic link between acute coronary infarction and the occurrence of arrhythmias. The electrophysiology of the epicardial surface of isolated hearts was examined using the voltage sensitive dye RH237 and optical action potentials were recorded from a 13x13mm area of left ventricular epicardium using a 16x16 element Hamamatsu photodiode array. Contraction motion artefacts were practically eliminated with blebbistatin (5µM). An alternative mechanical uncoupler, BDM, was found to be not suitable for the study of arrhythmic behaviour associated with ischaemia. After occlusion of the left coronary artery, a progressive reduction in action potential duration (APD), and slowing of upstroke was observed in an area of the left ventricle anterior surface, accompanied by ECG S-T segment elevation. These effects were reversed when the coronary artery occlusion was released. Ligation (duration 12-15mins) caused a decrease in APD50 (APD at 50% repolarisation), in the zone of reduced perfusion, from 141±5.2ms to 53.3±9.3ms (mean±SEM, n=10 hearts, P<0.001). After ligation was reversed and full perfusion restored, APD50 returned to normal values (149±7.0ms, n.s.). Trise (action potential rise time from 10-90% depolarisation) increased from 7.2±1.0ms to 15.8±2.8ms (P<0.01). In the non-infarcted area of myocardium, no significant changes in APD50 (147±7.0ms vs. 147±8.1ms) or Trise (6.4±0.4ms vs 8.8±1.4ms) were observed during occlusion. T-wave alternans behaviour was observed frequently during local ischaemia and associated with alternans of optical action potentials (OAPs) in the ischaemic border zone (BZ) and in ischaemic zone (IZ). T-wave alternans amplitude was not maintained during local ischaemia but OAPs continued to show alternating behaviour. Arrhythmias (VT and VF) were common when conduction block occurred at the interface between the normal and ischaemic zone, but arrhythmias were absent when conduction into the IZ was retained. This observation suggests that the conduction block was the crucial precipitating event for the generation of arrhythmias. Acute local ischaemia was also imposed in a heart with an existing infarct scar to examine the effects of pre-existing ischaemic damage. The incidence of arrhythmias was similar to that observed in the absence of an infarct scar indicating that pre-existing damage did not predispose the heart to arrhythmias. Global ischaemic challenges, both low flow and zero flow produced similar reductions in APD and rise time and were followed by arrhythmias, but the associated changes in the ECG were complex and could not be easily interpreted. Significant temporal variability in electrophysiology was observed in global ischaemia, but absent in the local ischaemic challenge. The underlying mechanisms of these temporal flucuations in cardiac electrophysiology may be dictated by either cellular metabolism or fluctuations in coronary flow. Long-term local ischaemia (~60mins) did not reveal a second phase of arrhythmias after 40-45mins as observed in other animal models, and nor were there signs of significant further electrophysiological changes as a consequence of the additional period of local ischaemia

Autonomic Nervous System and Neurocardiac Physiopathology

The autonomic nervous system regulates multiple physiological functions; how distinct neurons in peripheral autonomic and intrathoracic ganglia communicate remains to be established. Increasing focus is being paid to functionality of the neurocardiac axis and crosstalk between the intrinsic nervous system and diverse organ systems. Current findings indicate that progression of cardiovascular disease comprises peripheral and central aspects of the cardiac nervous system hierarchy. Indeed, autonomic neuronal dysfunction is known to participate in arrhythmogenesis and sudden cardiac death; diverse interventions (pharmacological, non-pharmacological) that affect neuronal remodeling in the heart following injury caused by cardiovascular disease (congestive heart failure, etc.) or acute myocardial infarction are being investigated. Herein we examine recent findings from clinical and animal studies on the role of the intrinsic cardiac nervous system on regulation of myocardial perfusion and the consequences of cardiac injury. We also discuss different interventions that target the autonomic nervous system, stimulate neuronal remodeling and adaptation, and thereby optimize patient outcomes

Evolution of spiral and scroll waves of excitation in a mathematical model of ischaemic border zone

Abnormal electrical activity from the boundaries of ischemic cardiac tissue is recognized as one of the major causes in generation of ischemia-reperfusion arrhythmias. Here we present theoretical analysis of the waves of electrical activity that can rise on the boundary of cardiac cell network upon its recovery from ischaemia-like conditions. The main factors included in our analysis are macroscopic gradients of the cell-to-cell coupling and cell excitability and microscopic heterogeneity of individual cells. The interplay between these factors allows one to explain how spirals form, drift together with the moving boundary, get transiently pinned to local inhomogeneities, and finally penetrate into the bulk of the well-coupled tissue where they reach macroscopic scale. The asymptotic theory of the drift of spiral and scroll waves based on response functions provides explanation of the drifts involved in this mechanism, with the exception of effects due to the discreteness of cardiac tissue. In particular, this asymptotic theory allows an extrapolation of 2D events into 3D, which has shown that cells within the border zone can give rise to 3D analogues of spirals, the scroll waves. When and if such scroll waves escape into a better coupled tissue, they are likely to collapse due to the positive filament tension. However, our simulations have shown that such collapse of newly generated scrolls is not inevitable and that under certain conditions filament tension becomes negative, leading to scroll filaments to expand and multiply leading to a fibrillation-like state within small areas of cardiac tissue.Comment: 26 pages, 13 figures, appendix and 2 movies, as accepted to PLoS ONE 2011/08/0

Highly trabeculated structure of the human endocardium underlies asymmetrical response to low-energy monophasic shocks

Novel low-energy defibrillation therapies are thought to be driven by virtual-electrodes (VEs), due to the interaction of applied monophasic electric shocks with fine-scale anatomical structures within the heart. Significant inter-species differences in the cardiac (micro)-anatomy exist, however, particularly with respect to the degree of endocardial trabeculations, which may underlie important differences in response to low-energy defibrillation protocols. Understanding the interaction of monophasic electric fields with the specific human micro-anatomy is therefore imperative in facilitating the translation and optimisation of these promising experimental therapies to the clinic. In this study, we sought to investigate how electric fields from implanted devices interact with the highly trabeculated human endocardial surface to better understand shock success in order to help optimise future clinical protocols. A bi-ventricular human computational model was constructed from high resolution (350 μm) ex-vivo MR data, including anatomically accurate endocardial structures. Monophasic shocks were applied between a basal right ventricular catheter and an exterior ground. Shocks of varying strengths were applied with both anodal [positive right ventricle (RV) electrode] and cathodal (negative RV electrode) polarities at different states of tissue refractoriness and during induced arrhythmias. Anodal shocks induced isolated positive VEs at the distal side of “detached” trabeculations, which rapidly spread into hyperpolarised tissue on the surrounding endocardial surfaces following the shock. Anodal shocks thus depolarised more tissue 10 ms after the shock than cathodal shocks where the propagation of activation from VEs induced on the proximal side of “detached” trabeculations was prevented due to refractory endocardium. Anodal shocks increased arrhythmia complexity more than cathodal shocks during failed anti-arrhythmia shocks. In conclusion, multiple detached trabeculations in the human ventricle interact with anodal stimuli to induce multiple secondary sources from VEs, facilitating more rapid shock-induced ventricular excitation compared to cathodal shocks. Such a mechanism may help explain inter-species differences in response to shocks and help to develop novel defibrillation strategies