DMD DiLab, Inc.

Home Help [Feedback] [For Subscribers] [Archive] [Search] [Contents]
 QUICK SEARCH:   [advanced]


     


Drug Metabolism and Disposition Fast Forward
First published on September 8, 2006; DOI: 10.1124/dmd.106.011262


0090-9556/06/3412-2083-2090$20.00
DMD 34:2083-2090, 2006

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
dmd.106.011262v1
34/12/2083    most recent
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Kienhuis, A. S.
Right arrow Articles by Stierum, R. H.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Kienhuis, A. S.
Right arrow Articles by Stierum, R. H.

Comparison of Coumarin-Induced Toxicity between Sandwich-Cultured Primary Rat Hepatocytes and Rats in Vivo: A Toxicogenomics Approach

Anne S. Kienhuis, Heleen M. Wortelboer, Jean-Christophe Hoflack, Edwin J. Moonen, Jos C.S. Kleinjans, Ben van Ommen, Joost H. M. van Delft, and Rob H. Stierum

Business Unit Biosciences, TNO Quality of Life (Netherlands Organization for Applied Scientific Research), Zeist, The Netherlands (A.S.K., H.M.W., B.v.O., R.H.S.); Department of Health Risk Analysis and Toxicology, University of Maastricht, Maastricht, The Netherlands (A.S.K., E.J.M., J.C.S.K., J.H.M.v.D.); Independent author (J.C.H.)

(Received June 1, 2006; accepted September 7, 2006)


    Abstract
 Top
 Abstract
 Materials and Methods
 Results
 Discussion
 References
 
Sandwich-cultured primary rat hepatocytes are often used as an in vitro model in toxicology and pharmacology. However, loss of liver-specific functions, in particular, the decline of cytochrome P450 (P450) enzyme activity, limits the value of this model for prediction of in vivo toxicity. In this study, we investigated whether a hepatic in vitro system with improved metabolic competence enhances the predictability for coumarin-induced in vivo toxicity by using a toxicogenomics approach. Therefore, primary rat hepatocytes were cultured in sandwich configuration in medium containing a mixture of low concentrations of P450 inducers, phenobarbital, dexamethasone, and ß-naphthoflavone. The toxicogenomics approach used enabled comparison of similar mechanistic end-points at the molecular level between in vitro and in vivo conditions, namely, compound-induced changes in multiple genes and signaling pathways. Toxicant-induced cytotoxic effects and gene expression profiles observed in hepatocytes cultured in modified medium and hepatocytes cultured in standard medium (without inducers) were compared with results from a rat in vivo study. Coumarin was used as a model compound because its toxicity depends on bioactivation by P450 enzymes. Metabolism of coumarin toward active metabolites, coumarin-induced cytotoxicity, and gene expression modulation were more pronounced in hepatocytes cultured in modified medium compared with hepatocytes cultured in standard medium. In addition, more genes and biological pathways were similarly affected by coumarin in hepatocytes cultured in modified medium and in vivo. In conclusion, these experiments showed that for coumarin-induced toxicity, sandwich-cultured hepatocytes maintained in modified medium better represent the situation in vivo compared with hepatocytes cultured in standard medium.


To assess possible hepatotoxicity, conventional studies rely on the use of animal model systems to examine tissue toxin levels, changes in serum levels of hepatic enzymes, and histopathological changes (Nuwaysir et al., 1999Go; Waring and Ulrich, 2000Go). Simple, well established in vitro assays, such as primary hepatocytes, precision-cut liver slices, and hepatic cell lines, are increasingly in demand for identifying potential hepatotoxicity in early stages of investigative toxicology and for decreasing attrition rates of drugs during lead optimization (Dambach et al., 2005Go). However, extrapolation of in vitro results to the in vivo situation remains a scientific challenge (Guillouzo, 1998Go).

Toxicogenomics, the application of the genomics technologies in toxicology, would be particularly useful in the extrapolation from in vitro experiments to the in vivo situation. Extrapolations can be made at the molecular level, comparing similar mechanistic endpoints, namely compound-induced changes in multiple genes and signaling pathways (Hamadeh et al., 2002Go; Stierum et al., 2005Go). Several toxicogenomics-based studies have already been performed comparing rat hepatic in vitro models with the situation in vivo (Waring et al., 2001Go; Boess et al., 2003Go; Jessen et al., 2003Go). These studies concluded that, to date, no toxicogenomics-based in vitro system allowed for prediction of hepatotoxic responses in vivo. The main limitation of hepatic in vitro assays used in these toxicogenomics-based studies is the loss of liver-specific functions, in particular, cytochrome P450 (P450) monooxygenase activities (Balls et al., 2002Go; Boess et al., 2003Go). Extrapolation from these in vitro models to the in vivo situation is therefore hampered when examining compounds for which toxicity depends on bioactivation by the P450 enzyme system.

In the aim to develop an in vitro toxicogenomics-based system, the relevance of a hepatocyte sandwich culture with improved metabolic competence toward prediction of in vivo toxicity was assessed. Sandwich-cultured hepatocytes were used because, compared with other hepatocyte-based in vitro models, hepatocyte longevity is increased and a polarized cell and membrane architecture resembling that in vivo is maintained for several weeks (Dunn et al., 1989Go; LeCluyse et al., 1996Go). Furthermore, rat hepatocytes maintained in sandwich configuration display a more optimal P450 inducibility (LeCluyse et al., 1999Go; Richert et al., 2002Go). In a separate study performed by this laboratory, it was shown that in sandwich-cultured hepatocytes maintained in modified medium, enriched with low concentrations of the known P450 inducers phenobarbital (PB), dexamethasone (DEX), and ß-naphthoflavone (ß-NF), the metabolic competence, as reflected by retention of P450 enzyme activities and gene expression levels of several phase I and phase II enzymes, was enhanced (A. S. Kienhuis, H. M. Wortelboer, W. J. Maas, M. van Herwijnen, J. C. S. Kleinjans, J. H. M. Delft, and R. H. Stierum, manuscript submitted for publication). In the present study, the effect of coumarin (cis-o-coumaric acid lactone), a compound for which toxicity depends on metabolic activation by the P450 enzyme system (Born et al., 2002Go; Lake et al., 2002Go), on toxicity and gene expression profiles was studied in sandwich-cultured hepatocytes maintained in standard medium, the standard model (without inducers), and in sandwich-cultured hepatocytes maintained in modified medium, the modified model (containing enzyme inducers). Coumarin-induced cytotoxicity data and gene expression profiles generated in both in vitro models were compared with toxicity and gene expression profiles in rat liver in the in vivo study in which rats were exposed to coumarin, since in vivo verification of in vitro toxicogenomics data has been proposed to be a necessity to judge the value of the in vitro model for in vivo toxicity prediction (Jaeschke, 2003Go). The expected differences of bioactivation of coumarin in the standard model versus the modified model were verified by measurement of coumarin and one of its metabolites in medium.


    Materials and Methods
 Top
 Abstract
 Materials and Methods
 Results
 Discussion
 References
 
Chemicals. Collagenase type B, Roche reagent kits, primer p(DT)15, and EasyHyb were obtained from Roche (Mannheim, Germany). Dulbecco's modified Eagle's medium (DMEM), fetal calf serum, gentamicin, phosphate-buffered saline, TRIzol, yeast tRNA, and human Cot-1 DNA were obtained from Invitrogen (Breda, The Netherlands). Insulin, glucagon, hydrocortisone, PB, DEX, ß-NF, 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT), and amino-allyl dUTP were purchased from Sigma-Aldrich (Zwijndrecht, The Netherlands). The RNeasy mini kit, the RNase-free DNase kit, and the QIAquick PCR purification kit were obtained from QIAGEN (Westburg B.V., Leusden, The Netherlands). Cy3 and Cy5 monofunctional reactive dyes, and poly(dA) · poly(dT) were purchased from Amersham Biosciences (Roosendaal, The Netherlands). Coumarin, CAS-no 91-64-5, purity by HPLC minimum 99% according to the manufacturer, was obtained from Sigma-Aldrich. All other chemicals were of analytical grade.

Animals. Male Wistar rats (Crl:(WI)WU BR), 9 to 12 weeks of age, 180 to 250 g, were obtained from Charles River GmbH (Sulzfeld, Germany). During the acclimatization period and until sacrifice, animals were housed individually in Macrolon cages with wire tops and sawdust bedding at 22°C and 50 to 60% humidity. The light cycle was 12-h light/12-h dark. Feed and tap water were available ad libitum.

Animal Treatment in Vivo. Wistar rats, housed under conditions as described above, were injected i.p. with 17.5 mg/kg b.wt. (low dose), 75 mg/kg b.wt. (mid dose), and 200 mg/kg b.wt. (high dose) coumarin dissolved in corn oil. Doses were defined in a range finding study. As a solvent control, only corn oil was injected (vehicle-treated control). In each dose group and vehicle-treated group, n = 5 rats were included. Injection volume for each treatment was 10 ml/kg b.wt. Body weight was recorded on day 0 and at 24 h, just before sacrifice. Rats were anesthetized by inhalation of CO2/O2. Twenty-four hours after i.p. injection, animals were sacrificed by bleeding through the aorta abdominalis; the blood was collected in a heparin tube, from which plasma was isolated for clinical chemistry. Thereafter, livers were immediately dissected, frozen in liquid nitrogen, and stored at –80°C until further processing. A section of the liver was kept aside in formalin for pathological examination.

Clinical Chemistry. Serum alanine aminotransferase (ALAT) activity, serum aspartate aminotransferase (ASAT) activity, lactate dehydrogenase (LDH) leakage, serum alkaline phosphatase activity, glucose, cholesterol, {gamma}-glutamyl transferase (GGT) activity, and phospholipid levels were analyzed on a Hitachi 911 centrifugal analyzer using Roche reagent kits. Differences between plasma levels in treated and nontreated animals were defined as statistically significant at a p value below 0.01, determined by one-way analysis of variance followed by Dunnett's test.

Culture of Rat Hepatocytes. Male Wistar rats similar to those used in the in vivo study and housed under identical conditions were used for hepatocyte isolation. Hepatocytes were isolated from livers from three individual rats according to a two-step collagenase perfusion technique as described by Seglen (Seglen, 1976Go; Paine et al., 1979Go) with minor modifications (Paine et al., 1979Go). Hepatocyte preparations with viability greater than 85% as determined by trypan blue exclusion were used and cultured on collagen gelprecoated six-well plates at a density of 1.25 x 106 cells per well. Sandwich cultures were essentially prepared according to the method of Beken et al. (2004Go). Hepatocytes were allowed to attach for 4 h in DMEM supplemented with 10% fetal calf serum, insulin (0.5 U/ml), glucagon (7 ng/ml), and gentamicin (50 µg/ml). After attachment, dead cells were removed by washing, and the upper collagen layer was applied. Cells were kept in standard medium consisting of DMEM containing 25 mM HEPES and 4.5 g/l D-glucose supplemented with insulin (0.5 U/ml), glucagon (7 ng/ml), hydrocortisone (7.5 µg/ml), and gentamicin (50 µg/ml). In the modified model, standard culture medium was modified by supplementation of an inducer mix that consisted of 1 mM PB, 10 µM DEX, and 5 µM ß-NF (A. S. Kienhuis, H. M. Wortelboer, W. J. Maas, M. van Herwijnen, J. C. S. Kleinjans, J. H. M. Delft, and R. H. Stierum, manuscript submitted for publication). PB was added as a concentrated stock solution in phosphate-buffered saline. DEX and ß-NF were added as concentrated stock solutions in dimethyl sulfoxide (DMSO). The final concentration of DMSO was equalized in all culture media and did not exceed 0.2% (v/v). Standard and modified media were applied 4 h after seeding the cells. Cultures were incubated at 37°C in a humidified incubator, O2 95% and CO2 5%. Medium was changed on a daily basis during a period of 72 h.

Hepatocyte Treatment and Cytotoxicity Analysis. Hepatocytes cultured in either standard or modified medium for 72 h were exposed to coumarin in a concentration range of 0 to 1000 µM dissolved in DMSO for 24 h. The final concentration of DMSO was 0.5% (v/v). Cytotoxicity was determined using the MTT reduction method (Mosmann, 1983Go). LDH leakage was determined spectrophotometrically in pooled medium obtained from three wells per rat per dose group, on a Hitachi 911 centrifugal analyzer using Roche reagent kits. The final coumarin concentrations for the gene expression study selected were 0 µM (control), 70 µM (low dose), 200 µM (mid dose), and 600 µM (high dose). Doses corresponded to 0%, 10%, 20%, and 50% cytotoxicity as determined by the MTT reduction assay, respectively (see results).

Gas Chromatography to Study Metabolism of Coumarin. Coumarin and the o-hydroxyphenylacetic acid (o-HPAA) metabolite were measured in standard and modified culture media of hepatocytes exposed to coumarin for 24 h using gas chromatography with flame ionization detection as described previously (Born et al., 2000aGo, 2003Go). Samples (0.1 ml) were diluted with water to a final volume of 1 ml. To increase the extraction efficiency of o-HPAA, 0.5 g of sodium chloride was added to each sample, after extraction with ethyl acetate. The extraction efficiencies of coumarin and o-HPAA were 30% and 65%, respectively. Coumarin and o-HPAA were separated and quantitated using a Fisons HRGC 8650 gas chromatograph (Fisons, Loughborough, UK) and a Varian VG-5ms column, 50 m x 0.25 mm (Varian, Inc., Palo Alto, CA). Quantification was accomplished by calculation of peak-area ratios relative to the n-tridecane internal standard, and use of the external standard curves. The same extraction procedure was applied to the external standards.

Total RNA Isolation. RNA was extracted using TRIzol (Invitrogen) according to the manufacturer's protocol. TRIzol was added to frozen liver samples obtained from the in vivo study which subsequently were pulverized with mortar and pestle under liquid nitrogen before extraction. To obtain RNA from sandwich cultures, TRIzol was added on the upper collagen layer, and cells were collected. RNA was purified using the RNeasy mini kit including an additional DNase digestion step. RNA concentration was calculated from the absorbance at 260 nm as measured spectrophotometrically. RNA quality was assessed by agarose gel electrophoresis.

Microarray Design. A local reference design was used for microarray hybridization (Table 1). For the in vivo samples, tester RNA labeled with Cy5 from individual rat livers was combined with Cy3-labeled reference RNA consisting of a pool of RNA extracted from livers obtained from five vehicle-treated control rats. For the in vitro studies, Cy5-labeled tester RNA from one experiment was combined with Cy3-labeled reference RNA from the same experiment; hepatocytes in one experiment were obtained from one rat.


View this table:
[in this window]
[in a new window]

 
TABLE 1 Experimental design

 

Microarray Labeling. RNA samples were indirectly labeled according to the amino-allyl labeling procedure for microarrays from The Institute for Genomic Research (TIGR) (http://www.tigr.org/tdb/microarray/protocolsTIGR.shtml). For each labeling reaction, 25 µg of total RNA was used as the starting amount for reverse transcription of mRNA. In brief, for reverse transcription, mRNA was selected by oligo(dT) priming with primer p(DT)15. The reverse transcription reaction, in which amino-allyl dUTP is incorporated, was conducted for 3 h at 42°C. Nontranscribed RNA was degraded by alkaline hydrolysis in a final concentration of 0.25 M NaOH for 30 min at 37°C. Thereafter, the mixture was neutralized with an equimolar amount of acetic acid. The cDNA was purified using the QIAGEN QIAquick PCR purification kit. Columns were washed with 10 mM sodium-borate in 80% ethanol (pH 8.5). Column-bound cDNA was eluted two times in 30 µl of MilliQ H2O. Samples were labeled with either Cy3 (reference) or Cy5 (tester) monofunctional reactive dyes and afterward cleaned from unincorporated Cy dyes using AutoseqG-50 Sephadex chromatography columns.

Microarray Hybridization. Samples were hybridized on rat oligonucleotide microarrays containing approximately 5800 different oligonucleotide fragments of 70 nucleotides in length (QIAGEN Operon, Westburg B.V., Leusden, The Netherlands), spotted in duplicate at the Frank Holstege group (Utrecht Genomics Laboratory, Utrecht, The Netherlands) as described previously (van de Peppel et al., 2003Go; Heijne et al., 2005Go). In brief, Cy3- and Cy5-labeled samples were combined according to the experimental design. To avoid nonspecific hybridization, yeast tRNA, poly(dA) · poly(dT), and human Cot-1 DNA were added. Samples were dissolved in 110 µl of EasyHyb hybridization buffer. After cDNA denaturation at 100°C, human Cot-1 DNA was allowed to anneal for half an hour at 42°C. The hybridization mixture was pipetted directly in the center of the hybridization chamber (Corning Life Sciences B.V., Schiphol-Rijk, The Netherlands). Slides were prehybridized according to the TIGR protocol and placed carefully on top. To keep the chamber moisturized during hybridization, a precut paper drenched in MilliQ H2O was put on top. Samples were hybridized overnight at 42°C in a water bath. Slides were washed with sodium chloride/sodium citrate buffers decreasing in stringency. Microarrays were scanned using a Packard Scanarray confocal laser scanner (PerkinElmer, Boston, MA). Resulting TIFF images were loaded into Imagene 5.0 (Biodiscovery Inc., El Segundo, CA) and saved to further process and analyze the data.

Microarray Quality Criteria. Criteria for microarrays to be further analyzed considered the homogeneity of the spot signal intensities, the effect of bleaching of the fluorescence, the number of manually flagged (excluded) spots, the spatial distribution of the signals over the slide surface, the balance between Cy3 and Cy5 signal intensity, the number of saturated spots, and the quality of the slide with respect to other slides of the experiment, as described by Heijne et al. (2005Go).

Microarray Data Analysis. Flagged spots and controls were excluded from further analyses. Only spots on qualified microarrays with intensities higher than 1.5 times the intensity of their local background were included in data analysis. Ratios of the background-corrected intensities of tester over reference were calculated for each slide. To account for technical variations introduced during labeling or hybridization, data were normalized using the LOWESS algorithm (Yang et al., 2002Go). Normalized ratios were log-transformed with base two in SAS Enterprise guide V2 (SAS Institute Inc., Cary, NC). The resulting data set was loaded into Excel (Microsoft Corporation, Redmond, WA). Replicate genes per array were averaged. For a gene to be included in the analysis for both in vivo and in vitro studies, a maximum of 30% missing values per gene was accepted. Two-sample t tests (with random variance model) were performed using BRB ArrayTools Version 3.3.1, developed by Dr. R. Simon and A. Peng Lam (Simon and Peng, 2002Go). To reduce false positives discovery, the nominal significance level of every univariate test was set at 0.001. To study whether similarities between biological processes could be observed between the in vitro models and in vivo, pathway analysis was performed using genes significantly expressed as determined by BRB ArrayTools, and Genelists were uploaded on the website of the Database for Annotation, Visualization and Integrated Discovery (DAVID) (http://apps1.niaid.nih.gov/david/). Pathway analyses were performed using the Kyoto Encyclopedia of Genes and Genomes (KEGG) (http://www.genome.ad.jp/kegg). A pathway was considered "triggered" when at least two genes within the pathway were significantly modulated by coumarin as assessed by means of Fisher's exact test. Therefore, for each pathway and per system (in vitro and in vivo), the total number of genes that gave a signal on the microarray was divided into four groups: the number of genes not significantly modulated by coumarin and representing the pathway, those not significantly modulated by coumarin and not representing the pathway, those significantly modulated by coumarin and representing the pathway, and finally, those significantly modulated by coumarin and not representing the pathway. Pathways were considered significant at a p value below 0.05. Furthermore, principal component analysis (PCA) was performed using Matlab software Version 6.5 (The MathWorks Inc., Natick, MA).


Figure 1
View larger version (34K):
[in this window]
[in a new window]

 
FIG. 1. Cytotoxicity of coumarin in the standard model (A) and the modified model (B). MTT reduction is expressed as a percentage of the MTT determination of control hepatocyte cultures incubated with vehicle only. LDH retention is expressed as an inverted percentage of the LDH leakage in Triton X-treated hepatocyte cultures. Data are means ± S.D. (n = 3).

 

    Results
 Top
 Abstract
 Materials and Methods
 Results
 Discussion
 References
 
In Vivo Toxicity. Hepatotoxicity in vivo was determined in rats treated with coumarin up to 200 mg/kg b.wt. Clinical chemistry parameters showed that coumarin significantly increased the plasma levels of cholesterol at a dose level of 75 mg/kg b.wt., and ASAT, ALAT, and GGT at a dose of 200 mg/kg b.wt.. The phospholipid level measured in plasma of rats exposed to 200 mg/kg b.wt. coumarin was significantly decreased, as presented in Table 2 (data represent mean ± S.D.). Histopathological findings in rat liver showed that at necropsy, the liver of one animal treated with 75 mg/kg b.wt. coumarin was focally discolored, probably by focal congestion. Upon microscopy, only animals treated with 200 mg/kg showed slight to severe single cell necrosis and, minimally, centrilobular and mononuclear cell infiltrate.


View this table:
[in this window]
[in a new window]

 
TABLE 2 Clinical chemistry parameters from n = 5 coumarin-treated rats per treatment group in the in vivo study

 


Figure 2
View larger version (25K):
[in this window]
[in a new window]

 
FIG. 2. Circles in the de Venn diagram represent the number of genes significantly changed in vivo, in the standard model (ST), and in the modified model (MOD). Overlaps contain the number of significantly modulated genes similar between systems.

 
In Vitro Toxicity. Cytotoxicity of coumarin was determined in the sandwich-cultured hepatocyte model containing either standard medium or modified medium (Fig. 1). No cytotoxicity of coumarin (up to 1 mM) was detected in the standard system. However, in the modified system, coumarin appeared to be cytotoxic in a dose-response manner. Based on these data, for the gene expression study, hepatocytes were exposed to coumarin at 0 µM (control), 70 µM (low dose), 200 µM (mid dose), and 600 µM (high dose) for 24 h. In the modified system, the low, mid, and high dose resulted in 100 to 90%, 90 to 80%, and <50% viability, respectively. In the standard system, the viability was 100% at all concentrations of coumarin.

Coumarin Metabolism. To verify that the improved metabolic competence of the modified model indeed affects the metabolism of coumarin and thus is responsible for differences in cytotoxicity of coumarin in standard culture medium and modified culture medium, both coumarin and the major coumarin metabolite reported to be detected in rat urine (Lake, 1999Go; Born et al., 2000aGo,bGo) were measured in both standard and modified culture media 24 h after exposure. A dose-dependent increase in the formation of o-HPAA was exclusively measured in the modified culture medium, but no o-HPAA was present in the standard culture medium. Furthermore, for all dose groups, coumarin concentration was significantly lower in the modified medium compared with the standard medium (the reduction varied from 15% for the high dose group to 59% for the low dose group; data not presented).

Gene Expression Profiling. Analysis of significantly modulated genes. In total, the criteria settings for microarray analysis used in this study resulted in data sets containing gene expression ratios of 1621 genes in vivo, 2536 genes in the standard system, and 2368 genes in the modified system. The large volume of data in the data set was reduced by performing robust statistical analysis that separated genes in which expression was actually changed by coumarin from genes that remained unchanged. Resulting data sets contained, regardless of dose dependence, 321, 13, and 92 genes that were significantly changed by coumarin in vivo, in the standard system, and in the modified system, respectively. Overlap of the changed genes between both in vitro models and in vivo is presented by the Venn diagram in Fig. 2. Twenty-three genes were altered in both the modified model and in vivo, and only three genes in both the standard system and in vivo. One gene was altered in all systems, in vitro and in vivo. These 27 genes form the subset of genes significantly changed by coumarin which are listed in Table 3. Both in vivo and in the in vitro models, most genes showed a dose-dependent down-regulation (for 52–85% of the genes, the Pearson correlation analysis coefficient is less than –0.8; data not shown). When comparing the different systems with each other, the direction of modulation was always comparable between the standard and the modified in vitro models and, in most of the cases, also between the in vivo model and one of the in vitro models.


View this table:
[in this window]
[in a new window]

 
TABLE 3 Genes significantly changed by coumarin in at least two systems: in vivo, in the standard, or in the modified model

Values are log-transformed ratios corrected for the control; significant values (p < 0.001) are marked in bold.

 

Pathway analysis. Pathways triggered by coumarin in at least two systems and significant in at least one are presented in Table 4. In the standard model, the modified model, and in vivo, one, eight, and six pathways, respectively, were significantly triggered by coumarin. Four pathways, i.e., methionine metabolism, fatty acid metabolism, {gamma}-hexachlorocyclohexane degradation, and complement and coagulation cascades, were significant in both the modified system and in vivo. Only one pathway was significantly triggered in the standard system, but not in vivo.


View this table:
[in this window]
[in a new window]

 
TABLE 4 Pathways triggered by coumarin in vivo, in the standard model, or in the modified model

 


Figure 3
View larger version (18K):
[in this window]
[in a new window]

 
FIG. 3. Principal component analysis of the complete data set of expression profiles resulting from in vitro and in vivo hepatotoxicant treatment of hepatocytes and rats. Two principal components were generated and plotted for each system and dose group. Dose groups are distinguished by circles and identified within (in vivo) or outside (ST, the standard model; MOD, the modified model) each circle.

 
Principal component analysis of whole data set. The complete data sets of the in vivo and in vitro studies were used for PCA (Fig. 3). The two major components within the total variation between in vivo and both in vitro models, principal component (PC) 1 and PC 2, explain 40% of the variance in the data set. The results show that genomics responses induced by the low and mid dose in the standard system could not be distinguished from controls, whereas a dose-response effect of coumarin in vivo and in the modified system was evident along the PC 2 axis, explaining 14% of the variance in the whole data set. Along the PC 1 axis, 26% of the variance in the data was explained by genes reacting differently in vivo and in vitro.


    Discussion
 Top
 Abstract
 Materials and Methods
 Results
 Discussion
 References
 
To investigate the relevance of hepatocyte sandwich cultures toward mimicking aspects of in vivo toxicity, cytotoxicity measures and gene expression profiles were compared with data from an in vivo rat study. Coumarin, a compound for which toxicity depends on metabolism by the P450 enzyme system, was used as a model compound. Specifically, P450 enzymes of the 1A and 2E subfamilies convert coumarin to the toxic metabolite coumarin 3,4-epoxide (CE) (Born et al., 2002Go; Lake et al., 2002Go). CE rearranges spontaneously to the more stable o-hydroxyphenylacetaldehyde (o-HPA). Both CE and o-HPA are assumed to contribute to the hepatotoxicity of coumarin in rats as they conjugate with critical cellular macromolecules (Lake et al., 1989Go, 1994Go; Fentem and Fry, 1993Go; Born et al., 1997Go,2000bGo, 2003Go).

Before engaging in gene expression studies, it was hypothesized that an increased metabolic competence of sandwich cultures would improve metabolic conversion of coumarin to toxic metabolites, thereby increasing the resemblance with aspects of in vivo coumarin-induced hepatotoxicity. Therefore, sandwich-cultured hepatocytes were maintained either in a standard medium or in a modified medium that was enriched with a mixture of low concentrations of known P450 inducers to enhance metabolic capacity (A. S. Kienhuis, H. M. Wortelboer, W. J. Maas, M. van Herwijnen, J. C. S. Kleinjans, J. H. M. Delft, and R. H. Stierum, manuscript submitted for publication). Findings in this study showed that both traditional measures of toxicity and coumarin-induced gene expression in the modified model were more pronounced and closer to in vivo compared with the standard model.

In vivo, coumarin toxicity was indicated by increased plasma ALAT, ASAT, and GGT levels, primarily at the highest dose level of 200 mg/kg b.wt. Moreover, hepatotoxicity did manifest itself as severe single cell necrosis and minimal centrilobular necrosis, as observed by histopathological observations in livers of rats treated with 200 mg/kg coumarin. These findings agree with other studies which report that administration of coumarin in doses ranging from 125 to 500 mg/kg results in severe centrilobular necrosis after 24 h (Lake et al., 1989Go, 1994Go; Lake, 1999Go).

In vitro, coumarin appeared to be cytotoxic in hepatocytes cultured in modified medium, whereas no cytotoxicity was measured in the standard model. In our laboratory, it has been shown that gene expression of CYP1A2 and CYP2E1 in the modified model is closer to in vivo compared with the standard model, in which gene expression of these genes is highly down-regulated. Furthermore, upon addition of the inducer mixture, enzyme activities of the 1A subfamily were increased to levels comparable to in vivo (A. S. Kienhuis, H. M. Wortelboer, W. J. Maas, M. van Herwijnen, J. C. S. Kleinjans, J. H. M. Delft, and R. H. Stierum, manuscript submitted for publication). Since coumarin is converted by CYP1A and CYP2E to o-HPA (Born et al., 2002Go; Lake et al., 2002Go), this metabolite possibly causes the cytotoxicity in the modified model. It was shown earlier that the compounds in the inducer mixture, PB, DEX, and ß-NF, increase the formation of o-HPA and other coumarin metabolites, enhancing coumarin toxicity (Peters et al., 1991Go; Fentem and Fry, 1992Go). In the present study, these findings were confirmed by measurement of the coumarin metabolite o-HPAA, which is formed after oxidation of o-HPA. After 24 h of coumarin exposure, a dose-dependent formation of o-HPAA was exclusively observed in hepatocytes cultured in the modified model, whereas no o-HPAA was formed in the standard model. Furthermore, coumarin clearance was significantly higher in the modified culture medium compared with the standard culture medium in all dose groups.

Analysis of gene expression revealed that both in vivo and in vitro, genes are affected by coumarin at dose levels at which no toxicity occurs as determined by traditional toxicology measures. This finding suggests that gene expression changes may very well be more sensitive indicators of potential adverse effects, as was indicated before by Heinloth et al. (2004Go). Considerably more genes were significantly altered and similarly affected by coumarin in both the modified system and in vivo compared with the standard system. In this study, not all affected genes are discussed in detail. It is of more relevance to investigate analogies in the biological pathways these genes trigger. The four pathways that were statistically significantly triggered both in vivo and in the modified model, and the significant differentially expressed genes within these pathways will be discussed.

Betaine-homocysteine methyltransferase (Bhmt) and CTL target antigen (Cth), a putative cystathione {gamma} lyase enzyme, both down-regulated in vivo and in the modified system, appear in the methionine metabolism pathway. It has been shown that down-regulation of Bhmt and Cth impairs conversion of homocysteine to methionine and cysteine, respectively, resulting in increased homocysteine levels in the cell (Torres et al., 1999Go; Forestier et al., 2003Go). Homocysteine accumulation is associated with impaired liver function, necrosis, liver cirrhosis, and fibrogenesis (Torres et al., 1999Go; Finkelstein, 2003Go).

Expression of CYP1A2, a gene which represents both the {gamma}-hexachlorocyclohexane degradation pathway and the fatty acid metabolism pathway, is significantly down-regulated in vivo and in the modified model. This P450 belongs to the CYP1A subfamily that converts coumarin to toxic metabolites. Therefore, down-regulation of gene expression of this enzyme could be interpreted as a classical negative feedback loop, turning out in the present case to prevent additional liver injury.

Fibrinogen {gamma}-polypeptide (fgg) is significantly up-regulated in vivo and significantly down-regulated in the modified model. This gene is part of the complement and coagulation cascade, to which coumarin is associated since coumarin derivatives, such as warfarin, are well known for their anticoagulant properties (Egan et al., 1990Go). One possible explanation for the significant, but different response of particular genes to coumarin in vivo and in vitro can be illustrated by fgg. fgg is mainly expressed in hepatocytes. Up-regulation of fgg, as occurred in vivo, would be expected in terms of toxicity since fibrinogen levels are increased after liver injury (Redman and Xia, 2001Go). Fibrinogen synthesis has been shown to be stimulated in hepatocytes by factors excreted by extrahepatic tissues or nonparenchymal cells (Otto et al., 1987Go). In vitro, these factors are absent. Therefore, stimulation of fibrinogen synthesis in hepatocyte cultures in vitro may only occur when these factors are added to the culture medium.

In addition to the statistical analysis of microarray data which resulted in a list of significantly modulated genes, PCA was performed. PCA allows inclusion of all genes into analysis, preserving their interrelationships. Results showed that even though differences between in vitro and in vivo data remain evident, the dose-response effect with respect to gene expression of coumarin in the modified system is more similar to in vivo data compared with the standard system.

Dose levels in vivo were compared with the coumarin concentrations applied in vitro. Using pharmacokinetic data from literature, a distribution coefficient of coumarin of 3.33 l/kg was estimated (Ritschel and Hussain, 1988Go; Lake, 1999Go; Born et al., 2003Go). Assuming that in vivo bioavailability of coumarin is 100% after i.p. injection, one can estimate a maximum in vivo plasma concentration of 36 µM, 154 µM, and 411 µM for the low, mid, and high in vivo doses, respectively. These concentrations are in the same range as the in vitro concentrations of 70 µM, 200 µM, and 600 µM. Although plasma levels of coumarin were not determined in vivo, the range of estimated maximum plasma levels in vivo and, therefore, potential target organ concentrations were comparable to the dose level range used in the culture media in vitro.

It is important to note that in addition to coumarin, the inducer mixture may also influence gene expression and that these effects interact with each other. In the design of the present study, however, the effect of the inducer mixture on gene expression was filtered out because matching controls were used; i.e., both test and control samples are from hepatocytes cultured in the modified system. Nonetheless, it remains possible that the inducer mixture affected the expression of some genes in the modified system in such a way that no additional effect of coumarin could have been detected.

Furthermore, effects on gene expression in the highest dose group of the modified model may not be solely attributed to coumarin. The highest coumarin dose group in this model resulted in a 50% loss of cell viability as measured by the MTT reduction assay. In a sandwich configuration, dead cells remain between the two layers of collagen. As a consequence, necrosis caused by coumarin exposure in one cell can affect gene expression in viable neighbor cells (Fielden and Zacharewski, 2001Go). Excluding the high dose group of the modified model from analysis, still 80% of the genes in Table 2 are retained. Moreover, the methionine metabolism and {gamma}-hexachlorocyclohexane degradation pathways are still significantly triggered by coumarin.

In summary, our experiments have shown that the metabolism of coumarin toward active metabolites, coumarin-induced toxicity, gene expression profiles, and, consequently, biological pathways in the modified system containing sandwich-cultured hepatocytes with enhanced metabolic capacity better represent the situation in vivo compared with conventional sandwich-cultured hepatocytes. This finding highlights the need for a metabolically competent, toxicogenomics-based, hepatocyte in vitro system.


    Acknowledgments
 
We thank M. Schut, Dr. N. Treijtel, and W. Maas for help in hepatocyte isolation and culture techniques; M. van den Wijngaard for sample isolation of the in vivo study; Dr. W. Heijne for help with the experimental design; A. de Kat Angelino-Bart, M. Havekes, and Dr. F. Schuren at the microarray facility; M. Dansen for microarray data quality analysis; A. Freidig for help with pharmacokinetic calculation; and Dr. S. Bijlsma for assistance in principal component analysis.


    Footnotes
 
Financial support was provided by The Netherlands Organisation for Health Research and Development, program Alternatives to Animal Experiments (3170.0049) and the Dutch Ministry of Economic Affairs. Financial support provided by Servier Nederland B.V. is greatly appreciated.

Article, publication date, and citation information can be found at http://dmd.aspetjournals.org.

doi:10.1124/dmd.106.011262.

ABBREVIATIONS: P450, cytochrome P450; PB, phenobarbital; DEX, dexamethasone; ß-NF, ß-naphthoflavone; DMEM, Dulbecco's modified Eagle's medium; MTT, 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide; ALAT, alanine aminotransferase; ASAT, aspartate aminotransferase; LDH, lactate dehydrogenase; GGT, {gamma}-glutamyl transferase; DMSO, dimethyl sulfoxide; o-HPAA, o-hydroxyphenylacetic acid; TIGR, The Institute for Genomic Research; BRB, Biometric Research Branch; PCA, principal component analysis; PC, principal component; CE, coumarin 3,4-epoxide; o-HPA, o-hydroxyphenylacetaldehyde.

Address correspondence to: Rob H. Stierum, TNO Quality of Life, Business Unit Biosciences, Physiological Genomics, P.O. Box 360 (PP 8), 3700 AJ, Zeist, The Netherlands. E-mail: rob.stierum{at}voeding.tno.nl.


    References
 Top
 Abstract
 Materials and Methods
 Results
 Discussion
 References
 


Balls M, Bogni A, Bremer S, Casati S, Coecke S, Eskes C, Prieto P, Sabbioni E, Worth A, Zuang V, et al. (2002) Alternative (non-animal) methods for chemicals testing: currents status and future prospects. Altern Lab Anim 30 (Suppl 2): 1–125.[Medline]

Beken S, Vanhaecke T, De Smet K, Pauwels M, Vercruysse A, and Rogiers V (2004) Collagengel cultures of rat hepatocytes: collagen-gel sandwich and immobilization cultures, in Cytochrome P450 Protocols (Shephard EA and Phillips IR eds), pp 303–310, Humana Press Inc., Totowa, NJ.

Boess F, Kamber M, Romer S, Gasser R, Muller D, Albertini S, and Suter L (2003) Gene expression in two hepatic cell lines, cultured primary hepatocytes, and liver slices compared to the in vivo liver gene expression in rats: possible implications for toxicogenomics use of in vitro systems. Toxicol Sci 73: 386–402.[Abstract/Free Full Text]

Born SL, Api AM, Ford RA, Lefever FR, and Hawkins DR (2003) Comparative metabolism and kinetics of coumarin in mice and rats. Food Chem Toxicol 41: 247–258.[CrossRef][Medline]

Born SL, Caudill D, Fliter KL, and Purdon MP (2002) Identification of the cytochromes P450 that catalyze coumarin 3,4-epoxidation and 3-hydroxylation. Drug Metab Dispos 30: 483–487.[Abstract/Free Full Text]

Born SL, Caudill D, Smith BJ, and Lehman-McKeeman LD (2000a) In vitro kinetics of coumarin 3,4-epoxidation: application to species differences in toxicity and carcinogenicity. Toxicol Sci 58: 23–31.[Abstract/Free Full Text]

Born SL, Hu JK, and Lehman-McKeeman LD (2000b) o-Hydroxyphenylacetaldehyde is a hepatotoxic metabolite of coumarin. Drug Metab Dispos 28: 218–223.[Abstract/Free Full Text]

Born SL, Rodriguez PA, Eddy CL, and Lehman-McKeeman LD (1997) Synthesis and reactivity of coumarin 3,4-epoxide. Drug Metab Dispos 25: 1318–1324.[Abstract/Free Full Text]

Dambach DM, Andrews BA, and Moulin F (2005) New technologies and screening strategies for hepatotoxicity: use of in vitro models. Toxicol Pathol 33: 17–26.[Abstract/Free Full Text]

Dunn JC, Yarmush ML, Koebe HG, and Tompkins RG (1989) Hepatocyte function and extracellular matrix geometry: long-term culture in a sandwich configuration. FASEB J 3: 174–177.[Abstract]

Egan D, O'Kennedy R, Moran E, Cox D, Prosser E, and Thornes RD (1990) The pharmacology, metabolism, analysis, and applications of coumarin and coumarin-related compounds. Drug Metab Rev 22: 503–529.[Medline]

Fentem JH and Fry JR (1992) Metabolism of coumarin by rat, gerbil and human liver microsomes. Xenobiotica 22: 357–367.[Medline]

Fentem JH and Fry JR (1993) Species differences in the metabolism and hepatotoxicity of coumarin. Comp Biochem Physiol C 104: 1–8.[Medline]

Fielden MR and Zacharewski TR (2001) Challenges and limitations of gene expression profiling in mechanistic and predictive toxicology. Toxicol Sci 60: 6–10.[Abstract/Free Full Text]

Finkelstein JD (2003) Methionine metabolism in liver diseases. Am J Clin Nutr 77: 1094–1095.[Free Full Text]

Forestier M, Banninger R, Reichen J, and Solioz M (2003) Betaine homocysteine methyltransferase: gene cloning and expression analysis in rat liver cirrhosis. Biochim Biophys Acta 1638: 29–34.[Medline]

Guillouzo A (1998) Liver cell models in in vitro toxicology. Environ Health Perspect 106 (Suppl 2): 511–532.[CrossRef][Medline]

Hamadeh HK, Bushel PR, Jayadev S, Martin K, DiSorbo O, Sieber S, Bennett L, Tennant R, Stoll R, Barrett JC, et al. (2002) Gene expression analysis reveals chemical-specific profiles. Toxicol Sci 67: 219–231.[Abstract/Free Full Text]

Heijne WH, Jonker D, Stierum RH, van Ommen B, and Groten JP (2005) Toxicogenomic analysis of gene expression changes in rat liver after a 28-day oral benzene exposure. Mutat Res 575: 85–101.[Medline]

Heinloth AN, Irwin RD, Boorman GA, Nettesheim P, Fannin RD, Sieber SO, Snell ML, Tucker CJ, Li L, Travlos GS, et al. (2004) Gene expression profiling of rat livers reveals indicators of potential adverse effects. Toxicol Sci 80: 193–202.[Abstract/Free Full Text]

Jaeschke H (2003) Are cultured liver cells the right tool to investigate mechanisms of liver disease or hepatotoxicity? Hepatology 38: 1053–1055.[Medline]

Jessen BA, Mullins JS, De Peyster A, and Stevens GJ (2003) Assessment of hepatocytes and liver slices as in vitro test systems to predict in vivo gene expression. Toxicol Sci 5: 208–222.

Lake BG (1999) Coumarin metabolism, toxicity and carcinogenicity: relevance for human risk assessment. Food Chem Toxicol 37: 423–453.[CrossRef][Medline]

Lake BG, Evans JG, Chapuis F, Walters DG, and Price RJ (2002) Studies on the disposition, metabolism and hepatotoxicity of coumarin in the rat and Syrian hamster. Food Chem Toxicol 40: 809–823.[CrossRef][Medline]

Lake BG, Evans JG, Lewis DF, and Price RJ (1994) Studies on the acute effects of coumarin and some coumarin derivatives in the rat. Food Chem Toxicol 32: 357–363.[CrossRef][Medline]

Lake BG, Gray TJ, Evans JG, Lewis DF, Beamand JA, and Hue KL (1989) Studies on the mechanism of coumarin-induced toxicity in rat hepatocytes: comparison with dihydrocoumarin and other coumarin metabolites. Toxicol Appl Pharmacol 97: 311–323.[CrossRef][Medline]

LeCluyse E, Bullock P, Madan A, Carroll K, and Parkinson A (1999) Influence of extracellular matrix overlay and medium formulation on the induction of cytochrome P-450 2B enzymes in primary cultures of rat hepatocytes. Drug Metab Dispos 27: 909–915.[Abstract/Free Full Text]

LeCluyse E, Bullock P, and Parkinson A (1996) Strategies for restoration and maintenance of normal hepatic structure and function in long-term cultures of rat hepatocytes. Adv Drug Delivery Rev 22: 133–186.[CrossRef]

Mosmann T (1983) Rapid colorimetric assay for cellular growth and survival: application to proliferation and cytotoxicity assays. J Immunol Methods 65: 55–63.[CrossRef][Medline]

Nuwaysir EF, Bittner M, Trent J, Barrett JC, and Afshari CA (1999) Microarrays and toxicology: the advent of toxicogenomics. Mol Carcinog 24: 153–159.[CrossRef][Medline]

Otto JM, Grenett HE, and Fuller GM (1987) The coordinated regulation of fibrinogen gene transcription by hepatocyte-stimulating factor and dexamethasone. J Cell Biol 105: 1067–1072.[Abstract/Free Full Text]

Paine AJ, Williams LJ, and Legg RF (1979) Determinants of cytochrome P-450 in liver cell cultures, in The Liver: Quantitative Aspects of Structure and Function (Preisig R and Bircher J eds), pp 99–109, Editio Cantor, Aulendorf.

Peters MM, Walters DG, van Ommen B, van Bladeren PJ, and Lake BG (1991) Effect of inducers of cytochrome P-450 on the metabolism of [3–14C]coumarin by rat hepatic microsomes. Xenobiotica 21: 499–514.[Medline]

Redman CM and Xia H (2001) Fibrinogen biosynthesis. Assembly, intracellular degradation, and association with lipid synthesis and secretion. Ann NY Acad Sci 936: 480–495.[Abstract/Free Full Text]

Richert L, Binda D, Hamilton G, Viollon-Abadie C, Alexandre E, Bigot-Lasserre D, Bars R, Coassolo P, and LeCluyse E (2002) Evaluation of the effect of culture configuration on morphology, survival time, antioxidant status and metabolic capacities of cultured rat hepatocytes. Toxicol In Vitro 16: 89–99.[CrossRef][Medline]

Ritschel WA and Hussain SA (1988) Transdermal absorption and topical bioavailability of coumarin. Methods Find Exp Clin Pharmacol 10: 165–169.[Medline]

Seglen PO (1976) Preparation of isolated rat liver cells. Methods Cell Biol 13: 29–83.[Medline]

Simon R and Peng A (2002) BRB ArrayTools, Users Guide. National Cancer Institute, Biometric Research Branch Technical Report 008, National Cancer Institute, Bethesda.

Stierum R, Heijne W, Kienhuis A, van Ommen B, and Groten J (2005) Toxicogenomics concepts and applications to study hepatic effects of food additives and chemicals. Toxicol Appl Pharmacol 207: 179–188.[Medline]

Torres L, Garcia-Trevijano ER, Rodriguez JA, Carretero MV, Bustos M, Fernandez E, Eguinoa E, Mato JM, and Avila MA (1999) Induction of TIMP-1 expression in rat hepatic stellate cells and hepatocytes: a new role for homocysteine in liver fibrosis. Biochim Biophys Acta 1455: 12–22.[Medline]

van de Peppel J, Kemmeren P, van Bakel H, Radonjic M, van Leenen D, and Holstege FC (2003) Monitoring global messenger RNA changes in externally controlled microarray experiments. EMBO Rep 4: 387–393.[CrossRef][Medline]

Waring JF, Ciurlionis R, Jolly RA, Heindel M, and Ulrich RG (2001) Microarray analysis of hepatotoxins in vitro reveals a correlation between gene expression profiles and mechanisms of toxicity. Toxicol Lett 1: 359–368.

Waring JF and Ulrich RG (2000) The impact of genomics-based technologies on drug safety evaluation. Annu Rev Pharmacol Toxicol 40: 335–352.[CrossRef][Medline]

Yang YH, Dudoit S, Luu P, Lin DM, Peng V, Ngai J, and Speed TP (2002) Normalization for cDNA microarray data: a robust composite method addressing single and multiple slide systematic variation. Nucleic Acids Res 30: e15.[Abstract/Free Full Text]



This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
dmd.106.011262v1
34/12/2083    most recent
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Kienhuis, A. S.
Right arrow Articles by Stierum, R. H.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Kienhuis, A. S.
Right arrow Articles by Stierum, R. H.


Home Help [Feedback] [For Subscribers] [Archive] [Search] [Contents]
All ASPET Journals Molecular Pharmacology Pharmacological Reviews
 Molecular Interventions Drug Metabolism and Disposition