A personalized computational model of edema formation in myocarditis based on long-axis biventricular MRI images

Myocarditis is defined as the inflammation of the myocardium, i.e. the cardiac muscle. Among the reasons that lead to this disease, we may include infections caused by a virus, bacteria, protozoa, fungus, and others. One of the signs of the inflammation is the formation of edema, which may be a consequence of the interaction between interstitial fluid dynamics and immune response. This complex physiological process was mathematically modeled using a nonlinear system of partial differential equations (PDE) based on porous media approach. By combing a model based on Biot’s poroelasticity theory with a model for the immune response we developed a new hydro-mechanical model for inflammatory edema. To verify this new computational model, T2 parametric mapping obtained by Magnetic Resonance (MR) imaging was used to identify the region of edema in a patient diagnosed with unspecific myocarditis. A patient-specific geometrical model was created using MRI images from the patient with myocarditis. With this model, edema formation was simulated using the proposed hydro-mechanical mathematical model in a two-dimensional domain. The computer simulations allowed us to correlate spatiotemporal dynamics of representative cells of the immune systems, such as leucocytes and the pathogen, with fluid accumulation and cardiac tissue deformation. This study demonstrates that the proposed mathematical model is a very promising tool to better understand edema formation in myocarditis. Simulations obtained from a patient-specific model reproduced important aspects related to the formation of cardiac edema, its area, position, and shape, and how these features are related to immune response.

[1]  Alan S. Perelson,et al.  Agent-based modeling of host–pathogen systems: The successes and challenges , 2008, Information Sciences.

[2]  Rodrigo Weber dos Santos,et al.  On the use of Gillespie stochastic simulation algorithm in a model of the human immune system response to the Yellow Fever vaccine , 2017, 2017 IEEE International Conference on Bioinformatics and Biomedicine (BIBM).

[3]  Manuel E. Cornely,et al.  Lipedema and Lymphatic Edema , 2006 .

[4]  Ruy Freitas Reis,et al.  Influence of the Immune System on the Biological Dynamics of the Interstitial Fluid Pressure , 2017 .

[5]  Colin Phipps,et al.  Mathematical Model of the Effect of Interstitial Fluid Pressure on Angiogenic Behavior in Solid Tumors , 2011, Comput. Math. Methods Medicine.

[6]  Marcelo Lobosco,et al.  On the Coupling of Two Models of the Human Immune Response to an Antigen , 2014, BioMed research international.

[7]  Justin Guinney,et al.  GSVA: gene set variation analysis for microarray and RNA-Seq data , 2013, BMC Bioinformatics.

[8]  Anders Logg,et al.  Automated Solution of Differential Equations by the Finite Element Method: The FEniCS Book , 2012 .

[9]  Rodrigo Weber dos Santos,et al.  On the computational modeling of the innate immune system , 2013, BMC Bioinformatics.

[10]  Helen M Byrne,et al.  Wound healing angiogenesis: the clinical implications of a simple mathematical model. , 2012, Journal of theoretical biology.

[11]  R. W. dos Santos,et al.  Ectopic beats arise from micro-reentries near infarct regions in simulations of a patient-specific heart model , 2018, Scientific Reports.

[12]  P E Seiden,et al.  A computer model of cellular interactions in the immune system. , 1992, Immunology today.

[13]  P. Lancellotti,et al.  [Update on myocarditis]. , 2018, Revue medicale de Liege.

[14]  Rodrigo Weber dos Santos,et al.  Computational Modeling of Microabscess Formation , 2012, Comput. Math. Methods Medicine.

[15]  IEEE International Conference on Bioinformatics and Biomedicine, BIBM 2018, Madrid, Spain, December 3-6, 2018 , 2018, BIBM.

[16]  Vicente Grau,et al.  A poroelastic model coupled to a fluid network with applications in lung modelling , 2014, International journal for numerical methods in biomedical engineering.

[17]  Alan S. Perelson,et al.  Stochastic Theory of Early Viral Infection: Continuous versus Burst Production of Virions , 2011, PLoS Comput. Biol..

[18]  Eike Nagel,et al.  T1 and T2 mapping in myocarditis: seeing beyond the horizon of Lake Louise criteria and histopathology , 2018, Expert review of cardiovascular therapy.

[19]  Rodrigo Weber dos Santos,et al.  An Hydro-Mechanical Model of Edema Formation Applied to Bacterial Myocarditis , 2018, 2018 IEEE International Conference on Bioinformatics and Biomedicine (BIBM).

[20]  Timoteo Carletti,et al.  The Stochastic Evolution of a Protocell: The Gillespie Algorithm in a Dynamically Varying Volume , 2011, Comput. Math. Methods Medicine.

[21]  Katharine Best,et al.  Mathematical modeling of within‐host Zika virus dynamics , 2018, Immunological reviews.

[22]  M. Biot General Theory of Three‐Dimensional Consolidation , 1941 .

[23]  Ingo Eitel,et al.  T2-weighted cardiovascular magnetic resonance in acute cardiac disease , 2011, Journal of cardiovascular magnetic resonance : official journal of the Society for Cardiovascular Magnetic Resonance.

[24]  P. Zunino,et al.  A computational model of drug delivery through microcirculation to compare different tumor treatments , 2014, International journal for numerical methods in biomedical engineering.

[25]  Virginia H. Huxley,et al.  Capillary Fluid Exchange: Regulation, Functions, and Pathology , 2010 .

[26]  G Thiene,et al.  Postmortem diagnosis in sudden cardiac death victims: macroscopic, microscopic and molecular findings , 2001, Cardiovascular research.

[27]  Renato Natal Jorge,et al.  Modelling skin wound healing angiogenesis: A review. , 2018, Journal of theoretical biology.

[28]  Amber Denomme,et al.  Pulmonary edema. , 1988, Journal of thoracic imaging.

[29]  Anders Logg,et al.  Solving PDEs in Python: The FEniCS Tutorial I , 2017 .

[30]  Lucia Russo,et al.  Mathematical modeling of infectious disease dynamics , 2013, Virulence.

[31]  Hilde van der Togt,et al.  Publisher's Note , 2003, J. Netw. Comput. Appl..

[32]  Malte Kelm,et al.  Abnormal T2 mapping cardiovascular magnetic resonance correlates with adverse clinical outcome in patients with suspected acute myocarditis , 2017, Journal of Cardiovascular Magnetic Resonance.

[33]  A.P.S. Selvadurai,et al.  Coupled hydro-mechanical effects in a poro-hyperelastic material , 2016 .

[34]  S. Goodreau,et al.  Mathematical Modeling of Infectious Disease , 2015 .

[35]  B Schnackenburg,et al.  T2 mapping cardiovascular magnetic resonance identifies the presence of myocardial inflammation in patients with dilated cardiomyopathy as compared to endomyocardial biopsy , 2018, European heart journal cardiovascular Imaging.

[36]  Richard B Thompson,et al.  Subclinical Pulmonary Edema Is Associated With Reduced Exercise Capacity in HFpEF and HFrEF. , 2017, Journal of the American College of Cardiology.

[37]  Denise E Kirschner,et al.  Multiple mechanisms allow Mycobacterium tuberculosis to continuously inhibit MHC class II-mediated antigen presentation by macrophages. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[38]  Triantafyllos Stylianopoulos,et al.  The role of mechanical forces in tumor growth and therapy. , 2014, Annual review of biomedical engineering.

[39]  E. Walshaw,et al.  A systematic approach , 2018, BDJ.

[40]  M. B E R N A S C H I 1 A N D F C A S T I G L I O N E Selection of escape mutants from immune recognition during HIV infection , 2002 .

[41]  Ruy Freitas Reis,et al.  Interstitial Pressure Dynamics Due to Bacterial Infection , 2016 .

[42]  Matthias Gutberlet,et al.  Cardiovascular Magnetic Resonance in Myocarditis: A JACC White Paper , 2009 .

[43]  E. Starling On the Absorption of Fluids from the Connective Tissue Spaces , 1896, The Journal of physiology.

[44]  John Samuel,et al.  A simple cellular automaton model for influenza A viral infections. , 2004, Journal of theoretical biology.

[45]  Marcelo Lobosco,et al.  A qualitatively validated mathematical-computational model of the immune response to the yellow fever vaccine , 2018, BMC Immunology.

[46]  Jing Li,et al.  The Failure of R 0 , 2011, Computational and mathematical methods in medicine.

[47]  P. Brigidi,et al.  The Three Genetics (Nuclear DNA, Mitochondrial DNA, and Gut Microbiome) of Longevity in Humans Considered as Metaorganisms , 2014, BioMed research international.

[48]  S. Sell,et al.  How the immune system works. , 1980, Medical times.

[49]  B Maisch,et al.  Report of the 1995 World Health Organization/International Society and Federation of Cardiology Task Force on the Definition and Classification of cardiomyopathies. , 1996, Circulation.

[50]  E. McVeigh,et al.  T2‐prepared SSFP improves diagnostic confidence in edema imaging in acute myocardial infarction compared to turbo spin echo , 2007, Magnetic resonance in medicine.

[51]  Marcelo Lobosco,et al.  Development of a Computational Model of Abscess Formation , 2018, Front. Microbiol..

[52]  Tiina Heliö,et al.  Current state of knowledge on aetiology, diagnosis, management, and therapy of myocarditis: a position statement of the European Society of Cardiology Working Group on Myocardial and Pericardial Diseases. , 2013, European heart journal.

[53]  R Puzone,et al.  A systematic approach to vaccine complexity using an automaton model of the cellular and humoral immune system. I. Viral characteristics and polarized responses. , 2000, Vaccine.

[54]  A.P.S. Selvadurai,et al.  On poro-hyperelastic shear , 2016 .

[55]  Ruy Freitas Reis,et al.  A Plasma Flow Model in the Interstitial Tissue Due to Bacterial Infection , 2016, IWBBIO.

[56]  H. Lee,et al.  Myocardial T1 and T2 Mapping: Techniques and Clinical Applications , 2017, Korean journal of radiology.

[57]  A. Huerta,et al.  Finite Element Methods for Flow Problems , 2003 .

[58]  Rodrigo Weber dos Santos,et al.  On the mathematical modeling of inflammatory edema formation , 2019, Comput. Math. Appl..

[59]  Christophe Geuzaine,et al.  Gmsh: A 3‐D finite element mesh generator with built‐in pre‐ and post‐processing facilities , 2009 .

[60]  L. Dubuquoy,et al.  Increased lymphatic vessel density and lymphangiogenesis in inflammatory bowel disease , 2011, Alimentary pharmacology & therapeutics.

[61]  B. Sua,et al.  Mathematical modelling of immune response in tissues , 2009 .

[62]  Francesco Migliavacca,et al.  An interactive simulation tool for patient-specific clinical decision support in single-ventricle physiology. , 2018, The Journal of thoracic and cardiovascular surgery.

[63]  Raymond J Kim,et al.  Late Gadolinium Enhancement Magnetic Resonance Imaging in the Diagnosis and Prognosis of Endomyocardial Fibrosis Patients , 2011, Circulation. Cardiovascular imaging.

[64]  A. Taghian,et al.  Hand Edema in Patients at Risk of Breast Cancer‐Related Lymphedema: Health Professionals Should Take Notice , 2018, Physical therapy.

[65]  Gilles Clermont,et al.  IN SILICO MODELS OF ACUTE INFLAMMATION IN ANIMALS , 2006, Shock.

[66]  S J Allen,et al.  Left ventricular myocardial edema. Lymph flow, interstitial fibrosis, and cardiac function. , 1991, Circulation research.