Electrophysiology Model for a Human Heart with Ischemic Scar and Realistic Purkinje Network

The role of Purkinje fibres in the onset of arrhythmias is controversial and computer simulations may shed light on possible arrhythmic mechanisms involving the Purkinje fibres. However, few computational modelling studies currently include a detailed Purkinje network as part of the model. We present a coupled Purkinje-myocardium electrophysiology model that includes an explicit model for the ischemic scar plus a detailed Purkinje network, and compare simulated activation times to those obtained by electro-anatomical mapping in vivo during sinus rhythm pacing. The results illustrate the importance of using sufficiently dense Purkinje networks in patient-specific studies to capture correctly the myocardial early activation that may be influenced by surviving Purkinje fibres in the infarct region

Ectopic beats arise from micro-reentries near infarct regions in simulations of a patient-specific heart model

Ectopic beats are known to be involved in the initiation of a variety of cardiac arrhythmias. Although their location may vary, ectopic excitations have been found to originate from infarct areas, regions of micro-fibrosis and other heterogeneous tissues. However, the underlying mechanisms that link ectopic foci to heterogeneous tissues have yet to be fully understood. In this work, we investigate the mechanism of micro-reentry that leads to the generation of ectopic beats near infarct areas using a patient-specific heart model. The patient-specific geometrical model of the heart, including scar and peri-infarct zones, is obtained through magnetic resonance imaging (MRI). The infarct region is composed of ischemic myocytes and non-conducting cells (fibrosis, for instance). Electrophysiology is captured using an established cardiac myocyte model of the human ventricle modified to describe ischemia. The simulation results clearly reveal that ectopic beats emerge from micro-reentries that are sustained by the heterogeneous structure of the infarct regions. Because microscopic information about the heterogeneous structure of the infarct regions is not available, Monte-Carlo simulations are used to identify the probabilities of an infarct region to behave as an ectopic focus for different levels of ischemia and different percentages of non-conducting cells. From the proposed model, it is observed that ectopic beats are generated when a percentage of non-conducting cells is near a topological metric known as the percolation threshold. Although the mechanism for micro-reentries was proposed half a century ago to be a source of ectopic beats or premature ventricular contractions during myocardial infarction, the present study is the first to reproduce this mechanism in-silico using patient-specific data.Peer ReviewedPostprint (published version

A rule-based method for predicting the electrical activation of the heart with cardiac resynchronization therapy from non-invasive clinical data

Item does not contain fulltex

Improved Hybrid/GPU Algorithm for Solving Cardiac Electrophysiology Problems on Purkinje Networks

Cardiac Purkinje fibres provide an important pathway to the coordinated contraction of the heart. We present a numerical algorithm for the solution of electrophysiology problems across the Purkinje network that is efficient enough to be used in in-silico studies on realistic Purkinje networks with physiologically detailed models of ion exchange at the cell membrane. The algorithm is based on operator splitting and is provided with three different implementations: pure CPU, hybrid CPU/GPU, and pure GPU. Compared to our previous work, we modify the explicit gap junction term at network bifurcations in order to improve its mathematical consistency. Due to this improved consistency of the model, we are able to perform an empirical convergence study against analytical solutions. The study verified that all three implementations produce equivalent convergence rates, which shows that the algorithm produces equivalent result across different hardware platforms. Finally, we compare the efficiency of all three implementations on Purkinje networks of increasing spatial resolution using membrane models of increasing complexity. Both hybrid and pure-GPU implementations outperform the pure-CPU implementation, but their relative performance difference depends on the size of the Purkinje network and the complexity of the membrane model used

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

Computational modelling of the human heart and multiscale simulation of its electrophysiological activity aimed at the treatment of cardiac arrhythmias related to ischaemia and Infarction

[ES] Las enfermedades cardiovasculares constituyen la principal causa de morbilidad y mortalidad a nivel mundial, causando en torno a 18 millones de muertes cada año. De entre ellas, la más común es la enfermedad isquémica cardíaca, habitualmente denominada como infarto de miocardio (IM). Tras superar un IM, un considerable número de pacientes desarrollan taquicardias ventriculares (TV) potencialmente mortales durante la fase crónica del IM, es decir, semanas, meses o incluso años después la fase aguda inicial. Este tipo concreto de TV normalmente se origina por una reentrada a través de canales de conducción (CC), filamentos de miocardio superviviente que atraviesan la cicatriz del infarto fibrosa y no conductora. Cuando los fármacos anti-arrítmicos resultan incapaces de evitar episodios recurrentes de TV, la ablación por radiofrecuencia (ARF), un procedimiento mínimamente invasivo realizado mediante cateterismo en el laboratorio de electrofisiología (EF), se usa habitualmente para interrumpir de manera permanente la propagación eléctrica a través de los CCs responsables de la TV. Sin embargo, además de ser invasivo, arriesgado y requerir mucho tiempo, en casos de TVs relacionadas con IM crónico, hasta un 50% de los pacientes continúa padeciendo episodios recurrentes de TV tras el procedimiento de ARF. Por tanto, existe la necesidad de desarrollar nuevas estrategias pre-procedimiento para mejorar la planificación de la ARF y, de ese modo, aumentar esta tasa de éxito relativamente baja. En primer lugar, realizamos una revisión exhaustiva de la literatura referente a los modelos cardiacos 3D existentes, con el fin de obtener un profundo conocimiento de sus principales características y los métodos usados en su construcción, con especial atención sobre los modelos orientados a simulación de EF cardíaca. Luego, usando datos clínicos de un paciente con historial de TV relacionada con infarto, diseñamos e implementamos una serie de estrategias y metodologías para (1) generar modelos computacionales 3D específicos de paciente de ventrículos infartados que puedan usarse para realizar simulaciones de EF cardíaca a nivel de órgano, incluyendo la cicatriz del infarto y la región circundante conocida como zona de borde (ZB); (2) construir modelos 3D de torso que permitan la obtención del ECG simulado; y (3) llevar a cabo estudios in-silico de EF personalizados y pre-procedimiento, tratando de replicar los verdaderos estudios de EF realizados en el laboratorio de EF antes de la ablación. La finalidad de estas metodologías es la de localizar los CCs en el modelo ventricular 3D para ayudar a definir los objetivos de ablación óptimos para el procedimiento de ARF. Por último, realizamos el estudio retrospectivo por simulación de un caso, en el que logramos inducir la TV reentrante relacionada con el infarto usando diferentes configuraciones de modelado para la ZB. Validamos nuestros resultados mediante la reproducción, con una precisión razonable, del ECG del paciente en TV, así como en ritmo sinusal a partir de los mapas de activación endocárdica obtenidos invasivamente mediante sistemas de mapeado electroanatómico en este último caso. Esto permitió encontrar la ubicación y analizar las características del CC responsable de la TV clínica. Cabe destacar que dicho estudio in-silico de EF podría haberse efectuado antes del procedimiento de ARF, puesto que nuestro planteamiento está completamente basado en datos clínicos no invasivos adquiridos antes de la intervención real. Estos resultados confirman la viabilidad de la realización de estudios in-silico de EF personalizados y pre-procedimiento de utilidad, así como el potencial del abordaje propuesto para llegar a ser en un futuro una herramienta de apoyo para la planificación de la ARF en casos de TVs reentrantes relacionadas con infarto. No obstante, la metodología propuesta requiere de notables mejoras y validación por medio de es[CA] Les malalties cardiovasculars constitueixen la principal causa de morbiditat i mortalitat a nivell mundial, causant entorn a 18 milions de morts cada any. De elles, la més comuna és la malaltia isquèmica cardíaca, habitualment denominada infart de miocardi (IM). Després de superar un IM, un considerable nombre de pacients desenvolupen taquicàrdies ventriculars (TV) potencialment mortals durant la fase crònica de l'IM, és a dir, setmanes, mesos i fins i tot anys després de la fase aguda inicial. Aquest tipus concret de TV normalment s'origina per una reentrada a través dels canals de conducció (CC), filaments de miocardi supervivent que travessen la cicatriu de l'infart fibrosa i no conductora. Quan els fàrmacs anti-arítmics resulten incapaços d'evitar episodis recurrents de TV, l'ablació per radiofreqüència (ARF), un procediment mínimament invasiu realitzat mitjançant cateterisme en el laboratori de electrofisiologia (EF), s'usa habitualment per a interrompre de manera permanent la propagació elèctrica a través dels CCs responsables de la TV. No obstant això, a més de ser invasiu, arriscat i requerir molt de temps, en casos de TVs relacionades amb IM crònic fins a un 50% dels pacients continua patint episodis recurrents de TV després del procediment d'ARF. Per tant, existeix la necessitat de desenvolupar noves estratègies pre-procediment per a millorar la planificació de l'ARF i, d'aquesta manera, augmentar la taxa d'èxit, que es relativament baixa. En primer lloc, realitzem una revisió exhaustiva de la literatura referent als models cardíacs 3D existents, amb la finalitat d'obtindre un profund coneixement de les seues principals característiques i els mètodes usats en la seua construcció, amb especial atenció sobre els models orientats a simulació de EF cardíaca. Posteriorment, usant dades clíniques d'un pacient amb historial de TV relacionada amb infart, dissenyem i implementem una sèrie d'estratègies i metodologies per a (1) generar models computacionals 3D específics de pacient de ventricles infartats capaços de realitzar simulacions de EF cardíaca a nivell d'òrgan, incloent la cicatriu de l'infart i la regió circumdant coneguda com a zona de vora (ZV); (2) construir models 3D de tors que permeten l'obtenció del ECG simulat; i (3) dur a terme estudis in-silico de EF personalitzats i pre-procediment, tractant de replicar els vertaders estudis de EF realitzats en el laboratori de EF abans de l'ablació. La finalitat d'aquestes metodologies és la de localitzar els CCs en el model ventricular 3D per a ajudar a definir els objectius d'ablació òptims per al procediment d'ARF. Finalment, a manera de prova de concepte, realitzem l'estudi retrospectiu per simulació d'un cas, en el qual aconseguim induir la TV reentrant relacionada amb l'infart usant diferents configuracions de modelatge per a la ZV. Validem els nostres resultats mitjançant la reproducció, amb una precisió raonable, del ECG del pacient en TV, així com en ritme sinusal a partir dels mapes d'activació endocardíac obtinguts invasivament mitjançant sistemes de mapatge electro-anatòmic en aquest últim cas. Això va permetre trobar la ubicació i analitzar les característiques del CC responsable de la TV clínica. Cal destacar que aquest estudi in-silico de EF podria haver-se efectuat abans del procediment d'ARF, ja que el nostre plantejament està completament basat en dades clíniques no invasius adquirits abans de la intervenció real. Aquests resultats confirmen la viabilitat de la realització d'estudis in-silico de EF personalitzats i pre-procediment d'utilitat, així com el potencial de l'abordatge proposat per a arribar a ser en un futur una eina de suport per a la planificació de l'ARF en casos de TVs reentrants relacionades amb infart. No obstant això, la metodologia proposada requereix de notables millores i validació per mitjà d'estudis de simulació amb grans cohorts de pacients.[EN] Cardiovascular diseases represent the main cause of morbidity and mortality worldwide, causing around 18 million deaths every year. Among these diseases, the most common one is the ischaemic heart disease, usually referred to as myocardial infarction (MI). After surviving to a MI, a considerable number of patients develop life-threatening ventricular tachycardias (VT) during the chronic stage of the MI, that is, weeks, months or even years after the initial acute phase. This particular type of VT is typically sustained by reentry through slow conducting channels (CC), which are filaments of surviving myocardium that cross the non-conducting fibrotic infarct scar. When anti-arrhythmic drugs are unable to prevent recurrent VT episodes, radiofrequency ablation (RFA), a minimally invasive procedure performed by catheterization in the electrophysiology (EP) laboratory, is commonly used to interrupt the electrical conduction through the CCs responsible for the VT permanently. However, besides being invasive, risky and time-consuming, in the cases of VTs related to chronic MI, up to 50% of patients continue suffering from recurrent VT episodes after the RFA procedure. Therefore, there exists a need to develop novel pre-procedural strategies to improve RFA planning and, thereby, increase this relatively low success rate. First, we conducted an exhaustive review of the literature associated with the existing 3D cardiac models in order to gain a deep knowledge about their main features and the methods used for their construction, with special focus on those models oriented to simulation of cardiac EP. Later, using a clinical dataset of a chronically infarcted patient with a history of infarct-related VT, we designed and implemented a number of strategies and methodologies to (1) build patient-specific 3D computational models of infarcted ventricles that can be used to perform simulations of cardiac EP at the organ level, including the infarct scar and the surrounding region known as border zone (BZ); (2) construct 3D torso models that enable to compute the simulated ECG; and (3) carry out pre-procedural personalized in-silico EP studies, trying to replicate the actual EP studies conducted in the EP laboratory prior to the ablation. The goal of these methodologies is to allow locating the CCs into the 3D ventricular model in order to help in defining the optimal ablation targets for the RFA procedure. Lastly, as a proof-of-concept, we performed a retrospective simulation case study, in which we were able to induce an infarct-related reentrant VT using different modelling configurations for the BZ. We validated our results by reproducing with a reasonable accuracy the patient's ECG during VT, as well as in sinus rhythm from the endocardial activation maps invasively recorded via electroanatomical mapping systems in this latter case. This allowed us to find the location and analyse the features of the CC responsible for the clinical VT. Importantly, such in-silico EP study might have been conducted prior to the RFA procedure, since our approach is completely based on non-invasive clinical data acquired before the real intervention. These results confirm the feasibility of performing useful pre-procedural personalized in-silico EP studies, as well as the potential of the proposed approach to become a helpful tool for RFA planning in cases of infarct-related reentrant VTs in the future. Nevertheless, the developed methodology requires further improvements and validation by means of simulation studies including large cohorts of patients.During the carrying out of this doctoral thesis, the author Alejandro Daniel López Pérez was financially supported by the Ministerio de Economía, Industria y Competitividad of Spain through the program Ayudas para contratos predoctorales para la formación de doctores, with the grant number BES-2013-064089.López Pérez, AD. (2019). Computational modelling of the human heart and multiscale simulation of its electrophysiological activity aimed at the treatment of cardiac arrhythmias related to ischaemia and Infarction [Tesis doctoral]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/124973TESI

lifex-ep: a robust and efficient software for cardiac electrophysiology simulations

Background: Simulating the cardiac function requires the numerical solution of multi-physics and multi-scale mathematical models. This underscores the need for streamlined, accurate, and high-performance computational tools. Despite the dedicated endeavors of various research teams, comprehensive and user-friendly software programs for cardiac simulations, capable of accurately replicating both normal and pathological conditions, are still in the process of achieving full maturity within the scientific community. Results: This work introduces lifex-ep, a publicly available software for numerical simulations of the electrophysiology activity of the cardiac muscle, under both normal and pathological conditions. lifex-ep employs the monodomain equation to model the heart's electrical activity. It incorporates both phenomenological and second-generation ionic models. These models are discretized using the Finite Element method on tetrahedral or hexahedral meshes. Additionally, lifex-ep integrates the generation of myocardial fibers based on Laplace-Dirichlet Rule-Based Methods, previously released in Africa et al., 2023, within lifex-fiber. As an alternative, users can also choose to import myofibers from a file. This paper provides a concise overview of the mathematical models and numerical methods underlying lifex-ep, along with comprehensive implementation details and instructions for users. lifex-ep features exceptional parallel speedup, scaling efficiently when using up to thousands of cores, and its implementation has been verified against an established benchmark problem for computational electrophysiology. We showcase the key features of lifex-ep through various idealized and realistic simulations conducted in both normal and pathological scenarios. Furthermore, the software offers a user-friendly and flexible interface, simplifying the setup of simulations using self-documenting parameter files. Conclusions: lifex-ep provides easy access to cardiac electrophysiology simulations for a wide user community. It offers a computational tool that integrates models and accurate methods for simulating cardiac electrophysiology within a high-performance framework, while maintaining a user-friendly interface. lifex-ep represents a valuable tool for conducting in silico patient-specific simulations

Contribution to the improvement of electrical therapies and to the comprehension of electrophysiological mechanisms in heart failure and acute ischemia using computational simulation

[ES] Una mejor comprensión de los mecanismos subyacentes a las arritmias ventriculares, así como una mejora de las terapias eléctricas y farmacológicas asociadas, son un factor clave para prevenir la muerte súbita cardíaca en pacientes con cardiopatías estructurales y eléctricas. Una miocardiopatía importante que puede provocar arritmias ventriculares potencialmente mortales es la insuficiencia cardíaca (HF). Los pacientes con HF a menudo sufren también de bloqueo de rama izquierda (LBBB) que deteriora su condición. Actualmente, el tratamiento más eficaz para estos pacientes es la terapia de resincronización cardíaca (CRT). Sin embargo, no se alcanza una respuesta positiva en todos los casos, por lo que es necesario un mayor estudio para mejorar este tratamiento. Una segunda patología cardíaca que también produce arritmias letales es la isquemia miocárdica. Evidencia experimental ha demostrado que las alteraciones electrofisiológicas en el miocardio ventricular constituyen un sustrato para la generación de arritmias durante la fase aguda de isquemia. Estas alteraciones son inducidas por los tres componentes isquémicos principales: hipercalemia, hipoxia y acidosis. Sin embargo, la influencia de cada componente en los mecanismos de inicio y mantenimiento de las arritmias no se comprende aún con claridad. Una primera parte de esta tesis doctoral, se centra en la optimización de la CRT durante su aplicación en un corazón que padece HF y LBBB. Para esto, se modificó el modelo de potencial de acción (AP) de O'Hara para simular una velocidad de conducción realista tanto en condiciones sanas como patológicas. Además, se estimó e incorporó un sistema de His-Purkinje (HPS) dentro de un modelo biventricular/torso humano 3D para simular un LBBB realista. A continuación, se desarrolló un conjunto de simulaciones computacionales para diferentes configuraciones de la CRT a fin de determinar la posición y el instante de estimulación óptimo que conducen a la duración más corta del QRS. Posteriormente, los resultados se compararon con otros criterios de optimización. Los principales hallazgos de este estudio mostraron la necesidad de definir criterios de optimización mejores o complementarios, como un índice basado en el tiempo hasta alcanzar el 90% del área del QRS sugerido en este trabajo, para alcanzar la mejor sincronía eléctrica ventricular durante la aplicación de la CRT. Además, nuestros resultados también muestran que el septo superior cercano al tracto de salida es un sitio alternativo para la estimulación del ventrículo derecho, lo cual evita los problemas de perforación de la pared apical durante el procedimiento típico de la CRT. Por último, para obtener mejores resultados de la CRT se deben considerar protocolos de estimulación endocárdica en el ventrículo izquierdo. En la segunda parte de esta tesis se investigó los efectos de los tres componentes principales de la isquemia sobre la vulnerabilidad a una reentrada, así como el papel del HPS y sus mecanismos de acción en la generación y mantenimiento de arritmias ventriculares. Para lograr este objetivo, en primer lugar, se modificó el modelo AP ventricular para simular de forma realista las principales alteraciones provocadas por la isquemia miocárdica aguda. Las simulaciones se realizaron en un modelo biventricular humano 3D, acoplado en un torso virtual, que incluye una geometría realista de las zonas isquémicas central y de borde, así como un HPS detallado. Se simularon cuatro escenarios de severidad isquémica correspondientes a diferentes minutos de oclusión de la arteria coronaria para evaluar los efectos de la evolución de la isquemia en el tiempo. Luego, se evaluó la influencia individual de la hipercalemia, hipoxia y acidosis en el ancho de la ventana vulnerable (VW) a reentradas durante siete escenarios de isquemia aguda. Finalmente, se repitió este último conjunto de simulaciones isquémicas utilizando el modelo anatómico sin el HPS para evaluar el efecto de este último en la VW. Los resultados muestran que una condición isquémica moderada es el peor escenario para la generación de una reentrada. La hipoxia es el componente isquémico con el efecto más significativo en el ancho de la VW. Además, el flujo de corriente retrógrado desde el miocardio hacia el HPS en la región isquémica, los bloqueos de conducción en secciones discretas del HPS y el grado de hiperkalemia que afecta a las células de Purkinje, son sugeridos como mecanismos que podrían favorecer la aparición de arritmias ventriculares.[EN] A better understanding of the mechanisms underlying ventricular arrhythmias, as well as an improvement of the associated electrical and pharmacological therapies, are a key factor to prevent sudden cardiac death in patients with structural and electrical heart diseases. An important cardiomyopathy that can lead to life-threatening ventricular arrhythmias is heart failure (HF). Patients with HF also often suffer from left bundle branch block (LBBB), which worsens their condition. Currently, the most effective treatment to these patients is cardiac resynchronization therapy (CRT). However, many patients are non-responders, so further studies are needed to improve this treatment. A second cardiac pathology that also produces lethal arrhythmias is myocardial ischemia. Substantial experimental evidence has shown that electrophysiological alterations in the ventricular myocardium constitute a substrate for the generation of arrhythmias during the acute phase of ischemia. These alterations are induced by the three main ischemic components: hyperkalemia, hypoxia and acidosis. However, the influence of each component in the mechanisms of arrhythmia initiation and maintenance is still not completely understood. In the first section of this doctoral thesis, we focus on the optimization of CRT during its application in a heart suffering from HF and LBBB. For this purpose, we modified the O'Hara action potential (AP) model to simulate a realistic conduction velocity both in healthy and pathological conditions. In addition, a His-Purkinje system (HPS) was generated and incorporated into a 3D human biventricular/torso model to simulate realistic LBBB. A set of computational simulations were performed for different CRT configurations to determine the optimal pacing leads location and delay values leading to the shortest QRS duration. Subsequently, results were compared with other optimization criteria. The main findings of this study showed the need of better or complementary optimization criteria, such as an index based on the time to reach the 90% of the QRS area suggested in this work, to reach the best ventricular electrical synchrony during the CRT application. In addition, our results also show that the upper septum close to the outflow tract is an alternative site for the right ventricle (RV) stimulation, which avoids the perforation problems of the RV apical wall during the typical CRT procedure. Finally, protocols of left ventricle endocardial pacing should be considered to obtain better CRT results. In the second section of this thesis, we investigated the effects of the three main components of ischemia on the vulnerability to reentry, as well as the role of the HPS and its mechanisms of action in the generation and maintenance of ventricular arrhythmias. In order to achieve our goal, we first modified the ventricular AP model to realistically simulate the major alterations caused by acute myocardial ischemia. Simulations were performed in a 3D human biventricular model, embedded in a virtual torso, which includes a realistic geometry of the central and border ischemic zones, as well as a detailed HPS. Four scenarios of ischemic severity corresponding to different minutes after coronary artery occlusion were simulated to evaluate the effects of the evolution of ischemia over time. Then, the individual influence of hyperkalemia, hypoxia and acidosis in the width of the vulnerable window (VW) for reentry was assessed during seven scenarios of acute ischemia. Finally, this last set of ischemic simulations was repeated using the anatomical model without the HPS to evaluate the effect of the latter in the VW. Results show that a moderate ischemic condition is the worst scenario for reentry generation. Hypoxia is the ischemic component with the most significant effect on the width of the VW. Furthermore, the retrograde current flow from the myocardium to the HPS in the ischemic region, conduction blocks in discrete sections of the HPS, and the degree of hyperkalemia affecting the Purkinje cells, are suggested as HPS mechanisms that could favor the triggering of ventricular arrhythmias.[CA] Una millor comprensió dels mecanismes subjacents a les arrítmies ventriculars, així com una millora de les teràpies elèctriques i farmacològiques associades, són un factor clau per a previndre la mort sobtada cardíaca en pacients amb cardiopaties estructurals i elèctriques. Una miocardiopatia important que pot provocar arrítmies ventriculars potencialment mortals és la insuficiència cardíaca (HF). Els pacients amb HF sovint pateixen també de bloqueig de branca esquerra (LBBB) que deteriora la seua condició. Actualment, el tractament més eficaç per a aquests pacients és la teràpia de resincronització cardíaca (CRT). No obstant això, no s'aconsegueix una resposta positiva en tots els casos, per la qual cosa és necessari un major estudi per a millorar aquest tractament. Una segona patologia cardíaca que també produeix arrítmies letals és la isquèmia miocàrdica. Evidència experimental ha demostrat que les alteracions electrofisiològiques en el miocardi ventricular constitueixen un substrat per a la generació d'arrítmies durant la fase aguda d'isquèmia. Aquestes alteracions són induïdes pels tres components isquèmics principals: hipercalèmia, hipòxia i acidosi. No obstant això, la influència de cada component en els mecanismes d'inici i manteniment de les arrítmies no es comprén encara amb claredat. Una primera part d'aquesta tesi doctoral, se centra en l'optimització de la CRT durant la seua aplicació en un cor que pateix HF i LBBB. Per a això, es va modificar el model de potencial d'acció (AP) de O'Hara per a simular una velocitat de conducció realista tant en condicions sanes com patològiques. A més, es va estimar i es va incorporar un sistema de His-Purkinje (HPS) dins d'un model biventricular/tors humà 3D per a simular un LBBB realista. A continuació, es va desenvolupar un conjunt de simulacions computacionals per a diferents configuracions de la CRT a fi de determinar la posició i l'instant d'estimulació òptim que condueixen a la duració més curta del QRS. Posteriorment, els resultats es van comparar amb altres criteris d'optimització. Les principals troballes d'aquest estudi van mostrar la necessitat de definir millors o complementaris criteris d'optimització, com un índex basat en el temps fins a aconseguir el 90% de l'àrea del QRS suggerida en aquest treball, per a aconseguir la millor sincronia elèctrica ventricular durant l'aplicació de la CRT. A més, els nostres resultats també mostren que el septe superior pròxim al tracte d'eixida és un lloc alternatiu per a l'estimulació del ventricle dret, la cual cosa evita els problemes de perforació de la paret apical durant el procediment típic de la CRT. Finalment, per a obtindre millors resultats de la CRT s'han de considerar protocols d'estimulació endocárdica en el ventricle esquerre. En la segona part d'aquesta tesi es va investigar els efectes dels tres components principals de la isquèmia sobre la vulnerabilitat a una reentrada, així com el paper del HPS i els seus mecanismes d'acció en la generació i manteniment d'arrítmies ventriculars. Per a aconseguir aquest objectiu, en primer lloc es va modificar el model AP ventricular per a simular de manera realista les principals alteracions provocades per la isquèmia miocàrdica aguda. Les simulacions es van realitzar en un model biventricular humà 3D, acoblat en un tors virtual, que inclou una geometria realista de les zones isquèmiques central i de vora, així com un HPS detallat. Es van simular quatre escenaris de severitat isquèmica corresponents a diferents minuts d'oclusió de l'artèria coronària per a avaluar els efectes de l'evolució de la isquèmia en el temps. Després, es va avaluar la influència individual de la hipercalèmia, hipòxia i acidosi en l'ample de la finestra vulnerable (VW) a reentradas durant set escenaris d'isquèmia aguda. Finalment, es va repetir aquest últim conjunt de simulacions isquèmiques utilitzant el model anatòmic sense el HPS per a avaluar l'efecte d'aquest últim en la VW. Els resultats mostren que una condició isquèmica moderada és el pitjor escenari per a la generació d'una reentrada. La hipòxia és el component isquèmic amb l'efecte més significatiu en l'ample de la VW. A més, el flux de corrent retrògrad des del miocardi cap al HPS a la regió isquèmica, els bloquejos de conducció en seccions discretes del HPS i el grau d'hiperkalèmia que afecta les cèl·lules de Purkinje, són suggerits com a mecanismes que podrien afavorir l'aparició d'arrítmies ventriculars.Carpio Garay, EF. (2021). Contribution to the improvement of electrical therapies and to the comprehension of electrophysiological mechanisms in heart failure and acute ischemia using computational simulation [Tesis doctoral]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/163041TESI

Personalized Electromechanical Modeling of the Human Heart : Challenges and Opportunities for the Simulation of Pathophysiological Scenarios

Mathematische Modelle des menschlichen Herzens entwickeln sich zu einem Eckpfeiler der personalisierten Medizin. Sie sind ein nützliches Instrument und helfen klinischen Entscheidungsträgern die zugrundeliegenden Mechanismen von Herzkrankheiten zu erforschen und zu verstehen. Aufgrund der Komplexität des Herzens benötigen derartige Modelle allerdings eine detaillierte Beschreibung der physikalischen Prozesse, welche auf verschiedenen räumlichen und zeitlichen Skalen miteinander interagieren. Aus mathematischer Perspektive stellen vor allem die Entwicklung robuster numerischer Methoden für die Lösung des Modells in Raum und Zeit sowie die Identifizierung von Parametern aus patientenspezifischen Messungen eine Herausforderung dar. In dieser Arbeit wird ein detailliertes mathematisches Modell vorgestellt, welches ein vollgekoppeltes Multiskalenmodell des menschlichen Herzens beschreibt. Das Modell beinhaltet unter anderem die Ausbreitung des elektrischen Signals und die mechanische Verformung des Herzmuskels sowie eine Beschreibung des Herz-Kreislauf-Systems. Basierend auf dem neusten Stand der Technik wurden Modelle der Membrankinetik sowie der Entwicklung der aktiven Kraft zu einem einheitlichen Modell einer Herzmuskelzelle zusammengeführt. Dieses beschreibt die elektromechanische Kopplung in Herzmuskelzellen der Vorhöfe und der Herzkammern basierend auf der Physiologie im Menschen und wurde mit Hilfe von experimentellen Daten aus einzelnen Zellen neu parametrisiert. Um das elektromechanisch gekoppelte Modell des menschlichen Herzens lösen zu können, wurde ein gestaffeltes Lösungsverfahren entwickelt, welches auf bereits existierenden Softwarelösungen der Elektrophysiologie und Mechanik aufbaut. Das neue Modell wurde verwendet, um den Einfluss elektromechanischer Rückkopplungseffekte auf das Herz im Sinusrhythmus zu untersuchen. Die Simulationsergebnisse zeigten, dass elektromechanische Rückkopplungseffekte auf zellulärer Ebene einen wesentlichen Einfluss auf das mechanische Verhalten des Herzens haben. Dahingegen hatte die Verformung des Herzens nur einen geringen Einfluss auf den Diffusionskoeffizienten des elektrischen Signals. Um die verschiedenen Komponenten der Simulationssoftware zu verifizieren, wurden spezielle Probleme definiert, welche die wichtigsten Aspekte der Elektrophysiologie und der Mechanik abdecken. Zusätzlich wurden diese Probleme dazu verwendet, den Einfluss von räumlicher und zeitlicher Diskretisierung auf die numerische Lösung zu bewerten. Die Ergebnisse zeigten, dass Raum- und Zeitdiskretisierung vor allem für das elektrophysiologische Problem die limitierenden Faktoren sind, während die Mechanik hauptsächlich anfällig für volumenversteifende Effekte ist. Weiterhin wurde das Modell verwendet, um zu untersuchen, wie sich eine Verteilung der Faserspannung auf den gesamten Herzmuskel auf die Funktion der linken Herzkammer auswirkt. Hierzu wurde zusätzlich eine Spannung in die Normalenrichtungen der Fasern einer idealisierten linken Herzkammer angewandt. Es zeigte sich, dass insbesondere eine Spannung senkrecht zu den Faserschichten zu einer physiologischeren Kontraktion der Kammer führte. Allerdings konnten diese Ergebnisse auf einem ganzen Herzen nicht vollständig bestätigt werden. In einem zweiten Projekt wurde mit Hilfe eines Modells der linken Herzkammer untersucht, wie sich das Rotationsmuster der Kammer unter Modifikation der lokalen elektromechanischen Eigenschaften verändert. Hierzu wurden in vivo Daten elektromechanischer Parameter von 30 Patienten mit Herzversagen und Linksschenkelblock in das Modell integriert, simuliert und ausgewertet. Die Ergebnisse konnten die klinisch aufgestellte Hypothese nicht bestätigen und es zeigte sich keine Korrelation zwischen den elektromechanischen Parametern und dem Rotationsverhalten. Die Auswirkungen von standardisierten Ablationsstrategien zur Behandlung von Vorhofflimmern in Bezug auf die kardiovaskuläre Leistung wurde in einem Modell des ganzen Herzens untersucht. Aufgrund der Narben im linken Vorhof wurde die elektrische Aktivierung und die Steifigkeit des Herzmuskels verändert. Dies führte zu einem reduzierten Auswurfvolumen, welches in direktem Zusammenhang mit dem inaktiven Gewebe steht. Abhängig von der Steifigkeit der Narben hat sich zusätzlich der Druck im linken Vorhof erhöht. Die linke Herzkammer war nur wenig beeinflusst. Zu guter Letzt wurden schrittweise pathologische Mechanismen in das Herzmodell integriert, welche in Zusammenhang mit Herzversagen stehen und in Patienten mit dilatativer Kardiomyopathie zu beobachten sind. Die Simulationen zeigten, dass vor allem zelluläre Veränderungen bezüglich der elektrophysiologischen Eigenschaften für die schlechte mechanische Aktivtät des Herzens verantwortlich sind. Weiterhin zeigte sich, dass strukturelle Veränderungen der Anatomie und die erhöhte Steifigkeit des Herzmuskels und die damit einhergehenden Anpassungen des Herz-Kreislauf-Systems nötig sind, um in vivo Messungen zu reproduzieren. In dieser Arbeit wurde eine Simulationsumgebung vorgestellt, welche die Berechnung der elektromechanischen Aktivität des Herzens und des Herz-Kreislauf-Systems ermöglicht. Die Simulationsumgebung wurde mit Hilfe von einfachen Beispielen verifiziert und unter Einbeziehung von Daten aus der Magnetresonanztomographie validiert. Zu guter Letzt wurde die Simulationsumgebung genutzt, um klinische Fragen zu beantworten, welche andernfalls im Dunkeln blieben