The evolution of the OATP hepatic uptake transport protein family in DMPK sciences: from obscure liver transporters to key determinants of hepatobiliary clearance

Over the last two decades the impact on drug pharmacokinetics of the organic anion transporting polypeptides (OATPs: OATP-1B1, 1B3 and 2B1), expressed on the sinusoidal membrane of the hepatocyte, has been increasingly recognized. OATP-mediated uptake into the hepatocyte coupled with subsequent excretion into bile via efflux proteins, such as MRP2, is often referred to as hepatobiliary excretion. OATP transporter proteins can impact some drugs in several ways including pharmacokinetic variability, pharmacodynamic response and drug-drug interactions (DDIs). The impact of transporter mediated hepatic clearance is illustrated with case examples, from the literature and also from the Pfizer portfolio. The currently available in vitro techniques to study the hepatic transporter proteins involved in the hepatobiliary clearance of drugs are reviewed herein along with recent advances in using these in vitro data to predict the human clearance of compounds recognized by hepatic uptake transporters.

Metabolism of myclobutanil and triadimefon by human and rat cytochrome P450 enzymes and liver microsomes

Metabolism of two triazole-containing antifungal azoles was studied using expressed human and rat cytochrome P450s (CYP) and liver microsomes. Substrate depletion methods were used due to the complex array of metabolites produced from myclobutanil and triadimefon. Myclobutanil was metabolized more rapidly than triadimefon, which is consistent with metabolism of the n-butyl side-chain in the former and the t-butyl group in the latter compound. Human and rat CYP2C and CYP3A enzymes were the most active. Metabolism was similar in microsomes prepared from livers of control and low-dose rats. High-dose (115 mg kg-1 day-1 of triadimefon or 150 mg kg-1 day-1 of myclobutanil) rats showed increased liver weight, induction of total CYP, and increased metabolism of the two triazoles, though the apparent Km appeared unchanged relative to the control. These data identify CYP enzymes important for the metabolization of these two triazoles. Estimated hepatic clearances suggest that CYP induction may have limited impact in vivo.

Database for physiologically based pharmacokinetic (PBPK) modeling: physiological data for healthy and health-impaired elderly

Physiologically based pharmacokinetic (PBPK) models have increasingly been employed in chemical health risk assessments. By incorporating individual variability conferred by genetic polymorphisms, health conditions, and physiological changes during development and aging, PBPK models are ideal for predicting chemical disposition in various subpopulations of interest. In order to improve the parameterization of PBPK models for healthy and health-impaired elderly (herein defined as those aged 65 yr and older), physiological parameter values were obtained from the peer-reviewed literature, evaluated, and entered into a Microsoft ACCESS database. Database records include values for key age-specific model inputs such as ventilation rates, organ volumes and blood flows, glomerular filtration rates, and other clearance-related processes. In total, 528 publications were screened for relevant data, resulting in the inclusion of 155 publications comprising 1051 data records for healthy elderly adults and 115 data records for elderly with conditions such as diabetes, chronic obstructive pulmonary disease (COPD), obesity, heart disease, and renal disease. There are no consistent trends across parameters or their associated variance with age; the gross variance in body weight decreased with advancing age, whereas there was no change in variance for brain weight. The database contains some information to inform ethnic and gender differences in parameters; however, the majority of the published data pertain to Asian (mostly Japanese) and Caucasian males. As expected, the number of records tends to decrease with advancing age. In addition to a general lack of data for parameters in the elderly with various health conditions, there is also a dearth of information on blood and tissue composition in all elderly groups. Importantly, there are relatively few records for alveolar ventilation rate; therefore, the relationship between this parameter and cardiac output (usually assumed to be 1:1) in the elderly is not well informed by the database. Despite these limitations, the database represents a potentially useful resource for parameterizing PBPK models for the elderly to facilitate the prediction of dose metrics in older populations for application in risk assessment.

Modeling single and repeated dose pharmacokinetics of PFOA in mice

Perfluorooctanoic acid (PFOA) displays complicated pharmacokinetics in that serum concentrations indicate long half-lives despite which steady state appears to be achieved rapidly. In this study, serum and tissue concentration time-courses were obtained for male and female CD1 mice after single, oral doses of 1 and 10 mg/kg of PFOA. When using one- and two-compartment models, the pharmacokinetics for these two dosages are not consistent with serum time-course data from female CD1 mice administered 60 mg/kg, or with serum concentrations following repeated daily doses of 20 mg/kg PFOA. Some consistency between dose regimens could be achieved using the saturable resorption model of Andersen et al. In this model PFOA is cleared from the serum into a filtrate compartment from which it is either excreted or resorbed into the serum by a process presumed transporter mediated with a Michaelis-Menten form. Maximum likelihood estimation found a transport maximum of T(m) = 860.9 (1298.3) mg/l/h and half-maximum concentration of K(T) = 0.0015 (0.0022) mg/l where the estimated standard errors (in parentheses) indicated large uncertainty. The estimated rate of flow into and out of the filtrate compartment, 0.6830 (1.0131) l/h was too large to be consistent with a biological interpretation. For these model parameters a single dose greater than 40 mg/kg, or a daily dose in excess of 5 mg/kg were necessary to observe nonlinear pharmacokinetics for PFOA in female CD1 mice. These data and modeling analyses more fully characterize PFOA in mice for purposes of estimating internal exposure for use in risk assessment.

Approaches for applications of physiologically based pharmacokinetic models in risk assessment

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

Comparing single and repeated dosimetry data for perfluorooctane sulfonate in rats

Perfluorooctane sulfonate (PFOS) is a member of a class of perfluorinated chemicals used in a variety of consumer and industrial applications because of their oleophobic and hydrophobic properties. It has been shown to cause toxicity in adult and developing laboratory animals. Because PFOS has also been shown to be widely distributed throughout the environment, there have been concerns about its potential health risk to humans. Limited pharmacokinetic data for PFOS are available in rodents and humans, while epidemiological studies of workers and extensive toxicity studies in rodents have been performed. The existing pharmacokinetic and toxicity database in rodents can be useful in the cross-species extrapolations needed to evaluate and interpret internal dosimetry in humans. A mathematical model that describes the disposition of PFOS in adult rats following intravenous, oral, and chronic dietary exposures was developed to gain a better understanding of the pharmacokinetics of PFOS and to determine whether single-dose kinetics are predictive of repeated-dose kinetics. In order to characterize existing time-course data, time-dependent and concentration-dependent changes in the pharmacokinetic parameters for urinary and biliary clearance and liver distribution were needed. Whether these time-dependent changes represent inconsistencies across experiments, effects of aging in the rats, or chemically induced changes in pharmacokinetics remains to be determined.

Predicting maternal rat and pup exposures: how different are they?

Risk and safety assessments for early life exposures to environmental chemicals or pharmaceuticals based on cross-species extrapolation would greatly benefit from information on chemical dosimetry in the young. Although relevant toxicity studies involve exposures during multiple life stages, the mother's exposure dose is frequently used for extrapolation of rodent toxicity findings to humans and represents a substantial source of uncertainty. A compartmental pharmacokinetic model augmented with biological information on factors changing during lactation and early postweaning was developed. The model uses adult pharmacokinetics, milk distribution, and relevant postnatal biology to predict dosimetry in the young for chemicals. The model addressed three dosing strategies employed in toxicity studies (gavage, constant ppm diet, and adjusted ppm diet) and the impact of different pharmacokinetic properties such as rates of clearance, milk distribution, and volume of distribution on the pup exposure doses and internal dosimetry. Developmental delays in clearance and recirculation of chemical in excreta from the pup to mother were evaluated. Following comparison with data for two chemicals, predictions were made for theoretical chemicals with a range of characteristics. Pup exposure was generally lower than the mother's with a shorter half-life, lower milk transfer, larger volume of distribution, and gavage dosing, while higher with longer half-life, higher milk transfer, smaller volume of distribution, and dietary exposures. The present model demonstrated pup exposures do not always parallel the mother's. The model predictions can be used to help design early life toxicity and pharmacokinetic studies and better interpret study findings.

Serologic evaluation of human microcystin exposure

Microcystins are among the most commonly detected toxins associated with cyanobacteria blooms worldwide. Two episodes of intravenous microcystin exposures occurred among kidney dialysis patients during 1996 and 2001. Analysis of serum samples collected during these episodes suggests that microcystins are detectable as free and bound forms in human serum. Our goal was to characterize the biochemical evidence for human exposure to microcystins, to identify uncertainties associated with interpretation of these observed results, and to identify research needs. We analyzed serum samples using enzyme-linked immunosorbent assay (ELISA) methods to detect free microcystins, and gas chromatography/mass spectrometry (GC/MS) to detect 2-methyl-3-methoxy-4-phenylbutyric acid (MMPB). MMPB is derived from both free and protein-bound microcystins by chemical oxidation, and it appears to represent total microcystins present in serum. We found evidence of free microcystins in patient serum for more than 50 days after the last documented exposure. Serum concentrations of free microcystins were consistently lower than MMPB quantification of total microcystins: free microcystins as measured by ELISA were only 8-51% of total microcystin concentrations as detected by the GC/MS method. After intravenous exposure episodes, we found evidence of microcystins in human serum in free and protein-bound forms, though the nature of the protein-bound forms is uncertain. Free microcystins appear to be a small but variable subset of total microcystins present in human serum. Research is needed to elucidate the human toxicokinetics of microcystins, in part to determine how observed serum concentrations can be used to estimate previous microcystin exposure.

Age-dependent partition coefficients for a mixture of volatile organic solvents in Sprague-Dawley rats and humans

The absorption, distribution, metabolism, and excretion of volatile organic compounds (VOCs) are critically determined by a few chemical-specific factors, notably their blood and tissue partition coefficients (PC) and metabolism. Age-specific values for PCs in rats have rarely been reported or utilized in pharmacokinetic modeling for predicting dosimetry in toxicity studies with rats progressing through their lifestages. A mixture of six VOCs (benzene, chloroform, methyl ethyl ketone, methylene chloride, trichloroethylene, and perchloroethylene) was used to determine blood:air and tissue:air PCs in rats at three different ages (postnatal d 10, 60 d, and 2 yr) and blood:air PCs in pediatric and adult human blood. No differences with age in human blood:air PCs for the six compounds were observed. Rat blood:air PCs increased with age varying with compound. Tissue:air PCs showed tissue-specific changes with age. Water-soluble methyl ethyl ketone showed no age-dependent differences. Partition coefficients, particularly the blood:air PC, are key determinants of the rodent and human blood concentrations; age-appropriate values improve the accuracy of pharmacokinetic model predictions of population variability and age-specific exposures.

Extrahepatic metabolism by CYP2E1 in PBPK modeling of lipophilic volatile organic chemicals: impacts on metabolic parameter estimation and prediction of dose metrics

Physiologically based pharmacokinetic (PBPK) models are increasingly available for environmental chemicals and applied in risk assessments. Volatile organic compounds (VOCs) are important pollutants in air, soil, and water. CYP2E1 metabolically activates many VOCs in animals and humans. Despite its presence in extrahepatic tissues, the metabolism by CYP2E1 is often described as restricted to the liver in PBPK models, unless target tissue dose metrics in extrahepatic tissues are needed for the model application, including risk assessment. The impact of accounting for extrahepatic metabolism by CYP2E1 on the estimation of metabolic parameters and the prediction of dose metrics was evaluated for three lipophilic VOCs: vinyl chloride, trichloroethylene, and carbon tetrachloride. Metabolic parameters estimated from fitting gas uptake data with and without extrahepatic metabolism were similar. The impact of extrahepatic metabolism on PBPK predictions was evaluated using inhalation exposure scenarios relevant for animal toxicity studies and human risk assessment. Although small, the relative role of extrahepatic metabolism and the differences in the predicted dose metrics were greater at low exposure concentrations. The impact was species dependent and influenced by Km for CYP2E1. The current study indicates that inhalation modeling for several representative VOCs that are CYP2E1 substrates is not affected by the inclusion of extrahepatic metabolism, implying that liver-only metabolism may be a reasonable simplification for PBPK modeling of lipophilic VOCs. The PBPK predictions using this assumption can be applied confidently for risk assessment, but this conclusion should not necessarily be applied to VOCs that are metabolized by other enzymes.

Characterizing Uncertainty and Variability in Physiologically Based Pharmacokinetic Models: State of the Science and Needs for Research and Implementation

Physiologically based pharmacokinetic (PBPK) models are used in mode-of-action based risk and safety assessments to estimate internal dosimetry in animals and humans. When used in risk assessment, these models can provide a basis for extrapolating between species, doses, and exposure routes or for justifying nondefault values for uncertainty factors. Characterization of uncertainty and variability is increasingly recognized as important for risk assessment; this represents a continuing challenge for both PBPK modelers and users. Current practices show significant progress in specifying deterministic biological models and nondeterministic (often statistical) models, estimating parameters using diverse data sets from multiple sources, using them to make predictions, and characterizing uncertainty and variability of model parameters and predictions. The International Workshop on Uncertainty and Variability in PBPK Models, held 31 Oct-2 Nov 2006, identified the state-of-the-science, needed changes in practice and implementation, and research priorities. For the short term, these include (1) multidisciplinary teams to integrate deterministic and nondeterministic/statistical models; (2) broader use of sensitivity analyses, including for structural and global (rather than local) parameter changes; and (3) enhanced transparency and reproducibility through improved documentation of model structure(s), parameter values, sensitivity and other analyses, and supporting, discrepant, or excluded data. Longer-term needs include (1) theoretical and practical methodological improvements for nondeterministic/statistical modeling; (2) better methods for evaluating alternative model structures; (3) peer-reviewed databases of parameters and covariates, and their distributions; (4) expanded coverage of PBPK models across chemicals with different properties; and (5) training and reference materials, such as cases studies, bibliographies/glossaries, model repositories, and enhanced software. The multidisciplinary dialogue initiated by this Workshop will foster the collaboration, research, data collection, and training necessary to make characterizing uncertainty and variability a standard practice in PBPK modeling and risk assessment.

Predicting Age-Appropriate Pharmacokinetics of Six Volatile Organic Compounds in the Rat Utilizing Physiologically Based Pharmacokinetic Modeling

The capability of physiologically based pharmacokinetic models to incorporate age-appropriate physiological and chemical-specific parameters was utilized to predict changes in internal dosimetry for six volatile organic compounds (VOCs) across different ages of rats. Typical 6-h animal inhalation exposures to 50 and 500 ppm perchloroethylene, trichloroethylene, benzene, chloroform, methylene chloride, or methyl ethyl ketone (MEK) were simulated for postnatal day 10 (PND10), 2-month-old (adult), and 2-year-old (aged) rats. With the exception of MEK, predicted venous blood concentrations of VOCs in the aged rat were equal or up to 1.5-fold higher when compared to the adult rat at both exposure levels, whereas levels were predicted to be up to 3.8-fold higher in the case of PND10 at 50 ppm. Predicted blood levels of MEK were similar in the adult and aged rat, but were more than 5-fold and 30-fold greater for PND10 rats at 500 and 50 ppm, respectively, reflecting high water solubility along with lower metabolic capability and faster ventilation rate per unit body weight (BW) of PND10 animals. Steady-state blood levels of VOCs, simulated by modeling constant exposure, were predicted to be achieved in the order PND10 > adult > aged, largely due to increasing fat volume. The dose metric, total amount metabolized per unit liver volume was generally much lower in PND10 than in adult rats. The blood:air partition coefficient, fat volume, and fat blood flow were identified as critical determinants for the predicted differences in venous blood concentrations between the adult and aged. The lower metabolic capability, largely due to a smaller liver size, and faster ventilation rate per unit BW of PND10 animals contribute the most to the differences between PND10 and adult rats. This study highlights the pharmacokinetic differences and the relevant parameters that may contribute to differential susceptibility to the toxic effects of VOCs across life stages of the rat.

Evaluation of physiologically based pharmacokinetic models for use in risk assessment

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

Mathematical model for the androgenic regulation of the prostate in intact and castrated adult male rats

The testicular-hypothalamic-pituitary axis regulates male reproductive system functions. Understanding these regulatory mechanisms is important for assessing the reproductive effects of environmental and pharmaceutical androgenic and antiandrogenic compounds. A mathematical model for the dynamics of androgenic synthesis, transport, metabolism, and regulation of the adult rodent ventral prostate was developed on the basis of a model by Barton and Anderson (1997). The model describes the systemic and local kinetics of testosterone (T), 5alpha-dihydrotestosterone (DHT), and luteinizing hormone (LH), with metabolism of T to DHT by 5alpha-reductase in liver and prostate. Also included are feedback loops for the positive regulation of T synthesis by LH and negative regulation of LH by T and DHT. The model simulates maintenance of the prostate as a function of hormone concentrations and androgen receptor (AR)-mediated signal transduction. The regulatory processes involved in prostate size and function include cell proliferation, apoptosis, fluid production, and 5alpha-reductase activity. Each process is controlled through the occupancy of a representative gene by androgen-AR dimers. The model simulates prostate dynamics for intact, castrated, and intravenous T-injected rats. After calibration, the model accurately captures the castration-induced regression of the prostate compared with experimental data that show that the prostate regresses to approximately 17 and 5% of its intact weight at 14 and 30 days postcastration, respectively. The model also accurately predicts serum T and AR levels following castration compared with data. This model provides a framework for quantifying the kinetics and effects of environmental and pharmaceutical endocrine active compounds on the prostate.

Agricultural chemical safety assessment: A multisector approach to the modernization of human safety requirements

Better understanding of toxicological mechanisms, enhanced testing capabilities, and demands for more sophisticated data for safety and health risk assessment have generated international interest in improving the current testing paradigm for agricultural chemicals. To address this need, the ILSI Health and Environmental Sciences Institute convened a large and diverse group of international experts to develop a credible and viable testing approach that includes scientifically appropriate studies that are necessary without being redundant, and that emphasize toxicological endpoints and exposure durations that are relevant for risk assessment. Benefits of the proposed approach include improved data for risk assessment, greater efficiency, use of fewer animals, and better use of resources. From the outset of this endeavor, it was unanimously agreed that a tiered approach should be designed to incorporate existing knowledge on the chemistry, toxicology, and actual human exposure scenarios of the compound, with integration of studies on metabolism/kinetics, life stages, and systemic toxicities. Three international task forces were charged with designing study types and endpoints on metabolism/kinetics, life stages, and systemic toxicities to be used in the tiered approach. This tiered testing proposal departs from the current standardized list of hazard studies used by many national authorities, and represents the first comprehensive effort of its kind to scientifically redesign the testing framework for agricultural chemicals.

The acquisition and application of absorption, distribution, metabolism, and excretion (ADME) data in agricultural chemical safety assessments

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. Incorporation of pharmacokinetic studies describing absorption, distribution, metabolism, and excretion is an essential tool for improving the design and interpretation of toxicity studies and their application for safety assessment. A tiered approach is described in which basic pharmacokinetic studies, similar to those for pharmaceuticals, are conducted for regulatory submission. Subsequent tiers provide additional information in an iterative manner, depending on pharmacokinetic properties, toxicity study results, and the intended uses of the compound.