A Simple Toxicokinetic Model Exhibiting Complex Dynamics and Nonlinear Exposure Response

Uncertainty in model predictions of exposure response at low exposures is a problem for risk assessment. A particular interest is the internal concentration of an agent in biological systems as a function of external exposure concentrations. Physiologically based pharmacokinetic (PBPK) models permit estimation of internal exposure concentrations in target tissues but most assume that model parameters are either fixed or instantaneously dose-dependent. Taking into account response times for biological regulatory mechanisms introduces new dynamic behaviors that have implications for low-dose exposure response in chronic exposure. A simple one-compartment simulation model is described in which internal concentrations summed over time exhibit significant nonlinearity and nonmonotonicity in relation to external concentrations due to delayed up- or downregulation of a metabolic pathway. These behaviors could be the mechanistic basis for homeostasis and for some apparent hormetic effects.     

Cancer Dose-Response Modeling of Epidemiological Data on Worker Exposures to Aldrin and Dieldrin

The paper applies classical statistical principles to yield new tools for risk assessment and makes new use of epidemiological data for human risk assessment. An extensive clinical and epidemiological study of workers engaged in the manufacturing and formulation of aldrin and dieldrin provides occupational hygiene and biological monitoring data on individual exposures over the years of employment and provides unusually accurate measures of individual lifetime average daily doses. In the cancer dose-response modeling, each worker is treated as a separate experimental unit with his own unique dose. Maximum likelihood estimates of added cancer risk are calculated for multistage, multistage-Weibull, and proportional hazards models. Distributional characterizations of added cancer risk are based on bootstrap and relative likelihood techniques. The cancer mortality data on these male workers suggest that low-dose exposures to aldrin and dieldrin do not significantly increase human cancer risk and may even decrease the human hazard rate for all types of cancer combined at low doses (e.g., 1 g/kg/day). The apparent hormetic effect in the best fitting dose-response models for this data set is statistically significant. The decrease in cancer risk at low doses of aldrin and dieldrin is in sharp contrast to the U.S. Environmental Protection Agency's upper bound on cancer potency based on mouse liver tumors. The EPA's upper bound implies that lifetime average daily doses of 0.0000625 and 0.00625 g/kg body weight/day would correspond to increased cancer risks of 0.000001 and 0.0001, respectively. However, the best estimate from the Pernis epidemiological data is that there is no increase in cancer risk in these workers at these doses or even at doses as large as 2 g/kg/day.     

Logistic Regression Models with Unspecified Low Dose–Response Relationships and Experimental Designs for Hormesis Studies

Hormesis refers to a nonmonotonic (biphasic) dose–response relationship in toxicology, environmental science, and related fields. In the presence of hormesis, a low dose of a toxic agent may have a lower risk than the risk at the control dose, and the risk may increase at high doses. When the sample size is small due to practical, logistic, and ethical considerations, a parametric model may provide an efficient approach to hypothesis testing at the cost of adopting a strong assumption, which is not guaranteed to be true. In this article, we first consider alternative parameterizations based on the traditional three‐parameter logistic regression. The new parameterizations attempt to provide robustness to model misspecification by allowing an unspecified dose–response relationship between the control dose and the first nonzero experimental dose. We then consider experimental designs including the uniform design (the same sample size per dose group) and the $c$‐optimal design (minimizing the standard error of an estimator for a parameter of interest). Our simulation studies showed that (1) the $c$‐optimal design under the traditional three‐parameter logistic regression does not help reducing an inflated Type I error rate due to model misspecification, (2) it is helpful under the new parameterization with three parameters (Type I error rate is close to a fixed significance level), and (3) the new parameterization with four parameters and the $c$‐optimal design does not reduce statistical power much while preserving the Type I error rate at a fixed significance level.     

Estimation of a Benchmark Dose in the Presence or Absence of Hormesis Using Posterior Averaging

U.S. Environment Protection Agency benchmark doses for dichotomous cancer responses are often estimated using a multistage model based on a monotonic dose‐response assumption. To account for model uncertainty in the estimation process, several model averaging methods have been proposed for risk assessment. In this article, we extend the usual parameter space in the multistage model for monotonicity to allow for the possibility of a hormetic dose‐response relationship. Bayesian model averaging is used to estimate the benchmark dose and to provide posterior probabilities for monotonicity versus hormesis. Simulation studies show that the newly proposed method provides robust point and interval estimation of a benchmark dose in the presence or absence of hormesis. We also apply the method to two data sets on carcinogenic response of rats to 2,3,7,8‐tetrachlorodibenzo‐p‐dioxin.     

Hormesis Without Cell Killing

Stochastic two-stage clonal expansion (TSCE) models of carcinogenesis offer the following clear theoretical explanation for U-shaped cancer dose-response relations. Low doses that kill initiated (premalignant) cells thereby create a protective effect. At higher doses, this effect is overwhelmed by an increase in the net number of initiated cells. The sum of these two effects, from cell killing and cell proliferation, gives a U-shaped or J-shaped dose-response relation. This article shows that exposures that do not kill, repair, or decrease cell populations, but that only hasten transitions that lead to cancer, can also generate U-shaped and J-shaped dose-response relations in a competing-risk (modified TSCE) framework where exposures disproportionately hasten transitions into carcinogenic pathways with relatively long times to tumor. Quantitative modeling of the competing effects of more transitions toward carcinogenesis (risk increasing) and a higher proportion of transitions into the slower pathway (risk reducing) shows that a J-shaped dose-response relation can occur even if exposure increases the number of initiated cells at every positive dose level. This suggests a possible new explanation for hormetic dose-response relations in response to carcinogenic exposures that do not have protective (cell-killing) effects. In addition, the examples presented emphasize the role of time in hormesis: exposures that monotonically increase risks at younger ages may nonetheless produce a U-shaped or J-shaped dose-response relation for lifetime risk of cancer.     

A Parametric Model for Detecting Hormetic Effects in Developmental Toxicity Studies

Hormetic effects have been observed at low exposure levels based on the dose-response pattern of data from developmental toxicity studies. This indicates that there might actually be a reduced risk of exhibiting toxic effects at low exposure levels. Hormesis implies the existence of a threshold dose level and there are dose-response models that include parameters that account for the threshold. We propose a function that introduces a parameter to account for hormesis. This function is a subset of the set of all functions that could represent a hormetic dose-response relationship at low exposure levels to toxic agents. We characterize the overall dose-response relationship with a piecewise function that consists of a hormetic u-shape curve at low dose levels and a logistic curve at high dose levels. We apply our model to a data set from an experiment conducted at the National Toxicology Program (NTP). We also use the beta-binomial distribution to model the litter response data. It can be seen by observing the structure of these data that current experimental designs for developmental studies employ a limited number of dose groups. These designs may not be satisfactory when the goal is to illustrate the existence of hormesis. In particular, increasing the number of low-level doses improves the power for detecting hormetic effects. Therefore, we also provide the results of simulations that were done to characterize the power of current designs in detecting hormesis and to demonstrate how this power can be improved upon by altering these designs with the addition of only a few low exposure levels.     

Research Integrity and Conflicts of Interest: The Case of Unethical Research-Misconduct Charges Filed by Edward Calabrese

Special-interest polluters often file research-misconduct (RM) charges against scientists whose research suggests needed pollutant regulation. This article argues that U.S. RM regulations are flawed in requiring RM assessors/experts/accused, but not accusers, to reveal possible conflicts of interest (COI) that could affect RM allegations. It (1) summarizes U.S. RM regulatory history; (2) uses a case study about 2011 RM allegations, filed by chemical-industry-funded toxicologist Edward Calabrese, to illustrate problems with RM regulations; and (3) offers 4 arguments in favor of revising RM regulations so as to require RM-accuser revelation of possible COI and who funded preparation of the RM allegations.