A porous circulation model of the human brain for in silico clinical trials in ischaemic stroke

The advancement of ischaemic stroke treatment relies on resource-intensive experiments and clinical trials. In order to improve ischaemic stroke treatments, such as thrombolysis and thrombectomy, we target the development of computational tools for in silico trials which can partially replace these animal and human experiments with fast simulations. This study proposes a model that will serve as part of a predictive unit within an in silico clinical trial estimating patient outcome as a function of treatment. In particular, the present work aims at the development and evaluation of an organ-scale microcirculation model of the human brain for perfusion prediction. The model relies on a three-compartment porous continuum approach. Firstly, a fast and robust method is established to compute the anisotropic permeability tensors representing arterioles and venules. Secondly, vessel encoded arterial spin labelling magnetic resonance imaging and clustering are employed to create an anatomically accurate mapping between the microcirculation and large arteries by identifying superficial perfusion territories. Thirdly, the parameter space of the problem is reduced by analysing the governing equations and experimental data. Fourthly, a parameter optimization is conducted. Finally, simulations are performed with the tuned model to obtain perfusion maps corresponding to an open and an occluded (ischaemic stroke) scenario. The perfusion map in the occluded vessel scenario shows promising qualitative agreement with computed tomography images of a patient with ischaemic stroke caused by large vessel occlusion. The results highlight that in the case of vessel occlusion (i) identifying perfusion territories is essential to capture the location and extent of underperfused regions and (ii) anisotropic permeability tensors are required to give quantitatively realistic estimation of perfusion change. In the future, the model will be thoroughly validated against experiments.

[1]  Lutz Jäncke,et al.  Brain size and grey matter volume in the healthy human brain , 2002, Neuroreport.

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

[3]  P. Vespa,et al.  Analysis of Thrombi Retrieved from Cerebral Arteries of Patients with Acute Ischemic Stroke. , 2005 .

[4]  P. Jezzard,et al.  Cerebral blood flow quantification using vessel-encoded arterial spin labeling , 2013, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[5]  Myriam Peyrounette,et al.  Multiscale modelling of blood flow in cerebral microcirculation: Details at capillary scale control accuracy at the level of the cortex , 2018, PloS one.

[6]  T. Józsa,et al.  Boundary conditions for flow simulations of abdominal aortic aneurysms , 2014 .

[7]  P. Tsui,et al.  Investigating cerebral oedema using poroelasticity. , 2016, Medical engineering & physics.

[8]  E. Arsava,et al.  Can Restoring Incomplete Microcirculatory Reperfusion Improve Stroke Outcome after Thrombolysis? , 2012, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[9]  H. Duvernoy,et al.  Cortical blood vessels of the human brain , 1981, Brain Research Bulletin.

[10]  J. Desilles,et al.  Modified Thrombolysis in Cerebral Infarction 2C/Thrombolysis in Cerebral Infarction 3 Reperfusion Should Be the Aim of Mechanical Thrombectomy: Insights From the ASTER Trial (Contact Aspiration Versus Stent Retriever for Successful Revascularization) , 2018, Stroke.

[11]  Liesbet Geris,et al.  UvA-DARE (Digital Academic Repository) Virtual physiological human 2016: translating the virtual physiological human to the clinic , 2022 .

[12]  Stephen J. Payne,et al.  Investigating the effects of a penetrating vessel occlusion with a multi-scale microvasculature model of the human cerebral cortex , 2018, NeuroImage.

[13]  S. Payne Cerebral Blood Flow and Metabolism: A Quantitative Approach , 2017 .

[14]  C. Iadecola The Neurovascular Unit Coming of Age: A Journey through Neurovascular Coupling in Health and Disease , 2017, Neuron.

[15]  Myriam Peyrounette,et al.  Brain Capillary Networks Across Species: A few Simple Organizational Requirements Are Sufficient to Reproduce Both Structure and Function , 2019, Front. Physiol..

[16]  Stephen Payne,et al.  The Influence of Network Structure on the Transport of Blood in the Human Cerebral Microvasculature , 2012, Microcirculation.

[17]  Stephen J Payne,et al.  A statistical model of the penetrating arterioles and venules in the human cerebral cortex , 2016, Microcirculation.

[18]  D. Liebeskind Reperfusion for acute ischemic stroke: arterial revascularization and collateral therapeutics , 2010, Current opinion in neurology.

[19]  Lucas O. Müller,et al.  Blood pressure gradients in cerebral arteries: a clue to pathogenesis of cerebral small vessel disease , 2017, Stroke and Vascular Neurology.

[20]  A. Pries,et al.  Resistance to blood flow in microvessels in vivo. , 1994, Circulation research.

[21]  J. Bogousslavsky,et al.  Arterial territories of the human brain , 1998, Neurology.

[22]  Henk A. van der Vorst,et al.  Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems , 1992, SIAM J. Sci. Comput..

[23]  L. Beenen,et al.  Automated Cerebral Infarct Volume Measurement in Follow-up Noncontrast CT Scans of Patients with Acute Ischemic Stroke , 2013, American Journal of Neuroradiology.

[24]  N. Voets,et al.  Cognition based bTBI mechanistic criteria; a tool for preventive and therapeutic innovations , 2018, Scientific Reports.

[25]  Hang Si,et al.  TetGen, a Delaunay-Based Quality Tetrahedral Mesh Generator , 2015, ACM Trans. Math. Softw..

[26]  P. Tofts,et al.  Normal cerebral perfusion measurements using arterial spin labeling: Reproducibility, stability, and age and gender effects , 2004, Magnetic resonance in medicine.

[27]  David J Gladstone,et al.  Toward Wisdom From Failure: Lessons From Neuroprotective Stroke Trials and New Therapeutic Directions , 2002, Stroke.

[28]  M. Krause,et al.  Endovascular therapy for ischemic stroke with perfusion-imaging selection. , 2015, The New England journal of medicine.

[29]  M. Goldacre,et al.  Determinants of the decline in mortality from acute stroke in England: linked national database study of 795 869 adults , 2019, BMJ.

[30]  Mariano Vázquez,et al.  Fully coupled fluid‐electro‐mechanical model of the human heart for supercomputers , 2018, International journal for numerical methods in biomedical engineering.

[31]  David Attwell,et al.  What is a pericyte? , 2016, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[32]  John A. Nelder,et al.  A Simplex Method for Function Minimization , 1965, Comput. J..

[33]  Alfons G. Hoekstra,et al.  Coupling one-dimensional arterial blood flow to three-dimensional tissue perfusion models for in silico trials of acute ischaemic stroke , 2020, Interface Focus.

[34]  S. Sherwin,et al.  Pulse wave propagation in a model human arterial network: assessment of 1-D numerical simulations against in vitro measurements. , 2007, Journal of biomechanics.

[35]  T. Sochi,et al.  A computationally efficient framework for the simulation of cardiac perfusion using a multi‐compartment Darcy porous‐media flow model , 2013, International journal for numerical methods in biomedical engineering.

[36]  Jorge Nocedal,et al.  Algorithm 778: L-BFGS-B: Fortran subroutines for large-scale bound-constrained optimization , 1997, TOMS.

[37]  S. Sherwin,et al.  Modelling the circle of Willis to assess the effects of anatomical variations and occlusions on cerebral flows. , 2007, Journal of biomechanics.

[38]  Hester F. Lingsma,et al.  MR CLEAN, a multicenter randomized clinical trial of endovascular treatment for acute ischemic stroke in the Netherlands: study protocol for a randomized controlled trial , 2014, Trials.

[39]  J A Reggia,et al.  Pathogenic mechanisms in ischemic damage: a computational study , 1999, Comput. Biol. Medicine.

[40]  Anders Logg,et al.  DOLFIN: Automated finite element computing , 2010, TOMS.

[41]  D. F. R. Heijtel,et al.  Cerebral Perfusion Measurements in Elderly with Hypertension Using Arterial Spin Labeling , 2015, PloS one.

[42]  Jaime Grutzendler,et al.  Regional Blood Flow in the Normal and Ischemic Brain Is Controlled by Arteriolar Smooth Muscle Cell Contractility and Not by Capillary Pericytes , 2015, Neuron.

[43]  Ayush Goyal,et al.  Multi-Scale Parameterisation of a Myocardial Perfusion Model Using Whole-Organ Arterial Networks , 2013, Annals of Biomedical Engineering.

[44]  Patrick Jenny,et al.  Depth-dependent flow and pressure characteristics in cortical microvascular networks , 2017, PLoS Comput. Biol..

[45]  A. Demchuk,et al.  Endovascular thrombectomy after large-vessel ischaemic stroke: a meta-analysis of individual patient data from five randomised trials , 2016, The Lancet.

[46]  Eric S Donkor,et al.  Stroke in the 21st Century: A Snapshot of the Burden, Epidemiology, and Quality of Life , 2018, Stroke research and treatment.

[47]  E. Mohammadi,et al.  Barriers and facilitators related to the implementation of a physiological track and trigger system: A systematic review of the qualitative evidence , 2017, International journal for quality in health care : journal of the International Society for Quality in Health Care.

[48]  Jack Lee,et al.  Parameterisation of multi-scale continuum perfusion models from discrete vascular networks , 2012, Medical & Biological Engineering & Computing.

[49]  P. Vespa,et al.  Analysis of Thrombi Retrieved From Cerebral Arteries of Patients With Acute Ischemic Stroke , 2006, Stroke.

[50]  Erlend Hodneland,et al.  A new framework for assessing subject-specific whole brain circulation and perfusion using MRI-based measurements and a multi-scale continuous flow model , 2019, PLoS Comput. Biol..

[51]  Olav Jansen,et al.  Diffusion-weighted imaging or computerized tomography perfusion assessment with clinical mismatch in the triage of wake up and late presenting strokes undergoing neurointervention with Trevo (DAWN) trial methods , 2017, International journal of stroke : official journal of the International Stroke Society.

[52]  M. Chappell,et al.  Relationship between haemodynamic impairment and collateral blood flow in carotid artery disease , 2017, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[53]  Anders Logg,et al.  The FEniCS Project Version 1.5 , 2015 .

[54]  D. Heistad,et al.  Effects of Chronic Hypertension and Sympathetic Nerves on the Cerebral Microvasculature of Stroke‐Prone Spontaneously Hypertensive Rats , 1984, Circulation research.

[55]  E. Ruppin,et al.  Modeling Brain Energy Metabolism and Function: A Multiparametric Monitoring Approach , 2006, Bulletin of mathematical biology.

[56]  Laurent Risser,et al.  A 3D-investigation shows that angiogenesis in primate cerebral cortex mainly occurs at capillary level , 2009, International Journal of Developmental Neuroscience.

[57]  V. L. Doyle,et al.  Middle Cerebral Artery Territory Infarction Spectroscopy and Initial Infarct Volume by MRI Predicts Outcome in Patients With Acetyl Aspartate Concentration by Magnetic Resonance N-Measurement of Initial , 1999 .

[58]  Hervé Delingette,et al.  Coupled personalization of cardiac electrophysiology models for prediction of ischaemic ventricular tachycardia , 2011, Interface Focus.

[59]  Alain Goriely,et al.  A physics-based model explains the prion-like features of neurodegeneration in Alzheimer’s disease, Parkinson’s disease, and amyotrophic lateral sclerosis , 2019, Journal of the Mechanics and Physics of Solids.

[60]  H. Bohlen,et al.  Microvascular adaptation in the cerebral cortex of adult spontaneously hypertensive rats. , 1984, Hypertension.

[61]  Rong Zhang,et al.  Distribution of cardiac output to the brain across the adult lifespan , 2017, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[62]  R. Gonzalez,et al.  In patients with suspected acute stroke, CT perfusion-based cerebral blood flow maps cannot substitute for DWI in measuring the ischemic core , 2017, PloS one.

[63]  Gergely Bárdossy,et al.  A “backward” calculation method for the estimation of central aortic pressure wave in a 1D arterial model network , 2013 .

[64]  Gary R. Mirams,et al.  mRNA Expression Levels in Failing Human Hearts Predict Cellular Electrophysiological Remodeling: A Population-Based Simulation Study , 2013, PloS one.

[65]  S. Sotiropoulos,et al.  On the mechanical behaviour of PEEK and HA cranial implants under impact loading. , 2017, Journal of the mechanical behavior of biomedical materials.

[66]  Guillaume Houzeaux,et al.  Large-scale CFD simulations of the transitional and turbulent regime for the large human airways during rapid inhalation , 2016, Comput. Biol. Medicine.

[67]  Alejandro F. Frangi,et al.  Subject-specific multi-poroelastic model for exploring the risk factors associated with the early stages of Alzheimer's disease , 2017, Interface Focus.

[68]  Guillaume Houzeaux,et al.  A massively parallel computational electrophysiology model of the heart , 2011 .

[69]  A. Hoekstra,et al.  In-Silico Trials for Treatment of Acute Ischemic Stroke , 2020, Frontiers in Neurology.

[70]  Jeroen Hendrikse,et al.  Mapping of cerebral perfusion territories using territorial arterial spin labeling: techniques and clinical application , 2013, NMR in biomedicine.

[71]  Hester F. Lingsma,et al.  A randomized trial of intraarterial treatment for acute ischemic stroke. , 2015, The New England journal of medicine.

[72]  Gert Kwakkel,et al.  Early Prediction of Outcome of Activities of Daily Living After Stroke: A Systematic Review , 2011, Stroke.

[73]  Hester F. Lingsma,et al.  Follow-up infarct volume as a mediator of endovascular treatment effect on functional outcome in ischaemic stroke , 2018, European Radiology.

[74]  I. Marshall,et al.  How does blood regulate cerebral temperatures during hypothermia? , 2018, Scientific Reports.

[75]  S. Payne,et al.  Multi-scale homogenization of blood flow in 3-dimensional human cerebral microvascular networks. , 2015, Journal of theoretical biology.

[76]  Céline Fouard,et al.  A Novel Three‐Dimensional Computer‐Assisted Method for a Quantitative Study of Microvascular Networks of the Human Cerebral Cortex , 2006, Microcirculation.

[77]  F. Bajd,et al.  Mathematical modeling of blood clot fragmentation during flow-mediated thrombolysis. , 2013, Biophysical journal.

[78]  K. S. Burrowes,et al.  Multi-scale computational models of the airways to unravel the pathophysiological mechanisms in asthma and chronic obstructive pulmonary disease (AirPROM) , 2013, Interface Focus.

[79]  Aart J. Nederveen,et al.  Accuracy and precision of pseudo-continuous arterial spin labeling perfusion during baseline and hypercapnia: A head-to-head comparison with 15O H2O positron emission tomography , 2014, NeuroImage.