Slow R-warfarin 7-hydroxylation mediated by P450 2C19 genetic variants in cynomolgus monkeys in vivo

Polymorphic cytochromes P450 in non-human primates

Cynomolgus macaques (Macaca fascicularis, an Old World monkey) are widely used in drug development because of their genetic and physiological similarities to humans, and this trend has continued with the use of common marmosets (Callithrix jacchus, a New World monkey). Information on the major drug-metabolizing cytochrome P450 (CYP, P450) enzymes of these primate species indicates that multiple forms of their P450 enzymes have generally similar substrate selectivities to those of human P450 enzymes; however, some differences in isoform, activity, and substrate specificity account for limited species differences in drug oxidative metabolism. This review provides information on the P450 enzymes of cynomolgus macaques and marmosets, including cDNA, tissue expression, substrate specificity, and genetic variants, along with age differences and induction. Typical examples of important P450s to be considered in drug metabolism studies include cynomolgus CYP2C19, which is expressed abundantly in liver and metabolizes numerous drugs. Moreover, genetic variants of cynomolgus CYP2C19 affect the individual pharmacokinetic data of drugs such as R-warfarin. These findings provide a foundation for understanding each P450 enzyme and the individual pharmacokinetic and toxicological results in cynomolgus macaques and marmosets as preclinical models. In addition, the effects of induction on some drug clearances mediated by P450 enzymes are also described. In summary, this review describes genetic and acquired individual differences in cynomolgus and marmoset P450 enzymes involved in drug oxidation that may be associated with pharmacological and/or toxicological effects.

Genetic variation in the Mauritian cynomolgus macaque population reflects variation in the human population

The cynomolgus macaque is an important species for preclinical research, however the extent of genetic variation in this population and its similarity to the human population is not well understood. Exome sequencing was conducted for 101 cynomolgus macaques to characterize genetic variation. The variant distribution frequency was 7.81 variants per kilobase across the sequenced regions, with a total of 2,770,009 single nucleotide variants identified from 2,996,041 loci. A large portion (85.6%) had minor allele frequencies greater than 5%. Enriched pathways for genes with high genetic diversity (≥10 variants per kilobase) were those involving signaling peptides and immune response. Compared to human, the variant distribution frequency and nucleotide diversity in the macaque exome was approximately 4 times greater; however the ratio of non-synonymous to synonymous variants was similar (0.735 and 0.831, respectively). Understanding genetic variability in cynomolgus macaques will enable better interpretation and human translation of phenotypic variability in this species.

Human total clearance values and volumes of distribution of typical human cytochrome P450 2C9/19 substrates predicted by single-species allometric scaling using pharmacokinetic data sets from common marmosets genotyped for P450 2C19

Common marmosets (Callithrix jacchus) are small non-human primates that genetically lack cytochrome P450 2C9 (CYP2C9). Polymorphic marmoset CYP2C19 compensates by mediating oxidations of typical human CYP2C9/19 substrates.Twenty-four probe substrates were intravenously administered in combinations to marmosets assigned to extensive or poor metaboliser (PM) groups by CYP2C19 genotyping. Eliminations from plasma of cilomilast, phenytoin, repaglinide, tolbutamide, and S-warfarin in the CYP2C19 PM group were significantly slow; these drugs are known substrates of human CYP2C8/9/19.Human total clearance values and volumes of distribution of the 24 test compounds were extrapolated using single-species allometric scaling with experimental data from marmosets and found to be mostly comparable with the reported values.Human total clearance values and volumes of distribution of 15 of the 24 test compounds similarly extrapolated using reported data sets from cynomolgus or rhesus monkeys were comparable to the present predicted results, especially to those based on data from PM marmosets.These results suggest that single-species allometric scaling using marmosets, being small, has advantages over multiple-species-based allometry and could be applicable for pharmacokinetic predictions at the discovery stage of drug development.

Cytochrome P450 1A1, 2C9, 2C19, and 3A4 Polymorphisms Account for Interindividual Variability of Toxicological Drug Metabolism in Cynomolgus Macaques

Cytochromes P450 (P450s) and their genetic variants in humans are important drug-metabolizing enzymes partly accounting for interindividual variations in drug metabolism and toxicity. However, these genetic variants in P450s have not been fully investigated in cynomolgus macaques, a nonhuman primate species widely used in toxicological studies. In this study, genetic variants found in cynomolgus CYP1A1, CYP2C9 (formerly CYP2C43), CYP2C19 (CYP2C75), and CYP3A4 (CYP3A8) were assessed on functional importance. Resequencing of CYP1A1 in cynomolgus macaques found 18 nonsynonymous variants, of which M121I and V382I were located in SRSs, domains potentially important for P450 function. By further analyzing these two variants, V382I was significantly associated with lower drug-metabolizing activities in the liver for the heterozygotes than the wild types. Similarly, the heterozygotes or homozygotes of CYP2C9 variants (A82V and H344R) and CYP2C19 variant (A490V) showed significantly lower drug-metabolizing activities in the liver than the wild types. Moreover, the homozygotes of CYP3A4 variant (S437N) showed significantly higher activities than the wild type in the liver. Kinetic analyses using recombinant proteins revealed that CYP2C9 variants (A82V and H344R) showed substantially lower  K_s values than the wild type, although CYP1A1 variant (V382I) showed kinetic parameters similar to the wild type. Likewise, CYP2C19 variant (A490V) showed substantially a lower  V_max/ K_m value than the wild type, whereas CYP3A4 variant (S437N) showed a higher V_max/ K_m value than the wild type. These results suggest the toxicologically functional importance of CYP2C9 variants (A82V and H344R), CYP2C19 variant (A490V), and CYP3A4 variant (S437N) for hepatic drug metabolism in cynomolgus macaques.

Development of genotyping method for functionally relevant variants of cytochromes P450 in cynomolgus macaques

In cynomolgus macaques (Macaca fascicularis), widely used in drug metabolism studies, CYP2C9, CYP2C76, CYP2D6, CYP3A4, and CYP3A5, important drug-metabolizing enzymes, are abundantly expressed in liver and metabolize cytochrome P450 substrates. CYP2C9 (c.334A>C), CYP2C76 (c.449TG>A), CYP2D6 (c.891A>G), CYP3A4 (IVS3 + 1G>del), and CYP3A5 (c.625A>T) substantially influence metabolic activity of enzymes, and thus are important variants in drug metabolism studies. In this study, a real-time PCR method was developed for genotyping these variants. The validity of the methods was verified by genotyping two wild type, two heterozygous, and two homozygous DNAs and was used to genotype 41 cynomolgus macaques (from Cambodia, Indonesia, the Philippines, or Vietnam) for the five variants, along with another important variant CYP2C19 (c.308C>T). The CYP2C9 and CYP2C19 variants were found only in Cambodian and Vietnamese animals, while the CYP2C76 and CYP2D6 variants were found only in Indonesian and Philippine animals. The CYP3A4 and CYP3A5 variants were not found in any of the animals analyzed. Mauritian animals, genotyped using next-generation sequencing data for comparison, possessed the CYP2C19 and CYP2D6 variants, but not the other variants. These results indicated differences in prevalence of these important variants among animal groups. Therefore, the genotyping tool developed is useful for drug metabolism studies using cynomolgus macaques.

Effects of aging and rifampicin pretreatment on the pharmacokinetics of human cytochrome P450 probes caffeine, warfarin, omeprazole, metoprolol and midazolam in common marmosets genotyped for cytochrome P450 2C19

1. The pharmacokinetics were investigated for human cytochrome P450 probes after single intravenous and oral administrations of 0.20 and 1.0 mg/kg, respectively, of caffeine, warfarin, omeprazole, metoprolol and midazolam to aged (10-14 years old, n = 4) or rifampicin-treated/young (3 years old, n = 3) male common marmosets all genotyped as heterozygous for a cytochrome P450 2C19 variant. 2. Slopes of the plasma concentration-time curves after intravenous administration of warfarin and midazolam were slightly, but significantly (two-way analysis of variance), decreased in aged marmosets compared with young marmosets. The mean hepatic clearances determined by in silico fitting for individual pharmacokinetic models of warfarin and midazolam in the aged group were, respectively, 23% and 56% smaller than those for the young group. 3. Significantly enhanced plasma clearances of caffeine, warfarin, omeprazole and midazolam were evident in young marmosets pretreated with rifampicin (25 mg/kg daily for 4 days). Two- to three-fold increases in hepatic intrinsic clearance values were observed in the individual pharmacokinetic models. 4. The in vivo dispositions of multiple simultaneously administered drugs in old, young and P450-enzyme-induced marmosets were elucidated. The results suggest that common marmosets could be experimental models for aged, induced or polymorphic P450 enzymes in P450-dependent drug metabolism studies.

In vivo individual variations in pharmacokinetics of efavirenz in cynomolgus monkeys genotyped for cytochrome P450 2C9

Cynomolgus monkeys are used frequently in preclinical studies for new drug development due to their evolutionary closeness to humans. An antiretroviral drug, efavirenz, is a typical probe substrate for human cytochrome P450 (P450) 2B6, but is mainly metabolized by cynomolgus monkey P450 2C9. In this study, plasma concentrations of efavirenz were assessed in six cynomolgus monkeys genotyped for P450 2C9 c.334 A > C (I112L) (three wild-type, one heterozygote and two homozygotes) by high performance liquid chromatography with tandem mass spectrometry. After intravenous administration at a dose of 1.0 mg/kg, biphasic plasma elimination curves of efavirenz were seen in these cynomolgus monkeys. The mean plasma concentration of the primary metabolite 8-hydroxyefavirenz (1 h after treatment, with hydrolysis by β-glucuronidase) in the wild-type group was significantly higher (4.0-fold) than the combined heterozygous and homozygous group mean. The area under the plasma concentration-time curve value of efavirenz in the homozygous group after oral administration at a dose of 2.0 mg/kg was significantly higher (2.0-fold) than the combined wild-type and heterozygous group. These results collectively indicated that P450 2C9 c.334 A > C (I112L) variation was associated with efavirenz metabolic clearance in vivo. Cynomolgus P450 2C9 polymorphism might account for interindividual variations of efavirenz metabolism in cynomolgus monkeys. Copyright © 2016 John Wiley & Sons, Ltd.

Differences in Toxicological and Pharmacological Responses Mediated by Polymorphic Cytochromes P450 and Related Drug-Metabolizing Enzymes

Research over the past 30 years has elucidated the roles of polymorphic human liver cytochrome P450 (P450) enzymes associated with toxicological and/or pharmacological actions. Thalidomide exerts its various pharmacological and toxic actions in primates through multiple mechanisms, including nonspecific modification of many protein networks after bioactivation by autoinduced human P450 enzymes. To overcome species differences between rodents, currently, nonhuman primates and/or mouse models with transplanted human hepatocytes are used. Interindividual variability of P450-dependent drug clearances in cynomolgus monkeys and common marmosets is partly accounted for by polymorphic P450 variants and/or aging, just as it is in humans with increased prevalence of polypharmacy. Genotyping of P450 genes in nonhuman primates would be beneficial before and/or after drug metabolism and toxicity testing and evaluation as well in humans. Genome-wide association studies in humans have been rapidly advanced; however, unique whole-gene deletion of P450 2A6 was subsequently developed to cover nicotine-related lung cancer risk study. Regarding polypharmacy, toxicological research should generally be aimed at identifying the risk of adverse drug events following specific potential drug exposures by examining single or multiple metabolic pathways involving single or multiple drug-metabolizing enzymes. Current and next-generation research of drug metabolism and disposition resulting in drug toxicity would be addressed under advanced knowledge of polymorphic and age-related intra- and/or interspecies differences of drug-metabolizing enzymes. In the near future, humanized animal models combining transplanted hepatocytes and a humanized immune system may be available to study human immune reactions caused by human-type drug metabolites. Such sophisticated models should provide preclinical predictions of human drug metabolism and potential toxicity.

Analysis of the chemical toxicity effects using the enrichment of Gene Ontology terms and KEGG pathways

BACKGROUND

Chemical toxicity is one of the major barriers for designing and detecting new chemical entities during drug discovery. Unexpected toxicity of an approved drug may lead to withdrawal from the market and significant loss of the associated costs. Better understanding of the mechanisms underlying various toxicity effects can help eliminate unqualified candidate drugs in early stages, allowing researchers to focus their attention on other more viable candidates.

METHODS

In this study, we aimed to understand the mechanisms underlying several toxicity effects using Gene Ontology (GO) terms and KEGG pathways. GO term and KEGG pathway enrichment theories were adopted to encode each chemical, and the minimum redundancy maximum relevance (mRMR) was used to analyze the GO terms and the KEGG pathways. Based on the feature list obtained by the mRMR method, the most related GO terms and KEGG pathways were extracted.

RESULTS

Some important GO terms and KEGG pathways were uncovered, which were concluded to be significant for determining chemical toxicity effects.

CONCLUSIONS

Several GO terms and KEGG pathways are highly related to all investigated toxicity effects, while some are specific to a certain toxicity effect.

GENERAL SIGNIFICANCE

The findings in this study have the potential to further our understanding of different chemical toxicity mechanisms and to assist scientists in developing new chemical toxicity prediction algorithms. This article is part of a Special Issue entitled "System Genetics" Guest Editor: Dr. Yudong Cai and Dr. Tao Huang.

Development of a genotyping tool for a functionally relevant CYP2C19 allele (Phe100Asn, Ala103Val and Ile112Leu) in cynomolgus macaques

In cynomolgus macaques, which are widely used in drug metabolism studies, CYP2C19 (formerly known as CYP2C75) is abundantly expressed in liver, metabolizes human CYP2C substrates and is thus an important drug-metabolizing enzyme. One of the cynomolgus CYP2C19 alleles (p.Phe100Asn, p.Ala103Val and p.Ile112Leu) results in substantially reduced metabolic activity and thus is an important allele in drug metabolism studies. For this allele, a genotyping tool was developed using allele-specific TaqMan probe. Genotyping 40 Cambodian cynomolgus macaques using this tool found 1 homozygote, 17 heterozygotes and 22 wild type animals, and the result was confirmed by direct-sequencing. Interestingly, this allele frequency was similar to that of Chinese cynomolgus macaques. The genotyping tool established is useful for drug metabolism studies using cynomolgus macaques.

Novel Marmoset Cytochrome P450 2C19 in Livers Efficiently Metabolizes Human P450 2C9 and 2C19 Substrates, S-Warfarin, Tolbutamide, Flurbiprofen, and Omeprazole

The common marmoset (Callithrix jacchus), a small New World monkey, has the potential for use in human drug development due to its evolutionary closeness to humans. Four novel cDNAs, encoding cytochrome P450 (P450) 2C18, 2C19, 2C58, and 2C76, were cloned from marmoset livers to characterize P450 2C molecular properties, including previously reported P450 2C8. The deduced amino acid sequence showed high sequence identities (>86%) with those of human P450 2Cs, except for marmoset P450 2C76, which has a low sequence identity (∼70%) with any human P450 2Cs. Phylogenetic analysis showed that marmoset P450 2Cs were more closely clustered with those of humans and macaques than other species investigated. Quantitative polymerase chain reaction analysis showed that all of the marmoset P450 2C mRNAs were predominantly expressed in liver as opposed to the other tissues tested. Marmoset P450 2C proteins were detected in liver by immunoblotting using antibodies against human P450 2Cs. Among marmoset P450 2Cs heterologously expressed in Escherichia coli, marmoset P450 2C19 efficiently catalyzed human P450 2C substrates, S-warfarin, diclofenac, tolbutamide, flurbiprofen, and omeprazole. Marmoset P450 2C19 had high Vmax and low Km values for S-warfarin 7-hydroxylation that were comparable to those in human liver microsomes, indicating warfarin stereoselectivity similar to findings in humans. Faster in vivo S-warfarin clearance than R-warfarin after intravenous administration of racemic warfarin (0.2 mg/kg) to marmosets was consistent with the in vitro kinetic parameters. These results indicated that marmoset P450 2C enzymes had functional characteristics similar to those of humans, and that P450 2C-dependent metabolic properties are likewise similar between marmosets and humans.