A numerical investigation of intrathecal isobaric drug dispersion within the cervical subarachnoid space

Magnetically Guided Microcatheter for Targeted Injection of Magnetic Particle Swarms

The initial delivery of small-scale magnetic devices such as microrobots is a key, but often overlooked, aspect for their use in clinical applications. The deployment of these devices within the dynamic environment of the human body presents significant challenges due to their dispersion caused by circulatory flows. Here, a method is introduced to effectively deliver a swarm of magnetic nanoparticles in fluidic flows. This approach integrates a magnetically navigated robotic microcatheter equipped with a reservoir for storing the magnetic nanoparticles. The microfluidic flow within the reservoir facilitates the injection of magnetic nanoparticles into the fluid stream, and a magnetic field gradient guides the swarm through the oscillatory flow to a target site. The microcatheter and reservoir are engineered to enable magnetic steering and injection of the magnetic nanoparticles. To demonstrate this approach, experiments are conducted utilizing a spinal cord phantom simulating intrathecal catheter delivery for applications in the central nervous system. These results demonstrate that the proposed microcatheter successfully concentrates nanoparticles near the desired location through the precise manipulation of magnetic field gradients, offering a promising solution for the controlled deployment of untethered magnetic micro-/nanodevices within the complex physiological circulatory systems of the human body.

Effects of buoyancy on the dispersion of drugs released intrathecally in the spinal canal

This paper investigates the transport of drugs delivered by direct injection into the cerebrospinal fluid (CSF) that fills the intrathecal space surrounding the spinal cord. Because of the small drug diffusivity, the dispersion of neutrally buoyant drugs has been shown in previous work to rely mainly on the mean Lagrangian flow associated with the CSF oscillatory motion. Attention is given here to effects of buoyancy, arising when the drug density differs from the CSF density. For the typical density differences found in applications, the associated Richardson number is shown to be of order unity, so that the Lagrangian drift includes a buoyancy-induced component that depends on the spatial distribution of the drug, resulting in a slowly evolving cycle-averaged flow problem that can be analysed with two-time scale methods. The asymptotic analysis leads to a nonlinear integro-differential equation for the spatiotemporal solute evolution that describes accurately drug dispersion at a fraction of the cost involved in direct numerical simulations of the oscillatory flow. The model equation is used to predict drug dispersion of positively and negatively buoyant drugs in an anatomically correct spinal canal, with separate attention given to drug delivery via bolus injection and constant infusion.

In vivo distribution of cerebrospinal fluid tracer in human upper spinal cord and brain stem

BACKGROUNDIntrathecal injection is an attractive route through which drugs can be administered and directed to the spinal cord, restricted by the blood-spinal cord barrier. However, in vivo data on the distribution of cerebrospinal fluid (CSF) substances in the human spinal cord are lacking. We conducted this study to assess the enrichment of a CSF tracer in the upper cervical spinal cord and the brain stem.METHODSAfter lumbar intrathecal injection of a magnetic resonance imaging (MRI) contrast agent, gadobutrol, repeated blood samples and MRI of the upper cervical spinal cord, brain stem, and adjacent subarachnoid spaces (SAS) were obtained through 48 hours. The MRI scans were then analyzed for tracer distribution in the different regions and correlated to age, disease, and amounts of tracer in the blood to determine CSF-to-blood clearance.RESULTSThe study included 26 reference individuals and 35 patients with the dementia subtype idiopathic normal pressure hydrocephalus (iNPH). The tracer enriched all analyzed regions. Moreover, tracer enrichment in parenchyma was associated with tracer enrichment in the adjacent SAS and with CSF-to-blood clearance. Clearance from the CSF was delayed in patients with iNPH compared with younger reference patients.CONCLUSIONA CSF tracer substance administered to the lumbar thecal sac can access the parenchyma of the upper cervical spinal cord and brain stem. Since CSF-to-blood clearance is highly individual and is associated with tracer level in CSF, clearance assessment may be used to tailor intrathecal treatment regimes.FUNDINGSouth-Eastern Norway Regional Health and Østfold Hospital Trust supported the research and publication of this work.

Determination of spinal tracer dispersion after intrathecal injection in a deformable CNS model

Background: Traditionally, there is a widely held belief that drug dispersion after intrathecal (IT) delivery is confined locally near the injection site. We posit that high-volume infusions can overcome this perceived limitation of IT administration. Methods: To test our hypothesis, subject-specific deformable phantom models of the human central nervous system were manufactured so that tracer infusion could be realistically replicated in vitro over the entire physiological range of pulsating cerebrospinal fluid (CSF) amplitudes and frequencies. The distribution of IT injected tracers was studied systematically with high-speed optical methods to determine its dependence on injection parameters (infusion volume, flow rate, and catheter configurations) and natural CSF oscillations in a deformable model of the central nervous system (CNS). Results: Optical imaging analysis of high-volume infusion experiments showed that tracers spread quickly throughout the spinal subarachnoid space, reaching the cervical region in less than 10 min. The experimentally observed biodispersion is much slower than suggested by the Taylor-Aris dispersion theory. Our experiments indicate that micro-mixing patterns induced by oscillatory CSF flow around microanatomical features such as nerve roots significantly accelerate solute transport. Strong micro-mixing effects due to anatomical features in the spinal subarachnoid space were found to be active in intrathecal drug administration but were not considered in prior dispersion theories. Their omission explains why prior models developed in the engineering community are poor predictors for IT delivery. Conclusion: Our experiments support the feasibility of targeting large sections of the neuroaxis or brain utilizing high-volume IT injection protocols. The experimental tracer dispersion profiles acquired with an anatomically accurate, deformable, and closed in vitro human CNS analog informed a new predictive model of tracer dispersion as a function of physiological CSF pulsations and adjustable infusion parameters. The ability to predict spatiotemporal dispersion patterns is an essential prerequisite for exploring new indications of IT drug delivery that targets specific regions in the CNS or the brain.

Intrathecal delivery and its applications in leptomeningeal disease

Intrathecal delivery (IT) of opiates into the cerebrospinal fluid (CSF) for anesthesia and pain relief has been used clinically for decades, but this relatively straightforward approach of bypassing the blood-brain barrier has been underutilized for other indications because of its lack of utility in delivering small lipid-soluble drugs. However, emerging evidence suggests that IT drug delivery be an efficacious strategy for the treatment of cancers in which there is leptomeningeal spread of disease. In this review, we discuss CSF flow dynamics and CSF clearance pathways in the context of intrathecal delivery. We discuss human and animal studies of several new classes of therapeutic agents-cellular, protein, nucleic acid, and nanoparticle-based small molecules-that may benefit from IT delivery. The complexity of the CSF compartment presents several key challenges in predicting biodistribution of IT-delivered drugs. New approaches and strategies are needed that can overcome the high rates of turnover in the CSF to reach specific tissues or cellular targets.

Investigation of Human Intrathecal Solute Transport Dynamics Using a Novel in vitro Cerebrospinal Fluid System Analog

BACKGROUND

Understanding the relationship between cerebrospinal fluid (CSF) dynamics and intrathecal drug delivery (ITDD) injection parameters is essential to improve treatment of central nervous system (CNS) disorders.

METHODS

An anatomically detailed in vitro model of the complete CSF system was constructed. Patient-specific cardiac- and respiratory-induced CSF oscillations were input to the model in the subarachnoid space and within the ventricles. CSF production was input at the lateral ventricles and CSF absorption at the superior sagittal sinus. A model small molecule simulated drug product containing fluorescein was imaged within the system over a period of 3-h post-lumbar ITDD injections and used to quantify the impact of (a) bolus injection volume and rate, (b) post-injection flush volume, rate, and timing, (c) injection location, and (d) type of injection device. For each experiment, neuraxial distribution of fluorescein in terms of spatial temporal concentration, area-under-the-curve (AUC), and percent of injected dose (%ID) to the brain was quantified at a time point 3-h post-injection.

RESULTS

For all experiments conducted with ITDD administration in the lumbar spine, %ID to the brain did not exceed 11.6% at a time point 3-h post-injection. Addition of a 12 mL flush slightly increased solute transport to the brain up to +3.9%ID compared to without a flush (p < 0.01). Implantation of a lumbar catheter with the tip at an equivalent location to the lumbar placed needle, but with rostral tip orientation, resulted in a small improvement of 1.5%ID to the brain (p < 0.05). An increase of bolus volume from 5 to 20 mL improved solute transport to the brain from 5.0 to 6.3%ID, but this improvement was not statistically significant. Increasing bolus injection rate from 5 to 13.3 mL/min lacked improvement of solute transport to the brain, with a value of 6.3 compared to 5.7%ID.

CONCLUSION

The in vitro modeling approach allowed precisely controlled and repeatable parametric investigation of ITDD injection protocols and devices. In combination, the results predict that parametric changes in lumbar spine ITDD-injection related parameters and devices can alter %ID to the brain and be tuned to optimize therapeutic benefit to CNS targets.

Intrathecal amyloid-beta oligomer administration increases tau phosphorylation in the medial temporal lobe in the African green monkey: A nonhuman primate model of Alzheimer's disease

AIMS

An obstacle to developing new treatment strategies for Alzheimer's disease (AD) has been the inadequate translation of findings in current AD transgenic rodent models to the prediction of clinical outcomes. By contrast, nonhuman primates (NHPs) share a close neurobiology with humans in virtually all aspects relevant to developing a translational AD model. The present investigation used African green monkeys (AGMs) to refine an inducible NHP model of AD based on the administration of amyloid-beta oligomers (AβOs), a key upstream initiator of AD pathology.

METHODS

AβOs or vehicle were repeatedly delivered over 4 weeks to age-matched young adult AGMs by intracerebroventricular (ICV) or intrathecal (IT) injections. Induction of AD-like pathology was assessed in subregions of the medial temporal lobe (MTL) by quantitative immunohistochemistry (IHC) using the AT8 antibody to detect hyperphosphorylated tau. Hippocampal volume was measured by magnetic resonance imaging (MRI) scans prior to, and after, intrathecal injections.

RESULTS

IT administration of AβOs in young adult AGMs revealed an elevation of tau phosphorylation in the MTL cortical memory circuit compared with controls. The largest increases were detected in the entorhinal cortex that persisted for at least 12 weeks after dosing. MRI scans showed a reduction in hippocampal volume following AβO injections.

CONCLUSIONS

Repeated IT delivery of AβOs in young adult AGMs led to an accelerated AD-like neuropathology in MTL, similar to human AD, supporting the value of this translational model to de-risk the clinical trial of diagnostic and therapeutic strategies.

Human in silico trials for parametric computational fluid dynamics investigation of cerebrospinal fluid drug delivery: impact of injection location, injection protocol, and physiology

Background

Intrathecal drug delivery has a significant role in pain management and central nervous system (CNS) disease therapeutics. A fluid-physics based tool to assist clinicians in choosing specific drug doses to the spine or brain may help improve treatment schedules.

Methods

This study applied computational fluid dynamics (CFD) and in vitro model verification to assess intrathecal drug delivery in an anatomically idealized model of the human CSF system with key anatomic features of the CNS. Key parameters analyzed included the role of (a) injection location including lumbar puncture (LP), cisterna magna (CM) and intracerebroventricular (ICV), (b) LP injection rate, injection volume, and flush volume, (c) physiologic factors including cardiac-induced and deep respiration-induced CSF stroke volume increase. Simulations were conducted for 3-h post-injection and used to quantify spatial–temporal tracer concentration, regional area under the curve (AUC), time to maximum concentration (T_max), and maximum concentration (C_max), for each case.

Results

CM and ICV increased AUC to brain regions by ~ 2 logs compared to all other simulations. A 3X increase in bolus volume and addition of a 5 mL flush both increased intracranial AUC to the brain up to 2X compared to a baseline 5 mL LP injection. In contrast, a 5X increase in bolus rate (25 mL/min) did not improve tracer exposure to the brain. An increase in cardiac and respiratory CSF movement improved tracer spread to the brain, basal cistern, and cerebellum up to ~ 2 logs compared to the baseline LP injection.

Conclusion

The computational modeling approach provides ability to conduct in silico trials representative of CSF injection protocols. Taken together, the findings indicate a strong potential for delivery protocols to be optimized to reach a target region(s) of the spine and/or brain with a needed therapeutic dose. Parametric modification of bolus rate/volume and flush volume was found to have impact on tracer distribution; albeit to a smaller degree than injection location, with CM and ICV injections resulting in greater therapeutic dose to brain regions compared to LP. CSF stroke volume and frequency both played an important role and may potentially have a greater impact than the modest changes in LP injection protocols analyzed such as bolus rate, volume, and flush.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12987-022-00304-4.

Translational approaches for brain delivery of biologics via cerebrospinal fluid

Delivery of biologics via cerebrospinal fluid (CSF) has demonstrated potential to access the tissues of the central nervous system (CNS) by circumventing the blood‐brain barrier and blood‐CSF barrier. Developing an effective CSF drug delivery strategy requires optimization of multiple parameters, including choice of CSF access point, delivery device technology, and delivery kinetics to achieve effective therapeutic concentrations in the target brain region, whereas also considering the biologic modality, mechanism of action, disease indication, and patient population. This review discusses key preclinical and clinical examples of CSF delivery for different biologic modalities (antibodies, nucleic acid‐based therapeutics, and gene therapy) to the brain via CSF or CNS access routes (intracerebroventricular, intrathecal‐cisterna magna, intrathecal‐lumbar, intraparenchymal, and intranasal), including the use of novel device technologies. This review also discusses quantitative models of CSF flow that provide insight into the effect of fluid dynamics in CSF on drug delivery and CNS distribution. Such models can facilitate delivery device design and pharmacokinetic/pharmacodynamic translation from preclinical species to humans in order to optimize CSF drug delivery to brain regions of interest.

Direction and magnitude of cerebrospinal fluid flow vary substantially across central nervous system diseases

Background

Several central nervous system diseases are associated with disturbed cerebrospinal fluid (CSF) flow patterns and have typically been characterized in vivo by phase-contrast magnetic resonance imaging (MRI). This technique is, however, limited by its applicability in space and time. Phase-contrast MRI has yet to be compared directly with CSF tracer enhanced imaging, which can be considered gold standard for assessing long-term CSF flow dynamics within the intracranial compartment.

Methods

Here, we studied patients with various CSF disorders and compared MRI biomarkers of CSF space anatomy and phase-contrast MRI at level of the aqueduct and cranio-cervical junction with dynamic intrathecal contrast-enhanced MRI using the contrast agent gadobutrol as CSF tracer. Tracer enrichment of cerebral ventricles was graded 0–4 by visual assessment. An intracranial pressure (ICP) score was used as surrogate marker of intracranial compliance.

Results

The study included 94 patients and disclosed marked variation of CSF flow measures across disease categories. The grade of supra-aqueductal reflux of tracer varied, with strong reflux (grades 3–4) in half of patients. Ventricular tracer reflux correlated with stroke volume and aqueductal CSF pressure gradient. CSF flow in the cerebral aqueduct was retrograde (from 4th to 3rd ventricle) in one third of patients, with estimated CSF net flow volume about 1.0 L/24 h. In the cranio-cervical junction, net flow was cranially directed in 78% patients, with estimated CSF net flow volume about 4.7 L/24 h.

Conclusions

The present observations provide in vivo quantitative evidence for substantial variation in direction and magnitude of CSF flow, with re-direction of aqueductal flow in communicating hydrocephalus, and significant extra-cranial CSF production. The grading of ventricular reflux of tracer shows promise as a clinical useful method to assess CSF flow pattern disturbances in patients.

Graphic abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s12987-021-00251-6.

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.

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.

Dispersion in porous media in oscillatory flow between flat plates: applications to intrathecal, periarterial and paraarterial solute transport in the central nervous system

Background

As an alternative to advection, solute transport by shear-augmented dispersion within oscillatory cerebrospinal fluid flow was investigated in small channels representing the basement membranes located between cerebral arterial smooth muscle cells, the paraarterial space surrounding the vessel wall and in large channels modeling the spinal subarachnoid space (SSS).

Methods

Geometries were modeled as two-dimensional. Fully developed flows in the channels were modeled by the Darcy–Brinkman momentum equation and dispersion by the passive transport equation. Scaling of the enhancement of axial dispersion relative to molecular diffusion was developed for regimes of flow including quasi-steady, porous and unsteady, and for regimes of dispersion including diffusive and unsteady.

Results

Maximum enhancement occurs when the characteristic time for lateral dispersion is matched to the cycle period. The Darcy–Brinkman model represents the porous media as a continuous flow resistance, and also imposes no-slip boundary conditions at the walls of the channel. Consequently, predicted dispersion is always reduced relative to that of a channel without porous media, except when the flow and dispersion are both unsteady.

Discussion/conclusions

In the basement membranes, flow and dispersion are both quasi-steady and enhancement of dispersion is small even if lateral dispersion is reduced by the porous media to achieve maximum enhancement. In the paraarterial space, maximum enhancement R_max = 73,200 has the potential to be significant. In the SSS, the dispersion is unsteady and the flow is in the transition zone between porous and unsteady. Enhancement is 5.8 times that of molecular diffusion, and grows to a maximum of 1.6E+6 when lateral dispersion is increased. The maximum enhancement produces rostral transport time in agreement with experiments.

Electronic supplementary material

The online version of this article (10.1186/s12987-019-0132-y) contains supplementary material, which is available to authorized users.

Comparison of phase-contrast MR and flow simulations for the study of CSF dynamics in the cervical spine

Background Investigators use phase-contrast magnetic resonance (PC-MR) and computational fluid dynamics (CFD) to assess cerebrospinal fluid dynamics. We compared qualitative and quantitative results from the two methods. Methods Four volunteers were imaged with a heavily T2-weighted volume gradient echo scan of the brain and cervical spine at 3T and with PC-MR. Velocities were calculated from PC-MR for each phase in the cardiac cycle. Mean pressure gradients in the PC-MR acquisition through the cardiac cycle were calculated with the Navier-Stokes equations. Volumetric MR images of the brain and upper spine were segmented and converted to meshes. Models of the subarachnoid space were created from volume images with the Vascular Modeling Toolkit. CFD simulations were performed with a previously verified flow solver. The flow patterns, velocities and pressures were compared in PC-MR and CFD flow images. Results PC-MR images consistently revealed more inhomogeneous flow patterns than CFD, especially in the anterolateral subarachnoid space where spinal nerve roots are located. On average, peak systolic and diastolic velocities in PC-MR exceeded those in CFD by 31% and 41%, respectively. On average, systolic and diastolic pressure gradients calculated from PC-MR exceeded those of CFD by 11% and 39%, respectively. Conclusions PC-MR shows local flow disturbances that are not evident in typical CFD. The velocities and pressure gradients calculated from PC-MR are systematically larger than those calculated from CFD.

An MRI-Compatible Hydrodynamic Simulator of Cerebrospinal Fluid Motion in the Cervical Spine

GOAL

Develop and test an MRI-compatible hydrodynamic simulator of cerebrospinal fluid (CSF) motion in the cervical spinal subarachnoid space. Four anatomically realistic subject-specific models were created based on a 22-year-old healthy volunteer and a five-year-old patient diagnosed with Chiari I malformation.

METHODS

The in vitro models were based on manual segmentation of high-resolution T2-weighted MRI of the cervical spine. Anatomically realistic dorsal and ventral spinal cord nerve rootlets (NR) were added. Models were three dimensional (3-D) printed by stereolithography with 50-μm layer thickness. A computer controlled pump system was used to replicate the shape of the subject specific in vivo CSF flow measured by phase-contrast MRI. Each model was then scanned by T2-weighted and 4-D phase contrast MRI (4D flow).

RESULTS

Cross-sectional area, wetted perimeter, and hydraulic diameter were quantified for each model. The oscillatory CSF velocity field (flow jets near NR, velocity profile shape, and magnitude) had similar characteristics to previously reported studies in the literature measured by in vivo MRI.

CONCLUSION

This study describes the first MRI-compatible hydrodynamic simulator of CSF motion in the cervical spine with anatomically realistic NR. NR were found to impact CSF velocity profiles to a great degree.

SIGNIFICANCE

CSF hydrodynamics are thought to be altered in craniospinal disorders such as Chiari I malformation. MRI scanning techniques and protocols can be used to quantify CSF flow alterations in disease states. The provided in vitro models can be used to test the reliability of these protocols across MRI scanner manufacturers and machines.