A model for transit time distributions through organs that accounts for fractal heterogeneity

Modeling across-trial variability in the Wald drift rate parameter

The shifted-Wald model is a popular analysis tool for one-choice reaction-time tasks. In its simplest version, the shifted-Wald model assumes a constant trial-independent drift rate parameter. However, the presence of endogenous processes—fluctuation in attention and motivation, fatigue and boredom—suggest that drift rate might vary across experimental trials. Here we show how across-trial variability in drift rate can be accounted for by assuming a trial-specific drift rate parameter that is governed by a positive-valued distribution. We consider two candidate distributions: the truncated normal distribution and the gamma distribution. For the resulting distributions of first-arrival times, we derive analytical and sampling-based solutions, and implement the models in a Bayesian framework. Recovery studies and an application to a data set comprised of 1469 participants suggest that (1) both mixture distributions yield similar results; (2) all model parameters can be recovered accurately except for the drift variance parameter; (3) despite poor recovery, the presence of the drift variance parameter facilitates accurate recovery of the remaining parameters; (4) shift, threshold, and drift mean parameters are correlated.

Usefulness of continuous probability distributions of rates for modelling radionuclide biokinetics in humans and animals

Modelling the biokinetics of radionuclide excretion or retention is important in nuclear medicine and following accidental/malicious radioactivity releases. Sums of discrete exponential decay rates are often used, but we hypothesized that continuous probability distributions (CPD) of decay rates can describe the data more parsimoniously and robustly. We tested this hypothesis on diverse human and animal data sets involving various radionuclides (including plutonium, strontium, caesium) measured in the laboratory and in regions contaminated by the Fukushima and Chernobyl nuclear accidents. We used four models on each data set: mono-exponential (ME) with one discrete decay rate, bi-exponential (BE) with two rates, gamma-exponential (GE) with a Gamma distribution of stretched-exponential rates, and power-decay (PD) with a Gamma distribution of power-decay rates. Information-theoretic model selection suggested that radionuclide biokinetics, e.g. for plutonium in humans, are often better described by CPD models like GE and PD, than by discrete rates (ME and BE). Extrapolation of models fitted to data at short times to longer times was frequently more robust for CPD formalisms. We suggest that using a set of several CPD and discrete-rate models, and comparing them by information-theoretic methods, is a promising strategy to enhance the analysis of radionuclide excretion and retention kinetics.

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.