HPC for Simulation and Modelling of Cardiac Infarction for the Development of Novel and Effective Treatment Thomas Franz Chris Barnard Division of Cardiothoracic Surgery, Programme for the Enhancement of Research Capacity - Research Office, Centre for Research in Computational and Applied Mechanics University of Cape Town; Centre for High Performance Computing www.freedigitalphotos.net
Acknowledgements Dr Neil Davies Prof Daya Reddy Dr Sebastian Skatulla Dr Greg Mitchell Dr Jun Liao Dr Bruce Spottiswoode Prof Ernesta Meintjes PhD Students Karen Kadner* Mazin Sirry Fulufhelo Masithulela Renee Miller Ritesh Rama Post-docs Dr Jeroen Kortsmit* Dr Dieter Legner Laura Dubuis MSc Students Renee Miller* James Mbewu* Ritesh Rama* Muhammad Saleh* Peter Wise Kevin Sack Mohammed Essack Jules Lancee (*completed) CHPC Kevin Colville Eric Mbele BSc Students Jesse Macadangdang* Gesant Abed* Devin Dollery* Tshenolo Tau* Dale Lewis* Aniel Prosper* Siddhartha Ram* Bruce Phelp* Kelebogile Motili* Mogoma Senthumule Wayde Herman Rimon Muchehiwa 1
WHO: World Health Statistics 2007. Geneva, 2007 2
Myocardial Infarct (MI) 3
Changes in the Heart Healthy Infarcted 4
Therapies for MI Acorn CSD Aim: Prevention of adverse remodelling of the heart Cell delivery therapy Injection of cells with carrier medium into infarcted region of the heart Positive outcomes Mechanisms unclear 5
Infarct Thickness & Cardiac Function 1 Month after Infarction * p<0.05 n=7-9 6
www.freedigitalphotos.net Questions on Injection Therapy Contribution to mechanical and biological properties of the heart Properties of injectable biomaterials Timing of delivery Injection site Injection volume 7
Canine Biventricular Cardiac Model Myofibre representation Constitutive models Passive strain energy function Active contraction model Coupling to Windkessel circulation model (Nielsen et al. 1991) Antero-apical infarct: 22% LV wall Infarct wall thinning Post-infarct LV dilation 100% (Kortsmit, Davies, Miller, Franz)
Hydrogel Injectates Gel inclusion in infarct region Volume: 9.4 ml (7.1% of LV wall) Transmural distribution 8 thin layers Single bulk layer Kadner et al, Biomaterials, 2012
Time after coronary occlusion Study Design and Ventricular Function Case C Functional Parameters V 0 (ml) SV (ml) EF (%) Healthy Control H 0.88 13.9 15.2 35.6 Ischemic Infarct II 0.88 17.8 6.61 15.5 + layered gel II + L 0.88 16.7 7.02 17.7 + bulk gel II + B 0.88 19.9 4.78 12.4 Scarred Dilated Infarct SDI 8.80 39.4 9.28 11.8 + layered gel SDI + L 8.80 27.9 13.0 19.0 + bulk gel SDI + B 8.80 27.2 13.2 19.2
Stress and Strain in Infarct Region End-diastolic Myofibre Strain End-diastolic Myofibre Stress End-systolic Myofibre Strain End-systolic Myofibre Stress
Main Findings Functional Improvement More pronounced in scarred dilated infarct (still dilated LV cavity volume volume overload) Layered injectate superior at early infarct stage Bulk injectate superior at late infarct stage Infarct Stress and Strain Layered injectates better than bulk for all infarct stages Greatest benefits Stresses at end-diastole (filling) Strains at end-systole (ejection)
Meshless parallelised code (SESKA) Geometry: Rat LV, truncated ellipsoid Vol.-fraction : 100.0% / 0.0% / 0.0% 84.2% / 15.8% / 0.0% 71.0% / 15.8% / 13.2% (Healthy / Infarct / B-Zone) Material: Healthy: Omens et al. (2002) Infarct: Stiffness x 0.5 Injection: Neo-Hookean B.C.: Endocardium fixed at the base Loading: Diastolic pressure (Legner, Skatulla, Mbewu, Rama, Sack, Mothibi, Reddy)
Fibre Strain Infarct + Injectate Infarct Healthy
Constitutive Properties of Infarcts (Sirry, Liao, Davies, Franz) 15
Stress (Pa) Infarct Stress-Strain Relationships 40 000 35 000 circumferential 28d longitudinal 28d 30 000 25 000 14d 7d 14d 7d 20 000 0d 0d 15 000 10 000 5 000 0-0.02 0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 Strain 16
Inverse FE Modelling 17
Cardiac Regeneration Mechanical Stimuli in Cardiac Environment Cell Differentiation Cell Proliferation Gene Expression
Cell-specific Geometrical Modelling Confocal microscopy images Slicing through sample 3D Cell geometry (Lancee, Davies, Dubuis, Franz)
Cell-specific Finite Element Model Microscope images Cytoskeleton 3D Reconstruction Cytoplasm Nucleus Cell membrane Stretching FEM Nucleus Focal adhesion Elastic substrate (Dubuis, Davies, Green, Franz)
Cell Deformation Bottom View 3D cell-specific FEM Cytoplasm, nucleus, membrane Focal adhesion Physiological stretch loading Top View
FSI Model Simulating blood flow around a cell: Fluid structure interaction
Tensile Strain in Fibroblast 3D view of nucleus 3D view on cell 3D Cross-sectional view
In Progress: Endothelial Cell Model