Modulation of Gene Expression by Cancer Chemopreventive Dithiolethiones through the Keap1-Nrf2 Pathway

IDENTIFICATION OF NOVEL GENE CLUSTERS FOR CELL SURVIVAL*

Mi-Kyoung KwakDagger , Nobunao WakabayashiDagger §, Ken Itoh§, Hozumi Motohashi§, Masayuki Yamamoto§, and Thomas W. KenslerDagger

From the Dagger  Department of Environmental Health Sciences, Johns Hopkins University Bloomberg School of Public Health, Baltimore, Maryland 21205 and the § Institute of Basic Medical Sciences and Center for Tsukuba Advanced Research Alliance, University of Tsukuba, Tennoudai, Tsukuba 305-8577, Japan

Received for publication, November 21, 2002, and in revised form, December 18, 2002

    ABSTRACT
TOP
ABSTRACT
INTRODUCTION
EXPERIMENTAL PROCEDURES
RESULTS
DISCUSSION
REFERENCES

Enzyme inducers such as 3H-1,2-dithiole-3-thione (D3T) enhance the detoxication of environmental carcinogens and protect against neoplasia. The putative molecular sensor for inducers is Keap1, a sulfhydryl-rich protein that sequesters the transcription factor Nrf2 in the cytoplasm. Expression of these detoxication enzymes is blunted in nrf2-deficient mice; moreover, these mice are more sensitive to carcinogenesis, and the protective actions of dithiolethiones are lost with nrf2 disruption. Hepatic gene expression profiles were examined by oligonucleotide microarray analysis in vehicle- or D3T-treated wild-type mice as well as in nrf2 single and keap1-nrf2 double knockout mice to identify those genes regulated by the Keap1-Nrf2 pathway. Transcript levels of 292 genes were elevated in wild-type mice 24 h after treatment with D3T; 79% of these genes were induced in wild-type, but not nrf2-deficient mice. These nrf2-dependent, D3T-inducible genes included known detoxication and antioxidative enzymes. Unexpected clusters included genes for chaperones, protein trafficking, ubiquitin/26 S proteasome subunits, and signaling molecules. Gene expression patterns in keap1-nrf2 double knockout mice were similar to those in nrf2-single knockout mice. D3T also led to nrf2-dependent repression of 31 genes at 24 h; principally genes related to cholesterol/lipid biosynthesis. Collectively, D3T increases the expression of genes through the Keap1-Nrf2 signaling pathway that directly detoxify toxins and generate essential cofactors such as glutathione and reducing equivalents. Induction of nrf2-dependent genes involved in the recognition and repair/removal of damaged proteins expands the role of this pathway beyond primary control of electrophilic and oxidative stresses into secondary protective actions that enhance cell survival.

    INTRODUCTION
TOP
ABSTRACT
INTRODUCTION
EXPERIMENTAL PROCEDURES
RESULTS
DISCUSSION
REFERENCES

Inducers of phase 2 and antioxidative enzymes are known to enhance the detoxication of environmental carcinogens in animals, often leading to protection against neoplasia (1, 2). Use of enzyme inducers as cancer chemopreventive agents in humans is currently under clinical investigation (3, 4). Regulation of both basal and inducible expression of cytoprotective genes is mediated in part by the antioxidant response element (ARE),1 a cis-acting sequence found in the 5'-flanking region of genes encoding many phase 2 enzymes including mouse glutathione S-transferase (GST) A1, human NAD(P)H quinone oxidoreductase (NQO1), and human gamma -glutamylcysteine ligase as well as murine Nrf2 itself (5-8). Recently, Maf and CNC-bZIP ("cap"n collar family of basic region leucine zipper proteins) transcription factors such as Nrf2 have been shown to be ARE-binding proteins (9, 10). Overexpression of Nrf2 in human hepatoma cells enhanced both basal and inducible activation of an ARE reporter gene (10). Increased nuclear accumulation of Nrf2 has been observed in the liver of mice treated with 3H-1,2-dithiole-3-thione (D3T), a potent chemopreventive agent (8, 11). Initially, this accumulation results from translocation of Nrf2 from the cytoplasm. Itoh et al. (12) have identified and characterized Keap1, an actin-binding protein localized to the cytoplasm that sequesters Nrf2 by specific binding to its N-terminal regulatory domain (13). Administration of sulfhydryl reactive compounds, such as diethylmaleate or D3T, abolishes Keap1 repression of Nrf2 activity in cells and facilitates the nuclear accumulation of Nrf2 (14, 15). Recently, Dinkova-Kostova et al. (16) has shown in a cell-free system that selective cysteine amino acids in Keap1 can react directly with a sufhydryl reactive inducer, dexamethasone mesylate, and trigger the release of Nrf2 from Keap1. These results support the hypothesis that Keap1 is a key regulatory molecule of Nrf2 and that the Keap1-Nrf2 complex acts as a sensor to oxidative or electrophilic stresses to induce protective genes promoting cell survival.

Studies with nrf2-disrupted mice indicated that Nrf2 was essential for the induction of GST and NQO1 activities in vivo by different classes of chemopreventive agents including dithiole-thiones, isothiocyanates, and phenolic antioxidants (9, 11, 17). A series of studies from several laboratories have been undertaken over the past decade to better understand the range of genes coordinately induced by these agents so as to elucidate the biochemical basis for protection. Genes now recognized as induced through the Nrf2-ARE signaling system include a panel of xenobiotic conjugating enzymes, enzymes that provide cofactors and reducing equivalents for these reactions, and antioxidative enzymes and proteins (11, 17-20). Collectively, these gene products facilitate the detoxication and elimination of toxic electrophilic and free radical metabolites of xenobiotics. The likely importance of these protective enzymes is highlighted by recent observations that nrf2-knockout mice were considerably more sensitive to the acute toxicities of acetaminophen, butylated hydroxytoluene, and hyperoxia (21-24). These mice also form higher levels of DNA adducts following exposure to carcinogens such as aflatoxin B1, diesel particulate matter, and benzo[a]pyrene (25-27). Moreover, nrf2-disrupted mice are substantially more sensitive to the tumorigenicity of benzo[a]pyrene (28).

Dithiolethiones such as D3T, anethole dithiolethione, and oltipraz inhibit the toxicity and carcinogenicity of many chemical carcinogens in multiple target organs (29-32) and are undergoing preclinical and clinical evaluations for use as cancer chemopreventive agents. Interestingly, the chemoprotective efficacy of dithiolethiones and other enzyme inducers is lost in nrf2-knockout mice (28). Gene-disrupted animals provide elegant models for assessing the molecular mechanisms underlying the pharmacodynamic actions of chemopreventive agents. In this study, we have used nrf2-knockout mice and keap1-nrf2 double knockout mice to probe the role of the Keap1-Nrf2-ARE signaling pathway in protecting cells against stress conditions. Through the use of oligonucleotide array technology it is apparent that not only are primary defense systems involving detoxication enzymes induced through this pathway, but that secondary defense mechanisms involving surveillance and repair of damaged biomolecules are also activated.

    EXPERIMENTAL PROCEDURES
TOP
ABSTRACT
INTRODUCTION
EXPERIMENTAL PROCEDURES
RESULTS
DISCUSSION
REFERENCES

Animals and Treatments-- Wild-type and nrf2-disrupted mice were generated from inbred nrf2 heterozygous mice (9, 11). Keap1-nrf2 double knockout mice were generated by mating nrf2 -/- mice with keap1 +/- mice to obtain the double homozygous mice (33). Genotypes were determined by PCR using primers TGGACGGGACTATTGAAGGCTG, GCCGCCTTTTCAGTAGATGGAGG, and GCGGATTGACCGTAATGGGATAGG for nrf2, and CGGGATCCCCATGGAAAGGCTTATTGAGTTC, GAAGTGCATGTAGATATACTCCC, and TCAGAGCAGCCGATTGTCTGTTGTGCCCAGTCA for keap1. For the microarray analyses three male mice were used per group. Mice were fed AIN-76A semipurified diet for at least 10 days before treatment. 10-12-week-old mice were treated with 0.5 mmol/kg of D3T by gavage in a suspension consisting of 1% Cremophor and 25% of glycerol. Mice were sacrificed 6 h or 24 h after treatment, and livers were removed and snap-frozen. Additional mice were fed basal AIN-76A diet or diet supplemented with 0.03% (w/w) D3T for 7 days and sacrificed at day 8. The selected doses of D3T are the respective maximum tolerated doses for gavage and chronic dietary administration and do not affect body weight or other clinical signs over longer periods. It is estimated that the dietary dose provided about half of the amount of D3T per day that was administered in the single oral dose. Keap1-nrf2 double knockout mice were treated with either vehicle or D3T (0.5 mmol/kg) by gavage and sacrificed 24 h after treatment for microarray analysis of hepatic gene expression.

Sample Preparation for Microarray Analysis-- Total RNA was isolated from frozen liver using the Totally RNA kit (Ambion, Austin, TX). Isolated RNA was then purified using the RNeasy Mini kit (Qiagen, Valencia, CA). cDNA was synthesized from total RNA using the Superscript Choice kit (Invitrogen, Carlsbad, CA) with a T7-(dT)24 primer. Biotin-labeled cRNA was prepared from cDNA by in vitro transcription (Enzo Biochemical, Farmingdale, NY) and fragmented by incubation at 94 °C for 45 min in 40 mM Tris acetate buffer, pH 8.1, with 100 mM potassium and 30 mM magnesium acetate.

Hybridization and Scanning-- Fragmented cRNA (12.5 µg) was hybridized at 45 °C for 16 h to a Murine Genome U74Av2 GeneChip® (Affymetrix, Santa Clara, CA) containing probes for ~12,400 genes and expressed sequence tags (ESTs). Genechips were washed and stained using a fluid station and scanned using an Affymetrix Genechip system confocal scanner.

Data Analysis-- Affymetrix Microarray Suite 4.1 was used for expression analysis and fluorescence intensity was measured for each probe array and normalized to the average fluorescence intensity for the entire probe array. Pairwise comparison between data from individual mice (n = 3), generating 9 comparisons, was performed. The change in gene expression for each gene was calculated as a ratio of average difference(D3T-treated)/average difference(vehicle-treated) and average fold changes from nine comparisons were obtained. To evaluate the reproducibility of paired comparisons, coefficients of variation (CV, S.D./mean) for the average difference change was applied, and 1.0 was used as a cutoff value (34). After evaluating reproducibility, genes having comparison numbers of >= 7 were selected to filter out the false positive responses. Genespring software (Silicon Genetics, San Carlos, CA) and the Affymetrix Analysis Center website were used for the annotation of genes. For analysis of the distributions of gene families, expression levels of each gene category were plotted as histograms as described by Bouton and Pevsner (35).

RT-PCR Analysis-- RT-PCR was performed on genes induced or repressed by D3T to confirm the microarray results. For the synthesis of cDNAs, 50 ng of total RNA was incubated for 20 min with 10 mM Tris (pH 8.4), 5 mM KCl, 5 mM MgCl2, 4 mM dNTPs, 0.125 µg of oligo(dT)12-18 and 30 units of M-MRV (Moloney Murine Leukemia Virus) reverse transcriptase (Invitrogen). PCR amplification for each gene was performed with a Fail Safe PCR kit (Epicentere, Madison, WI) using a DNA thermal cycler (MJ Research, Watertown, MA). Amplification conditions were 26-29 cycles for 5 min at 95 °C, 30 s at 56 °C and 40 s at 72 °C. PCR primers specific to NQO1, mEH, and albumin were used as described previously (11). Other primers were synthesized by Integrated DNA Technology (Coralville, IA) and were as follows: GST Mu3 (GenBankTM accession number J03953) CAAGTTATGGACACCCGCAT and AGGCACTTGGGCTCAAACAT; UDP-glucose dehydrogenase (X06358) GACATGAATGACTACCAGAG and GTACGGAATTCTCTTGGAAG; heat shock protein 40 homolog b11 (AW122551) TACGATGATATCACTTCCTC and GGCGAGATTTCTATAAGATC; p58 (U28423) GTACGATGAAGCCATTCAGG and CCTTGCCATGAGTTCCAACT; proteasome subunit beta type 3 (AW045339) TTCAGCGTCCTGGTGGTGAT and ACAGAGCCTGTCATTGCTGG; 26 S proteasome-associated pad 1 homolog (Y13071); TATCAACACTCAGCAGAGCT and AATCCTTCCATCCAACTCT; 26 S proteasome ATPase 3 (AA409481) CCCTCAGACCACGGA and TGTGTGAGCGGGTATGAT; mCAR2 (AF009328) TTGTGCAGTTCAAGCCTCCAG and GCGTCCTCCATCTTGTAGCAA; carbonic anhydrase 3 (AJ006474) AACAGGAGCAAGAAAGAG and AGGTCTTCCCATTGTTCA. PCR products were electrophoresed on 1.8% agarose gels and the gel image captured and quantified with a VersaDoc Imaging System (BioRad).

    RESULTS
TOP
ABSTRACT
INTRODUCTION
EXPERIMENTAL PROCEDURES
RESULTS
DISCUSSION
REFERENCES

D3T-altered Gene Expression in Mouse Liver-- Oligonucleotide microarray was used to analyze the gene expression profiles in mouse liver following D3T treatment with a primary goal of identifying those genes that are regulated through the Keap1/Nrf2 signaling pathway. Expression levels of 352 genes (2.84% of the total genes on the array) were changed in the liver of wild-type mice 24 h after treatment with the potent chemoprotective agent D3T. 129 and 292 genes were increased 6 and 24 h after D3T treatment, respectively, while only 34 were elevated after feeding D3T for 1 week. Most of the D3T-induced genes were dependent on Nrf2 as 82, 79, and 97% of D3T-inducible genes were increased in wild-type mice but not in nrf2-disrupted mice in the 6 h, 24 h, and dietary treatment groups, respectively (Fig. 1). A much smaller number of genes were down-regulated by D3T (60 genes at 24 h). Also shown in Fig. 1, D3T suppressed the expression of 31 genes 24 h after treatment and 30 genes following dietary administration in wild-type mice only. Expression of most of these Nrf2-dependent, D3T-modulated genes was not altered in the keap1-nrf2 double knockout mice 24 h after treatment with D3T. Only 15 out of the 231 Nrf2-dependent, D3T-inducible genes and 2 of the Nrf2-dependent, D3T-repressed genes were affected by D3T in these double knockout mice. Thus, D3T treatment altered the expression of only a modest subset of genes in the liver of wild-type mice and most of these genes were regulated through the Keap1/Nrf2 pathway. Comparisons of expression patterns at different time points indicated that maximal changes occurred 24 h after a single administration of D3T.


View larger version (19K):
[in this window]
[in a new window]
 
Fig. 1.   Modulation of Nrf2-dependent gene expression by D3T in mouse liver. Hepatic gene expression patterns were analyzed 6 or 24 h after a single administration of D3T and following addition of D3T to the diet for 7 days. The number of genes increased or decreased only in wild-type mice, not in nrf2-knockout mice, are indicated for each group.

Nrf2-dependent, D3T-inducible Genes at 24 h-- Genes up-regulated by D3T in wild-type, but not nrf2-disrupted mice, were regarded as Nrf2-dependent, D3T-inducible genes. The 231 Nrf2-dependent genes elevated 24 h after treatment with D3T can be classified in several functional categories including xenobiotic metabolizing enzymes, antioxidants, chaperone/stress response proteins, ubiquitin/proteasome systems, acute response/immunity proteins, cytoskeletal organization proteins, protein trafficking proteins, membrane transport proteins, signaling molecules, transcription factors and regulators, and RNA processing/translation-related factors (Table I). As expected from earlier studies, a substantial number (26) of genes related to xenobiotic metabolism were induced by D3T treatment in a nrf2-dependent manner. This gene category includes seven genes belonging to cytochrome P450 subfamilies. Ten genes related to glutathione synthesis and conjugation were induced by D3T only in wild-type mice and those included seven subunits of GSTs, glutathione reductase, and the regulatory and catalytic subunits of gamma -glutamylcysteine ligase. NQO1, mEH, and UDP-glucuronosyl transferase 2 were also induced. Flavin-containing monooxygenase-1, carbonyl reductase and aldo-keto reductase family 1 were previously unrecognized D3T-inducible genes in mouse liver regulated through Nrf2. These results confirmed our previous observations that xenobiotic detoxifying genes are a major category of genes affected by D3T and nrf2 genotype. Several antioxidative genes were also induced by D3T; thioredoxin, thioredoxin reductase, and peroxiredoxin transcripts were increased by 2.2-, 3.2-, and 1.7-fold, respectively, 24 h after treatment. General enzymes including nucleotide diphosphatase (3.9-fold), aldolase 1A (2.2-fold), and amino levulinate synthase 1 (3.9-fold) were increased only in wild-type mice. UDP-glucose dehydrogenase was increased by D3T only in wild-type mice 6 h after treatment and in both genotypes 24 h after treatment; however, the induction in wild-type (7.5-fold) was much higher than in nrf2-disrupted (1.7-fold) mice.

                              
View this table:
[in this window]
[in a new window]
 
Table I
Hepatic Keap1-Nrf2-dependent, D3T-inducible genes

Twenty genes related to proteasomes as well as genes involved in ubiquitination (ubiquitin-conjugating enzymes E2-32, Huntington-interacting protein-2 and valosin-containing protein) were identified as Nrf2-dependent, D3T-induced genes. The average fold changes of these proteasome subunits ranged from 1.8 to 3.5. Five additional proteasome subunits (proteasome subunit alpha types 2, 5, and 7; delta type 4; and ATPase subunit 2) were induced in both wild-type and nrf2-disrupted mice, but with a higher AFC in the wild-type mice (AFC of wild-type/knockout, subunit alpha 2, 2.5/1.4; alpha 5, 2.9/1.6; alpha 7, 2.7/1.4; delta 4, 3.5/1.8; and subunit ATPase 2, 2.9/1.5). Further, proteasome subunit delta 4 was induced equally by D3T in wild-type and nrf2-knockout mice. Therefore, in the aggregate, induction of 25 proteasome subunit genes by D3T was affected by nrf2 genotype.

Other major categories of genes induced by D3T through a Nrf2 pathway were chaperone-protein trafficking systems. These included heat shock proteins, ClpP protease, peptidylproryly isomerase B, and FK506-binding protein 5, which are involved in the repair of unfolded/misfolded proteins as well as genes for protein transport across the endoplasmic reticulum such as the vesicle docking proteins and Sec61. Acute response and inflammatory proteins such as serum amyloids and CD8b and cytoskeletal organization/binding proteins including profilin 2 and protein 4.1 were increased by D3T in wild-type mouse liver. Signaling molecules affected by D3T included NFkappa B essential modulator, protein inhibitor of nitric-oxide synthase and GTP-binding proteins (ARF-1, Ran, and Gna11). Some genes related to transcription were also categorized as Nrf2-dependent, D3T-induced genes in mice liver. The transcription factor ATF4 family ATF7 was increased 3.6-fold and CRE/ATF family protein X-box-binding protein 1 was increased 2.1-fold by D3T in wild-type mice only. Nuclear receptor mCAR2, transcription activator COP9 subunit 5, and tumor rejection antigen gp96 were also elevated. Lastly, genes affecting mRNA stability/processing-related (small nuclear RNA and poly(rC)-binding protein) and translation (eIF3-p44 and eIF2a) were induced by D3T.

Nrf2-dependent Genes at Other Time Points-- Transcripts for 106 genes were increased by D3T 6 h after treatment and remained elevated for 54 of these genes 24 h after treatment. Many genes related to xenobiotic metabolism were substantially induced severalfold or more by 6 h, e.g. Cyp2a4 (12.4-fold), NQO1 (5.2-fold), and gamma -glutamylcysteine ligase catalytic subunit (3.6-fold), and remained elevated at 24 h. Redox-regulating molecules such as thioredoxin and thioredoxin reductase also showed rapid induction. Expression of half of the rapidly induced genes returned to basal levels at 24 h after treatment. Crystalline alpha C (2.0-fold), aldehyde oxidase 1 (2.9-fold), antiapoptotic Bag3 (2.5-fold), mATF4 (1.7-fold), and gelsolin (2.5-fold) were increased only 6 h after treatment by D3T in an Nrf2-dependent manner.

Gene expression profiling was also examined following feeding D3T for 7 days for comparison to that of single administration of D3T. Expression levels of most genes, which were induced by single administration of D3T, returned to basal levels despite continuous administration of D3T in the diet for 7 days. Transcripts for 34 genes remained elevated in this group, 33 of which were for Nrf2-dependent genes. Among the xenobiotic metabolizing enzymes, only two genes (Cyp3a13 and Cyp4a10) were sustainably elevated by D3T in the diet and transcripts for other genes such as NQO1 and the GSTs were not induced over this longer time period. Similarly, no proteasome subunits were detected as D3T-inducible genes with D3T in the diet for 7 days. Deoxyribonuclease IIalpha , GRO1 oncogene, ICAM-1, and cytokine inducible SH2-containing protein 3 exhibited large increases in this exposure group. D3T in the diet for 7 days also repressed the expression of 30 genes in an Nrf2-dependent manner. Genes repressed at 24 h as well as this later time point included carbonic anhydrase 3, sterol regulatory element-binding protein 1, and Spot 14.

Nrf2-dependent, D3T-inducible Genes in keap1-nrf2 Double Knockout Mice-- Gene expression profiles were examined in the liver of keap1-nrf2 double knockout mice to determine whether loss of Keap1, the putative sensor target for D3T, altered inducer-mediated gene expression patterns beyond that affected by loss of Nrf2. Overall, 94% of genes induced by D3T in wild-type mice were not increased in keap1-nrf2 double knockout mice. Among the 231 genes classified as Nrf2-dependent, D3T-inducible from the wild-type to nrf2-disrupted comparisons, just 15 genes including Cyp4a14, Cyp2c39, mEH, and GST theta2 were increased in these double knockout mice (Table I). Since the keap1-nrf2 double knockout mice behave like nrf2 single knockout mice upon D3T treatment, it appears that most D3T-inducible genes are induced through the Keap1-Nrf2 complex.

Nrf2-dependent, D3T-repressed Genes-- Shown in Table II, D3T treatment repressed the expression of 31 genes at 24 h and 30 largely overlapping genes following dietary administration of D3T in a Nrf2-dependent manner. Several cytochrome P450 isozymes (Cyp4a12, Cyp2d9, Cyp7a1, and Cyp8b1) were repressed by D3T treatment although no phase 2 genes were so affected. Expression levels of carbonic anhydrase 3 and 14 were reduced by half. An intriguing category of genes repressed by D3T treatment was the cholesterol/lipid biosynthesis metabolism-related genes. The transcription factor, sterol regulatory element-binding protein-1, which regulates cholesterol and lipid biosynthesis, was repressed 50-70% in the single gavage and dietary treatment groups.

                              
View this table:
[in this window]
[in a new window]
 
Table II
Hepatic Keap1-Nrf2 complex-dependent, D3T-repressed genes

RT-PCR Confirmation of Microarray Results-- Induction and repression of genes by D3T treatment and the effects of nrf2 genotype on the microarray results were confirmed using RT-PCR analysis in the 24-h treatment group (Fig. 2A). Gene-specific PCR primers were designed and RT-PCR reactions carried out. Gel images were quantified using densitometry and normalized using amplified albumin gene intensities. Expression levels of 22 genes among the D3T-inducible genes were analyzed by RT-PCR. NQO1 transcript levels were increased 10-fold by D3T in wild-type mice compared with vehicle-treated counterparts and were not increased in nrf2-disrupted mice in the RT-PCR analysis. Microarray analysis had indicated that NQO1 was increased 6-fold by D3T in the wild-type mice only. Transcripts levels for mEH, GST Mu3, and Hsp40 homolog showed 4.6-, 3.7-, and 5.5-fold increases, respectively, by D3T treatment in wild-type mice according to RT-PCR analysis. Corresponding fold changes of these genes in the microarray analysis were 4.0, 4.6, and 4.1, respectively. Expression of proteasome subunits were also confirmed and showed comparable induction by RT-PCR analysis. Representative repressed genes (sterol regulatory element-binding protein-1 and carbonic anhydrase 3) detected in the microarray analysis also showed reduced mRNA levels by RT-PCR analysis.


View larger version (32K):
[in this window]
[in a new window]
 
Fig. 2.   RT-PCR analyses. A, RT-PCR analysis of Nrf2-dependent, D3T-inducible genes in mouse liver. Gene-specific PCR primers were used to detect mRNA levels of each gene. PCR amplification was carried out with serially diluted cDNAs from vehicle (Veh.) or D3T-treated wild-type (WT) and nrf2-knockout mice (nrf2-KO). Albumin levels were used for normalization. B, RT-PCR analyses for NQO1, mEH, GST Mu3, p58, proteasome subunit beta type 3, and 26 S proteasome-associated pad1 homolog were performed with cDNA from vehicle (Veh.) or D3T-treated wild-type (WT) and keap1-nrf2 double knockout mice (keap1-nrf2 KO).

RT-PCR analysis was also performed to confirm the results from the microarray analysis in the double knockout mice (Fig. 2B). Expression levels of NQO1, GST Mu3, p58, proteasome subunits beta type 3 and 26 S proteasome-associated pad1 homolog were not increased by D3T treatment in the liver of keap1-nrf2 double knockout mice, while, mEH was slightly increased.

    DISCUSSION
TOP
ABSTRACT
INTRODUCTION
EXPERIMENTAL PROCEDURES
RESULTS
DISCUSSION
REFERENCES

Nrf2 is a critical transcription factor for mediating amplification of the mammalian defense system against environmental stresses. To identify the genes regulated by Nrf2, global gene expression patterns were examined in mouse liver following D3T treatment of nrf2 single and keap1-nrf2 double knockout mice using oligonucleotide microarray analysis. For evaluation of the microarray data, altered gene expression levels having a CV of less than 1 and significant differences in at least 7 of 9 possible comparisons were selected so as to filter out false positive responses. Filtering out false positive genes enabled the detection of induced genes with average fold changes of 1.2 and greater and of repressed genes under 0.7. The fold changes confirmed by RT-PCR consistently reflected the average fold changes observed in the microarray analyses. For example, induction of genes showing less than 2-fold increases by microarray analysis (eIF3p-44 and Sec61) exhibited similar fold changes by RT-PCR analysis. These results indicate that filtering genes by considering reproducibility and biological variance in individual responses provided an effective means to identify genes affected to small, but potentially biologically important degrees.

Genes increased or decreased by D3T treatment in wild-type but not nrf2-disrupted animals were classified as Nrf2-dependent genes. While, some genes increased in both genotypes, they showed higher fold induction in wild-type mice compared with nrf2-disrupted mice. These included heat shock protein 86 kDa 1 (AFCD3T/AFCVehicle; 4.11/1.39), Cyp2a4 (37.01/2.16), Cyp3a13 (3.41/1.73), UDP-glucose dehydrogenase (24 h after treatment, 7.52/1.67), and lipocalin (15.46/2.87). Therefore, expression of these genes by D3T treatment was affected by the nrf2-genotype. Collectively, there were 231 Nrf2-dependent, D3T-induced genes detected 24 h after treatment, representing 79% of total D3T-induced genes in wild-type mice. Similar percentages of Nrf2-dependent genes among all D3T-induced genes were observed at 6 h after a single treatment and with dietary exposure to D3T for 1 week. Thus, the large majority of D3T-induced genes appeared to be regulated through the Nrf2 pathway in mouse liver; 50% of D3T-repressed genes were also Nrf2-dependent genes.

The products of genes induced by D3T through a Nrf2-dependent pathway can be classified as xenobiotic metabolizing enzymes, antioxidants, general enzymes, chaperone/stress response proteins, ubiquitin/proteasome systems, protein trafficking proteins, acute response/immunity proteins, cytoskeletal organization proteins, membrane transport proteins, signaling molecules, transcription-related genes, and RNA processing and translation-related factors. This profile of changes in gene expression evoked by D3T is similar to that of the cellular response to oxidative stress in mammalian cells (36). The gene categories most affected by nrf2 genotype were xenobiotic metabolizing genes, antioxidants and proteasome subunits. D3T treatment induced genes belonging to other categories such as cytoskeletal organization proteins, and acute response genes in nrf2-disrupted mice; however, very few xenobiotic metabolizing genes and proteasome subunits were induced in the nrf2-knockout mice. These Nrf2-specific inducible gene categories may be critical contributors to the cytoprotective and chemopreventive effects of dithiolethiones.

Induction of xenobiotic detoxifying genes and antioxidants by dithiolethiones in a Nrf2-dependent manner is consistent with our previous results. Seven isozymes of GST and glutathione-related genes (glutathione reductase and gamma -glutamylcysteine ligases) were induced as expected and expression levels of NQO1, mEH, and UDP-glucuronyl transferase 2 family, member 5 were increased by D3T only in wild-type mice. Novel xenobiotic-metabolizing genes regulated by Nrf2 pathway are carbonyl reductase, flavin-containing monooxygenase 1, aldo-keto reductase 1A1, and aldehyde oxidase. Induction of carbonyl reductase, aldo-keto reductase 1A1, and aldehyde oxidase may contribute to the protective effects of dithiolethiones by detoxifying cytotoxic reactive carbonyl compounds (37-39). Induction of the antioxidant genes, peroxiredoxin, thioredoxin reductase, and thioredoxin restore normal cellular redox status under situations of oxidative stress (40). Other reducing pathways, namely, NADPH-generating enzymes such as UDP-glucose dehydrogenase and malate oxidoreductase were increased by D3T, again in an Nrf2-dependent pathway.

A previously unrecognized category of Nrf2-dependent, D3T-inducible genes is the ubiquitin/proteasome system. Twenty-six genes from this system were elevated by D3T treatment in wild-type mice. As depicted in Fig. 3, the genechip used in this study contains 39 genes related to proteasomes, and 67% of these proteasome genes were significantly induced by D3T. A substantial shift in the distribution of expression ratios was seen in this class of genes relative to all of the genes on the array (p < 0.001). By comparison, the expression distribution of glutathione synthetic and utilizing genes, reflecting a large number of known Nrf2-ARE regulated genes, was trimodal whereas that of 38 cytochrome P450 genes indicated an overall lack of responsiveness of these genes to D3T. A small subset of highly inducible P450s is apparent, however. Among the induced 26 S proteasome subunits, six subunits showed induction in both genotypes; however, induced fold changes were considerably higher in wild-type mice than nrf2-knockout mice. D3T treatment also induced heat shock proteins. These chaperone proteins together with the ubiquitin/proteasome system play an essential role in response to stress (41, 42). Accumulation of unfolded polypeptides following oxidative stress can disturb normal cellular functions and trigger apoptosis through the JNK pathway (42, 43). Heat shock proteins including Hsp40, Hsp70, Hsp90, and small Hsp (alpha -cystallin and Hsp27) function as molecular chaperones, and bind selectively to denatured or partially unfolded polypeptides (44). Many studies indicate that accumulation of unfolded or misfolded proteins induces a response in which chaperone genes and components of the ubiquitin pathway are up-regulated to repair damaged proteins in mammalian cells (45, 46). These induced heat shock proteins serve protective roles in neurons and expression of Hsp70 and Hsp40 protected neuronal cells from apoptosis by abnormal protein accumulation (47, 48). Increased expression of 26 S proteasome subunits by oxidative stress has not been reported in mammalian systems, although, proteasome subunits can be induced by interferon-gamma in human cell lines during viral infection (49). The expression of 26 S proteasome subunits is coordinately regulated by the transcription factor Rpn4p in Saccharomyces cerevisiae, and their expression can be affected by stress (50, 51). However, there does not appear to be any homology between Rpn4p and Nrf2. The increased expression of proteasome subunits by D3T in the presence of Nrf2 indicates that Nrf2 may be a common transcriptional regulator for the expression of 26 S proteasomes upon stress conditions in mammalian cells. Recently, Sekhar et al. (52) have shown that Nrf2 is degraded by the ubiquitin-proteasome pathway and that Nrf2 can be co-immunoprecipitated with ubiquitin polymer. Induction of chaperone and proteasome genes by D3T may facilitate the repair or degradation of damaged or oxidized proteins as an additional means to respond to oxidative and electrophilic stresses (Fig. 4).


View larger version (25K):
[in this window]
[in a new window]
 
Fig. 3.   Distribution profiles of genes as a function of their expression ratios. Ratios are those for D3T-treated over vehicle-treated wild-type mice. The distributions are for total genes (n = 12,487), cytochrome P450s (n = 38), glutathione (GSH)-related genes (n = 31), and proteasomes (n = 39). The numbers indicate the number of genes in each functional category contained on the mouse genechip. The distribution of D3T-induced proteasome genes is significantly different from the distribution for total genes (p < 0.001, Mann-Whitney rank sum test).


View larger version (31K):
[in this window]
[in a new window]
 
Fig. 4.   A model for Keap1-Nrf2 complex dependent genes involved in cytoprotection against environmental stresses. Dithiolethiones may react with sulfhydryl groups of Keap1 and liberate Nrf2 to translocate into the nucleus. Nrf2 transactivates the ARE with other bZIP proteins to increase the expression of detoxifying and antioxidative enzymes, NADPH-generating enzymes, glutathione biosynthesis enzymes, chaperone proteins, and ubiquitin-proteasomes. These gene products promote protection against electrophiles and reactive oxygen species from environmental and endogenous sources by detoxifying these reactive intermediates and enhancing the repair and degradation of damaged proteins.

Gene expression patterns were studied at different times. Many xenobiotic metabolizing genes were increased rapidly and persistently following a single dose of D3T, while genes related to restoration of damaged biomolecules such as the ubiquitination-proteasome system were principally elevated at 24 h. General enzymes, protein trafficking proteins and cytoskeletal components were increased early and only transiently after treatment. Relatively few gene transcripts remained elevated after feeding a D3T-containing diet for 7 days. The phase 2 detoxifying genes and proteasomes, with the exception of mEH, did not maintain their elevated transcript levels during continuous dietary exposure to D3T. These results indicate that repetitive administration of an Nrf2 activator may induce a cellular adaptation to the continuous signaling of gene transcription. This refractory transcriptional response to repetitive administration of dithiolethiones is consistent with our earlier observations (53). Comparison of the levels of GST subunit proteins and catalytic activity, mRNA levels and transcriptional rates of gst a1 showed that transcript level and transcription rate returned to near basal levels after 3 days of feeding the dithiolethione oltipraz, while, protein levels and activity remained maximally elevated for at least 1 week with feeding. Presumably altered post-translational modifications to the induced proteins contribute to the sustained GST activity. Thimmulappa et al. (19) has shown Nrf2-dependent, sulforaphane-inducible genes in the small intestine of mice following 7 daily doses of sulforaphane by gavage. They identified 26 genes as Nrf2-dependent, sulforaphane-inducible genes, and this gene cluster included xenobiotic-metabolizing enzymes, glutathione biosynthesis enzymes, and NADPH-generating enzymes. Whereas, these results are in accord with our present findings, they indicate that dosing schedule and sampling time influence the dynamics of gene expression. Therefore, analysis of gene expression patterns 24 h after single administration of inducer appears optimal for identifying Keap1Nrf2-dependent pathways. Companion proteomic analyses, however, are likely to reveal additional contributors to inducible cytoprotective responses.

The cytoplasmic repressor of Nrf2, Keap1 has been proposed to be a cellular sensor for protection against oxidative and electrophilic stresses. Direct interaction of enzyme inducers with cysteine residues of Keap1 has been demonstrated recently using the irreversible sulfhydryl reactant dexamethasone mesylate (16). Several cysteines in the intervening region of Keap1 are modified by this inducer. The chemical basis of dithiolethione activation of Nrf2 signaling is not clear as yet; however, dithiolethiones are sulfhydryl reactive compounds (54, 55) and stimulate the translocation of Nrf2 into the nucleus (8, 11). The critical role of Keap1 in the regulation of Nrf2 function has been established in studies using keap1-knockout mice. keap1-disrupted mice are not viable after 3 weeks of age; however, expression levels of phase 2 detoxifying genes such as GST and NQO1 are much higher in neonatal keap1-disrupted animals than in age-matched wild-type mice (33). The lethal phenotype of the keap1-disrupted mouse can be rescued by conjoint disruption of Nrf2 expression. Gene expression profiles in keap1-nrf2 double knockout mice following D3T treatment clearly show that most D3T-inducible genes are under the control of the Keap1-Nrf2 complex. Most (94%) of genes induced by D3T in wild-type mice were not increased in these double knockout mice. Fifteen genes classified as Nrf2-dependent, D3T-induced genes in the wild-type to nrf2-disrupted mice comparisons were induced in these keap1-nrf2 double knockout mice. These genes included mEH and GST theta 2, suggesting a partial recovery of alternative regulatory pathway(s) to express these protective genes in the double knockout mice.

The Keap1-Nrf2 complex appears to be the primary molecular target of chemopreventive agents such as the dithiolethiones. D3T increases the expression of a broad range of genes in a Keap1-Nrf2-dependent manner that act to directly detoxify toxins as well as generate essential cofactors such as glutathione and reducing equivalents. Induction of genes involved in the recognition and repair of damaged proteins expands the role of nrf2-dependent genes beyond primary control of electrophilic and oxidative stresses into secondary protective actions. Because dithiolethiones lose protective efficacy in nrf2-disrupted mice, these D3T-inducible, Keap1-Nrf2-dependent genes are likely to be central contributors to cell survival.

    ACKNOWLEDGEMENTS

We thank Kim Mai of the Microarray Core of the NIEHS Center in Urban Environmental Health for assistance with the microarray analyses. We also thank Patrick Dolan for maintaining and genotyping mice.

    FOOTNOTES

* This work was supported by Grants CA39416 and CA94076 from the National Institutes of Health and Center Grant ES 03819.The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

To whom correspondence should be addressed: Dept. of Environmental Health Sciences, Johns Hopkins University Bloomberg School of Public Health, 615 N. Wolfe St., Baltimore, MD 21205. Tel.: 410-955-4712; Fax: 410-955-0116; E-mail: tkensler@jhsph.edu.

Published, JBC Papers in Press, December 27, 2002, DOI 10.1074/jbc.M211898200

    ABBREVIATIONS

The abbreviations used are: ARE, antioxidant response element; D3T, 3H-1,2-dithiole-3-thione; GST, glutathione S-transferase; NQO1, NAD(P)H:quinone oxidoreductase; mEH, microsomal epoxide hydrolase; NFkappa B, nuclear factor kappa B; Hsp, heat shock protein; CV, coefficients of variation; AFC, average fold change; RT-PCR, reverse transcriptase-polymerase chain reaction.

    REFERENCES
TOP
ABSTRACT
INTRODUCTION
EXPERIMENTAL PROCEDURES
RESULTS
DISCUSSION
REFERENCES

1. Kensler, T. W. (1997) Environ. Health Perspect. 105, 965-970[Medline] [Order article via Infotrieve]
2. Hayes, J. D., and McMahon, M. (2000) Cancer Lett. 174, 103-113
3. Wang, J. S., Shen, X., He, X., Zhu, Y. R., Zhang, B. C., Wang, J. B., Qian, G. S., Kuang, S. Y., Zarba, A., Egner, P. A., Jacobson, L. P., Munoz, A., Helzlsouer, K. J., Groopman, J. D., and Kensler, T. W. (1999) J. Natl. Cancer Inst. 91, 347-354[Abstract/Free Full Text]
4. Lam, S., MacAulay, C., Le Riche, J. C., Dyachkova, Y., Coldman, A., Guillaud, M., Hawk, E., Christen, M. O., and Gazdar, A. F. (2002) J. Natl. Cancer Inst. 94, 1001-1009[Abstract/Free Full Text]
5. Rushmore, T. H., King, R. G., Paulson, K. E., and Pickett, C. B. (1990) Proc. Natl. Acad. Sci. U. S. A. 87, 3826-3830[Abstract]
6. Favreau, L., and Pickett, C. B. (1995) J. Biol. Chem. 270, 24468-24474[Abstract/Free Full Text]
7. Moinova, H. R., and Mulcahy, R. T. (1998) J. Biol. Chem. 273, 14683-14689[Abstract/Free Full Text]
8. Kwak, M-K., Itoh, K., Yamamoto, M., and Kensler, T. W. (2002) Mol. Cell. Biol. 22, 2883-2892[Abstract/Free Full Text]
9. Venugopal, R., and Jaiswal, A. K. (1996) Proc. Natl. Acad. Sci. U. S. A. 93, 14960-14965[Abstract/Free Full Text]
10. Itoh, K., Chiba, T., Takahashi, S., Ishii, T., Igarashi, K., Katoh, Y., Oyake, T., Hayashi, N., Satoh, K., Iatayama, I., Yamamoto, M., and Nabeshima, Y. (1997) Biochem. Biophys. Res. Commun. 236, 313-322[CrossRef][Medline] [Order article via Infotrieve]
11. Kwak, M-K., Itoh, K., Yamamoto, M., Sutter, T. R., and Kensler, T. W. (2001) Mol. Med. 7, 135-145[Medline] [Order article via Infotrieve]
12. Itoh, K., Wakabayashi, N., Katoh, Y., Ishii, T., Igarashi, K., Engel, J. D., and Yamamoto, M. (1999) Genes Dev. 13, 76-86[Abstract/Free Full Text]
13. Dhakshinamoorthy, S., and Jaiswal, A. K. (2001) Oncogene 20, 3906-3917[CrossRef][Medline] [Order article via Infotrieve]
14. Zipper, L. M., and Mulcahy, R. T. (2002) J. Biol. Chem. 277, 36544-36552[Abstract/Free Full Text]
15. Sekhar, K. R., Spitz, D. R., Harris, S., Nguyen, T. T., Meredith, M. J., Holt, J. T., Guis, D., Marnett, L. J., Summar, M. L., and Freeman, M. L. (2002) Free Radic. Biol. Med. 32, 650-662[CrossRef][Medline] [Order article via Infotrieve]
16. Dinkova-Kostova, A. T., Holtzclaw, W. D., Cole, R. N., Itoh, K., Wakabayashi, N., Katoh, Y., Yamamoto, M., and Talalay, P. (2002) Proc. Natl. Acad. Sci. U. S. A. 99, 11908-11913[Abstract/Free Full Text]
17. MacMahon, M., Itoh, K., Yamamoto, M., Chanas, S. A., Henderson, C. J., McLellan, L. I., Wolf, C. R., Cavin, C., and Hayes, J. D. (2001) Cancer Res. 61, 3299-3307[Abstract/Free Full Text]
18. Primiano, T., Gastel, J. A., Kensler, T. W., and Sutter, T. R. (1996) Carcinogenesis 17, 2297-2303[Abstract]
19. Thimmulappa, R. K., Mai, K. H., Srisuma, S., Kensler, T. W., Yamamoto, M., and Biswal, S. (2002) Cancer Res. 62, 5196-5203[Abstract/Free Full Text]
20. Lee, J. M., Hanson, J. M., Chu, W. A., and Johnson, J. A. (2001) J. Biol. Chem. 276, 20011-20016[Abstract/Free Full Text]
21. Enomoto, A., Itoh, K., Nagayoshi, E., Haruta, J., Kimura, T., O'Conner, T., Harada, T., and Yamamoto, M. (2001) Toxicol. Sci. 59, 169-177[Abstract/Free Full Text]
22. Chan, K., and Kan, Y. W. (1999) Proc. Natl. Acad. Sci. U. S. A. 96, 12731-12736[Abstract/Free Full Text]
23. Chan, K., Han, X. D., and Kan, Y. W. (2001) Proc. Natl. Acad. Sci. U. S. A. 98, 4611-4616[Abstract/Free Full Text]
24. Cho, H-Y., Jedlicka, A. E., Sekhar, M. S., Reddy, P. M., Zhang, L. Y., Kensler, T. W., Yamamoto, M., and Kleeberger, S. R. (2001) Am. J. Respir. Cell Mol. Biol. 26, 42-51[Abstract/Free Full Text]
25. Kwak, M-K., Egner, P. A., Dolan, P. M., Ramos-Gomez, M., Groopman, J. D., Itoh, K., Yamamoto, M., and Kensler, T. W. (2001) Mutat. Res. 480-481, 305-315
26. Aoki, Y., Sato, H., Nishimura, N., Takahashi, S., Itoh, K., and Yamamoto, M. (2001) Toxicol. Appl. Pharmacol. 173, 154-160[CrossRef][Medline] [Order article via Infotrieve]
27. Ramos-Gomez, M., Dolan, P. M., Itoh, K., Yamamoto, M., Talalay, P., and Kensler, T. W. (2003) Carcinogenesis, in press
28. Ramos-Gomez, M., Kwak, M-K., Dolan, P. M., Itoh, K., Yamamoto, M., Talalay, P., and Kensler, T. W. (2001) Proc. Natl. Acad. Sci. U. S. A. 98, 3410-3415[Abstract/Free Full Text]
29. Wattenberg, L. W., and Bueding, E. (1986) Carcinogenesis 7, 1379-1388[Abstract]
30. Rao, C. V., Rivenson, A., Katiwalla, M., Kelloff, G. J., and Reddy, B. S. (1993) Cancer Res. 53, 2502-2506[Abstract]
31. Kensler, T. W., Egner, P. A., Dolan, P., Groopman, J. D., and Roebuck, B. D. (1987) Cancer Res. 47, 4271-4277[Abstract]
32. Kensler, T. W., Groopman, J. D., Eaton, D. L., Curphey, T. J., and Roebuck, B. D. (1992) Carcinogenesis 13, 95-100[Abstract]
33. Yamamoto, M., Wakabayashi, N., Ishii, T., Kobayashi, M., Kato, Y., and Itoh, K. (2002) Proceedings of the 14th International Symposium on Microsomes and Drug Oxidation , p. 48, Sapporo, Japan
34. Li, J., and Johnson, J. A. (2002) Physiol. Genomics 9, 137-144[Abstract/Free Full Text]
35. Bouton, C. M., and Pevsner, J. (2002) Bioinformatics 18, 323-324[Abstract/Free Full Text]
36. Halliwell, B., and Gutteridge, J. M. C. (1999) Free Radicals in Biology and Medicine , Oxford University Press, Inc., New York
37. Forrest, G. L., and Gonzalez, B. (2000) Chem. Biol. Interact. 129, 21-40[CrossRef][Medline] [Order article via Infotrieve]
38. Oppermann, U. C., and Maser, E. (2000) Toxicology 144, 71-81[CrossRef][Medline] [Order article via Infotrieve]
39. Yoshihara, S., and Ohta, S. (1998) Arch. Biochem. Biophys. 360, 93-98[CrossRef][Medline] [Order article via Infotrieve]
40. Tanaka, T., Nakamura, H., Nishiyama, A., Hosoi, F., Masutani, H., Wada, H., and Yodoi, J. (2000) Free Radic. Res. 33, 851-855[Medline] [Order article via Infotrieve]
41. Gabai, V. L., Meriin, A. B., Yaglom, J. A., Volloch, V. Z., and Sherman, M. Y. (1998) FEBS Lett. 438, 1-4[CrossRef][Medline] [Order article via Infotrieve]
42. Gabai, V., L., Meriin, A. B., Mosser, D. D., Caron, A. W., Rits, S., Shifrin, V. I., and Sherman, M. Y. (1997) J. Biol. Chem. 272, 18033-18037[Abstract/Free Full Text]
43. Samali, A., and Orrenius, S. (1998) Cell Stress. Chaperones 3, 228-236[CrossRef][Medline] [Order article via Infotrieve]
44. Johnson, J. L., and Craig, E. A. (1997) Cell 90, 201-204[Medline] [Order article via Infotrieve]
45. Watt, R., and Piper, P. W. (1997) Mol. General Genet. 253, 439-447[CrossRef]
46. Hershko, A., and Ciechanover, A. (1998) Annu. Rev. Biochem. 67, 425-479[CrossRef][Medline] [Order article via Infotrieve]
47. Kobayashi, Y., Kume, A., Li, M., Doyu, M., Hata, M., Ohtsuka, K., and Sobue, G. (2000) J. Biol. Chem. 275, 8772-8778[Abstract/Free Full Text]
48. Satyal, S. H., Schmidt, E., Kitagawa, K., Sondheimer, N., Lindquist, S., Kramer, J. M., and Morimoto, R. I. (2000) Proc. Natl. Acad. Sci. U. S. A. 97, 5750-5755[Abstract/Free Full Text]
49. Foss, G. S., Larsen, F., Solheim, J., and Prydz, H. (1998) Biochim. Biophys. Acta 1402, 17-28[CrossRef][Medline] [Order article via Infotrieve]
50. Xie, Y., and Varshavsky, A. (2001) Proc. Natl. Acad. Sci. U. S. A. 98, 3056-3061[Abstract/Free Full Text]
51. Ng, D. T., Spear, E. D., and Walter, P. (2000) J. Cell Biol. 150, 77-88[Abstract/Free Full Text]
52. Sekhar, K. R., Yan, X. X., and Freeman, M. L. (2002) Oncogene 21, 6829-6834[CrossRef][Medline] [Order article via Infotrieve]
53. Davidson, N. E., Egner, P. A., and Kensler, T. W. (1990) Cancer Res. 50, 2251-2255[Abstract]
54. Levron, B., Burgot, G., and Burgot, J. L. (2000) Arch. Biochem. Biophys. 382, 189-194[CrossRef][Medline] [Order article via Infotrieve]
55. Carey, K. A., Kensler, T. W., and Fishbein, J. C. (2001) Chem. Res. Toxicol. 14, 939-945[CrossRef][Medline] [Order article via Infotrieve]


Copyright © 2003 by The American Society for Biochemistry and Molecular Biology, Inc.