Performance Assessment and Translation of Physiologically Based Pharmacokinetic Models From acslX to Berkeley Madonna, MATLAB, and R Language: Oxytetracycline and Gold Nanoparticles As Case Examples.

Many physiologically based pharmacokinetic (PBPK) models for environmental chemicals, drugs, and nanomaterials have been developed to aid risk and safety assessments using acslX. However, acslX has been rendered sunset since November 2015. Alternative modeling tools and tutorials are needed for future PBPK applications. This forum article aimed to: (1) demonstrate the performance of 4 PBPK modeling software packages (acslX, Berkeley Madonna, MATLAB, and R language) tested using 2 existing models (oxytetracycline and gold nanoparticles); (2) provide a tutorial of PBPK model code conversion from acslX to Berkeley Madonna, MATLAB, and R language; (3) discuss the advantages and disadvantages of each software package in the implementation of PBPK models in toxicology, and (4) share our perspective about future direction in this field. Simulation results of plasma/tissue concentrations/amounts of oxytetracycline and gold from different models were compared visually and statistically with linear regression analyses. Simulation results from the original models were correlated well with results from the recoded models, with time-concentration/amount curves nearly superimposable and determination coefficients of 0.86-1.00. Step-by-step explanations of the recoding of the models in different software programs are provided in the Supplementary Data. In summary, this article presents a tutorial of PBPK model code conversion for a small molecule and a nanoparticle among 4 software packages, and a performance comparison of these software packages in PBPK model implementation. This tutorial helps beginners learn PBPK modeling, provides suggestions for selecting a suitable tool for future projects, and may lead to the transition from acslX to alternative modeling tools.

Improving the risk assessment of lipophilic persistent environmental chemicals in breast milk.

Lipophilic persistent environmental chemicals (LPECs) have the potential to accumulate within a woman's body lipids over the course of many years prior to pregnancy, to partition into human milk, and to transfer to infants upon breastfeeding. As a result of this accumulation and partitioning, a breastfeeding infant's intake of these LPECs may be much greater than his/her mother's average daily exposure. Because the developmental period sets the stage for lifelong health, it is important to be able to accurately assess chemical exposures in early life. In many cases, current human health risk assessment methods do not account for differences between maternal and infant exposures to LPECs or for lifestage-specific effects of exposure to these chemicals. Because of their persistence and accumulation in body lipids and partitioning into breast milk, LPECs present unique challenges for each component of the human health risk assessment process, including hazard identification, dose-response assessment, and exposure assessment. Specific biological modeling approaches are available to support both dose-response and exposure assessment for lactational exposures to LPECs. Yet, lack of data limits the application of these approaches. The goal of this review is to outline the available approaches and to identify key issues that, if addressed, could improve efforts to apply these approaches to risk assessment of lactational exposure to these chemicals.

Lifetime PCB 153 bioaccumulation and pharmacokinetics in pilot whales: Bayesian population PBPK modeling and Markov chain Monte Carlo simulations.

Physiologically based pharmacokinetic (PBPK) models for wild animal populations such as marine mammals typically have a high degree of model uncertainty and variability due to the scarcity of information and the embryonic nature of this field. Parameters values used in marine mammals models are usually taken from other mammalian species (e.g. rats or mice) and might not be entirely suitable to properly explain the kinetics of pollutants in marine mammals. Therefore, several parameters for a PBPK model for the bioaccumulation and pharmacokinetics of PCB 153 in long-finned pilot whales were estimated in the present study using the Bayesian approach executed with Markov chain Monte Carlo (MCMC) simulations. This method uses 'prior' information of the parameters, either from the literature or from previous model runs. The advantage is that this method uses such 'prior' parameters to calculate probability distributions to determine 'posterior' values that best explain the field observations. Those field observations or datasets were PCB 153 concentrations in blubber of long-finned pilot whales from Sandy Cape and Stanley, Tasmania, Australia. The model predictions showed an overall decrease in PCB 153 levels in blubber over the lifetime of the pilot whales. All parameters from the Sandy Cape model were updated using the Stanley dataset, except for the concentration of PCB 153 in the milk. The model presented here is a promising and preliminary start to PBPK modeling in long-finned pilot whales that would provide a basis for non-invasive studies in these protected marine mammals.

Application of Bayesian population physiologically based pharmacokinetic (PBPK) modeling and Markov chain Monte Carlo simulations to pesticide kinetics studies in protected marine mammals: DDT, DDE, and DDD in harbor porpoises.

Physiologically based pharmacokinetic (PBPK) modeling in marine mammals is a challenge because of the lack of parameter information and the ban on exposure experiments. To minimize uncertainty and variability, parameter estimation methods are required for the development of reliable PBPK models. The present study is the first to develop PBPK models for the lifetime bioaccumulation of p,p'-DDT, p,p'-DDE, and p,p'-DDD in harbor porpoises. In addition, this study is also the first to apply the Bayesian approach executed with Markov chain Monte Carlo simulations using two data sets of harbor porpoises from the Black and North Seas. Parameters from the literature were used as priors for the first "model update" using the Black Sea data set, the resulting posterior parameters were then used as priors for the second "model update" using the North Sea data set. As such, PBPK models with parameters specific for harbor porpoises could be strengthened with more robust probability distributions. As the science and biomonitoring effort progress in this area, more data sets will become available to further strengthen and update the parameters in the PBPK models for harbor porpoises as a species anywhere in the world. Further, such an approach could very well be extended to other protected marine mammals.

A physiologically based pharmacokinetic model of rifampin in mice.

One problem associated with regimen-based development of antituberculosis (anti-TB) drugs is the difficulty of a systematic and thorough in vivo evaluation of the large number of possible regimens that arise from consideration of multiple drugs tested together. A mathematical model capable of simulating the pharmacokinetics and pharmacodynamics of experimental combination chemotherapy of TB offers a way to mitigate this problem by extending the use of available data to investigate regimens that are not initially tested. In order to increase the available mathematical tools needed to support such a model for preclinical anti-TB drug development, we constructed a preliminary whole-body physiologically based pharmacokinetic (PBPK) model of rifampin in mice, using data from the literature. Interindividual variability was approximated using Monte Carlo (MC) simulation with assigned probability distributions for the model parameters. An MC sensitivity analysis was also performed to determine correlations between model parameters and plasma concentration to inform future model development. Model predictions for rifampin concentrations in plasma, liver, kidneys, and lungs, following oral administration, were generally in agreement with published experimental data from multiple studies. Sensitive model parameters included those descriptive of oral absorption, total clearance, and partitioning of rifampin between blood and muscle. This PBPK model can serve as a starting point for the integration of rifampin pharmacokinetics in mice into a larger mathematical framework, including the immune response to Mycobacterium tuberculosis infection, and pharmacokinetic models for other anti-TB drugs.

An in silico approach for evaluating a fraction-based, risk assessment method for total petroleum hydrocarbon mixtures.

Both the Massachusetts Department of Environmental Protection (MADEP) and the Total Petroleum Hydrocarbon Criteria Working Group (TPHCWG) developed fraction-based approaches for assessing human health risks posed by total petroleum hydrocarbon (TPH) mixtures in the environment. Both organizations defined TPH fractions based on their expected environmental fate and by analytical chemical methods. They derived toxicity values for selected compounds within each fraction and used these as surrogates to assess hazard or risk of exposure to the whole fractions. Membership in a TPH fraction is generally defined by the number of carbon atoms in a compound and by a compound's equivalent carbon (EC) number index, which can predict its environmental fate. Here, we systematically and objectively re-evaluate the assignment of TPH to specific fractions using comparative molecular field analysis and hierarchical clustering. The approach is transparent and reproducible, reducing inherent reliance on judgment when toxicity information is limited. Our evaluation of membership in these fractions is highly consistent (˜80% on average across various fractions) with the empirical approach of MADEP and TPHCWG. Furthermore, the results support the general methodology of mixture risk assessment to assess both cancer and noncancer risk values after the application of fractionation.

Computational toxicology of chloroform: reverse dosimetry using Bayesian inference, Markov chain Monte Carlo simulation, and human biomonitoring data.

BACKGROUND: One problem of interpreting population-based biomonitoring data is the reconstruction of corresponding external exposure in cases where no such data are available. OBJECTIVES: We demonstrate the use of a computational framework that integrates physiologically based pharmacokinetic (PBPK) modeling, Bayesian inference, and Markov chain Monte Carlo simulation to obtain a population estimate of environmental chloroform source concentrations consistent with human biomonitoring data. The biomonitoring data consist of chloroform blood concentrations measured as part of the Third National Health and Nutrition Examination Survey (NHANES III), and for which no corresponding exposure data were collected. METHODS: We used a combined PBPK and shower exposure model to consider several routes and sources of exposure: ingestion of tap water, inhalation of ambient household air, and inhalation and dermal absorption while showering. We determined posterior distributions for chloroform concentration in tap water and ambient household air using U.S. Environmental Protection Agency Total Exposure Assessment Methodology (TEAM) data as prior distributions for the Bayesian analysis. RESULTS: Posterior distributions for exposure indicate that 95% of the population represented by the NHANES III data had likely chloroform exposures < or = 67 microg/L [corrected] in tap water and < or = 0.02 microg/L in ambient household air. CONCLUSIONS: Our results demonstrate the application of computer simulation to aid in the interpretation of human biomonitoring data in the context of the exposure-health evaluation-risk assessment continuum. These results should be considered as a demonstration of the method and can be improved with the addition of more detailed data.

Persistent tissue kinetics and redistribution of nanoparticles, quantum dot 705, in mice: ICP-MS quantitative assessment.

BACKGROUND: Quantum dots (QDs) are autofluorescent semiconductor nanocrystals that can be used for in vivo biomedical imaging. However, we know little about their in vivo disposition and health consequences. OBJECTIVES: We assessed the tissue disposition and pharmacokinetics of QD705 in mice. METHODS: We determined quantitatively the blood and tissue kinetics of QD705 in mice after single intravenous (iv) injection at the dose of 40 pmol for up to 28 days. Inductively coupled plasma-mass spectrometry (ICP-MS) measurement of cadmium was the primary method of quantification of QD705. Fluorescence light microscopy revealed the localization of QD705 in tissues. RESULTS: Plasma half-life of QD705 in mice was short (18.5 hr), but ICP-MS analyses revealed QD705 persisted and even continued to increase in the spleen, liver, and kidney 28 days after an iv dose. Considerable time-dependent redistribution from body mass to liver and kidney was apparent between 1 and 28 days postdosing. The recoveries at both time points were near 100%; all QD705s reside in the body. Neither fecal nor urinary excretion of QD705 was detected appreciably in 28 days postdosing. Fluorescence microscopy demonstrated deposition of QD705 in the liver, spleen, and kidneys. CONCLUSION: Judging from the continued increase in the liver (29-42% of the administered dose), kidney (1.5-9.2%), and spleen (4.8-5.2%) between 1 and 28 days without any appreciable excretion, QD705 has a very long half-life, potentially weeks or even months, in the body and its health consequences deserve serious consideration.