A multiscale subvoxel perfusion model to estimate diffusive capillary wall conductivity in multiple sclerosis lesions from perfusion MRI data

We propose a new mathematical model to learn capillary leakage coefficients from dynamic susceptibility contrast MRI data. To this end, we derive an embedded mixed-dimension flow and transport model for brain tissue perfusion on a sub-voxel scale. This model is used to obtain the contrast agent concentration distribution in a single MRI voxel during a perfusion MRI sequence. We further present a magnetic resonance signal model for the considered sequence including a model for local susceptibility effects. This allows modeling MR signal-time curves that can be compared to clinical MRI data. The proposed model can be used as a forward model in the inverse modeling problem of inferring model parameters such as the diffusive capillary wall conductivity. Acute multiple sclerosis lesions are associated with a breach in the integrity of the blood brain barrier. Applying the model to perfusion MR data of a patient with acute multiple sclerosis lesions, we conclude that diffusive capillary wall conductivity is a good indicator for characterizing activity of lesions, even if other patient-specific model parameters are not well-known. This article is protected by copyright. All rights reserved.

[1]  Steven P Sourbron,et al.  Classic models for dynamic contrast‐enhanced MRI , 2013, NMR in biomedicine.

[2]  J. Slotboom,et al.  Characterization of Microcirculation in Multiple Sclerosis Lesions by Dynamic Texture Parameter Analysis (DTPA) , 2013, PloS one.

[3]  A. Katchalsky,et al.  Thermodynamic analysis of the permeability of biological membranes to non-electrolytes. , 1958, Biochimica et biophysica acta.

[4]  Miroslav Kuchta,et al.  Sub-voxel Perfusion Modeling in Terms of Coupled 3d-1d Problem , 2017, Lecture Notes in Computational Science and Engineering.

[5]  Rainer Helmig,et al.  A New Simulation Framework for Soil–Root Interaction, Evaporation, Root Growth, and Solute Transport , 2018, Vadose Zone Journal.

[6]  Paul W. Cleary,et al.  An Efficient Computational Approach to Characterize DSC-MRI Signals Arising from Three-Dimensional Heterogeneous Tissue Structures , 2014, PloS one.

[7]  Characterization of Enhancing MS Lesions by Dynamic Texture Parameter Analysis of Dynamic Susceptibility Perfusion Imaging , 2016, BioMed research international.

[8]  A. Quarteroni,et al.  On the coupling of 1D and 3D diffusion-reaction equations. Applications to tissue perfusion problems , 2008 .

[9]  C. Crone,et al.  THE PERMEABILITY OF CAPILLARIES IN VARIOUS ORGANS AS DETERMINED BY USE OF THE 'INDICATOR DIFFUSION' METHOD. , 1963, Acta physiologica Scandinavica.

[10]  Timothy W. Secomb,et al.  Green's Function Methods for Analysis of Oxygen Delivery to Tissue by Microvascular Networks , 2004, Annals of Biomedical Engineering.

[11]  E. Rostrup,et al.  Abnormal blood–brain barrier permeability in normal appearing white matter in multiple sclerosis investigated by MRI☆☆☆ , 2013, NeuroImage: Clinical.

[12]  C. Nicholson,et al.  Diffusion in brain extracellular space. , 2008, Physiological reviews.

[13]  A. Quarteroni.  Cardiovascular mathematics , 2000 .

[14]  Daniel Foreman-Mackey,et al.  emcee: The MCMC Hammer , 2012, 1202.3665.

[15]  Frauke Zipp,et al.  Changes in cerebral perfusion precede plaque formation in multiple sclerosis: a longitudinal perfusion MRI study. , 2004, Brain : a journal of neurology.

[16]  Alireza Minagar,et al.  Blood-brain barrier disruption in multiple sclerosis , 2003, Multiple sclerosis.

[17]  Andreas Dedner,et al.  A generic grid interface for parallel and adaptive scientific computing. Part I: abstract framework , 2008, Computing.

[18]  Rebecca J Shipley,et al.  Multiscale Modeling of Fluid Transport in Tumors , 2008, Bulletin of mathematical biology.

[19]  Anna Devor,et al.  Oxygen advection and diffusion in a three- dimensional vascular anatomical network. , 2008, Optics express.

[20]  T. Secomb,et al.  Theoretical Simulation of Oxygen Transport to Brain by Networks of Microvessels: Effects of Oxygen Supply and Demand on Tissue Hypoxia , 2000, Microcirculation.

[21]  Giustina Casagrande,et al.  A computational model for microcirculation including Fahraeus‐Lindqvist effect, plasma skimming and fluid exchange with the tissue interstitium , 2018, International journal for numerical methods in biomedical engineering.

[22]  Jeffrey A. Cohen,et al.  Diagnostic criteria for multiple sclerosis: 2010 Revisions to the McDonald criteria , 2011, Annals of neurology.

[23]  A. Powell,et al.  Characterisation and application of ultra-high spin clusters as magnetic resonance relaxation agents. , 2015, Dalton transactions.

[24]  V G Kiselev.  On the theoretical basis of perfusion measurements by dynamic susceptibility contrast MRI , 2001, Magnetic resonance in medicine.

[25]  Roland Wiest,et al.  Differentiating enhancing multiple sclerosis lesions, glioblastoma, and lymphoma with dynamic texture parameters analysis (DTPA): A feasibility study , 2017, Medical physics.

[26]  Daniel Foreman-Mackey,et al.  corner.py: Scatterplot matrices in Python , 2016, J. Open Source Softw..

[27]  Masako Sugihara-Seki,et al.  Blood flow and permeability in microvessels , 2005 .

[28]  D L Buckley,et al.  Tracer kinetic modelling in MRI: estimating perfusion and capillary permeability , 2012, Physics in medicine and biology.

[29]  F. Curry,et al.  Microvascular permeability. , 1999, Physiological reviews.

[30]  Carlo D'Angelo,et al.  Finite Element Approximation of Elliptic Problems with Dirac Measure Terms in Weighted Spaces: Applications to One- and Three-dimensional Coupled Problems , 2012, SIAM J. Numer. Anal..

[31]  Steffen Müthing,et al.  DuMux: DUNE for multi-{phase,component,scale,physics,…} flow and transport in porous media , 2011 .

[32]  Rebecca J. Shipley,et al.  Multiscale Modelling of Fluid and Drug Transport in Vascular Tumours , 2010, Bulletin of mathematical biology.

[33]  Paolo Decuzzi,et al.  On the near-wall accumulation of injectable particles in the microcirculation: smaller is not better , 2013, Scientific Reports.

[34]  E. Larsson,et al.  Calculation of cerebral perfusion parameters using regional arterial input functions identified by factor analysis , 2006, Journal of magnetic resonance imaging : JMRI.

[35]  A. Pries,et al.  Biophysical aspects of blood flow in the microvasculature. , 1996, Cardiovascular research.

[36]  Rakesh K. Jain,et al.  Transport of molecules across tumor vasculature , 2004, Cancer and Metastasis Reviews.

[37]  Mads Nybo,et al.  Measurement of plasma viscosity by free oscillation rheometry: imprecision, sample stability and establishment of a new reference range , 2014, Annals of clinical biochemistry.

[38]  James M. Provenzale,et al.  Principles of T2*‐weighted dynamic susceptibility contrast MRI technique in brain tumor imaging , 2015, Journal of magnetic resonance imaging : JMRI.

[39]  E. Motti,et al.  The terminal vascular bed in the superficial cortex of the rat. An SEM study of corrosion casts. , 1986, Journal of neurosurgery.

[40]  Barbara I. Wohlmuth,et al.  A Local Error Estimate for the Poisson Equation with a Line Source Term , 2015, ENUMATH.

[41]  P. Tofts,et al.  Measurement of the blood‐brain barrier permeability and leakage space using dynamic MR imaging. 1. Fundamental concepts , 1991, Magnetic resonance in medicine.

[42]  T E Yankeelov,et al.  A theoretical framework to model DSC-MRI data acquired in the presence of contrast agent extravasation , 2009, Physics in medicine and biology.

[43]  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.

[44]  A. Traboulsee,et al.  What Have We Learned from Perfusion MRI in Multiple Sclerosis? , 2018, American Journal of Neuroradiology.

[45]  B E Kendall,et al.  Breakdown of the blood-brain barrier precedes symptoms and other MRI signs of new lesions in multiple sclerosis. Pathogenetic and clinical implications. , 1990, Brain : a journal of neurology.

[46]  Jonathan R Goodman,et al.  Ensemble samplers with affine invariance , 2010 .

[47]  B. Murray,et al.  Sucrose diffusion in aqueous solution , 2016, Physical chemistry chemical physics : PCCP.

[48]  A. Quarteroni,et al.  One-dimensional models for blood flow in arteries , 2003 .

[49]  A. Linninger,et al.  Cerebral Microcirculation and Oxygen Tension in the Human Secondary Cortex , 2013, Annals of Biomedical Engineering.

[50]  Barbara I. Wohlmuth,et al.  Hybrid Models for Simulating Blood Flow in Microvascular Networks , 2018, Multiscale Model. Simul..

[51]  Joanna M. Wardlaw,et al.  Tracer kinetic modelling for DCE-MRI quantification of subtle blood–brain barrier permeability , 2016, NeuroImage.

[52]  S. McQuaid,et al.  Abnormal Endothelial Tight Junctions in Active Lesions and Normal‐appearing White Matter in Multiple Sclerosis , 2002, Brain pathology.

[53]  B. Haraldsson,et al.  Transport of macromolecules across microvascular walls: the two-pore theory. , 1994, Physiological reviews.

[54]  R M Weisskoff,et al.  Relative cerebral blood volume maps corrected for contrast agent extravasation significantly correlate with glioma tumor grade, whereas uncorrected maps do not. , 2006, AJNR. American journal of neuroradiology.

[55]  Daniel S. Reich,et al.  Direct MRI detection of impending plaque development in multiple sclerosis , 2015, Neurology: Neuroimmunology & Neuroinflammation.

[56]  V. Kiselev,et al.  Transverse NMR relaxation as a probe of mesoscopic structure. , 2002, Physical review letters.

[57]  V. Huxley,et al.  Permeability of single capillaries to intermediate-sized colored solutes. , 1983, The American journal of physiology.

[58]  Andreas Dedner,et al.  A generic grid interface for parallel and adaptive scientific computing. Part II: implementation and tests in DUNE , 2008, Computing.

[59]  Glyn Johnson,et al.  Dynamic susceptibility contrast perfusion MR imaging of multiple sclerosis lesions: characterizing hemodynamic impairment and inflammatory activity. , 2005, AJNR. American journal of neuroradiology.

[60]  C. Heesen,et al.  T1 Recovery Is Predominantly Found in Black Holes and Is Associated with Clinical Improvement in Patients with Multiple Sclerosis , 2017, American Journal of Neuroradiology.

[61]  B. Douglas Ward,et al.  A novel technique for modeling susceptibility-based contrast mechanisms for arbitrary microvascular geometries: The finite perturber method , 2008, NeuroImage.

[62]  B. Rosen,et al.  High resolution measurement of cerebral blood flow using intravascular tracer bolus passages. Part I: Mathematical approach and statistical analysis , 1996, Magnetic resonance in medicine.

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

[64]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[65]  P Zunino,et al.  Modelling mass and heat transfer in nano-based cancer hyperthermia , 2015, Royal Society Open Science.

[66]  R K Jain,et al.  Transport of fluid and macromolecules in tumors. I. Role of interstitial pressure and convection. , 1989, Microvascular research.

[67]  F. Calamante.  Perfusion MRI Using Dynamic-Susceptibility Contrast MRI: Quantification Issues in Patient Studies , 2010, Topics in magnetic resonance imaging : TMRI.

[68]  Barbara Wohlmuth,et al.  Mathematical modelling, analysis and numerical approximation of second order elliptic problems with inclusions , 2018 .

[69]  Axel R Pries,et al.  A hybrid discrete-continuum approach for modelling microcirculatory blood flow. , 2019, Mathematical medicine and biology : a journal of the IMA.

[70]  P. Villoslada,et al.  T2 hypointense rims and ring-enhancing lesions in MS , 2010, Multiple sclerosis.

[71]  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.

[72]  Mikael Mortensen,et al.  Preconditioning trace coupled 3d‐1d systems using fractional Laplacian , 2016, Numerical Methods for Partial Differential Equations.

[73]  F. Curry.  Permeability coefficients of the capillary wall to low molecular weight hydrophilic solutes measured in single perfused capillaries of frog mesentery. , 1979, Microvascular research.

[74]  J. Squire,et al.  Similar endothelial glycocalyx structures in microvessels from a range of mammalian tissues: evidence for a common filtering mechanism? , 2011, Biophysical journal.

[75]  V. Kiselev.  Transverse relaxation effect of MRI contrast agents: A crucial issue for quantitative measurements of cerebral perfusion , 2005, Journal of magnetic resonance imaging : JMRI.

[76]  Maximilian Reiser,et al.  Quantification of Perfusion and Permeability in Multiple Sclerosis: Dynamic Contrast-Enhanced MRI in 3D at 3T , 2012, Investigative radiology.

[77]  Levick.  Capillary filtration‐absorption balance reconsidered in light of dynamic extravascular factors , 1991, Experimental physiology.

[78]  J. Mintorovitch,et al.  Comparison of Magnetic Properties of MRI Contrast Media Solutions at Different Magnetic Field Strengths , 2005, Investigative radiology.

[79]  Laura Cattaneo,et al.  Computational models for fluid exchange between microcirculation and tissue interstitium , 2014, Networks Heterog. Media.