Search PubMed⌕ Search

Biomedical subjects

M E Andersen

Publications and source records attributed to M E Andersen.

At least 55 records · Page 3Linked to original sources

Regional hepatic CYP1A1 and CYP1A2 induction with 2,3,7,8-tetrachlorodibenzo-p-dioxin evaluated with a multicompartment geometric model of hepatic zonation.

A physiologically based pharmacokinetic (PBPK) model for 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) was combined with a five-compartment geometric model of hepatic zonation to predict both total and regional induction of CYP450 proteins within the liver. Three literature studies on TCDD pharmacokinetics and protein induction in female rats were analyzed. In simulating low-dose behavior for mRNA in whole liver and, particularly, in representing immunohistochemical observations, the five-compartment model was more successful than conventional homogeneous one-compartment liver models. The five-compartment liver model was used with the affinity of TCDD for the Ah receptor (AhR) held constant across all the liver (Kb = 0.2 nM). The presumed affinities of the AhR-TCDD complex for TCDD responsive elements in the CYP1A1 (Kd1) and CYP1A2 (Kd2) genes varied between adjacent compartments by a factor of 3. This parameterization leads to predicted 81-fold differences in affinities between the centrilobular and the periportal regions. The affinities used for AhR-TCDD complex binding to TCDD response elements for CYP1A2 in compartment 3 (the midzonal area) ranged from 0.08 to 1.0 nM in the three studies modeled. For CYP1A1 the corresponding dissociation constant in compartment 3 varied from 0.6 to 2.0 nM. In each compartment, the Hill coefficient for induction had to be 4 or greater to match the immunohistochemical results. This multi-compartment liver model is consistent with data on protein and mRNA induction throughout the liver and on the regional distribution of these proteins. No previous model has incorporated regional variations in induction. The PBPK analysis based on the multicompartment liver model suggests that the low-dose behavior for hepatic CYP1A1/CYP1A2 induction by TCDD is highly non-linear.

Animals↗

Glutathione S-transferase-mediated mutagenicity of trihalomethanes in Salmonella typhimurium: contrasting results with bromodichloromethane off chloroform.

Trihalomethanes (THMs) are the most prevalent disinfection by-products identified in chlorinated drinking water. Among the THMs, chloroform (CHCl3) generally occurs at the highest concentration in finished water, but the concentrations of each of the brominated THMs (CHBrCl2, CHBr2Cl, and CHBr3) can exceed that of CHCl3. Each of these four THMs was carcinogenic in rodents in chronic oral dosing studies. This study assessed THM mutagenicity in a strain of Salmonella typhimurium TA1535 that was transfected with rat theta-class glutathione S-transferase T1-1 (+GST). The +GST strain and its nontransfected parent strain (-GST) were employed in a plate-incorporation assay and exposed for 24 hr to the vapor of individual THMs at concentrations up to 25,600 ppm in sealed Tedlar bags. Base-substitution revertants were produced in the +GST strain in a dose-dependent fashion by CHBrCl2 but not by CHCl3. At 4800 ppm CHBrCl2, which produced a calculated agar concentration of 0.67 mM, there were 419 +/- 75 revertants per plate compared to a spontaneous level of 23 +/- 5. CHCl3 produced a doubling of revertants only at the two highest concentrations tested (19,200 and 25,600 ppm). These results indicate that bromination of THMs confers the capability for theta-class GST-mediated transformation to mutagenic intermediates at low substrate concentrations, suggesting the possibility of a similar activation route in humans. Further, the very low affinity of the GSH-dependent pathway for CHCl3 demonstrates that different THMs can induce adverse effects via different mechanisms, indicating that risk evaluations of THMs should not treat members of this class as if they shared a common mode of action.

Animals↗

Hepatic foci in rats after diethylnitrosamine initiation and 2,3,7,8-tetrachlorodibenzo-p-dioxin promotion: evaluation of a quantitative two-cell model and of CYP 1A1/1A2 as a dosimeter.

2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) is a potent hepatic tumor promoter in female rats. We used a quantitative, stochastic initiation-promotion model based on R. B. Conolly and J. S. Kimbell (Toxicol. Appl. Pharmacol. 124, 284-295, 1994) to analyze initiation-promotion results from a previously published study (H. C. Pitot et al., Carcinogenesis 8, 1491-1499, 1987) within the context of a negative selection model of tumor promotion. In this model, two types of initiated cells (called A and B cells) are produced by DEN initiation. Visually excellent correspondence between model predictions and data (i.e., foci/cm3 liver and percentage of liver occupied by foci) are obtained when TCDD is described as having dose-responsive effects on division and death (apoptotic) rates of these two cell types. For A cells, both the division and the death rates increase while the difference between division and apoptotic rates decreases. For B cells, the difference between division and apoptotic rates increases, primarily due to a decrease in the apoptotic rate. We also linked these alterations in cell kinetics to a pharmacokinetic model for TCDD incorporating a five subcompartment model of the liver acinus with induction of CYP1A1 and 1A2 genes in the subcompartments. Alterations in A cell kinetics correlate with effects of TCDD in the region most sensitive to induction (subcompartment 5-centrilobular region); B cell dynamics correlate with induction in subcompartments 3-5 (centrilobular and mid-zonal regions). In summary, these modeling exercises show that (1) the two-cell model, without presuming effects of TCDD on the mutation rate of normal hepatocytes, reproduces the data of Pitot et al. (1987) and (2) induction of CYP1A1/1A2 in different regions of the hepatic acinus can be used as a general correlate of these presumed changes in cell growth kinetics.

Animals↗

A physiologically based pharmacokinetic model for retinoic acid and its metabolites.

BACKGROUND: A physiologically based pharmacokinetic (PBPK) model for all-trans-retinoic acid (tretinoin) was developed to provide a coherent description of tretinoin absorption, distribution metabolism, and excretion across species and routes of administration. OBJECTIVE: The goal of developing such a model is to provide a measure of internal dose that would be a biologically relevant surrogate for administered dose in assessing human teratogenic risk from topically applied tretinoin emollient cream. METHODS: The developed PBPK model included compartments for plasma, liver, gut, intestinal lumen, fat, skin, richly and slowly perfused tissues, placenta, and embryo. Tretinoin metabolism to 13-cis retinoic acid, oxidation, and glucuronidation were incorporated. Dose surrogates, including the maximum plasma concentration (Cmax) and area under the concentration-versus-time curve were calculated from the model. RESULTS: The ability of the model to predict tretinoin pharmacokinetics and to extrapolate across species and routes of administration was tested and validated. Model-derived estimates of dose surrogates demonstrated that the internal exposure to retinoids after topical treatment with 0.05% tretinoin emollient cream is minimal in comparison to that for teratogenic oral doses. The ratio of areas under the curve for total active retinoids after teratogenic oral doses in monkeys versus therapeutic topical doses in human beings, for example, was greater than 450,000 to 1. CONCLUSION: For topical application of tretinoin in human beings, detoxification via the glucuronidation pathway predominates, resulting in a much lower internal exposure to active retinoids than was inferred from total radioactivity data. The model predicts that topical application of tretinoin results in an internal exposure that is four to six orders of magnitude lower than a minimally teratogenic dose.

Abnormalities, Drug-Induced↗

Use of physiologically based pharmacokinetic modeling to investigate individual versus population risk.

Because of the heterogeneity of the human population, it is generally expected that there will be a broad range of observed susceptibilities to the biological effects of exposure to chemicals or drugs. Often it is possible to distinguish specific classes of individuals, such as infants or the elderly, who appear to be more susceptible to a specific effect. Non-cancer risk assessment often address this variability by dividing the experimentally determined acceptable exposure level by an uncertainty factor of 10 to protect sensitive individuals; cancer risk assessments typically do not address this issue in any quantitative fashion. Physiologically based pharmacokinetic (PBPK) modeling provides the capability to quantitatively describe the potential impact of pharmacokinetic factors on the variability of individual risk. In particular, PBPK models can be used to determine the impact of differences in key metabolism enzymes, whether due to multiple genotypic expression, such as cytochrome P450 polymorphisms, or just due to normal variation in enzyme activities within the general population. Other potential modulators of sensitivity which can be addressed quantitatively with a PBPK model include physical condition, level of activity, disease states, age, hormonal status, and interactions with other chemicals and drugs. In each case, the PBPK model provides a quantitative structure for determining the effect of these various factors on the relationship between the external (environmental) exposure and the internal (biologically effective) target tissue exposure. When coupled with Monte Carlo analysis, the PBPK model provides a method to assess the quantitative impact of these sources of variability on individual risk (as opposed to average population risk) by comparing model-predicted risks over the distribution of individual parameter values.

Animals↗

Predicting cancer risk from vinyl chloride exposure with a physiologically based pharmacokinetic model.

A physiologically based pharmacokinetic (PBPK) model capable of describing the metabolism of vinyl chloride (VC) in rats, mice, and humans has been developed and validated by comparison with experimental data from experiments not used in model development. This PBPK model has been used to predict measures of delivered dose (reactive VC metabolites produced in the livers of the affected species) hypothesized to be involved in the induction of liver angiosarcoma in rats, mice, and human populations exposed to VC. Measures of delivered dose in rats were fit to an empirical dose-response model (the linearized multistage model of Crump et al.) and used to make predictions of liver angiosarcoma incidence in mice and human populations exposed to VC. This procedure gave a good prediction of angiosarcoma incidence in mice. Predictions of angiosarcoma incidence in humans were more than two orders of magnitude lower than risk estimations which did not utilize pharmacokinetic data, but were still almost an order of magnitude higher than actually observed in exposed human populations.

Animals↗

Physiologically based pharmacokinetic/pharmacodynamic modeling of the toxicologic interaction between carbon tetrachloride and Kepone.

Carbon tetrachloride (CCl4) lethality in Sprague-Dawley rats is greatly amplified by pretreatment of Kepone (decachlorooctahydro-1,3,2-metheno-2H-cyclobuta[cd] pentalen-2-one). The increase in lethality was attributed to the obstruction of liver regenerative processes. These processes are essential for restoring the liver to its full functional capacity following injury by CCl4. Based on the available mechanistic information on Kepone/CCl4 interaction, a physiologically based pharmacokinetic/pharmacodynamic (PBPK/PD) model was constructed where the following effects of Kepone on CCl4 toxicity are incorporated: (1) inhibition of mitosis; (2) reduction of repair mechanism of hepatocellular injury; (3) suppression of phagocytosis. The PBPK/PD model provided computer simulation consistent with previously published time-course results of hepatotoxicity (i.e., pyknotic, injured and mitotic cells) of CCl4 with or without Kepone. As a further verification of this model, the computer simulations were also consistent with exhalation kinetic data for rats injected with different intraperitoneal (i.p.) doses of CCl4 in our laboratory. Subsequently, the PBPK/PD model, coupled with Monte Carlo simulation, was used to predict lethalities of rats treated with CCl4 alone and CCl4 in combination with Kepone. The experimental lethality studies performed in our laboratories were as follows: Sprague-Dawley rats were given either control diet or diet containing 10 ppm Kepone for 15 days. On day 16, rats in the Kepone treated group were given i.p. doses of 0, 10, 50, and 100 microliters/kg CCl4 (n = 9) while control rats were exposed to 0, 100, 1000, 3000, and 6000 microliters/kg CCl4 (n = 9). Lethality was observed at the 1000 (1/9), 3000 (4/9), and 6000 (8/9) microliters/kg doses for the control group and at the 50 (4/9) and 100 (8/9) microliters/kg for the treated group. Based on Monte Carlo simulation, which was used to run electronically 1000 lethality experiments for each dosing situation, the LD50 estimates for CCl4 toxicity with and without Kepone pretreatment were 47 and 2890 microliters/kg, respectively. Monte Carlo simulation coupled with the PBPK/PD model produced lethality rates which were not significantly different from the observed mortality, with the exception of CCl4 at very high doses (e.g., 6000 microliters/kg, p = 0.014). Deviation at very high doses of the predicted mortality from the observed may be attributed to extrahepatic systemic toxicities of CCl4, or solvent effects on tissues at high concentrations, which were not presently included in the model. Our modeling and experimental results verified the earlier findings of Mehendale (1990) for the 67-fold amplification of CCl4 lethality in the presence of Kepone. However, much of this amplification of CCl4 lethality with Kepone pretreatment was probably due to pharmacokinetic factors, because when target tissue dose (i.e., model estimated amount of CCl4 metabolites) was used to evaluate lethality, this amplification was reduced to 4-fold.

Animals↗

A pharmacokinetic model of anaerobic in vitro carbon tetrachloride metabolism.

Carbon tetrachloride (CCl4) is a potent hepatotoxic agent whose toxicity is mediated through cytochome P450-dependent metabolism. Results from anaerobic in vitro experiments with hepatic microsomes isolated from male F-344 rats indicate that chlorofom (CHCl3) formation from CCl4 is nonlinear with dose. Dose is traditionally expressed as the amount of CCl4 added to the vial. In this study, a pharmacokinetic model has been developed to calculate the concentration of CCl4 in the microsomal suspension. Hepatic microsomes prepared from fed and fasted animals were incubated with CCl4 under anaerobic conditions and formation of CHCl3 over a 5-min incubation period was monitored by headspace gas chromatography. Dose-response curves, based on total amount of CCl4 added to the microsomes, revealed a nonlinear, biphasic appearance of CHCl3, with fasting slightly increasing CHCl3 production in microsomes prepared from fasted rats. Microsomes were also pretreated with the CYP2E1 inhibitor, diallyl sulfone (DAS), before addition of CCl4. In uninhibited microsomes, there appeared to be a high-affinity saturable phase of metabolism occurring at lower concentrations followed by a linear phase at higher CCl4 concentrations. Following DAS pretreatment, the saturable portion of the dose-response curve was inhibited more than the linear phase with the biphasic CHCl3 production becoming more linear. DAS inhibition eliminated the effect of fasting on CHCl3 formation. The best fit kinetic constants for the saturable phase resulted in an estimate of V(max) of 0.017 mg/h/mg protein (V(maxc) = 7.61 mg/h/kg) and Km of 2.3 mg/l (15 microM). The linear phase rate constant (kf) was determined to be 0.046 h-1) (kfc = 0.03 h-1). In conclusion, a pharmacokinetic model has been developed for anaerobic in vitro metabolism of CCl4 to CHCl3 that estimates metabolic rates based on CHCl3 formation and actual CCl4 concentration in the microsomal suspension.

Allyl Compounds↗

PBPK modeling/Monte Carlo simulation of methylene chloride kinetic changes in mice in relation to age and acute, subchronic, and chronic inhalation exposure.

During a 2-year chronic inhalation study on methylene chloride (2000 or 0 ppm; 6 hr/day, 5 days/week), gas-uptake pharmacokinetic studies and tissue partition coefficient determinations were conducted on female B6C3F1, mice after 1 day, 1 month, 1 year, and 2 years of exposure. Using physiologically based pharmacokinetic (PBPK) modeling coupled with Monte Carlo simulation and bootstrap resampling for data analyses, a significant induction in the mixed function oxidase (MFO) rate constant (Vmaxc) was observed at the 1-day and 1-month exposure points when compared to concurrent control mice while decreases in glutathione S-transferase (GST) rate constant (Kfc) were observed in the 1-day and 1-month exposed mice. Within exposure groups, the apparent Vmaxc maintained significant increases in the 1-month and 2-year control groups. Although the same initial increase exists in the exposed group, the 2-year Vmaxc is significantly smaller than the 1-month group (p < 0.001). Within group differences in median Kfc values show a significant decrease in both 1-month and 2-year groups among control and exposed mice (p < 0.001). Although no changes in methylene chloride solubility as a result of prior exposure were observed in blood, muscle, liver, or lung, a marginal decrease in the fat:air partition coefficient was found in the exposed mice at p = 0.053. Age related solubility differences were found in muscle:air, liver:air, lung:air, and fat:air partition coefficients at p < 0.001, while the solubility of methylene chloride in blood was not affected by age (p = 0.461). As a result of this study, we conclude that age and prior exposure to methylene chloride can produce notable changes in disposition and metabolism and may represent important factors in the interpretation for toxicologic data and its application to risk assessment.

Administration, Inhalation↗

Biologically based, quantitative risk assessment of neurotoxicants.

The need for biologically based, quantitative risk assessment procedures for noncancer endpoints such as neurotoxicity has been discussed in reports by the United States Congress (Office of Technology Assessment, OTA), National Research Council (NRC), and a federal coordinating council. According to OTA, current attention and resources allocated to health risk assessment research are inadequate and not commensurate with its impact on public health and the economy. Methods to include continuous rather than dichotomous data for neurotoxicity endpoints, biomarkers of exposure and effects, and pharmacokinetic and mechanistic data have been proposed for neurotoxicity risk assessment but require further review and validation before acceptance. The purpose of this symposium was to examine procedures to enhance the risk assessment process for neurotoxicants and to discuss techniques to make the process more quantitative. Accordingly, a review of the currently used safety factor risk assessment approach for neurotoxicants is provided along with specific examples of how this process may be enhanced with the use of the benchmark dose approach. The importance of including physiologically based pharmacokinetic data in the risk assessment process and specific examples of this approach is presented for neurotoxicants. The role of biomarkers of exposure and effect and mechanistic information in the risk assessment process are also addressed. Finally, quantitative approaches with the use of continuous neurotoxicity data are demonstrated and the outcomes compared to those generated by currently used risk assessment procedures.

Biomarkers↗

Negative selection in hepatic tumor promotion in relation to cancer risk assessment.

Mechanistic studies with phenobarbital (PB), 2,3,7,8,-tetrachlorodibenzo-p-dioxin (TCDD) and other liver tumor promoters support a general model of promotion involving negative selection where specifically-mutated cells derive a growth advantage in the presence of persistent mitosuppression. Exposure to these liver tumor promoters appears to transiently enhance hepatocyte replication, presumably via transcriptional activation of growth regulatory genes, leading to a homeostatic increase in mitoinhibitory growth factors in the liver to constrain proliferation. Transforming growth factor beta 1 (TGF-beta), a potent mitoinhibitory growth factor for hepatocytes, has been associated with the mitosuppression caused by PB and certain peroxisomal proliferators. Escape from TGF-beta mitosuppression may involve loss or alteration of function of the mannose 6-phosphate/insulin-like growth factor II (M6P/IGFII) receptor, which is required for TGF-beta 1 activation, or alterations of the TGF-beta types I, II and III signal transduction receptors. A risk assessment based on a negative selection mechanism could be conducted for tumor promotion endpoints with TCDD and compared with current approaches that implicitly regard TCDD as an initiator. Benchmark dose calculation using centrilobular induction of cytochromes P450 1A1 and 1A2 as a surrogate for periportal growth stimulation would provide a rational starting point for application of conventional safety factor approaches, similar to those used with non-cancer effects. In the future, tissue and plasma concentrations of specific growth factors, e.g. TGF-beta or hepatocyte growth factor, HGF, might be considered as more direct dose surrogates for tumor-promoting effects of xenobiotics. Uncertainty factor adjustments to a TCDD benchmark dose calculation should eventually rely on direct knowledge of regulation of specific growth regulatory genes and their receptors in relevant species and on species differences in TCDD pharmacokinetics, instead of application of default animal-to-human and interindividual uncertainty factors.

Animals↗

Development of a mechanism-based dosimetry model for 2,4,4-trimethyl-2-pentanol-induced alpha 2u-globulin nephropathy in male Fischer 344 rats.

A mechanism-based dosimetry model was developed to describe 2,4,4-trimethyl-2-pentanol (TMP-2-OH) dosimetry and renal alpha 2u-globulin (alpha 2u) nephropathy in the male Fischer 344 rat. Experimental data were collected to estimate the chemical-specific parameters (metabolic constants, tissue solubility, and oral absorption rate) necessary to describe TMP-2-OH dosimetry in male rats. The concentrations of alpha 2u and TMP-2-OH were measured in male rats up to 64 hr after a single oral dose of TMP-2-OH (6, 60, or 600 mg/kg). The model predicted the time course behavior of TMP-2-OH and alpha 2u in the kidney, but overestimated their renal concentrations by two or threefold. Simulations of renal alpha 2u concentration were sensitive to changes in TMP-2-OH-alpha 2u-binding affinity and degradation rate of the TMP-2-OH-protein complex. In contrast, simulation of the concentration of TMP-2-OH in the kidney was most sensitive to the amount of protein present. Oral absorption of TMP-2-OH was dose dependent. The model predicted that alpha 2u and TMP-2-OH concentration in the kidney is sensitive to changes in the rate of TMP-2-OH absorbed after oral administration. This model permitted a more rigorous evaluation than has previously been possible of the combination of protein characteristics and chemical dosimetry required for the accumulation of alpha 2u in the kidney of male rats. The behavior of the model is consistent with the qualitative aspects of the alpha 2u hypothesis. However, further characterization of alpha 2u distribution and renal hydrolysis will be required in order to fully characterize the hypothesis at the quantitative level.

Administration, Oral↗

Applying simulation modeling to problems in toxicology and risk assessment--a short perspective.

The goals of this perspective have been to examine areas where quantitative simulation models may be useful in toxicology and related risk assessment fields and to offer suggestions for preparing manuscripts that describe these models. If developments in other disciplines serve as a bell-wether, the use of mathematical models in toxicology will continue to increase, partly, at least, because the new generations of scientists are being trained in an electronic environment where computation of all kinds is learned at an early age. Undoubtedly, however, the utility of these models will be directly tied to the skills of investigators in accurately describing models in their research papers. These publications should convey descriptions of both the insights obtained and the opportunities provided by these models to integrate existing data bases and suggest new and useful experiments. We hope these comments serve to facilitate the expansion of good modeling practices as applied to toxicological problems.

Computer Simulation↗

A physiologically based dosimetry description of acrylonitrile and cyanoethylene oxide in the rat.

The cytochrome P450-mediated oxidation of acrylonitrile (ACN) to the mutagen 2-cyanoethylene oxide (CEO) is thought to be important for the carcinogenic effects of ACN in rats, while glutathione (GSH) conjugation of ACN and CEO is regarded as detoxication. A physiologically based dosimetry description for ACN and CEO in the male F-344 rat has been developed from in vitro data and studies of the iv pharmacokinetics of ACN and CEO. The dosimetry description includes tissue partition coefficients and in vitro estimates of the rates of reaction of ACN and CEO with hemoglobin and blood macromolecules and the reaction of CEO with tissue GSH. Metabolic parameters for ACN and CEO were estimated from iv pharmacokinetic studies. Rats were given bolus doses of 3.4, 47, 55, or 84 mg ACN/kg via the femoral vein and blood samples were collected at selected time points. ACN and CEO blood concentrations were determined by gas chromatography. The iv pharmacokinetics of CEO were also determined using 0.6 or 5.3 mg CEO/kg. ACN elimination from blood was described by saturable P450 epoxidation (Vmax of 6.5 mg/hr/kg and Km of 1.5 mg/liter) and first-order GSH conjugation (30 hr-1/kg). CEO elimination was described by first-order GSH conjugation (750 hr-1/kg). Calculation of hepatic clearance values shows first-pass hepatic extractions of 61 and 90% for ACN and CEO, respectively. The dosimetry description accurately simulated the dose-dependent urinary excretion of ACN metabolites derived from epoxidation to CEO and from direct GSH conjugation of ACN. The dose-dependent formation of hemoglobin adducts from ACN was also well simulated.

Acrylonitrile↗

Considering pharmacokinetic and mechanistic information in cancer risk assessments for environmental contaminants: examples with vinyl chloride and trichloroethylene.

Risk assessments for vinyl chloride (VC) and trichloroethylene (TCE) are presented as examples of approaches for incorporating chemical-specific pharmacokinetic and mechanistic information into a more scientifically plausible cancer risk assessment. For VC, the evidence regarding mode of action includes direct reaction of a metabolite with DNA, resulting in DNA adducts and mistranscription, and cross-species target-tissue correspondence of a rare tumor type. Risk estimates for human exposure to VC predicted with a physiologically-based pharmacokinetic (PBPK) model and the linearized multistage (LMS) model were lower than those currently used in environmental decision-making by a factor of 30 to 50, and were more consistent with human epidemiological data. For TCE, there is evidence of increased cell proliferation due to receptor interaction or cytotoxicity in every instance in which tumors are observed, and the tumors typically represent an increase in the incidence of a commonly observed, species-specific lesion. Virtually safe exposure estimates for human exposure to TCE predicted with a PBPK model and a margin of exposure (MOE) approach were higher than those obtained by the conventional LMS approach by roughly a factor of 100. The MOE approach is recommended as an alternative to the LMS approach for chemicals with a carcinogenic mode of action which entails increased cell proliferation, leading to the expectation of a highly nonlinear cancer dose-response.

Animals↗

Development of physiologically based pharmacokinetic and physiologically based pharmacodynamic models for applications in toxicology and risk assessment.

Pharmacokinetics (PK) involves the study of the rates of absorption, distribution, excretion, and biotransformation of chemicals and their metabolites. PK models can be used to reconstruct extensive time-course data sets based on a small number of kinetic parameters. These models can be used to predict the results of new experiments and integrate studies on kinetics, disposition and metabolism in various animal species [1]. The 2 main approaches that have been pursued in developing PK models are: (1) data-based compartmental modeling; and (2) physiologically based compartmental modeling. Data-based models rely on the collection of time-course concentration data and fitting these data with mathematical models. Compartments in these models do not necessarily reflect the anatomy and physiology of the animal, and the kinetic constants derived from these models do not have obvious physiological or biochemical counterparts. In physiologically based pharmacokinetic (PBPK) models, compartments correspond more closely to actual anatomical structures, defined with respect to their volumes, blood flows, chemical binding (partitioning) characteristics, and ability to metabolize or excrete the compounds of interest. Because the kinetic parameters of these models reflect tissue blood flows, partitioning, and biochemical constants, these models are more readily scaled from one animal species to another [2]. PBPK models have been used to understand the disposition of chemicals in the body for almost 70 years. Their more widespread application in toxicology dates back only 15 years or so to models developed for polychlorinated biphenyls and other persistent lipophilic compounds. Quantitative applications of PBPK models in risk assessment date to the development of a number of PBPK models for methylene chloride in the mid 1980s. The burgeoning use of PBPK models in toxicology research and chemical risk assessment today is primarily related to their ability to make more accurate predictions of target tissue dose for different exposure situations in different animal species, including humans. This overview includes a discussion of the development of these PBPK models in toxicology and speculates about future applications of PBPK and physiologically based pharmacodynamic (PBPD) models in chemical risk assessment.

Animals↗

Physiologically based pharmacokinetic (PB-PK) models in the study of the disposition and biological effects of xenobiotics and drugs.

Physiologically based pharmacokinetic (PB-PK) models have been used to study the mechanisms of disposition of drugs and xenobiotics for almost 70 years. Their widespread application in toxicology began 15 years ago with models for polychlorinated biphenyls and other persistent lipophilic compounds. Quantitative applications of PB-PK moels in carcinogen risk assessment date to the development of a number of PB-PK models for dichloromethane in the mid 1980s. The expanding use of these models is primarily related to their ability to make more accurate predictions of target tissue dose for different exposure situations in different animal species, including humans, and to evaluate quantitatively the mechanisms of disposition of chemicals within the body. This paper discusses contemporary uses of PB-PK modeling in the context of risk assessment with xenobiotics and of safety assessment with drugs.

Animals↗

An improved approximation to the exact solution of the two-stage clonal growth model of cancer.

Multistage clonal growth models are of interest for cancer risk assessment because they can explicitly incorporate data on cell replication. Both approximate and exact formulations of the two stage growth model have been described. The exact solution considers the conditional probability of tumors arising in previously tumor-free animals; the approximate solution estimates total probability of tumor formation. The exact solution is much more computationally intensive when time-dependent cell growth parameters are included. The approximate solution deviates from the exact solution at high incidences and probabilities of tumor. This report describes a computationally tractable, 'improved approximation' to the exact solution. Our improved approximation includes a correction term to adjust the unconditional expectation of intermediate cells based on the time history of formation of intermediate cells by mutation of normal cells (recruitment) or by cell division in the intermediate cell population (expansion). The improved approximation provided a much better match to the exact solution than the approximate solution for a wide range of parameter values. The correction term also appears to provide insight into the biological factors that contribute to the variance of the expectation for the number of intermediate cells over time.

Algorithms↗