Infusion Mechanisms in Brain White Matter and Their Dependence on Microstructure: An Experimental Study of Hydraulic Permeability

Cell adhesion affects the properties of interstitial fluid flow: A study using multiscale poroelastic composite modeling

In this study, we conduct a multiscale, multiphysics modeling of the brain gray matter as a poroelastic composite. We develop a customized representative volume element based on cytoarchitectural features that encompass important microscopic components of the tissue, namely the extracellular space, the capillaries, the pericapillary space, the interstitial fluid, cell-cell and cell-capillary junctions, and neuronal and glial cell bodies. Using asymptotic homogenization and direct numerical simulation, the effective properties at the tissue level are identified based on microscopic properties. To analyze the influence of various microscopic elements on the effective/macroscopic properties and tissue response, we perform sensitivity analyses on cell junction (cluster) stiffness, cell junction diameter (dimensions), and pericapillary space width. The results of this study suggest that changes in cell adhesion can greatly affect both mechanical and hydraulic (interstitial fluid flow and porosity) features of brain tissue, consistent with the effects of neurodegenerative diseases.

Porosity-permeability tensor relationship of closely and randomly packed fibrous biomaterials and biological tissues: Application to the brain white matter

The constitutive model for the porosity-permeability relationship is a powerful tool to estimate and design the transport properties of porous materials, which has attracted significant attention for the advancement of novel materials. However, in comparison with other materials, biomaterials, especially natural and artificial tissues, have more complex microstructures e.g. high anisotropy, high randomness of cell/fibre dimensions/position and very low porosity. Consequently, a reliable microstructure-permeability relationship of fibrous biomaterials has proven elusive. To fill this gap, we start a mathematical derivation from the fundamental brain white matter (WM) formed by nerve fibres. This is augmented by a numerical characterisation and experimental validations to obtain an anisotropic permeability tensor of the brain WM as a function of the tissue porosity. A versatile microstructure generation software (MicroFiM) for fibrous biomaterial with complex microstructure and low porosity was built accordingly and made freely accessible here. Moreover, we propose an anisotropic poro-hyperelastic model enhanced by the newly defined porosity-permeability tensor relationship which precisely captures the tissues macro-scale permeability changes due to the microstructural deformation in an infusion scenario. The constitutive model, theories and protocols established in this study will both provide improved design strategies to tailor the transport properties of fibrous biomaterials and enable the non-invasive characterisation of the transport properties of biological tissues. This will lead to the provision of better patient-specific medical treatments, such as drug delivery. STATEMENT OF SIGNIFICANCE: Due to the microstructural complexity, a reliable microstructure-permeability relationship of fibrous biomaterials has proven elusive, which hinders our way of tuning the fluid transport property of the biomaterials by directly programming their microstructure. The same problem hinders non-invasive characterisations of fluid transport properties in biological tissues, which can significantly improve the efficiency of treatments e.g. drug delivery, directly from the tissues accessible microstructural information, e.g. porosity. Here, we developed a validated mathematical formulation to link the random microstructure to a fibrous material's macroscale permeability tensor. This will advance our capability to design complex biomaterials and make it possible to non-invasively characterise the permeability of living tissues for precise treatment planning. The newly established theory and protocol can be easily adapted to various types of fibrous biomaterials.

Linking fluid-axons interactions to the macroscopic fluid transport properties of the brain

Many brain disorders, including Alzheimer's Disease and Parkinson's Disease, and drug delivery procedures are linked to fluid transport in the brain; yet, while neurons are extremely soft and can be easily deformed, how the microscale channel flow interacts with the neuronal structures (especially axons) deformation and how these interactions affect the macroscale tissue function and transport properties is poorly understood. Misrepresenting these relationships may lead to the erroneous prediction of e.g. disease spread, drug delivery, and nerve injury in the brain. However, understanding fluid-neuron interactions is an outstanding challenge because the behaviours of both phases are not only dynamic but also occur at an extremely small length scale (the width of the flow channel is ∼100 nm), which cannot be captured by state-of-the-art experimental techniques. Here, by explicitly simulating the dynamics of the flow and axons at the microstructural level, we, for the first time, establish the link between micromechanical tissue response to the physical laws governing the macroscopic transport property of the brain white matter. We found that interactions between axons and the interstitial flow are very strong, thus playing an essential role in the brain fluid/mass transport. Furthermore, we proposed the first anisotropic pressure-dependent permeability tensor informed by microstructural dynamics for more accurate brain modelling at the macroscale, and analysed the effect of the variation of the microstructural parameters that influence such tensor. These findings will shed light on some unsolved issues linked to brain functions and medical treatments relying on intracerebral transport, and the mathematical model provides a framework to more realistically model the brain and design brain-tissue-like biomaterials. STATEMENT OF SIGNIFICANCE: This study reveals how neurons interact with the fluid flowing around them and how these microscale interactions affect macroscale transport behaviour of the brain tissue. The findings provide unprecedented insights into some unsolved issues linked to brain functions and medical treatments relying on intracerebral fluid transport. Furthermore, we, for the first time, established a microstructure-informed permeability tensor as a function of local hydraulic pressure and pressure gradient for the brain tissue, which inherently captures the dynamic transport property of the brain. This study is a cornerstone to advance the predicting accuracy of brain tissue transport property and neural tissue engineering.

Poro-viscoelastic material parameter identification of brain tissue-mimicking hydrogels

Understanding and characterizing the mechanical and structural properties of brain tissue is essential for developing and calibrating reliable material models. Based on the Theory of Porous Media, a novel nonlinear poro-viscoelastic computational model was recently proposed to describe the mechanical response of the tissue under different loading conditions. The model contains parameters related to the time-dependent behavior arising from both the viscoelastic relaxation of the solid matrix and its interaction with the fluid phase. This study focuses on the characterization of these parameters through indentation experiments on a tailor-made polyvinyl alcohol-based hydrogel mimicking brain tissue. The material behavior is adjusted to ex vivo porcine brain tissue. An inverse parameter identification scheme using a trust region reflective algorithm is introduced and applied to match experimental data obtained from the indentation with the proposed computational model. By minimizing the error between experimental values and finite element simulation results, the optimal constitutive model parameters of the brain tissue-mimicking hydrogel are extracted. Finally, the model is validated using the derived material parameters in a finite element simulation.

Distinguishing poroelasticity and viscoelasticity of brain tissue with time scale

Brain tissue is considered to be biphasic, with approximately 80% liquid and 20% solid matrix, thus exhibiting viscoelasticity due to rearrangement of the solid matrix and poroelasticity due to fluid migration within the solid matrix. However, how to distinguish poroelastic and viscoelastic effects in brain tissue remains challenging. In this study, we proposed a method of unconfined compression-isometric hold to measure the force versus time relaxation curves of porcine brain tissue samples with systematically varied sample lengths. Upon scaling the measured relaxation force and relaxation time with different length-dependent physical quantities, we successfully distinguished the poroelasticity and viscoelasticity of the brain tissue. We demonstrated that during isometric hold, viscoelastic relaxation dominated the mechanical behavior of brain tissue in the short-time regime, while poroelastic relaxation dominated in the long-time regime. Furthermore, compared with poroelastic relaxation, viscoelastic relaxation was found to play a more dominant role in the mechanical response of porcine brain tissue. We then evaluated the differences between poroelastic and viscoelastic effects for both porcine and human brain tissue. Because of the draining of pore fluid, the Young's moduli in poroelastic relaxation were lower than those in viscoelastic relaxation; brain tissue changed from incompressible during viscoelastic relaxation to compressible during poroelastic relaxation, resulting in reduced Poisson ratios. This study provides new insights into the physical mechanisms underlying the roles of viscoelasticity and poroelasticity in brain tissue. STATEMENT OF SIGNIFICANCE: Although the poroviscoelastic model had been proposed to characterize brain tissue mechanical behavior, it is difficult to distinguish the poroelastic and viscoelastic behaviors of brain tissue. The study distinguished viscoelasticity and poroelasticity of brain tissue with time scales and then evaluated the differences between poroelastic and viscoelastic effects for both porcine and human brain tissue, which helps to accurate selection of constitutive models suitable for application in certain situations (e.g., pore-dominant and viscoelastic-dominant deformation).

Reconstruction of ovine axonal cytoarchitecture enables more accurate models of brain biomechanics

There is an increased need and focus to understand how local brain microstructure affects the transport of drug molecules directly administered to the brain tissue, for example in convection-enhanced delivery procedures. This study reports a systematic attempt to characterize the cytoarchitecture of commissural, long association and projection fibres, namely the corpus callosum, the fornix and the corona radiata, with the specific aim to map different regions of the tissue and provide essential information for the development of accurate models of brain biomechanics. Ovine samples are imaged using scanning electron microscopy combined with focused ion beam milling to generate 3D volume reconstructions of the tissue at subcellular spatial resolution. Focus is placed on the characteristic cytological feature of the white matter: the axons and their alignment in the tissue. For each tract, a 3D reconstruction of relatively large volumes, including a significant number of axons, is performed and outer axonal ellipticity, outer axonal cross-sectional area and their relative perimeter are measured. The study of well-resolved microstructural features provides useful insight into the fibrous organization of the tissue, whose micromechanical behaviour is that of a composite material presenting elliptical tortuous tubular axonal structures embedded in the extra-cellular matrix. Drug flow can be captured through microstructurally-based models using 3D volumes, either reconstructed directly from images or generated in silico using parameters extracted from the database of images, leading to a workflow to enable physically-accurate simulations of drug delivery to the targeted tissue.

Imaging and reconstruction of sheep brain axonal cytoarchitecture provides insight for brain biomechanics models that simulate drug delivery and other biological processes governed by interstitial fluid flow and transport.

On the microstructurally driven heterogeneous response of brain white matter to drug infusion pressure

Delivering therapeutic agents into the brain via convection-enhanced delivery (CED), a mechanically controlled infusion method, provides an efficient approach to bypass the blood–brain barrier and deliver drugs directly to the targeted focus in the brain. Mathematical methods based on Darcy’s law have been widely adopted to predict drug distribution in the brain to improve the accuracy and reduce the side effects of this technique. However, most of the current studies assume that the hydraulic permeability and porosity of brain tissue are homogeneous and constant during the infusion process, which is less accurate due to the deformability of the axonal structures and the extracellular matrix in brain white matter. To solve this problem, a multiscale model was established in this study, which takes into account the pressure-driven deformation of brain microstructure to quantify the change of local permeability and porosity. The simulation results were corroborated using experiments measuring hydraulic permeability in ovine brain samples. Results show that both hydraulic pressure and drug concentration in the brain would be significantly underestimated by classical Darcy’s law, thus highlighting the great importance of the present multiscale model in providing a better understanding of how drugs transport inside the brain and how brain tissue responds to the infusion pressure. This new method can assist the development of both new drugs for brain diseases and preoperative evaluation techniques for CED surgery, thus helping to improve the efficiency and precision of treatments for brain diseases.

Role of Tissue Hydraulic Permeability in Convection-Enhanced Delivery of Nanoparticle-Encapsulated Chemotherapy Drugs to Brain Tumour

PURPOSE

Tissue hydraulic permeability of brain tumours can vary considerably depending on the tissue microstructure, compositions in interstitium and tumour cells. Its effects on drug transport and accumulation remain poorly understood.

METHODS

Mathematical modelling is applied to predict the drug delivery outcomes in tumours with different tissue permeability upon convection-enhanced delivery. The modelling is based on a 3-D realistic tumour model that is extracted from patient magnetic resonance images.

RESULTS

Modelling results show that infusing drugs into a permeable tumour can facilitate a more favourable hydraulic environment for drug transport. The infused drugs will exhibit a relatively uniform distribution and cover a larger tumour volume for effective cell killing. Cross-comparisons show the delivery outcomes are more sensitive to the changes in tissue hydraulic permeability and blood pressure than the fluid flow from the brain ventricle. Quantitative analyses demonstrate that increasing the fluid gain from both the blood and brain ventricle can further improve the interstitial fluid flow, and thereby enhance the delivery outcomes. Furthermore, similar responses to the changes in tissue hydraulic permeability can be found for different types of drugs.

CONCLUSIONS

Tissue hydraulic permeability as an intrinsic property can influence drug accumulation and distribution. Results from this study can deepen the understanding of the interplays between drug and tissues that are involved in the drug delivery processes in chemotherapy.

Insights into Infusion-Based Targeted Drug Delivery in the Brain: Perspectives, Challenges and Opportunities

Targeted drug delivery in the brain is instrumental in the treatment of lethal brain diseases, such as glioblastoma multiforme, the most aggressive primary central nervous system tumour in adults. Infusion-based drug delivery techniques, which directly administer to the tissue for local treatment, as in convection-enhanced delivery (CED), provide an important opportunity; however, poor understanding of the pressure-driven drug transport mechanisms in the brain has hindered its ultimate success in clinical applications. In this review, we focus on the biomechanical and biochemical aspects of infusion-based targeted drug delivery in the brain and look into the underlying molecular level mechanisms. We discuss recent advances and challenges in the complementary field of medical robotics and its use in targeted drug delivery in the brain. A critical overview of current research in these areas and their clinical implications is provided. This review delivers new ideas and perspectives for further studies of targeted drug delivery in the brain.

Cortex tissue relaxation and slow to medium load rates dependency can be captured by a two-phase flow poroelastic model

This paper investigates the complex time-dependent behavior of cortex tissue, under adiabatic condition, using a two-phase flow poroelastic model. Motivated by experiments and Biot's consolidation theory, we tackle time-dependent uniaxial loading, confined and unconfined, with various geometries and loading rates from 1μm/s to 100μm/s. The cortex tissue is modeled as the porous solid saturated by two immiscible fluids, with dynamic viscosities separated by four orders, resulting in two different characteristic times. These are respectively associated to interstitial fluid and glial cells. The partial differential equations system is discretized in space by the finite element method and in time by Euler-implicit scheme. The solution is computed using a monolithic scheme within the open-source computational framework FEniCS. The parameters calibration is based on Sobol sensitivity analysis, which divides them into two groups: the tissue specific group, whose parameters represent general properties, and sample specific group, whose parameters have greater variations. Our results show that the experimental curves can be reproduced without the need to resort to viscous solid effects, by adding an additional fluid phase. Through this process, we aim to present multiphase poromechanics as a promising way to a unified brain tissue modeling framework in a variety of settings.

Microscale characterisation of the time-dependent mechanical behaviour of brain white matter

Brain mechanics is a topic of deep interest because of the significant role of mechanical cues in both brain function and form. Specifically, capturing the heterogeneous and anisotropic behaviour of cerebral white matter (WM) is extremely challenging and yet the data on WM at a spatial resolution relevant to tissue components are sparse. To investigate the time-dependent mechanical behaviour of WM, and its dependence on local microstructural features when subjected to small deformations, we conducted atomic force microscopy (AFM) stress relaxation experiments on corpus callosum (CC), corona radiata (CR) and fornix (FO) of fresh ovine brain. Our experimental results show a dependency of the tissue mechanical response on axons orientation, with e.g. the stiffness of perpendicular and parallel samples is different in all three regions of WM whereas the relaxation behaviour is different for the CC and FO regions. An inverse modelling approach was adopted to extract Prony series parameters of the tissue components, i.e. axons and extra cellular matrix with its accessory cells, from experimental data. Using a bottom-up approach, we developed analytical and FEA estimates that are in good agreement with our experimental results. Our systematic characterisation of sheep brain WM using a combination of AFM experiments and micromechanical models provide a significant contribution for predicting localised time-dependent mechanics of brain tissue. This information can lead to more accurate computational simulations, therefore aiding the development of surgical robotic solutions for drug delivery and accurate tissue mimics, as well as the determination of criteria for tissue injury and predict brain development and disease progression.

On the microstructural origin of brain white matter hydraulic permeability

Brain microstructure plays a key role in driving the transport of drug molecules directly administered to the brain tissue, as in Convection-Enhanced Delivery procedures. The proposed research analyzes the hydraulic permeability of two white matter (WM) areas (corpus callosum and fornix) whose three-dimensional microstructure was reconstructed starting from the acquisition of electron microscopy images. We cut the two volumes with 20 equally spaced planes distributed along two perpendicular directions, and, on each plane, we computed the corresponding permeability vector. Then, we considered that the WM structure is mainly composed of elongated and parallel axons, and, using a principal component analysis, we defined two principal directions, parallel and perpendicular, with respect to the axons' main direction. The latter were used to define a reference frame onto which the permeability vectors were projected to finally obtain the permeability along the parallel and perpendicular directions. The results show a statistically significant difference between parallel and perpendicular permeability, with a ratio of about two in both the WM structures analyzed, thus demonstrating their anisotropic behavior. Moreover, we find a significant difference between permeability in corpus callosum and fornix, which suggests that the WM heterogeneity should also be considered when modeling drug transport in the brain. Our findings, which demonstrate and quantify the anisotropic and heterogeneous character of the WM, represent a fundamental contribution not only for drug-delivery modeling, but also for shedding light on the interstitial transport mechanisms in the extracellular space.

Convection Enhanced Delivery in the Setting of High-Grade Gliomas

Development of effective treatments for high-grade glioma (HGG) is hampered by (1) the blood–brain barrier (BBB), (2) an infiltrative growth pattern, (3) rapid development of therapeutic resistance, and, in many cases, (4) dose-limiting toxicity due to systemic exposure. Convection-enhanced delivery (CED) has the potential to significantly limit systemic toxicity and increase therapeutic index by directly delivering homogenous drug concentrations to the site of disease. In this review, we present clinical experiences and preclinical developments of CED in the setting of high-grade gliomas.