An asymptotic description of a basic FcRn-regulated clearance mechanism and its implications for PBPK modelling of large antibodies.

A basic FcRn-regulated clearance mechanism is investigated using the method of matched asymptotic expansions. The broader aim of the work is to obtain further insight on the mechanism, thereby providing theoretical support for future pharmacologically-based pharmacokinetic modelling efforts. The corresponding governing equations are first non-dimensionalised and the order of magnitudes of the model parameters are assessed based on their values reported in the literature. Under the assumption of high FcRn-binding affinity, analytical approximations are derived that are valid over the characteristic phases of the problem. Additionally, relatively simple equations relating clearance and AUC to physiological model parameters are derived, which are valid over the longest characteristic time scale of the problem. For lower to moderate doses clearance is effectively linear, whereas for higher doses it is nonlinear. It is shown that for all doses sufficiently high the leading-order approximation for the IgG concentration in plasma, over the longest characteristic time scale, is independent of the initial dose. This is because IgG that is in 'excess' of FcRn is eliminated over a time scale much shorter than that of the terminal phase. In conclusion, analytical approximations of the basic FcRn mechanism have been derived using matched asymptotic expansions, leading to a simple equation relating clearance to FcRn binding affinity, the ratio of degradation and FcRn concentration, and the volumes of the system.

Fingerprinting of neurotoxic compounds using a mouse embryonic stem cell dual luminescence reporter assay.

Identification of neurotoxic drugs and environmental chemicals is an important challenge. However, only few tools to address this topic are available. The aim of this study was to develop a neurotoxicity/developmental neurotoxicity (DNT) test system, using the pluripotent mouse embryonic stem cell line CGR8 (ESCs). The test system uses ESCs at two differentiation stages: undifferentiated ESCs and ESC-derived neurons. Under each condition, concentration-response curves were obtained for three parameters: activity of the tubulin alpha 1 promoter (typically activated in early neurons), activity of the elongation factor 1 alpha promoter (active in all cells), and total DNA content (proportional to the number of surviving cells). We tested 37 compounds from the ESNATS test battery, which includes polypeptide hormones, environmental pollutants (including methylmercury), and clinically used drugs (including valproic acid and tyrosine kinase inhibitors). Different classes of compounds showed distinct concentration-response profiles. Plotting of the lowest observed adverse effect concentrations (LOAEL) of the neuronal promoter activity against the general promoter activity or against cytotoxicity, allowed the differentiation between neurotoxic/DNT substances and non-neurotoxic controls. Reporter activity responses in neurons were more susceptible to neurotoxic compounds than the reporter activities in ESCs from which they were derived. To relate the effective/toxic concentrations found in our study to relevant in vivo concentrations, we used a reverse pharmacokinetic modeling approach for three exemplary compounds (teriflunomide, geldanamycin, abiraterone). The dual luminescence reporter assay described in this study allows high-throughput, and should be particularly useful for the prioritization of the neurotoxic potential of a large number of compounds.

Prediction of in vivo developmental toxicity of all-trans-retinoic acid based on in vitro toxicity data and in silico physiologically based kinetic modeling.

The use of laboratory animals for toxicity testing in chemical safety assessment meets increasing ethical, economic and legislative constraints. The development, validation and application of reliable alternatives for in vivo toxicity testing are therefore urgently needed. In order to use toxicity data obtained from in vitro assays for risk assessment, in vitro concentration-response data need to be translated into in vivo dose-response data that are needed to obtain points of departure for risk assessment, like a benchmark dose (BMD). In the present study, we translated in vitro concentration-response data of the retinoid all-trans-retinoic acid (ATRA), obtained in the differentiation assay of the embryonic stem cell test, into in vivo dose-response data using a physiologically based kinetic model for rat and human that is mainly based on kinetic model parameter values derived using in vitro techniques. The predicted in vivo dose-response data were used for BMD modeling, and the obtained BMDL10 values [lower limit of the 95 % confidence interval on the BMD at which a benchmark response equivalent to a 10 % effect size (BMR10) is reached (BMD10)] for rat were compared with BMDL10 values derived from in vivo developmental toxicity data in rats reported in the literature. The results show that the BMDL10 values from predicted dose-response data differ about sixfold from the BMDL10 values obtained from in vivo data, pointing at the feasibility of using a combined in vitro-in silico approach for defining a point of departure for toxicological risk assessment.

Human embryonic stem cell-derived test systems for developmental neurotoxicity: a transcriptomics approach.

Developmental neurotoxicity (DNT) and many forms of reproductive toxicity (RT) often manifest themselves in functional deficits that are not necessarily based on cell death, but rather on minor changes relating to cell differentiation or communication. The fields of DNT/RT would greatly benefit from in vitro tests that allow the identification of toxicant-induced changes of the cellular proteostasis, or of its underlying transcriptome network. Therefore, the 'human embryonic stem cell (hESC)-derived novel alternative test systems (ESNATS)' European commission research project established RT tests based on defined differentiation protocols of hESC and their progeny. Valproic acid (VPA) and methylmercury (MeHg) were used as positive control compounds to address the following fundamental questions: (1) Does transcriptome analysis allow discrimination of the two compounds? (2) How does analysis of enriched transcription factor binding sites (TFBS) and of individual probe sets (PS) distinguish between test systems? (3) Can batch effects be controlled? (4) How many DNA microarrays are needed? (5) Is the highest non-cytotoxic concentration optimal and relevant for the study of transcriptome changes? VPA triggered vast transcriptional changes, whereas MeHg altered fewer transcripts. To attenuate batch effects, analysis has been focused on the 500 PS with highest variability. The test systems differed significantly in their responses (<20 % overlap). Moreover, within one test system, little overlap between the PS changed by the two compounds has been observed. However, using TFBS enrichment, a relatively large 'common response' to VPA and MeHg could be distinguished from 'compound-specific' responses. In conclusion, the ESNATS assay battery allows classification of human DNT/RT toxicants on the basis of their transcriptome profiles.

OSWG Recommended Approaches to the Nonclinical Pharmacokinetic (ADME) Characterization of Therapeutic Oligonucleotides.

This white paper summarizes the recommendations of the absorption, distribution, metabolism, and excretion (ADME) Subcommittee of the Oligonucleotide Safety Working Group for the characterization of absorption, distribution, metabolism, and excretion of oligonucleotide (ON) therapeutics in nonclinical studies. In general, the recommended approach is similar to that for small molecule drugs. However, some differences in timing and/or scope may be warranted due to the greater consistency of results across ON classes as compared with the diversity among small molecule classes. For some types of studies, a platform-based approach may be appropriate; once sufficient data are available for the platform, presentation of these data should be sufficient to support development of additional ONs of the same platform. These recommendations can serve as a starting point for nonclinical study design and foundation for discussions with regulatory agencies.

Preclinical anti-tumour activity of HexaBody-CD38, a next-generation CD38 antibody with superior complement-dependent cytotoxic activity.

BACKGROUND: HexaBody®-CD38 (GEN3014) is a hexamerization-enhanced human IgG1 that binds CD38 with high affinity. The E430G mutation in its Fc domain facilitates the natural process of antibody hexamer formation upon binding to the cell surface, resulting in increased binding of C1q and potentiated complement-dependent cytotoxicity (CDC). METHODS: Co-crystallization studies were performed to identify the binding interface of HexaBody-CD38 and CD38. HexaBody-CD38-induced CDC, antibody-dependent cellular cytotoxicity (ADCC), antibody-dependent cellular phagocytosis (ADCP), trogocytosis, and apoptosis were assessed using flow cytometry assays using tumour cell lines, and MM patient samples (CDC). CD38 enzymatic activity was measured using fluorescence spectroscopy. Anti-tumour activity of HexaBody-CD38 was assessed in patient-derived xenograft mouse models in vivo. FINDINGS: HexaBody-CD38 binds a unique epitope on CD38 and induced potent CDC in multiple myeloma (MM), acute myeloid leukaemia (AML), and B-cell non-Hodgkin lymphoma (B-NHL) cells. Anti-tumour activity was confirmed in patient-derived xenograft models in vivo. Sensitivity to HexaBody-CD38 correlated with CD38 expression level and was inversely correlated with expression of complement regulatory proteins. Compared to daratumumab, HexaBody-CD38 showed enhanced CDC in cell lines with lower levels of CD38 expression, without increasing lysis of healthy leukocytes. More effective CDC was also confirmed in primary MM cells. Furthermore, HexaBody-CD38 efficiently induced ADCC, ADCP, trogocytosis, and apoptosis after Fc-crosslinking. Moreover, HexaBody-CD38 strongly inhibited CD38 cyclase activity, which is hypothesized to relieve immune suppression in the tumour microenvironment. INTERPRETATION: Based on these preclinical studies, a clinical trial was initiated to assess the clinical safety of HexaBody-CD38 in patients with MM. FUNDING: Genmab.

Next Generation Exon 51 Skipping Antisense Oligonucleotides for Duchenne Muscular Dystrophy.

In the last two decades, antisense oligonucleotides (AONs) that induce corrective exon skipping have matured as promising therapies aimed at tackling the dystrophin deficiency that underlies the severe and progressive muscle fiber degeneration in Duchenne muscular dystrophy (DMD) patients. Pioneering first generation exon 51 skipping AONs like drisapersen and eteplirsen have more recently been followed up by AONs for exons 53 and 45, with, to date, a total of four exon skipping AON drugs having reached (conditional) regulatory US Food and Drug Administration (FDA) approval for DMD. Nonetheless, considering the limited efficacy of these drugs, there is room for improvement. The aim of this study was to develop more efficient [2'-O-methyl-modified phosphorothioate (2'OMePS) RNA] AONs for DMD exon 51 skipping by implementing precision chemistry as well as identifying a more potent target binding site. More than a hundred AONs were screened in muscle cell cultures, followed by a selective comparison in the hDMD and hDMDdel52/mdx mouse models. Incorporation of 5-methylcytosine and position-specific locked nucleic acids in AONs targeting the drisapersen/eteplirsen binding site resulted in 15-fold higher exon 51 skipping levels compared to drisapersen in hDMDdel52/mdx mice. However, with similarly modified AONs targeting an alternative site in exon 51, 65-fold higher skipping levels were obtained, restoring dystrophin up to 30% of healthy control. Targeting both sites in exon 51 with a single AON further increased exon skipping (100-fold over drisapersen) and dystrophin (up to 40%) levels. These dystrophin levels allowed for normalization of creatine kinase (CK) and lactate dehydrogenase (LDH) levels, and improved motor function in hDMDdel52/mdx mice. As no major safety observation was obtained, the improved therapeutic index of these next generation AONs is encouraging for further (pre)clinical development.

An improved model to predict physiologically based model parameters and their inter-individual variability from anthropometry.

We developed a population physiology model, physB, which provides a statistical description of the physiological characteristics in the human population, in terms of the physiological parameters that are needed in physiologically based pharmacokinetic modelling. The model predicts individual organ weights, blood flows and some respiratory parameters from anthropometric properties (body height and weight, age and gender). It draws on two existing models, PK-Pop and P(3)M, but various changes and improvements were made. The conceptual differences among the three models are discussed and they are quantitatively compared by running all three models for various specific combinations of anthropometric properties.

Toxicodynamic analysis of the combined cholinesterase inhibition by paraoxon and methamidophos in human whole blood.

Theoretical work has shown that the isobole method is not generally valid as a method for testing the absence or presence of interaction (in the biochemical sense) between chemicals. The present study illustrates how interaction can be tested by fitting a toxicodynamic model to the results of a mixture experiment. The inhibition of cholinesterases (ChE) in human whole blood by various dose combinations of paraoxon and methamidophos was measured in vitro. A toxicodynamic model describing the processes related to both OPs in inhibiting AChE activity was developed, and fit to the observed activities. This model, not containing any interaction between the two OPs, described the results from the mixture experiment well, and it was concluded that the OPs did not interact in the whole blood samples. While this approach of toxicodynamic modeling is the most appropriate method for predicting combined effects, it is not rapidly applicable. Therefore, we illustrate how toxicodynamic modeling can be used to explore under which conditions dose addition would give an acceptable approximation of the combined effects from various chemicals. In the specific case of paraoxon and methamidophos in whole blood samples, it was found that dose addition gave a reasonably accurate prediction of the combined effects, despite considerable difference in some of their rate constants, and mildly non-parallel dose-response curves. Other possibilities of validating dose-addition using toxicodynamic modeling are briefly discussed.

Toxicodynamic analysis of the inhibition of isolated human acetylcholinesterase by combinations of methamidophos and methomyl in vitro.

The applicability of dose addition to combinations of OP-esters and carbamates has been questioned based on theoretical considerations, but these have not been well supported by experimental findings. In the present study, the inhibition of AChE by combinations of methamidophos (an OP-ester) and methomyl (a carbamate) was examined in vitro. AChE inhibition was measured by the Ellman assay. We addressed the question of interaction between the OP-ester and carbamate by a toxicodynamic (TD) model reflecting the mechanism of action of the individual chemicals, without incorporating any interactions between them. The model was extended by including the experimental actions in the Ellman assay to correct for the difference in reactivation rates between phosphorylated and carbamylated AChE, which caused a bias in the observations from the assay. This zero-interactive TD model described the observations well, indicating that the OP-ester and carbamate did not interact. The applicability of dose addition was further explored by applying dose addition to the predicted inhibition by the TD model. Despite the differences in dynamics between methamidophos and methomyl, their dose-response curves were close to parallel, and dose addition gave a reasonably accurate prediction of the combined effects.

Dose addition and the isobole method as approaches for predicting the cumulative effect of non-interacting chemicals: a critical evaluation.

The prediction of the effect of cumulative exposure to similarly acting chemicals is commonly done by dose addition, such as in the relative potency factor approach. This can only be done under the assumption of zero interaction between the chemicals. The related, but not equivalent, isobole method is the most common criterion to judge whether interactions between similarly acting chemicals have taken place in a mixture experiment. Many who apply this latter method assume that it is applicable to any combination of substances, regardless of the shape of the dose-response curves of the individual substances or their underlying mechanism of action. Proponents commonly refer to the work of Berenbaum, who claimed to have proven the general applicability of the isobole method based on zero interaction. In this article, we argue that his argumentation is not generally valid. We further demonstrate that the isobole method, just like dose addition, has limited applicability. Using a physiologically based mathematical model, we provide a theoretical example of a combination of chemicals with zero interaction where the isobole method would result in the decision that they do interact. We discuss the implications for research focusing on detecting or defining interactions, and for the prediction of effects from combined exposures assuming zero interaction.

An integrated probabilistic framework for cumulative risk assessment of common mechanism chemicals in food: an example with organophosphorus pesticides.

This paper presents a framework for integrated probabilistic risk assessment of chemicals in the diet which accounts for the possibility of cumulative exposure to chemicals with a common mechanism of action. Variability between individuals in the population with respect to food consumption, concentrations of chemicals in the consumed foods, food processing habits and sensitivity towards the chemicals is addressed by Monte Carlo simulations. A large number of individuals are simulated, for which the individual exposure (iEXP), the individual critical effect dose (iCED) and the ratio between these values (the individual margin of exposure, iMoE) are calculated by drawing random values for all variable parameters from databases or specified distributions. This results in a population distribution of the iMoE, and the fraction of this distribution below 1 indicates the fraction of the population that may be at risk. Uncertainty in the assessment is treated as a separate dimension by repeating the Monte Carlo simulations many times, each time drawing random values for all uncertain parameters. In this framework, the cumulative exposure to common mechanism chemicals is addressed by incorporation of the relative potency factor (RPF) approach. The framework is demonstrated by the cumulative risk assessment of organophosphorus pesticides (OPs). By going through this example, the various choices and assumptions underlying the cumulative risk assessment are made explicit. The problems faced and the solutions chosen may be more generic than the present example with OPs. This demonstration may help to familiarize risk assessors and risk managers with the somewhat more complex output of probabilistic risk assessment.

Nonclinical Exon Skipping Studies with 2'-O-Methyl Phosphorothioate Antisense Oligonucleotides in mdx and mdx-utrn-/- Mice Inspired by Clinical Trial Results.

Duchenne muscular dystrophy is a severe, progressive muscle-wasting disease that is caused by mutations that abolish the production of functional dystrophin protein. The exon skipping approach aims to restore the disrupted dystrophin reading frame, to allow the production of partially functional dystrophins, such as found in the less severe Becker muscular dystrophy. Exon skipping is achieved by antisense oligonucleotides (AONs). Several chemical modifications have been tested in nonclinical and clinical trials. The morpholino phosphorodiamidate oligomer eteplirsen has been approved by the Food and Drug Administration, whereas clinical development with the 2'-O-methyl phosphorothioate (2OMePS) AON drisapersen was recently stopped. In this study, we aimed to study various aspects of 2OMePS AONs in nonclinical animal studies. We show that while efficiency of exon skipping restoration is comparable in young and older C57BL/10ScSn-Dmdmdx/J (mdx/BL10) mice, functional improvement was only observed for younger treated mice. Muscle quality did not affect exon skipping efficiency as exon skip and dystrophin levels were similar between mdx/BL10 and more severely affected, age-matched D2-mdx mice. We further report that treadmill running increases AON uptake and dystrophin levels in mdx/BL10 mice. Finally, we show that even low levels of exon skipping and dystrophin restoration are sufficient to significantly increase the survival of mdx-utrn-/- mice from 70 to 97 days.

The Pharmacokinetics of 2'-O-Methyl Phosphorothioate Antisense Oligonucleotides: Experiences from Developing Exon Skipping Therapies for Duchenne Muscular Dystrophy.

Delivery to the target site and adversities related to off-target exposure have made the road to clinical success and approval of antisense oligonucleotide (AON) therapies challenging. Various classes of AONs have distinct chemical features and pharmacological properties. Understanding the similarities and differences in pharmacokinetics (PKs) among AON classes is important to make future development more efficient and may facilitate regulatory guidance of AON development programs. For the class of 2'-O-methyl phosphorothioate (2OMe PS) RNA AONs, most nonclinical and clinical PK data available today are derived from development of exon skipping therapies for Duchenne muscular dystrophy (DMD). While some publications have featured PK aspects of these AONs, no comprehensive overview is available to date. This article presents a detailed review of absorption, distribution, metabolism, and excretion of 2OMe PS AONs, compiled from publicly available data and previously unpublished internal data on drisapersen and related exon skipping candidates in preclinical species and DMD patients. Considerations regarding drug-drug interactions, toxicokinetics, and pharmacodynamics are also discussed. From the data presented, the picture emerges of consistent PK properties within the 2OMe PS class, predictable behavior across species, and a considerable overlap with other single-stranded PS AONs. A level of detail on muscle as a target tissue is provided, which was not previously available. Furthermore, muscle biopsy samples taken in DMD clinical trials allowed confirmation of the applicability of interspecies scaling approaches commonly applied in the absence of clinical target tissue data.

Suppression of Mutant Protein Expression in SCA3 and SCA1 Mice Using a CAG Repeat-Targeting Antisense Oligonucleotide.

Spinocerebellar ataxia type 3 (SCA3) and type 1 (SCA1) are dominantly inherited neurodegenerative disorders that are currently incurable. Both diseases are caused by a CAG-repeat expansion in exon 10 of the Ataxin-3 and exon 8 of the Ataxin-1 gene, respectively, encoding an elongated polyglutamine tract that confers toxic properties to the resulting proteins. We have previously shown lowering of the pathogenic polyglutamine protein in Huntington's disease mouse models using (CUG)7, a CAG repeat-targeting antisense oligonucleotide. Here we evaluated the therapeutic capacity of (CUG)7 for SCA3 and SCA1, in vitro in patient-derived cell lines and in vivo in representative mouse models. Repeated intracerebroventricular (CUG)7 administration resulted in a significant reduction of mutant Ataxin-3 and Ataxin-1 proteins throughout the brain of SCA3 and SCA1 mouse models, respectively. Furthermore, in both a SCA3 patient cell line and the MJD84.2 mouse model, (CUG)7 induced formation of a truncated Ataxin-3 protein species lacking the polyglutamine stretch, likely arising from (CUG)7-mediated exon 10 skipping. In contrast, skipping of exon 8 of Ataxin-1 did not significantly contribute to the Ataxin-1 protein reduction observed in (CUG)7-treated SCA1154Q/2Q mice. These findings support the therapeutic potential of a single CAG repeat-targeting AON for the treatment of multiple polyglutamine disorders.

Assessment of Dermal Absorption of Aluminum from a Representative Antiperspirant Formulation Using a 26 Al Microtracer Approach.

A clinical pharmacokinetic study was performed in 12 healthy women to evaluate systemic exposure to aluminum following topical application of a representative antiperspirant formulation under real-life use conditions. A simple roll-on formulation containing an extremely rare isotope of aluminum (26 Al) chlorohydrate (ACH) was prepared to commercial specifications. A 26 Al radio-microtracer was used to distinguish dosed aluminum from natural background, using accelerated mass spectroscopy. The 26 Al citrate was administered intravenously (i.v.) to estimate fraction absorbed (Fabs ) following topical delivery. In blood samples after i.v. administration, 26 Al was readily detected (mean area under the curve (AUC) = 1,273 ± 466 hours×fg/mL). Conversely, all blood samples following topical application were below the lower limit of quantitation (LLOQ; 0.12 fg/mL), except two samples (0.13 and 0.14 fg/mL); a maximal AUC was based on LLOQs. The aluminum was above the LLOQ (61 ag/mL) in 31% of urine samples. From the urinary excretion data, a conservative estimated range for dermal Fabs of 0.002-0.06% was calculated, with a mean estimate of 0.0094%.

To Apply Microdosing or Not? Recommendations to Single Out Compounds with Non-Linear Pharmacokinetics.

Microdosing studies allow clinical investigation of pharmacokinetics earlier in drug development, before all high-dose safety concerns have been sorted out. Furthermore, microdosing allows inclusion of target groups that are inadmissible in high-dose phase I trials. A potential concern when considering a microdosing study is that a particular drug candidate may display non-linear pharmacokinetics. Saturation of, for example, membrane transport or metabolism at exposure levels between the microdose and therapeutic dose may limit the predictivity of high-dose pharmacokinetics from microdose observations. Guidance on the likelihood of appreciable non-linear pharmacokinetics based on preclinical information can be helpful in staging the clinical phase and the place of microdosing in it. We present a decision tree that evaluates concerns about non-linearities raised in the preclinical phase and their potential impact on the proportionality between microdose and intended therapeutic dose as predicted from preclinical information. The expected maximum concentrations at relevant sites are estimated by non-compartmental methods. These are compared with dissolution, Michaelis constants for active or enzymatic processes, and binding protein concentrations to assess the potential saturation of the processes below therapeutic doses. The decision tree was applied to ten published cases comparing microdose and therapeutic dose pharmacokinetics, for which concerns about non-linear pharmacokinetics were raised a priori. The decision tree was able to discriminate cases showing substantial non-linearities from cases displaying dose-proportional pharmacokinetics. The recommendations described in this paper may be useful in deciding whether a microdosing study is a sensible option to gain early insight in clinical pharmacokinetics of drug candidates.

Proceedings in uncovering the mechanism behind peroxisome proliferator-induced hepatocarcinogenesis.

The cancer risk assessment of peroxisome proliferators (PPs) has been a matter of debate for decades. PPs are non-genotoxic carcinogens that cause liver tumours in rodents. There is considerable evidence that humans are refractory to the carcinogenic effect of PPs. Still, some toxicologists argue that these chemicals should be considered carcinogenic until the opposite has been proven. To solve this matter, much work has been dedicated to uncovering the mode of action of PPs. The proceedings are reviewed and it is discussed whether these provide an explanation for the observed interspecies differences and shed a new light on human health risk assessment. The efforts in the past years have significantly increased our understanding of the pathways involved, but have also raised new questions. The question why humans seem to be refractory is yet to be answered. From this point of view, the safety assessment of PPs to humans therefore remains an issue of discussion.

Interpreting in vitro developmental toxicity test battery results: The consideration of toxicokinetics.

In the EU collaborative project ChemScreen an alternative, in vitro assay-based test strategy was developed to screen compounds for reproductive toxicity. A toxicokinetic modeling approach was used to allow quantitative comparison between effective concentrations in the in vitro test battery and observations of developmental toxicity in vivo. This modeling strategy is based on (1) the definition of relevant observations of toxicity in vivo, (2) simulation of the corresponding systemic concentrations in vivo by toxicokinetic modeling, and (3) correction for differences in protein binding and lipid partitioning between plasma and in vitro test media. The test results of a feasibility study with a number of known reproductive toxicants has been described previously (Piersma et al. [15]). In the present paper, we take a more detailed look at the toxicokinetics of these compounds, and add the analysis of some compounds from subsequent studies. We discuss how the consideration of toxicokinetics allowed comparison between test systems with differing test medium composition, has helped to interpret the in vitro findings in light of in vivo observations, and to gain confidence in the predictive value of the test battery outcomes. The same toxicokinetic modeling strategy, in reverse order, can now be used for risk assessment purposes to predict toxic doses in vivo from effective concentrations in vitro.

Toxicokinetic Triage for Environmental Chemicals.

Toxicokinetic (TK) models link administered doses to plasma, blood, and tissue concentrations. High-throughput TK (HTTK) performs in vitro to in vivo extrapolation to predict TK from rapid in vitro measurements and chemical structure-based properties. A significant toxicological application of HTTK has been "reverse dosimetry," in which bioactive concentrations from in vitro screening studies are converted into in vivo doses (mg/kg BW/day). These doses are predicted to produce steady-state plasma concentrations that are equivalent to in vitro bioactive concentrations. In this study, we evaluate the impact of the approximations and assumptions necessary for reverse dosimetry and develop methods to determine whether HTTK tools are appropriate or may lead to false conclusions for a particular chemical. Based on literature in vivo data for 87 chemicals, we identified specific properties (eg, in vitro HTTK data, physico-chemical descriptors, and predicted transporter affinities) that correlate with poor HTTK predictive ability. For 271 chemicals we developed a generic HT physiologically based TK (HTPBTK) model that predicts non-steady-state chemical concentration time-courses for a variety of exposure scenarios. We used this HTPBTK model to find that assumptions previously used for reverse dosimetry are usually appropriate, except most notably for highly bioaccumulative compounds. For the thousands of man-made chemicals in the environment that currently have no TK data, we propose a 4-element framework for chemical TK triage that can group chemicals into 7 different categories associated with varying levels of confidence in HTTK predictions. For 349 chemicals with literature HTTK data, we differentiated those chemicals for which HTTK approaches are likely to be sufficient, from those that may require additional data.