Spatio-Temporal Simulation of First Pass Drug Perfusion in the Liver

The liver is the central organ for detoxification of xenobiotics in the body. In pharmacokinetic modeling, hepatic metabolization capacity is typically quantified as hepatic clearance computed as degradation in well-stirred compartments. This is an accurate mechanistic description once a quasi-equilibrium between blood and surrounding tissue is established. However, this model structure cannot be used to simulate spatio-temporal distribution during the first instants after drug injection. In this paper, we introduce a new spatially resolved model to simulate first pass perfusion of compounds within the naive liver. The model is based on vascular structures obtained from computed tomography as well as physiologically based mass transfer descriptions obtained from pharmacokinetic modeling. The physiological architecture of hepatic tissue in our model is governed by both vascular geometry and the composition of the connecting hepatic tissue. In particular, we here consider locally distributed mass flow in liver tissue instead of considering well-stirred compartments. Experimentally, the model structure corresponds to an isolated perfused liver and provides an ideal platform to address first pass effects and questions of hepatic heterogeneity. The model was evaluated for three exemplary compounds covering key aspects of perfusion, distribution and metabolization within the liver. As pathophysiological states we considered the influence of steatosis and carbon tetrachloride-induced liver necrosis on total hepatic distribution and metabolic capacity. Notably, we found that our computational predictions are in qualitative agreement with previously published experimental data. The simulation results provide an unprecedented level of detail in compound concentration profiles during first pass perfusion, both spatio-temporally in liver tissue itself and temporally in the outflowing blood. We expect our model to be the foundation of further spatially resolved models of the liver in the future.

[1]  Gerhard Rakhorst,et al.  Hypothermic Machine Preservation in Liver Transplantation Revisited: Concepts and Criteria in the New Millennium , 2004, Annals of Biomedical Engineering.

[2]  B O Palsson,et al.  Mathematical modelling of dynamics and control in metabolic networks. I. On Michaelis-Menten kinetics. , 1984, Journal of theoretical biology.

[3]  Stephen G. Monismith,et al.  Evaluation of Advective Schemes for Estuarine Salinity Simulations , 2000 .

[4]  Michael S. Roberts,et al.  A dispersion model of hepatic elimination: 1. Formulation of the model and bolus considerations , 1986, Journal of Pharmacokinetics and Biopharmaceutics.

[5]  M. Weiss A model for transit time distributions through organs that accounts for fractal heterogeneity. , 2012, Journal of theoretical biology.

[6]  R. Courant,et al.  Über die partiellen Differenzengleichungen der mathematischen Physik , 1928 .

[7]  R. Chabiniok,et al.  A novel porous mechanical framework for modelling the interaction between coronary perfusion and myocardial mechanics , 2012, Journal of biomechanics.

[8]  S. Narimatsu,et al.  Separate evaluation of intestinal and hepatic metabolism of three benzodiazepines in rats with cannulated portal and jugular veins: comparison with the profile in non-cannulated mice , 2009, Xenobiotica; the fate of foreign compounds in biological systems.

[9]  Milton Abramowitz,et al.  Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables , 1964 .

[10]  Horst K. Hahn,et al.  Reconstruction of vessel structures from serial whole slide sections of murine liver samples , 2013, Medical Imaging.

[11]  Thomas J. Raub,et al.  Roles of P-Glycoprotein, Bcrp, and Mrp2 in Biliary Excretion of Spiramycin in Mice , 2007, Antimicrobial Agents and Chemotherapy.

[12]  Alexis B. Carter Review of methods in medical informatics: Fundamentals of healthcare programming in Perl, Python and Ruby by Jules J. Berman , 2011, Journal of Pathology Informatics.

[13]  M. Roberts,et al.  Metabolite Mean Transit Times in the Liver as Predicted By Various Models of Hepatic Elimination , 1997, Journal of Pharmacokinetics and Biopharmaceutics.

[14]  Malcolm Rowland,et al.  Hepatic clearance of drugs. I. Theoretical considerations of a “well-stirred” model and a “parallel tube” model. Influence of hepatic blood flow, plasma and blood cell binding, and the hepatocellular enzymatic activity on hepatic drug clearance , 1977, Journal of Pharmacokinetics and Biopharmaceutics.

[15]  Hannah M. Jones,et al.  Modelling and PBPK Simulation in Drug Discovery , 2009, The AAPS Journal.

[16]  C. Nickerson A note on a concordance correlation coefficient to evaluate reproducibility , 1997 .

[17]  John Wambaugh,et al.  Simulating Microdosimetry in a Virtual Hepatic Lobule , 2010, PLoS Comput. Biol..

[18]  T. F. Russell,et al.  An Eulerian-Lagrangian localized adjoint method for the advection-diffusion equation , 1990 .

[19]  H. Teutsch The modular microarchitecture of human liver , 2005, Hepatology.

[20]  Gjlles Aubert,et al.  Mathematical problems in image processing , 2001 .

[21]  A. Henney,et al.  The virtual liver: a multidisciplinary, multilevel challenge for systems biology , 2012, Wiley interdisciplinary reviews. Systems biology and medicine.

[22]  Aijie Cheng,et al.  A preliminary study on multiscale ELLAM schemes for transient advection‐diffusion equations , 2010 .

[23]  R. Siegel Pharmacokinetic transfer functions and generalized clearances , 1986, Journal of Pharmacokinetics and Biopharmaceutics.

[24]  Takuji Nishimura,et al.  Mersenne twister: a 623-dimensionally equidistributed uniform pseudo-random number generator , 1998, TOMC.

[25]  Mattia Quattrocelli,et al.  Development of a New Tool for 3D Modeling for Regenerative Medicine , 2011, Int. J. Biomed. Imaging.

[26]  T. Minor,et al.  The isolated perfused rat liver: standardization of a time-honoured model , 2006, Laboratory animals.

[27]  David L. Woods,et al.  Targeted drug delivery by high intensity focused ultrasound mediated hyperthermia combined with temperature-sensitive liposomes : Computational modelling and preliminary in vivo validation , 2012 .

[28]  Jörg Lippert,et al.  From physicochemistry to absorption and distribution: predictive mechanistic modelling and computational tools , 2005, Expert opinion on drug metabolism & toxicology.

[29]  L. Kuepfer Towards whole-body systems physiology , 2010, Molecular systems biology.

[30]  Klaas Nicolay,et al.  Block-copolymer-stabilized iodinated emulsions for use as CT contrast agents. , 2010, Biomaterials.

[31]  Horst Karl Hahn,et al.  Registration of histological whole slide images guided by vessel structures , 2013, Journal of pathology informatics.

[32]  A. Lonardo,et al.  Review article: hepatic steatosis and insulin resistance , 2005, Alimentary pharmacology & therapeutics.

[33]  J. Oberholzer,et al.  Pancreatic islet cell transplantation: an update for interventional radiologists. , 2012, Journal of vascular and interventional radiology : JVIR.

[34]  L. Kuepfer,et al.  Using Expression Data for Quantification of Active Processes in Physiologically Based Pharmacokinetic Modeling , 2012, Drug Metabolism and Disposition.

[35]  Rebecca Fahrig,et al.  Deconvolution-Based CT and MR Brain Perfusion Measurement: Theoretical Model Revisited and Practical Implementation Details , 2011, Int. J. Biomed. Imaging.

[36]  C. Buechler,et al.  Does hepatic steatosis affect drug metabolizing enzymes in the liver? , 2011, Current drug metabolism.

[37]  Paolo Crosetto,et al.  Implicit Coupling of One-Dimensional and Three-Dimensional Blood Flow Models with Compliant Vessels , 2013, Multiscale Model. Simul..

[38]  Markus P. Schlaich,et al.  Change in Sympathetic Nerve Firing Pattern Associated with Dietary Weight Loss in the Metabolic Syndrome , 2011, Front. Physio..

[39]  J. Calixto,et al.  ABCC transporter inhibition reduces zymosan‐induced peritonitis , 2007, Journal of leukocyte biology.

[40]  Nan Zheng,et al.  A Cell-based Computational Modeling Approach for Developing Site-Directed Molecular Probes , 2012, PLoS Comput. Biol..

[41]  Erwin Fehlberg,et al.  Klassische Runge-Kutta-Formeln vierter und niedrigerer Ordnung mit Schrittweiten-Kontrolle und ihre Anwendung auf Wärmeleitungsprobleme , 1970, Computing.

[42]  Hong Wang,et al.  Runge–Kutta characteristic methods for first‐order linear hyperbolic equations , 1997 .

[43]  M. Weiss,et al.  Advanced pharmacokinetic models based on organ clearance, circulatory, and fractal concepts , 2007, The AAPS Journal.

[44]  M. Weiss,et al.  Pharmacokinetic model based on circulatory transport , 1979, European Journal of Clinical Pharmacology.

[45]  T. Chang,et al.  Numerical investigation of non-Newtonian microcirculatory blood flow in hepatic lobule. , 2006, Journal of biomechanics.

[46]  J. Scoazec,et al.  Vascular Development and Differentiation During Human Liver Organogenesis , 2008, Anatomical record.

[47]  F. Hossler,et al.  Vascular Corrosion Casting: Review of Advantages and Limitations in the Application of Some Simple Quantitative Methods , 2001, Microscopy and Microanalysis.

[48]  J. Lippert,et al.  Integration of dissolution into physiologically‐based pharmacokinetic models III: PK‐Sim® , 2012, The Journal of pharmacy and pharmacology.

[49]  Thomas J. R. Hughes,et al.  A stabilized mixed finite element method for Darcy flow , 2002 .

[50]  Fabian Kiessling,et al.  Virtual elastic sphere processing enables reproducible quantification of vessel stenosis at CT and MR angiography. , 2011, Radiology.

[51]  A. Seifalian,et al.  Impairment of Hepatic Microcirculation in Fatty Liver , 2003, Microcirculation.

[52]  T. Sochi,et al.  A computationally efficient framework for the simulation of cardiac perfusion using a multi‐compartment Darcy porous‐media flow model , 2013, International journal for numerical methods in biomedical engineering.

[53]  G R Wilkinson,et al.  Clearance approaches in pharmacology. , 1987, Pharmacological reviews.

[54]  George Loizou,et al.  Development of good modelling practice for physiologically based pharmacokinetic models for use in risk assessment: the first steps. , 2008, Regulatory toxicology and pharmacology : RTP.

[55]  Richard E. Ewing,et al.  A Family of Eulerian-Lagrangian Localized Adjoint Methods for Multi-dimensional Advection-Reaction Equations , 1999 .

[56]  M. Roberts,et al.  Hepatic Pharmacokinetics of Cationic Drugs in a High-Fat Emulsion-Induced Rat Model of Nonalcoholic Steatohepatitis , 2011, Drug Metabolism and Disposition.

[57]  Y. Sai,et al.  Prediction of theophylline clearance in CCl4-treated rats using in vivo CYP1A2 and CYP3A2 contents assessed with the PKCYP test. , 2005, Drug metabolism and pharmacokinetics.

[58]  Peter Bastian,et al.  On a monotonicity preserving Eulerian–Lagrangian localized adjoint method for advection–diffusion equations , 2005 .

[59]  Robert Roskoski Michaelis-Menten Kinetics , 2007 .

[60]  B. Luxon,et al.  Hepatic transport kinetics and plasma disappearance curves: distributed modeling versus conventional approach. , 1978, The American journal of physiology.

[61]  D. Cutler A linear recirculation model for drug disposition , 1979, Journal of Pharmacokinetics and Biopharmaceutics.

[62]  G. Brandi,et al.  Prevalence of and Risk Factors for Hepatic Steatosis in Northern Italy , 2000, Annals of Internal Medicine.

[63]  H. Elias,et al.  Observations on the anatomy of the rat liver , 1954, The Anatomical record.

[64]  Richard E. Ewing,et al.  Numerical treatment of multiphase flows in porous media : proceedings of the international workshop, held at Beijing, China, 2-6 August 1999 , 2000 .

[65]  Joachim Georgii,et al.  Workshop on Virtual Reality Interaction and Physical Simulation Vriphys (2011) Focused Ultrasound - Efficient Gpu Simulation Methods for Therapy Planning , 2022 .

[66]  J. Hampe,et al.  A Mechanistic, Model-Based Approach to Safety Assessment in Clinical Development , 2012, CPT: pharmacometrics & systems pharmacology.

[67]  M. Jamei,et al.  PBPK modelling of inter-individual variability in the pharmacokinetics of environmental chemicals. , 2010, Toxicology.

[68]  L. Weber,et al.  Hepatotoxicity and Mechanism of Action of Haloalkanes: Carbon Tetrachloride as a Toxicological Model , 2003, Critical reviews in toxicology.

[69]  K. Kichikawa,et al.  Radiologic placement of side-hole catheter with tip fixation for hepatic arterial infusion chemotherapy. , 2003, Journal of vascular and interventional radiology : JVIR.

[70]  Heinz-Otto Peitgen,et al.  Clinical relevance of model based computer-assisted diagnosis and therapy , 2008, SPIE Medical Imaging.

[71]  Jing-yu Yu,et al.  Cell-Based Multiscale Computational Modeling of Small Molecule Absorption and Retention in the Lungs , 2010, Pharmaceutical Research.

[72]  Steffen Borchers,et al.  Integrating Cellular Metabolism into a Multiscale Whole-Body Model , 2012, PLoS Comput. Biol..

[73]  O. Dirsch,et al.  A biphasic model for sinusoidal liver perfusion remodeling after outflow obstruction , 2010, Biomechanics and modeling in mechanobiology.

[74]  P. Georgopoulos,et al.  Tissue‐Level Modeling of Xenobiotic Metabolism in Liver: An Emerging Tool for Enabling Clinical Translational Research , 2009, Clinical and translational science.

[75]  L. Lin,et al.  A concordance correlation coefficient to evaluate reproducibility. , 1989, Biometrics.

[76]  Bernhard Preim,et al.  Analysis of vasculature for liver surgical planning , 2002, IEEE Transactions on Medical Imaging.

[77]  M. Labarbera Principles of design of fluid transport systems in zoology. , 1990, Science.

[78]  Abu-Zahra Tn,et al.  EFFECT OF ZONAL TRANSPORT AND METABOLISM ON HEPATIC REMOVAL: ENALAPRIL HYDROLYSIS IN ZONAL, ISOLATED RAT HEPATOCYTES IN VITRO AND CORRELATION WITH PERFUSION DATA , 2000 .

[79]  Hong Wang,et al.  A summary of numerical methods for time-dependent advection-dominated partial differential equations , 2001 .

[80]  Tatsuhisa Takahashi,et al.  Quantitative assessments of morphological and functional properties of biological trees based on their fractal nature. , 2007, Journal of applied physiology.

[81]  Stefan Willmann,et al.  Physiology-Based Simulations of a Pathological Condition , 2008, Clinical pharmacokinetics.

[82]  Hong Wang,et al.  Eulerian-Lagrangian localized adjoint methods for systems of nonlinear advective-diffusive-reactive transport equations , 1996 .

[83]  J. Mak,et al.  Significant Reduction of Brain Cysts Caused by Toxoplasma gondii after Treatment with Spiramycin Coadministered with Metronidazole in a Mouse Model of Chronic Toxoplasmosis , 2012, Antimicrobial Agents and Chemotherapy.

[84]  T. Vree,et al.  Pharmacokinetics of spiramycin in the rhesus monkey: transplacental passage and distribution in tissue in the fetus , 1994, Antimicrobial Agents and Chemotherapy.

[85]  Dieter Haemmerich,et al.  Targeted drug delivery by high intensity focused ultrasound mediated hyperthermia combined with temperature-sensitive liposomes: Computational modelling and preliminary in vivovalidation , 2012, International journal of hyperthermia : the official journal of European Society for Hyperthermic Oncology, North American Hyperthermia Group.

[86]  Wolfgang Weiss,et al.  A Computational Systems Biology Software Platform for Multiscale Modeling and Simulation: Integrating Whole-Body Physiology, Disease Biology, and Molecular Reaction Networks , 2011, Front. Physio..

[87]  W W Lautt,et al.  The 1995 Ciba-Geigy Award Lecture. Intrinsic regulation of hepatic blood flow. , 1996, Canadian journal of physiology and pharmacology.

[88]  R. Koehler,et al.  Blood viscosity modulates tissue perfusion: sometimes and somewhere. , 2007, Transfusion alternatives in transfusion medicine : TATM.

[89]  Panos Macheras,et al.  A Model for Transport and Dispersion in the Circulatory System Based on the Vascular Fractal Tree , 2003, Annals of Biomedical Engineering.

[90]  Horst K. Hahn,et al.  A comparison of sampling strategies for histological image analysis , 2011, Journal of pathology informatics.

[91]  R. Gebhardt Metabolic zonation of the liver: regulation and implications for liver function. , 1992, Pharmacology & therapeutics.

[92]  W Schreiner,et al.  Shear stress distribution in arterial tree models, generated by constrained constructive optimization. , 1999, Journal of theoretical biology.

[93]  N. Furusawa Binding profile of spiramycin to oviducal proteins of laying hens. , 2000, Journal of veterinary medicine. A, Physiology, pathology, clinical medicine.

[94]  Matthias Hermes,et al.  Prediction and validation of cell alignment along microvessels as order principle to restore tissue architecture in liver regeneration , 2010, Proceedings of the National Academy of Sciences.

[95]  K. Pang,et al.  Effect of zonal transport and metabolism on hepatic removal: enalapril hydrolysis in zonal, isolated rat hepatocytes in vitro and correlation with perfusion data. , 2000, Drug metabolism and disposition: the biological fate of chemicals.

[96]  Tim Morris,et al.  Physiological Parameters in Laboratory Animals and Humans , 1993, Pharmaceutical Research.

[97]  V. Gebski,et al.  Randomised trial of SIR-Spheres plus chemotherapy vs. chemotherapy alone for treating patients with liver metastases from primary large bowel cancer. , 2001, Annals of oncology : official journal of the European Society for Medical Oncology.

[98]  Tobias Preusser,et al.  Analysis and Algorithmic Generation of Hepatic Vascular Systems , 2012, International journal of hepatology.

[99]  N. Michoux,et al.  Transvascular and interstitial transport in rat hepatocellular carcinomas: Dynamic contrast‐enhanced MRI assessment with low‐ and high‐molecular weight agents , 2008, Journal of magnetic resonance imaging : JMRI.

[100]  T. F. Russell,et al.  An overview of research on Eulerian-Lagrangian localized adjoint methods (ELLAM). , 2002 .

[101]  D. Madoff,et al.  Liver Anatomy: Microcirculation of the Liver , 2008, Seminars in interventional radiology.

[102]  Philippe Ackerer,et al.  Solving the advection-diffusion equation with the Eulerian-Lagrangian localized adjoint method on unstructured meshes and non uniform time stepping , 2005 .