parameterize_1comp {httk} | R Documentation |
Parameters for a one compartment (empirical) toxicokinetic model
Description
This function initializes the parameters needed in the function solve_1comp. Volume of distribution is estimated by using a modified Schmitt (2008) method to predict tissue particition coefficients (Pearce et al., 2017) and then lumping the compartments weighted by tissue volume:
Usage
parameterize_1comp(
chem.cas = NULL,
chem.name = NULL,
dtxsid = NULL,
species = "Human",
default.to.human = FALSE,
adjusted.Funbound.plasma = TRUE,
adjusted.Clint = TRUE,
regression = TRUE,
restrictive.clearance = TRUE,
well.stirred.correction = TRUE,
suppress.messages = FALSE,
clint.pvalue.threshold = 0.05,
minimum.Funbound.plasma = 1e-04,
class.exclude = TRUE,
physchem.exclude = TRUE,
Caco2.options = list(),
...
)
Arguments
chem.cas |
Chemical Abstract Services Registry Number (CAS-RN) – the chemical must be identified by either CAS, name, or DTXISD |
chem.name |
Chemical name (spaces and capitalization ignored) – the chemical must be identified by either CAS, name, or DTXISD |
dtxsid |
EPA's DSSTox Structure ID (https://comptox.epa.gov/dashboard) – the chemical must be identified by either CAS, name, or DTXSIDs |
species |
Species desired (either "Rat", "Rabbit", "Dog", "Mouse", or default "Human"). |
default.to.human |
Substitutes missing rat values with human values if true. |
adjusted.Funbound.plasma |
Uses Pearce et al. (2017) lipid binding adjustment for Funbound.plasma (which impacts volume of distribution) when set to TRUE (Default). |
adjusted.Clint |
Uses Kilford et al. (2008) hepatocyte incubation binding adjustment for Clint when set to TRUE (Default). |
regression |
Whether or not to use the regressions in calculating partition coefficients in volume of distribution calculation. |
restrictive.clearance |
In calculating elimination rate and hepatic bioavailability, protein binding is not taken into account (set to 1) in liver clearance if FALSE. |
well.stirred.correction |
Uses correction in calculation of hepatic clearance for well-stirred model if TRUE. This assumes clearance relative to amount unbound in whole blood instead of plasma, but converted to use with plasma concentration. |
suppress.messages |
Whether or not to suppress messages. |
clint.pvalue.threshold |
Hepatic clearance for chemicals where the in vitro clearance assay result has a p-value greater than the threshold are set to zero. |
minimum.Funbound.plasma |
Monte Carlo draws less than this value are set equal to this value (default is 0.0001 – half the lowest measured Fup in our dataset). |
class.exclude |
Exclude chemical classes identified as outside of domain of applicability by relevant modelinfo_[MODEL] file (default TRUE). |
physchem.exclude |
Exclude chemicals on the basis of physico-chemical properties (currently only Henry's law constant) as specified by the relevant modelinfo_[MODEL] file (default TRUE). |
Caco2.options |
A list of options to use when working with Caco2 apical
to basolateral data |
... |
Additional arguments, not currently used. |
Details
V_d,steady-state = Sum over all tissues (K_i * V_i) + V_plasma
where K_i is the tissue:unbound plasma concentration partition coefficient for tissue i.
Because this model does not simulate exhalation, inhalation, and other processes relevant to volatile chemicals, this model is by default restricted to chemicals with a logHenry's Law Constant less than that of Acetone, a known volatile chemical. That is, chemicals with logHLC > -4.5 (Log10 atm-m3/mole) are excluded. Volatility is not purely determined by the Henry's Law Constant, therefore this chemical exclusion may be turned off with the argument "physchem.exclude = FALSE". Similarly, per- and polyfluoroalkyl substances (PFAS) are excluded by default because the transporters that often drive PFAS toxicokinetics are not included in this model. However, PFAS chemicals can be included with the argument "class.exclude = FALSE".
Value
Vdist |
Volume of distribution, units of L/kg BW. |
Fabsgut |
Fraction of the oral dose absorbed and surviving gut metabolism, i.e. the fraction of the dose that enters the gutlumen. |
kelim |
Elimination rate, units of 1/h. |
hematocrit |
Percent volume of red blood cells in the blood. |
Fabsgut |
Fraction of the oral dose absorbed, i.e. the fraction of the dose that enters the gutlumen. |
Fhep.assay.correction |
The fraction of chemical unbound in hepatocyte assay using the method of Kilford et al. (2008) |
kelim |
Elimination rate, units of 1/h. |
hematocrit |
Percent volume of red blood cells in the blood. |
kgutabs |
Rate chemical is absorbed, 1/h. |
million.cells.per.gliver |
Millions cells per gram of liver tissue. |
MW |
Molecular Weight, g/mol. |
Rblood2plasma |
The ratio of the concentration of the chemical in the blood to the concentration in the plasma. Not used in calculations but included for the conversion of plasma outputs. |
hepatic.bioavailability |
Fraction of dose remaining after first pass clearance, calculated from the corrected well-stirred model. |
BW |
Body Weight, kg. |
Author(s)
John Wambaugh and Robert Pearce
References
Pearce RG, Setzer RW, Strope CL, Wambaugh JF, Sipes NS (2017). “Httk: R package for high-throughput toxicokinetics.” Journal of Statistical Software, 79(4), 1. doi:10.18637/jss.v079.i04.
Schmitt W (2008). “General approach for the calculation of tissue to plasma partition coefficients.” Toxicology in vitro, 22(2), 457–467. doi:10.1016/j.tiv.2007.09.010.
Pearce RG, Setzer RW, Davis JL, Wambaugh JF (2017). “Evaluation and calibration of high-throughput predictions of chemical distribution to tissues.” Journal of pharmacokinetics and pharmacodynamics, 44, 549–565. doi:10.1007/s10928-017-9548-7.
Kilford PJ, Gertz M, Houston JB, Galetin A (2008). “Hepatocellular binding of drugs: correction for unbound fraction in hepatocyte incubations using microsomal binding or drug lipophilicity data.” Drug Metabolism and Disposition, 36(7), 1194–1197. doi:10.1124/dmd.108.020834.
See Also
Examples
parameters1 <- parameterize_1comp(chem.name='Bisphenol-A',species='Rat')
parameters2 <- parameterize_1comp(chem.cas='80-05-7',
restrictive.clearance=FALSE,
species='rabbit',
default.to.human=TRUE)
# The following will not work because Diquat dibromide monohydrate's
# Henry's Law Constant (-3.912) is higher than that of Acetone (~-4.5):
try(parameters3 <- parameterize_1comp(chem.cas = "6385-62-2"))
# However, we can turn off checking for phys-chem properties, since we know
# that Diquat dibromide monohydrate is not too volatile:
parameters3 <- parameterize_1comp(chem.cas = "6385-62-2",
physchem.exclude = FALSE)
out <- solve_1comp(parameters=parameters1, days=1)