Nonuniform Moving Boundary Method for Computational Fluid Dynamics Simulation of Intrathecal Cerebrospinal Flow Distribution in a Cynomolgus Monkey

Harnessing cerebrospinal fluid circulation for drug delivery to brain tissues

Initially thought to be useful only to reach tissues in the immediate vicinity of the CSF circulatory system, CSF circulation is now increasingly viewed as a viable pathway to deliver certain therapeutics deeper into brain tissues. There is emerging evidence that this goal is achievable in the case of large therapeutic proteins, provided conditions are met that are described herein. We show how fluid dynamic modeling helps predict infusion rate and duration to overcome high CSF turnover. We posit that despite model limitations and controversies, fluid dynamic models, pharmacokinetic models, preclinical testing, and a qualitative understanding of the glymphatic system circulation can be used to estimate drug penetration in brain tissues. Lastly, in addition to highlighting landmark scientific and medical literature, we provide practical advice on formulation development, device selection, and pharmacokinetic modeling. Our review of clinical studies suggests a growing interest for intra-CSF delivery, particularly for targeted proteins.

Intrathecal catheter implantation decreases cerebrospinal fluid dynamics in cynomolgus monkeys

A detailed understanding of the CSF dynamics is needed for design and optimization of intrathecal drug delivery devices, drugs, and protocols. Preclinical research using large-animal models is important to help define drug pharmacokinetics-pharmacodynamics and safety. In this study, we investigated the impact of catheter implantation in the sub-dural space on CSF flow dynamics in Cynomolgus monkeys. Magnetic resonance imaging (MRI) was performed before and after catheter implantation to quantify the differences based on catheter placement location in the cervical compared to the lumbar spine. Several geometric and hydrodynamic parameters were calculated based on the 3D segmentation and flow analysis. Hagen-Poiseuille equation was used to investigate the impact of catheter implantation on flow reduction and hydraulic resistance. A linear mixed-effects model was used in this study to investigate if there was a statistically significant difference between cervical and lumbar implantation, or between two MRI time points. Results showed that geometric parameters did not change statistically across MRI measurement time points and did not depend on catheter location. However, catheter insertion did have a significant impact on the hydrodynamic parameters and the effect was greater with cervical implantation compared to lumbar implantation. CSF flow rate decreased up to 55% with the catheter located in the cervical region. The maximum flow rate reduction in the lumbar implantation group was 21%. Overall, lumbar catheter implantation disrupted CSF dynamics to a lesser degree than cervical catheter implantation and this effect remained up to two weeks post-catheter implantation in Cynomolgus monkeys.

In vitro and numerical simulation of blood removal from cerebrospinal fluid: comparison of lumbar drain to Neurapheresis therapy

BACKGROUND

Blood removal from cerebrospinal fluid (CSF) in post-subarachnoid hemorrhage patients may reduce the risk of related secondary brain injury. We formulated a computational fluid dynamics (CFD) model to investigate the impact of a dual-lumen catheter-based CSF filtration system, called Neurapheresis™ therapy, on blood removal from CSF compared to lumbar drain.

METHODS

A subject-specific multiphase CFD model of CSF system-wide solute transport was constructed based on MRI measurements. The Neurapheresis catheter geometry was added to the model within the spinal subarachnoid space (SAS). Neurapheresis flow aspiration and return rate was 2.0 and 1.8 mL/min, versus 0.2 mL/min drainage for lumbar drain. Blood was modeled as a bulk fluid phase within CSF with a 10% initial tracer concentration and identical viscosity and density as CSF. Subject-specific oscillatory CSF flow was applied at the model inlet. The dura and spinal cord geometry were considered to be stationary. Spatial-temporal tracer concentration was quantified based on time-average steady-streaming velocities throughout the domain under Neurapheresis therapy and lumbar drain. To help verify CFD results, an optically clear in vitro CSF model was constructed with fluorescein used as a blood surrogate. Quantitative comparison of numerical and in vitro results was performed by linear regression of spatial-temporal tracer concentration over 24-h.

RESULTS

After 24-h, tracer concentration was reduced to 4.9% under Neurapheresis therapy compared to 6.5% under lumbar drain. Tracer clearance was most rapid between the catheter aspiration and return ports. Neurapheresis therapy was found to have a greater impact on steady-streaming compared to lumbar drain. Steady-streaming in the cranial SAS was ~ 50× smaller than in the spinal SAS for both cases. CFD results were strongly correlated with the in vitro spatial-temporal tracer concentration under Neurapheresis therapy (R² = 0.89 with + 2.13% and - 1.93% tracer concentration confidence interval).

CONCLUSION

A subject-specific CFD model of CSF system-wide solute transport was used to investigate the impact of Neurapheresis therapy on tracer removal from CSF compared to lumbar drain over a 24-h period. Neurapheresis therapy was found to substantially increase tracer clearance compared to lumbar drain. The multiphase CFD results were verified by in vitro fluorescein tracer experiments.

Impact of Neurapheresis System on Intrathecal Cerebrospinal Fluid Dynamics: A Computational Fluid Dynamics Study

It has been hypothesized that early and rapid filtration of blood from cerebrospinal fluid (CSF) in postsubarachnoid hemorrhage patients may reduce hospital stay and related adverse events. In this study, we formulated a subject-specific computational fluid dynamics (CFD) model to parametrically investigate the impact of a novel dual-lumen catheter-based CSF filtration system, the Neurapheresis™ system (Minnetronix Neuro, Inc., St. Paul, MN), on intrathecal CSF dynamics. The operating principle of this system is to remove CSF from one location along the spine (aspiration port), externally filter the CSF routing the retentate to a waste bag, and return permeate (uncontaminated CSF) to another location along the spine (return port). The CFD model allowed parametric simulation of how the Neurapheresis system impacts intrathecal CSF velocities and steady-steady streaming under various Neurapheresis flow settings ranging from 0.5 to 2.0 ml/min and with a constant retentate removal rate of 0.2 ml/min simulation of the Neurapheresis system were compared to a lumbar drain simulation with a typical CSF removal rate setting of 0.2 ml/min. Results showed that the Neurapheresis system at a maximum flow of 2.0 ml/min increased average steady streaming CSF velocity 2× in comparison to lumbar drain (0.190 ± 0.133 versus 0.093 ± 0.107 mm/s, respectively). This affect was localized to the region within the Neurapheresis flow loop. The mean velocities introduced by the flow loop were relatively small in comparison to normal cardiac-induced CSF velocities.

Non-invasive MRI quantification of cerebrospinal fluid dynamics in amyotrophic lateral sclerosis patients

Background

Developing novel therapeutic agents to treat amyotrophic lateral sclerosis (ALS) has been difficult due to multifactorial pathophysiologic processes at work. Intrathecal drug administration shows promise due to close proximity of cerebrospinal fluid (CSF) to affected tissues. Development of effective intrathecal pharmaceuticals will rely on accurate models of how drugs are dispersed in the CSF. Therefore, a method to quantify these dynamics and a characterization of differences across disease states is needed.

Methods

Complete intrathecal 3D CSF geometry and CSF flow velocities at six axial locations in the spinal canal were collected by T2-weighted and phase-contrast MRI, respectively. Scans were completed for eight people with ALS and ten healthy controls. Manual segmentation of the spinal subarachnoid space was performed and coupled with an interpolated model of CSF flow within the spinal canal. Geometric and hydrodynamic parameters were then generated at 1 mm slice intervals along the entire spine. Temporal analysis of the waveform spectral content and feature points was also completed.

Results

Comparison of ALS and control groups revealed a reduction in CSF flow magnitude and increased flow propagation velocities in the ALS cohort. Other differences in spectral harmonic content and geometric comparisons may support an overall decrease in intrathecal compliance in the ALS group. Notably, there was a high degree of variability between cases, with one ALS patient displaying nearly zero CSF flow along the entire spinal canal.

Conclusion

While our sample size limits statistical confidence about the differences observed in this study, it was possible to measure and quantify inter-individual and cohort variability in a non-invasive manner. Our study also shows the potential for MRI based measurements of CSF geometry and flow to provide information about  the hydrodynamic environment of the spinal subarachnoid space. These dynamics may be studied further to understand the behavior of CSF solute transport in healthy and diseased states.

Anthropomorphic Model of Intrathecal Cerebrospinal Fluid Dynamics Within the Spinal Subarachnoid Space: Spinal Cord Nerve Roots Increase Steady-Streaming

Cerebrospinal fluid (CSF) dynamics are thought to play a vital role in central nervous system (CNS) physiology. The objective of this study was to investigate the impact of spinal cord (SC) nerve roots (NR) on CSF dynamics. A subject-specific computational fluid dynamics (CFD) model of the complete spinal subarachnoid space (SSS) with and without anatomically realistic NR and nonuniform moving dura wall deformation was constructed. This CFD model allowed detailed investigation of the impact of NR on CSF velocities that is not possible in vivo using magnetic resonance imaging (MRI) or other noninvasive imaging methods. Results showed that NR altered CSF dynamics in terms of velocity field, steady-streaming, and vortical structures. Vortices occurred in the cervical spine around NR during CSF flow reversal. The magnitude of steady-streaming CSF flow increased with NR, in particular within the cervical spine. This increase was located axially upstream and downstream of NR due to the interface of adjacent vortices that formed around NR.

Characterization of intrathecal cerebrospinal fluid geometry and dynamics in cynomolgus monkeys (macaca fascicularis) by magnetic resonance imaging

Recent advancements have been made toward understanding the diagnostic and therapeutic potential of cerebrospinal fluid (CSF) and related hydrodynamics. Increased understanding of CSF dynamics may lead to improved detection of central nervous system (CNS) diseases and optimized delivery of CSF based CNS therapeutics, with many proposed therapeutics hoping to successfully treat or cure debilitating neurological conditions. Before significant strides can be made toward the research and development of interventions designed for human use, additional research must be carried out with representative subjects such as non-human primates (NHP). This study presents a geometric and hydrodynamic characterization of CSF in eight cynomolgus monkeys (Macaca fascicularis) at baseline and two-week follow-up.

Results showed that CSF flow along the entire spine was laminar with a Reynolds number ranging up to 80 and average Womersley number ranging from 4.1–7.7. Maximum CSF flow rate occurred ~25 mm caudal to the foramen magnum. Peak CSF flow rate ranged from 0.3–0.6 ml/s at the C3-C4 level. Geometric analysis indicated that average intrathecal CSF volume below the foramen magnum was 7.4 ml. The average surface area of the spinal cord and dura was 44.7 and 66.7 cm² respectively. Subarachnoid space cross-sectional area and hydraulic diameter ranged from 7–75 mm² and 2–3.7 mm, respectively. Stroke volume had the greatest value of 0.14 ml at an axial location corresponding to C3-C4.

Computer simulation of syringomyelia in dogs

Background

Syringomyelia is a pathological condition in which fluid-filled cavities (syringes) form and expand in the spinal cord. Syringomyelia is often linked with obstruction of the craniocervical junction and a Chiari malformation, which is similar in both humans and animals. Some brachycephalic toy breed dogs such as Cavalier King Charles Spaniels (CKCS) are particularly predisposed. The exact mechanism of the formation of syringomyelia is undetermined and consequently with the lack of clinical explanation, engineers and mathematicians have resorted to computer models to identify possible physical mechanisms that can lead to syringes. We developed a computer model of the spinal cavity of a CKCS suffering from a large syrinx. The model was excited at the cranial end to simulate the movement of the cerebrospinal fluid (CSF) and the spinal cord due to the shift of blood volume in the cranium related to the cardiac cycle. To simulate the normal condition, the movement was prescribed to the CSF. To simulate the pathological condition, the movement of CSF was blocked.

Results

For normal conditions the pressure in the SAS was approximately 400 Pa and the same applied to all stress components in the spinal cord. The stress was uniformly distributed along the length of the spinal cord. When the blockage between the cranial and spinal CSF spaces forced the cord to move with the cardiac cycle, shear and axial normal stresses in the cord increased significantly. The sites where the elevated stress was most pronounced coincided with the axial locations where the syringes typically form, but they were at the perimeter rather than in the central portion of the cord. This elevated stress originated from the bending of the cord at the locations where its curvature was high.

Conclusions

The results suggest that it is possible that repetitive stressing of the spinal cord caused by its exaggerated movement could be a cause for the formation of initial syringes. Further consideration of factors such as cord tethering and the difference in mechanical properties of white and grey matter is needed to fully explore this possibility.

Electronic supplementary material

The online version of this article (10.1186/s12917-018-1410-7) contains supplementary material, which is available to authorized users.

A 3D subject-specific model of the spinal subarachnoid space with anatomically realistic ventral and dorsal spinal cord nerve rootlets

BACKGROUND

The spinal subarachnoid space (SSS) has a complex 3D fluid-filled geometry with multiple levels of anatomic complexity, the most salient features being the spinal cord and dorsal and ventral nerve rootlets. An accurate anthropomorphic representation of these features is needed for development of in vitro and numerical models of cerebrospinal fluid (CSF) dynamics that can be used to inform and optimize CSF-based therapeutics.

METHODS

A subject-specific 3D model of the SSS was constructed based on high-resolution anatomic MRI. An expert operator completed manual segmentation of the CSF space with detailed consideration of the anatomy. 31 pairs of semi-idealized dorsal and ventral nerve rootlets (NR) were added to the model based on anatomic reference to the magnetic resonance (MR) imaging and cadaveric measurements in the literature. Key design criteria for each NR pair included the radicular line, descending angle, number of NR, attachment location along the spinal cord and exit through the dura mater. Model simplification and smoothing was performed to produce a final model with minimum vertices while maintaining minimum error between the original segmentation and final design. Final model geometry and hydrodynamics were characterized in terms of axial distribution of Reynolds number, Womersley number, hydraulic diameter, cross-sectional area and perimeter.

RESULTS

The final model had a total of 139,901 vertices with a total CSF volume within the SSS of 97.3 cm³. Volume of the dura mater, spinal cord and NR was 123.1, 19.9 and 5.8 cm³. Surface area of these features was 318.52, 112.2 and 232.1 cm² respectively. Maximum Reynolds number was 174.9 and average Womersley number was 9.6, likely indicating presence of a laminar inertia-dominated oscillatory CSF flow field.

CONCLUSIONS

This study details an anatomically realistic anthropomorphic 3D model of the SSS based on high-resolution MR imaging of a healthy human adult female. The model is provided for re-use under the Creative Commons Attribution-ShareAlike 4.0 International license (CC BY-SA 4.0) and can be used as a tool for development of in vitro and numerical models of CSF dynamics for design and optimization of intrathecal therapeutics.