Development and Application of a Life-Stage Physiologically Based Pharmacokinetic (PBPK) Model to the Assessment of Internal Dose of Pyrethroids in Humans

To address concerns around age-related sensitivity to pyrethroids, a life-stage physiologically based pharmacokinetic (PBPK) model, supported by in vitro to in vivo extrapolation (IVIVE) was developed. The model was used to predict age-dependent changes in target tissue exposure of 8 pyrethroids; deltamethrin (DLM), cis-permethrin (CPM), trans-permethrin, esfenvalerate, cyphenothrin, cyhalothrin, cyfluthrin, and bifenthrin. A single model structure was used based on previous work in the rat. Intrinsic clearance (CLint) of each individual cytochrome P450 or carboxylesterase (CES) enzyme that are active for a given pyrethroid were measured in vitro, then biologically scaled to obtain in vivo age-specific total hepatic CLint. These IVIVE results indicate that, except for bifenthrin, CES enzymes are largely responsible for human hepatic metabolism (>50% contribution). Given the high efficiency and rapid maturation of CESs, clearance of the pyrethroids is very efficient across ages, leading to a blood flow-limited metabolism. Together with age-specific physiological parameters, in particular liver blood flow, the efficient metabolic clearance of pyrethroids across ages results in comparable to or even lower internal exposure in the target tissue (brain) in children than that in adults in response to the same level of exposure to a given pyrethroid (Cmax ratio in brain between 1- and 25-year old = 0.69, 0.93, and 0.94 for DLM, bifenthrin, and CPM, respectively). Our study demonstrated that a life-stage PBPK modeling approach, coupled with IVIVE, provides a robust framework for evaluating age-related differences in pharmacokinetics and internal target tissue exposure in humans for the pyrethroid class of chemicals.

Physiologically Based Pharmacokinetic Modeling in Risk Assessment: Case Study With Pyrethroids

The assessment of potentially sensitive populations is an important application of risk assessment. To address the concern for age-related sensitivity to pyrethroid insecticides, life-stage physiologically based pharmacokinetic (PBPK) modeling supported by in vitro to in vivo extrapolation was conducted to predict age-dependent changes in target tissue exposure to 8 pyrethroids. The purpose of this age-dependent dosimetry was to calculate a Data-derived Extrapolation Factor (DDEF) to address age-related pharmacokinetic differences for pyrethroids in humans. We developed a generic human PBPK model for pyrethroids based on our previously published rat model that was developed with in vivo rat data. The results demonstrated that the age-related differences in internal exposure to pyrethroids in the brain are largely determined by the differences in metabolic capacity and in physiology for pyrethroids between children and adults. The most important conclusion from our research is that, given an identical external exposure, the internal (target tissue) concentration is equal or lower in children than in adults in response to the same level of exposure to a pyrethroid. Our results show that, based on the use of the life-stage PBPK models with 8 pyrethroids, DDEF values are essentially close to 1, resulting in a DDEF for age-related pharmacokinetic differences of 1. For risk assessment purposes, this indicates that no additional adjustment factor is necessary to account for age-related pharmacokinetic differences for these pyrethroids.

Bringing Big Data to Bear in Environmental Public Health: Challenges and Recommendations

Understanding the role that the environment plays in influencing public health often involves collecting and studying large, complex data sets. There have been a number of private and public efforts to gather sufficient information and confront significant unknowns in the field of environmental public health, yet there is a persistent and largely unmet need for findable, accessible, interoperable, and reusable (FAIR) data. Even when data are readily available, the ability to create, analyze, and draw conclusions from these data using emerging computational tools, such as augmented and artificial inteligence (AI) and machine learning, requires technical skills not currently implemented on a programmatic level across research hubs and academic institutions. We argue that collaborative efforts in data curation and storage, scientific computing, and training are of paramount importance to empower researchers within environmental sciences and the broader public health community to apply AI approaches and fully realize their potential. Leaders in the field were asked to prioritize challenges in incorporating big data in environmental public health research: inconsistent implementation of FAIR principles in data collection and sharing, a lack of skilled data scientists and appropriate cyber-infrastructures, and limited understanding of possibilities and communication of benefits were among those identified. These issues are discussed, and actionable recommendations are provided.

A Time-Embedding Network Models the Ontogeny of 23 Hepatic Drug Metabolizing Enzymes

Pediatric patients are at elevated risk of adverse drug reactions, and there is insufficient information on drug safety in children. Complicating risk assessment in children, there are numerous age-dependent changes in the absorption, distribution, metabolism, and elimination of drugs. A key contributor to age-dependent drug toxicity risk is the ontogeny of drug metabolism enzymes, the changes in both abundance and type throughout development from the fetal period through adulthood. Critically, these changes affect not only the overall clearance of drugs but also exposure to individual metabolites. In this study, we introduce time-embedding neural networks in order to model population-level variation in metabolism enzyme expression as a function of age. We use a time-embedding network to model the ontogeny of 23 drug metabolism enzymes. The time-embedding network recapitulates known demographic factors impacting 3A5 expression. The time-embedding network also effectively models the nonlinear dynamics of 2D6 expression, enabling a better fit to clinical data than prior work. In contrast, a standard neural network fails to model these features of 3A5 and 2D6 expression. Finally, we combine the time-embedding model of ontogeny with additional information to estimate age-dependent changes in reactive metabolite exposure. This simple approach identifies age-dependent changes in exposure to valproic acid and dextromethorphan metabolites and suggests potential mechanisms of valproic acid toxicity. This approach may help researchers evaluate the risk of drug toxicity in pediatric populations.

Developmental Expression of the Cytosolic Sulfotransferases in Human Liver

The liver is the predominant organ of metabolism for many endogenous and foreign chemicals. Cytosolic sulfotransferases (SULTs) catalyze the sulfonation of drugs and other xenobiotics, as well as hormones, neurotransmitters, and sterols, with consequences that include enhanced drug elimination, hormone inactivation, and procarcinogen bioactivation. SULTs are classified into six gene families, but only SULT1 and SULT2 enzymes are expressed in human liver. We characterized the developmental expression patterns of SULT1 and SULT2 mRNAs and proteins in human liver samples using reverse transcription quantitative polymerase chain reaction (RT-qPCR), RNA sequencing, and targeted quantitative proteomics. Using a set of prenatal, infant, and adult liver specimens, RT-qPCR analysis demonstrated that SULT1A1 (transcript variant 1) expression did not vary appreciably during development; SULT1C2, 1C4, and 1E1 mRNA levels were highest in prenatal and/or infant liver, and 1A2, 1B1, and 2A1 mRNA levels were highest in infant and/or adult. Hepatic SULT1A1 (transcript variant 5), 1C3, and 2B1 mRNA levels were low regardless of developmental stage. Results obtained with RNA sequencing of a different set of liver specimens (prenatal and pediatric) were generally comparable results to those of the RT-qPCR analysis, with the additional finding that SULT1A3 expression was highest during gestation. Analysis of SULT protein content in a library of human liver cytosols demonstrated that protein levels generally corresponded to the mRNAs, with the major exception that SULT1C4 protein levels were much lower than expected based on mRNA levels. These findings further support the concept that hepatic SULTs play important metabolic roles throughout the human life course, including early development.

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.

A Chemical Category-Based Prioritization Approach for Selecting 75 Per- and Polyfluoroalkyl Substances (PFAS) for Tiered Toxicity and Toxicokinetic Testing

Per- and polyfluoroalkyl substances (PFASs) are a group of fluorinated substances of interest to researchers, regulators, and the public due to their widespread presence in the environment. A few PFASs have comparatively extensive amounts of human epidemiological, exposure, and experimental animal toxicity data (e.g., perfluorooctanoic acid), whereas little toxicity and exposure information exists for much of the broader set of PFASs. Given that traditional approaches to generate toxicity information are resource intensive, new approach methods, including in vitro high-throughput toxicity (HTT) testing, are being employed to inform PFAS hazard characterization and further (in vivo) testing. The U.S. Environmental Protection Agency (EPA) and the National Toxicology Program (NTP) are collaborating to develop a risk-based approach for conducting PFAS toxicity testing to facilitate PFAS human health assessments. This article describes the construction of a PFAS screening library and the process by which a targeted subset of 75 PFASs were selected. Multiple factors were considered, including interest to the U.S. EPA, compounds within targeted categories, structural diversity, exposure considerations, procurability and testability, and availability of existing toxicity data. Generating targeted HTT data for PFASs represents a new frontier for informing priority setting. https://doi.org/10.1289/EHP4555.

Application of epigenetic data in human health risk assessment

Despite the many recent advances in the field of epigenetics, application of this knowledge in environmental health risk assessment has been limited. In this paper, we identify opportunities for application of epigenetic data to support health risk assessment. We consider current applications and present a vision for the future.

Determination of Human Hepatic CYP2C8 and CYP1A2 Age-Dependent Expression to Support Human Health Risk Assessment for Early Ages

Predicting age-specific metabolism is important for evaluating age-related drug and chemical sensitivity. Multiple cytochrome P450s and carboxylesterase enzymes are responsible for human pyrethroid metabolism. Complete ontogeny data for each enzyme are needed to support in vitro to in vivo extrapolation (IVIVE). This study was designed to determine age-dependent human hepatic CYP2C8 expression, for which only limited ontogeny data are available, and to further define CYP1A2 ontogeny. CYP2C8 and 1A2 protein levels were measured by quantitative Western blotting using liver microsomal samples prepared from 222 subjects with ages ranging from 8 weeks gestation to 18 years after birth. The median CYP2C8 expression was significantly greater among samples from subjects older than 35 postnatal days (n = 122) compared with fetal samples and those from very young infants (fetal to 35 days postnatal, n = 100) (0.00 vs. 13.38 pmol/mg microsomal protein; p < 0.0001). In contrast, the median CYP1A2 expression was significantly greater after 15 months postnatal age (n = 55) than in fetal and younger postnatal samples (fetal to 15 months postnatal, n = 167) (0.0167 vs. 2.354 pmol/mg microsomal protein; p < 0.0001). CYP2C8, but not CYP1A2, protein levels significantly correlated with those of CYP2C9, CYP2C19, and CYP3A4 (p < 0.001), consistent with CYP2C8 and CYP1A2 ontogeny probably being controlled by different mechanisms. This study provides key data for the physiologically based pharmacokinetic model-based prediction of age-dependent pyrethroid metabolism, which will be used for IVIVE to support pyrethroid risk assessment for early life stages.

FutureTox III: Bridges for Translation

Future Tox III, a Society of Toxicology Contemporary Concepts in Toxicology workshop, was held in November 2015. Building upon Future Tox I and II, Future Tox III was focused on developing the high throughput risk assessment paradigm and taking the science of in vitro data and in silico models forward to explore the question-what progress is being made to address challenges in implementing the emerging big-data toolbox for risk assessment and regulatory decision-making. This article reports on the outcome of the workshop including 2 examples of where advancements in predictive toxicology approaches are being applied within Federal agencies, where opportunities remain within the exposome and AOP domains, and how collectively the toxicology community across multiple sectors can continue to bridge the translation from historical approaches to Tox21 implementation relative to risk assessment and regulatory decision-making.

Prediction of Warfarin Dose in Pediatric Patients: An Evaluation of the Predictive Performance of Several Models

OBJECTIVES

The objective of this study was to evaluate the performance of pediatric pharmacogenetic-based dose prediction models by using an independent cohort of pediatric patients from a multicenter trial.

METHODS

Clinical and genetic data (CYP2C9 [cytochrome P450 2C9] and VKORC1 [vitamin K epoxide reductase]) were collected from pediatric patients aged 3 months to 17 years who were receiving warfarin as part of standard care at 3 separate clinical sites. The accuracy of 8 previously published pediatric pharmacogenetic-based dose models was evaluated in the validation cohort by comparing predicted maintenance doses to actual stable warfarin doses. The predictive ability was assessed by using the proportion of variance (R(2)), mean prediction error (MPE), and the percentage of predictions that fell within 20% of the actual maintenance dose.

RESULTS

Thirty-two children reached a stable international normalized ratio and were included in the validation cohort. The pharmacogenetic-based warfarin dose models showed a proportion of variance ranging from 35% to 78% and an MPE ranging from -2.67 to 0.85 mg/day in the validation cohort. Overall, the model developed by Hamberg et al showed the best performance in the validation cohort (R(2) = 78%; MPE = 0.15 mg/day) with 38% of the predictions falling within 20% of observed doses.

CONCLUSIONS

Pharmacogenetic-based algorithms provide better predictions than a fixed-dose approach, although an optimal dose algorithm has not yet been developed.

Integration of Life-Stage Physiologically Based Pharmacokinetic Models with Adverse Outcome Pathways and Environmental Exposure Models to Screen for Environmental Hazards

A computational framework was developed to assist in screening and prioritizing chemicals based on their dosimetry, toxicity, and potential exposures. The overall strategy started with contextualizing chemical activity observed in high-throughput toxicity screening (HTS) by mapping these assays to biological events described in Adverse Outcome Pathways (AOPs). Next, in vitro to in vivo (IVIVE) extrapolation was used to convert an in vitro dose to an external exposure level, which was compared with potential exposure levels to derive an AOP-based margins of exposure (MOE). In this study, the framework was applied to estimate MOEs for chemicals that can potentially cause developmental toxicity following a putative AOP for fetal vasculogenesis/angiogenesis. A physiologically based pharmacokinetic (PBPK) model was developed to describe chemical disposition during pregnancy, fetal, neonatal, and infant to adulthood stages. Using this life-stage PBPK model, maternal exposures were estimated that would yield fetal blood levels equivalent to the chemical concentration that altered in vitro activity of selected HTS assays related to the most sensitive vasculogenesis/angiogenesis putative AOP. The resulting maternal exposure estimates were then compared with potential exposure levels using literature data or exposure models to derive AOP-based MOEs.

Expression Patterns of Organic Anion Transporting Polypeptides 1B1 and 1B3 Protein in Human Pediatric Liver

Determining appropriate pharmacotherapy in young children can be challenging due to uncertainties in the development of drug disposition pathways. With knowledge of the ontogeny of drug-metabolizing enzymes and an emerging focus on drug transporters, the developmental pattern of the uptake transporters organic anion transporting polypeptide (OATP) 1B1 and 1B3 was assessed by relative protein quantification using Western blotting in 80 human pediatric liver specimens covering an age range from 9 days to 12 years. OATP1B3 exhibited high expression at birth, which declined over the first months of life, and then increased again in the preadolescent period. In comparison with children 6-12 years of age, the relative protein expression of highly glycosylated (total) OATP1B3 was 235% (357%) in children <3 months of age, 33% (64%) in the age group from 3 months to 2 years, and 50% (59%) in children 2-6 years of age. The fraction of highly glycosylated to total OATP1B3 increased with age, indicating ontogenic processes not only at the transcriptional level but also at the post-translational level. Similar to OATP1B3, OATP1B1 showed high interindividual variability in relative protein expression but no statistically significant difference among the studied age groups.

Role of Chromatin Structural Changes in Regulating Human CYP3A Ontogeny

Variability in drug-metabolizing enzyme developmental trajectories contributes to interindividual differences in susceptibility to chemical toxicity and adverse drug reactions, particularly in the first years of life. Factors linked to these interindividual differences are largely unknown, but molecular mechanisms regulating ontogeny are likely involved. To evaluate chromatin structure dynamics as a likely contributing mechanism, age-dependent changes in modified and variant histone occupancy were evaluated within known CYP3A4 and 3A7 regulatory domains. Chromatin immunoprecipitation using fetal or postnatal human hepatocyte chromatin pools followed by quantitative polymerase chain reaction DNA amplification was used to determine relative chromatin occupancy by modified and variant histones. Chromatin structure representing a poised transcriptional state (bivalent chromatin), indicated by the occupancy by modified histones associated with both active and repressed transcription, was observed for CYP3A4 and most 3A7 regulatory regions in both postnatal and fetal livers. However, the CYP3A4 regulatory regions had significantly greater occupancy by modified histones associated with repressed transcription in the fetal liver. Conversely, some modified histones associated with active transcription exhibited greater occupancy in the postnatal liver. CYP3A7 regulatory regions also had significantly greater occupancy by modified histones associated with repressed transcription in the fetus. The observed occupancy by modified histones is consistent with chromatin structural dynamics contributing to CYP3A4 ontogeny, although the data are less conclusive regarding CYP3A7. Interpretation of the latter data may be confounded by cell-type heterogeneity in the fetal liver.

Age-Dependent Human Hepatic Carboxylesterase 1 (CES1) and Carboxylesterase 2 (CES2) Postnatal Ontogeny

Human hepatic carboxylesterase 1 and 2 (CES1 and CES2) are important for the disposition of ester- and amide-bond-containing pharmaceuticals and environmental chemicals. CES1 and CES2 ontogeny has not been well characterized, causing difficulty in addressing concerns regarding juvenile sensitivity to adverse outcomes associated with exposure to certain substrates. To characterize postnatal human hepatic CES1 and CES2 expression, microsomal and cytosolic fractions were prepared using liver samples from subjects without liver disease (N = 165, aged 1 day to 18 years). Proteins were fractionated, detected, and quantitated by Western blotting. Median microsomal CES1 was lower among samples from subjects younger than 3 weeks (n = 36) compared with the rest of the population (n = 126; 6.27 vs. 17.5 pmol/mg microsomal protein, respectively; P < 0.001; Kruskal-Wallis test). Median cytosolic CES1 expression was lowest among samples from individuals between birth and 3 weeks of age (n = 36), markedly greater among those aged 3 weeks to 6 years (n = 90), and modestly greater still among those older than 6 years (n = 36; median values = 4.7, 15.8, and 16.6 pmol/mg cytosolic protein, respectively; P values < 0.001 and 0.05, respectively; Kruskal-Wallis test). Median microsomal CES2 expression increased across the same three age groups with median values of 1.8, 2.9, and 4.2 pmol/mg microsomal protein, respectively (P < 0.001, both). For cytosolic CES2, only the youngest age group differed from the two older groups (P < 0.001; median values = 1.29, 1.93, 2.0, respectively). These data suggest that infants younger than 3 weeks of age would exhibit significantly lower CES1- and CES2-dependent metabolic clearance compared with older individuals.

Baseline Chromatin Modification Levels May Predict Interindividual Variability in Ozone-Induced Gene Expression

Traditional toxicological paradigms have relied on factors such as age, genotype, and disease status to explain variability in responsiveness to toxicant exposure; however, these are neither sufficient to faithfully identify differentially responsive individuals nor are they modifiable factors that can be leveraged to mitigate the exposure effects. Unlike these factors, the epigenome is dynamic and shaped by an individual's environment. We sought to determine whether baseline levels of specific chromatin modifications correlated with the interindividual variability in their ozone (O3)-mediated induction in an air-liquid interface model using primary human bronchial epithelial cells from a panel of 11 donors. We characterized the relationship between the baseline abundance of 6 epigenetic markers with established roles as key regulators of gene expression-histone H3 lysine 4 trimethylation (H3K4me3), H3K27 acetylation (H3K27ac), pan-acetyl H4 (H4ac), histone H3K27 di/trimethylation (H3K27me2/3), unmodified H3, and 5-hydroxymethylcytosine (5-hmC)-and the variability in the O3-induced expression of IL-8, IL-6, COX2, and HMOX1. Baseline levels of H3K4me3, H3K27me2/3, and 5-hmC, but not H3K27ac, H4ac, and total H3, correlated with the interindividual variability in O3-mediated induction of HMOX1 and COX2. In contrast, none of the chromatin modifications that we examined correlated with the induction of IL-8 and IL-6. From these findings, we propose an "epigenetic seed and soil" model in which chromatin modification states between individuals differ in the relative abundance of specific modifications (the "soil") that govern how receptive the gene is to toxicant-mediated cellular signals (the "seed") and thus regulate the magnitude of exposure-related gene induction.

Oxidative stress-responsive transcription factor NRF2 is not indispensable for the human hepatic Flavin-containing monooxygenase-3 (FMO3) gene expression in HepG2 cells

The flavin-containing monooxygenases (FMOs) are important for the oxidation of a variety of endogenous compounds and xenobiotics. The hepatic expression of FMO3 is highly variable and until recently, it was thought to be uninducible. In this study, human FMO3 gene regulation by the oxidative stress transcription factor, nuclear factor (erythroid-derived 2)-like 2 (NRF2) was examined. Constitutive FMO3 gene expression is repressed in HepG2 cells, thus this cell can be a good model for FMO3 gene regulation studies. Over-expression of NRF2 in HepG2 cells increased NRF2 target gene expression, heme oxygenase-1 (HMOX1) and

NAD(P)H

quinone oxidoreductase-1 (NQO1), but did not alter FMO3 gene expression. Co-transfection studies with NRF2 or its cytosolic regulatory protein, Kelch-like ECH-associated protein 1 (KEAP1), expression vectors, along with FMO3 promoter luciferase reporter constructs of various lengths (5kb or 6kb), did not change FMO3 reporter gene activity significantly. Furthermore, treatment with tert-butyl hydroperoxide (tBHP) and tert-butyl hydroquinone (tBHQ) did not alter FMO3 reporter construct activity. In summary, in vitro results suggest that the transcriptional regulation of FMO3 might not involve the NRF2-KEAP1 regulatory pathway.

Use of the RISK21 roadmap and matrix: human health risk assessment of the use of a pyrethroid in bed netting

The HESI-coordinated RISK21 roadmap and matrix are tools that provide a transparent method to compare exposure and toxicity information and assess whether additional refinement is required to obtain the necessary precision level for a decision regarding safety. A case study of the use of a pyrethroid, "pseudomethrin," in bed netting to control malaria is presented to demonstrate the application of the roadmap and matrix. The evaluation began with a problem formulation step. The first assessment utilized existing information pertaining to the use and the class of chemistry. At each stage of the step-wise approach, the precision of the toxicity and exposure estimates were refined as necessary by obtaining key data which enabled a decision on safety to be made efficiently and with confidence. The evaluation demonstrated the concept of using existing information within the RISK21 matrix to drive the generation of additional data using a value-of-information approach. The use of the matrix highlighted whether exposure or toxicity required further investigation and emphasized the need to address the default uncertainty factor of 100 at the highest tier of the evaluation. It also showed how new methodology such as the use of in vitro studies and assays could be used to answer the specific questions which arise through the use of the matrix. The matrix also serves as a useful means to communicate progress to stakeholders during an assessment of chemical use.

Risk assessment in the 21st century: roadmap and matrix

Abstract The RISK21 integrated evaluation strategy is a problem formulation-based exposure-driven risk assessment roadmap that takes advantage of existing information to graphically represent the intersection of exposure and toxicity data on a highly visual matrix. This paper describes in detail the process for using the roadmap and matrix. The purpose of this methodology is to optimize the use of prior information and testing resources (animals, time, facilities, and personnel) to efficiently and transparently reach a risk and/or safety determination. Based on the particular problem, exposure and toxicity data should have sufficient precision to make such a decision. Estimates of exposure and toxicity, bounded by variability and/or uncertainty, are plotted on the X- and Y-axes of the RISK21 matrix, respectively. The resulting intersection is a highly visual representation of estimated risk. Decisions can then be made to increase precision in the exposure or toxicity estimates or declare that the available information is sufficient. RISK21 represents a step forward in the goal to introduce new methodologies into 21st century risk assessment. Indeed, because of its transparent and visual process, RISK21 has the potential to widen the scope of risk communication beyond those with technical expertise.

A 21st century roadmap for human health risk assessment

The Health and Environmental Sciences Institute (HESI)-coordinated Risk Assessment in the 21st Century (RISK21) project was initiated to develop a scientific, transparent, and efficient approach to the evolving world of human health risk assessment, and involved over 120 participants from 12 countries, 15 government institutions, 20 universities, 2 non-governmental organizations, and 12 corporations. This paper provides a brief overview of the tiered RISK21 framework called the roadmap and risk visualization matrix, and articulates the core principles derived by RISK21 participants that guided its development. Subsequent papers describe the roadmap and matrix in greater detail. RISK21 principles include focusing on problem formulation, utilizing existing information, starting with exposure assessment (rather than toxicity), and using a tiered process for data development. Bringing estimates of exposure and toxicity together on a two-dimensional matrix provides a clear rendition of human safety and risk. The value of the roadmap is its capacity to chronicle the stepwise acquisition of scientific information and display it in a clear and concise fashion. Furthermore, the tiered approach and transparent display of information will contribute to greater efficiencies by calling for data only as needed (enough precision to make a decision), thus conserving animals and other resources.

Raising the bar for reproducible science at the U.S. Environmental Protection Agency Office of Research and Development

Considerable concern has been raised regarding research reproducibility both within and outside the scientific community. Several factors possibly contribute to a lack of reproducibility, including a failure to adequately employ statistical considerations during study design, bias in sample selection or subject recruitment, errors in developing data inclusion/exclusion criteria, and flawed statistical analysis. To address some of these issues, several publishers have developed checklists that authors must complete. Others have either enhanced statistical expertise on existing editorial boards, or formed distinct statistics editorial boards. Although the U.S. Environmental Protection Agency, Office of Research and Development, already has a strong Quality Assurance Program, an initiative was undertaken to further strengthen statistics consideration and other factors in study design and also to ensure these same factors are evaluated during the review and approval of study protocols. To raise awareness of the importance of statistical issues and provide a forum for robust discussion, a Community of Practice for Statistics was formed in January 2014. In addition, three working groups were established to develop a series of questions or criteria that should be considered when designing or reviewing experimental, observational, or modeling focused research. This article describes the process used to develop these study design guidance documents, their contents, how they are being employed by the Agency’s research enterprise, and expected benefits to Agency science. The process and guidance documents presented here may be of utility for any research enterprise interested in enhancing the reproducibility of its science.

Human hepatic UGT2B15 developmental expression

Human hepatic UGT2B15 developmental expression changes may alter the metabolism of important drugs and toxicants such as bisphenol A (BPA). Previously, UGT2B15 ontogeny knowledge consisted of transcript data, a dubious surrogate for protein expression. Herein, UGT2B15 protein content was determined in human hepatic microsomes (n = 236, 8 weeks gestation to 18 years). The impact of a common, functional single nucleotide polymorphism (g.253G>T), present in UGT2B15*2 and *5 alleles, was also tested. UGT2B15 expression began during late fetal life, at about 18% of mature values (medians = 48, 267 pmoles/mg of microsomal protein, respectively; p < 0.001). UGT2B15 neonatal (n = 39) and late fetal (≥28 weeks, n = 10) content was similar, but lower than that of infants between 3 and 15 weeks age (n = 46; medians = 38, 48, 404 pmoles/mg microsomal protein, respectively; p < 0.001). Values for the latter group were higher compared with the remaining age group (15 weeks to 18 years; n = 82, p < 0.001). UGT2B15 expression varied 31-fold across the entire sample, and within groups, ranged from 4- to 27-fold. Among postnatal samples, age group, the presence of g.253T and male gender were each significantly associated with greater UGT2B15 expression (p < 0.001, <0.01, and <0.05, respectively; stepwise linear regression). In summary, hepatic UGT2B15 protein onset begins in late gestation; however, the greatest rate of change occurs during the first few weeks after birth. We speculate that the fetus and neonate may have lower clearance of some UGT2B15 substrates, such as BPA, compared with older individuals.

Tolerance to acetaminophen hepatotoxicity in the mouse model of autoprotection is associated with induction of flavin-containing monooxygenase-3 (FMO3) in hepatocytes

Acetaminophen (APAP) pretreatment with a hepatotoxic dose (400 mg/kg) in mice results in resistance to a second, higher dose (600 mg/kg) of APAP (APAP autoprotection). Recent microarray work by our group showed a drastic induction of liver flavin containing monooxygenase-3 (Fmo3) mRNA expression in our mouse model of APAP autoprotection. The role of liver Fmo3, which detoxifies xenobiotics, in APAP autoprotection is unknown. The purpose of this study was to characterize the gene regulation and protein expression of liver Fmo3 during APAP hepatotoxicity. The functional consequences of Fmo3 induction were also investigated. Plasma and livers were collected from male C57BL/6J mice over a period of 72 h following a single dose of APAP (400 mg/kg) to measure Fmo3 mRNA and protein expression. Although Fmo3 mRNA levels increased significantly following APAP treatment, protein expression changed marginally. In contrast, both Fmo3 mRNA and protein expression were significantly higher in APAP autoprotected livers. Unlike male C57BL/6J mice, female mice have ∼80-times higher constitutive Fmo3 mRNA levels and are highly resistant to APAP hepatotoxicity. Coadministration of APAP with the FMO inhibitor methimazole rendered female mice susceptible to APAP hepatotoxicity, with no changes in susceptibility detected in male mice. Furthermore, a human hepatocyte cell line (HC-04) clone over-expressing human FMO3 showed enhanced resistance to APAP cytotoxicity. Taken together, these findings establish for the first time induction of Fmo3 protein expression and function by xenobiotic treatment. Our results also indicate that Fmo3 expression and function plays a role in protecting the liver from APAP-induced toxicity. Although the mechanism(s) of this protection remains to be elucidated, this work describes a novel protective function for this enzyme.

Toll-like receptor genetic variants are associated with Gram-negative infections in VLBW infants

OBJECTIVE

To test the hypothesis that single-nucleotide polymorphisms (SNPs) in Toll-like receptor (TLR) genes alter susceptibility to bacterial infections and modulate white blood cell (WBC) counts during infections in very low birth weight (VLBW) infants (birth weight <1500 g).

STUDY DESIGN

VLBW infants recruited in a multicenter study were genotyped for nine functional TLR SNPs and associations between SNPs and infection rates examined. WBC counts obtained during infections were compared among infants with and without SNPs.

RESULT

In our cohort (n=408), 90 infants developed bacterial infections. Presence of TLR4 (rs4986790 and rs4986791) variants were associated with Gram-negative (G-ve) infections. Female infants heterozygous for the X-linked IRAK1 (rs1059703) SNP had less G-ve infections. In regression models controlling for confounders, the TLR4 (rs4986790) SNP was associated with increased G-ve infections. The TLR5 (rs5744105) variant was associated with elevated WBC counts during infections.

CONCLUSION

TLR genetic variants can contribute to increased risk of bacterial infections and altered immune responses in VLBW infants.

Hepatobiliary disposition of 17-OHPC and taurocholate in fetal human hepatocytes: a comparison with adult human hepatocytes

Little information is available in the literature regarding the expression and activity of transporters in fetal human liver or cultured cells. A synthetic progesterone structural analog, 17α-hydroxyprogesterone caproate (17-OHPC), is used in the prevention of spontaneous abortion in women with a history of recurrent miscarriage (habitual abortion). 17-OHPC has been reported to traverse the placental barrier and gain access to fetal circulation. In this study, the role of transporters in the disposition of 17-OHPC in fetal and adult human hepatocytes was examined. Progesterone metabolites have been reported to induce trans-inhibition of bile acid transporter, ABCB11. Thus, we investigated the effect of 17-OHPC or its metabolites on [(3)H]taurocholic acid transport in sandwich-cultured human fetal and adult hepatocytes. 17-OHPC was taken up rapidly into the cells and transported out partially by an active efflux process that was significantly inhibited by cold temperature, cyclosporine, verapamil, and rifampin. The active efflux mechanism was observed in both adult and fetal hepatocyte cultures. 17-OHPC produced a concentration-dependent inhibition of taurocholate efflux into canaliculi in sandwich-cultured adult and fetal human hepatocytes. However, given the high concentrations required to cause inhibition of these transport processes, no adverse effects would be anticipated from therapeutic levels of 17-OHPC. We also evaluated the expression of various hepatic transporters (ABCB1, ABCB4, SLCO1B1, SLCO1B3, SLCO2B1, ABCB11, SLC10A1, ABCC2, ABCC3, ABCC4, and ABCG2) in fetal and adult hepatocytes. With the exception of ABCB4, all transporters examined were expressed, albeit at lower mRNA levels in fetal hepatocytes compared with adults.

A TLR5 (g.1174C > T) variant that encodes a stop codon (R392X) is associated with bronchopulmonary dysplasia

Current evidence supports a major role for inherited factors in determining bronchopulmonary dysplasia (BPD) susceptibility. The Toll-like receptor (TLR) family of proteins maintain pulmonary homeostasis in the developing lung by aiding pathogen recognition and clearance, regulating inflammation, and facilitating reparative tissue growth. We hypothesized that sequence variation in the TLR pathway genes would alter the susceptibility/severity of BPD in preterm infants. Very low birth-weight infants were recruited prospectively in a multi-center study involving collection of blood samples and clinical information. Nine TLR pathway single-nucleotide polymorphisms were genotyped using a multiplexed single-base extension assay. BPD outcomes were compared among infants with and without the variant allele using Chi-square or Fisher's exact tests. In our cohort (n = 289), 66 (23.6%) infants developed BPD, out of which 32 (11.2%) developed severe BPD. The TLR5 (g.1174C > T) variant was associated with BPD (P = 0.03) and severe BPD (P = 0.004). The TIRAP (g.2054C > T) variant was associated with BPD (P = 0.04). Infants heterozygous for the X-linked IRAK1 (g.6435T > C) variant had a lower incidence of BPD compared to infants homozygous for either the reference or variant allele (P = 0.03). In regression models that controlled for potential epidemiological confounders, the TIRAP variant was associated with BPD, and the TLR5 variant was associated with severe BPD. Our data support the hypothesis that aberrant pathogen recognition in premature infants arising from TLR pathway genetic variation can contribute to BPD pathogenesis.

Prenatal and postnatal expression of glutathione transferase ζ 1 in human liver and the roles of haplotype and subject age in determining activity with dichloroacetate

Glutathione transferase ζ 1 (GSTZ1), also known as maleylacetoacetate isomerase, catalyzes the penultimate step of tyrosine catabolism and metabolizes several α-halocarboxylic acids, including dichloroacetic acid (DCA), an investigational drug used for lactic acidosis and, recently, solid tumors. Age-related differences have been suggested in DCA pharmacotoxicology, but no information is available on GSTZ1 ontogeny in humans. Here, we investigated the cytosolic GSTZ1 developmental expression pattern and the influence of haplotype on GSTZ1 activity with DCA by using human livers from donors between 10 weeks gestation and 74 years. GSTZ1 expression was very low in fetal livers (<2 pmol of GSTZ1/mg cytosol). The expression began to increase after birth in an age-dependent manner until age 7 years. GSTZ1 was then sustained at stable, yet variable, levels (median, 20.0 pmol/mg cytosol; range, 4.8-47.3 pmol/mg cytosol) until age 74 years. GSTZ1 activity with DCA was strongly associated with haplotype and expression level. Samples homozygous or heterozygous for GSTZ1A exhibited ∼3-fold higher DCA dechlorinating activity than samples carrying other alleles at a given level of expression. The correlations (r²) between activity and expression were 0.90 and 0.68, respectively, for GSTZ1A carriers (n = 11) and noncarriers (n = 61). GSTZ1 is expressed in mitochondria in addition to cytosol. The GSTZ1A allele exhibited similar effects in the mitochondrial fraction by conferring a higher activity with DCA. In summary, we report a neonatal onset and an age-related increase in GSTZ1 protein expression during human liver development. Haplotype influenced GSTZ1 activity with DCA but not protein expression.

Approaches for assessing risks to sensitive populations: lessons learned from evaluating risks in the pediatric population

Assessing the risk profiles of potentially sensitive populations requires a "tool chest" of methodological approaches to adequately characterize and evaluate these populations. At present, there is an extensive body of literature on methodologies that apply to the evaluation of the pediatric population. The Health and Environmental Sciences Institute Subcommittee on Risk Assessment of Sensitive Populations evaluated key references in the area of pediatric risk to identify a spectrum of methodological approaches. These approaches are considered in this article for their potential to be extrapolated for the identification and assessment of other sensitive populations. Recommendations as to future research needs and/or alternate methodological considerations are also made.

A Tiered Approach to Life Stages Testing for Agricultural Chemical Safety Assessment

A proposal has been developed by the Agricultural Chemical Safety Assessment (ACSA) Technical Committee of the ILSI Health and Environmental Sciences Institute (HESI) for an improved approach to assessing the safety of crop protection chemicals. The goal is to ensure that studies are scientifically appropriate and necessary without being redundant, and that tests emphasize toxicological endpoints and exposure durations that are relevant for risk assessment. The ACSA Life Stages Task Force proposes a tiered approach to toxicity testing that assesses a compound's potential to cause adverse effects on reproduction, and that assesses the nature and severity of effects during development and adolescence, with consideration of the sensitivity of the elderly. While incorporating many features from current guideline studies, the proposed approach includes a novel rat reproduction and developmental study with enhanced endpoints and a rabbit development study. All available data, including toxicokinetics, ADME data, and systemic toxicity information, are considered in the design and interpretation of studies. Compared to existing testing strategies, the proposed approach uses fewer animals, provides information on the young animal, and includes an estimation of human exposure potential for making decisions about the extent of testing required.

The Role of Flavin-Containing Monooxygenase (FMO) in the Metabolism of Tamoxifen and Other Tertiary Amines

Tamoxifen is utilized in breast cancer therapy and in chemoprevention. Tamoxifen may enhance risk for other neoplasias, especially endometrial cancer. The risk:benefit depends on the rate of metabolic activation versus detoxication. Cytochrome P450-dependent alpha-hydroxylation, followed by sulfonation, represents a metabolic activation pathway, producing products capable of covalent DNA adduction. In contrast, tamoxifen N-oxygenation represents a detoxication pathway, with the caveat that N-oxides can be reduced back to the parent amines. The N-oxygenation pathway will be the focus for this review. Dr. David Kupfer pioneered studies on cytochrome P450 and flavin-containing monooxygenase (FMO) tamoxifen metabolism. We collaborated with Dr. Kupfer's laboratory and recently determined that the low level of tamoxifen N-oxide production in human liver microsomes may be explained by the kinetics of FMO1 versus FMO3.

Developmental changes in human liver CYP2D6 expression

Within the human cytochrome P450 family, specific forms show developmental expression patterns that can affect drug clearance, efficacy, and safety. The objective of this study was to use dextromethorphan O-demethylase activity and quantitative Western blotting to identify CYP2D6 developmental expression patterns in a large (n = 222) and developmentally diverse set of pediatric liver samples. Immunodetectable levels of CYP2D6 protein determined for selected samples across all age categories showed a significant correlation with the corresponding dextromethorphan O-demethylase activity. Of gender, ethnicity, postmortem interval, and genotype, only increasing gestational age was associated with CYP2D6 activity and protein content in prenatal samples. In contrast, both age and genotype were associated with CYP2D6 expression in postnatal samples. CYP2D6 expression in liver samples from neonates less than 7 days of age was higher than that observed in first and second trimester samples, but not significantly higher than third trimester fetal samples. In contrast, expression in postnatal samples greater than 7 days of age was substantially higher than that for any earlier age category. Higher CYP2D6 expression also was observed in liver samples from Caucasians versus African Americans. Finally, using phenotype categories inferred from genotype, CYP2D6 activity was higher in postnatal samples predicted to be extensive or intermediate metabolizers versus poor metabolizers. These results suggest that age and genetic determinants of CYP2D6 expression constitute significant determinants of interindividual variability in CYP2D6-dependent metabolism during ontogeny.

The ontogeny of drug metabolism enzymes and implications for adverse drug events

Profound changes in drug metabolizing enzyme (DME) expression occurs during development that impacts the risk of adverse drug events in the fetus and child. A review of our current knowledge suggests individual hepatic DME ontogeny can be categorized into one of three groups. Some enzymes, e.g., CYP3A7, are expressed at their highest level during the first trimester and either remain at high concentrations or decrease during gestation, but are silenced or expressed at low levels within one to two years after birth. SULT1A1 is an example of the second group of DME. These enzymes are expressed at relatively constant levels throughout gestation and minimal changes are observed postnatally. ADH1C is typical of the third DME group that are not expressed or are expressed at low levels in the fetus, usually during the second or third trimester. Substantial increases in enzyme levels are observed within the first one to two years after birth. Combined with our knowledge of other physiological factors during early life stages, knowledge regarding DME ontogeny has permitted the development of robust physiological based pharmacokinetic models and an improved capability to predict drug disposition in pediatric patients. This review will provide an overview of DME developmental expression patterns and discuss some implications of the data with regards to drug therapy. Common themes emerging from our current knowledge also will be discussed. Finally, the review will highlight gaps in knowledge that will be important to advance this field.

Novel CYP2C9 promoter variants and assessment of their impact on gene expression

There are a considerable number of reports identifying and characterizing genetic variants within the CYP2C9 coding region. Much less is known about polymorphic promoter sequences that also might contribute to interindividual differences in CYP2C9 expression. To address this problem, approximately 10,000 base pairs of CYP2C9 upstream information were resequenced using 24 DNA samples from the Coriell Polymorphism Discovery Resource. Thirty-one single-nucleotide polymorphisms (SNPs) were identified; nine SNPs were novel, whereas 22 were reported previously. Using both sequencing and multiplex single-base extension, individual SNP frequencies were determined in 193 DNA samples obtained from unrelated, self-reported Hispanic Americans of Mexican descent, and they were compared with similar data obtained from a non-Latino white cohort. Significant interethnic differences were observed in several SNP frequencies, some of which seemed unique to the Hispanic population. Analysis using PHASE 2.1 inferred nine common (>1%) variant haplotypes, two of which included the g.3608C>T (R144C) CYP2C9(*)2 and two the g.42614A>C (I359L) CYP2C9(*)3 SNPs. Haplotype variants were introduced into a CYP2C9/luciferase reporter plasmid using site-directed mutagenesis, and the impact of the variants on promoter activity assessed by transient expression in HepG2 cells. Both constitutive and pregnane X receptor-mediated inducible activities were measured. Haplotypes 1B, 3A, and 3B each exhibited a 65% decrease in constitutive promoter activity relative to the reference haplotype. Haplotypes 1D and 3B exhibited a 50% decrease and a 40% increase in induced promoter activity, respectively. These data suggest that genetic variation within CYP2C9 regulatory sequences is likely to contribute to differences in CYP2C9 phenotype both within and among different populations.

Ontogeny of human hepatic cytochromes P450

Significant changes in drug-metabolizing enzyme (DME) expression occur during ontogeny. Such changes can have a profound effect on therapeutic efficacy in the fetus and child, as well as the risk for adverse drug reactions. To gain a better understanding of DME ontogeny, enzyme contents for six key cytochromes P450 were measured in 240 human liver samples representing ages from 8 weeks gestation to 18 years. Where possible, both quantitative western blotting and activity assays with probe substrates were performed. Although oversimplified, the DME can be grouped into one of three categories. As typified by CYP3A7, some enzymes are expressed at their highest level during the first trimester and either remain at high concentrations or decrease during gestation and are silenced or expressed at low levels within 1-2 years after birth. These data cause one to query whether these enzymes have an important endogenous function. Representatives of a second group, CYP3A5 and CYP2C19, are expressed at relatively constant levels throughout gestation. Postnatal increases in CYP2C19 are observed within the first year, but not for CYP3A5. CYP2C9, 2E1, and 3A4 are more typical of a third group of enzymes that are not expressed or are expressed at low levels in the fetus with the onset of expression generally in either the second or third trimester. Substantial increases in expression are observed within the first 1-2 years after birth; however, considerable interindividual variability is observed in the immediate postnatal (1-6 months) onset or increase in expression of these enzymes, often resulting in a window of hypervariability.

Mechanisms regulating human FMO3 transcription

Flavin-containing monooxygenases (FMOs) are important oxidative drug metabolizing enzymes. FMO3 is the primary human adult liver FMO enzyme, but is developmentally regulated. FMO3 promoter characterization using in vitro DNA binding assays with HepG2 cell and fetal and adult liver nuclear protein, as well as FMO3/reporter construct transient expression in HepG2 cells, provided evidence for specific mechanisms contributing to both developmental and constitutive adult regulation. NFY, USF1, an unidentified GC box binding protein, and YY1 appear to play major roles regulating constitutive FMO3 transcription, while Pbx(2) as a heterodimer with an unidentified Hox isoform also may contribute to FMO3 developmental expression.

Identification and functional analysis of a novel human CYP2E1 far upstream enhancer

Both transcriptional and post-transcriptional CYP2E1 regulatory mechanisms are known, resulting in 20-fold or greater variation in CYP2E1 expression. To evaluate functional regulatory elements controlling transcription, CYP2E1 promoter constructs were used to make adenovirus vectors containing CYP2E1 promoter-driven luciferase reporters for analyses in both primary human hepatocytes and HepG2 cells. A 1.2-kilobase pair portion of the CYP2E1 promoter was associated with 5- to 10-fold greater luciferase activity. This upstream region contained five direct repeats of 59 base pairs (bp) that increased thymidine kinase-driven luciferase reporter activity in HepG2 cells more than 5-fold, regardless of orientation. Electrophoretic mobility shift assays (EMSAs) identified sequence-specific nuclear protein binding to the 59-bp repeats that was dependent on a 17-bp sequence containing a canonical GATA binding site (WGATAR). Competitive and supershift EMSA identified the participation of GATA4, another GATA family member or GATA-like factor, and a third factor unrelated to the GATA family. Involvement of the tricho-rhino-phalangeal syndrome-1 factor, which also binds a GATA sequence, was eliminated. Rather, competitive EMSA using known binding sequences for the orphan nuclear receptors, steroidogenic factor-1 (or NR5A1), and fetoprotein transcription factor (or NR5A2) implicated an NR5A member in binding a sequence overlapping the canonical GATA. Chromatin immunoprecipitation assay demonstrated in vivo binding of NR5A2 to the enhancer sequence in human hepatocytes. The enhancer sequence is conserved within the human population but seems species-specific. The identification of this novel enhancer and its putative mechanism adds to the complexities of human CYP2E1 regulation.

Identification and Functional Analysis of Common Human Flavin-Containing Monooxygenase 3 Genetic Variants

Flavin-containing monooxygenases (FMOs) are important for the disposition of many therapeutics, environmental toxicants, and nutrients. FMO3, the major adult hepatic FMO enzyme, exhibits significant interindividual variation. Eighteen FMO3 single-nucleotide polymorphism (SNP) frequencies were determined in 202 Hispanics (Mexican descent), 201 African Americans, and 200 non-Latino whites. Using expressed recombinant enzyme with methimazole, trimethylamine, sulindac, and ethylenethiourea, the novel structural variants FMO3 E24D and K416N were shown to cause modest changes in catalytic efficiency, whereas a third novel variant, FMO3 N61K, was essentially devoid of activity. The latter variant was present at an allelic frequency of 5.2% in non-Latino whites and 3.5% in African Americans, but it was absent in Hispanics. Inferring haplotypes using PHASE, version 2.1, the greatest haplotype diversity was observed in African Americans followed by non-Latino whites and Hispanics. Haplotype 2A and 2B, consisting of a hypermorphic promoter SNP cluster (-2650C>G, -2543T>A, and -2177G>C) in linkage with synonymous structural variants was inferred at a frequency of 27% in the Hispanic population, but only 5% in non-Latino whites and African Americans. This same promoter SNP cluster in linkage with one or more hypomorphic structural variant also was inferred in multiple haplotypes at a total frequency of 5.6% in the African-American study group but less than 1% in the other two groups. The sum frequencies of the hypomorphic haplotypes H3 [15,167G>A (E158K)], H5B [-2650C>G, 15,167G>A (E158K), 21,375C>T (N285N), 21,443A>G (E308G)], and H6 [15,167G>A (E158K), 21,375C>T (N285N)] was 28% in Hispanics, 23% in non-Latino whites, and 24% in African Americans.

Epirubicin Glucuronidation and UGT2B7 Developmental Expression

The usefulness of epirubicin in the treatment of adult and childhood malignant diseases is related in part to the potential reduction in cardiac toxicity compared with that of other anthracyclines given at equivalent doses. An important pathway for epirubicin detoxification is UGT2B7-dependent glucuronidation. This study was implemented to provide a preclinical evaluation of the metabolism of epirubicin with respect to age-related changes in epirubicin glucuronidation in pediatric liver microsomes. Rates of epirubicin glucuronidation and levels of UGT2B7 were determined for liver microsomes from four pediatric age categories (n = 32) and one adult age category (n = 8). Both sets of data showed an increase in UGT2B7 activity and content with increasing age. Epirubicin glucuronidation activity in the adult group was statistically higher compared with all pediatric age groups (p < or = 0.01). UGT2B7 expression also was statistically higher in adults compared with children below 11 years of age, with evidence of significant differences in protein levels among the pediatric age categories. A positive correlation (r = 0.68) between UGT2B7 levels and postnatal age was observed, suggesting a progressive increase in UGT2B7 protein expression with increasing age. However, allometric scaling using the (3/4) power rule suggested no difference in activity between any of the pediatric age categories and the adult, although only a single neonatal sample was included in the analysis. In summary, these in vitro data show differences in epirubicin glucuronidation and UGT2B7 content within pediatric age groups and support the use of epirubicin in pediatric patients at least 6 months of age.

Pharmacogenomics and the future of drug therapy

With the completion of the human genome project, many investigators are striving to translate the resulting wealth of new information into new and improved clinical practices. Pharmacogenomics represents one of the most promising of these applications for adult- and pediatric-based therapies. This article provides a historical perspective, but most importantly, uses this background to illustrate important principles of the field. The application of pharmacogenomics to asthma therapy is presented as an example of the current status of pharmacogenomics as it is being applied to an important pediatric health problem. Finally, a discussion of future promises and challenges to the application of pharmacogenomics is presented, including economic and ethical issues.

Enzyme-Mediated Protein Haptenation of Dapsone and Sulfamethoxazole in Human Keratinocytes: I. Expression and Role of Cytochromes P450

Cutaneous drug reactions (CDRs) are among the most common adverse drug reactions and are responsible for numerous minor to life-threatening complications. Several arylamine drugs, such as sulfamethoxazole (SMX) and dapsone (DDS), undergo bioactivation, resulting in adduction to cellular proteins. These adducted proteins may initiate the immune response that ultimately results in a CDR. Recent studies have demonstrated that normal human epidermal keratinocytes (NHEKs) can bioactivate these drugs, resulting in protein haptenation. We sought to identify the enzyme(s) responsible for this bioactivation in NHEKs. Using immunofluorescence confocal microscopy and an adduct-specific enzyme-linked immunosorbent assay (ELISA), we found that N-acetylation of the primary amine of SMX and DDS markedly reduced the level of protein haptenation in NHEKs. Detection of mRNA and/or protein confirmed the presence of CYP3A4, CYP3A5, and CYP2E1 in NHEKs. In contrast, although a faint band suggestive of CYP2C9 protein was detected in one NHEK sample, a CYP2C9 message was not detectable. We also examined the ability of chemical inhibitors of cytochromes P450 (aminobenzotriazole and 1-dichloroethylene) and cyclooxygenase (indomethacin) to reduce protein haptenation when NHEKs were incubated with SMX or DDS by either confocal microscopy or ELISA. These inhibitors did not significantly attenuate protein adduction with either SMX or DDS, indicating that cytochromes P450 and cyclooxygenase do not play important roles in the bioactivation of these xenobiotics in NHEKs and thus suggesting the importance of other enzymes in these cells.

Enzyme-mediated protein haptenation of dapsone and sulfamethoxazole in human keratinocytes: II. Expression and role of flavin-containing monooxygenases and peroxidases

Arylamine compounds, such as sulfamethoxazole (SMX) and dapsone (DDS), are metabolized in epidermal keratinocytes to arylhydroxylamine metabolites that auto-oxidize to arylnitroso derivatives, which in turn bind to cellular proteins and can act as antigens/immunogens. Previous studies have demonstrated that neither cytochromes P450 nor cyclooxygenases mediate this bioactivation in normal human epidermal keratinocytes (NHEKs). In this investigation, we demonstrated that methimazole (MMZ), a prototypical substrate of the flavin-containing monooxygenases (FMOs), attenuated the protein haptenation observed in NHEKs exposed to SMX or DDS. In addition, recombinant FMO1 and FMO3 were able to bioactivate both SMX and DDS, resulting in covalent adduct formation. Western blot analysis confirmed the presence of FMO3 in NHEKs, whereas FMO1 was not detectable. In addition to MMZ, 4-aminobenzoic acid hydrazide (ABH) also attenuated SMX- and DDS-dependent protein haptenation in NHEKs. ABH did not alter the bioactivation of these drugs by recombinant FMO3, suggesting its inhibitory effect in NHEKs was due to its known ability to inhibit peroxidases. Studies confirmed the presence of peroxidase activity in NHEKs; however, immunoblot analysis and reverse transcription-polymerase chain reaction indicated that myeloperoxidase, lactoperoxidase, and thyroid peroxidase were absent. Thus, our results suggest an important role for FMO3 and yet-to-be identified peroxidases in the bioactivation of sulfonamides in NHEKs.

Modeling interchild differences in pharmacokinetics on the basis of subject-specific data on physiology and hepatic CYP2E1 levels: A case study with toluene

The objective of the present study was to evaluate the magnitude of interindividual variability in the internal dose of toluene in children of various age groups, on the basis of subject-specific hepatic CYP2E1 content and physiology. The methodology involved the use of a previously validated physiologically based pharmacokinetic (PBPK) model, in which the intrinsic clearance for hepatic metabolism (CL(int)) was expressed in terms of the CYP2E1 content. The adult toluene PBPK model, with enzyme content-normalized CL(int), facilitated the calculation of child-specific CL(int) based on knowledge of hepatic CYP2E1 protein levels. The child-specific physiological parameters, except liver volume, were computed with knowledge of age and body weight, whereas physicochemical parameters for toluene were kept age-invariant based on available data. The actual individual-specific liver volume (autopsy data) was also included in the model. The resulting model was used to simulate the blood concentration profiles in children exposed by inhalation, to 1 ppm toluene for 24 h. For this exposure scenario, the area under the venous blood concentration vs. time curve (AUC) ranged from 0.30 to 1.01 microg/ml x h in neonates with low CYP2E1 concentration (<3.69 pmol/mg protein). The simulations indicated that neonates with higher levels of CYP2E1 (4.33 to 55.93 pmol/mg protein) as well as older children would have lower AUC (0.16 to 0.43 microg/ml x h). The latter values were closer to those simulated for adults. Similar results were also obtained for 7 h exposure to 17 ppm toluene, a scenario previously evaluated in human volunteers. The interindividual variability factor for each subgroup of children and adults, calculated as the ratio of the 95th and 50th percentile values of AUC, was within a factor of 2. The 95th percentile value of the low metabolizing neonate group, however, was greater than the mean adult AUC by a factor of 3.9. This study demonstrates the feasibility of incorporating subject-specific data on hepatic CYP2E1 content and physiology within PBPK models for evaluating the age, interchild and population variability of internal dose for use in risk assessment of inhaled volatile organics.