Understanding pharmacokinetics using realistic computational models of fluid dynamics: biosimulation of drug distribution within the CSF space for intrathecal drugs

We introduce how biophysical modeling in pharmaceutical research and development, combining physiological observations at the tissue, organ and system level with selected drug physiochemical properties, may contribute to a greater and non-intuitive understanding of drug pharmacokinetics and therapeutic design. Based on rich first-principle knowledge combined with experimental data at both conception and calibration stages, and leveraging our insights on disease processes and drug pharmacology, biophysical modeling may provide a novel and unique opportunity to interactively characterize detailed drug transport, distribution, and subsequent therapeutic effects. This innovative approach is exemplified through a three-dimensional (3D) computational fluid dynamics model of the spinal canal motivated by questions arising during pharmaceutical development of one molecular therapy for spinal cord injury. The model was based on actual geometry reconstructed from magnetic resonance imaging data subsequently transformed in a parametric 3D geometry and a corresponding finite-volume representation. With dynamics controlled by transient Navier–Stokes equations, the model was implemented in a commercial multi-physics software environment established in the automotive and aerospace industries. While predictions were performed in silico, the underlying biophysical models relied on multiple sources of experimental data and knowledge from scientific literature. The results have provided insights into the primary factors that can influence the intrathecal distribution of drug after lumbar administration. This example illustrates how the approach connects the causal chain underlying drug distribution, starting with the technical aspect of drug delivery systems, through physiology-driven drug transport, then eventually linking to tissue penetration, binding, residence, and ultimately clearance. Currently supporting our drug development projects with an improved understanding of systems physiology, biophysical models are being increasingly used to characterize drug transport and distribution in human tissues where pharmacokinetic measurements are difficult or impossible to perform. Importantly, biophysical models can describe emergent properties of a system, i.e. properties not identifiable through the study of the system’s components taken in isolation.

[1]  U Klose,et al.  The Influence of Pulse and Respiration on Spinal Cerebrospinal Fluid Pulsation , 2004, Investigative radiology.

[2]  M. Seeberger,et al.  Bupivacaine Concentrations in the Lumbar Cerebrospinal Fluid of Patients During Spinal Anesthesia , 2010 .

[3]  J. Wildsmith,et al.  Failed Spinal Anesthesia: Mechanisms, Management, and Prevention , 2010 .

[4]  Anusha Sivaramakrishnan,et al.  Magnetic resonance imaging-based measurements of cerebrospinal fluid and blood flow as indicators of intracranial compliance in patients with Chiari malformation. , 2005, Journal of neurosurgery.

[5]  H. Stockman,et al.  Effect of anatomical fine structure on the flow of cerebrospinal fluid in the spinal subarachnoid space. , 2006, Journal of biomechanical engineering.

[6]  K. Drasner Spinal Anesthesia: A Century of Refinement, and Failure Is Still an Option , 2010 .

[7]  N. Alperin,et al.  MR-Intracranial pressure (ICP): a method to measure intracranial elastance and pressure noninvasively by means of MR imaging: baboon and human study. , 2000, Radiology.

[8]  M. Schwab,et al.  Inhibition of Nogo: A key strategy to increase regeneration, plasticity and functional recovery of the lesioned central nervous system , 2005, Annals of medicine.

[9]  P. White,et al.  Neural Blockade: In Clinical Anesthesia and Pain Medicine, 4th ed. , 2010 .

[10]  K. Drasner Spinal anaesthesia: a century of refinement, and failure is still an option. , 2009, British journal of anaesthesia.

[11]  I. Russell Neural Blockade in Clinical Anesthesia and Pain Medicine , 2010 .

[12]  J. Wildsmith,et al.  Failed spinal anaesthesia: mechanisms, management, and prevention. , 2009, British journal of anaesthesia.

[13]  M. Seeberger,et al.  Bupivacaine concentrations in the lumbar cerebrospinal fluid of patients during spinal anaesthesia. , 2009, British journal of anaesthesia.

[14]  S. McJames,et al.  The Effect of Prior Dural Puncture on Cerebrospinal Fluid Sufentanil Concentrations in Sheep After the Administration of Epidural Sufentanil , 1998, Anesthesia and analgesia.

[15]  R K Jain,et al.  Transport of molecules in the tumor interstitium: a review. , 1987, Cancer research.

[16]  A. Møller,et al.  [Spinal and epidural anesthesia]. , 2007, Ugeskrift for laeger.

[17]  T. Takarada,et al.  Riluzole, a Glutamate Release Inhibitor, Induces Loss of Righting Reflex, Antinociception, and Immobility in Response to Noxious Stimulation in Mice , 2007, Anesthesia and analgesia.

[18]  D. Gaver,et al.  An experimental investigation of oscillating flow in a tapered channel , 1986, Journal of Fluid Mechanics.

[19]  C. Bernards Cerebrospinal Fluid and Spinal Cord Distribution of Baclofen and Bupivacaine during Slow Intrathecal Infusion in Pigs , 2006, Anesthesiology.

[20]  D J Crouch,et al.  Effects of intrathecal morphine on the ventilatory response to hypoxia. , 2000, The New England journal of medicine.

[21]  B. Launder,et al.  The numerical computation of turbulent flows , 1990 .

[22]  B. Mokri The Monro–Kellie hypothesis , 2001, Neurology.

[23]  R. R. Lee,et al.  Dynamic Physiologic Changes in Lumbar CSF Volume Quantitatively Measured By Three-Dimensional Fast Spin-Echo MRI , 2001, Spine.

[24]  G. di Chiro,et al.  Spinal descent of cerebrospinal fluid in man , 1976, Neurology.

[25]  N Alperin,et al.  Hydrodynamic modeling of cerebrospinal fluid motion within the spinal cavity. , 2000, Journal of biomechanical engineering.

[26]  Steven Edward Kern,et al.  The Pharmacokinetics of the Conopeptide Contulakin-G (CGX-1160) After Intrathecal Administration: An Analysis of Data from Studies in Beagles , 2007, Anesthesia and Analgesia.

[27]  M. Seeberger,et al.  Bupivacaine concentrations in lumbar cerebrospinal fluid in patients with failed spinal anaesthesia. , 2009, British journal of anaesthesia.

[28]  A. Martins,et al.  Dynamics of the cerebrospinal fluid and the spinal dura mater1 , 1972, Journal of neurology, neurosurgery, and psychiatry.

[29]  Fred H. Gage,et al.  Therapeutic interventions after spinal cord injury , 2006, Nature Reviews Neuroscience.

[30]  G R Harrison,et al.  The topographical anatomy of the lumbar epidural space. , 1985, Journal of anatomy.