Characterization of xenobiotic metabolizing of a reconstructed human epidermal model from adult hair follicles Daniel Bacqueville, Carine Jacques, Laure Duprat, Emilien L. Jamin, Beatrice Guiraud, Elisabeth Perdu, Sandrine Bessou-Touya, Daniel Zalko, Hélène Duplan

To cite this version:

Daniel Bacqueville, Carine Jacques, Laure Duprat, Emilien L. Jamin, Beatrice Guiraud, et al.. Characterization of xenobiotic metabolizing enzymes of a reconstructed human epidermal model from adult hair follicles. Toxicology and Applied Pharmacology, Elsevier, 2017, 329, pp.190-201. ￿10.1016/j.taap.2017.05.040￿. ￿hal-01608588￿

HAL Id: hal-01608588 https://hal.archives-ouvertes.fr/hal-01608588 Submitted on 26 May 2020

HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied journal pertain. be it manuscript a This 10.1016/j.taap.2017.05.040 captured epidermal Duplan Jamin, Please Accepted date: Revised date: Received date: T Reference: DOI: PII: T Jamin, Daniel reconstructed humanepidermalmodelfromadulthairfollicles Characterization Accepted Manuscript o appearin: ouya, DanielZalko,HélèneDuplan is service discovered published is Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040 cite Beatrice a Bacqueville, , Beatrice Characterization PDF as to model this will af our in which filiation file article under Guiraud, customers its from Guiraud, of of Comment citer cedocument: final Carine could go an xenobiotic as: adult 30 May2017 10 May2017 7 March2017 T YT doi: S0041-008X(17)30258-2 for oxicology andAppliedPharmacology copyediting, unedited form. Elisabeth Daniel AAP 13975 of 10.1016/j.taap.2017.05.040 all af Elisabeth we hair Jacques, xenobiotic fect authors.

Please are follicles. Bacqueville, the manuscript metabolizing Perdu, providing typesetting, content, Laure note Perdu, metabolizing Please Sandrine The that Duprat, that this Sandrine and Carine address during check enzymes and has early all Bessou-T review Emilien enzymes Jacques, legal been if the for Bessou- version appropriate. production the of disclaimers accepted of ouya, L. corresponding Laure of the a of a resulting Daniel reconstructed the for Duprat, process Ytaap(2017), that manuscript. publication. Zalko, proof apply author Emilien errors Hélène human before to may doi: was The the As L. Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied Manuscript correspon  UPS, Toulouse, France. 2 Département Pharmacologie, 1 affiliations:Complete Perdu DanielBacqueville adulthairfollicles Characterizationofxenobiotic metabolizing enzymes of a reconstructed humanepidermal model from E Fax: Tel.:+33 534506427 DanielBacqueville *

Toxalim (Research ToxalimCentre (Research in Food Toxicology), UniversitéToulouse,de INRA, ENVT,INP Pierre Fabre Dermo - Equalcontribution C mail address: orrespo

+33 2 , , SandrineBessou Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040 534503424 nding author:

daniel

1,*

ACCEPTED MANUSCRIPT , -  . cosmétique, , bacqueville

Carine Jacques

dence ACCEPTED MANUSCRIPT - Comment citer cedocument: Touya Centre R&D Pierre Fabre, 3avenue Hubert Curien,Toulouse, France.

1 ,

@pierre DanielZalko ServicePharmacologie Division 2

1 ,  , - Laure Duprat fabre.com 2 , Hélène Duplan

1 , Emilien L. Jamin 1

et Pharmacocinétique et Cutané, 2 , Beatrice Guiraud 1 - , Purpan, Elisabeth

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied andnon metabolized modelswith7 including NAT1. generallywere higher than those of Phase after topicalandsystemic application m expressionlevel nativehuman skin alsowere present at thegene levelintheORS was consistentwith native human skin. fromtheouter root sheath (ORS) ofhair follicles ( gene expressionandenzyme functionality Inthis study, Abstract functional activity words: Key cosmeticingredients. ofbe interestwhen evaluating application routes,supporting their use incutaneousmetabolism studies.Such a reproducible modelwill induced Inconclusion, the ORS odel was - by saturable

Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040

responsiveto AhRligands

within 24 h after topicalapplication. reconstructed a a comprehensivecharacterization of - ethoxycoumarin, testosterone, benzo(a)py The presence of functional CYPs, UGTsand SULTs was confirmedbyincubating the

ACCEPTED metabolism byMANUSCRIPTthe XMEsandthe on residence timewithin themodel.

s

were - Functional XME . activities also were demonstrated using systemicor topical RHE modelexpress ACCEPTED MANUSCRIPT the Comment citer cedocument: human

observed for alcohol/aldehyde dehydrogenasesand (carboxyl) CYP inducers,CYP the cutaneous metabolism and potential toxicity ofinnovative dermo

epiderm s ,

evident at thegene andenzyme activity level is, I was was made

XMEs,thehighest levelswere 3 es - hairfollicle, cutaneous metabolism, gene expression, methylcholanthrene (

a a numberof

The extentof metabolism was dependent on ORS xenobiotic metabolizing enzymes ( in P450( - a a reconstructed skinepidermal model derived RHE rene and rene 3 Phase I Phase CYPs ) . The - RHE model. )

- 3 and consistentlyreported todetected be in ORS MC - MC II , ofallwhich were -

) and RHE model XMEs, GSTs Thehighest  -

and , naphthoflavone someof .

XMEs Phase II Phase XME

whichmay be esterases. The esterases. I Phase rapidly ) gene based onbased

saturable XME

profile ( β

XME - NF)

levels

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied metabolizing enzymes(XMEs). carboxyl ubiquitin reversetranscription (RT PCR epiderm human reconstructed Economic for Organization (3 methylcholanthrene 10 S ethoxycoumarin(7 (CYP); ( (AhR); Abbreviations: β - (GST);transferaseHighSp Resolution Mass (CES1/2); 2 and 1 carboxylesteraseNF); ( HSD17B10 pxd hdoae EH) etrs D ED; S suc i te eaie oe (NI mode negative the in source ESI (ESD); D esterase (EPHX); Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040 APCI source in the positiveinsourceAPCIthe mode(PI loo dhdoeae (D) adhd dhdoeae (ALDH); dehydrogenases aldehyde (ADH); dehydrogenases alcohol ; yrxcuai (C; long (HC); hydroxycoumarin ); - - EC); ethoxyresorufinO EC); ACCEPTED MANUSCRIPT xenobiotic (UGT); glucuronosyltransferase UDP (UCHL3); L3 esterase terminal - C; N(alpha) MC); ACCEPTED MANUSCRIPT Comment citer cedocument:

- is q

oprto ad eeomn (ED; ue ro sheath root outer (OECD); Development and Cooperation

PCR); PCR); spermine (ORS

- H) prostaglandinRHE); - ctlrnfrs 2 (A2) N (NAA20); 20 acetyltransferase - deethylase (EROD); flavinglutathionedeethylase(EROD); (FMO);monooxidase catechol ectrometry (HRMS);ectrometry

N1 - acetyltransferase1 (SAT1); (SULT); - - - chain APCI);benzo(a)pyrene (B(a)P); O - methyltransferase - - fatty noeoie ytae PG2; el time real (PTGS2); synthase endoperoxide - acid hydroxysteroid(17 - o lgs 1 AS1; 3 (ACSL1); 1 CoA ( - COMT ctl rnfrs (NAT); transferase acetyl ); - P450  hydrocarbon aryl  -

naphthoflavone dehydrogenase - - 7 ESI); derived - -

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied modeltometabolism studies genoprotection bipyrimidine photoproduct production modelbeenhas shown to a be gooda differentiation markers similarwere tothose of native epidermis (Guiraudal., et 2014). alternativetonative human epidermis, suchthat thearchitecture, generalstratification andlocalization commonly usedin commercial models. healthyordiseased) ofease use of epidermalmodel follicles. Thismodel humanskinequivalent model The 1 the in interprettoxicity endpoints,important itis to characterize the metabolic capacityof theskin model metabolites(Bronaughand Stewart, 1985 metabolism, which can detoxify them to more water (Lockley et al., 2 The rate ofpercutaneous absorption ofsome compounds isknown tobe influenced byskin metabolism in these models (Eilsteinal., et 2015; characterizing Phenion

bio ORS outersheath root assays. assays. ORS the rapidandnon Introduction ® -

RHE model Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040 Full

derivedkeratinocytes affords a numberofadvantages over otherRHE models Inthis study, we aimedto characterize as precisely as possible the metabolic capabilities of -

hcns skinmodelThickness andEpiDerm™

the function (Bacquevilleal., et 2015) 004; Jacques Jacques et 004;al.,2010 that

ACCEPTED MANUSCRIPT and

was developedin orderto establish a cost - could be customized according to specific assay needs Tothis end,we derivedreco ,

- ACCEPTED MANUSCRIPT moreover, invasive Comment citer cedocument: al

activitiesofseveral I Phase in engineered the same thesame way as nstructed humanepidermal samplecol the lternativeto human skinfor Jäckh al., et 2011 , firstinvestigated

repair The select . , Oesch , a fromadult undifferentiated keratinocytes from Here, Here, o ), anda ), ORS

lection, theability toselect a particular donor ionof andinduction ofapoptosis

previousreports

- urintention et et al. RHE model s s compoundspenetrate the skin, theymay be subject to

adult donors, as opposed to neonatal donors that are - , solublemetabolites orbioactivate them to reactive

that , 2007, ) the .

and

support overall

and II was was ( was was shown be to - ORS

effective andreproducible xenobiotic metabolizing enzymes (XMEs)

us 2 studying of theeffects

basal to extend 014 their ing - RHE) model ). Therefore, ). inorderto help RHE models applicability tothis purpose XMEgene expression profile andinvestigating ) the theuse ofthe ORS

(Guiraud et al.,2014)

a a good

is , a a three

such as suchas andrepro This UV UV radiation(e.g. ORS ,

in in vitro including the - (gender,age, dimensional ORS SkinEthic,

cells - sunscreen RHE ducible

- . human s

RHE ofhair

The

used by of of

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied usings pathway complexities functional XMEactivities by measuring the Inaddition toanalyzing theXME profile ofthe ORS CYP1A1/2 andCYP1B1 selective ,ethoxyresorufin. addition, expression located and3 2013). suchmodels couldbe altered during thecourse of usedto investigate topically applied chemicals. t the absorption as one ofthemodel compounds to develop skin alternative models. corresponding technicalguidance docume EconomicCooperation andDevelopment (OECD) adopted guideline 28 functional reportedto be good metabolized through reductivepathways. Sincethereductive capacity oftheskin has previously been steroidand well a metabolize steroids both testosterone, studies,tocapture thekinetics ofmetabolism. he mosthe relevant route of ORS P - methylchol h imilar Therefore, we measured the ofeffect two

ase I ase belowthe model) andtopical (added the tou

- we we confirmedinduction effect RHE Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040 oxido of

and B(a)P

timepoint a model XMEgene does not necessarily indicate thatthe - eutv capacity reductiveofthemodel a II - n established diagnostic substrate for activities.CYP

, and

hee ( threne sequentialpathways ( ( ACCEPTED MANUSCRIPT Haag Haag et al.,2012 vanEijl . using s The twoendpoints

3 (24 ACCEPTED MANUSCRIPT - Comment citer cedocument: exposureto theskin istopical MC 3 RT

h) - et et al MC . . 7 . . Additional timepoints - qPCR - ., 2012., ) EC represents represents EC a chemicalwith a ,XMEon gene )

; therefore, ; we selected t analysis

s (De (De Kanter et al.,1999) )

by incubating controlandinducer nt testosterone represents an idea (gene expressionand

metabolism ( OECD,2004b

The test chemicals 7 were of Thisisalso of importance sincethe metabolic capacity of

. knownAhR

morethan Testosterone isalso specified inthe Organization for expression topical pper surface ppersurface of themodel) -

RHE model, we compl (48 and(48 h) added72were for functional metabolic itimportantis to show, thatRHE models can be offour test chemicalswi treatment with a test compound (Hewitt et )

9

todescribe methodsassessing drug after after systemic (addedto theculture medium ligands, namely 0

estosterone Phase I Phase . functional corresponding Theskin knownis to ablebe to

simple metaboliteprofile Thissubstrate can alsobe

and 7 - lsubstrate toconfirm the ( ethoxycoumarin (7 OECD, 2004a - metabolism) as treated treated e 0

mented thismentedby a 

Phase Phase B(a)P XMEfunctionalis application - model naphthoflavone ( th varyingth skin

II is ofan

) XME a a knowngenotoxi comparedwere modelswith the and a .

T endogenous - involving metabolic EC), investigating he  In . NF .

al., Since

)

c ,

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied evaluated, as well as theapplication of test compoundsby thesystemic and topical routes. and RHE modelanalyzingby bothgeneexpression profiles These studiesaimed toprovide a comprehensive measure itsmetabolism was Itused toinvestigat models(Aardema et al.,2013) humanskin andUGTs, soit metabolism model, Epiderm metabolite compound hydroxybenzo(a)pyrene, 9 4,5 Benzo(a)pyrene fromSigma MoravekBiochemicals Amersham Biosciences(Buckinghamshire, UK). testosterone with 7 2.1 2 -

ethoxy[3 -

dihydrodiol, benzo(a)pyrene test test chemicalmetabolism) Chemicals andMaterials methods Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040 s - which , especially, tohelp interpret negative outcomes in genotoxicity assays usin - 14 ( andthat th lrc (StQuentinAldrich Fallavier, France) had a purity greater than 98.5%. Conney et al., 1994 C] - ™ - isimportant tounderstand thebalance ofmetabolism ofsuch genotoxi a coumarin, [U

r specificactivities of2,2.074, - , toxicity 7 24h

,

t ACCEPTED MANUSCRIPT - 8 is

, e theinducibility e ofXMEs in ORS t (Brea, United States). alongside after - ACCEPTED MANUSCRIPT model 9 Comment citer cedocument: relieson - , hydroxybenzo(a)pyrene, 3 t - 10 topicalapplication - . .

14 - Thesuitability ofthemodel forcompound metabolism assays was 3 tetrahydrotetrol respon ) - C] - . 1,6

MC its inductionits effects. B(a)P - CYP450 a - bisphenolA,[5,6 dione

is ds

a a directgenotoxin

has been to CYP1 induction ,

benzo(a)pyrene - Unlabeled7 depend

, benz,

2 (Götzal., et 2012), suggesting itinduces its own

characterizationof theXMEs present inthe and 3 also - M - o(a)pyrene - e 2.03GBq/mmol, respectively, purchasedwere from 14 ethylcholanthrene hydroxybenzo(a)pyrene shownto induce EROD activities intheepidermal nt bioactivation, into

C] b C] (PCR)

- - EC and - . RHEs RHEs and therefore it enzo(a)pyrene[ 3,6 B(a)P ,

and functional activities - is metabolizedis B(a)P in, benzo(a)pyrenedione, - cis

- is detoxifiedalso byGSTs, SULTs 7,8 , 3 - - - dihydrodiol [ MC 3 H] was was purchased H] from 14

and testosterone,purchased

C 7 , ultimate DNA reactive via - B(a)P], B(a)P], [4,7 -

hydroxybenzo(a)pyrene was P benz , hase hase I

gskin epidermal

- ofinterest to 7,8 ( c compound c EROD EROD activity

and o(a)pyrene -

in, 8 dione, - 14 novel ORS II C]

pathways - -

- in cis - - .

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied Yvelines, Merck (Briare andethylfrom acetate Scharlau chemie S.A. (Barcelona,Spain); ethanoland acetic acidfrom acetate, sodium acetate, phosphate buffer(0.01 Other chemicalsand solvents (analytical grade) purchasedwere from thefollowing sources: ammonium dione. an 11 hydroxytestosterone,11 hydroxytestosterone,6 Testosterone metabolitesreference fromwere Sigma (Kansas City,MO,USA) and obtainedwere thefrom National Cancer Institut Briefly, cultureof the houseaccording three Humanepidermis was reconstructed 2.2 purchased fromSigma (streptomycin/penicillin, gentamycin)andtheantifongic (fungizone (DMEM) andL drostan 

- ol healthy

- 3,17 Preparation humanofreconstituted epidermis pluckedhairfollicles explanted were - Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040 France) France) was used for 3 -  in, 4 dione, - Caucasian ORS - ol Le - -

g 17 to lutamine obtained were from Gibco - Canal, France). Ultrapure water from a Milli

- - the RHE tissues andthestructure of thetissues at thetime of incubation (histology). ACCEPTED MANUSCRIPT 5 one, - androsten - 

Aldrich. protocol - - volunteerswhoha ketotestosterone,1 hydroxytestosterone, 7 ACCEPTED MANUSCRIPT Comment citer cedocument:  - androstan were were

- HPLC HPLC mobilephase 16 describedby Guiraud et al.

 used compounds.as reference - ol using - -

3 3,17  - ol d given theirinformed consent. T ORS - - - dione dehydrotestosterone,4 17

on a a on

-  Chemical e CarcinogenStandardsReference Repository M)an

n, 5 one, keratinocytesisolated hairfromfollicles obtainedfrom - epitestosterone , hydroxytestosterone microporous membrane of a cellcultureinsert (Costar, preparation - TM Aldrich (StQuentin Fallavier, France): 2 d sodium hydroxide, from Sigma  (Cergy Pontoise, France). The antibiotics - androstan

( 2014 - system (Millipore, Saint QuentinQ en s . Dulbecco’s . Modified EagleMedium )

. , androsterone,, epiandrosterone, 5 - - usedin the culture media) were Figure1 androstene 3,17 , 11 - dione and5  - A hydroxytestosterone, 16 he t

providesoverviewan of the - 3,17 issue - in, 4 dione,  s

- were engineeredwere - androstan Aldrich; - androsten  - - 3,17   - - -

- in -

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied ortopical). For qPCR analysis, Figure1 Measurement 2.3 of XME expressiongene differentiation, thetissues treated were with the orderto differentiate thetissue keratinocytes againwith post recons iodothyronine, insulin, amphotericin andantibiotics and was refreshed th (3:1)supplemented with10%fetal calfserum, epidermal growthfactor, hydrocortisone, adenine,tri incubation fo to3 weeks C theon underside. Corning, New York,withUS) post cDNA bycDNA using the RT quantif Total RNA was extracted usingthe RNeasy Mini kitandRNAse The panelofthe metabolic genes isavailable inFigure2 I Phase and II Phase enzymes r RHE DMSOwere (systemic application) or acetone (topicalapplication). RNA w nmoles/cm medium) were ealtime reverse transcription PCR (RT models24 h after truction ied using a using a NanoDropied ND B Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040

2 providesan overview ofthe cultureofthe RHE tissuesand the incubation conditions (systemic

before ORS cellharvesting and long

3 r - reached , the , MC EROD EROD activity and - treated treated with mitotic humandermal fibroblasts

in acetone acetone into thesurface of the model (topicalapplication). Thecontrol treated tissues ORS

ACCEPTED MANUSCRIPT Follicles

confluence, at which 2 xenobiotic

first strandfirst kitandAppliedan Biosystems2720 ThermalCycler. The I Phase and cells were cellst were ACCEPTED MANUSCRIPT Comment citer cedocument: 25

were submergedwere inthemedium after 6 daysand then cultured for a further were  and ORS M - 1000spectrophotometer. Then,1µg RNA was reverse transcripted into qPCR treatment andthe XMEgeneexpression profiles measured were using - 

toform mitotic humandermal fibroblasts used accordingto thesupplier’s instructions Biosciences/Qiagen). (SA hen plated a at density of0.3 x10 NF inDMSO to themedium (systemicapplication) orof 100 - RHE

- / q metabolismassays, respectively) contained DMEM/Ham’s F12 PCR). ThePCR). RT models(1.12cm a a functional timethecultures

-

em storage term test test compoun as as describedpreviously

stratumcorneum 2 and

Profiler™ PCR Arrays humandrug metabolis 2

were liftedwere totheair culturedwith1ml supplemented DMEM . . The growthmedium (0.5ml and1 ml per d Supplementary .

- free free DNAse I set (Qiagen), and

6 sub

cells/cm After . ree ree times a week. - ( lethally treated withmitomycin insert as as isolatedfromthe Table1 2

on cellon culture inserts, picture) 2 – - mediuminterface 3days, .

. After After 16. days of

the

F or tissue ORS - in - m 2

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied of100 nmoles/cm toxicconcentrations of ORS 2.4.1 2.4 reconstructed epider 32or more was considered to be absent,and fold  genes includedwere the in arrays for normalization of data,and theresul A dissociation curve was performedwith a 10min preventingnon RT II arrays runwere on a iCycleriQ™real timePCR dete pmoles). The resorufin was quantified directlyfrom culturetissue medium accordingto a standardcurve (1.5to 100 (Bertholdtechnologies) witha out ethoxyresorufin ORS approximately 5 met application). The Figure(see 1 C 2

for1 h at 37°C, after which hod

SYBR GreenMastermix containing HotStart DNA Taq Polymerase enzymeand other components T - -

RHE model RHE methodbyusing the web Induction ofethoxyresorufin O Measurement of functional XME activities previously Finally, t Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040 models(0.33cm B ). The ). controltissues were -

pcfc amplificationspecific and dyesreference was mixed with andcDNA amplification was - and µM33dic od higherfoldthansystemically applied doses were placedwere 24in a concentration 2 he epidermis was lyzedbyusing

describedby 3

ACCEPTED MANUSCRIPTmis(n=3 - MC  NF ACCEPTED MANUSCRIPT

hot start at 95°C,followed by 40 cycles of 95°C for andsec 15 60°Cfor 1 min. performed Comment citer cedocument: or 2  (25

- culturedwith0.5ml supplemented DMEM medium) were 8) depending8) the xenobiotic treatments. NF in acetone tothesurface of the model(topical application, 10 n o s -

fluorescence umarol fr e based RT Harri  ofthe xcitation M) after after runningcheck tosingle amplicons. Assays for 5housekeeping - or wellplatecontaining 500 EROD µLb s s et al.

inducers 3 - - eshlyprepared inculture medium. The 2 treated treated with deethylase deethylase (EROD) activity MC

wavelength Profiler™ PCR Array Data Analysis Software. A C

( was was measured 2002

(10

andthe - changes of changes geneexpression obtainedwere fromdifferent )  0.5M such that,, a on M)

DMSO (systemicapplication) or acetone (topical of measurement of EROD activity w inDMSO to themedium (systemic application) or

530andnm emission NaOH andNaOH totalproteins assayed were with a

using . . After incubation and PBS washing, the ctionsystem (Bio a multiplate reader Mithras  g/cm

ts were analyzedwere ts according to the 2

basis,topical dosesare uffer containing 3.4µM wavelengthof incubation was - Rad laboratories).Rad The treated treated with ere LB940

based on the 590n T

 valueof carried l/cm m.

non 2 ) -

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied Radioactivity inculture media, studied. washed twice withacetonitrile or methanol, accordingto Jacques et al(2010 a,b and 201 Compartmentalanalysis 72h testosterone Test compounds were 4 ORS 2.4.2 donors expressed microBC Assay proteinquantification kit(Up (PerkinElmer metabolites separated by HPLC monitoredwere by on B,55 5 The temperature ofthecolumn was maintained at 30°Cand the following gradient was used: m KromasilC18guard precolumn (Interchim,Montluçon, France). Mobilephases used a at flow rate of1 250×4.6 testosterone. spectrometry counting

nmoles/cm min 10% l /m

n consistedofin, A: ammonium (20 acetate and -

RHE min 62.5% , each , Test compound Test metabolism

mm,(5 ORS were storedwere at toobtain themass balance Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040

inpmole/mg/min. Experiments performedwere from

models( B,15 , B(a) 2

, For 3 assayed

- , France). France). The, scintillation cocktailused forFlo as as of test compoundin acetone (10 RHE

µm) ZorbaxSb

previously described Jacques by et al

P P and B,65 min18% - MC,

1

ACCEPTED modelsMANUSCRIPT . intriplicate (n=9) 12 radiolabelled

min100% – 3

the HPLC system for B(a)P, 7 ACCEPTED MANUSCRIPT - cm 20°Cuntil theanalysis. incubationsAll carried were out in triplicate Comment citer cedocument: MC

B,28 2 were were

culturedwith to reach toreach the requireddose RHE 4 -

C18column (Agilent, Interchim, Montluçon, France) protected by a Briefly, ). min 18% dissolved B, 70B,

. 7

- Metabolitea extracts extracts andwashing solutions - EC andtestosterone

EC,

.

min100%

testosterone, B(a) consistedofAgilent an 1100(

B,33 1 ORS in

mlsupplemented DMEM medium) were

4 tima 

mM,pH3.5): ml l/cm

- min28% RHE nalyses

)

ofsolvable

B,71 inorder tomeasure EROD activity 2 ) to the surface tothesurface of ) the model(topical application). - .

i e radioactivityline detection usin surfaces, tissue cultureinserts and wells were (2010a .

the Culturemedia were min0% was B, P P and physicochemical propertiesof

- 40 acetonitrile(95:5; v/v) and acetonitrile.B: one detection onewas Ultima FlowM (Packard

performed

performed byradio overnight and epidermis reconstructed from

min 28%

B, 75B, 3 , - band MC was was quantified byscintillation

min0% were were , France) France) associated with a

(PerkinElmer, France).

adjustedwith unlabeled 2014

as as previously described B,45 collected at 24 )

for B. min40% - HPLC HPLC andmass

3 7 H - EC,B(a) - g a Flo a g incubatedwith , which was whichwas , 3 - MC andrelated

B,50 the

h - 0 (3 , 48,h 3batch one A500one

P P and min 0% compound

min 50% batch

7 by -

EC, and

1 e

s to B, ) .

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied 3.1 3 25 follows: vaporizertempera temperature state as follow: were mobile ph (NI spectrometer (Thermo Fisher, Les Ulis, France) operatingeither with an ESIsource inthe negative mode Mass spectrometry analyses chromatographicpeaks Instruments).3 Ingeneral, 3.1.1 <25Ct assignedto samples for which the was Ct>32; whereas, strongexpression was assigned to samples with according application markedly affected the basal expressionofXMEs. after topicalapplication expressionof treated systemically withDMSO the inmedium (1% v/vfinal concentration) The basalgeneexpression of

au,

- ESI) oran sourceAPCI inthepositive mode(PI

optimized and and Phase I Phase expressiongene Basal of XMEs Results ase to the ase source. ion Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040 to capillary temperature 150°C the variety of genes expressed as well as theirwellgenesvarietyexpressed astheas of “ 36 weak the Ct P 0°C.Optimal sourceAPCI conditions - hase hase MC andmetabolitesits quantifiedwere integratingby from reference frommetabolitesreference of

sprayvoltage XMEs and discussion ”

(which are to ACCEPTED MANUSCRIPT I XMEs after systemic application ofDMSO was very similar tothat in “ . good

of10 ACCEPTED MANUSCRIPT Comment citer cedocument: ture 30 ” of3

a a panelof  expression between were these values - also 2.5 l/cm Since reference Sincemetabolitesreference of3

- MCmetabolites 0°C

kV, sheathkV, gas (N2) shownin Supplementary Table1) 2

acetone, suggesting that neither , Corona, discharge4

.

gene targets for TheOrbitrap mass

benzo(a)pyrenewhich display similar structures, and were w ere ere performed using an LTQ

- were obtainedwere from3 APCI). A APCI). post columnsplit diverted 25% of the 30 P hase hase

a µ analyzer u , sheath gas A,(N ,auxiliary gas (N Thelevels ofexpression assignedwere levels of expression lowerlevels wereglobally of I - and MC,operating

. II was was used a at resolution of whereby;

the

the ofarea theradio metabolism in solventthenor method of - MC analysis is 2 ) 2 ) 60 shownin Figure2. “ ESIsource parameters - 1 no expressionno Orbitrap 0

au,auxiliary gas (N

au ORS

and were set were as - - - XLmass capillary RHE ORS ” 60 was was models - RHEs ,000. Basal 2 )

a

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied carboxyl expressed inORS tohelp toprotect the epidermis agains dehydrogenasethat reduces quinoneshydroquinones.to is It present in nativehuman andskin thoughtis findingconsistent is with itsexclusively dermallocation (Luu notpresent these inORS contributing tothegene expressionprofile, itappeared not toso do since selectively dermalXMEs were feeder layer theon underside ofthe andconsequently itwas notinvestigated toxicological assays, especially isoformin (Cheungskin et al. 1999) was expressed intheRHE tissues,namely ADH5. It has beenreported thatADH3 isthemost abundant S theexpression levels of10 ALDH isoforms goodwere (Ct CYP thanfornonfor CES2is responsible forthe majority ofhydrolase activity in thesmall intestine bothbutare responsible expressed at goodlevelsin ORS a theinactivation of ofpolycyclic aromatichydrocarbons e.g. 3 The gene was also dermisbutis presentalso in theepidermis ofnative human skin(Luu expressed inthe ORS Maibach,2007)and D, ESD (esterase involved inthe recycling of sialic acids) alsowere strongly expressed lso upplementary Table

expression playan - terminal esterase L3, esterase terminal a protease), HSD17B10(hydroxysteroid 17 Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040

inthe basement epidermis andinvolved inestrogen formation/metabolism (Hikima and

important expressed

ofthis XMEis important sincethis enzyme involved is inthe detoxification viahydrolysis -

many ACCEPTED RHE tissues. MANUSCRIPTIn addition tothegood expression of dehydrogenases, UCHL3 (ubiquitin - 1 RHE tissues. Theepoxid - role ). ). Int CYP phaseI XMEs CYP

ACCEPTED MANUSCRIPT - reactive epo at low(DMSO treated) and good (acetone treated) levels in ORS Comment citer cedocument: RHE models.For example, ADH1B was absent intheRHE modelsa inthehydrolysis large ofa number of structurally diversedrugs erestingly, onlyoneof theseven alcohol dehydrogenases (ADH) measured for - RHE models.CES1 is themajor hydrolasein human liver organ:thereference liver, ADH3 transwellinsert,which couldthe affect overall XMEprofile by . Sincethe . genepanel xides(Benhamou et al.,1998). Carboxylesterase 1and 2 (CES1/2) t oxidativet stress (Hu et al.,2010). TheNQO1 gene was alsowell

inour - MCand B(a)P especially,, aldehyde dehydrogenases (ALDH), forwhich e hydrolase e gene, EPHX1, preferentiallyis expressed in the stu dy . Although .

, viahydrolysis, andplays an important rolein used <30)orstrong (Ct - The et Theet al.,2009). NQO1 a iscytoplasmic was was developedtheby our skin

was was notincluded inth - The et Theet al.,2009).Likewise, this model includes a fibroblast - 

dehydrogenase 10 <25)(Figure2A supplier for - ; whereas, whereas, ; RHE t . They . is

genepanel

and nd this ) , issues.

were were

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied (Oesch et al.,2014) whereby only halftheskin samples from individuals et al.,2009). In addition, there is a largevariation in thepresence of FMOs inskinfrom different donors, affected by theculture me 2009, Hu al., et 2010). Th FMO5 lacking fromtheRHE models because they are selectivelyexpressed in dermis(Luu isof interest because Of theFMOs, only FMO4 was found tobe present inmeasureable amountsinthe ORS Accordingly reportedto be present innative forthe metabolism of multiple drugs(Imai, 2006 humanskinfro CYP1A2 Wi strongly expressed isabsent nativewhole skin tested by two other laboratories (Luu andCYP27B1 CYP2C18,CYP2R1, CYP2S1, genes foundtobe The e gand et et gandal., 2014 expressionofthe were were have beenreported to in and Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040

adulthuman liver (Edwards et al.,1998)) present inmeasureable levels (weak orgood expression). These absent orpresent onlylow at levels inORS the CES3 was not ,expressed inORS - The al., et 200

CYP2B6; ho . mone or more donors

In previous studies,a . c

). ACCEPTED MANUSCRIPT ompared to other CYPs

itreportedis thatthe expressionofthis XME is differently intheepidermis and

majorhepaticphase No ACCEPTED MANUSCRIPT Comment citer cedocument: e e expression tableXMEs thatnot were expressed dium andthe presence of fibroblasts which inhibitorysecrete factors (Luu w , Wiegand9, ever, these be presentbe inthe epidermis orwhole ofskin nativeskin (Luu

human whereas, skin; the expression of CES3 is lacking (Huet al.,20 CYP3A5, CYP4B1, CYP4F2,

llof these genes, except CYP4F11, (Luu

of FMOgenes inRHE models

et et al.,2014). While it were were - I The et Theet al.,2009, Wi XMEs, innative human skin all reportedto be

andOe - RHE models.

namely , havebeen shown CYP4B1 andCYP26B1 - - The et Theet al.,2009, Huet al.,2010). RHE modelsand sch

CYP1A2,CYP2C9/18/19 and CYP3A4 were

et et al. i s expected that s FMO2and FMO3 would be CYP4F3, inORS expressed e

gand et al.,2014, Hu al., et 2010)

,2007 (Luu -

ha - RHE models toexpress ). Bothof these ). isoforms are also The et Theet al.,2009 thedermis, FMO1,FMO4 and CYP4F11,CYP7B1, were were s

at at verylow levels only few ofa the

previously been genes foundto be were were

FMO1,3,and 4 mRNA

were CYP1B1,were were CYP1A1,were also reportedto be - , Hu et Huet ,al.,2010, RHE tissues. This CYP1B1 - expressed in

reportedto be Theet al., 42 in CYP26B1 native target CYP

,and

(which - The 10 ). ).

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied mostpredomina There a were number of Phase 3.1.2 namely CYP3A7. skin, CYP3A5 was present (Hu al., et 2010) presentbut in native skin three from dependentsince itwas absent in the exhibits anyfunctional a substrates (e.g. bufuralol ordextromethorphan) should be conducted to determine whether thisCYP polymorphic nature ofthis CYP. modelsand al.,2014) 2014) reportedwere notto expressedbe functional activity (NATs) isknown to organbe modelsderivedfrom thesame donors(Wiegand al., et 2014) modelsuchthat it absentis in monolayers ofkeratinocytes but present 3D inepidermal andfull theUGT1A family genes ( skinbyothers (Hu al., et 20 microsomal reactive oxygenspecies such as that inducedU by expressed inthe skin ( ,TST)

. targeted undetectable.were CYP2D6 Phase Phase . In our. study, it Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040 notall at GSTisoforms alsowere present at goodlevelsandhave shownalsoto be present innative II n l expressed GSTtly was GSTP1 is

XMEs were allwere expres

notably, UGT1A10,,which is repor intheskin ACCEPTED MANUSCRIPT reported to a be majordetected CYP at the mRNA levelinnative human skin (Oesch et

intopically applied acetone models (Supplemen Luu ctivity ACCEPTED MANUSCRIPT

Comment citer cedocument: was was - The et Theet al.,2009). 10 at at the gene level - thanin the ). specificsuch thatNAT1 ispreferentially expressed foundtobe expressed II

intheORS SULT2B1,SULT1A3, sed ingood XMEs expressed at significant levelsinRHE tissues (Figure 2B). To further check this, additional s Allthe inhuman keratinocytes (Saeki et 2002) al ORS

liver - UGT isoforms RHE tissues measured here andin native from skin one donor - RHE models inthe AllGSTsplay importantan role levelsin ( Hu al., et 2010

(Ct =20.5) VB VB radiation(Black et al.,2008). ORS

only at verylow levels inDMSO other ORS ted to be affected tedto affected beby the complexity oftheskin - SULT1E1

. RHE modelsnotbut thefetal form of this CYP, foundbe to CYP3A4 expression appears to donorbe - donors(Wiegand et al., , whichis, known tobe themajor RHE tissues;however, other . , The expressionofN Kawakubo et al.,1990) and tudiesusing CYP2D6

expressed ay Table1) tary t hiosulfatesulfurtransferase

andhuman skin (Oesch et al., inthedetoxification of inourmodel belongdo to andhas a higher . . This - acetyl transferases 2014). Other - . In keeping. with treated treated control

mayreflect the - selective cytosolic and As As withnative GST - thickness isoform

The

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied application A comparisonofthee 3.2 intracellular concentrationof polyamines andtheir transportout of cells. SAT1catalyzes theacetylati also noteda strong expression of thisenzyme neurotransmittersand catechol hormones ( been suggestedto a be surprising occurrence since acetyltransferase 1) involved inlipid metabolism), COMT Other strongly expressed this, the presence of NAT1 thebut absence of NAT2 inthe verylow higherthanthat of CYP1B1 (13.9 inducedby both significant (P<0.0 notably CYP1A and CYP1B1 isoforms (Liet al.1998). Inline withthis, there than a ofCt 32 significantly (P<0.05)increased by Supplementary Table1 The valuesCt

 NF and Induction ofPhase leve Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040  of3 NF (133NF of l

andthestatistical analysis 3 .

- -  MC MC NF and expressionofCYP1A1 5 )

and NAA20 (N(alpha)

ACCEPTED MANUSCRIPT induction ofCYP1A1 and CYP1B1 orsystemic application ( - Figure3 od andfold) xpression . A biologically A relevantinduction ACCEPTED MANUSCRIPT 3 Phase II Phase Comment citer cedocument: - MC on ofonspermidine and spermine, and involvedis intheregulation ofthe I ) XME XME expression . . DespiteCYP1A1 beinge

3 are are ligandsforAhR and of - - MC XME foldby phase more

(25 (catechol

genes ACSL1were (

for control andinducer I than2 insolvent treated tissues 0

 Luu - - XMEs acetyltransferase acetyltransferase 20) NF andNF 11.9 fold) ofβ , suggesting, its presence isimportant in metabolismskin - The et Theet al.,2009) - -

foldandwhen the inducedexpressed level was greater O - to NF in - thisenzyme m clearlymeasurable ORS ethyltransferase), isshown in Figure3 - AhR od byfold expressionin - xpressed at verylow levels response was called when a gene was RHE ORS Long

regulates a numberofPhase

; however, ;we andothers (Huet al.,20 is - models . - 3 treated models are shownin

RHE The . involved in - -

MC chain expression ), which was most duelikely th to was was confirmedin ourstudies. ORS SAT1 levels withthat24 h -

fatty .

- RHE

. Theresponse . (spermineN1 the metabolism of -

was was a selective acid ofCOMT inskinhas models - CoA ligase1,CoA

it was readily after after treated treated with I - topical

was was much XMEs, and

. 10 e )

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied A A numberof comparedto whichcouldbe due to a higher was higherintopical applications) thatless ofthedose penetrated the topicawere was higher more oracetone, respectively).By contrast, EROD activitywas undetectablein ethoxyresorufin incubating controland treated The induction of CYP1A1 andCYP1B1 at the genelevel was confirmed the at functional by level however,the increase notdid reach statistical significance There twowere othergenes that inducedwere bymore than 2 32.34). fold), the induced level was still belowthecut 32.By contrast, although the be systemic MC(2.2 CYP4F3 biologically batches

- ,ALDH3A1 or fold, P<0.05). Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040

lly applied. and moreconsistent across donors topical application

 additional genes inducedwere eitherby treated treated topically andsystemically with NF (25 relevant after the systemicapplication of , ORS

ACCEPTED MANUSCRIPT h24after the start ofinduction -

RHE and NQO1  M)). Thiscould ALDH3A1 ACCEPTED MANUSCRIPT Comment citer cedocument: models

topical

ORS of affinity by both .

likely be treated treated topically orsystemically withvehicle control solvents(DMSO ORS Inaddition, the induction response was higher using - was was significantly induced(P<0.001)

the RHE application of

inducer

of - systemicapplication of RHE tissues with 3 -

dueto difference inthebioavailability of theinducers such MC

- s when applied topically insystemically treated offvalue considered to biologicallybe relevant(Ct was

(Figure3) with forAhR 3 -  MC   theCYP1A1/2 andCYP1B1 selective substrate, NF or NF orNF NF or3 

(despitethelower dose of caused a significant

; however,; theinduction was onlyconsidered to NF measurable in (P (P 3 3 values: - - since fold(ALDH1A3 andCY4F12) - - MC MC MC  NF and . . Theextent ofinduction inducers by

theinduced level was below Ct a of (Figure4) (but

0.323 ORS (even thoughthe by>2 never ORS PTGS2by topical - and 0. and RHE

- . - induction ofthis gene

RHE fold bothbyinducers) EROD EROD activity was tissues than whenthey 245 following the 3 tissues onefrom or , respectively, - MC theoretical 3 -

(10 MC ly

applied 3 by  than M) ,namely,  ) dose NF and/or 

(2.64 NF ; - , -

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied change was namelyINMT T 22nmoles/cm byothers using hepatocytes (Soars al., et 2004). CYP1A1 and CYP1B1 induction (UGT1A3/5/ only controlandinducer 3 A comparisonoftheexpression 3.3 theinduced level was not biologically relevant was (Ct~32 for CYP4F12) was was (whichw h);however, in orderto compare the functional profileswiththe gene expression analysis described above and throughtheskin ( knowntobe metabolized byphase whether Havingchar 3. - here were here were MC 4

2014 also onetype of orsystemic Induction ofPhase Metabolism ofMetabolism substratesmodel after topical application this ) considered as as measured after 24h), we have expressed thedose ) Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040

and thereforewere appliedtopically.

two phase notsignificant fromcontrol treated acterized the Phase 6/7 correlated with 2

andSULT1A1; however, thiswas effect considered notbe to biologically  8 and 10),/8 NF Phase fromthestudies presented here andfrom previous rep application -

ACCEPTED MANUSCRIPT models treated are shown in tobe (3.8 II II - XMEs thatdown were ACCEPTED MANUSCRIPT

4.6

XMEgene Comment citer cedocument: a

and

suitable functional XMEactivities by incubating them with four test compounds - II

fold)than of β XME XME expression theextent ofinduction was alsorelatively low I of . and The lowerinduction ofUGTs - phase I NF

timepoint s and/or

II that is shownisin Figure

aft XMEgeneexpression II was er XMEsin II

topicalapplication of100 nmoles/cm

XMEs.

inducedby sinceall ofthe 4model substrates foundwere tobe -

ORS regulatedby more than 2 T Supplementary Table1 The incubations included three timepoints (24,48 and 72

he responses were - ORS RHE Allfour test chemicals are known to 5 -  RHEs models . NF and TheCt values -

dependent in compared to

ORS withthat h24 .

3

higher - MC - RHE orts( .

effects usingeffects h24incubations. - Unlike Phase ,namely UGTs fold .

and thestatistical analysis

after after systemicapplication of tissues, we then

Jacques Jacques et al.,2010 CYPs has CYPs alsobeen reported

( by topicallyby applied 2 after after topicalapplication -

to 2

3 5 - MC - fold I

XMEs,there relevant

(2.6 ) penetrate

comparedto

investigated - 2.9 a and ba

because because - fold) 3 - MC for was

.

It of ,

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied couldhavecontributed tothe (shownto be 1999), whichalsois levels biotransformation basalCYP2E1 expression was absent in these producedin metabolite into hydroxycoumarin CYP The mainmetabolites of7 3. extensivelymetabolized after 24h (SupplementaryFigure Interesti 24 hthereafter,decreased but most likely dueto further conjugation to HC 7 4nmoles/model culturemedium after 24h. theonly oneof t epidermallayer metabolized 7 theapplied dose The amountof - 4 EC concentrationEC used. .1 1

A1 (Figure2A). 7 ngly, the production ofHC - and CYP2E1 Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040 Ethoxycoumarin the expressed - EC, indicatingEC, radioactivity s HC ORS he 4he test chemicals that was present as the parent chemicalinsignificant amountsinthe

(Figure6A) (forall mediasamples into the medium (HC) (HC) of7

- The glucuronidation of HC isreported to

ACCEPTED MANUSCRIPTglucuronide expressed - RHE -

EC tohydroxycoumarin was likely toCYP1B1,be which was expressed (and a minor amountvia 1 inthe ( A ACCEPTED MANUSCRIPT Shimada . Comment citer cedocument:

- Duringthecourse of the h72incubation, t ) EC The , model

, applied (7 which we attributedtothe shortresidence timeof7 suggestingthattheenzymes inv sulfonation that ORS producedin - re re was a at

EC) and

, with, goodlevelsin (an observationmadealso thischemical does not require metabolism inorder to pass through the et al.,199 (SupplementaryFigure - - s RHE

that was HC ulfate

anddosesof themajor metabolite being HC dose –

of HC. models sulfate andHC recombinanthuman liver microsomes 7 - ), dependent recovered inthemedium ranged between 60%and 82% ORS ),with a subsequent conjugation oftheintermediate

ORS

CYP1A2 ( in acetone treated control models De De Kanter et al.,1999 - - g RHE

7 luc - - RHE EC

uronide linearproduction of three all metabolites models

and 1 over72 h) models. )

. previously be mediatedbe by UGT1A6 (Lampe et al., olved were not were olved saturated

CYP1B1

reached reached a plateau after 48h , the ,responsibleCYP for There were There were he productionhe of HC was maximal at andof this, 54

- ) ) g (Jacques (Jacques et al.,2010 .

luc mediated The - s uronide ulf se ate

- six three metabolites were EC inEC the RHE models, are are

andHC O

(Figure2B SULT isoforms - (Figure 73% was was non73% - producedvia deethylation , evenat thehighest - the g

6) b luc ) ) . )

. uronide which Since at at good This was

over1 to - to

of .

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied 2000;Desilleal., et 1999)) and6 inORS neither of which includedwere ingeneour was androsterone (peak XIII),the production of whichiscatalyzed CYP2C18/19by (Yengi, 2003), hydroxy (peak X), incubatedinourstudies, namely producedhumanby liver microsomes human Testosterone ishydroxylated multiple by and CYPs reduced byreductases 3. formation compared to ofthe suchthat onceit had reached the culturemedium, nof reduced metaboliteswould have beenexpected (van Eijl et al.,2012). expressionpresent reconstructedin skin models compared to reductases capacity, higherproportions of suggesting metabolites there was a more extensiveproduction ofhydroxylated forthe purposeof determining thedose Due tothemultiple metabolites formed, we grouped them as 3,17 andCYP3A43 (Soucek al., et 2001; Waxman et al.,1991)), as well as a reduced metabolite, 5α 4 .2 - in (peak dioneXIII’,catalysed by 5

O - liver Testosterone testosterone (peak testosteroneIII) - RHE tissues,7 deethylated Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040 hydroxy

that . microsomes This was Thiswas alsoo thisobservation was notspecific to - 

ACCEPTED MANUSCRIPT (peak 4.dione metabolite was detected in themedium, whichmay duebe to its

 the faster thefaster rate of its - ACCEPTED MANUSCRIPT andfresh humanskin explants hydroxytestosterone(peak producedI, by CYP2A6 and2A13 (Bogaards et al., Comment citer cedocument: bserved

(Figure7B)  epiandrosterone(peakXII), androstenedione (peak XI), - s hydroxytestosterone(peak II, produced CYP1A1/2,by CYP3A4/5/7

IV in

 (Jacques (Jacques et al.,2 , V’ andVIII

ex ex vivo - andtime dependence reductase reductase ( . conjugation

analysis.Two other hydroxylated metabolites producedwere One ofthemajor metabolites formedin theORS

humanskin our ), Liuand Yamauchi, 2008)) 4 urther

- 0 ORS androsten metabolitesoftestosterone thanreduced in 14) were alsowere 14)produced by the (Jacques (Jacques et al.,2 o the to

andpig ear skin (Jacques et al, 2014) - metabolism was possible. RHE either (Figure7A) HC - 16α model. - hydroxyl glucuronide - Possiblereasons for the ol - 0 3,17 . 14 Considering thelower CYP , producing manymetabolites in An interesting finding was that ated ). Th ). - . dione

andHC or e e majormetabolites

reduced metabolites (peak s

low ORS Comparatively less -  sulfate

6 VI) VI) and16α rate of - - tes RHE - RHE models relative - t androstan . osterone

models -

-

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied saturation was also Non metabolites as decreased the dose increased,suggesting these enzymes saturatedwere at higherdoses. metabolite formation was linear (and saturated)not over thisdoserange, therate ofproduction ofreduced dose modellong enough tosubjectbe to extensive meta present as non theapplied dose The amountof follicle sample although thishas beennot confirmed withexperimental data dueto thelow amount ofcells in thehair theRHE models are derived may exhibit XMEprofiles withpredominant composition ofthemedium and/orthe differences inthe contribution ofhydroxy metabolitesofthis chemical penetrate into the receptor fluid penetration B 3. application. respectively hydroxylated metabolitesincreased 34%from after 24h to 48% and 54% after 48h and 72 h, furtherwithin the incubation duration 24h to 17% and 13% after 48h andh, 72 respectively (Jacques et al.,2014 ( 4 a . ) 3 - P P was selected as a test chemical - linearproduction oftestosterone metabolites

andtime Benzo(a) Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040

ofthis lipophilic compound suggesting, - - s metabolized testosterone dependent radioactivity . (forall mediasamples

pyrene (B previously

ACCEPTED) MANUSCRIPT ) . T ain andreduction lationpathways totestosterone metabolism he that the ACCEPTED MANUSCRIPT

Comment citer cedocument: (Fig

percentage of ( applied a ure ) observedusing pig a skin P) rate rate

7

andSupplementary Figure1 that was of metabolism

because because itsmetabolism ismore complex (Fig

anddosesof testosterone ,indicating thatthis ishighly

reduced metabolites ure recovered inthemedium ranged between 62%and 78% 7

andSupplementary Figure1

dependenton its (including bolism decreased after decreased , suggesting, that these metabolites are metabolized t higherdoses skin .

The production of metabolitesall was

modelmolecule ( hydroxylated metabolites) models J also acques acques et al.,2010 B) over 72h) and of this, only metabolism . decreased decreased

48 W . T

(400 hiletheincrease in hydroxy h following testosterone he hairfollicle cellswhichfrom hydrox –

B) residedin the

800nmoles/skin over time such thatonly the than7 .

T a ylation he could be dueto ).

Likewise, inthe percentage of - EC and because andEC because the

( andenzyme from25% after pathways ORS 2 - explant 4 - % was was % RHE both the - ,

of

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied been ind inducing to effects incubation.is It likely thatCYP1A1 isalso induced sinceB(a)P by this toresponsiblebe forthe CYP foundtobe inpartby CYP2C9, andCYP3A4 (Gautier et al.,1996). produced same metabolites are alsoproduced humanby liver microsomes (Jacques et al.,2010 metabolized B(a)P 1 doses ofB(a)P non current study nmolesincetheamount of these conjugatedmetabolites produced reached a plateau this at dose. (Figure The production ofphase models. was present, alongw conjugating SULT forms(SULT1A1, 1B1, 1A3)and SULT1E1 (Meinl al., et 2008). Of these,UGT1A6 (1A9,1A6 ,1A1) (Bock al., et 1999; Bockand Kohle, 2005;Girard et al.,2008) andthe phenol producedin mono The mainphase C ). ). - metabolized B A -

sulfate conjugatessulfateof B(a)P metaboliteproduced to be ableidentifyto themby lth

8A uced. uh thesyntheough by CYP1B1,by CYP1A1 and microsomal epoxidehydrolase Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040 ). ). Bycontrast, expressed at the gene level solventincontrol RHE tissues ex ex vivo by , only 1 ,

over 72h ORS II metabolites ofproducedB(a)P inORS ( 

a

ACCEPTEDpigskin (Jacques et al., 2010a). B( MANUSCRIPT NF and3 - - ithUGT1A10, 1A7and 1A8, andSULT1A3 and1E1 presentwere inORS ) 3% RHE P . . B , the , majority ofbeingB(a)P metabolized inthefirst h24( ACCEPTED MANUSCRIPT I t of theradioactivity Comment citer cedocument: i sulfation and glucuronidation c standards c etween 9% and 21% of theapplied dose w metabolites of models - - MC;therefore, thismayCYP contribute metabolismB(a)P to once it mediated metabolism ofB( - hydroxide(Figure8B), which consistentis with themetabolites tohydroxide, diol, catecholand dione metabolites

forsome metabolites not were avail

B(a)P

present intheculture medium

was was dose

a)P phase a)P

pathways saturatedwere at doses higherthan1 Of Of these XMEs,only - - mass spectrometryand a)P, a)P, at leastin the few first hours of the dependent and overlinear thedose rangetested RHE models are the mono I metabolites are conjugated viaUGTs s as as recovered

(EPHX) . . T is herefore able an AhRligand withsimilar was was foundto correspond to

there was enough, ofeach CYP1B1 ( Hvastkovs et al.,2007) , these isoformsare NMR. in SupplementaryFigure mediasamplesfor all - a

glucuronide and (Figure8B) ) and )

B(a)P and EPHX

are are

was - mainly - RHE

. The .

we has

likely

re and

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied sulfo intheORS minor role) primarily responsible forthe activation of2 a thebioactivation depend are reported tomorebe mutagenic than thehydroxides (Flesher al., et 1998) andtheSULTs involved 1982 including ( metabolized 19 The metabolic activation of 3. the appliedthe(for all samplesdosemedia of78%55% recoveredand medium inbetweenamount radioactivity theranged ofwas The that applied detected after 24 hincubation of UGTs chemicals,UGTswe RHE modelslikelyis to be viaEP “hydroxylated” position of produced(<10 functiona which are then ctivation of 1 4 92 . 4 - The hydroxides ). The primary ). metabolites are hydroxides,

conjugation of

at at the genelevel after topical application Figure (see 5B) 3

activities.l - Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040 trans Methylcholanthrene ( - hydroxyl functionsnorconjugation positions,

(Meinl et al.,2013). RHE models(Figure 2B);which sugge to -

hydroxy  1 - or hydrolyzed 1,2 - g) keto “ - i notdidallow NMR identification reduced ACCEPTED MANUSCRIPTdihydroxy T 3 re re involved inthe detoxification of - - he s he 3 MC - are are subsequentlyconjugat - 3 MC ACCEPTED MANUSCRIPT - onthehydroxide Comment citer cedocument: MC tandard syntheticmetabolites of

metabolitesor ” (Stoming et al.,1977) 3 , andcholanthrene, (Myers et al.,1989, Myers andFlesher, 1990)

- metabolites -

MC C and MC (withSULT1C3 and SULT2A1playing a minor role) Interestingly, noneof these SULT 3 3

- H to DNA toDNA reactive metabolitesinvolves CYP1 XMEs(Guengerich, - MC MC X1,which was expressed in these tissues.

with 11,12 ) and dosesandof

(Figure

that conjugated - hydroxyl ORS -

the genelevel of these XMEs does correlatenot withthe dihydroxy

1 ed tosulfates ed andglucuronides. - - 9A ) hydroxy RHE

sts either sts and number a of multiple T . ) - SULT1B1:is primarily responsibleforthe . 3 3

Detoxification viaepoxide in the - - herefore, MC MC models 3 - - 11,12 anddenoted these metabolites MC - 3 3

- - over 72 h) and ofoverthis,h)72and 2 only that (withSULT1A2and SULT2A1 playing a MC MC . Interestingly, . we foundthat3 - . dihydro and these conjugates readilywere

anotherisoform was responsible

we we were notwere available and2 isoforms presentwere at the gene level were unablewere to - - MC

hydroxy -

yrxltd metaboliteshydroxylated (ClaytonandClayton, 1981

As As with - Thesulfated conjugates 3 - MC determinethe andthelow amount andSULT1A1 is

theother 3

, which are whichare further , either as - 4% was presentwas4% - , MC induced

forthe exact ORS in ,

ir

-

-

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied nativehuman skin formetabolism assays. functional activity ofXMEs in ORS The aimsofthe currentstudies towere makea comprehensivecharacterization ofthe geneexpression and 4 be expressionUGT of and ORS inlevelthe gene thepresent at were ( SULT1A1 and SULT1A2SULT2A1, 2005 3 lower was than theformation metabolitesofII phase (glucuronidated and sulfated conjugates) metabolism. extensive to subject non as models, made undercontrolled conditions, are available andare suited to flexibility inthe selection ofthedonor, especially ofadultan skin source. available epidermalmodels metabolism capabilities. Therefore, it Our modified tomimic topicalapplication. test systems metabolism underrealistic conversely,detoxification applied chemicals -

MC en

UGT data showeddata thatt and ) reported as UGTenzymes involved inthe 3

phase - metabolized Conclusions 1A7 were present at the genelevelpresentwereORS 1A7thein at the Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040 h SLs novd n h frain f h slae cnuae ae SULT1B1, SULT1C3, are conjugates sulfated the of formation the in involved SULTs the ,which notdo havethe same XMEprofile I metabolites are conjugated viaUGTs ( by ACCEPTED MANUSCRIPTand1A8 3

enabling the measurement of - he ORS MC ACCEPTED MANUSCRIPT bothof which impactlocal ,skin andsystemic toxicity Comment citer cedocument: idctn ta ti chemica this indicating that , exposure . The use . oftheORS ofthehair follicle a isnew concept UGT - RHE model 1A10 were induced after topicalinduced 1A10wereapplication after 3 of In

-

scenarios andis RHE modelsand to characterize them as a suitablesurrogate for contrast could

Meinl Such models - is a a isnovel alternativeapproach RHE models2B).(FigureRHE represent a to t l, 03. Interes 2013). al., et

- B(a)P MCmetabolism intheliterature.

potential metabolic activation to toxic metabolites or, more predictive 1A5, s

,

complementary are importantin thesafety assessment of topically as as skin(van Eijlal., et 2012) andcannot be h frain f hydroxylated formationthe of - be to enough models long RHE the in resided l RHE modelsRHE2B)(Figure 1A7,1A6 igy nn o tee UT isoforms SULT these of tingly, none

, 1A1 , than

Of the UGTs,UGT the Of modelto a number a of different assay Although several commercial ) ( ) usingextrapolations from withextensive xenobiotic Gramsal, 2000et Thismodel .

the commercially

. Interestin. - MC but they have nothave theybut MC that 3

allowsfor - MC 1A gly, the genethegly, refl 5 , Finelal, et metabolites

, (Figure6) UGT ects ects skin liver 1 A6 A6

.

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied I keeping withthe proteomic profileofXMEs in nativehuman reportedskin by van Eijl et al.(2012). Alsoin XMEs was observed foralcohol/aldehyde dehydrogenases and(carboxyl)esterases, whichconsistentis humannative skin at byleast oneother laboratory (Table presence of many of1 thePhase and 2XMEs the at genelevel which havebeen reportedto be present in The ORS wanted to themost relevant exposure routefor cosmetics ingredients therefore excludedifferences in variations.Therefore, a was affected by tr notedthat endpoints (e.g. genotoxicity (Aardema et al.,2013), skincorrosion/irritation (Alépée et al.,2015)), we there et al.,2012), making them unsuitable for skin penetration studies idealbut for metabolism studies(since sufficient.RHE models exhibit a reducedbarrier function compared with thatof native skin (Thakoersing application toepidermal models, and the epidermal XMEs There islittle information onthecomparison ofmetabolites produced after systemicand topical CYP27B1. includ detected innative human skin anddehydrogenases ( Although theabundance ofenzymes CYP reportedis tol be XMEs,thehighest levelsobservedwere for GSTs(notably GSTP1)and transferases, including NAT1.

ishigher a bioavailability totheepidermal XMEs).Our study demonstrate thatthe ORS : CYP1B1, CYP2C18,e:CYP2D6, CYP2S1,CYP3A5, CYP4B1, CYP4F3,CYP7B1, CYP26B1 and

withothe - Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040 RHE model a has XME profile that consistentis with nativehuman skin, supportedtheby investigatewhether XME

s

profiles morewere sensitivethanother endpoints such thatthereproducibility ofdata ACCEPTED MANUSCRIPTreports,r the levels ofII Phase ansportationconditions e.g. the duration oftransport our tolaboratory andseasonal vanEijl et al., ni ACCEPTED MANUSCRIPT mportant aspect aspect mportantof thismodel that is iti Comment citer cedocument: . These . topical application enzymesfunctionality 2012),t

were

also present at thegene levelintheORS here here are a numberof that CYPs are consistentlyreportedto be

XMEs presentgenerallywere higher than those of Phase of test compounds

dueto t

1

andfor environmental pollutants ). ). bioavailabilityoften is assumedto be ow incomparisonow with reductases, esterases The hi ransport effects. ransporteffects. s generated ghest gene ghestexpression levelofP couldemployedbe insuchassays.

withinlaboratory our and Since topical application is - RHE mode w , - RHE e e also l

and hase hase

I

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied metabolism studies;therefore, itisimporta CYP1A1 and CYP1B1 (133 a numberofUGT isoforms, although theextent ofinduction was not as high (2 notnecessarily evident at thegene andenzyme acti bothpotently induced bytheAhR ligands, Importantly, regulation ofandrogens in theskin which in regulatesturn hairgrowth thepresence of reductases such as 5 testosterone, andB(a)P 3 presence of functional was confirmed CYPs theby extensive CYP model, whichisin keeping withtheinternal exposuremodel described by Manwaring et al., (2015). ofmetabolism was dependent theon metabolism bythe XMEsa models are ableto rapidly metabolize four test compounds within 24h of topical application. The extent potential toxicity appl metabolism studies demonstratedusing eithersystemic or topical application routes, supporting their use incutaneous inducedinresponse totheincu Inconclusion, the ORS compound itselfi.e. i c ation Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040 Such a reproducible model . will of be interest when evaluating

the model isresponsive prototypical to inducers, CYP such that CYP1A1 and CYP1B1 are correspondto functional of innovative

and theprediction ofthesystemic exposure toparent

ACCEPTEDautoinduction MANUSCRIPT - RHE model ACCEPTED MANUSCRIPT - Comment citer cedocument: MC. TheMC. metabolitesproduced after incubation withtestosterone also indicates -

and25 ain with bationinducers.CYP Functional XMEactivities also were active

during the incubation period tested vity level, whichofissignificance because gene expression  0 was was shownto express a numberof important XMEs,which -

- ingredients od respectively).fold, reductase, the presence of which isreported to pivotalbe inthe

XMEactivities ( t toconfirm ntwhether XMEs can alteredbe b  NF and .

3 - MC Hewitt et al.,2013) TheORS . Theinduction ofCYP1A1 and CYP1B1 was ndtheon residence timewithin the - mediated metabolism of7 - RHE models are intended foruse in .

(LiuandYamauchi, 2008). and/ the . ormetabolites af  cutaneous - andNF

to5 - fold) as thatfor 3 y the y test - metabolism and MC -

also inducedalso EC, ter ter topical levelsdo could

The be

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied microsomal epoxidehydrolase genotypes. BenhamouReinikainenS, Bouchardy M, C, Dayer P, Hirvonen A. Association between lung andcancer genoprotection. J S,Duplan new H. A hair follicle Bacqueville D,Douki T, Duprat L, Rebelo 29(7):1779 andirritation assessment usingSkinEthic™ Reconstructed Human Epidermis. A assay. requiring metabolic activation in the EpiDerm™ 3Dhuman reconstructed skinmicronucleus (RSMN) Aardema MJ,Barnett BB, MunGC, EL,Dahl Curren RD,Hewitt NJ, PfuhlerEvaluationS. ofchemicals References Toxalim, Toulouse, France), MetaboHUB MetaToul the by achieved were analyses spectrometry Mass during thepreparation of themanuscript. The authorswould liketothank Acknowledgement aldehydedehydrogenase enzymesin skin.Biochem Biophys Res Commun. 1999 CheungC, Smith CK,Hoog diffusion cell.Pharm J Sci.1985 BronaughRL, Stewart RF. Methods forin vitro percutaneous absorption studiesIV: Theflow 10.1080/00498250010021684 rabbit, dog,micropig, monke Determining the animalbest model for humancytochrome P450 activities: a comparisonofmouse, rat, Bogaards JJ,Bertrand M,Jackson OudshoornP, MJ, Weaver RJ, vanBladeren Walther PJ, B. MethodsEnzymol. 2005, 400:57 BockKW, Köhle C. UDP 1999, 31(2):411 transcriptional regulation of PAH BockKW, Gschaidmeier H,Heel H,Lehmköster 2008, onantioxidant expression undifferentiatedin anddifferentiated mouse keratinocytes. Carcinogenesi BlackAT, Gray J lépée N,Grandidier MH,Tornier C,Cotovio J. Anintegrated testing strategy forin vitroskin corrosion 29(1):219 Mutat Res. 2013, Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040 - 92.

- - doi: 22. 25. P, ShakarjianMP,Laskin DL, Heck DE,Laskin JD. Distinct ofeffects ultraviolet B light Photochem

ACCEPTED MANUSCRIPT10.1016/j.tiv.2015.07.012. 750(1 ACCEPTED MANUSCRIPT Comment citer cedocument: - glucuronosyltransferase 1A6:structural,functional, andregulatory aspec

JO,Hotchkiss SA. Expressionandlocalization ofhuman alcohol and y andman. yXenobiotica. 2000,30(12):1131

Photobiol B. 2015, - 2):40

- - , NicolaHewitt forproviding writing supportand for helpful discussions eie humanepidermal derivedmodel for theevaluation ofsunscreen 75. - 74(1):64 inducible humanUDP -

9.

doi:

- - - Cancer Cancer Res. 1998, 7. ANR MoreiraGuiraudS, B, Dromigny H, Perier V, Bessou 10.1016/j.mrgentox.2012.08.009.

- 151:31 11 T,Münzel PA, Bock - INBS - 8. - glucuronosyltransferases. Drug Metab Rev. - 0010 doi: 10.1016/j.jphotobiol.2015.06.015. 58(23):5291 .

- XO pafr (NA UMR1331 (INRA, platform AXIOM - HennigBS. Functions and - 2 doi: 52. - 3.

Toxicol , 261(1):100

In Vitro.2015, - 7. - through

- Touya s. s. ts. ts.

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied 1998, 243: sulf FlesherJW 315(3):1143 glucuronosyltransferase1A5:identification, expression, and activity. J PharmacolExp Ther. 2005, Finel M,LiX, Gardner 10.1111/exd.12694. hum Eilstein Pharmacol. 1998, 56: antibodies againstthemajor xenobiotic metabolizing forms ofcytochrome P450 in Edwards RJ,Adams DA, Watts PS, DaviesD 27;68(10):1437 activity inpiglivers and hepatocytes intendedforxenotherapy. Transplantation. Nov1999 DesilleM1,Corcos L, test system fordrug metabolism inhumanliver, lung and kidney. De Kanter R, Olinga De P, Jager MH, Merema MT,Meijer DK, Groothius GM. 1994, dependentdifferences inthe mutagenic profile ofits ultimate carcinogenic m ConneyAH,Chang RL, Jerina DM, Wei SJ. Studiesonthemetabolism ofbenzo[a]pyrene and dose Toxic Clayton, GD,FE 36:485 keratinocytesand comparison withtwo commercially modelsandnative H,Bessou Guiraud B,Hernandez Guengerich FP. tract. Arch BiochemBiophys. 2000, 377(2):255 inducibility3 by GramsB, Harms A,Braun StrassburgS, MannsCP, MP,Obermayer EpiDerm™models. Xenobiotica. 2012 Effects of thegenotoxic compounds, benzo[a]pyreneand cyclophosphamide onphase 1and 2 ac GötzC,Hewitt NJ, Jermann E,Tigges J,Kohne Z, Hübenthal U,Krutmann J, Merk HF,Fritsche E. MutatRes. 2008, 644(1 UGT1A9functional variants,meat intake, andcolon cancer, among Caucasians and African GirardH,Butler LM,Villeneuve L,Millikan RC, Sinha R,Sandler RS, Guillemette C. UGT1A1and fromtheanalysis ofindividual kineticsteps in Gautier ate ester ester 1 ate an skinandanSkinEthic reconstructed human models.skin Dermatol.Exp lg. 3rdology. ed. New York:John Wiley Sons 26(1 –

9 JC Léreaux J,G, Arbey E, Daronnat E,Wilkinson Duché S, D. Xenobiotic metabolizing enzymes in 3. Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040 - 30 - ,Urban Touya S. ,Horn 2):125 - doi: 10.1111/ics.12150. . doi: 10.1124/jpet.105.0919009. – - 35. sulfooxy - 43. Metabolicactivation ofcarcinogens. Pharmacol Ther -

methylcholanthrene offamily 1UDP Clayton (eds.). Patty's Industrial Hygiene and Toxicology: Volume2B,2A,2C: - J,Lehner AF. 63.

P

ACCEPTED, Beaune MANUSCRIPT Characterisation of a humanepidermis modelreconstructed fromhairfollicle 377

- L'Helgoualc'hA,Frémond B,Campion JP, Guillouzo A,Clément B. Detoxifying - - - PigeonH, Ceruti I, Mas S,

StephenD, Bratton MackenzieS, PI, Radominska 2):56 – 3 ACCEPTED MANUSCRIPT Comment citer cedocument: - 387. methylcholanthrene

- P 3 doi: 10.1016/j.mrfmmm.2008.07.002.63.

,Pompon

Carcinogenicity of1

,

42(6):526

D.

Simulation ofhuman benzo[a]pyrene metabolism deduced

S,B

recom - inSprague 5 doi: 10.1006/abbi.2000.177765. ,1981 - oobisAR. Palvadeau Y,Saint 37. binant yeast, ChemRes -

hydroxy - doi: 10.3109/00498254.2011.643255. - glucuronosyltrans 1982 - Dawley rats. Biochem

p.3367. , Developmentof - 3 -

ehlhlnhee andmethylcholanthreneits electrophilic Toxicol InVitro. 1999, - - . Straub P.Distribution and

MartoryCastex C, 1992

skin.Int J Cosmet intheferases gastrointestinal rat - Pandya A.HumanPandya ,

2015, etabolite. 54:17 a a comprehensivepanel of Toxicol.

Organ slices as in an vitro

Biophys - 61. 24(7):547 humans.Biochem

1996, 9(2): DrugMetab R - 13(4 RizziN,Duplan Res

Sci.2014, - -

Americans. 9. UDP -

5):737 Commun.

tivities in doi: 418 -

- - - ev. ev. 425. 44.

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied Pharmacokinet. 2006 Imai T. Human carboxylesterase isoz 80(6):2272. frombenzo[a]pyrene metabolites. AnalChem. 2007 Electrochemiluminescent arrays for HvastkovsEG, SoM, Krishnan S,Bajrami B,Tarun M, Jansson I, Schenkman 24(5):1450 EpiDermin vitro3Dhuman epidermis modelcomparedto human skin. Toxicol InVitro. 2010 BB,Osborne R,Reichling TD,Nemeth ST,Aardema MJ. HuT, Khambatta ZS,Hayden PJ, BolmarcichJ,Binder RL, Robinson MK, Carr GJ, Tiesman JP, Jarrold hydroxysteroid dehydrogenase inhuman skin invitr Hikima T,Maibach HI. Genderdifferences of enzymaticactivity anddistribution of 17. models forandaccurate ethicalrisk assessment: metabolicconsid EilsteinLatil J,A, Kenny J, Moore KuehnlC, Barroso J, J,Fautz R, Pfuhler S. Use o HewittNJ, Edwards RJ,Fritsche E, Goebel AebyC, P, ScheelJ, Reisinger K,Ouédraogo G, Duche D, reconstructe dependenton glutathione S Harris IR,Siefken W,Beck hormonemeta relateddifferences insteroid metabolism inskinmalepoint a keytorole of sebocytes incutaneous Haag M,Hamann KulleT,AE, Riepe FG, BlattT, WenckH, HolterhusPM,Peirano RI. Age andskin site Lampe JW,Bigler ofarylamines. Skin Pharmacol Kawakubo Y,Yamazoe Y,Kato R,Nishikawa T:H Pharmacol Physiol. 2014, dermaldelivery oftestosterone: qualitative assessment using a new short Jacques C, Perdu JaminE, CravediEL, JP, Mav 10.1016/j.tiv.2010.04.006. of[14C] Jacques C, P approaches. Toxicol Lett. 2010a,199(1):22 Disposition andbiotransformation of14C Jacques C, Perdu DuplanE, H, Jami Toxicol InVitro. andUDP enzymeactivities ofCytochrome P450 enzymes, Jäckh C,Blatz V,Fabian GuthE, K,van Ravenzwaay B,Reisinger K,Landsiedel R.. Characterization of

doi: 10.1093/toxsci/kft080. - ethoxycoumarin in a ear pigski - Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040 glucuronyltransferases inhuman reconstructed epidermis and full -

63. d epidermalmodels.Skin Pharm erdu DorioE,C, Bacqueville D,Mavon A,Zalko D. Percutaneous absorptionandmetabolism bolism. Dermato doi: 10.1016/j.tiv.2010.03.013. 2011,

J,Horner NK, Potter JD. UDP ACCEPTED MANUSCRIPT , 21(3):173 25(6):1209 ACCEPTED MANUSCRIPT Comment citer cedocument: 27(4):188. - - transferase andtransferase cytochrome P Oldach K,Brandt M, Wittern KP, Pollet D. Comparison ofactivities

. . 1990, 3:180 endocrinol.2012, - 85. - 14.

n EL,nCanlet Debrauwer L,C, Cravedi JP, MavonA, Zalko D.

ms catalyticymes:properties and rational drugdesign. Drug Metab doi: 10.1159/000351683.

doi: 10.1016/j.tiv.2011.03.012. nmodel. Toxicol InVitro. 2010b - Benzo(a)pyrenein a pig ear skinmodel: vivo ex andin vitro - – 3 doi: 10.1016/j.toxlet.2010.08.001.33. acolAppl Skin Physiol. 2002,15( 5.

onA, Duplan H,Zalko D. Effect of skinmetabolism on 4(1):58 Flavin - glucuronosyltransferase (UGT1A1*28andUGT1A6*2) o. Skin Pharmacolo. Physiol. 2007, , ighcapacity ofhuman skinforN 79(5):1897 - activatedgenotoxicity screening. DNA damage - - eedn ,dependent N 4 doi: 10.4161/derm.19201.64. - Xenobiotic metabolism geneexpression in the 450IA1 inculturedkeratinocytes and

- 0. Erratumin: AnalChem. 906. 2008 erations.Toxicol Sci.2013,

- termskin model. , - 24(5):1426 hcns skinmodels.thickness 1

) JB,Rusling JF. :59

- 67. 17beta f humanf vitroin skin - acetylation - 20(4):168

34. - acetyltransferases

Skin doi: - , 133(2):209

- 74.

,

-

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied hydroxy Meinl W 2008 TC7cells varying at stages of differentiation androle benzo[a]pyrenein metabolism. DrugMetab Dispos. Meinl W,Ebert B, Glatt H,Lampen A. 10.1016/j.taap.2015.05.016. ofaromatic amine hairdyes vitro. inToxicol ApplPharmacol. 2015 Extrapolation of systemicbioavailability assessing skinabsorption and epidermal and hepatic metabolism ManwaringJ, fromEpiskin. Steroid J Biochem Biol.Mol 2009, 116(3 metabolizing enzymesin humanskin andthe reconstructed skinmodels Episkin andfull thicknessmodel Luu Toxicol. 2004 LockleyDJ, Howes D 368(4):858 testosterone metabolism inhumanfollicular dermalpapilla cells.B LiuYa S, 2,3,7,8 receptor inhuman cells: CYP1A2 expression in theLS180 colon carcinoma cell line after treatment with LiW Pharmacogenetics. 1999, polymorphisms inCaucasians andAsians: relationships toserum bilirubin concent SaekiM 10.1007/s00 pig, guineapig, man,and inhuman skin models. Arch Toxicol. 2014Dec;88(12):2135 Oesch F,Fabian GuthE, K,Landsiedel R. Xenobiotic skinofman, rat, andpig. DrugMetab Rev. Oesch F,Fabian Oesch E, AdoptedApril,13 2004. OECD, assessment (No. 28) 5March, 2004. OECD, 2004a.Guidance documentforthe conductof skin absorption studies. OECD series ontesting and preparations.Drug Metab Dispos.1990 Myers SR,Flesher JW. Metabolism of Interact. 1989, Myers SR,Blake JW, Flesher JW. Metabolism of3 expressed - 36(2):276, TheV, Duche D,Ferraris C,Meunier JR, LeclaireJ, Labrie , Harper PA, Tang BK,Okey AB. Regulation ofcy - tetrachlorodibenzo 2004b. OECD guideline forthetesting of chemicals. Skin absorption: invitro method. 428. - , Saito, Y, Nagano M, Teshima R, Ozawa S, Sawada J. mRNA expressionof multiple cytochrome 3 ,Tsoi C,Swedmark TibbsS,ZE, Falany CN, Glatt H. Highly selective bioactivation of1 mauchi H. Different patterns of5alpha Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040 - inSalmonella typhimurium. Mutagenesis. 2013 - methylcholanthrene tomutagens individualby human and other mammalian sulphotransferases 64. 204

, RotheH, Obringer C,Foltz DJ, Baker TR, Troutman JA, Hewitt NJ, Goebel C. 71(4):393 78:617 - - doi: 10.1016/j.bbrc.2008.01.130. 014 3 doi: 10.1124/dmd.107.01803683.

- ACCEPTED MANUSCRIPT 1382 – ,Williams FM.Percutaneous penetration andmetabolism of2 28.

- ACCEPTED MANUSCRIPT 401. - - 9(3):341 Comment citer cedocument:

- p 8. BartlomowiczB, Werner C,Landsiedel R. Drug

-

ixn or 3 dioxin

- 9.

the carcinogen3 , Sulfotransferase formsexpressed inhuman intestinal Caco

- 18(5):664 methylcholanthrene. BiochemPharmacol. 1998

2007

- , reductase reductase expression,cellular distribution, and 39(4):659

- 9. - methylcholanthrene in rat liver cy -

eaoiig enzymesmetabolizingin theskin ofrat, mouse, - - tochromeP450 enzymes by arylhydrocarbon methylcholanthrene inhumanbone marrow 5):178 , - 28(5):609 8 doi: 10.1080/0360253070169036698. - 6 doi: 10.1016/j.jsbmb.2009.05.011.86.

F. ExpressionF. profiles ofphases 1 and 2 iochemBiophys Res Commun. 2008, , - 287(2): 19. doi: 10.1093/mutage/get039. - metabolizing enzymes the in 139 - 48. rations.

- doi: butoxyethanol. Arch - tosol.Chem Biol 0 doi:90. ,

56(5):599 - -

- 2 and and2 612.

-

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied ofcytochrome P450 mRNA levelsinhuman skin. AnalBiochem. 2003316(1):103 YengiLG, Xiang Q,Pan J,Scatina Kao J,J,Ball SE,Fruncillo R,Ferron G, RolandWolf Quantitation C. 7(7):e41721. doi: 10.1371/journal.pone.0041721. metabolism pathways inhuman skin andhuman skin models by prot vanEijl S,Z,Zhu Cupitt J,Gierula M, Götz FritscheC, E, Edwards RJ. Elucidation ofxenobiotic 70. doeshuman skin maintain stratum itscorneum lipid properties vitro?in Exp ThakoersingVS,Danso MO,Mulder A,Gooris G, ElGhalbzouri A,Bouwstra JA. doi: 10.1124/dmd.108.024794. Svensson CK. Biotransformation ofdrugs in human sk - Stoming TA,Bronstein W, Bresnick E. 32(1):140 glucuronosyltransferaseinduction usingprimary human hepat Soars MG,Petullo DM, EcksteinJA, Kasper SC,Wrighton SA. P4503A, 2Aand 2C enzymeshavesimilar propertieshumanto analogs. BMCPharmacol 2001, SoucekP, ZuberR, Anzenbacherova AnzenbacherE, P, Guengerich FP. Minipig cytochr byrecombinant humancytochrome P450 1B1. ShimadaT p450isozymes infourtypes of cultured skincells. Int Arch Allergy Immunol. 2002 27(5) nativehuman skin, fourin vitro skin test systems and a liversyste WiegandC,Hewitt NJ,Merk ReisingerHF, K. Dermalxenobiotic metabolism: 1991, 290:160 hydroxylasespecificities ofeleven cDNA WaxmanDJ, LapensonD

a a reinvestigation. BiochemBiophys Res Commun. doi: 10.1111/exd.12031. :263 - - Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040 5 doi: 10.1159/00035827275. 8. ,Gillam EM,Sutter TR, Strickland PT, Guengerich FP,Yamazaki H.. Oxidation ofxenobiotics

– 166.

ACCEPTED MANUSCRIPT ACCEPTED MANUSCRIPT Comment citer cedocument: P, Aoyama GelboinT, HV, GonzalezKorzekwa FJ, K. Steroid hormone

.

Themetabolism of3 - expressed human cytochrome P450s. Arch BiochemBiophys. Dr

ug MetabDispos. 1977, in. DrugMetab Dispos.2009 Feb;37(2):247 21 - 79(2):461 methylcholanthrene rat byliver microsomes ocytes. Drug An assessment of udp m.Skin PharmacolPhysiol. 2014, 1997, eomicprofiling. PLoS One. 2012, - 9.

25(5):617 Metab Dispos.2004, Dermatol. 2012

a a comparisonbetween - Nature Nature versus nurture: , 10. - 22. 127(4):333 -

ome ,

21(11):865 -

6. 1:

- 11. 53.

-

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied inred indic bymore than 2 comparedto control. The lines represent unity ± 2  XMEgeneprofile ofORS Figure 3 – assignedwere accordingthe as Ct follows: noexpression ( systemically withDMSO (1% v/vfinal concentrati Figure 2. metabolism assays. Figure 1 Figure legends 0.001. statistically significant fromp <0.05. notNS significant p>0.05, *p<0.05, ** p < 0.01 and ***p< a expressed as the mean ±standard error inpmole application) orafter application of100 nmoles/cm 3 Figure 4 statistically significant. ndstatistical analysis was calculated using a two - M 32;Good expression ( MC. EROD activity was measured 24 haft  NF inNF the medium orafter topicalapplication of100 nmoles

. . An overviewof thegeneration of theORS

Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040 EROD EROD activityin ORS Effect of systemic application of B t the change was statisticallyate significant andXMEs in black indicate thechange was not asal gene expression ofI Phase (A) and (B) II XME - od andfold symbols in bluerepresent genes that bydecreased were more than2

ACCEPTED MANUSCRIPT

ACCEPTED MANUSCRIPT Comment citer cedocument: - RHE models.Thegene profilesmeasured were 24 hafter treatment with25 ) =Ct 25 - RHE models treated systemically (A) ortopically with (B) –

30 and30strong expression (  NF (A) (A) NF and topical application of3 er er treatment with25µM s /mg/min. Data were /mg/min. Data obtainedwere from 3 donors intriplicate - - 2 tailedStudent’s t od Symbolsfold. inred represent genes that inducedwere on)in themedium. Thedifferent levels ofexpression

of - RHE modeland theincubation conditions for  NF orNF 3 ) =Ct) >weak 32; expression ( - target genes MC (topical application). EROD was /cm - test. test. Differences consideredwere 2  ) 3 NF orNF 10 µM 3 =Ct <25. - MC. Valuesare fold increase inORS - MC(B) tothe I Phase

- RHE tissues treated - MC(systemic -  fold.XMEs NF orNF ) =30Ct

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied andHC correlationcoefficients (R glucuronide ( metabolites(n=3 donors, mean± after topicalapplication of7 Figure 6. significant. inred indicate the change was significant bymore than 2 comparedto control. The lines represent unity ± 2  XMEgeneprofile ofORS Figure 5. testosterone); XIPeak =androstenedione; XIIPeak =epiandrosterone; XIII’Peak = 5α hydroxytestosterone(standard notavailable); X Peak =17β 3,17 hydroxy hydroxytestosterone; IV Peak = hydroxy follows: I Peak = 7 hydroxylated ( metabolites(n=3 after topicalapplication oftestosterone. The doses applied 1, were 2 and 4 nmoles/cm Figure 7. M -  dione); Peak VIII Peak =dione); hydroxy NF inNF the medium orafter topicalapplication of100 nmoles/cm - - s 4 ulfate,respectively.

Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040 - Dose Dose Effect of systemic applicatio

androste - - ,Pea - od andfold symbols in bluerepresent genes that bydecreased were more than2 dependency (A) and radiochromatogram(B) ofmetabolite production inthe medium h24 dependency (A) and radiochromatogram(B) ) metabolites.) The metabolites oftestost

donors,mean±  ne ACCEPTED MANUSCRIPT - k I inradiochromatogram) andHC

- hydroxytestosterone; II Peak =6 3,17 ACCEPTED MANUSCRIPT Comment citer cedocument: 2 - - ) forthegoodness of fit 0.94,were 0.98and 0.86 for theHC, HC RHE models.Thegene profilesmeasured were 24 hafter treatment with25 dione; Peak VI Peak =4 dione; - EC. Thedoses applied 1, were 2 and 4 nmoles/cm - 4 SEM SEM - androstene of n

) shownin (A) are combined valuesfor reduced( ) formed HC were ( - 4 and XMEsin blackindicate thechange was notstatistically - androstene  NF (A) (A) NF and topical application of3 - - androsten 3,17 - od Symbolsfold. inred represent genes that inducedwere - in (standard available);notdione IXPeak = -  3,17 - - sulfate ( hydroxytestosterone; III Peak =16 erone shownin theradiochromatogram (B) are as - 16α - - in (standardnot available); dione V’ Peak = hydroxyandrosta

ofmetabolite production inthe medium h24 , - ol peak visiblenot inchromatogram), HC - 3,17 , Peak ,II Peak inradiochromatogram). The 2 - in (16 dione 3 - MC. Valuesare fold incr - 4,6 2 . The . amountof -  - MC(B) tothe II Phase die - hydroxy n - 2 3 . Theamount . of - - androstan n ( one -  ) and 4 - - - glucuronide - fold.XMEs androstene  6 - ease ease - 3,17 - - -

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied metabolites(n=3 donors, mean± after topicalapplication of3 Figure 9. 1.00,1.00and 0.61for I, thePhase HC B(a)P III =B(a)P radiochromatogram are as (B) follows: I Peak =B(a)P glucuronidated ( metabolites(n=3 after topicalapplication The of B(a)P. doses applied 1,were 2 and 4 nmoles/cm Figure 8. 0.98for thereduced andhydroxylat dione andXIIIPeak =androsterone. Thecorrelation coefficients (R were 0.80,were 0.97and 0.98forthe I, Phase HC =2 hydroxy 3MC sulfate hydroxy hydroxy radiochromatogram are as (B) follows: I Peaks and =3 II glucuronidated ( - hydroxy -

- sulfate; XPeak =3 catechol; VII Peak = B(a)P ; Peak VIII Peak ; = 2 - - - 3MC 3MC 3MC

Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040 Dose Dose - - diol 3MC; Peak XVII Peak 3MC;=hydroxy - - - IV Peak = sulfate;2

- - - glucuronide; Peak glucuronide; XIV Peak =2 glucuronide; Peak glucuronide; IV Peak = B(a)P dependency (A) and radiochromatogram(B) ofmetabolite production inthe medium h24 dependency (A) and radiochromatogram(B) ofmetabolite production inthe medium h24

ACCEPTED) andsulfated ( MANUSCRIPT) andsulfated ( donors,mean± - hydroxy ACCEPTED MANUSCRIPT Comment citer cedocument: - hydroxy - MC. Thedoses applied 1, were 2 and 4 nmoles/cm - 3MC - - VI =2 SEM hydroxide.The correlation coefficients(R SEM - 3MC metabolites, respectively.ed hydroxy metabolites.) Th metabolites.) The metabolites ofshownB(a)P inthe -

glucuronide +2

- ) shownin (A) are combined values glucuronide andHC ) shownin (A) are combined valuesfor I Phase ( - - -

hydroxy 3MC. The 3MC.correlation coefficients(R sulfate; Peaks sulfate;XIPeaks and =XII 3 - hydroxy - - 3MC glucuronide andHC - O - - - glucuronide; 3MC - glucuronide; V Peak = 3 3MC e e metabolitesof3 - - diol hydroxy -

- glucuronide; VII Peak = 2 hydroxy ; Peak XVPeak ; - sulfate; II Peak =B(a)P - sulfate,respectively. Peak V Peak = B(a)P -

3MC - - sulfat 3MC 2 ) forthe goodnessoffit 0.97were and

- =Hydroxy

- sulfate; IXPeak =2 - MCshown in the e,respectively. - hydroxy glucuronide; Peak glucuronide;III Peak =3

- for I Phase ( 2 hydroxy ) for thegoodness of fitwere - 2 2 2 O . Theamount . of ) for thegoodness of fit - - . Theamount . of 3MC XIII Peak = 3MC;

- sulfate; VI Peak = - - hydroxy 7,8,9,10 - 3MC - sulfate ), - ), - uft and2 sulfate - hydroxy - tetrol; Peak Peak tetrol; 3MC

; XVIPeak - -

- -

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied hydroxylated sulfated( formed from B(a)P shownin (B) are combined valuesfor reduced ( HC were ( produced±SD, n=3(triplicate ORS nmole/cm topical application of (A) 7 Supplementary Figure 1 2 Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040

(0.33cm ) metabolites.) The metabolites 3 from ), HC ( ),glucuronidated ( - glucuronide ( (

ACCEPTED MANUSCRIPT 2

discarea). ) ACCEPTED MANUSCRIPT shownin (C) Comment citer cedocument: .

- Time EC, (B) EC,testosterone, (C)

- The amountof parent andmetabolites was expressed as the mean nmole dependency ofmetabolite production inthe medium over h72after ) and) HC - RHE models). ) andsulfated (

are are comb

- sulfate ( inedvalues forI Phase - ) and) hydroxylated ( MC

Themetabolites formedfrom 7  ( ) metabolites.) ). ). The metabolites formedfrom testosterone ▲ a)P ) shownin (D) are combined valuesfor

and (D)

3 - MC

( ) metabolites.Themetabolites . . The dose applied was 1 ), glucuronidated ( - EC ( ) shownin (A) ) and) ( )

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied Fig.1

Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040 ACCEPTED MANUSCRIPT ACCEPTED MANUSCRIPT Comment citer cedocument:

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied Fig.2

Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040 ACCEPTED MANUSCRIPT ACCEPTED MANUSCRIPT Comment citer cedocument:

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied Fig.3

Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040 ACCEPTED MANUSCRIPT ACCEPTED MANUSCRIPT Comment citer cedocument:

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied Fig.4

Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040 ACCEPTED MANUSCRIPT ACCEPTED MANUSCRIPT Comment citer cedocument:

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied Fig.5

Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040 ACCEPTED MANUSCRIPT ACCEPTED MANUSCRIPT Comment citer cedocument:

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied Fig.6

Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040 ACCEPTED MANUSCRIPT ACCEPTED MANUSCRIPT Comment citer cedocument:

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied Fig.7

Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040 ACCEPTED MANUSCRIPT ACCEPTED MANUSCRIPT Comment citer cedocument:

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied Fig.8

Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040 ACCEPTED MANUSCRIPT ACCEPTED MANUSCRIPT Comment citer cedocument:

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied Fig.9

Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040 ACCEPTED MANUSCRIPT ACCEPTED MANUSCRIPT Comment citer cedocument:

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied Hu et Huet al.(2010). =absent; NR =not reported; = ND Notdetermined; =reported a by Luu (2009),al. Hu(2010),et Wiegand (2014) al. and et Svensso comparisonnative humanto skin. Genes shown to presentbe in native human fromskin are Luu 1. Table CYP4F11 CYP4F3 CYP4F2 CYP4B1 CYP3A5 CYP2S1 CYP2R1 CYP2C18 CYP1B1 CYP1A1 1 Phase UCHL1 ALDH7A1/9A1 ALDH4A1/5A1/6A1 ALDH3B1/2 ALDH3A2 ALDH2 ALDH1A1/2/3 ADH5 MAOA FMO4 5 CYP27B1 CYP26B1 CYP19A1 CYP7B1  - Reductase

XME

SummaryidentifiedXMEscontrolofORS in Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040

ACCEPTED MANUSCRIPT

ACCEPTED MANUSCRIPT Comment citer cedocument: Native human Native skin NR NR NR + + + + + + + + + + + + + + + + + + + - - -

Control ORS Control RHE model RHE ND

+ + + + + + + + + + + + + + + + + + + + + + - - - levelthegenein present at modelsare thatRHE

(2009)al. (2014).et and Oesch present, +=

- - The et Theet al.(2009); b =reported by - The et Theet al. -

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied SULT1A3 SULT1A1 MGST1/2/3 GSTT2 GSTT1 GSTP1 GSTO1/2 GSTM2/3/4 GSTK1 GST1A4 2 Phase XHD DPYD HSD17B10 PTGS1/2 NQO1/2 CES1/2/5A EPX1 ESD UCHL3 TPMT NNMT GAMT COMT AS3MT SAT1 NAA20 NAT1 UGT1A6/8/10 UGT1A3/5/7 TST SULT2B1 SULT1E1 SULT1B1 SULT1A5/8

Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040

ACCEPTED MANUSCRIPT ACCEPTED MANUSCRIPT Comment citer cedocument: + NR NR NR NR NR NR NR NR NR NR NR NR NR NR NR NR NR a + + + + + + + + + + + + + + + + ,

- b

+ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + - -

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied

CCBL1 DDOST

Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040 ACCEPTED MANUSCRIPT ACCEPTED MANUSCRIPT Comment citer cedocument: NR NR

+ +

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied Graphical abstract Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040

ACCEPTED MANUSCRIPT ACCEPTED MANUSCRIPT Comment citer cedocument:

Version postprint Bessou-Touya, Zalko, D.,Duplan(2017). Characterization ofxenobiotic metabolizing enzymesof Bacqueville, D.(Auteur decorrespondance), Jacques, Duprat,Jamin, E.L.,Guiraud, Perdu,E., a reconstructed human epidermal model fromadult hair follicles. Toxicology and Applied HIGHLIGHTS     

ORS enzymes are inXenobiotic metabolizing inducibleORS ORS in functional are enzymes metabolizing Xenobiotic ORS in expressed are enzymes metabolizing xenobiotic II and I Phase ORS rtpcl applicationor topical derived from follicles adult hair Pharmacology, 329, 190-201. DOI : 10.1016/j.taap.2017.05.040 - - RHE andforreliable toolevaluationxenobioticsathe convenientis of RHE metabolismstudying skin modelfora suitableis ACCEPTED MANUSCRIPT ACCEPTED MANUSCRIPT Comment citer cedocument:

- RHE - RHE following either systemic systemic either following RHE

- H model RHE