In vitro measurements of metabolism for application in pharmacokinetic modeling

The potential of enteroids derived from children and adults to study age-dependent differences in intestinal CYP3A4/5 metabolism

Drug metabolism in the intestinal wall affects bioavailability of orally administered drugs and is influenced by age. Hence, it is important to fully understand the drug metabolizing capacity of the gut to predict systemic exposure. The aim of this study was to investigate the potential of enteroids as a tool to study CYP3A4/5 -mediated metabolism in both children and adults. Bioconversion of midazolam, a CYP3A4/5 model substrate, was studied using enteroid monolayers as well as tissue explants in the Ussing chamber, both derived from pediatric [median (range age): 54 weeks (2 days - 13 years), n = 21] and adult (n = 5) tissue. Caco-2 cellular monolayers were employed as controls. In addition, mRNA expression of CYP3A4 was determined in enteroid monolayers (n = 11), tissue (n = 23) and Caco-2 using RT-qPCR. Midazolam metabolism was successfully detected in all enteroid monolayers, as well as in all tissue explants studied in the Ussing chamber, whereas Caco-2 showed no significant metabolite formation. The extracted fraction of midazolam was similar between enteroid monolayers and tissue. The fraction of midazolam extracted increased with age in enteroid monolayers derived from 0 to 70 week old donors. No statistically significant correlation was observed in tissue likely due to high variability observed and the smaller donor numbers included in the study. At the level of gene expression, CYP3A4 increased with age in tissues (n = 32), while this was not reflected in enteroid monolayers (n = 16). Notably, asymmetric metabolite formation was observed in enteroids and tissue, with higher metabolite formation on the luminal side of the barrier. In summary, we demonstrated that enteroids can be used to measure CYP3A4/5 midazolam metabolism, which we show is similar as observed in fresh isolated tissue. This was the case both in children and adults, indicating the potential of enteroids to predict intestinal metabolism. This study provides promising data to further develop enteroids to study drug metabolism in vitro and potentially predict oral absorption for special populations as an alternative to using fresh tissue.

Exploring Kp,uu,BBB values smaller than unity in remoxipride: A physiologically-based CNS model approach highlighting brain metabolism in drugs with passive blood-brain barrier transport

(AIM)

Kp,uu,BBB values are crucial indicators of drug distribution into the brain, representing the steady-state relationship between unbound concentrations in plasma and in brain extracellular fluid (brainECF). Kp,uu,BBB values < 1 are often interpreted as indicators of dominant active efflux transport processes at the blood-brain barrier (BBB). However, the potential impact of brain metabolism on this value is typically not addressed. In this study, we investigated the brain distribution of remoxipride, as a paradigm compound for passive BBB transport with yet unexplained brain elimination that was hypothesized to represent brain metabolism.

(METHODS)

The physiologically-based LeiCNS pharmacokinetic predictor (LeiCNS-PK model) was used to compare brain distribution of remoxipride with and without Michaelis-Menten kinetics at the BBB and/or brain cell organelle levels. To that end, multiple in-house (IV 0.7, 3.5, 4, 5.2, 7, 8, 14 and 16 mg kg-1) and external (IV 4 and 8 mg kg-1) rat microdialysis studies plasma and brainECF data were analysed.

(RESULTS)

The incorporation of active elimination through presumed brain metabolism of remoxipride in the LeiCNS-PK model significantly improved the prediction accuracy of experimentally observed brainECF profiles of this drug. The model integrated with brain metabolism in both barriers and organelles levels is named LeiCNS-PK3.5.

(CONCLUSION)

For drugs with Kp,uu,BBB values < 1, not only the current interpretation of dominant BBB efflux transport, but also potential brain metabolism needs to be considered, especially because these may be concentration dependent. This will improve the mechanistic understanding of the processes that determine brain PK profiles.

Interspecies comparison of metabolism of two novel prototype PFAS

As a result of proposed global restrictions and regulations on current-use per-and polyfluoroalkyl substances (PFAS), research on possible alternatives is highly required. In this study, phase I in vitro metabolism of two novel prototype PFAS in human and rat was investigated. These prototype chemicals are intended to be safer-by-design and expected to mineralize completely, and thus be less persistent in the environment compared to the PFAS available on the market. Following incubation with rat liver S9 (RL-S9) fractions, two main metabolites per initial substance were produced, namely an alcohol and a short-chain carboxylic acid. While with human liver S9 (HL-S9) fractions, only the short-chain carboxylic acid was detected. Beyond these major metabolites, two and five additional metabolites were identified at very low levels by non-targeted screening for the ether- and thioether-linked prototype chemicals, respectively. Overall, complete mineralization during the in vitro hepatic metabolism of these novel PFAS by HL-S9 and RL-S9 fractions was not observed. The reaction kinetics of the surfactants was determined by using the metabolite formation, rather than the substrate depletion approach. With rat liver enzymes, the formation rates of primary metabolite alcohols were at least two orders of magnitude higher than those of secondary metabolite carboxylic acids. When incubating with human liver enzymes, the formation rates of single metabolite carboxylic acids, were similar or smaller than those experienced in rat. It also indicates that the overall metabolic rate and clearance of surfactants are significantly higher in rat liver than in human liver. The maximum formation rate of the thioether congener exceeded 10-fold that of the ether in humans but were similar in rats. Overall, the results suggest that metabolism of the prototype chemicals followed a similar trend to those reported in studies of fluorotelomer alcohols.

A physiologically-based pharmacokinetic/pharmacodynamic (PBPK/PD) model for the insecticide dimethoate

1. Dimethoate is an organophosphate insecticide that is converted in vivo to omethoate, the active toxic moiety. Omethoate inhibits acetylcholinesterase (AChE) in the brain and red blood cells (RBCs). This paper describes the development of rat and human physiologically-based pharmacokinetic/pharmacodynamic (PBPK/PD) models for dimethoate.2. The model simulates the absorption and distribution of dimethoate and omethoate, the conversion of dimethoate to omethoate and to other metabolites, the metabolism and excretion of omethoate, and the inhibition of RBC and brain AChE. An extensive data collection program to estimate metabolism and inhibition parameters is described.3. The suite of models includes an adult rat, post-natal rat, and human model. The rat models were evaluated by comparing model predictions of dimethoate and omethoate to measured blood time course data, and with RBC and brain AChE inhibition estimates from an extensive database of in vivo AChE measurements.4. After the demonstration of adequately fitted rat models that were robust to sensitivity analysis, the human model was applied for estimation of points-of-departure (PODs) for risk assessment using the human-specific parameters in the human PBPK/PD model. Thus, the standard interspecies uncertainty factor can be reduced from 10X to 1X.

The Role of "Physiologically Based Pharmacokinetic Model (PBPK)" New Approach Methodology (NAM) in Pharmaceuticals and Environmental Chemical Risk Assessment

Physiologically Based Pharmacokinetic (PBPK) models are mechanistic tools generally employed in the pharmaceutical industry and environmental health risk assessment. These models are recognized by regulatory authorities for predicting organ concentration-time profiles, pharmacokinetics and daily intake dose of xenobiotics. The extension of PBPK models to capture sensitive populations such as pediatric, geriatric, pregnant females, fetus, etc., and diseased populations such as those with renal impairment, liver cirrhosis, etc., is a must. However, the current modelling practices and existing models are not mature enough to confidently predict the risk in these populations. A multidisciplinary collaboration between clinicians, experimental and modeler scientist is vital to improve the physiology and calculation of biochemical parameters for integrating knowledge and refining existing PBPK models. Specific PBPK covering compartments such as cerebrospinal fluid and the hippocampus are required to gain mechanistic understanding about xenobiotic disposition in these sub-parts. The PBPK model assists in building quantitative adverse outcome pathways (qAOPs) for several endpoints such as developmental neurotoxicity (DNT), hepatotoxicity and cardiotoxicity. Machine learning algorithms can predict physicochemical parameters required to develop in silico models where experimental data are unavailable. Integrating machine learning with PBPK carries the potential to revolutionize the field of drug discovery and development and environmental risk. Overall, this review tried to summarize the recent developments in the in-silico models, building of qAOPs and use of machine learning for improving existing models, along with a regulatory perspective. This review can act as a guide for toxicologists who wish to build their careers in kinetic modeling.

IVIVE: Facilitating the Use of In Vitro Toxicity Data in Risk Assessment and Decision Making

During the past few decades, the science of toxicology has been undergoing a transformation from observational to predictive science. New approach methodologies (NAMs), including in vitro assays, in silico models, read-across, and in vitro to in vivo extrapolation (IVIVE), are being developed to reduce, refine, or replace whole animal testing, encouraging the judicious use of time and resources. Some of these methods have advanced past the exploratory research stage and are beginning to gain acceptance for the risk assessment of chemicals. A review of the recent literature reveals a burst of IVIVE publications over the past decade. In this review, we propose operational definitions for IVIVE, present literature examples for several common toxicity endpoints, and highlight their implications in decision-making processes across various federal agencies, as well as international organizations, including those in the European Union (EU). The current challenges and future needs are also summarized for IVIVE. In addition to refining and reducing the number of animals in traditional toxicity testing protocols and being used for prioritizing chemical testing, the goal to use IVIVE to facilitate the replacement of animal models can be achieved through their continued evolution and development, including a strategic plan to qualify IVIVE methods for regulatory acceptance.

Considerations for Improving Metabolism Predictions for In Vitro to In Vivo Extrapolation

High-throughput (HT) in vitro to in vivo extrapolation (IVIVE) is an integral component in new approach method (NAM)-based risk assessment paradigms, for rapidly translating in vitro toxicity assay results into the context of in vivo exposure. When coupled with rapid exposure predictions, HT-IVIVE supports the use of HT in vitro assays for risk-based chemical prioritization. However, the reliability of prioritization based on HT bioactivity data and HT-IVIVE can be limited as the domain of applicability of current HT-IVIVE is generally restricted to intrinsic clearance measured primarily in pharmaceutical compounds. Further, current approaches only consider parent chemical toxicity. These limitations occur because current state-of-the-art HT prediction tools for clearance and metabolite kinetics do not provide reliable data to support HT-IVIVE. This paper discusses current challenges in implementation of IVIVE for prioritization and risk assessment and recommends a path forward for addressing the most pressing needs and expanding the utility of IVIVE.

Characterization of furathiocarb metabolism in in vitro human liver microsomes and recombinant cytochrome P450 enzymes

Furathiocarb is a carbamate insecticide detected in ecosystems. Its main metabolite carbofuran has been alluded to affect birth outcomes and disturb hormone levels in humans. The metabolism of furathiocarb in humans has not been characterized. The metabolism studies were performed using hepatic microsomes from ten donors and fifteen human cDNA-expressed CYPs. The initial screening and identification of the metabolites were performed by LC-TOF. Quantifications and fragmentations were performed by LC/MS-MS. Furathiocarb was metabolized to eight phase I metabolites via two general pathways, carbofuran metabolic pathway and furathiocarb oxidation pathway. Six metabolites in the carbofuran metabolic pathway (carbofuran, 3-hydroxycarbofuran, 3-ketocarbofuran, 3-keto-7-phenolcarbofuran, 3-hydroxy-7-phenolcarbofuran, and 7-phenolcarbofuran) were identified with the help of authentic standards. The two unidentified metabolites in the furathiocarb oxidation pathway are probably hydroxylated and sulfoxidated derivatives of furathiocarb. The carbofuran metabolic pathway was more predominant than the furathiocarb oxidation pathway, ratios ranged from 24- to 115-fold in a 10-donor panel of hepatic microsomes. On the basis of recombinant CYP studies, the carbofuran pathway was dominated by CYP3A4 (95.9%); contributions by CYP1A2 (1.3%) and CYP2B6 (2.0%) were minor. The minor furathiocarb oxidation pathway was catalyzed by CYP2C19 and CYP2D6 (hydroxylated/sulfoxidated metabolite A) and by CYP3A5, CYP3A4 and CYP2A6 (metabolite B). High and significant correlation between carbofuran metabolic pathway and CYP3A4 marker activities (midazolam-1'-hydroxylation and omeprazole-sulfoxidation) were observed. Ketoconazole, a CYP3A4-inhibitor, inhibited the carbofuran pathway by 32–86% and hydroxylated/sulfoxidated metabolite-B formations by 41–62%. The data suggest that in humans, the carbofuran metabolic pathway is dominant, and CYP3A4 is the major enzyme involved. These results provide useful scientific information for furathiocarb risk assessment in humans.

•

Eight Phase I metabolites were detected by LC-TOF-MS/MS.

•

The carbofuran pathway was more rapid than the furathiocarb oxidation pathway

•

The carbofuran pathway was dominated by CYP3A4 (96%).

•

Ketoconazole inhibited the carbofuran pathway by 32–86%.

•

The findings provide useful scientific information for furathiocarb risk assessment in humans.

Quaternary Ammonium Compounds: Bioaccumulation Potentials in Humans and Levels in Blood before and during the Covid-19 Pandemic

Quaternary ammonium compounds (QACs) are commonly used in a variety of consumer, pharmaceutical, and medical products. In this study, bioaccumulation potentials of 18 QACs with alkyl chain lengths of C8-C18 were determined in the in vitro-in vivo extrapolation (IVIVE) model using the results of human hepatic metabolism and serum protein binding experiments. The slowest in vivo clearance rates were estimated for C12-QACs, suggesting that these compounds may preferentially build up in blood. The bioaccumulation of QACs was further confirmed by the analysis of human blood (sera) samples (n = 222). Fifteen out of the 18 targeted QACs were detected in blood with the ΣQAC concentrations reaching up to 68.6 ng/mL. The blood samples were collected during two distinct time periods: before the outbreak of the COVID-19 pandemic (2019; n = 111) and during the pandemic (2020, n = 111). The ΣQAC concentrations were significantly higher in samples collected during the pandemic (median 6.04 ng/mL) than in those collected before (median 3.41 ng/mL). This is the first comprehensive study on the bioaccumulation and biomonitoring of the three major QAC groups and our results provide valuable information for future epidemiological, toxicological, and risk assessment studies targeting these chemicals.

Human internal exposure to organophosphate esters: A short review of urinary monitoring on the basis of biological metabolism research

As alternatives to traditional brominated flame retardants, organophosphate flame retardants (OPFRs), especially for organophosphate esters (OPEs) -- the most widely used and investigated OPFRs, have raised people's concern on their environmental and health-related risks over the years. Considering their extensive environmental occurrence and potential adverse effects, precise estimation on the human body burden of OPEs will be conducive to the restrictions on the usage of these compounds scientifically. Biomonitoring research can provide precise information on human exposure to OPEs as it reveals the degree of external exposure from all exposure routes. Knowledge on biotransformation and metabolism of OPEs in the biosystems is of great significance for our understanding of the internal exposure to these compounds. In this study, the biological metabolic processes of nine OPEs prevalent in the environment, involving tris(2-chloroethyl) phosphate (TCEP), tris(1-chloro-2-propyl) phosphate (TCIPP), tris(1,3-dichloro-2-propyl) phosphate (TDCIPP), tripropyl phosphate (TPrP), tri-n-butyl phosphate (TnBP), tris(2-butoxyethyl) phosphate (TBOEP), triphenyl phosphate (TPhP), 2-ethylhexyl diphenyl phosphate (EHDPP), and tricresyl phosphate (TCrP), are comprehensively reviewed. Specifically, the metabolic pathway, kinetics and mechanism of OPEs are depicted in detail. Under this context, the advances and limitations on biomonitoring of OPE metabolites in human urine are summarized. The requirements of specificity, quantitative stability, high detection frequency/concentration are needed for OPE metabolites to be considered and validated as biomarkers. Thus far, deeper elucidations on the metabolic processes and identification of biomarkers of OPEs are urgently required, given that some OPEs have no suitable biomarkers in human biomonitoring. For better assessment of the body burden of OPEs in humans, reliable and effective methodologies for urine sampling and estimation on internal exposure to OPEs need to be further developed in the future.

Arsenic toxicokinetic modeling and risk analysis: Progress, needs and applications

Arsenic (As) poses unique challenges in PBTK model development and risk analysis applications. Arsenic metabolism is complex, adequate information to attribute specific metabolites to particular adverse effects in humans is sparse, and measurement of relevant metabolites in biological media can be difficult. Multiple As PBTK models have been published and used or adapted for use in various exposure and risk analysis applications. These applications illustrate the broad utility of PBTK models for exposure and dose-response analysis, particularly for arsenic where multi-pathway, multi-route exposures and multiple toxic effects are of concern. Arsenic PBTK models have been used together with exposure reconstruction and dose-response functions to estimate risk of specific adverse health effects due to drinking water exposure and consumption of specific foodstuffs (e.g. rice, seafood), as well as to derive safe exposure levels and develop consumption advisories. Future refinements to arsenic PBTK models can enhance the confidence in such analyses. Improved estimates for methylation biotransformation parameters based on in vitro to in vivo extrapolation (IVIVE) methods and estimation of interindividual variability in key model parameters for specific toxicologically relevant metabolites are two important areas for consideration.

Estimation of toluene exposure in air from BMA (S-benzylmercapturic acid) urinary measures using a reverse dosimetry approach based on physiologically pharmacokinetic modeling

This study aimed to use a reverse dosimetry PBPK modeling approach to estimate toluene atmospheric exposure from urinary measurements of S-benzylmercapturic acid (BMA) in a small group of individuals and to evaluate the uncertainty associated to urinary spot-sampling compared to 24-h collected urine samples. Each exposure assessment technique was developed namely to estimate toluene air exposure from BMA measurements in 24-h urine samples (24-h-BMA) and from distributions of daily urinary BMA spot measurements (DUBSM). Model physiological parameters were described based upon age, weight, size and sex. Monte Carlo simulations with the PBPK model allowed converting DUBSM distribution (and 24-h-BMA) into toluene air levels. For the approach relying on DUBSM distribution, the ratio between the 95% probability of predicted toluene concentration and its 50% probability in each individual varied between 1.2 and 1.4, while that based on 24-h-BMA varied between 1.0 and 1.1. This suggests more variability in estimated exposure from spot measurements. Thus, estimating toluene exposure based on DUBSM distribution generated about 20% more uncertainty. Toluene levels estimated (0.0078-0.0138 ppm) are well below Health Canada's maximum chronic air guidelines. PBPK modeling and reverse dosimetry may be combined to interpret urinary metabolites data of VOCs and assess related uncertainties.

Comparison of in vivo derived and scaled in vitro metabolic rate constants for several volatile organic compounds (VOCs)

Metabolic rate parameters estimation using in vitro data is necessary due to numbers of chemicals for which data are needed, trend towards minimizing laboratory animal use, and limited opportunity to collect data in human subjects. We evaluated how well metabolic rate parameters derived from in vitro data predict overall in vivo metabolism for a set of environmental chemicals for which well validated and established methods exist. We compared values of VmaxC derived from in vivo vapor uptake studies with estimates of VmaxC scaled up from in vitro hepatic microsomal metabolism studies for VOCs for which data were available in male F344 rats. For 6 of 7 VOCs, differences between the in vivo and scaled up in vitro VmaxC estimates were less than 2.6-fold. For bromodichloromethane (BDCM), the in vivo derived VmaxC was approximately 4.4-fold higher than the in vitro derived and scaled up VmaxC. The more rapid rate of BDCM metabolism estimated based in vivo studies suggests other factors such as extrahepatic metabolism, binding or other non-specific losses making a significant contribution to overall clearance. Systematic and reliable utilization of scaled up in vitro biotransformation rate parameters in PBPK models will require development of methods to predict cases in which extrahepatic metabolism and binding as well as other factors are likely to be significant contributors.

Estimating renal and hepatic clearance rates of organophosphate esters in humans: Impacts of intrinsic metabolism and binding affinity with plasma proteins

The renal and hepatic clearance rates of organophosphate esters (OPEs) in humans were estimated. Six OPEs and their corresponding diester metabolites (mOPEs) were quantified respectively in 30 paired human plasma and urine samples collected in Hengshui, Hebei province, China. The renal clearance rate (CL_renal) of triphenyl phosphate (TPHP), tris(chloroethyl) phosphate (TCEP) and tris(1,3-dichloro-isopropyl) phosphate (TDCIPP) was estimated to be 68.9, 50.9 and 33.3 mL/kg/day, respectively, while it was not calculated for other three OPEs due to the low detection frequencies in human samples. To estimate the clearance rates of the target OPEs, hepatic clearance rates (CL_h) of OPEs were extrapolated from their in vitro intrinsic clearance data in human liver microsomes (CL_HLM). The calculated CL_h values of TCEP and TDCIPP were comparable to their CL_renal, indicating that the in vitro extrapolation method was suitable for estimating the clearance rates of OPEs. The higher binding affinity of TDCIPP with plasma proteins could reduce its renal clearance. The estimated half-lives of Cl-OPEs in human were longer than those of the aryl- and alkyl-OPEs. This study provided a feasible in vitro method to predict the clearance and half-lives of OPEs in human, which is significant for their accurate health risk assessment.

First insight into human extrahepatic metabolism of flame retardants: Biotransformation of EH-TBB and Firemaster-550 components by human skin subcellular fractions

2-ethylhexyl-2,3,4,5-tetrabromobenzoate (EH-TBB) and a mixture of EH-TBB, Bis(2-ethylhexyl)tetrabromphthalate (BEH-TEBP) and Triphenyl phosphate (TPhP), prepared in a ratio similar to the Firemaster-550™ (FM550) flame retardant formulation, were exposed to human skin subcellular fractions (S9) to evaluate their dermal in vitro metabolism for the first time. After 60 min of incubation, tetrabromobenzoic acid (TBBA) and diphenyl phosphate (DPhP) were identified as the major metabolites of EH-TBB and TPhP, respectively using UPLC-Q-Exactive Orbitrap™-MS analysis. Dermal biotransformation of EH-TBB and TPhP was catalyzed by skin carboxylesterases rather than CYP450 enzymes, while no stable metabolites could be identified for BEH-TEBP. Metabolite formation rates of EH-TBB as individual compound and as a component of FM550 fitted the Michaelis-Menten model, while no steady state could be reached for TPhP under experimental conditions. Estimated maximum metabolic rate (V_max) for TBBA formation upon exposure to FM550 was lower than V_max for EH-TBB (1.08 and 15.2 pmol min^-1 mg protein^-1, respectively). This indicates dermal metabolism would contribute less to the clearance of EH-TBB body burden than hepatic metabolism (V_max = 644 pmol min^-1 mg protein^-1). Implications for human exposure include EH-TBB accumulation in skin tissue and human exposure to dermal metabolic products, which may have different toxicokinetic and toxicodynamic parameters than parent flame retardants.

Relevance of In Vitro Metabolism Models to PET Radiotracer Development: Prediction of In Vivo Clearance in Rats from Microsomal Stability Data

The prediction of in vivo clearance from in vitro metabolism models such as liver microsomes is an established procedure in drug discovery. The potentials and limitations of this approach have been extensively evaluated in the pharmaceutical sector; however, this is not the case for the field of positron emission tomography (PET) radiotracer development. The application of PET radiotracers and classical drugs differs greatly with regard to the amount of substance administered. In typical PET imaging sessions, subnanomolar quantities of the radiotracer are injected, resulting in body concentrations that cannot be readily simulated in analytical assays. This raises concerns regarding the predictability of radiotracer clearance from in vitro data. We assessed the accuracy of clearance prediction for three prototypical PET radiotracers developed for imaging the A₁ adenosine receptor (A₁AR). Using the half-life (t_1/2) approach and physiologically based scaling, in vivo clearance in the rat model was predicted from microsomal stability data. Actual clearance could be accurately predicted with an average fold error (AFE) of 0.78 and a root mean square error (RMSE) of 1.6. The observed slight underprediction (1.3-fold) is in accordance with the prediction accuracy reported for classical drugs. This result indicates that the prediction of radiotracer clearance is possible despite concentration differences of more than three orders of magnitude between in vitro and in vivo conditions. Consequently, in vitro metabolism models represent a valuable tool for PET radiotracer development.

Relevance of CYP2C9 Function in Valproate Therapy

BACKGROUND

Genetic polymorphisms of drug metabolizing enzymes can substantially modify the pharmacokinetics of a drug and eventually its efficacy or toxicity; however, inferring a patient's drug metabolizing capacity merely from his or her genotype can lead to false prediction. Non-genetic host factors (age, sex, disease states) and environmental factors (nutrition, comedication) can transiently alter the enzyme expression and activities resulting in genotypephenotype mismatch. Although valproic acid is a well-tolerated anticonvulsant, pediatric patients are particularly vulnerable to valproate injury that can be partly attributed to the age-related differences in metabolic pathways.

METHODS

CYP2C9 mediated oxidation of valproate, which is the minor metabolic pathway in adults, appears to become the principal route in children. Genetic and non-genetic variations in CYP2C9 activity can result in significant inter- and intra-individual differences in valproate pharmacokinetics and valproate induced adverse reactions.

RESULTS

The loss-of-function alleles, CYP2C9*2 or CYP2C9*3, display significant reduction in valproate metabolism in children; furthermore, low CYP2C9 expression in patients with CYP2C9*1/*1 genotype also leads to a decrease in valproate metabolizing capacity. Due to phenoconversion, the homozygous wild genotype, expected to be translated to CYP2C9 enzyme with normal activity, is transiently switched into poor (or extensive) metabolizer phenotype.

CONCLUSION

Novel strategy for valproate therapy adjusted to CYP2C9-status (CYP2C9 genotype and CYP2C9 expression) is strongly recommended in childhood. The early knowledge of pediatric patients' CYP2C9-status facilitates the optimization of valproate dosing which contributes to the avoidance of misdosing induced adverse reactions, such as abnormal blood levels of ammonia and alkaline phosphatase, and improves the safety of children's anticonvulsant therapy.

The Impact of Scaling Factor Variability on Risk-Relevant Pharmacokinetic Outcomes in Children: A Case Study Using Bromodichloromethane (BDCM)

Biotransformation rates extrapolated from in vitro data are used increasingly in human physiologically based pharmacokinetic (PBPK) models. This practice requires use of scaling factors, including microsomal content (mg of microsomal protein/g liver, MPPGL), enzyme specific content, and liver mass as a fraction of body weight (FVL). Previous analyses indicated that scaling factor variability impacts pharmacokinetic (PK) outcomes used in adult population dose-response studies. This analysis was extended to pediatric populations because large inter-individual differences in enzyme ontogeny likely would further contribute to scaling factor variability. An adult bromodichloromethane (BDCM) model (Kenyon, E. M., Eklund, C., Leavens, T. L., and Pegram, R. A. (2016a). Development and application of a human PBPK model for bromodichloromethane (BDCM) to investigate impacts of multi-route exposure. J. Appl. Toxicol. 36, 1095-1111) was re-parameterized for neonates, infants, and toddlers. Monte Carlo analysis was used to assess the impact of pediatric scaling factor variation on model-derived PK outcomes compared with adult findings. BDCM dose metrics were estimated following a single 0.05-liter drink of water or a 20-min bath, under typical (5 µg/l) and plausible higher (20 µg/l) BDCM concentrations. MPPGL, CYP2E1, and FVL values reflected the distribution of reported pediatric population values. The impact of scaling factor variability on PK outcome variation was different for each exposure scenario, but similar for each BDCM water concentration. The higher CYP2E1 expression variability during early childhood was reflected in greater variability in predicted PK outcomes in younger age groups, particularly for the oral exposure route. Sensitivity analysis confirmed the most influential parameter for this variability was CYP2E1, particularly in neonates. These findings demonstrate the importance of age-dependent scaling factor variation used for in vitro to in vivo extrapolation of biotransformation rates.

Reverse dosimetry modeling of toluene exposure concentrations based on biomonitoring levels from the Canadian health measures survey

Biomonitoring might provide useful estimates of population exposure to environmental chemicals. However, data uncertainties stemming from interindividual variability are common in large population biomonitoring surveys. Physiologically based pharmacokinetic (PBPK) models might be used to account for age- and gender-related variability in internal dose. The objective of this study was to reconstruct air concentrations consistent with blood toluene measures reported in the third Canadian Health Measures Survey using reverse dosimetry PBPK modeling techniques. Population distributions of model's physiological parameters were described based upon age, weight, and size for four subpopulations (12-19, 20-39, 40-59, and 60-79 years old). Monte Carlo simulations applied to PBPK modeling allowed converting the distributions of venous blood measures of toluene obtained from CHMS into related air levels. Based upon blood levels observed at the 50^th, 90^th and 95^th percentiles, corresponding air toluene concentrations were estimated for teenagers aged 12-19 years as being, respectively, 0.009, 0.04 and 0.06 ppm. Similarly, values were computed for adults aged 20-39 years (0.007, 0.036, and 0.06 ppm), 40-59 years (0.007, 0.036 and 0.06 ppm) and 60-79 years (0.006, 0.022 and 0.04 ppm). These estimations are well below Health Canada's maximum recommended chronic air guidelines for toluene. In conclusion, PBPK modeling and reverse dosimetry may be combined to help interpret biomonitoring data for chemical exposure in large population surveys and estimate the associated toxicological health risk.

Microbial biotransformation of bioactive and clinically useful steroids and some salient features of steroids and biotransformation

Steroids are perhaps one of the most widely used group of drugs in present day. Beside the established utilization as immunosuppressive, anti-inflammatory, anti-rheumatic, progestational, diuretic, sedative, anabolic and contraceptive agents, recent applications of steroid compounds include the treatment of some forms of cancer, osteoporosis, HIV infections and treatment of declared AIDS. Steroids isolated are often available in minute amounts. So biotransformation of natural products provides a powerful means in solving supply problems in clinical trials and marketing of the drug for obtaining natural products in bulk amounts. If the structure is complex, it is often an impossible task to isolate enough of the natural products for clinical trials. The microbial biotransformation of steroids yielded several novel metabolites, exhibiting different activities. The metabolites produced from pregnenolone acetate by Cunning hamella elegans and Rhizopus stolonifer were screened against tyrosinase and cholinesterase showed significant inhibitory activities than the parent compound. Diosgenin and its transformed sarsasapogenin were screened for their acetyl cholinesterase and butyryl cholinesterase inhibitory activities. Sarsasapogenin was screened for phytotoxicity, and was found to be more active than the parent compound. Diosgenin, prednisone and their derivatives were screened for their anti-leishmanial activity. All derivatives were found to be more active than the parent compound. The biotransformation of steroids have been reviewed to a little extent. This review focuses on the biotransformation and functions of selected steroids, the classification, advantages and agents of enzymatic biotransformation and examines the potential role of new enzymatically transformed steroids and their derivatives in the chemoprevention and treatment of other diseases. tyrosinase and cholinesterase inhibitory activities, severe asthma, rheumatic disorders, renal disorders and diseases of inflammatory bowel, skin, gastrointestinal tract.

Organochlorine pesticide acetofenate and its hydrolytic metabolite in rabbits: Enantioselective metabolism and cytotoxicity

Acetofenate (AF) is a chiral organochlorine pesticide used for controlling hygiene pests. In this study, the metabolism of AF in rabbits in vivo and in vitro was investigated and the primary chiral metabolite acetofenate-alcohol (AF-A) was analyzed. The cytotoxicity of AF and AF-A was also determined. AF in rabbits in vivo was eliminated so rapidly that AF could not be detected within 10min after intravenous administration at 20mg/kg (body weight), and AF-A was quickly formed. In vitro metabolism assay, using plasma and liver microsomes, showed that AF was also quickly metabolized to AF-A and the metabolic process was significantly enantioselective with preferential degradation of (-)-AF and formation of (-)-AF-A. The cytotoxicity of AF and AF-A were investigated by assessing cell proliferation, apoptosis and generation of reactive oxygen species. The results showed that AF and AF-A induce enantioselective cytotoxicity. This study will be helpful for improving knowledge about the metabolism and toxicity of AF on an enantiomeric level and providing evidence to understand the potential environmental risk.

Biotransformation of the Flame Retardant 1,2-Dibromo-4-(1,2-dibromoethyl)cyclohexane (TBECH) in Vitro by Human Liver Microsomes

The technical mixture of 1,2-dibromo-4-(1,2-dibromoethyl)cyclohexane (TBECH or DBE-DBCH) and the pure β-TBECH isomer were subjected to in vitro biotransformation by human liver microsomes (HLM). After 60 min of incubation, 5 potential metabolites of TBECH were identified in microsomal assays of both the TBECH mixture and β-TBECH using ultraperformance liquid chromatography-Q-Exactive Orbitrap mass spectrometry. These include mono- and dihydroxylated TBECH and mono- and dihydroxylated TriBECH as well as an α-oxidation metabolite bromo-(1,2-dibromocyclohexyl)-acetic acid. Our results indicate potential hepatic biotransformation of TBECH via cyctochrome P450-catalyzed hydroxylation, debromination, and α-oxidation. Kinetic studies revealed that the formation of monohydroxy-TBECH, dihydroxy-TBECH, and monohydroxy-TriBECH were best fitted to a Michaelis-Menten enzyme kinetic model. Respective estimated V_max values (maximum metabolic rate) for these metabolites were 11.8 ± 4, 0.6 ± 0.1, and 10.1 ± 0.8 pmol min^-1 mg protein^-1 in TBECH mixture and 4992 ± 1340, 14.1 ± 4.9, and 66.1 ± 7.3 pmol min^-1 mg protein^-1 in β-TBECH. This indicates monohydroxy-TBECH as the major metabolite of TBECH by in vitro HLM-based assay. The estimated in vitro intrinsic clearance (Cl_int) of TBECH mixture was slower (P < 0.05) than that of pure β-TBECH. While the formation of monohydroxy-TBECH may reduce the bioaccumulation potential and provide a useful biomarker for monitoring TBECH exposure, further studies are required to fully understand the levels and toxicological implications of the identified metabolites.

Etoxazole is Metabolized Enantioselectively in Liver Microsomes of Rat and Human in Vitro

Acaricide etoxazole belongs to the ovicides/miticides diphenyloxazole class, affecting adults to lay sterile eggs by inhibiting chitin biosynthesis possibly. The reverse-phase HPLC-MS/MS method was used to determine the etoxazole enantiomers. The enantioselective degradation behavior of rac-etoxazole in liver microsomes of rat and human in vitro with NADPH was dramatically different. The t1/2 of (R)-etoxazole was 15.23 min in rat liver microsomes and 30.54 min in human liver microsomes, while 21.73 and 23.50 min were obtained for (S)-etoxazole, respectively. The Vmax of (R)-etoxazole was almost 5-fold of (S)-etoxazole in liver microsomes of rat in vitro. However, the Vmax of (S)-etoxazole was almost 2-fold of (R)-etoxazole in liver microsomes of human in vitro. The CLint of etoxazole was also shown the enantioselectivity on the contrary in liver microsomes of rat and human. These results indicated that the metabolism of two etoxazole enantiomers was selective in liver microsomes of rat and human in vitro, and enantioselectivity in the two kinds of liver microsomes was in the difference in degradation performance. The reason might be related to the composition and content involved in the enzyme system.

Animal-Free Chemical Safety Assessment

The exponential growth of the Internet of Things and the global popularity and remarkable decline in cost of the mobile phone is driving the digital transformation of medical practice. The rapidly maturing digital, non-medical world of mobile (wireless) devices, cloud computing and social networking is coalescing with the emerging digital medical world of omics data, biosensors and advanced imaging which offers the increasingly realistic prospect of personalized medicine. Described as a potential “seismic” shift from the current “healthcare” model to a “wellness” paradigm that is predictive, preventative, personalized and participatory, this change is based on the development of increasingly sophisticated biosensors which can track and measure key biochemical variables in people. Additional key drivers in this shift are metabolomic and proteomic signatures, which are increasingly being reported as pre-symptomatic, diagnostic and prognostic of toxicity and disease. These advancements also have profound implications for toxicological evaluation and safety assessment of pharmaceuticals and environmental chemicals. An approach based primarily on human in vivo and high-throughput in vitro human cell-line data is a distinct possibility. This would transform current chemical safety assessment practice which operates in a human “data poor” to a human “data rich” environment. This could also lead to a seismic shift from the current animal-based to an animal-free chemical safety assessment paradigm.

Comparison of the inhibitory effects of tolcapone and entacapone against human UDP-glucuronosyltransferases

Tolcapone and entacapone are two potent catechol-O-methyltransferase (COMT) inhibitors with a similar skeleton and displaying similar pharmacological activities. However, entacapone is a very safe drug used widely in the treatment of Parkinson's disease, while tolcapone is only in limited use for Parkinson's patients and needs careful monitoring of hepatic functions due to hepatotoxicity. This study aims to investigate and compare the inhibitory effects of entacapone and tolcapone on human UDP-glucosyltransferases (UGTs), as well as to evaluate the potential risks from the view of drug-drug interactions (DDI). The results demonstrated that both tolcapone and entacapone exhibited inhibitory effects on UGT1A1, UGT1A7, UGT1A9 and UGT1A10. In contrast to entacapone, tolcapone exhibited more potent inhibitory effects on UGT1A1, UGT1A7, and UGT1A10, while their inhibitory potentials against UGT1A9 were comparable. It is noteworthy that the inhibition constants (Ki) of tolcapone and entacapone against bilirubin-O-glucuronidation in human liver microsomes (HLM) are determined as 0.68μM and 30.82μM, respectively, which means that the inhibition potency of tolcapone on UGT1A1 mediated bilirubin-O-glucuronidation in HLM is much higher than that of entacapone. Furthermore, the potential risks of tolcapone or entacapone via inhibition of human UGT1A1 were quantitatively predicted by the ratio of the areas under the plasma drug concentration-time curve (AUC). The results indicate that tolcapone may result in significant increase in AUC of bilirubin or the drugs primarily metabolized by UGT1A1, while entacapone is unlikely to cause a significant DDI through inhibition of UGT1A1.

Enantioselective Metabolism of Flufiprole in Rat and Human Liver Microsomes

The enantioselective metabolism of flufiprole in rat and human liver microsomes in vitro was investigated in this study. The separation and determination were performed using a liquid chromatography system equipped with a triple-quadrupole mass spectrometer and a Lux Cellulose-2 chiral column. The enantioselective metabolism of rac-flufiprole was dramatically different in rat and human liver microsomes in the presence of the β-nicotinamide adenine dinucleotide phosphate regenerating system. The half-lives (t1/2) of flufiprole in rat and human liver microsomes were 7.22 and 21.00 min, respectively, for R-(+)-flufiprole, whereas the values were 11.75 and 17.75 min, respectively, for S-(-)-flufiprole. In addition, the Vmax of R-(+)-flufiprole was about 3-fold that of S-(-)-flufiprole in rat liver microsomes, whereas its value in the case of S-(-)-flufiprole was about 2-fold that of R-(+)-flufiprole in human liver microsomes. The CLint of rac-flufiprole also showed opposite enantioselectivy in rat and human liver microsomes. The different compositions and contents of metabolizing enzyme in the two liver microsomes might be the reasons for the difference in the metabolic behavior of the two enantiomers.

QSARs for estimating intrinsic hepatic clearance of organic chemicals in humans

Quantitative structure-activity relationships (QSARs) were developed to predict the in vitro clearance (CLINT) of xenobiotics metabolised in human hepatocytes (118 compounds) and microsomes (115 compounds). Clearance values were gathered from the scientific literature and multiple linear models were built and validated selecting at most 6 predictors from a pool of over 2000 potential molecular descriptors. For the hepatocytes QSAR, the explained variance (Radj(2)) was 67% and the predictive ability (Rext(2)) was 62%. For the microsomes QSAR, Radj(2) was 50% and Rext(2) 30%. For both liver assays, the most important descriptor relates to electronic properties of the compound. Functional groups of fragments were useful to identify specific compounds that have a deviating reaction rate compared to the others, such as polychlorobiphenyls (PCBs) and organic amides which were poorly metabolised by hepatocytes and microsomes, respectively. For hepatocytes, clearance was predominantly determined by electronic characteristics, while size and shape characteristics were less important and partitioning properties were absent. This may suggest that uptake across the membrane and enzyme binding are not rate-limiting steps. Particularly for hepatocytes the QSAR statistics are encouraging, allowing application of the outcomes in in vitro to in vivo extrapolation.

Kinetics of tris (1-chloro-2-propyl) phosphate (TCIPP) metabolism in human liver microsomes and serum

Tris(1-chloro-2-propyl) phosphate (TCIPP) is an emerging contaminant which is ubiquitous in the indoor and outdoor environment. Moreover, its presence in human body fluids and biota has been evidenced. Since no quantitative data exist on the biotransformation or stability of TCIPP in the human body, we performed an in vitro incubation of TCIPP with human liver microsomes (HLM) and human serum (HS). Two metabolites, namely bis(2-chloro-isopropyl) phosphate (BCIPP) and bis(1-chloro-2-propyl) 1-hydroxy-2-propyl phosphate (BCIPHIPP), were quantified in a kinetic study using HLM or HS (only BCIPP, the hydrolysis product) and LC-MS. The Michaelis-Menten model fitted best the NADPH-dependent formation of BCIPHIPP and BCIPP in HLM, with respective V(MAX) of 154 ± 4 and 1470 ± 110 pmol/min/mg protein and respective apparent K(m) of 80.2 ± 4.4 and 96.1 ± 14.5 μM. Hydrolases, which are naturally present in HLM, were also involved in the production of BCIPP. A HS paraoxonase assay could not detect any BCIPP formation above 38.6 ± 10.8 pmol/min/μL serum. Our data indicate that BCIPP is the major metabolite of TCIPP formed in the liver. To our knowledge, this is the first quantitative assessment of the stability of TCIPP in tissues of humans or any other species. Further research is needed to confirm whether these biotransformation reactions are associated with a decrease or increase in toxicity.

Development and application of a human PBPK model for bromodichloromethane to investigate the impacts of multi-route exposure

As a result of its presence in water as a volatile disinfection byproduct, bromodichloromethane (BDCM), which is mutagenic, poses a potential health risk from exposure via oral, dermal and inhalation routes. We developed a refined human physiologically based pharmacokinetic (PBPK) model for BDCM (including new chemical-specific human parameters) to evaluate the impact of BDCM exposure during showering and bathing on important measures of internal dose compared with oral exposure. The refined model adequately predicted data from the published literature for oral, dermal and bathing/showering exposures. A liter equivalency approach (L-eq) was used to estimate BDCM concentration in a liter of water consumed by the oral route that would be required to produce the same internal dose of BDCM resulting from a 20-min bath or a 10-min shower in water containing 10 µg l(-1) BDCM. The oral liter equivalent concentrations for the bathing scenario were 605, 803 and 5 µg l(-1) BDCM for maximum venous blood concentration (Cmax), the area under the curve (AUCv) and the amount metabolized in the liver per hour (MBDCM), respectively. For a 10-min showering exposure, the oral L-eq concentrations were 282, 312 and 2.1 µg l(-1) for Cmax, AUC and MBDCM, respectively. These results demonstrate large contributions of dermal and inhalation exposure routes to the internal dose of parent chemical reaching the systemic circulation, which could be transformed to mutagenic metabolites in extrahepatic target tissues. Thus, consideration of the contribution of multiple routes of exposure when evaluating risks from water-borne BDCM is needed, and this refined human model will facilitate improved assessment of internal doses from real-world exposures. Published 2015. This article has been contributed to by US Government employees and their work is in the public domain in the USA.

Toxicokinetic Model Development for the Insensitive Munitions Component 2,4-Dinitroanisole

The Armed Forces are developing new explosives that are less susceptible to unintentional detonation (insensitive munitions [IMX]). 2,4-Dinitroanisole (DNAN) is a component of IMX. Toxicokinetic data for DNAN are required to support interpretation of toxicology studies and refinement of dose estimates for human risk assessment. Male Sprague-Dawley rats were dosed by gavage (5, 20, or 80 mg DNAN/kg), and blood and tissue samples were analyzed to determine the levels of DNAN and its metabolite 2,4-dinitrophenol (DNP). These data and data from the literature were used to develop preliminary physiologically based pharmacokinetic (PBPK) models. The model simulations indicated saturable metabolism of DNAN in rats at higher tested doses. The PBPK model was extrapolated to estimate the toxicokinetics of DNAN and DNP in humans, allowing the estimation of human-equivalent no-effect levels of DNAN exposure from no-observed adverse effect levels determined in laboratory animals, which may guide the selection of exposure limits for DNAN.

The application of global sensitivity analysis in the development of a physiologically based pharmacokinetic model for m-xylene and ethanol co-exposure in humans

Global sensitivity analysis (SA) was used during the development phase of a binary chemical physiologically based pharmacokinetic (PBPK) model used for the analysis of m-xylene and ethanol co-exposure in humans. SA was used to identify those parameters which had the most significant impact on variability of venous blood and exhaled m-xylene and urinary excretion of the major metabolite of m-xylene metabolism, 3-methyl hippuric acid. This analysis informed the selection of parameters for estimation/calibration by fitting to measured biological monitoring (BM) data in a Bayesian framework using Markov chain Monte Carlo (MCMC) simulation. Data generated in controlled human studies were shown to be useful for investigating the structure and quantitative outputs of PBPK models as well as the biological plausibility and variability of parameters for which measured values were not available. This approach ensured that a priori knowledge in the form of prior distributions was ascribed only to those parameters that were identified as having the greatest impact on variability. This is an efficient approach which helps reduce computational cost.

In vitro biotransformation of tris(2-butoxyethyl) phosphate (TBOEP) in human liver and serum

Tris(2-butoxyethyl) phosphate (TBOEP) is a plasticizer present in indoor dust, reaching levels of several micrograms per gram. Such levels could lead to significant daily exposure of adults and children. Currently, no toxicokinetic data are available to estimate TBOEP clearance in humans after uptake and therefore, one objective of this study was to investigate intrinsic clearance of TBOEP by human liver microsome (HLM) and serum enzymes. Another objective was to generate information to identify and prioritize several metabolites of TBOEP for investigation of human exposure by biomonitoring. 1D and 2D-NMR methodologies were successfully applied on a mixture of the metabolites to confirm the structure of 3-HO-TBOEP (bis(2-butoxyethyl) 3-hydroxyl-2-butoxyethyl phosphate) and to tentatively assign structures to 1-HO-TBOEP and 2-HO-TBOEP. HO-TBOEP isomers and bis(2-butoxyethyl) phosphate (BBOEP), bis(2-butoxyethyl) hydroxyethyl phosphate (BBOEHEP) were further monitored by liquid chromatography-tandem mass spectrometry. Rates of formation of BBOEHEP and HO-TBOEP metabolites by liver enzymes were best described by the Michaelis-Menten model. Apparent Km values for BBOEHEP, 3-HO-TBOEP, and sum of 1- and 2-HO-TBOEP isomer formation were 152, 197 and 148μM, respectively. Apparent Vmax values for the formation of BBOEHEP, 3-HO-TBOEP, and the sum of 1- and 2-HO-TBOEP isomers were 2560, 643, and 254pmol/min/mg protein, respectively. No detectable formation of BBOEP occurred with liver or serum enzymes. Our findings indicate that intrinsic clearance of TBOEP is mainly catalyzed by oxidative enzymes in the liver and that its major in vitro metabolite is BBOEHEP. These findings can be applied in human biomonitoring studies and risk assessment.

The utilisation of structural descriptors to predict metabolic constants of xenobiotics in mammals

Quantitative structure-activity relationships (QSARs) were developed to predict the Michaelis-Menten constant (Km) and the maximum reaction rate (Vmax) of xenobiotics metabolised by four enzyme classes in mammalian livers: alcohol dehydrogenase (ADH), aldehyde dehydrogenase (ALDH), flavin-containing monooxygenase (FMO), and cytochrome P450 (CYP). Metabolic constants were gathered from the literature and a genetic algorithm was employed to select at most six predictors from a pool of over 2000 potential molecular descriptors using two-thirds of the xenobiotics in each enzyme class. The resulting multiple linear models were cross-validated using the remaining one-third of the compounds. The explained variances (R(2)adj) of the QSARs were between 50% and 80% and the predictive abilities (R(2)ext) between 50% and 60%, except for the Vmax QSAR of FMO with both R(2)adj and R(2)ext less than 30%. The Vmax values of FMO were independent of substrate chemical structure because the rate-limiting step of its catalytic cycle occurs before compound oxidation. For the other enzymes, Vmax was predominantly determined by functional groups or fragments and electronic properties because of the strong and chemical-specific interactions involved in the metabolic reactions. The most relevant predictors for Km were functional groups or fragments for the enzymes metabolising specific compounds (ADH, ALDH and FMO) and size and shape properties for CYP, likely because of the broad substrate specificity of CYP enzymes. The present study can be helpful to predict the Km and Vmax of four important oxidising enzymes in mammals and better understand the underlying principles of chemical transformation by liver enzymes.

Assessing human variability in kinetics for exposures to multiple environmental chemicals: a physiologically based pharmacokinetic modeling case study with dichloromethane, benzene, toluene, ethylbenzene, and m-xylene

The objective of this study was to compare the magnitude of interindividual variability in internal dose for inhalation exposure to single versus multiple chemicals. Physiologically based pharmacokinetic models for adults (AD), neonates (NEO), toddlers (TODD), and pregnant women (PW) were used to simulate inhalation exposure to "low" (RfC-like) or "high" (AEGL-like) air concentrations of benzene (Bz) or dichloromethane (DCM), along with various levels of toluene alone or toluene with ethylbenzene and xylene. Monte Carlo simulations were performed and distributions of relevant internal dose metrics of either Bz or DCM were computed. Area under the blood concentration of parent compound versus time curve (AUC)-based variability in AD, TODD, and PW rose for Bz when concomitant "low" exposure to mixtures of increasing complexities occurred (coefficient of variation (CV) = 16-24%, vs. 12-15% for Bz alone), but remained unchanged considering DCM. Conversely, AUC-based CV in NEO fell (15 to 5% for Bz; 12 to 6% for DCM). Comparable trends were observed considering production of metabolites (AMET), except for NEO's CYP2E1-mediated metabolites of Bz, where an increased CV was observed (20 to 71%). For "high" exposure scenarios, Cmax-based variability of Bz and DCM remained unchanged in AD and PW, but decreased in NEO (CV= 11-16% to 2-6%) and TODD (CV= 12-13% to 7-9%). Conversely, AMET-based variability for both substrates rose in every subpopulation. This study analyzed for the first time the impact of multiple exposures on interindividual variability in toxicokinetics. Evidence indicates that this impact depends upon chemical concentrations and biochemical properties, as well as the subpopulation and internal dose metrics considered.

Incorporating population variability and susceptible subpopulations into dosimetry for high-throughput toxicity testing

Momentum is growing worldwide to use in vitro high-throughput screening (HTS) to evaluate human health effects of chemicals. However, the integration of dosimetry into HTS assays and incorporation of population variability will be essential before its application in a risk assessment context. Previously, we employed in vitro hepatic metabolic clearance and plasma protein binding data with in vitro in vivo extrapolation (IVIVE) modeling to estimate oral equivalent doses, or daily oral chemical doses required to achieve steady-state blood concentrations (Css) equivalent to media concentrations having a defined effect in an in vitro HTS assay. In this study, hepatic clearance rates of selected ToxCast chemicals were measured in vitro for 13 cytochrome P450 and five uridine 5'-diphospho-glucuronysyltransferase isozymes using recombinantly expressed enzymes. The isozyme-specific clearance rates were then incorporated into an IVIVE model that captures known differences in isozyme expression across several life stages and ethnic populations. Comparison of the median Css for a healthy population against the median or the upper 95th percentile for more sensitive populations revealed differences of 1.3- to 4.3-fold or 3.1- to 13.1-fold, respectively. Such values may be used to derive chemical-specific human toxicokinetic adjustment factors. The IVIVE model was also used to estimate subpopulation-specific oral equivalent doses that were directly compared with subpopulation-specific exposure estimates. This study successfully combines isozyme and physiologic differences to quantitate subpopulation pharmacokinetic variability. Incorporation of these values with dosimetry and in vitro bioactivities provides a viable approach that could be employed within a high-throughput risk assessment framework.

Use of in vitro data in developing a physiologically based pharmacokinetic model: Carbaryl as a case study

In vitro-derived information has been increasingly used to support and improve human health risk assessment for exposure to chemicals. Physiologically based pharmacokinetic (PBPK) modeling is a key component in the movement toward in vitro-based risk assessment, providing a tool to integrate diverse experimental data and mechanistic information to relate in vitro effective concentrations to equivalent human exposures. One of the challenges, however, in the use of PBPK models for this purpose has been the need for extensive chemical-specific parameters. With the remarkable advances in in vitro methodologies in recent years, in vitro-derived parameters can now be easily incorporated into PBPK models. In this study we demonstrate an in vitro data based parameterization approach to develop a physiologically based pharmacokinetic and pharmacodynamic (PBPK/PD) model, using carbaryl as a case study. In vitro experiments were performed to provide the chemical-specific pharmacokinetic (PK) and pharmacodynamic (PD) parameters for carbaryl in the PBPK model for this compound. Metabolic clearance and cholinesterase (ChE) interaction parameters for carbaryl were measured in rat and human tissues. These in vitro PK and PD data were extrapolated to parameters in the whole body PBPK model using biologically appropriate scaling. The PBPK model was then used to predict the kinetics and ChE inhibition dynamics of carbaryl in vivo. This case study with carbaryl provides a reasonably successful example of utilizing the in vitro to in vivo extrapolation (IVIVE) approach for PBPK model development. This approach can be applied to other carbamates with an anticholinesterase mode of action as well as to environmental chemicals in general with further refinement of the current shortcomings in the approach. It will contribute to minimizing the need for in vivo human data for PBPK model parameterization and evaluation in human risk assessments.