In Vivo Intrathecal Tracer Dispersion in Cynomolgus Monkey Validates Wide Biodistribution Along Neuraxis.

OBJECTIVE: It is commonly believed that in intrathecal (IT) drug delivery, agent distribution is confined to a narrow region close to the injection site, thereby undermining the efficacy of the method. METHODS: To test the claim, multimodal in vivo imaging was used to experimentally observe the effects of IT infusion in cynomolgus monkey, looking at cerebrospinal fluid flow, anatomy, and dispersion of a radiolabeled tracer. RESULTS: At high infusion rates, the tracer reached the cervical region after only 2 h, demonstrating rapid and wide distribution. The same in vivo nonhuman primate imaging data also provided evidence in support of a computational fluid dynamic model for the prediction of drug distribution following IT injection. Tracer dispersion was predicted in two specimens matching the distribution acquired with positron emission tomography (PET). For the third specimen, tracer dispersion simulations were conducted as a blind study: predictions were made before in vivo biodistribution data was known. In all cases, the computational fluid dynamics (CFD) predictions of drug dispersion after IT administration showed close spatio-temporal agreement with tracer biodistribution in vivo. CONCLUSION: Validation by in vivo nonhuman primate data confirms our ability to predict the biodistribution of intrathecally administered agents in subject-specific models of the central nervous system from first principles. SIGNIFICANCE: The experiments reinstate IT delivery as a viable administration method when targeting molecules to the whole spine or the brain. The proposed computational methodology enables rational design of novel therapies for neurological diseases that require reliable, efficient, and safe delivery of therapeutic agents to specific target sites in the central nervous system.

Cellular Obstruction Clearance in Proximal Ventricular Catheters Using Low-Voltage Joule Heating.

OBJECTIVE: Proximal obstruction due to cellular material is a major cause of shunt failure in hydrocephalus management. The standard approach to treat such cases involves surgical intervention which unfortunately is accompanied by inherent surgical risks and a likelihood of future malfunction. We report a prototype design of a proximal ventricular catheter capable of noninvasively clearing cellular obstruction. Methods: In-vitro cell-culture methods show that low-intensity ac signals successfully destroy a cellular layer in a localized manner by means of Joule heating induced hyperthermia. A detailed electrochemical model for determining the temperature distribution and ionic current density for an implanted ventricular catheter supports our experimental observations. RESULTS: In-vitro experiments with cells cultured in a plate as well as cells seeded in mock ventricular catheters demonstrated that localized heating between 43 °C and 48 °C caused cell death. This temperature range is consistent with hyperthermia. The electrochemical model verified that Joule heating due to ionic motion is the primary contributor to heat generation. CONCLUSION: Hyperthermia induced by Joule heating can clear cellular material in a localized manner. This approach is feasible to design a noninvasive self-clearing ventricular catheter system. SIGNIFICANCE: A shunt system capable of clearing cellular obstruction could significantly reduce the need for future surgical interventions, lower the cost of disease management, and improve the quality of life for patients suffering from hydrocephalus.

Modeling the diffusion of D-2-hydroxyglutarate from IDH1 mutant gliomas in the central nervous system.

Background: Among diffusely infiltrative gliomas in adults, 20%-30% contain a point mutation in isocitrate dehydrogenase 1 (IDH1mut), which increases production of D-2-hydroxyglutarate (D2HG). This is so efficient that D2HG often reaches 30 mM within IDH1mut gliomas. Yet, while up to 100 µM D2HG can be detected in the circulating cerebrospinal fluid of IDH1mut glioma patients, the exposure of nonneoplastic cells within and surrounding an IDH1mut glioma to D2HG is unknown and difficult to measure directly. Methods: Conditioned medium from patient-derived wild type IDH1 (IDH1wt) and IDH1mut glioma cells was analyzed for D2HG by liquid chromatography-mass spectrometry (LC-MS). Mathematical models of D2HG release and diffusion around an IDH1mut glioma were independently generated based on fluid dynamics within the brain and on previously reported intratumoral and cerebrospinal D2HG concentrations. Results: LC-MS analysis indicates that patient-derived IDH1mut glioma cells release 3.7-97.0 pg D2HG per cell per week. Extrapolating this to an average-sized tumor (30 mL glioma volume and 1 × 108 cells/mL tumor), the rate of D2HG release by an IDH1mut glioma (SA) is estimated at 3.2-83.0 × 10-12 mol/mL/sec. Mathematical models estimate an SA of 2.9-12.9 × 10-12 mol/mL/sec, within the range of the in vitro LC-MS data. In even the most conservative of these models, the extracellular concentration of D2HG exceeds 3 mM within a 2 cm radius from the center of an IDH1mut glioma. Conclusions: The microenvironment of an IDH1mut glioma is likely being exposed to high concentrations of D2HG, in the low millimolar range. This has implications for understanding how D2HG affects nonneoplastic cells in an IDH1mut glioma.

Large-scale subject-specific cerebral arterial tree modeling using automated parametric mesh generation for blood flow simulation.

In this paper, we present a novel technique for automatic parametric mesh generation of subject-specific cerebral arterial trees. This technique generates high-quality and anatomically accurate computational meshes for fast blood flow simulations extending the scope of 3D vascular modeling to a large portion of cerebral arterial trees. For this purpose, a parametric meshing procedure was developed to automatically decompose the vascular skeleton, extract geometric features and generate hexahedral meshes using a body-fitted coordinate system that optimally follows the vascular network topology. To validate the anatomical accuracy of the reconstructed vasculature, we performed statistical analysis to quantify the alignment between parametric meshes and raw vascular images using receiver operating characteristic curve. Geometric accuracy evaluation showed an agreement with area under the curves value of 0.87 between the constructed mesh and raw MRA data sets. Parametric meshing yielded on-average, 36.6% and 21.7% orthogonal and equiangular skew quality improvement over the unstructured tetrahedral meshes. The parametric meshing and processing pipeline constitutes an automated technique to reconstruct and simulate blood flow throughout a large portion of the cerebral arterial tree down to the level of pial vessels. This study is the first step towards fast large-scale subject-specific hemodynamic analysis for clinical applications.

Backflow-free catheters for efficient and safe convection-enhanced delivery of therapeutics.

Convection-enhanced delivery (CED) is an invasive drug delivery technique used to target specific regions of the brain for the treatment of cancer and neurodegenerative diseases while bypassing the blood-brain barrier. In order to prevent the possibility of backflow, low volumetric flow rates are applied which limit the achievable drug distribution volumes from CED. This can render CED treatment ineffective since a small convective flow produces narrow drug distribution inside the treatment region. Novel catheter designs and CED protocols are needed to improve the drug distribution inside the treatment region. This is especially important when administering toxic chemotherapeutics which could adversely affect other organs if backflow occurred and these drugs entered the circulating blood stream. In order to help elucidate the causes of backflow and to design backflow-free catheters, we have studied the impact that microfluid flow has on deformable brain phantom gels experimentally as well as numerically. We found that fluid injections into porous media have considerable effects on local transport properties such as porosity and hydraulic conductivity. These phenomena not only alter the bulk flow velocity distribution of the microfluid flow due to the changing porosity, but significantly modify flow direction and even volumetric flow distribution due to induced local hydraulic conductivity anisotropy. These studies led us to the development of novel backflow-free catheters with safe volumetric flow rates up to 10 µL/min. The catheter designs, numerical simulations and experimental results are described throughout this article.

Computational and In Vitro Experimental Investigation of Intrathecal Drug Distribution: Parametric Study of the Effect of Injection Volume, Cerebrospinal Fluid Pulsatility, and Drug Uptake.

BACKGROUND: Intrathecal drug delivery is an attractive option to circumvent the blood-brain barrier for pain management through its increased efficacy of pain relief, reduction in adverse side effects, and cost-effectiveness. Unfortunately, there are limited guidelines for physicians to choose infusion or drug pump settings to administer therapeutic doses to specific regions of the spine or the brain. Although empiric trialing of intrathecal drugs is critical to determine the sustained side effects, currently there is no inexpensive in vitro method to guide the selection of spinal drug delivery parameters. The goal of this study is to demonstrate current computational capabilities to predict drug biodistribution while varying 3 parameters: (1) infusion settings, (2) drug chemistry, and (3) subject-specific anatomy and cerebrospinal fluid dynamics. We will discuss strategies to systematically optimize these 3 parameters to administer drug molecules to targeted tissue locations in the central nervous system. METHODS: We acquired anatomical data from magnetic resonance imaging (MRI) and velocity measurements in the spinal cerebrospinal fluid with CINE-MRI for 2 subjects. A bench-top surrogate of the subject-specific central nervous system was constructed to match measured anatomical dimensions and volumes. We generated a computational mesh for the bench-top model. Idealized simulations of tracer distribution were compared with bench-top measurements for validation. Using reconstructions from MRI data, we also introduced a subject-specific computer model for predicting drug spread for the human volunteer. RESULTS: MRI velocity measurements at 3 spinal regions of interest reasonably matched the simulated flow fields in a subject-specific computer mesh. Comparison between the idealized spine computations and bench-top tracer distribution experiments demonstrate agreement of our drug transport predictions to this physical model. Simulated multibolus drug infusion theoretically localizes drug to the cervical and thoracic region. Continuous drug pump and single bolus injection were successful to target the lumbar spine in the simulations. The parenchyma might be targeted suitably by multiple boluses followed by a flush infusion. We present potential guidelines that take into account drug specific kinetics for tissue uptake, which influence the speed of drug dispersion in the model and potentially influence tissue targeting. CONCLUSIONS: We present potential guidelines considering drug-specific kinetics of tissue uptake, which determine the speed of drug dispersion and influence tissue targeting. However, there are limitations to this analysis in that the parameters were obtained from an idealized healthy patient in a supine position. The proposed methodology could assist physicians to select clinical infusion parameters for their patients and provide guidance to optimize treatment algorithms. In silico optimization of intrathecal drug delivery therapies presents the first steps toward a possible care paradigm in the future that is specific to personalized patient anatomy and diseases.

Starling forces drive intracranial water exchange during normal and pathological states.

AIM: To quantify the exchange of water between cerebral compartments, specifically blood, tissue, perivascular pathways, and cerebrospinal fluid-filled spaces, on the basis of experimental data and to propose a dynamic global model of water flux through the entire brain to elucidate functionally relevant fluid exchange phenomena. METHODS: The mechanistic computer model to predict brain water shifts is discretized by cerebral compartments into nodes. Water and species flux is calculated between these nodes across a network of arcs driven by Hagen-Poiseuille flow (blood), Darcy flow (interstitial fluid transport), and Starling's Law (transmembrane fluid exchange). Compartment compliance is accounted for using a pressure-volume relationship to enforce the Monro-Kellie doctrine. This nonlinear system of differential equations is solved implicitly using MATLAB software. RESULTS: The model predictions of intraventricular osmotic injection caused a pressure rise from 10 to 22 mmHg, followed by a taper to 14 mmHg over 100 minutes. The computational results are compared to experimental data with R2=0.929. Moreover, simulated osmotic therapy of systemic (blood) injection reduced intracranial pressure from 25 to 10 mmHg. The modeled volume and intracranial pressure changes following cerebral edema agree with experimental trends observed in animal models with R2=0.997. CONCLUSION: The model successfully predicted time course and the efficacy of osmotic therapy for clearing cerebral edema. Furthermore, the mathematical model implicated the perivascular pathways as a possible conduit for water and solute exchange. This was a first step to quantify fluid exchange throughout the brain.

CNS wide simulation of flow resistance and drug transport due to spinal microanatomy.

Spinal microstructures are known to substantially affect cerebrospinal fluid patterns, yet their actual impact on flow resistance has not been quantified. Because the length scale of microanatomical aspects is below medical image resolution, their effect on flow is difficult to observe experimentally. Using a computational fluid mechanics approach, we were able to quantify the contribution of micro-anatomical aspects on cerebrospinal fluid (CSF) flow patterns and flow resistance within the entire central nervous system (CNS). Cranial and spinal CSF filled compartments were reconstructed from human imaging data; microscopic trabeculae below the image detection threshold were added artificially. Nerve roots and trabeculae were found to induce regions of microcirculation, whose location, size and vorticity along the spine were characterized. Our CFD simulations based on volumetric flow rates acquired with Cine Phase Contrast MRI in a normal human subject suggest a 2-2.5 fold increase in pressure drop mainly due to arachnoid trabeculae. The timing and phase lag of the CSF pressure and velocity waves along the spinal canal were also computed, and a complete spatio-temporal map encoding CSF volumetric flow rates and pressure was created. Micro-anatomy induced fluid patterns were found responsible for the rapid caudo-cranial spread of an intrathecally administered drug. The speed of rostral drug dispersion is drastically accelerated through pulsatile flow around microanatomy induced vortices. Exploring massive parallelization on a supercomputer, the feasibility of computational drug transport studies was demonstrated. CNS-wide simulations of intrathecal drugs administration can become a practical tool for in silico design, interspecies scaling and optimization of experimental drug trials.

Impedance Changes Indicate Proximal Ventriculoperitoneal Shunt Obstruction In Vitro.

Extracranial cerebrospinal fluid (CSF) shunt obstruction is one of the most important problems in hydrocephalus patient management. Despite ongoing research into better shunt design, robust and reliable detection of shunt malfunction remains elusive. The authors present a novel method of correlating degree of tissue ingrowth into ventricular CSF drainage catheters with internal electrical impedance. The impedance based sensor is able to continuously monitor shunt patency using intraluminal electrodes. Prototype obstruction sensors were fabricated for in-vitro analysis of cellular ingrowth into a shunt under static and dynamic flow conditions. Primary astrocyte cell lines and C6 glioma cells were allowed to proliferate up to 7 days within a shunt catheter and the impedance waveform was observed. During cell ingrowth a significant change in the peak-to-peak voltage signal as well as the root-mean-square voltage level was observed, allowing the impedance sensor to potentially anticipate shunt malfunction long before it affects fluid drainage. Finite element modeling was employed to demonstrate that the electrical signal used to monitor tissue ingrowth is contained inside the catheter lumen and does not endanger tissue surrounding the shunt. These results may herald the development of "next generation" shunt technology that allows prediction of malfunction before it affects patient outcome.