Abstract
Monkeys are widely used as a primate model to study drug metabolism because they generally show a metabolic pattern similar to humans. However, the paucity of information on cytochrome P450 (P450) genes has hampered a deep understanding of drug metabolism in the monkey. In this study, we report identification of the CYP2C76 cDNA newly identified in cynomolgus monkey and characterization of this CYP2C along with cynomolgus CYP2C20, CYP2C43, and CYP2C75. The CYP2C76 cDNA contains the open reading frame encoding a protein of 489 amino acids that are only approximately 80% identical to any human or monkey P450 cDNAs. Gene and protein expression of CYP2C76 was confirmed in the liver of cynomolgus and rhesus monkeys but not in humans or the great apes. Moreover, CYP2C76 is located at the end of the CYP2C gene cluster in the monkey genome, the region of which corresponds to the intergenic region adjacent to the CYP2C cluster in the human genome, strongly indicating that this gene does not have the ortholog in humans. Among the four CYP2C genes expressing predominantly in the liver, the expression level of CYP2C76 was the greatest, suggesting that CYP2C76 is a major CYP2C in the monkey liver. Assays for the capacity of CYP2C76 to metabolize drugs using several substrates typical for human CYP2Cs revealed that CYP2C76 showed unique metabolic activity. These results suggest that CYP2C76 contributes to overall drug-metabolizing activity in the monkey liver and might account for species difference occasionally seen in drug metabolism between monkeys and humans.
Cytochrome P450s (P450s) are one of the most important drug-metabolizing enzymes and form a superfamily consisting of a large number of subfamilies (Nelson et al., 1996, 2004). The cDNA sequences encoding P450s have been reported for many species of not only mammals but also birds, insects, plants, bacteria, and others (see http://drnelson.utmem.edu/CytochromeP450.html). In humans, 57 functional genes have been identified to date (Nelson et al., 2004). The human CYP2C subfamily, comprising CYP2C8, CYP2C9, CYP2C18, and CYP2C19, is essential in metabolizing approximately 20% of all prescribed drugs, including tolbutamide, phenytoin, warfarin, and ibuprofen (Goldstein, 2001). The CYP2C subfamily consists of multiple members in each mammalian species, including 15 in mice, 12 in rats, and 9 in rabbits (for the latest information, see http://drnelson.utmem.edu/CytochromeP450.html). Between humans and rodents, the number of the subfamily members is different, and none of the CYP2Cs seems to show a clear orthologous relationship between the two species, suggesting that the data from rodents must be cautiously interpreted and extrapolated to humans (Nelson et al., 2004).
For monkeys, which generally mean Old or New World monkeys, three CYP2C cDNAs have been identified in the macaque and cynomolgus (Macaca fascicularis) and rhesus (Macaca mulatta) monkeys. Two sequences have been published including cynomolgus CYP2C20 (Komori et al., 1992) and rhesus CYP2C43 (Matsunaga et al., 2002), whereas rhesus CYP2C75 has been reported to GenBank but unpublished. CYP2C20 shows ∼95% homology to human CYP2C8, whereas CYP2C43 and CYP2C75 have ∼95% identity to both human CYP2C9 and CYP2C19. Among these monkey CYP2Cs, only CYP2C43 has been analyzed for drug-metabolizing capacity using the recombinant protein, showing activities toward S-mephenytoin but not tolbutamide, similar to metabolic properties for human CYP2C19 (Matsunaga et al., 2002).
The monkey, the animal evolutionarily close to humans, shows not only phenotypic but also physiological similarities to humans in such biological circumstances as aging (Roth et al., 2004), reproduction (Bellino and Wise, 2003), and neurological disease, including Parkinson's disease (Takagi et al., 2005). This resemblance can be partly explained by a relatively high similarity (∼97%) of protein sequences generally seen between monkeys and humans (Magness et al., 2005). For preclinical trials during drug development, monkeys, especially macaques, are used as a large and nonrodent species to evaluate drug effects because they generally show pharmacokinetics more similar to humans than any nonprimate models. However, it has become apparent that monkeys are not always similar to humans in drug metabolism (Stevens et al., 1993; Sharer et al., 1995; Guengerich, 1997; Weaver et al., 1999; Bogaards et al., 2000; Narimatsu et al., 2000). We hypothesized that this species difference was due to the difference in the genetic components essential for drug metabolism such as P450s between monkeys and humans. In monkeys, information on the characteristics of drug-metabolizing enzymes is largely scarce, especially at molecular level, preventing a deep understanding of drug metabolism. Therefore, to test our hypothesis, we have identified a number of the cynomolgus monkey cDNA clones for P450s; the characterization of these clones is currently ongoing. Successful drug development requires accuracy in the extrapolation of drug metabolism and toxicity data from experimental animals to humans. Then characterization of these cynomolgus P450s should help to better understand drug metabolism in monkeys.
In this article, we report the isolation and characterization of cynomolgus CYP2C76, a novel P450 with a low homology to any human or monkey CYP2C cDNAs. RT-PCR and immunoblotting indicated the expression of mRNA and protein homologous to CYP2C76 in cynomolgus and rhesus monkeys but not in humans. Moreover, the genomic analysis indicated that the CYP2C76 gene was located at the end of the CYP2C cluster in the macaque genome, the location of which corresponded to an intergenic region in the human genome, suggesting that the P450 homologous to CYP2C76 does not exist in humans. The hepatic expression of CYP2C76 was higher than any other CYP2Cs analyzed, indicating that this P450 was a major CYP2C in the monkey liver. CYP2C76 protein was active in the metabolism of tolbutamide and testosterone. Because of its species specificity and functional importance, CYP2C76 might account for the species difference occasionally seen in drug metabolism between monkeys and humans.
Materials and Methods
Chemicals and Reagents. 4-Hydroxytolbutamide, 6α-hydroxypaclitaxel, 3-hydroxypaclitaxel, 4-hydroxy-S-mephenytoin, and 6β-hydroxytestosterone were purchased from Ultrafine Chemicals (Manchester, UK). Pooled hepatic microsomes from human subjects and male cynomolgus monkeys were both purchased from BD Gentest (Woburn, MA). Labeled [2-benzoyl ring-U-14C]paclitaxel (2.23 MBq/mg) was obtained from Sigma-Aldrich (St. Louis, MO), whereas [ring-U-14C]tolbutamide (2.26 GBq/mmol), S-[4-14C]mephenytoin (9.42 MBq/mg), and [4-14C]testosterone (2.11 GBq/mmol) were from GE Healthcare (Little Chalfont, Buckinghamshire, UK). The radiochemical purities of these 14C-labeled chemicals were >99%. Oligonucleotides were synthesized by Sigma-Genosys (Ishikari, Japan). All other reagents were purchased from Sigma-Aldrich unless otherwise specified.
Tissue Samples and RNA Extraction. Tissue samples were collected from individual monkeys, including six cynomolgus monkeys (three male and three female) and two male rhesus monkeys, which were kept under the established guidelines and standard procedures at Shin Nippon Biomedical Laboratories (Tokyo, Japan). The study was approved by the local ethics committee. The tissue samples from brain, lung, heart, liver, kidney, adrenal gland, small intestine, testis, ovary, and uterus were frozen in liquid nitrogen right after removal from animals to prevent potential RNA degradation. Orangutan and chimpanzee liver samples were kindly provided by GAIN (Great Ape Information Network, Japan). The frozen tissues were first ground with mortar and pestle and then processed in TRIzol (Invitrogen, Carlsbad, CA) with a Polytron homogenizer (Kinematica, Basel, Switzerland), followed by extraction of total RNA according to the manufacturer's instruction. After treatment with DNase I (Takara, Tokyo, Japan), the RNA was purified using GenElute Mammalian Total RNA Mini Kit (Sigma-Aldrich).
Cell Culture and RNA Extraction. COS1 cells from the American Type Culture Collection (Manassas, VA) was cultured as described previously (Saito et al., 2001). RNA extraction from the cell and the subsequent DNase I treatment were performed using RNeasy Mini Kit (QIAGEN, Valencia, CA) according to the manufacturer's protocol.
Cloning of CYP2C76 Homologous cDNA in Other Primate Species. RNA was extracted from liver (for rhesus monkey, orangutan, and chimpanzee) and COS1 cell (for African green monkey). The first-strand cDNA was generated in a mixture containing 1 μg of total RNA, oligo (dT) or random primers, and Moloney murine leukemia virus reverse transcriptase (Toyobo, Osaka, Japan) at 37°C for 1 h. The resultant cDNA was diluted 25-fold and used as a template for the subsequent PCR. For human cDNA, the liver cDNA available from BD Biosciences (San Jose, CA) was used. The amplification was carried out using KOD Plus DNA polymerase (Toyobo) according to the manufacturer's protocol with the MJ Research thermal cycler (MJ Research, Watertown, CA). PCR conditions include an initial denaturation at 95°C for 2 min and 30 cycles of 95°C for 20 s, 55°C for 20 s, and 72°C for 2 min, followed by a final extension at 72°C for 10 min. Among several different primer-pairs tested, a CYP2C76 homologous sequence for rhesus and African green monkeys was successfully amplified by PCR using the following primer pairs: mf27B9 (5rt2), 5′-CCCAGCAATGGATCTCTTCA-3′, and mf27B9 (3polyA2a), 5′-TGCCTAGACAGGTAGATAGGAGTG-3′, for the rhesus monkey; mf27B9 (5rt2) and mf27B9 (3ex4b), 5′-GAAAAGTGGGATCACAGGGA-3′, for the African green monkey. After the addition of 3′ A-overhangs, the PCR products were cloned into vectors using TOPO TA Cloning Kit (Invitrogen). The inserts were then sequenced using ABI Prism BigDye Terminator v3.0 Ready Reaction Cycle Sequencing Kit (Applied Biosystems, Foster City, CA), followed by electrophoresis with the ABI Prism 3730 DNA Analyzer (Applied Biosystems).
Sequence Analysis. Raw sequence data were imported into DNASIS Pro (Hitachi Software, Tokyo, Japan) for most sequence analyses. After vector sequences and regions of sequence with unacceptable quality were removed, the trimmed sequences were assembled to the full-length sequence. A homology search was conducted by the BLAST program (National Center for Biotechnology Information). Multiple alignment of amino acid sequences was performed with the ClustalW program; the resultant alignment was used to create a phylogenetic tree using the PHYLIP program by default parameter. The human and chimpanzee genome data were searched for the sequence homologous to CYP2C76 by BLAT search (UCSC Genome Bioinformatics). Likewise, the macaque genome data (Baylor College of Medicine Human Genome Sequencing Center) were used to identify and analyze the genome sequences corresponding to macaque CYP2Cs.
Amplification of CYP2C76 Introns. All the introns were amplified from cynomolgus monkey genomic DNA by PCR with 5 pmol each of forward and reverse gene-specific primers, 0.5 mM dNTPs, 2 mM MgCl2, and 1 unit of LA Taq polymerase (Takara) in a total volume of 20 μl. The primers used to amplify introns 1 to 8 were mf27B9 (5rt2), mf27B9 (5ex2a), mf27B9 (5ex3a), mf27B9 (5ex4a), mf27B9 (5ex5a), mf27B9 (5ex6a), mf27B9 (5qrt1), or mf27B9 (5gen1) as a forward primer and mf27B9 (3ex2a), mf27B9 (3ex3a), mf27B9 (3ex4a), mf27B9 (3ex5a), mf27B9 (3ex6a), mf27B9 (3ex7a), mf27B9 (3qrt1), or mf27B9 (3rt1) as a reverse primer, respectively. The nucleotide sequence for each primer is listed below. Thermal cycler conditions were as follows: 95°C for 2 min; 35 cycles of 95°C for 20 s, 55°C for 30 s, and 72°C for 5 min; and a final extension at 72°C for 20 min. After electrophoresis in 0.8% agarose gels, the PCR products were gel-purified, cloned into vectors using TOPO XL Cloning Kit (Invitrogen), and sequenced. Sequencing and sequence analysis were performed as described above to determine the entire sequence of each intron.
The sequences of oligonucleotide primers used are as follows: mf27B9 (5ex2a), 5′-GTATTTTCTGGCCGAGGGAG-3′; mf27B9 (5ex3a), 5′-CGGCGTTTCTCTCTCATGGT-3′; mf27B9 (5ex4a), 5′-GGGTTGTGTTCCCTGTAATGTC-3′; mf27B9 (5ex5a), 5′-CATCAGGAATCTCTGGACATC-3′; mf27B9 (5ex6a), 5′-CAGAGACAACAAGCACCACAA-3′; mf27B9 (5qrt1), 5′-CCCATGCAGTGACTCAAGAC-3′; mf27B9 (5gen1), 5′-GCCACTTCCTGGACGAAAG-3′; mf27B9 (3ex2a), 5′-CTCCCTCGGCCAGAAAATAC-3′; mf27B9 (3ex3a), 5′-ATGCTTCCACCAGACACAAG-3′; mf27B9 (3ex4a), 5′-ACAGGGAACACAACCCAGAA-3′; mf27B9 (3ex5a), 5′-CGAGGGTTATTGATGTCCAGAG-3′; mf27B9 (3ex6a), 5′-GGAGCATCAGTCCATATCTCATT-3′; mf27B9 (3ex7a), 5′-ATTGGTGGGGATGAGGTCAATA-3′; mf27B9 (3qrt1), 5′-AAGTGGCCAGGGTCAAACT-3′; and mf27B9 (3rt1), 5′-ACAGCCTTGCTCTGCAATC-3′.
Isolation and Analysis of Monkey BAC Clones. The BAC clone containing the CYP2C genes was isolated by screening a rhesus monkey BAC library (BACPAC, Oakland, CA) using the CYP2C75 and CYP2C76 cDNAs as probes, because the BAC library was not available for cynomolgus monkeys. Hybridization with the library filters was carried out as recommended by the manufacturer using the probes synthesized in the presence of [α-32P]dCTP (GE Healthcare) with the RadPrime DNA labeling system (Invitrogen). The identified BAC clones were obtained from the BACPAC. The BAC DNA was purified using DNA PhasePrep BAC DNA Kit (Sigma-Aldrich). To identify the CYP2C genes contained in each BAC DNA, the purified DNA was used as a template for the PCR with specific primers for 5′ or 3′ of each gene, 0.5 mM dNTPs, 2 mM MgCl2, and 1 unit of AmpliTaq Gold DNA polymerase (Applied Biosystems) in a total volume of 20 μl. PCR conditions were as follows: 95°C for 10 min; 30 cycles of 95°C for 20 s, 55°C for 20 s, and 72°C for 1 min; and final extension at 72°C for 10 min. The different primer pairs were initially designed at exons 1 and 9 of each gene, the location of which was determined by comparing each cynomolgus cDNA to the human CYP2C genes. For the genes highly homologous to CYP2C43 and CYP2C75, the designed primers did not show a gene-specific amplification pattern. Therefore, the gene-specific indels were identified in introns 1 and 8 of these genes by searching the macaque genome data and the primers recognizing these indels were designed. The sequences of primers used for the PCR are listed in Table 1. The amplification pattern was examined to determine an arrangement of the CYP2C genes in the genome as described previously (Gray et al., 1995). For the same purpose, the DNA was also used for the BAC end sequencing and for a restriction enzyme mapping with BamHI or EcoRI as recommended by the BACPAC.
Real-Time RT-PCR. The primers and TaqMan MGB probes specific for each gene (Table 2) were designed using Primer Express software (Applied Biosystems). The 5′ end of the probes was labeled with 5-carboxyfluorescein fluorescence reporter dye. Reverse transcriptase reaction was carried out using random primers as described above. A twenty-fifth volume of the reaction mixture was then used for the subsequent PCR that was carried out in a total volume of 25 μl using TaqMan Universal PCR Master Mix (Applied Biosystems) with the ABI Prism 7700 sequence detection system (Applied Biosystems) following the manufacturer's protocol. Final concentration of each primer set was 0.3 μM for CYP2C20 and CYP2C43, 0.9 μM for CYP2C75, and 0.1 μM for CYP2C76. The final concentration of the probes was 0.25 μM for all CYP2Cs. Thermal cycler conditions for all reactions were 2 min at 50°C and 10 min at 95°C, followed by 40 cycles of 15 s at 95°C and 1 min at 60°C. Standard curves were generated by serial 10-fold dilutions of a plasmid for the corresponding cDNA. The specificity of assays for all CYP2C genes was confirmed by sequencing a single DNA band with the expected size in agarose gels and by performing the highly efficient amplification of cDNA plasmid for the target gene over that of the other CYP2C genes. Relative expression level of each gene was normalized to the 18S ribosomal RNA level measured using a predeveloped kit available from Applied Biosystems. At least three amplifications were performed for each gene.
Heterologous Expression of Four Recombinant CYP2Cs in Escherichia coli. Protein expression of the four CYP2Cs was carried out as described previously (Iwata et al., 1998). To enhance protein expression, the eight residues of the N terminus were replaced with the corresponding ones of the modified bovine CYP17, MALLLAVF (Barnes et al., 1991) by amplifying the open reading frame of each cDNA with the primers listed in Table 3 using the KOD Plus DNA polymerase. The forward and reverse primers contained NdeI and XbaI sites, respectively, so that after restriction enzyme digestion, the PCR products can be easily cloned into the pCW vector (Barnes, 1996), in which the human reductase cDNA has been already accommodated (Iwata et al., 1998). The resultant construct for CYP2C20, CYP2C43, CYP2C75, or CYP2C76 was used to transform DH5α competent cells (Invitrogen). Thereafter, sequence and orientation of each insert were confirmed by sequencing. To express proteins, the bacteria cells grown overnight in Luria-Bertani broth were diluted 100-fold and cultured in the presence of 200 μg/ml ampicillin for 6 to 12 h at 30°C in the modified Terrific broth (Iwata et al., 1998) until the optical density at 600 nm reached approximately 0.6 to 0.8. Isopropyl-β-d-thiogalactoside was then added to the culture at a final concentration of 1.5 mM. After 16 to 20 h, the cultured cells were harvested and cell membrane fraction was prepared as described previously (Daigo et al., 2002). The content of each CYP2C protein in the membrane preparation was determined by Fe2+ ·CO versus Fe2+ difference spectra, according to the method described by Omura and Sato (1964) using the U-3000 spectrophotometer (Shimadzu, Kyoto, Japan). The concentration of NAPDHP450 reductase was also measured as reported previously (Iwata et al., 1998; Daigo et al., 2002).
Drug-Metabolizing Activity of the Partially Purified Preparations of Monkey CYP2Cs. All the recombinant CYP2C proteins were analyzed for their activities to metabolize drugs with four substrates prototypical for human P450s, including paclitaxel, tolbutamide, S-mephenytoin, and testosterone. To prepare the reaction mixture, [14C]paclitaxel (6 μM), [14C]tolbutamide (100 μM), [14C]S-mephenytoin (50 μM), or [14C]testosterone (50 μM) was preincubated in a 100 mM sodium phosphate buffer solution, pH 7.4, with hepatic microsomes (1 mg/ml) or the partially purified recombinant CYP2Cs (200 pmol/ml) at 37°C for 5 min. This was followed by the addition of a NADPH regenerating system containing 1.3 mM NADP+, 3.3 mM glucose 6-phosphate, 0.4 U/ml glucose-6-phosphate dehydrogenase, and 3.3 mM magnesium chloride in 100 mM sodium phosphate buffer, pH 7.4, so that the metabolic reaction was initiated at a final concentration of 1 mg of protein/ml or 200 pmol of P450/ml. After incubation at 37°C for 15 min (paclitaxel and testosterone), 45 min (S-mephenytoin), or 60 min (tolbutamide), the reaction was quenched by adding an equal volume of 100% methanol solution. For the recombinant CYP2Cs, the incubation was carried out for 30 min with all the tested substrates. The reaction-terminated samples were centrifuged, the aliquots of the supernatant were evaporated to dryness, and the residue was dissolved in 15 μl of methanol. The analysis of 6α-hydroxypaclitaxel and 3-hydroxypaclitaxel was performed as reported previously (Fujino et al., 2001). Aliquots (∼2 μl) of the supernatant were spotted onto TLC plates (Silicagel 60F254, 20 × 20 cm; Merck, Darmstadt, Germany) and developed with toluene-acetone-formic acid [60:39:1 (v/v/v)] to 12 cm in a horizontal TLC chamber that was saturated with solvent vapor. The analysis of 6β-hydroxytestosterone was carried out as follows: supernatant was applied to TLC plates and developed with dichloromethane-acetone [4:1 (v/v)] to 16 cm. The measurement of 4-hydroxytolbutamide was performed as reported previously (Ludwig et al., 1998). In brief, the supernatant was spotted and developed with toluene-acetone-formic acid [60:39:1 (v/v/v)] to 10 cm. The assay for S-mephenytoin 4-hydroxylase activity was also performed according to a previous report (Shimada et al., 1985). The spotted supernatant was developed with chloroform-methanol-28% ammonium [90:10:1 (v/v/v)] to 12 cm. The TLC plates were dried and placed in contact with a phosphor imaging plate for 12 h. The amounts of unchanged drug and metabolites were determined using the BAS-2500 (Fuji Photo Film Co., Tokyo, Japan). The radioactive metabolites were positively identified by a comparison of Rf values using authentic unlabeled standard.
Immunoblotting. Polyclonal antibodies raised to CYP2C76 were produced by NeoMPS (San Diego, CA) using specific peptides for this protein. In brief, a peptide specific for CYP2C76 (Fig. 1) was designed, synthesized, purified, coupled through the terminal cysteine thiol with a keyhole limpet hemocyanin, and used to immunize New Zealand white rabbits. The recombinant P450 proteins (1.0 pmol each) were run in 10% SDS polyacrylamide gels and transferred to Hybond-P filters (GE Healthcare). The filters were immunoblotted with the rabbit anti-CYP2C76 (1:250) and the donkey anti-rabbit IgG conjugated with horseradish peroxidase (SantaCruz Biotechnology, Santa Cruz, CA). To detect protein disulfide isomerase (PDI) as a loading control, rabbit anti-PDI (1:200, SantaCruz Biotechnology) was also used. A specific band was visualized using an ECL Western blotting detection reagent (GE Healthcare) according to the manufacturer's instructions.
Immunohistochemistry. The sections of cynomolgus monkey liver were immunostained with the anti-CYP2C76 antibody after the standard procedure. In brief, the primary antibodies were diluted 50-fold and applied to the sections at 4°C overnight. The bound antibodies were detected using the EnVision+ System (Dako North America, Inc., Carpinteria, CA) and liquid diaminobenzidine (Dako North America) according to the manufacturer's instruction. Slides were counterstained with Harris hematoxylin. As a negative control, rabbit preimmune serum was used instead of primary antibodies. To validate immunohistochemical specificity, the antibodies were preincubated at 4°C overnight with excess amount of the CYP2C76 specific peptide (0.05 mg/ml) and this mixture was used in place of primary antibodies for immunohistochemistry.
Results
Identification of Cynomolgus CYP2C cDNAs. CYP2C76, along with CYP2C20, CYP2C43, and CYP2C75, was originally identified as a unique cDNA clone by searching our in-house EST database that was established using a full-length cDNA library prepared from the cynomolgus monkey liver (Y. Uno, Y. Suzuki, Y. Sakamoto, H. Sano, K. Hashimoto, S. Sugano, and I. Inoue, unpublished data). Among these, the sequences newly identified in cynomolgus monkey, CYP2C43, CYP2C75, and CYP2C76, have been deposited to GenBank under accession numbers of DQ074806, DQ074805, and DQ074807, respectively. The CYP2C76 cDNA contained the open reading frame of 489 amino acids (Fig. 1). The deduced amino acid sequence showed primary sequence structures common to CYP2C molecules, including a highly hydrophobic N terminus, heme-binding region, and six potential substrate recognition sites (Gotoh, 1992). Blast analysis using the deduced amino acid sequences showed that CYP2C76 had only ∼71% identity to any human CYP2Cs, whereas CYP2C20, CYP2C43, and CYP2C75 were ∼92% homologous to human CYP2C (Table 4). This, together with a phylogenetic comparison of CYP2C amino acid sequences among mammalian species (Fig. 2), indicated possibilities that the CYP2C76 might be monkey-specific or that the human ortholog has not been isolated. To examine the latter possibility, we attempted to identify the sequence homologous to CYP2C76 in humans and other primate species by RT-PCR using gene-specific primer pairs. RNA samples used were from humans, the great apes (chimpanzee and orangutan), and Old World monkeys (cynomolgus, rhesus, and African green monkeys). The amplification was seen for Old World monkeys, and sequences of the PCR products were determined (data not shown). These CYP2C76 homologous sequences were ≥99% identical to each other, reflecting evolutionary closeness of these species. In contrast, the amplification with the human, chimpanzee, and orangutan samples showed no detectable bands in agarose gels (data not shown). Moreover, searching human and chimpanzee genome databases by BLAT (UCSC Genome Bioinformatics) showed no potential CYP2C genome sequence ≥90% homologous to the CYP2C76, raising the possibility that the CYP2C76 is monkey-specific.
Genomic Organization of the Monkey CYP2C Locus. To confirm the species specificity for CYP2C76, the location of the CYP2C76 in the genome was determined by analyzing rhesus monkey CYP2C BAC clones. The CYP2C-positive clone was used for PCR as a template with gene-specific primers that were assigned at the 5′ and 3′ ends of each cynomolgus CYP2C cDNA. During the course of this study, the genome sequence data of the rhesus monkey became available and was used to confirm that the designed primers could be used for the rhesus monkey. The analysis of the data identified the genome sequence highly homologous to human CYP2C18. Therefore, the primers specific for this CYP2C18-like gene were also designed and used for amplification with the BAC clones. The amplification pattern of each gene, together with the end-sequencing and the restriction enzyme mapping of the BAC clones, indicated that the five CYP2C genes together form a gene cluster in the monkey genome similar to humans (Fig. 3). Moreover, CYP2C76 was located at the end of the cluster, corresponding to the intergenic region adjacent to the CYP2C cluster in the human genome. These results strongly support the idea that CYP2C76 is expressed in monkeys but not in humans.
Gene Structure of CYP2C76. To determine the gene structure of CYP2C76, long PCR amplification was performed with the genomic DNA of cynomolgus monkey as a template. Gene-specific primers were designed on each exon to amplify each intron. The sequences of the PCR products were determined as described above and assembled into the entire intron. The CYP2C76 gene spanned approximately 19.6 kb and contained nine exons, as has been described for all human CYP2C genes. Sizes of exons and introns ranged from 142 to 693 base pairs and from 937 to 4307 base pairs, respectively (Table 5). All exons were flanked by GU and AG dinucleotides consistent with the consensus sequences for splice junctions in eukaryotic genes, with the exception of the 5′ splice site for intron 8, where GU was replaced by GC.
Tissue Distribution of Gene Expression. To analyze the expression of cynomolgus CYP2C20, CYP2C43, CYP2C75, and CYP2C76, real-time RT-PCR was performed with gene-specific primers and TaqMan MGB probes using RNAs prepared from brain, lung, heart, liver, kidney, adrenal gland, small intestine, testis, ovary, and uterus. All the four CYP2C genes were expressed predominantly in the liver with some extrahepatic expression (Fig. 4). Among these CYP2Cs, the expression level of CYP2C76 was the highest, indicating that CYP2C76 is a major CYP2C in the monkey liver.
Drug-Metabolizing Activities of Monkey Recombinant CYP2Cs. The activity of monkey CYP2Cs to metabolize drugs was characterized by incubating partially purified recombinant CYP2Cs with NADPH-regenerating system in the presence of the radiolabeled substrates typical for human CYP2Cs, paclitaxel, tolbutamide, S-mephenytoin, and testosterone. The results showed that CYP2C20 was involved only in the metabolism of paclitaxel among the four substrates examined, similar to the metabolic properties observed for human CYP2C8 (Fig. 5A). Tolbutamide was metabolized by CYP2C75 and CYP2C76 (Fig. 5B), whereas S-mephenytoin was only weakly metabolized by CYP2C43 and CYP2C75 (Fig. 5C). Testosterone was efficiently metabolized by all the four CYP2Cs except for CYP2C20, in which the metabolites generated were different (Fig. 5D). It is noteworthy that all three P450s capable of metabolizing testosterone generated one common but unknown metabolite. In addition, the 6β-hydroxylation of testosterone and the 3-hydroxylation of paclitaxel, both of which are mediated by CYP3A4 in humans, were not observed in the presence of CYP2Cs. Our results indicate that each monkey CYP2C has the characteristic substrate specificity.
Immunoblotting and Immunohistochemistry. The peptide specific for CYP2C76 was synthesized and used to raise anti-CYP2C76 antibodies. To investigate the specificity of the antibodies, immunoblotting was performed using the recombinant proteins, including monkey CYP2C20, CYP2C43, CYP2C75, and CYP2C76, in addition to human CYP2C8, CYP2C9, CYP2C18, and CYP2C19. Among these proteins, the antibodies detected a single ∼50-kDa band only in CYP2C76 (Fig. 6A), indicating the immunospecificity of the antibodies. The immunoblotting for CYP2C76 was also carried out with liver microsomes prepared from five primate species: human, chimpanzee, orangutan, and cynomolgus and rhesus monkeys. A signal of the expected size was clearly detected only with liver microsomes from cynomolgus and rhesus monkeys (Fig. 6B), coincided well with the results of the gene expression pattern of CYP2C76.
To determine the cellular localization of the CYP2C76 protein, the cryosections prepared from the cynomolgus monkey liver were stained with the anti-CYP2C76 antibodies. Strong staining was seen in the cytoplasm of hepatocytes but not in the cells lining the bile duct or the vein (Fig. 7A). Little or no staining was seen after peptide blocking (Fig. 7B) or with preimmune serum (Fig. 7C), indicating that the staining was specific for CYP2C76.
Discussion
Monkeys have been employed in the preclinical studies of drug metabolism because they are believed to show a metabolic profile similar to that of humans. However, a different pattern is occasionally seen in drug metabolism between monkeys and humans (Stevens et al., 1993; Sharer et al., 1995; Guengerich, 1997; Weaver et al., 1999; Bogaards et al., 2000; Narimatsu et al., 2000). The molecular mechanism(s) behind this phenomenon remains unclear, partly as a result of the lack of detailed information on the genes and molecules responsible for drug metabolism in monkeys. In this study, to understand a possible cause responsible for this species difference, we identified a cDNA for CYP2C76 and characterized along with other CYP2C cDNAs encoding CYP2C20, CYP2C43, and CYP2C75 in cynomolgus monkeys. We investigated species specificity and tissue distribution of gene expression, in addition to genomic organization and metabolic properties.
The CYP2C subfamily has been known to be diverged in each species during evolution (Nelson et al., 1996, 2004), indicating that species specificity and orthologous relationship must be determined cautiously for the CYP2C subfamily. In this report, gene and protein expression pattern, genomic organization, and metabolic properties suggest that CYP2C76 does not seem to have the ortholog in humans. During the preparation of this article, a partial cDNA sequence for rhesus CYP2C76 was reported under GenBank accession number CX078602 as one of the clones identified by EST sequencing (Magness et al., 2005). In their report, PCR using gene-specific primers with the genomic DNA from several primate species showed the amplification in macaques but not in humans, further supporting our results. A definite conclusion could be made when a complete sequence of the monkey genome will be available.
The highly homologous CYP2C genes tend to be located near each other within the gene cluster, whereas CYP2C44, least homologous to other subfamily members, is located ∼4 × 106 bases away from the mouse CYP2C gene cluster (Nelson et al., 2004). In contrast, CYP2C76 (also least identical to any other CYP2C genes and thus placed outside the CYP2C group in the phylogenetic tree) is located within but on the edge of the CYP2C cluster in the monkey genome, similar to the mouse CYP2C70 gene. The outer location probably did not allow CYP2C76 for efficient crossover, resulting in the lower identity to other CYP2C subfamily members as has been proposed for CYP2C44 and CYP2C70 (Nelson et al., 2004). It is of great interest to know how CYP2C76 has arisen after human and Old World monkeys diverged from a common ancestor around 25 million years ago (Kumar and Hedges 1998). Sequencing the CYP2C cluster of the closely related primate species should give an insight into this question.
Transcript variants influence the function of the P450 genes. Two transcript variants of CYP4F3 contain either exon 3 or exon 4 generated through alternative splicing, leading to the high varieties of the synthesized proteins to accommodate different substrates as well as tissue specificity in gene expression (Christmas et al., 2001). In contrast, we identified a transcript variant for CYP2C76 that lacks one of the exons, probably because of alternative splicing, which does not seem to have functional importance because of the premature termination codon (PTC) generated (data not shown). The PTC mRNAs can be subjected to a rapid degradation by nonsense-mediated decay (NMD) when PTCs are located more than 50 nucleotides before the last exon-exon junction (Holbrook et al., 2004). NMD is responsible for the degradation of the transcript variants generated by the CYP3A5*3 allele, which can explain the difference in the expression levels of mRNAs between this and other genotype groups (Kuehl et al., 2001; Busi and Cresteil, 2005). Moreover, at least one third of alternative transcripts in humans were identified as PTC mRNAs, potential targets for RNA decay pathway through NMD (Lewis et al., 2003). Because gene expression is not completely diminished by NMD, expression levels remained would vary considerably between RNA isoforms, cell types, and even individuals (Holbrook et al., 2004). Therefore, together with the effect of SNPs on alternative splicing, the transcript variant of CYP2C76 we identified, if subjected to NMD, could increase functional complexity of this gene.
CYP2C76 showed characteristic metabolic properties compared with CYP2C20, CYP2C43, and CYP2C75. CYP2C76 was critically involved in the metabolism of tolbutamide and testosterone, the metabolic properties of which are different from the other monkey CYP2Cs analyzed. CYP2C76, as a species-specific CYP2C enzyme, certainly adds the complexity to drug metabolism in monkeys and might account for the species difference in drug metabolism between monkeys and humans. Indeed, we have revealed that the CYP2C76 is at least partly responsible for the difference between monkeys and humans in the metabolism of a currently prescribed drug (Y. Uno, Y. Suzuki, Y. Sakamoto, H. Sano, K. Hashimoto, S. Sugano, and I. Inoue, unpublished observations). Further investigation of CYP2C76 on metabolic activity for a variety of substrates will help to better understand drug metabolism in monkeys and species difference between monkeys and humans.
The species difference of drug metabolism has been a major issue in drug development, because the results obtained with experimental animals need to be extrapolated to humans. The analysis for the capacity of hepatic microsomes from several animal species to metabolize drugs revealed that some differences occur in monkeys for the metabolism of marker substrates compared with humans (Sharer et al., 1995; Weaver et al., 1999; Bogaards et al., 2000). Tolbutamide hydroxylase activity is at least 3-fold lower in cynomolgus monkeys than humans (Sharer et al., 1995; Weaver et al., 1999), in contrary to our expectation because of monkey CYP2C75 and CYP2C76 exhibiting activities toward this substrate. This discrepancy might be due to the lower enzymatic activities of these enzymes compared with CYP2C9. However, other factors also need to be considered, such as genetic polymorphisms in individuals, from which the cDNA sequence used for protein expression was derived, because as in humans, monkeys have a diverse genetic background. Considering that even a single amino acid substitution can alter the enzymatic activity of P450s (Guengerich, 1997), nonsynonymous single nucleotide polymorphisms might have reduced the activity of CYP2C75 and CYP2C76. To examine this possibility, the genetic polymorphisms in CYP2C76 and other CYP2C genes must be identified and characterized.
Understanding the species difference and the mechanisms behind it is an inevitable task to improve the accuracy in extrapolating the animal data to humans. To accomplish this goal, advancing genomic techniques such as EST or genome sequencing, microarray, and comparative genomics should be helpful because they can identify genomic components specific for each animal species including genes, transcripts, and regulatory elements. The ESTs specific for cynomolgus and rhesus monkeys have been identified by our and other groups (Magness et al., 2005), including CYP2C76. Further identification and characterization of monkey-specific ESTs will help to better understand the species uniqueness of monkeys in drug metabolism.
In conclusion, we have identified cynomolgus CYP2C76, which does not have the corresponding gene in the human genome. CYP2C76 contains nine exons and is located in a single CYP2C cluster in the monkey genome, similar to the human CYP2C genes. Our data show that CYP2C76 is predominantly expressed in the liver, and its expression level is the greatest among the four CYP2C genes analyzed. Moreover, CYP2C76 has a characteristic metabolic profile different from the other CYP2Cs. From these observations, we conclude that cynomolgus CYP2C76 is a major CYP2C contributing substantially to overall drug-metabolizing activity in the liver.
Acknowledgments
We thank the GAIN for providing us invaluable tissue samples from orangutan and chimpanzee, and Ms. Makiko Hase (Applied Biosystems) for assistance in designing gene-specific primers and probes for real-time RT-PCR. We also would like to acknowledge Dr. Aleksandar Milosavljevic at the Baylor College of Medicine for allowing us to access invaluable information on the rhesus monkey genome.
Footnotes
-
ABBREVIATIONS: P450, cytochrome P450; PCR, polymerase chain reaction; BAC, bacterial artificial chromosome; RT, reverse transcription; TLC, thin-layer chromatography; PDI, protein disulfide isomerase; EST, expressed sequence tag; PTC, premature termination codon; NMD, nonsense-mediated decay.
- Received January 20, 2006.
- Accepted April 21, 2006.
- The American Society for Pharmacology and Experimental Therapeutics