Use of human liver and EpiSkin™ S9 subcellular fractions as a screening assays to compare the in vitro hepatic and dermal metabolism of 47 cosmetics‐relevant chemicals

In vitro to in vivo extrapolation to derive a metabolism factor for estimating the aggregate exposure to salicylic acid after dermal exposure of its esters

As part of the safety assessment of salicylate esters in cosmetics, we developed a metabolism factor based on in vitro to in vivo extrapolation (IVIVE) to provide a better estimation of the aggregate internal exposure to the common metabolite, salicylic acid. Optimal incubation conditions using human liver S9 were identified before measuring salicylic acid formation from 31 substances. Four control substances, not defined as salicylic esters but which could be mistaken as such due to their nomenclature, did not form salicylic acid. For the remaining substances, higher in vitro intrinsic clearance (CLint, in vitro) values generally correlated with lower LogP values. A "High-Throughput Pharmacokinetic" (HTPK) model was used to extrapolate CLint, in vitro values to human in vivo clearance and half-lives. The latter were used to calculate the percentage of substance metabolised to salicylic acid in 24 h in vivo following human exposure to the ester, i.e. the "metabolism factor". The IVIVE model correctly reproduced the observed elimination rate of 3 substances using in silico or in vitro input parameters. For other substances, in silico only-based predictions generally resulted in lower metabolism factors than when in vitro values for plasma binding and liver S9 CLint, in vitro were used. Therefore, in vitro data input provides the more conservative metabolism factors compared to those derived using on in silico input. In conclusion, these results indicate that not all substances contribute equally (or at all) to the systemic exposure to salicylic acid. Therefore, we propose a realistic metabolism correction factor by which the potential contribution of salicylate esters to the aggregate consumer exposure to salicylic acid from cosmetic use can be estimated.

Evaluation of new alternative methods for the identification of estrogenic, androgenic and steroidogenic effects: a comparative in vitro/in silico study

A suite of in vitro assays and in silico models were evaluated to identify which best detected the endocrine-disrupting (ED) potential of 10 test chemicals according to their estrogenic, androgenic and steroidogenic (EAS) potential compared to the outcomes from ToxCast. In vitro methods included receptor-binding, CALUX transactivation, H295R steroidogenesis, aromatase activity inhibition and the Yeast oestrogen (YES) and Yeast androgen screen (YAS) assays. The impact of metabolism was also evaluated. The YES/YAS assays exhibited a high sensitivity for ER effects and, despite some challenges in predicting AR effects, is a good initial screening assay. Results from receptor-binding and CALUX assays generally correlated and were in accordance with classifications based on ToxCast assays. ER agonism and AR antagonism of benzyl butyl phthalate were abolished when CALUX assays included liver S9. In silico final calls were mostly in agreement with the in vitro assays, and predicted ER and AR effects well. The efficiency of the in silico models (reflecting applicability domains or inconclusive results) was 43-100%. The percentage of correct calls for ER (50-100%), AR (57-100%) and aromatase (33-100%) effects when compared to the final ToxCast call covered a wide range from highly reliable to less reliable models. In conclusion, Danish (Q)SAR, Opera, ADMET Lab LBD and ProToxII models demonstrated the best overall performance for ER and AR effects. These can be combined with the YES/YAS assays in an initial screen of chemicals in the early tiers of an NGRA to inform on the MoA and the design of mechanistic in vitro assays used later in the assessment. Inhibition of aromatase was best predicted by the Vega, AdmetLab and ProToxII models. Other mechanisms and exposure should be considered when making a conclusion with respect to ED effects.

Limits of the regulatory evaluation of resorcinol as a thyroid disruptor: When limited experimental data challenge established effects in humans

Severe hypothyroidism has been reported in humans during resorcinol therapeutic use. However, available data highlight differences in the severity of resorcinol-induced thyroid effects between humans and rodents, leading to a debate on the relevance of human data for its classification as a thyroid disruptor. The aim of this review is to illustrate some of the limitations of the evaluation framework for thyroid disrupters using resorcinol as a case study of a chemical with clear thyroid-disrupting properties in humans that could not have been identified solely from regulatory studies on animals. The reliability of human data has been called into question due to the specific exposure patterns in humans and the paucity of robust toxicokinetic data. In humans, therapeutic use of resorcinol induces severe hypothyroidism, but in rodents, thyroid disruption is limited to decreased thyroxine concentrations and histological changes in the thyroid. The adverse effects of thyroid disruption, such as impaired neurodevelopment, have not been sufficiently investigated, and experimental neurobehavioral data for resorcinol remain scarce and inconclusive. Although regulatory toxicological evaluations have not included in-depth investigations of thyroid regulation and related adverse effects, they have been used to challenge the relevance of human data. Resorcinol is an emblematic example of how the framework for regulatory evaluations of thyroid disruptors relies almost exclusively on animal studies which may not be suitable for assessing thyroid disruption. This review highlights the need to revise regulatory guidelines and to adopt strategies based on up-to-date, scientifically sound approaches to identify thyroid disruptors. The limits of the current regulatory framework for identifying thyroid disruptors can lead to opposing positions between regulatory bodies. The French Agency for Food, Environmental and Occupational Health & Safety (ANSES)'proposal to identify resorcinol as a "substance of very high concern" due to its ED properties has not been adopted by the European instances.

Impact of chemical structure on the in vitro hydrolysis of fatty esters of 2-ethylhexanoic acid or 2-ethylhexanol and extrapolation to the in vivo situation

Fatty esters of 2-ethylhexanoic acid (EHA) and 2-ethylhexanol (EH) are commonly used in cosmetics. Human liver and skin S9 and human plasma were used to determine the in vitro rates of clearance (CL_int) of a series of compounds, with a range of 2-11 carbons on the acid or alcohol moiety and branching at the C2 position. The impact of carbon chain length on in vitro CL_int was most prominent for the liver metabolism of esters of EH, while for in vitro skin metabolism it was greater for esters of EHA. The position of the branching also impacted the liver hydrolysis rates, especially for the C3, C4, and C5 esters with lower CL_{int in vitro} rates for esters of EHA relative to those of EH. When the in vitro intrinsic clearance rates were scaled to in vivo rates of hepatic clearance, all compounds approximated the rate for hepatic blood flow, mitigating this dependence of metabolism on structure. This work shows how structural changes to the molecule can affect in vitro metabolism and, furthermore, allows for an estimation of the in vivo metabolism.

Use of in vitro metabolism and biokinetics assays to refine predicted in vivo and in vitro internal exposure to the cosmetic ingredient, phenoxyethanol, for use in risk assessment

A novel approach was developed to help characterize the biokinetics of the cosmetic ingredient, phenoxyethanol, to help assess the safety of the parent and its major stable metabolite. In the first step of this non-animal tiered approach, primary human hepatocytes were used to confirm or refute in silico predicted metabolites, and elucidate the intrinsic clearance of phenoxyethanol. A key result was the identification of the major metabolite, phenoxyacetic acid (PAA), the exposure to which in the kidney was subsequently predicted to far exceed that of phenoxyethanol in blood or other tissues. Therefore, a novel aspect of this approach was to measure in the subsequent step the formation of PAA in the cells dosed with phenoxyethanol that were used to provide points of departure (PoDs) and express the intracellular exposure as the Cmax and AUC24. This enabled the calculation of the intracellular concentrations of parent and metabolite at the PoD in the cells used to derive this value. These concentrations can be compared with in vivo tissue levels to conclude on the safety margin. The lessons from this case study will help to inform the design of other non-animal safety assessments.

Use of Physiologically-Based Kinetics Modelling to Reliably Predict Internal Concentrations of the UV Filter, Homosalate, After Repeated Oral and Topical Application

Ethical and legal considerations have led to increased use of non-animal methods to evaluate the safety of chemicals for human use. We describe the development and qualification of a physiologically-based kinetics (PBK) model for the cosmetic UV filter ingredient, homosalate, to support its safety without the need of generating further animal data. The intravenous (IV) rat PBK model, using PK-Sim^®, was developed and validated using legacy in vivo data generated prior to the 2013 EU animal-testing ban. Input data included literature or predicted physicochemical and pharmacokinetic properties. The refined IV rat PBK model was subject to sensitivity analysis to identify homosalate-specific sensitive parameters impacting the prediction of C_max (more sensitive than AUC_(0-∞)). These were then considered, together with population modeling, to calculate the confidence interval (CI) 95% C_max and AUC_(0-∞). Final model parameters were established by visual inspection of the simulations and biological plausibility. The IV rat model was extrapolated to oral administration, and used to estimate internal exposures to doses tested in an oral repeated dose toxicity study. Next, a human PBK dermal model was developed using measured human in vitro ADME data and a module to represent the dermal route. Model performance was confirmed by comparing predicted and measured values from a US-FDA clinical trial (Identifier: NCT03582215, https://clinicaltrials.gov/). Final exposure estimations were obtained in a virtual population and considering the in vitro and input parameter uncertainty. This model was then used to estimate the C_max and AUC_(0–24 h) of homosalate according to consumer use in a sunscreen. The developed rat and human PBK models had a good biological basis and reproduced in vivo legacy rat and human clinical kinetics data. They also complied with the most recent WHO and OECD recommendations for assessing the confidence level. In conclusion, we have developed a PBK model which predicted reasonably well the internal exposure of homosalate according to different exposure scenarios with a medium to high level of confidence. In the absence of in vivo data, such human PBK models will be the heart of future completely non-animal risk assessments; therefore, valid approaches will be key in gaining their regulatory acceptance.

Clinical Trial Registration: https://clinicaltrials.gov/, identifier, NCT03582215

Demonstration of the first-pass metabolism in the skin of the hair dye, 4-amino-2-hydroxytoluene, using the Chip2 skin-liver microphysiological model

We used TissUse's HUMIMIC Chip2 microfluidic model, incorporating reconstructed skin models and liver spheroids, to investigate the impact of consumer-relevant application scenarios on the metabolic fate of the hair dye, 4-amino-2-hydroxytoluene (AHT). After a single topical or systemic application of AHT to Chip2 models, medium was analysed for parent and metabolites over 5 days. The metabolic profile of a high dose (resulting in a circuit concentration of 100 μM based on 100% bioavailability) of AHT was the same after systemic and topical application to 96-well EpiDerm™ models. Additional experiments indicated that metabolic capacity of EpiDerm™ models were saturated at this dose. At 2.5 μM, concentrations of AHT and several of its metabolites differed between application routes. Topical application resulted in a higher C_max and a 327% higher area under the curve (AUC) of N-acetyl-AHT, indicating a first-pass effect in the EpiDerm™ models. In accordance with in vivo observations, there was a concomitant decrease in the C_max and AUC of AHT-O-sulphate after topical, compared with systemic application. A similar alteration in metabolite ratios was observed using a 24-well full-thickness skin model, EpiDermFT™, indicating that a first-pass effect was also possible to detect in a more complex model. In addition, washing the EpiDermFT™ after 30 min, thus reflecting consumer use, decreased the systemic exposure to AHT and its metabolites. In conclusion, the skin-liver Chip2 model can be used to (a) recapitulate the first-pass effect of the skin and alterations in the metabolite profile of AHT observed in vivo and (b) provide consumer-relevant data regarding leave-on/rinse-off products.

Characterization of application scenario-dependent pharmacokinetics and pharmacodynamic properties of permethrin and hyperforin in a dynamic skin and liver multi-organ-chip model

Microphysiological systems (MPS) aim to mimic the dynamic microenvironment and the interaction between tissues. While MPS exist for investigating pharmaceuticals, the applicability of MPS for cosmetics ingredients is yet to be evaluated. The HUMIMIC Chip2 ("Chip2″), is the first multi-organ chip technology to incorporate skin models, allowing for the topical route to be tested. Therefore, we have used this model to analyze the impact of different exposure scenarios on the pharmacokinetics and pharmacodynamics of two topically exposed chemicals, hyperforin and permethrin. The Chip2 incorporated reconstructed human epidermis models (EpiDerm™) and HepaRG-stellate spheroids. Initial experiments using static incubations of single organoids helped determine the optimal dose. In the Chip2 studies, parent and metabolites were analyzed in the circuit over 5 days after application of single and repeated topical or systemic doses. The gene expression of relevant xenobiotic metabolizing enzymes in liver spheroids was measured to reflect toxicodynamics effects of the compounds in liver. The results show that 1) metabolic capacities of EpiDerm™ and liver spheroids were maintained over five days; 2) EpiDerm™ model barrier function remained intact; 3) repeated application of compounds resulted in higher concentrations of parent chemicals and most metabolites compared to single application; 4) compound-specific gene induction e.g. induction of CYP3A4 by hyperforin depended on the application route and frequency; 5) different routes of application influenced the systemic concentrations of both parents and metabolites in the chip over the course of the experiment; 6) there was excellent intra- and inter-lab reproducibility. For permethrin, a process similar to the excretion in a human in vivo study could be simulated which was remarkably comparable to the in vivo situation. These results support the use of the Chip2 model to provide information on parent and metabolite disposition that may be relevant to risk assessment of topically applied cosmetics ingredients.

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

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