CYP2C19 polymorphisms account for inter-individual variability of drug metabolism in cynomolgus macaques

Key regulators in the architecture of substrate access/egress channels in mammalian cytochromes P450 governing flexibility in substrate oxyfunctionalization

Cytochrome P450s (CYP) represent a superfamily of b-type hemoproteins catalyzing oxifunctionalization of a vast array of endogenous and exogenous compounds. The present review focuses on assessment of the topology of prospective determinants in substrate entry and product release channels of mammalian P450s, steering the conformational dynamics of substrate accessibility and productive ligand orientation toward the iron-oxene core. Based on a generalized, CYP3A4-related construct, the sum of critical elements from diverse target enzymes was found to cluster within the known substrate recognition sites. The majority of prevalent substrate access/egress tunnels revealed to be of fairly balanced functional importance. The hydrophobicity profile of the candidates revealed to be the most salient feature in functional interaction throughout the conduits, while bulkiness of the residues imposes steric restrictions on substrate traveling. Thus, small amino acids such as prolines and glycines serve as hinges, driving conformational flexibility in ligand passage. Similarly, bottlenecks in the tunnel architecture, being narrowest encounter points within the CYP3A4 model, have a vital function in substrate selectivity along with clusters of aromatic amino acids acting as gatekeepers. In addition, peripheral patches in conduits may house determinants modulating allosteric cooperativity between remote and central domains in the P450 structure. Remarkably, the bulk critical residues lining tunnels in the various isozymes reside in helices B'/C and F/G inclusive of their interhelical turns as well as in helix I. This suggests these regions to represent hotspots for targeted genetic engineering to tailor more sophisticated mammalian P450s exploitable in industrial, biotechnological and medicinal areas.

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.

Recent advances in computational metabolite structure predictions and altered metabolic pathways assessment to inform drug development processes

Many drug candidates fail during preclinical and clinical trials due to variable or unexpected metabolism which may lead to variability in drug efficacy or adverse drug reactions. The drug metabolism field aims to address this important issue from many angles which range from the study of drug-drug interactions, pharmacogenomics, computational metabolic modeling, and others. This manuscript aims to provide brief but comprehensive manuscript summaries highlighting the conclusions and scientific importance of seven exceptional manuscripts published in recent years within the field of drug metabolism. Two main topics within the field are reviewed: novel computational metabolic modeling approaches which provide complex outputs beyond site of metabolism predictions, and experimental approaches designed to discern the impacts of interindividual variability and species differences on drug metabolism. The computational approaches discussed provide novel outputs in metabolite structure and formation likelihood and/or extend beyond the saturated field of drug phase I metabolism, while the experimental metabolic pathways assessments aim to highlight the impacts of genetic polymorphisms and clinical animal model metabolic differences on human metabolism and subsequent health outcomes.

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.

RNA sequencing of kidney and liver transcriptome obtained from wild cynomolgus macaque (Macaca fascicularis) originating from Peninsular Malaysia

Objective

Using high-throughput RNA sequencing technology, this study aimed to sequence the transcriptome of kidney and liver tissues harvested from Peninsular Malaysia cynomolgus macaque (Macaca fascicularis). M. fascicularis are significant nonhuman primate models in the biomedical field, owing to the macaque’s biological similarities with humans. The additional transcriptomic dataset will supplement the previously described Peninsular Malaysia M. fascicularis transcriptomes obtained in a past endeavour.

Results

A total of 75,350,240 sequence reads were obtained via Hi-seq 2500 sequencing technology. A total of 5473 significant differentially expressed genes were called. Gene ontology functional categorisation showed that cellular process, catalytic activity, and cell part categories had the highest number of expressed genes, while the metabolic pathways category possessed the highest number of expressed genes in the KEGG pathway analysis. The additional sequence dataset will further enrich existing M. fascicularis transcriptome assemblies, and provide a dataset for further downstream studies.

Electronic supplementary material

The online version of this article (10.1186/s13104-018-4014-1) contains supplementary material, which is available to authorized users.

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.

Molecular and Functional Characterization of N-Acetyltransferases NAT1 and NAT2 in Cynomolgus Macaque

Arylamine N-acetyltransferases (NATs) are drug-metabolizing enzymes essential for the metabolism of endogenous substrates and xenobiotics, and their molecular characteristics have been extensively investigated in humans, but not in cynomolgus macaques, nonhuman primate species important for drug metabolism studies. In this study, cynomolgus NAT1 and NAT2 cDNAs were isolated from livers. NAT1 and NAT2 were characterized by molecular analyses and drug-metabolizing assays. A total of 9 transcript variants were found for cynomolgus NAT1, similar to human NAT1, and contained 1-4 exons with the coding region largely conserved with human NAT1. Genomic organization was similar between cynomolgus macaques and humans. Cynomolgus NAT1 and NAT2 amino acid sequences showed high sequence homology (95% and 89%, respectively) and showed close relationships with human NAT1 and NAT2 in a phylogenetic tree. Cynomolgus NAT2 mRNA was predominantly expressed in liver among the 10 different tissues analyzed, followed by kidney and jejunum. In contrast, cynomolgus NAT1 mRNA showed more ubiquitous expression with relatively more abundant expression in liver, kidney, and jejunum, along with testis. Metabolic assays using recombinant proteins showed that cynomolgus NAT1 and NAT2 metabolized human NAT substrates, including p-aminobenzoic acid, sulfamethazine, isoniazid, and 2-aminofluorene. Interestingly, p-aminobenzoic acid and isoniazid were largely metabolized by NAT1 and NAT2, respectively, in cynomolgus macaques and humans; sulfamethazine, a human NAT2 substrate, was metabolized by both NAT enzymes in cynomolgus macaques. These results suggest molecular and enzymatic similarities of NAT1 and NAT2 between cynomolgus macaques and humans, despite some small differences in substrate specificity of the enzymes.

Genetic Variants of Glutathione S-Transferase GSTT1 and GSTT2 in Cynomolgus Macaques: Identification of GSTT Substrates and Functionally Relevant Alleles

Glutathione S-transferase (GST) is a family of important drug-metabolizing enzymes, conjugating endogenous and exogenous compounds. Genetic polymorphisms result in the inter-individual variability of GST activity in humans. Especially, human GSTT1 and GSTT2 null alleles are associated with toxicity and various cancers derived from chemicals. Cynomolgus macaque, a nonhuman primate species widely used in drug metabolism studies, has molecular and enzymatic similarities of GSTs to the human orthologs; however, genetic polymorphisms have not been investigated in this species. In this study, resequencing of GSTT1 and GSTT2 in 64 cynomolgus and 32 rhesus macaques found 15 nonsynonymous variants and 1 nonsense variant for GSTT1 and 15 nonsynonymous variants for GSTT2. Some of these GSTT variants were distributed differently in Indochinese and Indonesian cynomolgus macaques and rhesus macaques. For analysis of functional relevance of the GSTT variants, 1-iodohexane and dibromomethane were determined to be suitable substrates for cynomolgus GSTT1 and GSTT2. However, the conjugation activities were roughly correlated with GSTT protein levels immunochemically quantified in cynomolgus liver samples with no statistical significances, implying the contributions of the GST genetic variants. Among the GSTT1 variants identified, the animals carrying R76C and D125G mutations showed lower conjugation activities toward dibromomethane than those of the wild-type in liver cytosolic fractions. Moreover, the recombinant R76C/D125G and D125G GSTT variant proteins showed significantly lower 1-iodohexane or dibromomethane conjugation activities than those of the wild-type protein. Therefore, inter-animal variability of GSTT-dependent drug metabolism is at least partly accounted for by GSTT1 and possibly GSTT2 variants in cynomolgus and rhesus 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.

Single nucleotide polymorphisms in the FcγR3A and TAP1 genes impact ADCC in cynomolgus monkey PBMCs

Phenotypic variability is often observed in cynomolgus monkeys on preclinical studies and may, in part, be driven by genetic variability. However, the role of monkey genetic variation remains largely unexplored in the context of drug response. This study evaluated genetic variation in cynomolgus monkey FcγR3A and TAP1 genes and the potential impact of identified polymorphisms on antibody-dependent cell-mediated cytotoxicity (ADCC) in vitro. Studies in humans have demonstrated that a single nucleotide polymorphism (SNP), F158V, in FcγR3A can influence response to rituximab through altered ADCC and that SNPs in TAP1/2 decrease natural killer (NK) cell activity against major histocompatibility complex (MHC) class I deficient cells, potentially through altered ADCC. Monkeys were genotyped for FcγR3A and TAP1 SNPs, and ADCC was assessed in vitro using peripheral blood mononuclear cells (PBMCs) treated with trastuzumab in the presence of NCI-N87 cells. FcγR3A g.1134A>C (exonic S42R), FcγR3A g.5027A>G (intronic), and TAP1 g.1A>G (start codon loss) SNPs were all significantly associated with decreased ADCC for at least one trastuzumab concentration ≥0.0001 μM when compared with wild type (WT). Regression analysis demonstrated significant association of the SNP-SNP pairs FcγR3A g.1134A>C/TAP1 g.1A>G and FcγR3A g.5027A>G/TAP1 g.1A>G with a combinatorial decrease on ADCC. Mechanisms underlying the decreased ADCC were investigated by measuring FcγR3A/IgG binding affinity and expression of FcγR3A and TAP1 in PBMCs; however, no functional associations were observed. These data demonstrate that genetic variation in cynomolgus monkeys is reflective of known human genetic variation and may potentially contribute to variable drug response in preclinical 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.

The Effects of CYP2C19 genotype on the susceptibility for nephrosis in cardio-cerebral vascular disease treated by anticoagulation

In recent years, the genetic factor has become one of the important predisposing factors of nephropathy susceptibility. There is a high incidence of nephropathy in CCVd. The CYP2C19 enzyme metabolizes most the drugs, including proton pump inhibitors commonly used medicines to treat CCVd, CYP2C19 genetic polymorphisms is association with multi-pathogenesis factors of nephropathy. The purpose of the study is to reveal the association between CYP2C19 genotype and the susceptibility of nephropathy in the CCVd patients. The study is composed of 623 samples from CCVd treated by anticoagulation. The patients were studied, including CCVd with hyperuricemia, coronary heart disease, diabetes, and other complication. Biochemical tests and CYP2C19 variants measurements were performed by the gene chip method. The association among CYP2C19 variants, complications, and nephropathy was analyzed in the CCVd. There is no correlation between nephropathy and complications in CCVd. In hyperuricemia, coronary heart disease and diabetes groups, the differences of renal function tests were significant between CYP2C19 mutant (P < 0.05). The nephropathy risk of wild genotype is 3.288 times higher than of mutation genotype in hyperuricemic group, 1.928 times higher than mutation genotype in coronary heart disease group, and 5.248 times higher than CYP2C19 mutation genotype in the diabetic group. There was significant correlation between the CYP2C19 wild type and the nephropathy susceptibility in CCVd patients. The CYP2C19 gene plays a potential maker to evaluate nephropathy in CCVd patients. We deduced that identification of CYP2C19 gene type may benefit for reducing and avoiding nephropathy caused by abnormal metabolism function in CCVd patients.

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.

Enhanced methamphetamine metabolism in rhesus macaque as compared with human: an analysis using a novel method of liquid chromatography with tandem mass spectrometry, kinetic study, and substrate docking

Methamphetamine (MA), which remains one of the widely used drugs of abuse, is metabolized by the cytochrome P450 (P450) family of enzymes in humans. However, metabolism of methamphetamine in macaques is poorly understood. Therefore, we first developed and validated a very sensitive liquid chromatography with tandem mass spectrometry (LC-MS/MS) method using solid phase extraction of rhesus plasma with a lower limit of quantitation at 1.09 ng/ml for MA and its metabolites, 4-hydroxy methamphetamine (4-OH MA), amphetamine (AM), 4-OH amphetamine (4-OH AM), and norephedrine. We then analyzed plasma samples of MA-treated rhesus, which showed >10-fold higher concentrations of AM (∼29 ng/ml) and 4-OH AM (∼28 ng/ml) than MA (∼2 ng/ml). Because the plasma levels of MA metabolites in rhesus were much higher than in human samples, we examined MA metabolism in human and rhesus microsomes. Interestingly, the results showed that AM and 4-OH AM were formed more rapidly and that the catalytic efficiency (Vmax/Km) for the formation of AM was ∼8-fold higher in rhesus than in human microsomes. We further examined the differences in these kinetic characteristics using three selective inhibitors of each human CYP2D6 and CYP3A4 enzymes. The results showed that each of these inhibitors inhibited both d- and l-MA metabolism by 20%-60% in human microsomes but not in rhesus microsomes. The differences between human and rhesus CYP2D6 and CYP3A4 enzymes were further assessed by docking studies for both d and l-MA. In conclusion, our results demonstrated an enhanced MA metabolism in rhesus compared with humans, which is likely to be caused by differences in MA-metabolizing P450 enzymes between these species.