U.S. flag

An official website of the United States government

Dot gov

The .gov means it's official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you're on a federal government site.


The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Share This:

httk: R Package for High-Throughput Toxicokinetics

Robert G. Pearce, R. Woodrow Setzer, Cory L. Strope, Nisha S. Sipes, John F. Wambaugh.
Journal of Statistical Software (2017) DOI: https://doi.org/10.18637/jss.v079.i04 PMID: Not available



Thousands of chemicals have been profiled by high-throughput screening programs such as ToxCast and Tox21; these chemicals are tested in part because most of them have limited or no data on hazard, exposure, or toxicokinetics. Toxicokinetic models aid in predicting tissue concentrations resulting from chemical exposure, and a "reverse dosimetry" approach can be used to predict exposure doses sufficient to cause tissue concentrations that have been identified as bioactive by high-throughput screening. We have created four toxicokinetic models within the R software package httk. These models are designed to be parameterized using high-throughput in vitro data (plasma protein binding and hepatic clearance), as well as structure-derived physicochemical properties and species-specific physiological data. The package contains tools for Monte Carlo sampling and reverse dosimetry along with functions for the analysis of concentration vs. time simulations. The package can currently use human in vitro data to make predictions for 553 chemicals in humans, rats, mice, dogs, and rabbits, including 94 pharmaceuticals and 415 ToxCast chemicals. For 67 of these chemicals, the package includes rat-specific in vitro data. This package is structured to be augmented with additional chemical data as they become available. Package httk enables the inclusion of toxicokinetics in the statistical analysis of chemicals undergoing high-throughput screening.


Figure 1. Models (A) 1compartment, (B) 3compartment, and (C) pbtk.

In order to preserve mass-balance, Qrest is defined as the difference between Qcardiac and the flow to the liver, kidney, and gut. Variable names are defined in Table 1.

Figure 2. Css at 3 doses per day, 1 mg/kg BW/day.

Figure 3. Days to steady state histogram.

Figure 4. Average vs. maximum concentration at steady state for 1 dose per day, 1 mg/kg BW/day.

Figure 5. Sampling distribution of zoxomide Css from the model pbtk.


Table 1. List of abbreviations.

Table 2. Model parameter and prediction comparison.

*Partition coefficients are needed in calculating Vdist . Clearances and fub are needed in calculating kelim.

Table 3. List of data tables in the package.

In Ring et al. (2017), a series of tables for generating populations based on variation in human physiology were added. They are described in that manuscript and vignettes.

Table 4. List of functions in the package – Part I.

Models are described in Table 2. Parameters are defined in Table 1. Jarnac and SBML are external languages for systems biology models.

Table 5. List of functions in the package – Part II.

Models are described in Table 2. Parameters are defined in Table 1. Jarnac and SBML are external languages for systems biology models.

Supplemental Materials

Supplemental Data