Mechanistic modeling of rodent liver tumor promotion at low levels of exposure: an example related to dose-response relationships for 2,3,7,8-tetrachlorodibenzo-p-dioxin

Computational experiments reveal plausible mechanisms for changing patterns of hepatic zonation of xenobiotic clearance and hepatotoxicity

No concrete, causal, mechanistic theory is available to explain how different hepatic zonation patterns of P450 isozyme levels and hepatotoxicity emerge following dosing with different compounds. We used the synthetic method of modeling and simulation to discover, explore, and experimentally challenge concrete mechanisms that show how and why biomimetic zonation patterns can emerge and change within agent-based analogues, expecting that those mechanisms may have counterparts in rats. Mobile objects map to compounds. One analogue represents a cross-section through a lobule. It is comprised of 460 identical, quasi-autonomous functional units called sinusoidal segments (SSs). SSs detect and respond to compound-generated response signals and the local level of an endogenous gradient. Each SS adapts by using those signals to adjust (or not) the probability that it will clear a detected compound during the next simulation cycle. The adjustment decision is based on the value of a biomimetic algorithm that is based on an assumed, evolution imposed, genetic mandate that normal hepatocytes resist increasing the cost of their actions. The algorithm estimates the long-term, discounted cost to a given SS of continuing to use its current clearance effort. Upon compound exposure, lobular analogues developed a variety of clearance and hepatotoxicity patterns that were strikingly similar to those reported in the literature. A degree of quantitative validation was achieved against data on hepatic zonation of CYP1A2 mRNA expression caused by three different doses of TCDD (2,3,7,8-tetracholorodibenzo-p-dioxone).

Hormesis and adaptive cellular control systems

Hormetic dose response occurs for many endpoints associated with exposures of biological organisms to environmental stressors. Cell-based U- or inverted U-shaped responses may derive from common processes involved in activation of adaptive responses required to protect cells from stressful environments. These adaptive pathways extend the region of cellular homeostasis and are protective against ultimate cell, organ, and system toxicity. However, the activation of stress responses carries a significant energetic cost to the cell, leading to alterations of a variety of basal cellular functions in adapted or stressed cells. This tradeoff of resources between the unstressed and adapted states may lead to U-or inverted U-shaped dose response curves for some precursor endpoints. We examine this general hypothesis with chlorine, a prototype oxidative stressor, using a combination of cellular studies with gene expression analysis of response pathways and with computational modeling of activation of control networks. Discrete cellular states are expected as a function of exposure concentration and duration. These cellular states include normal functioning state, adaptive and stressed states at mild to intermediate exposures, and overt toxicity in the presence of an overwhelming concentration of stressors. These transitions can be used to refine default risk assessment practices that do not currently accommodate adaptive responses.

Dioxin cancer risk--example of hormesis?

A recent case-control study implied an inverse correlation between the measured body burden of dioxins (polychlorinated dibenzo-p-dioxins and polychlorinated dibenzofurans, PCDD/F) and the risk of soft tissue sarcoma in normal population exposed to dioxins mainly via food. The surprising result could not be explained by biases or confounding. There is no a priori confounding by occupational chemicals in a random sample from general population, but exposures to other lipid soluble chemicals with similar sources might be expected to associate with that of dioxins. One such group is polychlorinated biphenyls (PCB). Therefore three most relevant dioxin-like PCB compounds PCB 77, PCB 126, and PCB 169 were now analyzed from the same patients. Cases were 110 soft-tissue sarcoma patients undergoing surgery for their disease, and referents were 227 patients operated for appendicitis. Dioxin and PCB concentrations were analyzed from subcutaneous fat samples by high-resolution gas chromatography-mass spectrometry and TCDD equivalent concentrations (WHO-TEq) were calculated by using toxicity equivalency factors of WHO. The highest risk of sarcoma was found in the septile with the lowest body burden of sum WHO-TEq, and the differences of septiles 2 and 6 from septile 1 were statistically significant. If soft sarcoma risk is true at high occupational levels of dioxins, the provocative result suggests that a possibility of a J-shaped dose-response curve should be taken into consideration and studied further. This is also supported by the similar J-shaped dose responses in animal studies.

The use of quantitative histological and molecular data for risk assessment and biologically based model development

In organs with diverse cell populations, it is not uncommon for one type of cell to respond while others are spared. Even in an organ with common cell types, such as hepatocytes within the liver, the population of cells may respond with different sensitivities for injury or for biochemical responses to toxicants. In the liver, many tumor promoters induce cytochrome P450 enzymes and other proteins in centrilobular cells at much lower doses than required to cause induction in periportal cells. In addition, these induction responses appear to occur at the level of individual cells--a 50% response of the liver for induction does not represent 50% induction in all cells. Instead, half of the cells are fully induced and half are unaffected. Cells "switch" from one phenotypic state to another. Over the past 10 years, several attempts have been made to model these cellular switches and to understand their relevance for hepatic tumor promotion and risk assessment. The data used for analyzing these switches include responses of the entire liver (total induction), responses of individual cells in the liver (regional induction), and cellular responses such as proliferation and apoptosis. This brief overview describes the development of biologically based, dose-response (BBDR) models for protein induction and tumor promotion in liver by 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) with emphasis on the role of specific types of histological and molecular data in providing insights about mechanisms for cellular switches and their implications for tumor promotion. As the biological basis of these switches become unraveled and incorporated into the models, these BBDR models should eventually serve to improve risk assessments with a variety of liver tumor promoters with receptor-based modes of action.