In vitro to in vivo extrapolation for trichloroethylene metabolism in humans

The Application of a Physiologically Based Toxicokinetic Model in Health Risk Assessment

Toxicokinetics plays a crucial role in the health risk assessments of xenobiotics. Classical compartmental models are limited in their ability to determine chemical concentrations in specific organs or tissues, particularly target organs or tissues, and their limited interspecific and exposure route extrapolation hinders satisfactory health risk assessment. In contrast, physiologically based toxicokinetic (PBTK) models quantitatively describe the absorption, distribution, metabolism, and excretion of chemicals across various exposure routes and doses in organisms, establishing correlations with toxic effects. Consequently, PBTK models serve as potent tools for extrapolation and provide a theoretical foundation for health risk assessment and management. This review outlines the construction and application of PBTK models in health risk assessment while analyzing their limitations and future perspectives.

Convergent somatic mutations in metabolism genes in chronic liver disease

The progression of chronic liver disease to hepatocellular carcinoma is caused by the acquisition of somatic mutations that affect 20-30 cancer genes1-8. Burdens of somatic mutations are higher and clonal expansions larger in chronic liver disease9-13 than in normal liver13-16, which enables positive selection to shape the genomic landscape9-13. Here we analysed somatic mutations from 1,590 genomes across 34 liver samples, including healthy controls, alcohol-related liver disease and non-alcoholic fatty liver disease. Seven of the 29 patients with liver disease had mutations in FOXO1, the major transcription factor in insulin signalling. These mutations affected a single hotspot within the gene, impairing the insulin-mediated nuclear export of FOXO1. Notably, six of the seven patients with FOXO1S22W hotspot mutations showed convergent evolution, with variants acquired independently by up to nine distinct hepatocyte clones per patient. CIDEB, which regulates lipid droplet metabolism in hepatocytes17-19, and GPAM, which produces storage triacylglycerol from free fatty acids20,21, also had a significant excess of mutations. We again observed frequent convergent evolution: up to fourteen independent clones per patient with CIDEB mutations and up to seven clones per patient with GPAM mutations. Mutations in metabolism genes were distributed across multiple anatomical segments of the liver, increased clone size and were seen in both alcohol-related liver disease and non-alcoholic fatty liver disease, but rarely in hepatocellular carcinoma. Master regulators of metabolic pathways are a frequent target of convergent somatic mutation in alcohol-related and non-alcoholic fatty liver disease.

Hepatic Scaling Factors for In Vitro-In Vivo Extrapolation of Metabolic Drug Clearance in Patients with Colorectal Cancer with Liver Metastasis

In vitro-in vivo extrapolation (IVIVE) linked with physiologically based pharmacokinetics (PBPK) modeling is used to predict the fates of drugs in patients. Ideally, the IVIVE-PBPK models should incorporate systems information accounting for characteristics of the specific target population. There is a paucity of such scaling factors in cancer, particularly microsomal protein per gram of liver (MPPGL) and cytosolic protein per gram of liver (CPPGL). In this study, cancerous and histologically normal liver tissue from 16 patients with colorectal liver metastasis were fractionated to microsomes and cytosol. Protein content was measured in homogenates, microsomes, and cytosol. The loss of microsomal protein during fractionation was accounted for using corrections based on NADPH cytochrome P450 reductase activity in different matrices. MPPGL was significantly lower in cancerous tissue (24.8 ± 9.8 mg/g) than histologically normal tissue (39.0 ± 13.8 mg/g). CPPGL in cancerous tissue was 42.1 ± 12.9 mg/g compared with 56.2 ± 16.9 mg/g in normal tissue. No correlations between demographics (sex, age, and body mass index) and MPPGL or CPPGL were apparent in the data. The generated scaling factors together with assumptions regarding the relative volumes of cancerous versus noncancerous tissue were used to simulate plasma exposure of drugs with different extraction ratios. The PBPK simulations revealed a substantial difference in drug exposure (area under the curve), up to 3.3-fold, when using typical scaling factors (healthy population) instead of disease-related parameters in cancer population. These indicate the importance of using population-specific scalars in IVIVE-PBPK for different disease states. SIGNIFICANCE STATEMENT: Accuracy in predicting the fate of drugs from in vitro data using IVIVE-PBPK depends on using correct scaling factors. The values for two of such scalars, namely microsomal and cytosolic protein per gram of liver, is not known in patients with cancer. This study presents, for the first time, scaling factors from cancerous and matched histologically normal livers. PBPK simulations of various metabolically cleared drugs demonstrate the necessity of population-specific scaling for model-informed precision dosing in oncology.

Influence of Proteome Profiles and Intracellular Drug Exposure on Differences in CYP Activity in Donor-Matched Human Liver Microsomes and Hepatocytes

Human liver microsomes (HLM) and human hepatocytes (HH) are important in vitro systems for studies of intrinsic drug clearance (CLint) in the liver. However, the CLint values are often in disagreement for these two systems. Here, we investigated these differences in a side-by-side comparison of drug metabolism in HLM and HH prepared from 15 matched donors. Protein expression and intracellular unbound drug concentration (Kpuu) effects on the CLint were investigated for five prototypical probe substrates (bupropion-CYP2B6, diclofenac-CYP2C9, omeprazole-CYP2C19, bufuralol-CYP2D6, and midazolam-CYP3A4). The samples were donor-matched to compensate for inter-individual variability but still showed systematic differences in CLint. Global proteomics analysis outlined differences in HLM from HH and homogenates of human liver (HL), indicating variable enrichment of ER-localized cytochrome P450 (CYP) enzymes in the HLM preparation. This suggests that the HLM may not equally and accurately capture metabolic capacity for all CYPs. Scaling CLint with CYP amounts and Kpuu could only partly explain the discordance in absolute values of CLint for the five substrates. Nevertheless, scaling with CYP amounts improved the agreement in rank order for the majority of the substrates. Other factors, such as contribution of additional enzymes and variability in the proportions of active and inactive CYP enzymes in HLM and HH, may have to be considered to avoid the use of empirical scaling factors for prediction of drug metabolism.

Evaluation of Quantitative Structure Property Relationship Algorithms for Predicting Plasma Protein Binding in Humans

The extent of plasma protein binding is an important compound-specific property that influences a compound's pharmacokinetic behavior and is a critical input parameter for predicting exposure in physiologically based pharmacokinetic (PBPK) modeling. When experimentally determined fraction unbound in plasma (fup) data are not available, quantitative structure-property relationship (QSPR) models can be used for prediction. Because available QSPR models were developed based on training sets containing pharmaceutical-like compounds, we compared their prediction accuracy for environmentally relevant and pharmaceutical compounds. Fup values were calculated using Ingle et al., Watanabe et al. and ADMET Predictor (Simulation Plus). The test set included 818 pharmaceutical and environmentally relevant compounds with fup values ranging from 0.01 to 1. Overall, the three QSPR models resulted in over-prediction of fup for highly binding compounds and under-prediction for low or moderately binding compounds. For highly binding compounds (0.01≤ fup ≤ 0.25), Watanabe et al. performed better with a lower mean absolute error (MAE) of 6.7% and a lower mean absolute relative prediction error (RPE) of 171.7 % than other methods. For low to moderately binding compounds, both Ingle et al. and ADMET Predictor performed better than Watanabe et al. with superior MAE and RPE values. The positive polar surface area, the number of basic functional groups and lipophilicity were the most important chemical descriptors for predicting fup. This study demonstrated that the prediction of fup was the most uncertain for highly binding compounds. This suggested that QSPR-predicted fup values should be used with caution in PBPK modeling.

Prediction of Drug Clearance from Enzyme and Transporter Kinetics

Accurate estimation of in vivo clearance in human is pivotal to determine the dose and dosing regimen for drug development. In vitro-in vivo extrapolation (IVIVE) has been performed to predict drug clearance using empirical and physiological scalars. Multiple in vitro systems and mathematical modeling techniques have been employed to estimate in vivo clearance. The models for predicting clearance have significantly improved and have evolved to become more complex by integrating multiple processes such as drug metabolism and transport as well as passive diffusion. This chapter covers the use of conventional as well as recently developed methods to predict metabolic and transporter-mediated clearance along with the advantages and disadvantages of using these methods and the associated experimental considerations. The general approaches to improve IVIVE by use of appropriate scalars, incorporation of extrahepatic metabolism and transport and application of physiologically based pharmacokinetic (PBPK) models with proteomics data are also discussed. The chapter also provides an overview of the advantages of using such dynamic mechanistic models over static models for clearance predictions to improve IVIVE.

Parameterization of Microsomal and Cytosolic Scaling Factors: Methodological and Biological Considerations for Scalar Derivation and Validation

Mathematical models that can predict the kinetics of compounds have been increasingly adopted for drug development and risk assessment. Data for these models may be generated from in vitro experimental systems containing enzymes contributing to metabolic clearance, such as subcellular tissue fractions including microsomes and cytosol. Extrapolation from these systems is facilitated by common scaling factors, known as microsomal protein per gram (MPPG) and cytosolic protein per gram (CPPG). Historically, parameterization of MPPG and CPPG has employed the use of recovery factors, commonly benchmarked to cytochromes P450 which work well in some contexts, but could be problematic for other enzymes. Here, we propose absolute quantification of protein content and supplementary assays to evaluate microsomal/cytosolic purity that should be employed. Examples include calculation of microsomal latency by mannose-6-phosphatase activity and immunoblotting of subcellular fractions with fraction-specific markers. Further considerations include tissue source, as disease states can affect enzyme expression and activity, and the methodology used for scalar parameterization. Regional- and organ-specific expression of enzymes, in addition to differences in organ physiology, is another important consideration. Because most efforts have focused on the liver that is, for the most part, homogeneous, derived scalars may not capture the heterogeneity of other major tissues contributing to xenobiotic metabolism including the kidneys and small intestine. Better understanding of these scalars, and how to appropriately derive them from extrahepatic tissues can provide support to the inferences made with physiologically based pharmacokinetic modeling, increase its accuracy in characterizing in vivo drug pharmacokinetics, and improve confidence in go-no-go decisions for clinical trials.

Metabolism and plasma protein binding of 16 straight- and branched-chain parabens in in vitro liver and skin models

Parabens are alkyl esters of 4-hydroxybenzoic acid (4-HBA), with short-chain parabens used as antimicrobials in cosmetics. We investigated the impact of chain structure on skin and liver metabolism. Incubations with primary human hepatocytes and human liver S9 indicated that methyl-, ethyl-, propyl- and butylparaben were rapidly metabolized to similar metabolites, including 4-HBA plus the corresponding alcohols. Liver and EpiSkin™ S9 were used to investigate the metabolism of 16 short and long straight- and branched-chain parabens. The rate of hydrolysis generally decreased with increasing chain length in liver S9, whereas the reverse was true for EpiSkin™ S9. Chain length also correlated with the number of metabolites, with more oxidized metabolites detected from longer chain parabens. The identity of the alcohol group impacted metabolism the most, in terms of the rate of metabolism and the contribution of cofactors. The majority of parabens (13/16) exhibited high plasma protein binding (PPB) (>90%); whereas, 4-HBA PPB was 38%. PPB was related to the LogP of the parabens. In conclusion, the major and common paraben metabolite in PHH, liver S9 and EpiSkin™ S9 was 4-HBA. The rate of metabolism, type of metabolite and contribution of hydrolysis was tissue-specific (liver, skin) and was influenced by the chain length (and hence LogP), structural isomeric form (straight vs branched), and/or the identity of the alkyl group. SHORT ABSTRACT: We investigated how the chain structure of parabens affects their metabolism by liver and EpiSkin™ S9. The major and common metabolite in primary human hepatocytes, liver S9 and EpiSkin™ S9 was 4-HBA plus the corresponding alcohols. The rate of metabolism, type of metabolite and contribution of hydrolysis was tissue-specific and influenced by the chain length, structural isomeric form (straight vs branched), and/or the identity of the alkyl group. Most parabens exhibited high PPB (>90%), whereas the PPB of 4-HBA was 38%.

Using in vitro derived enzymatic reaction rates of metabolism to inform pesticide body burdens in amphibians

Understanding how pesticide exposure to non-target species influences toxicity is necessary to accurately assess the ecological risks these compounds pose. To assess the potential metabolic activation of broad use pesticides in amphibians, in vitro and in vivo metabolic rate constants were derived from toad (Anaxyrus terrestris) livers in experiments measuring the depletion of atrazine (ATZ), triadimefon (TDN), and fipronil (FIP) as well as formation of their metabolites. To determine the predictability of these in vitro derived rate constants, Fowler's toads (Anaxyrus fowleri) were exposed to soil contaminated with each of the pesticides at maximum application rate. Desethyl atrazine (DEA) and deisopropyl atrazine (DIA), both metabolites of ATZ, exhibited similar velocities (V_max) while the K_M constant for DIA was two times higher than DEA. TDN was metabolized into two diastereomers of triadimenol (TDL A and TDL B), where TDL B had a V_max around two times higher than TDL A. The metabolite fipronil sulfone's V_max and K_M were 150 pmol min^-1 mg^-1 and 29 μM, respectively. While intrinsic clearance rates for the pesticides ranged from 0.54 to 38.31 mL min^-1 kg^-1. Thus, gaining knowledge on differences in metabolism of pesticides within amphibians is important in estimating risk to these non-target species since the inherent toxicity of metabolites can differ from the parent compound.

In vitro to in vivo extrapolation for high throughput prioritization and decision making

In vitro chemical safety testing methods offer the potential for efficient and economical tools to provide relevant assessments of human health risk. To realize this potential, methods are needed to relate in vitro effects to in vivo responses, i.e., in vitro to in vivo extrapolation (IVIVE). Currently available IVIVE approaches need to be refined before they can be utilized for regulatory decision-making. To explore the capabilities and limitations of IVIVE within this context, the U.S. Environmental Protection Agency Office of Research and Development and the National Toxicology Program Interagency Center for the Evaluation of Alternative Toxicological Methods co-organized a workshop and webinar series. Here, we integrate content from the webinars and workshop to discuss activities and resources that would promote inclusion of IVIVE in regulatory decision-making. We discuss properties of models that successfully generate predictions of in vivo doses from effective in vitro concentration, including the experimental systems that provide input parameters for these models, areas of success, and areas for improvement to reduce model uncertainty. Finally, we provide case studies on the uses of IVIVE in safety assessments, which highlight the respective differences, information requirements, and outcomes across various approaches when applied for decision-making.

A framework and case studies for evaluation of enzyme ontogeny in children's health risk evaluation

Knowledge of the ontogeny of Phase I and Phase II metabolizing enzymes may be used to inform children's vulnerability based upon likely differences in internal dose from xenobiotic exposure. This might provide a qualitative assessment of toxicokinetic (TK) variability and uncertainty pertinent to early lifestages and help scope a more quantitative physiologically based toxicokinetic (PBTK) assessment. Although much is known regarding the ontogeny of metabolizing systems, this is not commonly utilized in scoping and problem formulation stage of human health risk evaluation. A framework is proposed for introducing this information into problem formulation which combines data on enzyme ontogeny and chemical-specific TK to explore potential child/adult differences in internal dose and whether such metabolic differences may be important factors in risk evaluation. The framework is illustrated with five case study chemicals, including some which are data rich and provide proof of concept, while others are data poor. Case studies for toluene and chlorpyrifos indicate potentially important child/adult TK differences while scoping for acetaminophen suggests enzyme ontogeny is unlikely to increase early-life risks. Scoping for trichloroethylene and aromatic amines indicates numerous ways that enzyme ontogeny may affect internal dose which necessitates further evaluation. PBTK modeling is a critical and feasible next step to further evaluate child-adult differences in internal dose for a number of these chemicals.

Content and activity of human liver microsomal protein and prediction of individual hepatic clearance in vivo

The lack of information concerning individual variation in content and activity of human liver microsomal protein is one of the most important obstacles for designing personalized medicines. We demonstrated that the mean value of microsomal protein per gram of liver (MPPGL) was 39.46 mg/g in 128 human livers and up to 19-fold individual variations existed. Meanwhile, the metabolic activities of 10 cytochrome P450 (CYPs) were detected in microsomes and liver tissues, respectively, which showed huge individual variations (200-fold). Compared with microsomes, the activities of liver tissues were much suitable to express the individual variations of CYP activities. Furthermore, individual variations in the in vivo clearance of tolbutamide were successfully predicted with the individual parameter values. In conclusion, we offer the values for MPPGL contents in normal liver tissues and build a new method to assess the in vitro CYP activities. In addition, large individual variations exist in predicted hepatic clearance of tolbutamide. These findings provide important physiological parameters for physiologically-based pharmacokinetics models and thus, establish a solid foundation for future development of personalized medicines.

Trichloroethylene biotransformation and its role in mutagenicity, carcinogenicity and target organ toxicity

Metabolism is critical for the mutagenicity, carcinogenicity, and other adverse health effects of trichloroethylene (TCE). Despite the relatively small size and simple chemical structure of TCE, its metabolism is quite complex, yielding multiple intermediates and end-products. Experimental animal and human data indicate that TCE metabolism occurs through two major pathways: cytochrome P450 (CYP)-dependent oxidation and glutathione (GSH) conjugation catalyzed by GSH S-transferases (GSTs). Herein we review recent data characterizing TCE processing and flux through these pathways. We describe the catalytic enzymes, their regulation and tissue localization, as well as the evidence for transport and inter-organ processing of metabolites. We address the chemical reactivity of TCE metabolites, highlighting data on mutagenicity of these end-products. Identification in urine of key metabolites, particularly trichloroacetate (TCA), dichloroacetate (DCA), trichloroethanol and its glucuronide (TCOH and TCOG), and N-acetyl-S-(1,2-dichlorovinyl)-L-cysteine (NAcDCVC), in exposed humans and other species (mostly rats and mice) demonstrates function of the two metabolic pathways in vivo. The CYP pathway primarily yields chemically stable end-products. However, the GST pathway conjugate S-(1,2-dichlorovinyl)glutathione (DCVG) is further processed to multiple highly reactive species that are known to be mutagenic, especially in kidney where in situ metabolism occurs. TCE metabolism is highly variable across sexes, species, tissues and individuals. Genetic polymorphisms in several of the key enzymes metabolizing TCE and its intermediates contribute to variability in metabolic profiles and rates. In all, the evidence characterizing the complex metabolism of TCE can inform predictions of adverse responses including mutagenesis, carcinogenesis, and acute and chronic organ-specific toxicity.

Modulation of trichloroethylene in vitro metabolism by different drugs in human

Toxicological interactions with drugs have the potential to modulate the toxicity of trichloroethylene (TCE). Our objective is to identify metabolic interactions between TCE and 14 widely used drugs in human suspended hepatocytes and characterize the strongest using microsomal assays. Changes in concentrations of TCE and its metabolites were measured by headspace GC-MS. Results with hepatocytes show that amoxicillin, cimetidine, ibuprofen, mefenamic acid and ranitidine caused no significant interactions. Naproxen and salicylic acid showed to increase both TCE metabolites levels, whereas acetaminophen, carbamazepine and erythromycin rather decreased them. Finally, diclofenac, gliclazide, sulphasalazine and valproic acid had an impact on the levels of only one metabolite. Among the 14 tested drugs, 5 presented the most potent interactions and were selected for confirmation with microsomes, namely naproxen, salicylic acid, acetaminophen, carbamazepine and valproic acid. Characterization in human microsomes confirmed interaction with naproxen by competitively inhibiting trichloroethanol (TCOH) glucuronidation (Ki=2.329 mM). Inhibition of TCOH formation was also confirmed for carbamazepine (partial non-competitive with Ki=70 μM). Interactions with human microsomes were not observed with salicylic acid and acetaminophen, similar to prior results in rat material. For valproic acid, interactions with microsomes were observed in rat but not in human. Inhibition patterns were shown to be similar in human and rat hepatocytes, but some differences in mechanisms were noted in microsomal material between species. Next research efforts will focus on determining the adequacy between in vitro observations and the in vivo situation.

Evaluation of simple in vitro to in vivo extrapolation approaches for environmental compounds

As part of an effort to support in silico/in vitro based risk assessment, we evaluated the accuracy associated with conducting simple in vitro to in vivo extrapolation (IVIVE) for environmental compounds using available in vitro human metabolism data. The IVIVE approach was applied to a number of compounds with a wide range of properties spanning the diversity of characteristics of environmental compounds, and where possible the resulting estimates of the in vivo steady-state blood concentration were compared with estimates derived on the basis of human in vivo kinetic data. There appears to be a systematic bias in the estimation of intrinsic clearance (Clint) from in vitro versus in vivo data, with in vitro based estimates underestimating in vivo clearance for small values of Clint but with the opposite relationship at large values of Clint. Nevertheless, the resulting estimates of Css were in good agreement. The chief drawback of the simple approach used in this study, which performs the IVIVE prediction for the parent compound only, is that it is not applicable for toxicity associated with a metabolite.

Cutting Edge PBPK Models and Analyses: Providing the Basis for Future Modeling Efforts and Bridges to Emerging Toxicology Paradigms

Physiologically based Pharmacokinetic (PBPK) models are used for predictions of internal or target dose from environmental and pharmacologic chemical exposures. Their use in human risk assessment is dependent on the nature of databases (animal or human) used to develop and test them, and includes extrapolations across species, experimental paradigms, and determination of variability of response within human populations. Integration of state-of-the science PBPK modeling with emerging computational toxicology models is critical for extrapolation between in vitro exposures, in vivo physiologic exposure, whole organism responses, and long-term health outcomes. This special issue contains papers that can provide the basis for future modeling efforts and provide bridges to emerging toxicology paradigms. In this overview paper, we present an overview of the field and introduction for these papers that includes discussions of model development, best practices, risk-assessment applications of PBPK models, and limitations and bridges of modeling approaches for future applications. Specifically, issues addressed include: (a) increased understanding of human variability of pharmacokinetics and pharmacodynamics in the population, (b) exploration of mode of action hypotheses (MOA), (c) application of biological modeling in the risk assessment of individual chemicals and chemical mixtures, and (d) identification and discussion of uncertainties in the modeling process.

Physiologically-based pharmacokinetic (PBPK) models in toxicity testing and risk assessment

Physiologically-based pharmacokinetic (PBPK) modeling offers a scientifically-sound framework for integrating mechanistic data on absorption, distribution, metabolism and elimination to predict the time-course of parent chemical, metabolite(s) or biomarkers in the exposed organism. A major advantage of PBPK models is their ability to forecast the impact of specific mechanistic processes and determinants on the tissue dose. In this regard, they facilitate integration of data obtained with in vitro and in silico methods, for making predictions of the tissue dosimetry in the whole animal, thus reducing and/or refining the use of animals in pharmacokinetic and toxicity studies. This chapter presents the principles and practice of PBPK modeling, as well as the application of these models in toxicity testing and health risk assessments.

Parameters for pyrethroid insecticide QSAR and PBPK/PD models for human risk assessment

In this review we have examined the status of parameters required by pyrethroid QSAR-PBPK/PD models for assessing health risks. In lieu of the chemical,biological, biochemical, and toxicological information developed on the pyrethroids since 1968, the finding of suitable parameters for QSAR and PBPK/PD model development was a monumental task. The most useful information obtained came from rat toxicokinetic studies (i.e., absorption, distribution, and excretion), metabolism studies with 14C-cyclopropane- and alcohol-labeled pyrethroids, the use of known chiral isomers in the metabolism studies and their relation to commercial products. In this review we identify the individual chiralisomers that have been used in published studies and the chiral HPLC columns available for separating them. Chiral HPLC columns are necessary for isomer identification and for developing kinetic values (Vm,, and Kin) for pyrethroid hydroxylation. Early investigators synthesized analytical standards for key pyrethroid metabolites, and these were used to confirm the identity of urinary etabolites, by using TLC. These analytical standards no longer exist, and muste resynthesized if further studies on the kinetics of the metabolism of pyrethroids are to be undertaken.In an attempt to circumvent the availability of analytical standards, several CYP450 studies were carried out using the substrate depletion method. This approach does not provide information on the products formed downstream, and may be of limited use in developing human environmental exposure PBPK/PD models that require extensive urinary metabolite data. Hydrolytic standards (i.e., alcohols and acids) were available to investigators who studied the carboxylesterase-catalyzed hydrolysis of several pyrethroid insecticides. The data generated in these studies are suitable for use in developing human exposure PBPK/PD models.Tissue:blood partition coefficients were developed for the parent pyrethroids and their metabolites, by using a published mechanistic model introduced by Poulin and Thiele (2002a; b) and log DpH 7.4 values. The estimated coefficients, especially those of adipose tissue, were too high and had to be corrected by using a procedure in which the proportion of parent or metabolite residues that are unbound to plasma albumin is considered, as described in the GastroPlus model (Simulations Plus, Inc.,Lancaster, CA). The literature suggested that Km values be adjusted by multiplying Km by the substrate (decimal amount) that is unbound to microsomal or CYPprotein. Mirfazaelian et al. (2006) used flow- and diffusion-limited compartments in their deltamethrin model. The addition of permeability areas (PA) having diffusion limits, such as the fat and slowly perfused compartments, enabled the investigators to bring model predictions in line with in vivo data.There appears to be large differences in the manner and rate of absorption of the pyrethroids from the gastrointestinal tract, implying that GI advanced compartmental transit models (ACAT) need to be included in PBPK models. This is especially true of the absorption of an oral dose of tefluthrin in male rats, in which 3.0-6.9%,41.3-46.3%, and 5.2-15.5% of the dose is eliminated in urine, feces, and bile,respectively (0-48 h after administration). Several percutaneous studies with the pyrethroids strongly support the belief that these insecticides are not readily absorbed, but remain on the surface of the skin until they are washed off. In one articular study (Sidon et al. 1988) the high levels of permethrin absorption through the forehead skin (24-28%) of the monkey was reported over a 7- to 14-days period.Wester et al. (1994) reported an absorption of 1.9% of pyrethrin that had been applied to the forearm of human volunteers over a 7-days period.SAR models capable of predicting the binding of the pyrethroids to plasma and hepatic proteins were developed by Yamazaki and Kanaoka (2004), Saiakhov et al. (2000), Colmenarejo et al. (2001), and Colmenarejo (2003). QikProp(Schrodinger, LLC) was used to obtain Fu values for calculating partition coefficients and for calculating permeation constants (Caco-2, MDCK, and logBBB). ADMET Predictor (Simulations Plus Inc.) provided Vm~,x and Km values for the hydroxylation of drugs/pyrethroids by human liver recombinant cytochrome P450 enzymes making the values available for possible use in PBPK/PD models.The Caco-2 permeability constants and CYP3A4 Vmax and Km values are needed in PBPK/PD models with GI ACAT sub models. Modeling work by Chang et al.(2009) produced rate constants (kcat) for the hydrolysis of pyrethroids by rat serumcarboxylesterases. The skin permeation model of Potts and Guy (1992) was used topredict K, values for the dermal absorption of the 15 pyrethroids.The electrophysiological studies by Narahashi (1971) and others (Breckenridgeet al. 2009; Shafer et al. 2005; Soderlund et al. 2002; Wolansky and Harrill 2008)demonstrated that the mode of action of pyrethroids on nerves is to interfere with the changes in sodium and potassium ion currents. The pyrethroids, being highly lipid soluble, are bound or distributed in lipid bilayers of the nerve cell membrane and exert their action on sodium channel proteins. The rising phase of the action potential is caused by sodium influx (sodium activation), while the falling phase is caused by sodium activation being turned off, and an increase in potassium efflux(potassium activation). The action of allethrin and other pyrethroids is caused by an inhibition or block of the normal currents. An equation by Tatebayashi and Narahashi (1994) that describes the action of pyrethroids on sodium channels was found in the literature. This equation, or some variation of it, may be suitable for use in the PD portion of pyrethroid PBPK models.

Comparison and evaluation of pharmacokinetics of PFOA and PFOS in the adult rat using a physiologically based pharmacokinetic model

Perfluoroalkyl acid carboxylates and sulfonates (PFAAs) have many consumer and industrial applications. The persistence and widespread distribution of PFAAs have brought them under intense scrutiny. Limited PK data for PFAAs is available for humans; however, toxicological and pharmacokinetic data exist for rats, which can be useful for cross-species extrapolation. In this work, PBPK models were developed for adult male and female rats to describe the pharmacokinetics of PFOA and PFOS. The models contain a description of saturable renal resorption, free fraction of chemical in plasma, and saturable binding in liver. Both male and female rat models for each chemical were consistent with available PK data resulting from IV, oral, and dietary dosing regimens. Predicted plasma concentration curves followed trends observed in experimental data, and model predictions were within a factor of two of experimental values. PFOA and PFOS rat model output is sensitive to parameters governing renal resorption, indicating that renal resorption is responsible for the long-half life. These models, along with the PFAA gestation and lactation models published in this issue, will help address concerns about possible health effects due to PFAA exposure in the fetus and neonate and will be useful in comparing PK across life stages.

Cytochrome P450-specific human PBPK/PD models for the organophosphorus pesticides: chlorpyrifos and parathion

Organophosphorus pesticides (OPs) remain a potential concern to human health because of their continuing use worldwide. Phosphororthioate OPs like chlorpyrifos and parathion are directly activated and detoxified by various cytochrome P450s (CYPs), with the primary CYPs involved being CYP2B6 and CYP2C19. The goal of the current study was to convert a previously reported human pharmacokinetic and pharmacodynamic (PBPK/PD) model for chlorpyrifos, that used chlorpyrifos metabolism parameters from rat liver, into a human CYP based/age-specific model using recombinant human CYP kinetic parameters (V(max), K(m)), hepatic CYP content and plasma binding measurements to estimate new values for acetylcholinesterase (AChE) and butyrylcholinesterase (BuChE) inhibition and to use the model as a template for the development of a comparable parathion PBPK/PD model. The human CYP based/age-specific PBPK/PD models were used to simulate single oral exposures of adults (19 year old) and infants (1 year) to chlorpyrifos (10,000, 1000 and 100 μg/kg) or parathion (100, 25 and 5 μg/kg). Model simulations showed that there is an age dependency in the amount of blood cholinesterase inhibition observed, however additional age-dependent data are needed to further optimize age-specific human PBPK/PD modeling for these OP compounds. PBPK/PD model simulations estimated that a 4-fold increase or decrease in relative CYP2B6 and CYP2C19 content would produce a 9-22% inhibition in blood AChE activity following exposure of an adult to chlorpyrifos (1000 μg/kg). Similar model simulation produced an 18-22% inhibition in blood AChE activity following exposure of an adult to parathion (25 μg/kg). Individuals with greater CYP2B6 content and lower CYP2C19 content were predicted to be most sensitive to both OPs. Changes in hepatic CYP2B6 and CYP2C19 content had more of an influence on cholinesterase inhibition for exposures to chlorpyrifos than parathion, which agrees with previously reported literature that these CYPs are more reaction biased for desulfurization (activation) and dearylation (detoxification) of chlorpyrifos compared to parathion. The data presented here illustrate how PBPK/PD models with human enzyme-specific parameters can assist ongoing risk assessment efforts and aid in the identification of sensitive individuals and populations.

Development and evaluation of a harmonized physiologically based pharmacokinetic (PBPK) model for perchloroethylene toxicokinetics in mice, rats, and humans

This article reports on the development of a "harmonized" PBPK model for the toxicokinetics of perchloroethylene (tetrachloroethylene or perc) in mice, rats, and humans that includes both oxidation and glutathione (GSH) conjugation of perc, the internal kinetics of the oxidative metabolite trichloroacetic acid (TCA), and the urinary excretion kinetics of the GSH conjugation metabolites N-Acetylated trichlorovinyl cysteine and dichloroacetic acid. The model utilizes a wider range of in vitro and in vivo data than any previous analysis alone, with in vitro data used for initial, or "baseline," parameter estimates, and in vivo datasets separated into those used for "calibration" and those used for "evaluation." Parameter calibration utilizes a limited Bayesian analysis involving flat priors and making inferences only using posterior modes obtained via Markov chain Monte Carlo (MCMC). As expected, the major route of elimination of absorbed perc is predicted to be exhalation as parent compound, with metabolism accounting for less than 20% of intake except in the case of mice exposed orally, in which metabolism is predicted to be slightly over 50% at lower exposures. In all three species, the concentration of perc in blood, the extent of perc oxidation, and the amount of TCA production is well-estimated, with residual uncertainties of ~2-fold. However, the resulting range of estimates for the amount of GSH conjugation is quite wide in humans (~3000-fold) and mice (~60-fold). While even high-end estimates of GSH conjugation in mice are lower than estimates of oxidation, in humans the estimated rates range from much lower to much higher than rates for perc oxidation. It is unclear to what extent this range reflects uncertainty, variability, or a combination. Importantly, by separating total perc metabolism into separate oxidative and conjugative pathways, an approach also recommended in a recent National Research Council review, this analysis reconciles the disparity between those previously published PBPK models that concluded low perc metabolism in humans and those that predicted high perc metabolism in humans. In essence, both conclusions are consistent with the data if augmented with some additional qualifications: in humans, oxidative metabolism is low, while GSH conjugation metabolism may be high or low, with uncertainty and/or interindividual variability spanning three orders of magnitude. More direct data on the internal kinetics of perc GSH conjugation, such as trichlorovinyl glutathione or tricholorvinyl cysteine in blood and/or tissues, would be needed to better characterize the uncertainty and variability in GSH conjugation in humans.

Prediction of modified release pharmacokinetics and pharmacodynamics from in vitro, immediate release, and intravenous data

The aim of this study was to demonstrate the value of mechanistic simulations in gaining insight into the behaviors of modified release (MR) formulations in vivo and to use the properly calibrated models for prediction of pharmacokinetics (PK) and pharmacodynamics (PD). GastroPlus (Simulations Plus, Inc.) was used to fit mechanistic models for adinazolam and metoprolol that describe the absorption, PK, and PD after intravenous (i.v.) and immediate release (IR) oral (p.o.) administration. The fitted model for adinazolam was then used to predict the PD profile for a MR formulation and to design a new formulation with desired onset and duration of action. The fitted metoprolol model was used to gain insight and to explain the in vivo behaviors of MR formulations. For each drug, a single absorption/PK model was fitted that provided simulated plasma concentration-time profiles closely matching observed in vivo profiles across several different i.v. and p.o doses. Sedation score profiles of adinazolam were fitted with an indirect PD model. For metoprolol, the fitted absorption/PK model for IR p.o. doses was used to select in vitro dissolution conditions that best matched the in vivo release of MR doses. This model also explained differences in exposure after administration of MR formulations with different release rates. Mechanistic absorption/PK models allow for detailed descriptions of all processes affecting the two drugs' bioavailability, including release/dissolution, absorption, and intestinal and hepatic first pass extraction. The insights gained can be used to design formulations that more effectively overcome identified problems.

Determination of microsome and hepatocyte scaling factors for in vitro/in vivo extrapolation in the rat and dog

1. In vivo clearance predictions from in vitro assays require scaling factors to relate the concentrations of hepatocytes or microsomal protein to the intact liver. 2. The aims were to measure the variability in scaling factors for Wistar rat and beagle dog for which the literature is particularly scarce and determine any sex differences. 3. Scaling factors were determined by comparing the cytochrome P450 (P450) content in hepatocytes or microsomes against the P450 content of fresh liver homogenate. The use of fresh homogenate is recommended as freezing can increase contamination and affect the P450 assay. 4. Mean(geo) hepatic microsomal concentrations in Wistar rats were 61 mg g(-1) liver (95% confidence interval (CI); 47-75 mg g(-1) liver) and in beagle dogs 55 mg g(-1) liver (95% CI = 48-62 mg g(-1) liver). Mean(geo) hepatocellularity was 163 x 10(6) cells g(-1) liver for Wistar rats (95% CI = 127-199 x 10(6) cells g(-1) liver) and 169 x 10(6) cells g(-1) liver (95% CI = 131-207 x 10(6) cells g(-1) liver) for beagle dogs. The data generated in this study indicate a consistency in scaling factors between rat and dog. No sex differences were observed.

Approaches for applications of physiologically based pharmacokinetic models in risk assessment

Physiologically based pharmacokinetic (PBPK) models are particularly useful for simulating exposures to environmental toxicants for which, unlike pharmaceuticals, there is often little or no human data available to estimate the internal dose of a putative toxic moiety in a target tissue or an appropriate surrogate. This article reviews the current state of knowledge and approaches for application of PBPK models in the process of deriving reference dose, reference concentration, and cancer risk estimates. Examples drawn from previous U.S. Environmental Protection Agency (EPA) risk assessments and human health risk assessments in peer-reviewed literature illustrate the ways and means of using PBPK models to quantify the pharmacokinetic component of the interspecies and intraspecies uncertainty factors as well as to conduct route to route, high dose to low dose and duration extrapolations. The choice of the appropriate dose metric is key to the use of the PBPK models for the various applications in risk assessment. Issues related to whether uncertainty factors are most appropriately applied before or after derivation of human equivalent dose (or concentration) continue to be explored. Scientific progress in the understanding of life stage and genetic differences in dosimetry and their impacts on variability in susceptibility, as well as ongoing development of analytical methods to characterize uncertainty in PBPK models, will make their use in risk assessment increasingly likely. As such, it is anticipated that when PBPK models are used to express adverse tissue responses in terms of the internal target tissue dose of the toxic moiety rather than the external concentration, the scientific basis of, and confidence in, risk assessments will be enhanced.

Enantiomeric Metabolic Interactions and Stereoselective Human Methadone Metabolism

Methadone is administered as a racemate, although opioid activity resides in the R-enantiomer. Methadone disposition is stereoselective, with considerable unexplained variability in clearance and plasma R/S ratios. N-Demethylation of methadone in vitro is predominantly mediated by cytochrome P450 CYP3A4 and CYP2B6 and somewhat by CYP2C19. This investigation evaluated stereoselectivity, models, and kinetic parameters for methadone N-demethylation by recombinant CYP2B6, CYP3A4, and CYP2C19, and the potential for interactions between enantiomers during racemate metabolism. CYP2B6 metabolism was stereoselective. CYP2C19 was less active, and stereoselectivity was opposite that for CYP2B6. CYP3A4 was not stereoselective. With all three isoforms, enantiomer N-dealkylation rates in the racemate were lower than those of (R)-(6-dimethyamino-4,4-diphenyl-heptan-3-one) hydrochloride (R-methadone) or (S)-(6-dimethyamino-4,4-diphenyl-heptan-3-one) hydrochloride (S-methadone) alone, suggesting an enantiomeric interaction and mutual metabolic inhibition. For CYP2B6, the interaction between enantiomers was stereoselective, with S-methadone as a more potent inhibitor of R-methadone N-demethylation than R-of S-methadone. In contrast, enantiomer interactions were not stereoselective with CYP2C19 or CYP3A4. For all three cytochromes P450, methadone N-demethylation was best described by two-site enzyme models with competitive inhibition. There were minor model differences between cytochromes P450 to account for stereoselectivity of metabolism and enantiomeric interactions. Changes in plasma R/S methadone ratios observed after rifampin or troleandomycin pretreatment in humans in vivo were successfully predicted by CYP2B6- but not CYP3A4-catalyzed methadone N-demethylation. CYP2B6 is a predominant catalyst of stereoselective methadone metabolism in vitro. In vivo, CYP2B6 may be a major determinant of methadone metabolism and disposition, and CYP2B6 activity and stereoselective metabolic interactions may confer variability in methadone disposition.

Evaluation of physiologically based pharmacokinetic models for use in risk assessment

Physiologically based pharmacokinetic (PBPK) models are sophisticated dosimetry models that offer great flexibility in modeling exposure scenarios for which there are limited data. This is particularly of relevance to assessing human exposure to environmental toxicants, which often requires a number of extrapolations across species, route, or dose levels. The continued development of PBPK models ensures that regulatory agencies will increasingly experience the need to evaluate available models for their application in risk assessment. To date, there are few published criteria or well-defined standards for evaluating these models. Herein, important considerations for evaluating such models are described. The evaluation of PBPK models intended for risk assessment applications should include a consideration of: model purpose, model structure, mathematical representation, parameter estimation, computer implementation, predictive capacity and statistical analyses. Model purpose and structure require qualitative checks on the biological plausibility of a model. Mathematical representation, parameter estimation, computer implementation involve an assessment of the coding of the model, as well as the selection and justification of the physical, physicochemical and biochemical parameters chosen to represent a biological organism. Finally, the predictive capacity and sensitivity, variability and uncertainty of the model are analysed so that the applicability of a model for risk assessment can be determined. Published in 2007 by John Wiley & Sons, Ltd.

Toxicokinetics of Inhaled Trichloroethylene and Tetrachloroethylene in Humans at 1 ppm: Empirical Results and Comparisons with Previous Studies

Trichloroethylene (TRI) and tetrachloroethylene (TETRA) are solvents that have been widely used in a variety of industries, and both are widespread environmental contaminants. In order to provide a better basis for understanding their toxicokinetics at environmental exposures, seven human volunteers were exposed by inhalation to 1 ppm of TRI or TETRA for 6 h, with biological samples collected for analysis during exposure and up to 6-days postexposure. Concentrations of TRI, TETRA, free trichloroethanol (TCOH), total TCOH (free TCOH plus glucuronidated TCOH), and trichloroacetic acid (TCA) were determined in blood and urine; TRI and TETRA concentrations were measured in alveolar breath. Toxicokinetic time courses and empirical analyses of classical toxicokinetic parameters were compared with those reported in previous human volunteer studies, most of which involved exposures that were at least 10-fold higher. Qualitatively, TRI and TETRA toxicokinetics were consistent with previous human studies. Quantitatively, alveolar retention and clearance by exhalation were similar to those found previously but blood and urine data suggest a number of possible toxicokinetic differences. For TRI, data from the current study support lower apparent blood-air partition coefficients, greater apparent metabolic clearance, less TCA production, and greater glucuronidation of TCOH as compared to previous studies. For TETRA, the current data suggest TCA formation that is similar or slightly lower than that of previous studies. Variability and uncertainty in empirical estimates of total TETRA metabolism are substantial, with confidence intervals among different studies substantially overlapping. Relative contributions to observed differences from concentration-dependent toxicokinetics and interindividual and interoccasion variability remain to be determined.

Bayesian population analysis of a harmonized physiologically based pharmacokinetic model of trichloroethylene and its metabolites

Bayesian population analysis of a harmonized physiologically based pharmacokinetic (PBPK) model for trichloroethylene (TCE) and its metabolites was performed. In the Bayesian framework, prior information about the PBPK model parameters is updated using experimental kinetic data to obtain posterior parameter estimates. Experimental kinetic data measured in mice, rats, and humans were available for this analysis, and the resulting posterior model predictions were in better agreement with the kinetic data than prior model predictions. Uncertainty in the prediction of the kinetics of TCE, trichloroacetic acid (TCA), and trichloroethanol (TCOH) was reduced, while the kinetics of other key metabolites dichloroacetic acid (DCA), chloral hydrate (CHL), and dichlorovinyl mercaptan (DCVSH) remain relatively uncertain due to sparse kinetic data for use in this analysis. To help focus future research to further reduce uncertainty in model predictions, a sensitivity analysis was conducted to help identify the parameters that have the greatest impact on various internal dose metric predictions. For application to a risk assessment for TCE, the model provides accurate estimates of TCE, TCA, and TCOH kinetics. This analysis provides an important step toward estimating uncertainty of dose-response relationships in noncancer and cancer risk assessment, improving the extrapolation of toxic TCE doses from experimental animals to humans.

Issues in the pharmacokinetics of trichloroethylene and its metabolites

Much progress has been made in understanding the complex pharmacokinetics of trichloroethylene (TCE) . Qualitatively, it is clear that TCE is metabolized to multiple metabolites either locally or into systemic circulation. Many of these metabolites are thought to have toxicologic importance. In addition, efforts to develop physiologically based pharmacokinetic (PBPK) models have led to a better quantitative assessment of the dosimetry of TCE and several of its metabolites. As part of a mini-monograph on key issues in the health risk assessment of TCE, this article is a review of a number of the current scientific issues in TCE pharmacokinetics and recent PBPK modeling efforts with a focus on literature published since 2000. Particular attention is paid to factors affecting PBPK modeling for application to risk assessment. Recent TCE PBPK modeling efforts, coupled with methodologic advances in characterizing uncertainty and variability, suggest that rigorous application of PBPK modeling to TCE risk assessment appears feasible at least for TCE and its major oxidative metabolites trichloroacetic acid and trichloroethanol. However, a number of basic structural hypotheses such as enterohepatic recirculation, plasma binding, and flow- or diffusion-limited treatment of tissue distribution require additional evaluation and analysis. Moreover, there are a number of metabolites of potential toxicologic interest, such as chloral, dichloroacetic acid, and those derived from glutathione conjugation, for which reliable pharmacokinetic data is sparse because of analytical difficulties or low concentrations in systemic circulation. It will be a challenge to develop reliable dosimetry for such cases.

Metabolism and tissue distribution of orally administered trichloroethylene in male and female rats: identification of glutathione- and cytochrome P-450-derived metabolites in liver, kidney, blood, and urine

Male and female Fischer 344 rats were administered trichloroethylene (TRI) (2, 5, or 15 mmol/kg body weight) in corn oil by oral gavage, and TRI and its metabolites were measured at times up to 48 h in liver, kidneys, blood, and urine. Studies tested the hypothesis that gender-dependent differences in distribution and metabolism of TRI could help explain differences in toxicity. Higher levels of TRI were generally observed in tissues of males at lower doses. Complex patterns of TRI concentration, sometimes with multiple peaks, were observed in liver, kidneys, and blood of both males and females, consistent with enterohepatic recirculation. Higher concentrations of cytochrome P-450 (P450)-derived metabolites were observed in livers of males than in females, whereas the opposite pattern was observed in kidneys. Trichloroacetate was the primary P450-derived metabolite in blood and urine, although it generally appeared at later times than chloral hydrate. Trichloroethanol was also a significant metabolite in urine. S-(1,2-Dichlorovinyl)glutathione (DCVG) was recovered in liver and kidneys of female rats only and in blood of both males and females, with generally higher amounts found in females. S-(1,2-Dichlorovinyl)-L-cysteine (DCVC), the penultimate nephrotoxic metabolite, was recovered in male and female liver, female kidneys, male blood, and in urine of both males and females. The relationship between gender-dependent differences in distribution and metabolism of TRI and susceptibility to TRI-induced toxicity is discussed.

Determination of a Human Hepatic Microsomal Scaling Factor for Predicting in Vivo Drug Clearance

PURPOSE

To determine a microsomal scaling factor for human liver suitable for prediction of in vivo drug clearance from in vitro data and to explore the role of inter-liver variability in this factor on the reported underprediction from microsomal parameters.

METHODS

Cytochrome P450 (henceforth P450) content in whole homogenates and microsomes from 38 donor livers was used to determine a microsomal scaling factor. In a subset (n = 20) of these preparations, individual P450 enzymes were examined by Western blotting and selective probe activities were determined.

RESULTS

The scaling factor from 38 livers averaged 40 mg microsomal protein per gram liver with a coefficient of variation of 31%. Western blotting experiments indicated that there was no P450 enzyme-specific trend in the distribution of individual P450 enzymes in liver microsomes relative to whole homogenate. Predictions based on an average scaling factor resulted in a satisfactory prediction of intrinsic clearance of three benzodiazepines similar to that obtained using individual factors for the same livers.

CONCLUSION

A value for human liver microsomal scaling of 40 mg microsomal protein per gram liver has been established. The reason for underprediction previously reported for 52 different drug substrates was not the use of an incorrect value for the scaling factor.

Under the skin: Biotransformation of para-aminophenol and para-phenylenediamine in reconstructed human epidermis and human hepatocytes

We investigated the biotransformation of the oxidative arylamine (AA) hair dye ingredients [14C]-para-aminophenol (PAP) and [14C]-para-phenylenediamine (PPD) in reconstructed human epidermis and human hepatocytes. Human epidermis quantitatively transformed PAP to its N-acetylated derivative (APAP), whereas hepatocytes transformed PAP to sulfate or glucuronic acid conjugates of APAP or PAP as well as free APAP. Epidermis and hepatocytes converted PPD to N-mono- (MAPPD) and N,N'-di-acetylated (DAPPD) derivatives. At higher concentrations of PPD (250-1000 microM), epidermis or hepatocytes produced more of the MAPPD, whereas concentrations below 250 microM and lower favoured formation of the DAPPD metabolite. When compared with epidermis, human hepatocytes had a three-fold or eight-fold greater capacity for generation of MAPPD or DAPPD, respectively. No evidence of transformation of PAP or PPD to N-hydroxylated derivatives was found in epidermis or hepatocytes. Our results suggest that (i) after dermal absorption of PAP or PPD, humans are systemically exposed to acetylated derivatives; (ii) current in vitro skin absorption studies may be inadapated for determination of human systemic exposure to AAs due to reduced or absent metabolic capacity of non-viable skin; (iii) due to qualitative differences between dermal and hepatic metabolism, oral toxicity studies may be unsuited for the hazard assessment of dermal exposure to AAs; and (iv) use of induced rodent liver S9 metabolic activation systems for in vitro genotoxicity studies may produce misleading results on the hazard of human dermal exposure to AAs. In conclusion, our data support the growing evidence that AAs are transformed in human skin and suggest that current practices of safety assessment of AAs should take these findings into account.

Moving from external exposure concentration to internal dose: duration extrapolation based on physiologically based pharmacokinetic derived estimates of internal dose

The potential human health risk(s) from chemical exposure must frequently be assessed under conditions for which adequate human or animal data are not available. The default method for exposure-duration adjustment, based on Haber's rule, C (external exposure concentration) or C(n) (the ten Berge modification) x t (exposure duration) = K (a constant toxic effect), has been criticized for prediction errors. A promising alternative approach to duration adjustment is based on equivalence of internal dose, that is, target-tissue dose levels, across different exposure durations. A proposed methodology for dose-duration adjustments for acute exposure guideline levels (AEGLs) based on physiologically based pharmacokinetic (PBPK) estimates of dose is illustrated with trichloroethylene (TCE). Steps in this methodology include: (1) selection and evaluation, or development and evaluation, of an appropriate PBPK model; (2) determination of an appropriate measure of internal dose; (3) estimation with the PBPK model of the tissue dose (the target tissue dose) resulting from the external exposure conditions (concentration, duration) of the critical effect; (4) estimation of the external exposure concentrations required to achieve tissue doses equivalent to the target tissue dose at exposure durations of interest; and (5) evaluation of sources of variability and uncertainty. For TCE, this PBPK modeling approach has allowed determination of dose metrics predictive of the acute neurotoxic effects of TCE and dose-duration adjustments based on estimates of internal dose.

Duration adjustment of acute exposure guideline level values for trichloroethylene using a physiologically-based pharmacokinetic model

Acute Exposure Guideline Level (AEGL) recommendations are developed for 10-minute, 30-minute, 1-hour, 4-hours, and 8-hours exposure durations and are designated for three levels of severity: AEGL-1 represents concentrations above which acute exposures may cause noticeable discomfort including irritation; AEGL-2 represents concentrations above which acute exposure may cause irreversible health effects or impaired ability to escape; and AEGL-3 represents concentrations above which exposure may cause life-threatening health effects or death. The default procedure for setting AEGL values across durations when applicable data are unavailable involves estimation based on Haber's rule, which has an underlying assumption that cumulative exposure is the determinant of toxicity. For acute exposure to trichloroethylene (TCE), however, experimental data indicate that momentary tissue concentration, and not the cumulative amount of exposure, is important. We employed an alternative approach to duration adjustments in which a physiologically-based pharmacokinetic (PBPK) model was used to predict the arterial blood concentrations [TCE(a)] associated with adverse outcomes appropriate for AEGL-1, -2, or -3-level effects. The PBPK model was then used to estimate the atmospheric concentration that produces equivalent [TCE(a)] at each of the AEGL-specific exposure durations. This approach yielded [TCE(a)] values of 4.89 mg/l for AEGL-1, 18.7 mg/l for AEGL-2, and 310 mg/l for AEGL-3. Duration adjustments based on equivalent target tissue doses should provide similar degrees of toxicity protection at different exposure durations.

The metabolic rate constants and specific activity of human and rat hepatic cytochrome P-450 2E1 toward toluene and chloroform

Chloroform (CHCl3) is a near-ubiquitous environmental contaminant, a by-product of the disinfection of drinking water sources and a commercially important compound. Standards for safe exposure have been established based on information defining its toxicity, which is mediated by metabolites. The metabolism of CHCl3 is via cytochrome P-450 2E1 (CYP2E1)-mediated oxidation to phosgene, which is known to obey a saturable mechanism. CYP2E1 is a highly conserved form, expressed in all mammalian systems studied, and is responsible for the metabolism of a great many low-molecular-weight (halogenated) compounds. However, the Michaelis-Menten rate constants for CHCl3 oxidation have not been derived in vitro, and the specific activity of CYP2E1 toward CHCl3 has not been reported. In this investigation with microsomal protein (MSP), apparent Vmax values of 27.6 and 28.3 nmol/h/mg MSP and apparent K(m) values of 1 and 0.15 microM in rats and human organ donors, respectively, were demonstrated. The specific activity of CYP2E1 toward CHCl3 in rats and humans was 5.29 and 5.24 pmol/min/pmol CYP2E1, respectively. Toluene metabolism to benzyl alcohol (BA), another CYP2E1-dependent reaction, was also highly dependent on CYP2E1 content in humans, and was more efficient than was CHCl3 metabolism. The specific activity of human CYP2E1 toward toluene metabolism in human MSP was 23 pmol/min/pmol CYP2E1. These results demonstrate that differences in CYP2E1 content of MSP among individuals and between species are largely responsible for observed differences in toluene and CHCl3 metabolism in vitro.

Physiologically based pharmacokinetic (PBPK) modeling of caffeine and theophylline in neonates and adults: implications for assessing children's risks from environmental agents

Children's risks can differ from those in adults for numerous reasons, one being differences in the pharmacokinetic handling of chemicals. Immature metabolism and a variety of other factors in neonates can affect chemical disposition and clearance. These factors can be incorporated into physiologically based pharmacokinetic (PBPK) models that simulate the fate of environmental toxicants in both children and adults. PBPK models are most informative when supported by empirical data, but typically pediatric pharmacokinetic data for toxicants are not available. In contrast, pharmacokinetic data in children are readily available for therapeutic drugs. The current analysis utilizes data for caffeine and theophylline, closely related xanthines that are both cytochrome P-450 (CYP) 1A2 substrates, in developing PBPK models for neonates and adults. Model development involved scale-up of in vitro metabolic parameters to whole liver and adjusting metabolic function for the ontological pattern of CYP1A2 and other CYPs. Model runs were able to simulate the large differences in half-life and clearance between neonates and adults. Further, the models were able to reproduce the faster metabolic clearance of theophylline relative to caffeine in neonates. This differential between xanthines was found to be due primarily to an extra metabolic pathway available to theophylline, back-methylation to caffeine, that is not available to caffeine itself. This pathway is not observed in adults exemplifying the importance of secondary or novel routes of metabolism in the immature liver. Greater CYP2E1 metabolism of theophylline relative to caffeine in neonates also occurs. Neonatal PBPK models developed for these drugs may be adapted to other CYP1A2 substrates (e.g., arylamine toxicants). A stepwise approach for modeling environmental toxicants in children is proposed.

Whole body pharmacokinetic models

The aim of the current review is to summarise the present status of physiologically based pharmacokinetic (PBPK) modelling and its applications in drug research, and thus serve as a reference point to people interested in the methodology. The review is structured into three major sections. The first discusses the existing methodologies and techniques of PBPK model development. The second describes some of the most interesting PBPK model implementations published. The final section is devoted to a discussion of the current limitations and the possible future developments of the PBPK modelling approach. The current review is focused on papers dealing with the pharmacokinetics and/or toxicokinetics of medicinal compounds; references discussing PBPK models of environmental compounds are mentioned only if they represent considerable methodological developments or reveal interesting interpretations and/or applications.The major conclusion of the review is that, despite its significant potential, PBPK modelling has not seen the development and implementation it deserves, especially in the drug discovery, research and development processes. The main reason for this is that the successful development and implementation of a PBPK model is seen to require the investment of significant experience, effort, time and resources. Yet, a substantial body of PBPK-related research has been accumulated that can facilitate the PBPK modelling and implementation process. What is probably lagging behind is the expertise component, where the demand for appropriately qualified staff far outreaches availability.

The impact of cytochrome P450 2E1-dependent metabolic variance on a risk-relevant pharmacokinetic outcome in humans

Risk assessments include assumptions about sensitive subpopulations, such as the fraction of the general population that is sensitive and the extent that biochemical or physiological attributes influence sensitivity. Uncertainty factors (UF) account for both pharmacokinetic (PK) and pharmacodynamic (PD) components, allowing the inclusion of risk-relevant information to replace default assumptions about PK and PD variance (uncertainty). Large numbers of human organ donor samples and recent advances in methods to extrapolate in vitro enzyme expression and activity data to the intact human enable the investigation of the impact of PK variability on human susceptibility. The hepatotoxicity of trichloroethylene (TCE) is mediated by acid metabolites formed by cytochrome P450 2E1 (CYP2E1) oxidation, and differences in the CYP2E1 expression are hypothesized to affect susceptibility to TCE's liver injury. This study was designed specifically to examine the contribution of statistically quantified variance in enzyme content and activity on the risk of hepatotoxic injury among adult humans. We combined data sets describing (1) the microsomal protein content of human liver, (2) the CYP2E1 content of human liver microsomal protein, and (3) the in vitro Vmax for TCE oxidation by humans. The 5th and 95th percentiles of the resulting distribution (TCE oxidized per minute per gram liver) differed by approximately sixfold. These values were converted to mg TCE oxidized/h/kg body mass and incorporated in a human PBPK model. Simulations of 8-hour inhalation exposure to 50 ppm and oral exposure to 5 micro g TCE/L in 2 L drinking water showed that the amount of TCE oxidized in the liver differs by 2% or less under extreme values of CYP2E1 expression and activity (here, selected as the 5th and 95th percentiles of the resulting distribution). This indicates that differences in enzyme expression and TCE oxidation among the central 90% of the adult human population account for approximately 2% of the difference in production of the risk-relevant PK outcome for TCE-mediated liver injury. Integration of in vitro metabolism information into physiological models may reduce the uncertainties associated with risk contributions of differences in enzyme expression and the UF that represent PK variability.

Inter-individual variability in levels of human microsomal protein and hepatocellularity per gram of liver

AIMS

To determine levels of microsomal protein (MPPGL) and hepatocellularity (HPGL) per gram of human liver and their interindividual variability.

METHODS

Triplicate liver samples were used to determine values of MPPGL (n = 20) and HPGL (n = 7) after accounting for the fractional loss of microsomal protein or hepatocytes during processing. Repeated measurements from each liver sample allowed the estimation of true interindividual variability in MPPGL and HPGL using ANOVA.

RESULTS

The value of MPPGL ranged from 26 to 54 mg g(-1) (mean(geo)= 33 mg g(-1)). The value of HPGL ranged from 65 to 185 x 10(6) cells g(-1) (mean(geo)= 10(7) x 10(6) cells g(-1)).

CONCLUSIONS

There is significant interindividual variability in MPPGL, which has implications for the accurate extrapolation of in vitro data on drug metabolism to predict in vivo metabolic clearance.

Identification of trichloroethylene and its metabolites in human seminal fluid of workers exposed to trichloroethylene

We have investigated the potential of the male reproductive tract to accumulate trichloroethylene (TCE) and its metabolites, including chloral, trichloroethanol (TCOH), trichloroacetic acid (TCA), and dichloroacetic acid (DCA). Human seminal fluid and urine samples from eight mechanics diagnosed with clinical infertility and exposed to TCE occupationally were analyzed. In in vivo experimental studies, TCE and its metabolites were determined in epididymis and testis of mice exposed to TCE (1000 ppm) by inhalation for 1 to 4 weeks. In other studies, incubations of monkey epididymal microsomes were performed in the presence of TCE and NADPH. Our results showed that seminal fluid from all eight subjects contained TCE, chloral, and TCOH. DCA was present in samples from two subjects, and only one contained TCA. TCA and/or TCOH were also identified in urine samples from only two subjects. TCE, chloral, and TCOH were detected in murine epididymis after inhalation exposure with TCE for 1 to 4 weeks. Levels of TCE and chloral were similar throughout the entire exposure period. TCOH levels were similar at 1 and 2 weeks but increased significantly after 4 weeks of TCE exposure. Chloral was identified in microsomal incubations with TCE in monkey epididymis. CYP2E1, a P450 that metabolizes TCE, was localized in human and monkey epididymal epithelium and testicular Leydig cells. These results indicated that TCE is metabolized in the reproductive tract of the mouse and monkey. Furthermore, TCE and its metabolites accumulated in seminal fluid, and suggested associations between production of TCE metabolites, reproductive toxicity, and impaired fertility.

Toxicokinetics and physiologically based toxicokinetics in toxicology and risk assessment

Toxicokinetics is the study of kinetics of absorption, distribution, metabolism, and excretion of a xenobiotic under the conditions of toxicity evaluation. Conventional toxicokinetics uses the hypothetical compartments, and the model is composed of rate equations that describe the time course of drug and chemical disposition. The utility of toxicokinetics in toxicity evaluation and interpretation of animal toxicology data is emerging as an important tool in product discovery and development. With implementation of the International Conference on Harmonization (ICH) guidelines on systemic exposure and dose selection, toxicokinetics have been integrated in routine toxicity evaluations. Although traditional compartmental/noncompartmental models are generally adequate for assessing systemic exposure, they are unable to the predict time course of drug disposition in target tissues and often fail to relate systemic drug levels to a biological response. Physiologically based toxicokinetic (PB-TK) models address this deficiency of traditional compartmental models. PB-TK models are the kinetic models of the uptake and disposition of chemicals based on rates of biochemical reactions, physiological and anatomical characteristics. These models, when developed appropriately, can predict target organ drug distribution in different species under variety of conditions. This minireview discusses the basic principles, and applications of traditional compartmental toxicokinetic and physiologically based toxicokinetics (PB-TK) models in drug development and risk assessment. Special emphasis will be placed on discussion related to interpretation of the ICH guidelines related to toxicokinetics and the utility of toxicokinetics data in dose selection for toxicity and carcinogenicity studies. The utility of PB-TK models in risk assessment of methylene chloride, vinyl chloride, retinoic acid, dioxin, and inhaled organic esters is discussed.

Incorporating human interindividual biotransformation variance in health risk assessment

The protection of sensitive individuals within a population dictates that measures other than central tendencies be employed to estimate risk. The refinement of human health risk assessments for chemicals metabolized by the liver to reflect data on human variability can be accomplished through (1) the characterization of enzyme expression in large banks of human liver samples, (2) the employment of appropriate techniques for the quantification and extrapolation of metabolic rates derived in vitro, and (3) the judicious application of physiologically based pharmacokinetic (PBPK) modeling. While in vitro measurements of specific biochemical reactions from multiple human samples can yield qualitatively valuable data on human variance, such measures must be put into the perspective of the intact human to yield the most valuable predictions of metabolic differences among humans. For quantitative metabolism data to be the most valuable in risk assessment, they must be tied to human anatomy and physiology, and the impact of their variance evaluated under real exposure scenarios. For chemicals metabolized in the liver, the concentration of parent chemical in the liver represents the substrate concentration in the Michaelis Menten description of metabolism. Metabolic constants derived in vitro may be extrapolated to the intact liver, when appropriate conditions are met. Metabolic capacity Vmax; the maximal rate of the reaction) can be scaled directly to the concentration of enzyme (or enzyme fraction) contained in the liver. Several environmental, genetic and lifestyle factors can influence the concentration of cytochrome P450 forms (CYP) in the liver by affecting either (1) the extent to which the CYP forms are expressed in the endoplasmic reticulum of the cell (isolated as the microsomal fraction from tissue homogenates), or (2) the expression of microsomal protein in intact liver tissue. Biochemically sound measures of the hepatic distribution of xenobiotic metabolizing enzymes among humans, based on expression of the enzymes within microsomal protein and the distribution of microsomal protein among intact livers, can be combined with metabolic constants derived in vitro to generate values consistent with those employed in PBPK models. When completed, the distribution (and bounds) of Vmax values can be estimated and included in PBPK models. Exercising such models under plausible exposure scenarios will demonstrate the extent to which human interindividual enzyme variance can influence parameters (i.e., the detoxication of a toxic chemical through metabolism) that may influence risk. In this article, we describe a methodology and conditions which must exist for such an approach to be successful.