Search

Pigs are sometimes used in preclinical drug metabolism studies, with growing interest, and thus their drug-metabolizing enzymes, including the cytochromes P450 (P450 or CYP; EC 1.14.14.1), need to be examined. In the present study, novel CYP4A cDNAs were isolated and characterized, namely, pig CYP4A23 and CYP4A90; cat CYP4A37 and CYP4A106; and tree shrew CYP4A11a, CYP4A11d, CYP4A11e, CYP4A11f, and CYP4A11g. For comparison, the following known CYP4A cDNAs were also analyzed: pig CYP4A21 and dog CYP4A37, CYP4A38, and CYP4A39. These CYP4A cDNAs all contained open reading frames of 504-513 amino acids and had high amino acid sequence identity (74%-80%) with human CYP4As. Phylogenetic analysis of amino acid sequences revealed that these CYP4As were clustered in each species. All CYP4A genes contained 12 coding exons and formed a gene cluster in the corresponding genomic regions. A range of tissue types were analyzed, and these CYP4A mRNAs were preferentially expressed in liver and/or kidney, except for pig CYP4A90, which showed preferential expression in lung and duodenum. CYP4A enzymes, heterologously expressed in Escherichia coli , preferentially catalyzed lauric acid 12-hydroxylation and arachidonic acid 20-hydroxylation, just as human CYP4A11 does, with the same regioselectivity (i.e., at the ω -position of fatty acids). These results imply that dog, cat, pig, and tree shrew CYP4As have functional characteristics similar to those of human CYP4A11, with minor differences in lauric acid 12-hydroxylation. SIGNIFICANCE STATEMENT: Cytochrome P450 (P450, CYP) 4As are important P450s in human biological processes because of their fatty acid-metabolizing ability. Pig CYP4A21, CYP4A23, and CYP4A90; cat CYP4A37 and CYP4A106; tree shrew CYP4A11a, CYP4A11d, CYP4A11e, CYP4A11f, and CYP4A11g; and dog CYP4A37, CYP4A38, and CYP4A39 cDNAs were isolated and analyzed. These CYP4A cDNAs shared relatively high sequence identities with human CYP4A11 and CYP4A22. Pig, cat, tree shrew, and dog CYP4As in the liver and kidneys are likely to catalyze the ω -hydroxylation of fatty acids., (Copyright © 2024 by The American Society for Pharmacology and Experimental Therapeutics.)

Omeprazole, a gastric acid pump inhibitor, is repeatedly administered and is oxidatively metabolized mainly by polymorphic cytochrome P450 2C19. The prescribed dosage of omeprazole was discontinued or reduced in 47 of the 135 patients who received omeprazole alone in this survey, as recorded in the Japanese Adverse Drug Event Report database. The days to onset of omeprazole-related disorders were 3-4 d (median) and 16 d for intravenous 20-40 mg and oral 20 mg daily doses, respectively, in 34 patients for whom relevant data were available. The maximum plasma concentration of omeprazole was pharmacokinetically modeled after a single oral 40-mg dose in P450 2C19-defective poor metabolizers and was 2.4-fold higher than that in extensive metabolizers. The modeled area under the hepatic concentration curves of omeprazole in P450 2C19 poor metabolizers after virtual daily 40-mg doses for 7 d was 5.2-fold higher than that in the extensive metabolizers. Omeprazole-induced P450 2C19 (approx. 2-fold), resulting in increased hepatic intrinsic clearance in repeated doses, was considered after the second day. Virtual plasma/hepatic exposure estimated using pharmacokinetic modeling in subjects with P450 2C19 poor metabolizers indicated that these exposure levels virtually estimated could be one of causal factors for unexpected hepatic disorders induced by prescribed omeprazole, such as those resulting from drug interactions with repeatedly co-administered medicines.

The hepatic elimination of chemical substances in pharmacokinetic models requires hepatic intrinsic clearance (CL h,int ) parameters for unbound drug in the liver, and these are regulated by the liver-to-plasma partition coefficients (K p,h ). Both Poulin and Theil and Rodgers and Rowland have proposed in silico expressions for K p,h for a variety of chemicals. In this study, two sets of in silico K p,h values for 14 model substances were assessed using experimentally reported in vivo steady-state K p,h data and time-dependent virtual internal exposures in the liver and plasma modeled by forward dosimetry in rats. The K p,h values for 14 chemicals independently calculated using the primary Poulin and Theil method in this study were significantly correlated with those obtained using the updated Rodgers and Rowland method and with reported in vivo steady-state K p,h data in rats. When pharmacokinetic parameters were derived based on individual in vivo time-dependent data for diazepam, phenytoin, and nicotine in rats, the modeled liver and plasma concentrations after intravenous administration of the selected substrates in rats using two sets of in silico K p,h values were mostly similar to the reported time-dependent in vivo internal exposures. Similar results for modeled liver and plasma concentrations were observed with input parameters estimated by machine-learning systems for hexobarbital, fingolimod, and pentazocine, with no reference to experimental pharmacokinetic data. These results suggest that the output values from rat pharmacokinetic models based on in silico K p,h values derived from the primary Poulin and Theil model would be applicable for estimating toxicokinetics or internal exposure to substances.

The estimation of health risks of chemical substances was historically investigated using animal studies; however, current research focuses on reducing the number of animal experiments. The toxicity of chemicals in fish screening systems is reportedly correlated with their hydrophobicity. The inverse relationship between absorption rates (intestinal cell permeability) and virtual hepatic/plasma pharmacokinetics of diverse chemicals has been previously evaluated by modeling oral administration in rats. In the current study, internal exposures, i.e., virtual maximum plasma concentrations (C max ) and areas under the concentration-time curves (AUC), of 56 food chemicals with reported hepatic lowest-observed-effect levels (LOELs) of ≤1000 mg/kg/d in rats were pharmacokinetically modeled using in silico estimated input pharmacokinetic parameters. After a virtual single oral dose of 1.0 mg/kg of 56 food chemicals, the output C max and AUC values in rat plasma generated by modeling using the corresponding in silico estimated input parameters were not significantly correlated with the reported hepatic LOEL values. However, significant inverse relationships between hepatic/plasma concentrations of selected lipophilic food chemicals (i.e., octanol-water partition coefficient log P > 1) using forward dosimetry and reported LOEL values (≤300 mg/kg/d) were observed (n = 14, r=-0.52-0.66, p ≤ 0.05). This simple modeling approach, which uses no experimental pharmacokinetic data, has the potential to play a significant role in reducing the use of animals to estimate toxicokinetics or internal exposures of lipophilic food components after oral doses. Therefore, these methods are valuable for estimating hepatic toxicity by using forward dosimetry in animal toxicity experiments.

Primary human hepatocytes (PHHs) have been commonly used as the gold standard in many drug metabolism studies, regardless of having large inter-individual variation. These inter-individual variations in PHHs arise primarily from genetic polymorphisms, as well as from donor health conditions and storage conditions prior to cell processing. To equalize the effects of the latter two factors, PHHs were transplanted to quality-controlled mice providing human hepatocyte proliferation niches, and engrafted livers were generated. Cells that were harvested from engrafted livers, call this as experimental human hepatocytes (EHH; termed HepaSH cells), were stably and reproducibly produced from 1014 chimeric mice produced by using 17 different PHHs. Expression levels of acute phase reactant (APR) genes as indicators of a systemic reaction to the environmental/inflammatory insults of liver donors varied widely among PHHs. In contrast to PHHs, the expression of APR genes in HepaSH cells was found to converge within a narrower range than in donor PHHs. Further, large individual differences in the expression levels of drug metabolism-related genes (28 genes) observed in PHHs were greatly reduced among HepaSH cells produced in a unified in vivo environment, and none deviated from the range of gene expression levels in the PHHs. The HepaSH cells displayed a similar level of drug-metabolizing enzyme activity and gene expression as the average PHHs but retained their characteristics for drug-metabolizing enzyme gene polymorphisms. Furthermore, long-term 2D culture was possible and HBV infection was confirmed. These results suggest that the stably and reproducibly providable HepaSH cells with lesser inter-individual differences in drug-metabolizing properties, may have a potential to substitution for PHH as practical standardized human hepatocytes in drug discovery research., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

Chimeric TK-NOG mice with a humanized liver (normal Hu-liver) are a unique animal model for predicting drug metabolism in humans. However, residual mouse hepatocytes occasionally prevent the precise evaluation of human drug metabolism. Herein, we developed a novel humanized liver TK-NOG mouse with a conditional knockout of liver-specific cytochrome P450 oxidoreductase (POR cKO Hu-liver). Immunohistochemical analysis revealed only a few POR-expressing cells around the portal vein in POR cKO mouse livers. NADPH-cytochrome c reductase and cytochrome P450 (P450)-mediated drug oxidation activity in liver microsomes from POR cKO mice was negligible. After the intravenous administration of S-warfarin, high circulating and urinary levels of S-7-hydroxywarfarin (a major human metabolite) were observed in POR cKO Hu-liver mice. Notably, the circulating and urinary levels of S-4'-hydroxywarfarin (a major warfarin metabolite in mice) were much lower in POR cKO Hu-liver mice than in normal Hu-liver mice. POR cKO Hu-liver mice with minimal interference from mouse hepatic P450 oxidation activity are a valuable model for predicting human drug metabolism., (© 2022. The Author(s).)

Pigs are an important species used in drug metabolism studies; however, the cytochromes P450 (P450s or CYPs) have not been fully investigated in pigs.In this study, pig CYP2C32, CYP2C33, CYP2C34, CYP2C36, CYP2C42, and CYP2C49 cDNAs were isolated and found to contain open reading frames of 490 or 494 amino acids that shared 64-82% sequence identity with human CYP2C8/9/18/19.Pig CYP2C genes formed a gene cluster in a genomic region that corresponded to that of the human CYP2C cluster; an additional gene cluster was formed by pig CYP2C33a and CYP2C33b distant from the first cluster but located in the same chromosome.Among the tissues analysed, these pig CYP2C mRNAs were preferentially expressed in liver, small intestine, and/or kidney; pig CYP2C49, CYP2C32, CYP2C34, and CYP2C33 mRNAs were the most abundant CYP2C mRNAs in liver, jejunum, ileum, and kidney, respectively.Metabolic assays showed that pig CYP2C proteins (heterologously expressed in Escherichia coli ) metabolised typical human CYP2C substrates diclofenac, warfarin, and/or omeprazole.The results suggest that these pig CYP2Cs are functional enzymes able to metabolise human CYP2C substrates in liver and small intestine, just as human CYP2Cs do.

Cynomolgus macaques are used in preclinical studies in part because of their evolutionary closeness to humans. However, drug transporters [including solute carrier (SLC) transporters] essential for the absorption and excretion of drugs have not been fully investigated at the molecular level in cynomolgus macaques. We identified and characterized cynomolgus macaque SLC15A1, SLC15A2, SLC22A1, SLC22A2, SLC22A6, SLC22A8, SLC47A1, and SLC47A2, along with SLCO (formerly SLC21A) transporters SLCO1A2, SLCO1B1, SLCO1B3, and SLCO2B1. These cynomolgus SLC transporters had high amino acid sequence identities (92-97%) with their human orthologs and contained sequence motifs characteristic of SLC transporters. Phylogenetic analysis showed that these cynomolgus SLC transporters were more closely clustered with their human orthologs than with those of dogs, rats, or mice. Gene structure and genomic organization were similar in macaques and humans. Cynomolgus SLC transporter mRNAs showed distinct tissue expression patterns, being most abundantly expressed in jejunum (SLC15A1), liver (SLC22A1, SLCO1B1, and SLCO2B1), and kidney (SLC15A2, SLC22A2, SLC22A6, SLC22A8, SLC47A1, SLC47A2, and SLCO1A2). In contrast, cynomolgus SLCO2B1 mRNA was more ubiquitously expressed. Among these SLC mRNAs, the most abundant in liver was SLCO1B1, in jejunum SLC15A1, and in kidney SLC22A2. These results suggest similar characteristics of SLC transporters in cynomolgus macaques and humans., Competing Interests: Declaration of competing interest The authors declare that there are no conflicts of interest., (Copyright © 2021 The Japanese Society for the Study of Xenobiotics. Published by Elsevier Ltd. All rights reserved.)

We developed a novel immunodeficient NOG mouse expressing HSVtk mutant clone 30 cDNA under the control of mouse transthyretin gene enhancer/promoter (NOG-TKm30) to acquire fertility in males and high inducibility of liver injury in females. Maximum human albumin levels (approx. 15 mg/mL plasma) in both male and female NOG-TKm30 mice engrafted with human hepatocytes (humanized liver mice) were observed 8-12 weeks after transplantation. Immunohistochemical analyses revealed abundant expression of major human cytochrome P450 (CYP) enzymes (CYP1A2, CYP2C9, CYP2D6, CYP2E1, and CYP3A4) in reconstituted liver with original zonal distribution. In vivo drug-drug interactions were observed in humanized liver mice as decreased area under the curve of midazolam (CYP3A4/5 substrate) and omeprazole (CYP3A4/5 and CYP2C19 substrate) after oral administration of rifampicin. Furthermore, we developed a pregnant model for evaluating prenatal exposure to drugs. The detection of thalidomide metabolites in the fetuses of pregnant humanized liver mice indicates that the novel TK model can be used for developmental toxicity studies requiring the assessment of human drug metabolism. These results suggest that the limitations of traditional TK-NOG mice can be addressed using NOG-TKm30 mice, which constitute a novel platform for humanized liver for both in vivo and in vitro studies., Competing Interests: Declarations of competing interest All authors declare that they have no conflicts of interest., (Copyright © 2021 The Japanese Society for the Study of Xenobiotics. Published by Elsevier Ltd. All rights reserved.)

Although human urinary aniline and 2,6-dimethylaniline were unexpectedly detected in biomonitoring data, little is known about the daily intake doses of aniline and 2,6-dimethylaniline in the living environment or their relation to tolerable daily intake (TDI) values in humans. In the current study, to evaluate the daily oral intake of aniline and 2,6-dimethylaniline in humans, forward and reverse dosimetry was carried out using simplified in silico physiologically based pharmacokinetic (PBPK) modeling established using in vivo experimental pharmacokinetic data. These data were from humanized-liver mice after single oral doses of 100 mg/kg aniline (previously determined) and 116 mg/kg 2,6-dimethylanine (currently investigated). The in vivo elimination rates of 2,6-dimethylaniline from plasma in humanized-liver mice were generally slow compared with those of aniline. Faster in vitro metabolic elimination rates of aniline mediated by liver 9000 × g supernatant fractions from rats than those from humans may suggest the existence of higher first-pass effects in rats than in humanized-liver mice. In silico aniline and 2,6-dimethylaniline concentration curves in human urine after virtual oral administrations were estimated by human PBPK models created with data from humanized-liver mice. Reverse dosimetry analysis using human PBPK models estimated the daily intake of aniline, based on reported human urinary concentrations in biomonitoring data, to be roughly similar to the aniline TDI level. These results suggest that forward and reverse dosimetry using simplified human PBPK models founded on data from humanized-liver mice can be used to evaluate possible higher than expected exposures of aniline and 2,6-dimethylaniline in humans.

Chimeric mice with humanized livers (humanized liver mice) are attractive experimental animal models for drug metabolism and pharmacokinetic studies. The "humanized liver" is a mature and functional liver with zonal position-specific expressions of human cytochrome P450 (P450) enzymes and a global gene expression pattern consistent with that of the mature human liver. Most P450-dependent drug oxidation activities were comparable between microsomes from livers of human and humanized liver mice based on similar expression levels of human P450 enzymes; however, some differences were observed between the two species, including considerable variations in activities of bufuralol 1'-hydroxylation and propafenone 4'-hydroxylation. Human disproportionate and/or unique metabolites of P450 substrate drugs were produced in humanized liver mice. Plasma concentration profiles of typical P450 substrate drugs in humans could be extrapolated from the corresponding data in humanized liver mice using simplified physiologically based pharmacokinetic modeling. Drug-drug interaction-mediated hepatic human CYP3A/2C induction by rifampicin (a human pregnane X receptor agonist) was observed in humanized liver mice. The major role of human CYP2C9 in in vivo diclofenac 4'-hydroxylation were determined using human CYP2C9-inactivated chimeric mice using a mechanism-based inhibitor, tienilic acid. Overall, based on the functional characteristics of hepatic human P450 enzymes, humanized liver mice are valuable experimental animals for studying metabolite profiling, pharmacokinetics, and drug interactions., Competing Interests: Conflict of interest The authors declare no conflict of interest., (Copyright © 2022 Elsevier Inc. All rights reserved.)

For drug development, species differences in drug-metabolism reactions present obstacles for predicting pharmacokinetics in humans. We characterized the species differences in hydrolases among humans and mice, rats, dogs, and cynomolgus monkeys. In this study, to expand the series of such studies, we attempted to characterize marmoset hydrolases. We measured hydrolase activities for 24 compounds using marmoset liver and intestinal microsomes, as well as recombinant marmoset carboxylesterase (CES) 1, CES2, and arylacetamide deacetylase (AADAC). The contributions of CES1, CES2, and AADAC to hydrolysis in marmoset liver microsomes were estimated by correcting the activities by using the ratios of hydrolase protein levels in the liver microsomes and those in recombinant systems. For six out of eight human CES1 substrates, the activities in marmoset liver microsomes were lower than those in human liver microsomes. For two human CES2 substrates and three out of seven human AADAC substrates, the activities in marmoset liver microsomes were higher than those in human liver microsomes. Notably, among the three rifamycins, only rifabutin was hydrolyzed by marmoset tissue microsomes and recombinant AADAC. The activities for all substrates in marmoset intestinal microsomes tended to be lower than those in liver microsomes, which suggests that the first-pass effects of the CES and AADAC substrates are due to hepatic hydrolysis. In most cases, the sums of the values of the contributions of CES1, CES2, and AADAC were below 100%, which indicated the involvement of other hydrolases in marmosets. In conclusion, we clarified the substrate preferences of hydrolases in marmosets. SIGNIFICANCE STATEMENT: This study confirmed that there are large differences in hydrolase activities between humans and marmosets by characterizing marmoset hydrolase activities for compounds that are substrates of human CES1, CES2, or arylacetamide deacetylase. The data obtained in this study may be useful for considering whether marmosets are appropriate for examining the pharmacokinetics and efficacies of new chemical entities in preclinical studies., (Copyright © 2021 by The American Society for Pharmacology and Experimental Therapeutics.)

Paraoxonase (PON) plays roles in the metabolism of organophosphate xenobiotics and drugs. Despite the importance of marmosets for research into drug metabolism and pharmacokinetics, marmoset paraoxonase has not yet been fully characterized. Consequently, we identified the PON1 gene in the marmoset genome by sequence homology analysis. Marmoset PON1 cDNA containing an open reading frame (1065 bp) was successfully cloned from marmoset liver by reverse transcription-polymerase chain reaction. The deduced amino acid sequence (355 amino acids) has approximately 93% identity with the human ortholog and contains important amino acid residues for substrate binding and calcium ion coordination. According to a phylogenetic tree of PON1 amino acid sequences constructed using data from seven animal species, marmoset PON1 is closer to human PON1 than it is to the PON1 orthologs of experimental animals such as pigs, rabbits, rats, and mice. Marmoset PON1 mRNA was predominantly expressed in liver among the five tissues examined. Marmoset PON1 protein secreted into plasma was detected by immunoblotting. The paraoxon-hydrolyzing activity in plasma was higher in marmosets than in humans. Based on these data, we concluded that marmoset and human PON1 have similar characteristics with regard to genomic structure, amino acid sequences, and tissue distribution., Competing Interests: Declaration of competing interest The authors declare that there are no conflicts of interest., (Copyright © 2021 The Japanese Society for the Study of Xenobiotics. Published by Elsevier Ltd. All rights reserved.)

p -Toluic acid, a metabolite of organic solvent xylene, has a high reported no-observed-effect level (NOEL, 1000 mg/kg) in rats, possibly because of direct glycine conjugation to methylhippuric acid. In this study, plasma levels of p -toluic acid and its glycine conjugate in mice and humanised-liver mice were evaluated after oral administrations.Although rapid conversion of p -toluic acid to its glycine conjugate was evident from mouse plasma concentrations, the biotransformation of p -toluic acid was slower in humanised-liver mice. The input parameters for physiologically based pharmacokinetic (PBPK) models were determined using fitting procedures to create PBPK-generated plasma concentration curves.The PBPK-modelled hepatic concentrations of p -toluic acid in humanised-liver mice were higher than those observed in plasma. PBPK-modelled hepatic and plasma concentrations of p -toluic acid also indicated slow elimination in humans.These results suggest that rapid conjugations of p -toluic acid reportedly observed in rats could result in overestimation of NOELs for conjugatable chemicals when extrapolated to humanised-liver mice or humans.

Flavin-containing monooxygenase (FMO) catalyzes the oxygenation of a wide variety of medicines and dietary-derived compounds. However, little information is available regarding drug interactions mediated by FMO3 in vivo. Consequently, we investigated interactions between FMO substrates in humanized-liver mice. Trimethylamine-d 9 and itopride were, respectively, intravenously and orally administered to humanized-liver mice (n = 5-7). The pharmacokinetic profiles of itopride (the victim drug) in the presence of trimethylamine (the perpetrator drug) were determined for 24 h after co-administration using liquid chromatography/tandem mass spectrometry. Itopride (10 mg/kg) was extensively oxygenated in humanized-liver mice to its N-oxide. The plasma concentrations of itopride N-oxide after co-administration of itopride and trimethylamine (10 and 100 mg/kg) were significantly suppressed in a dose-dependent manner, but only during the early phase, i.e., up to 2 h after co-administration. With the higher dose of trimethylamine, the areas under the concentration-time curves of itopride and its N-oxide significantly increased (1.6-fold) and decreased (to 60%), respectively; modeling suggested that these modified pharmacokinetics resulted from suppression of the in vivo hepatic intrinsic clearance (to 67%). These results suggest that food-derived trimethylamine may result in interactions with FMO drug substrates immediately after administration; however, the potential for this to occur in vivo may be limited., Competing Interests: Declaration of competing interest All authors declare that they have no conflicts of interest., (Copyright © 2020 The Japanese Society for the Study of Xenobiotics. Published by Elsevier Ltd. All rights reserved.)

Plasma concentrations of acetaminophen, its glucuronide and sulfate conjugates, and cysteinyl acetaminophen were experimentally determined after oral administrations of 10 mg/kg in humanised-liver mice, control mice, rats, common marmosets, cynomolgus monkeys, and minipigs; the results were compared with reported human pharmacokinetic data. Among the animals tested, only rats predominantly converted acetaminophen to sulfate conjugates, rather than glucuronide conjugates. In contrast, the values of area under the plasma concentration curves of acetaminophen, its glucuronide and sulfate conjugates, and cysteinyl acetaminophen after oral administration of acetaminophen in marmosets and minipigs were consistent with those reported in humans under the present conditions. Physiologically based pharmacokinetic (PBPK) models (consisting of the gut, liver, and central compartments) for acetaminophen and its primary metabolite could reproduce and estimate, respectively, the plasma and hepatic concentrations of acetaminophen in experimental animals and humans after single virtual oral doses. The values of area under the curves of hepatic concentrations of acetaminophen estimated using PBPK models were correlated with the measured levels of cysteinyl acetaminophen (a deactivated metabolite) in plasma fractions in these species. Consequently, using simple PBPK models and plasma data to predict hepatic chemical concentrations after oral doses could be helpful as an indicator of in vivo possible hepatotoxicity of chemicals such as acetaminophen.

Recently developed computational models can estimate plasma, hepatic, and renal concentrations of industrial chemicals in rats. Typically, the input parameter values (i.e., the absorption rate constant, volume of systemic circulation, and hepatic intrinsic clearance) for simplified physiologically based pharmacokinetic (PBPK) model systems are calculated to give the best fit to measured or reported in vivo blood substance concentration values in animals. The purpose of the present study was to estimate in silico these three input pharmacokinetic parameters using a machine learning algorithm applied to a broad range of chemical properties obtained from several cheminformatics software tools. These in silico estimated parameters were then incorporated into PBPK models for predicting internal exposures in rats. Following this approach, simplified PBPK models were set up for 246 drugs, food components, and industrial chemicals with a broad range of chemical structures. We had previously generated PBPK models for 158 of these substances, whereas 88 for which concentration series data were available in the literature were newly modeled. The values for the absorption rate constant, volume of systemic circulation, and hepatic intrinsic clearance could be generated in silico by equations containing between 14 and 26 physicochemical properties. After virtual oral dosing, the output concentration values of the 246 compounds in plasma, liver, and kidney from rat PBPK models using traditionally determined and in silico estimated input parameters were well correlated ( r ≥ 0.83). In summary, by using PBPK models consisting of chemical receptor (gut), metabolizing (liver), excreting (kidney), and central (main) compartments with in silico -derived input parameters, the forward dosimetry of new chemicals could provide the plasma/tissue concentrations of drugs and chemicals after oral dosing, thereby facilitating estimates of hematotoxic, hepatotoxic, or nephrotoxic potential as a part of risk assessment.

Tetrabromobisphenol A, a brominated flame retardant, is increasingly prevalent worldwide and presents a potential health risk. Adjusted animal biomonitoring equivalents of tetrabromobisphenol A after orally administered doses in humanized-liver mice were scaled up to humans using known species allometric scaling factors to set up simplified physiologically based pharmacokinetic (PBPK) models. Absorbed tetrabromobisphenol A was slightly, moderately, and extensively metabolized in vivo to its glucuronide in rats, control mice, and humanized-liver mice tested, respectively. In silico estimated hepatic exposures of tetrabromobisphenol A and its glucuronide generated using the rat PBPK model-generated plasma concentration profiles were consistent with the reported values. The extent of hepatic injury in humanized-liver mice caused by tetrabromobisphenol A was evaluated by detecting human albumin mRNA in mouse plasma after oral administration of a high dose of tetrabromobisphenol A (1000 mg/kg). Reverse dosimetry analyses were carried out using two human PBPK models (set up based on the humanized-liver-mouse model and by optimizing the input parameters for reported human plasma concentrations of tetrabromobisphenol A glucuronide) to estimate the tetrabromobisphenol A daily intake based on reported human serum concentrations of total tetrabromobisphenol A from biomonitoring data. Within the predictability of the forward and reverse dosimetry estimations, the calculated daily intake was found to be far below established health benchmark levels (i.e., the suggested daily reported reference dose) with a wide (4 orders of magnitude) safety margin. These results suggest that the simplified PBPK models can be successfully applied to forward and reverse dosimetry estimations of tissue and/or blood exposures of tetrabromobisphenol A in humans after oral doses.

The teratogenicity of the chemotherapeutic drug thalidomide is species-specific and affects humans, non-human primates, and rabbits. The primary oxidation of thalidomide in previously investigated rodents predominantly resulted in the formation of deactivated 5'-hydroxythalidomide. In the current study, similar in vivo biotransformations to 5-hydroxythalidomide and 5'-hydroxythalidomide were confirmed by the analysis of blood plasma from male rabbits, a thalidomide-sensitive species, after oral administration of thalidomide (2.0 mg/kg). Similar levels of thalidomide in seminal plasma and in blood plasma were detected using liquid chromatography-tandem mass spectrometry at 4 hr and 7 hr after oral doses in male rabbits. Seminal plasma concentrations of 5-hydroxythalidomide and 5'-hydroxythalidomide were also seen in male rabbits in a roughly similar time-dependent manner to those in the blood plasma after oral doses of thalidomide (2.0 mg/kg). Furthermore, the values generated by a simplified physiologically based pharmacokinetic rabbit model were in agreement with the measured in vivo blood plasma data under metabolic ratios of 0.01 for the hepatic intrinsic clearance of thalidomide to both unconjugated 5-hydroxythalidomide and 5'-hydroxythalidomide. These results suggest that metabolic activation of thalidomide may be dependent on rabbit liver enzymes just it was for cytochrome P450 enzymes in humanized-liver mice; in contrast, rodent livers predominantly mediate biotransformation of thalidomide to 5'-hydroxythalidomide. A developmental toxicity test system with experimental animals that involves intravaginal exposures to the chemotherapeutic drug thalidomide via semen should be considered in the future.

Troglitazone and its major metabolite troglitazone sulfate were intravenously administered to chimeric mice with different ratios of liver replacement by human hepatocytes. Total clearances were converted to hepatic intrinsic clearances normalized to their liver weight, with the assumption that extra-hepatic elimination of these compounds was negligible. These values were plotted against the replacement indices, and postulated values for virtual 100% chimeric mice were assumed to be equivalent to those in humans. Metabolic formation ratio was estimated by comparing AUCs of troglitazone sulfate after separate administration of troglitazone and troglitazone sulfate. Liver to plasma concentration ratios were obtained from direct measurement. These parameters were extrapolated to 100% chimeric mice and subjected to semi-physiological pharmacokinetic modeling using pharmacokinetic parameters for oral administration taken from literature. Our simulated plasma concentration-time profile of troglitazone agreed well with observed values obtained in clinical study. However, the profile of troglitazone sulfate was far below the reported values. Although the possible reasons for this discrepancy remains unsolved, the combination of chimeric mice with semi-physiological PK modeling proved to be a useful tool in understanding the function of each PK parameter in human pharmacokinetics of troglitazone and its conjugated metabolite., Competing Interests: Declaration of competing interest The authors declare no conflict of interest., (Copyright © 2020 The Japanese Society for the Study of Xenobiotics. Published by Elsevier Ltd. All rights reserved.)

Polymorphic human and cynomolgus macaque flavin-containing monooxygenases (FMO) 3 are important oxygenation enzymes for nitrogen-containing drugs. Inter-animal variability of FMO3-dependent drug oxygenations in vivo is suspected in cynomolgus macaques because such variability is evident in humans. Therefore, this follow-up study was performed to investigate the pharmacokinetics of orally administered deuterium-labeled trimethylamine in three cynomolgus macaques genotyped for FMO3. Trimethylamine-d 9 was rapidly absorbed and attained plasma concentrations greater than the background levels of non-labeled trimethylamine. Trimethylamine-d 9 was then converted to trimethylamine-d 9 N-oxide. The half-lives, maximum plasma concentrations, and areas under the curve for trimethylamine-d 9 and its N-oxygenated metabolite and the total clearance for orally administered trimethylamine-d 9 were not different among the heterozygote for Q506K FMO3, the heterozygote for V325I FMO3, and the heterozygote for both S99N and F510S FMO3. Trimethylamine N-oxygenation activities mediated by liver microsomes prepared from the same three animals were not substantially different. However, recombinant proteins of the corresponding cynomolgus FMO3 variants showed apparent reduced trimethylamine N-oxygenation activities compared with the wild-type proteins. This study suggests only limited polymorphic effects on the in vivo catalytic function of cynomolgus FMO3. These findings yield important insights in terms of both quantitative and qualitative variations of polymorphic FMO3 in cynomolgus liver., Competing Interests: Declaration of competing interest All authors declare that they have no conflicts of interest., (Copyright © 2020 The Japanese Society for the Study of Xenobiotics. Published by Elsevier Ltd. All rights reserved.)

The cynomolgus macaque is a nonhuman primate species that is often used in drug metabolism studies during drug development. However, the localization of UDP-glucuronosyltransferases (UGTs), essential drug-metabolizing enzymes, has not been fully investigated in the liver and small intestine of cynomolgus macaques. In this study, UGT activities were analyzed in liver (five lobes) and small intestine (the duodenum and six sections from the proximal jejunum to the distal ileum) using typical probe substrates of human UGTs: 7-hydroxycoumarin, estradiol, serotonin, propofol, and zidovudine. In liver, UGT activities with respect to all substrates were detected, and the activity levels were similar in all liver lobes of the cynomolgus macaques tested. In contrast, in the small intestine, UGT activities toward all substrates were detected, but their levels generally decreased from jejunum to ileum in cynomolgus macaques. The localization of estradiol 3-O-glucuronosyltransferases and serotonin O-glucuronosyltransferases (which are mainly UGT1A enzymes) appear to be different in liver and small intestine. These results collectively suggest that, in cynomolgus macaques, UGT1As are differentially localized in the small intestine but are relatively homogeneously distributed in the liver., Competing Interests: Declaration of competing interest The authors declare that there are no conflicts of interest., (Copyright © 2020 The Japanese Society for the Study of Xenobiotics. Published by Elsevier Ltd. All rights reserved.)

Prediction of human pharmacokinetics is important in the preclinical stage. Values for total clearance of compounds from plasma should be one of the most important pharmacokinetic parameters for predictions. Although several physiological and empirical methods including single-species allometry for prediction of values for human clearance of compounds using humanized-liver mice have been reported, further improvement of prediction accuracies would be still expected. To optimize these approaches, we proposed methods for unbound intrinsic clearance in virtually 100% humanized-liver mouse by incorporating unbound plasma fractions of compounds in differently humanized-liver mice. Comparisons of prediction accuracies of values for human clearance of 15 model compounds were performed among our current physiological and previously reported models and single-species allometry using humanized-liver mice. Incorporation of the actual unbound plasma fractions of compounds and correction of residual mice hepatocyte in humanized-liver mice showed comparable prediction accuracy to that by single-species allometry. After exclusion of 3 compounds with large species differences in values of clearance and unbound plasma fractions between mice and humans out of 15 compounds, prediction accuracies were improved in the methods investigated. The previously and present reported physiological methods could show the good prediction accuracy of values for clearance of drugs from plasma., Competing Interests: Declaration of competing interest The authors declare that there are no conflicts of interest., (Copyright © 2020 The Japanese Society for the Study of Xenobiotics. Published by Elsevier Ltd. All rights reserved.)

Recently developed high-throughput in vitro assays in combination with computational models could provide alternatives to animal testing. The purpose of the present study was to model the plasma, hepatic, and renal pharmacokinetics of approximately 150 structurally varied types of drugs, food components, and industrial chemicals after virtual external oral dosing in rats and to determine the relationship between the simulated internal concentrations in tissue/plasma and their lowest-observed-effect levels. The model parameters were based on rat plasma data from the literature and empirically determined pharmacokinetics measured after oral administrations to rats carried out to evaluate hepatotoxic or nephrotic potentials. To ensure that the analyzed substances exhibited a broad diversity of chemical structures, their structure-based location in the chemical space underwent projection onto a two-dimensional plane, as reported previously, using generative topographic mapping. A high-throughput in silico one-compartment model and a physiologically based pharmacokinetic (PBPK) model consisting of chemical receptor (gut), metabolizing (liver), central (main), and excreting (kidney) compartments were developed in parallel. For 159 disparate chemicals, the maximum plasma concentrations and the areas under the concentration-time curves obtained by one-compartment models and modified simple PBPK models were closely correlated. However, there were differences between the PBPK modeled and empirically obtained hepatic/renal concentrations and plasma maximal concentrations/areas under the concentration-time curves of the 159 chemicals. For a few compounds, the lowest-observed-effect levels were available for hepatotoxicity and nephrotoxicity in the Hazard Evaluation Support System Integrated Platform in Japan. The areas under the renal or hepatic concentration-time curves estimated using PBPK modeling were inversely associated with these lowest-observed-effect levels. Using PBPK forward dosimetry could provide the plasma/tissue concentrations of drugs and chemicals after oral dosing, thereby facilitating estimates of nephrotoxic or hepatotoxic potential as a part of the risk assessment.

Common marmoset (Callithrix jacchus) is an attractive animal model primate species for potential use in drug metabolism and pharmacokinetic studies. In this study, marmoset cytochrome P450 (P450) 2S1, 4V2, 7A1, 7B1, 8B1, 24A1, 26A1, 26C1, 27A1, 39A1, and 51A1 cDNAs were isolated from marmoset tissues (brains, lungs, livers, kidneys, and jejunums). Deduced amino acid sequences (89-98% homologous) of the marmoset P450 gene suggested similarity of molecular characteristics of marmoset P450s to human counterparts, compared with those of pig, rabbit, and rodents. Phylogenetic analysis using amino acid sequences indicated 11 marmoset P450 forms clustered with those of human and other primate counterparts, suggesting marmoset P450s have an evolutionary close relationship to human and other primate counterparts. Tissue expression patterns of these P450 mRNAs except for P450 7B1 mRNA were generally similar to those of human P450s in the five tissue types analyzed. These results suggest similarity of molecular characteristics for P450 2S1, 4V2, 7A1, 7B1, 8B1, 24A1, 26A1, 26C1, 27A1, 39A1, and 51A1 between marmosets and humans, in addition to the orthologs of human P450 1, 2, 3, and 4 families previously identified and characterized in marmosets., Competing Interests: Declaration of competing interest The authors declare that there are no conflicts of interest., (Copyright © 2019 The Japanese Society for the Study of Xenobiotics. Published by Elsevier Ltd. All rights reserved.)

1. In order to investigate the pharmacokinetics of diisononyl phthalate (DINP) in humans, we administered [phenyl-U- 14 C]DINP at a dose of 50.0 mg/kg orally to chimeric mice (humanized-liver mice) in which the liver of TK-NOG mice (control mice) was replaced with human hepatocytes. 2. The plasma radioactivity concentrations peaked (18.0 and 59.9 µg equivalent of DINP/mL, respectively) at 2 h after administration in control and humanized-liver mice. Concentrations rose again at 8 h in controls, but not in humanized-liver mice. 3. The cumulative excretion rates in urine and feces, respectively, were 58.1% and 37.3% of the doses in controls up to 48 h, but were 86.0% and 7.7% in humanized-liver mice. 4. The main circulating metabolites in control and humanized-liver mice were monoisononyl phthalate (MINP) and the glucuronide of oxidized MINP, respectively. The urinary excretion ratio of the glucuronide of oxidized MINP in control mice was one-third of that in humanized-liver mice. 5. The present results suggested that the oxidation rates of the primary metabolite of DINP and their excretion routes to urine/feces were different for control and humanized-liver mice. Species differences in liver activities could be a determinant factor in the in vivo metabolism and disposition of diallyl phthalates such as DINP.

Evaluation of uptake of lipophilic acid compounds into hepatocytes was an unresolved drug development issue because of their adsorption to cells and materials and low analytical sensitivity and accuracy in assessment of protein bindings. Uptake assays of compounds using hepatocytes suspended in serum were expected to solve these problems for prediction of in vivo hepatic clearance. Here, for compounds with high protein binding (>99%), diflunisal, montelukast, cerivastatin, telmisartan, fluvastatin and six new drug candidates, in vivo hepatic clearance predicted based on hepatic depletion and uptake (CL h, uptake, predicted ) data using hepatocytes in the absence and presence of sera was investigated. In vitro hepatic uptake results with hepatocytes suspended in serum improved prediction of human hepatic clearance values for highly lipophilic montelukast and telmisartan. In vivo CL h, uptake, predicted values of six new highly lipophilic acid drug candidates (protein binding >99.97%) and diflunisal, montelukast and cerivastatin predicted based on hepatocytes suspended in serum were within threefold differences of their total clearance in vivo in rats, guinea pigs or monkeys, except for montelukast in monkeys (5.8-fold). These results suggest that the human hepatic uptake in hepatocytes suspended in serum is useful for prediction of CL h, uptake, predicted , especially for highly lipophilic/protein binding acid compounds.

Diisononyl phthalate (DINP) used as a plasticizer is a mixture of compounds consisting of isononyl esters of phthalic acid. There are concerns about the bioaccumulation of such esters in humans. A [phenyl-U- 14 C]DINP mixture was synthesized and orally administered (50 mg/kg body weight) to control and humanized-liver mice and their pharmacokinetics were determined. Monoisononyl phthalate (MINP, a primary metabolite of DINP), oxidized MINP (isomers with hydroxy, carbonyl, and carboxy functional groups), and their glucuronides were detected in plasma from control and humanized-liver mice. Biphasic plasma concentration-time curves of MINP and its glucuronide were seen in control mice. In contrast, no such biphasic relationship was seen in humanized-liver mice, in which MINP and oxidized MINP were extensively excreted in the urine within 48 h. Animal biomonitoring equivalents of MINP and oxidized MINP from humanized-liver mice studies were scaled to human equivalents using known species allometric scaling factors with a simple physiologically based pharmacokinetic (PBPK) model. Estimated urinary oxidized MINP concentrations in humans were roughly consistent with reported concentrations of MINP (with a different side chain). The simplified PBPK model could estimate human urinary concentrations of MINP after ingestion of DINP and was capable of both forward and reverse dosimetry.

In this study, total body clearance (CL t ), volume of distribution at steady state (V ss ) and plasma concentration-time profiles in humans of model compounds were predicted using chimeric mice with humanized livers. On the basis of assumption that unbound intrinsic clearance (CL Uint ) per liver weight in chimeric mice was equal to those in humans, CL t were predicted by substituting human liver blood flow and liver weights in well-stirred model. V ss were predicted by Rodgers equation using scaling factors of tissue-plasma concentration ratios (SF Kp ) in chimeric mice estimated from a difference between the observed and predicted V ss . These physiological approaches showed high prediction accuracy for CL t and V ss values in humans. We compared the predictability of CL t and V ss determined by the physiologically based predictive approach using chimeric mice with those from predictive methods reported by Pharmaceutical Research Manufacturers of America. The physiological approach using chimeric mice indicated the best prediction accuracy in each predictive method. Simulation of human plasma concentration-time profiles were generally successful with physiologically based pharmacokinetic (PBPK) model incorporating CL Uint and SF Kp obtained from chimeric mice. Combined application of chimeric mice and PBPK modeling is effective for prediction of human PK in various compounds.

The common marmoset (Callithrix jacchus) is a useful experimental animal to evaluate the pharmacokinetics of drug candidates. Cytochrome P450 (P450) 2B enzyme in marmoset livers has been identified; however, only limited information on the enzymatic properties and distribution has been available. Marmoset P450 2B6 amino acids showed high sequence identities (>86%) with those of primates including humans and cynomolgus monkeys. Phylogenetic analysis using amino acid sequences indicated that marmoset P450 2B6 was closer to human and cynomolgus monkey P450 2B6 than to P450 2B orthologs of other species, including pigs, dogs, rabbits and rodents. Quantitative polymerase chain reaction analysis using specific primers showed P450 2B6 mRNA predominantly expressed in livers among the five marmoset tissues, similar to those of humans and cynomolgus monkeys. Marmoset P450 2B6 heterologously expressed in Escherichia coli membranes oxidized 7-ethoxycoumarin, pentoxyresorufin, propofol and testosterone, at roughly similar rates to those of humans and/or cynomolgus monkeys. A high capacity of marmoset P450 2B6 with propofol 4-hydroxylation (at low ionic strength conditions) with a low K m value was relatively comparable to that for marmoset livers. These results collectively indicated a high propofol 4-hydroxylation activity of P450 2B6 expressed in marmoset livers.

Dibutyl phthalate (DBP) was widely used as a plasticizer but it has been recently replaced with other kinds of phthalates such as di(2-ethylhexyl)phthalate and diisononyl phthalate because of its toxicity. To evaluate the human risk of DBP, forward and reverse dosimetry was conducted using in silico simplified physiologically based pharmacokinetic (PBPK) modeling based on in vivo experimental pharmacokinetic data in humanized-liver mice (HL-mice) obtained after an oral dose of 100 mg/kg. Absorbed DBP was converted to monobutyl phthalate (MBP) and its glucuronide extensively in vivo. HL-mice had higher concentrations of MBP glucuronide in plasma than did the control mice. Concentrations of MBP glucuronide in 0-7 h accumulated urine samples from HL-mice were significantly higher than those in control mice. Similarly, in vitro MBP glucuronidation rates mediated by pooled microsomes from rat or mouse livers were lower than those mediated by human liver microsomes. Liver damage by MBP to humanized liver was detected by measuring human albumin mRNA in HL-mouse plasma. By simple PBPK modeling, in silico concentration curves in plasma, liver, or urine following virtual oral administration of DBP were created for rats, control mice, and HL-mice. A human PBPK model for MBP was established based on the HL-mouse PBPK model using allometric scaling without consideration of interspecies factors in terms of liver metabolism. Human PBPK models were used to estimate urinary and plasma concentrations of MBP and its glucuronide throughout 14 days of oral DBP administration (1.2 and 13 μg/kg/day). Reverse dosimetry PBPK modeling found that reported 50th and 95th percentile MBP urine and plasma concentrations of the general population could potentially imply exposures similar to or exceeding tolerable daily intake levels (5-10 μg/kg/day) recommended by the European and Japanese authorities. Further in-depth assessment of DBP is needed to assess the validity of assumptions made based on human biomonitoring data.

Long-term exposure to certain volatile organic compounds is a significant public health concern. A variety of food containers and drinking cups prepared from polystyrene or polystyrene-related plastics could contain styrene monomer. In the current study, the concentrations of styrene in plasma and liver were surveyed and determined after oral doses of 25 mg/kg to rats and 200 mg/kg to control and humanized-liver mice. Plasma concentrations of styrene in rats were still detected 2 hr after 10-25 mg/kg oral doses. In contrast, after an order of magnitude higher oral dose of styrene (200 mg/kg) to mice, styrene in mouse plasma was rapidly cleared within 15 min to the limit-of-detection level. However, unmetabolized styrene was detected in mouse liver 24 hr after oral treatment. A simple physiologically based pharmacokinetic (PBPK) model capable of estimating blood and liver concentrations of styrene was established for rats. A human PBPK model was then set up for styrene by using the same intrinsic hepatic clearances in rats and humans and by applying allometric scaling to rat parameters obtained from the plasma concentrations of styrene in rats. By reverse dosimetry analysis (from concentrations to doses), we found that the 95th percentile values of styrene concentrations (0.132 ng/mL) reported in United States biomonitoring data of more than 1000 human blood samples may imply exposure to repeated oral doses of styrene of 2.89 µg/kg/day. These results suggest that styrene biomonitoring data in human blood samples imply exposures roughly similar to or lower than the established tolerable daily intake level of 7.7 μg/kg/day.

Background: Common marmosets (Callithrix jacchus) are potentially useful nonhuman primate models for preclinical studies. Information for major drug-metabolizing cytochrome P450 (P450) enzymes is now available that supports the use of this primate species as an animal model for drug development. Here, we collect and provide an overview of information on the activities of common marmoset hepatic and intestinal microsomes with respect to 28 typical human P450 probe oxidations., Results: Marmoset P450 2D6/8-dependent R-metoprolol O-demethylation activities in hepatic microsomes were significantly correlated with those of midazolam 1'- and 4-hydroxylations, testosterone 6β-hydroxylation, and progesterone 6β-hydroxylation, which are probe reactions for marmoset P450 3A4/5/90. In marmosets, the oxidation activities of hepatic microsomes and intestinal microsomes were roughly comparable for midazolam and terfenadine. Overall, multiple forms of marmoset P450 enzymes in livers and intestines had generally similar substrate recognition functionalities to those of human and/or cynomolgus monkey P450 enzymes., Conclusion: The marmoset could be a model animal for humans with respect to the first-pass extraction of terfenadine and related substrates. These findings provide a foundation for understanding individual pharmacokinetic and toxicological results in nonhuman primates as preclinical models and will help to further support understanding of the molecular mechanisms of human P450 function., (Copyright© Bentham Science Publishers; For any queries, please email at epub@benthamscience.net.)

1. Simulated clearances of R-warfarin and efavirenz were recently reported for individual cynomolgus monkeys genotyped for cytochrome P450 2C19 and 2C9, respectively. To expand and verify this modeling procedure, simulations of R/S-omeprazole and R/S-warfarin clearances after oral administrations in individual marmosets were performed using individual simplified physiologically based pharmacokinetic (PBPK) modeling consisting of gut, liver and central compartments. 2. Pharmacokinetics of R/S-omeprazole were chirally determined using the previously reported plasma microsamples in this study. The areas under the plasma concentration/time curves (AUC) of R-omeprazole and S-warfarin, but not S-omeprazole and R-warfarin, after oral administrations in the P450 2C19 homozygous mutant group were significantly higher than those in the wild-type group. These modeled hepatic intrinsic clearances were also significantly associated with the marmoset P450 2C19 genotypes. Other parameter values, e.g. absorption rate constants or systemic circulation volumes, were not likely determining factors. 3. The reported individual AUC values measured in 4-6 marmosets after oral R-omeprazole and S-warfarin administrations were significantly correlated with the AUC values predicted using the PBPK models after virtual administrations. 4. This study indicates that clearances of R-omeprazole, S-warfarin and related medicines associated with polymorphic P450 2C19 in individual marmosets can be simulated using simplified individual PBPK models.

1. A potentially useful animal model for preclinical studies is the common marmoset (Callithrix jacchus). In this study, using reverse-transcription polymerase chain reaction from marmoset livers, we identified a novel cytochrome P450 (P450) 2F1 cDNA with an open reading frame of 1473 bp. 2. High sequence identities of 92-94% with primate P450 2 F amino acid sequences were indicated by deduced amino acid sequences of P450 2F1 cDNA. Phylogenetic analysis indicates that marmoset P450 2F1 is more congruent with primate P450 2 F forms than those of other species such as rodents. 3. Among five tissue types examined, abundant expression of marmoset P450 2F1 mRNA and P450 2F1 protein in lungs was shown. Cynomolgus monkey P450 2F1 mRNA was abundantly expressed in lungs as well as testes and ovaries in 10 tissue types. 4. Similar to those of humans and cynomolgus monkeys, marmoset P450 2F1 heterologously expressed in Escherichia coli membranes efficiently catalyzed 7-ethoxycoumarin O-deethylation and biphenyl hydroxylation, however unlike human P450 2F1, marmoset P450 2F1 exhibited hydroxylation activity toward coumarin and chlorzoxazone. 5. These findings indicated that P450 2F1 enzyme expressed in marmoset lungs and also catalyzed metabolism of xenobiotics, suggesting the importance of P450 2 F-dependent drug metabolism in marmoset lungs.

1. Roles of human cytochrome P450 (P450) 3A4 in oxidation of an antihistaminic drug terfenadine have been previously investigated in association with terfenadine-ketoconazole interaction. Several antihistamine drugs have been recently identified as substrates for multiple P450 enzymes. In this study, overall roles of P450 3A4, 2J2, and 4F12 enzymes in terfenadine t-butyl hydroxylation were investigated in small intestines and livers from humans, marmosets, and/or cynomolgus monkeys. 2. Human liver microsomes and liver and small intestine microsomes from marmosets and cynomolgus monkeys effectively mediated terfenadine t-butyl hydroxylation. Ketoconazole and N-hydroxy-N'-(4-butyl-2-methylphenyl)-formamidine (a P450 4A/F inhibitor) almost completely and moderately inhibited these activities, respectively, in human liver microsomes; however, these chemicals did not show substantially suppression in marmoset liver. Anti-human P450 3A and 4F antibodies showed the roughly supportive inhibitory effects. 3. Recombinant P450 3A4/90 and 4F12 showed high terfenadine t-butyl hydroxylation activities with substrate inhibition constants of 84-144 μM (under 26-76 μM of K m values), in similar manners to liver and intestine microsomes. 4. These results suggest that human and marmoset P450 3A4/90 and 4F12 in livers or small intestines played important roles in terfenadine t-butyl hydroxylation. Marmosets could be a model for humans during first pass extraction of terfenadine and related substrates.

Lenalidomide has been shown to be potentially teratogenic in thalidomide-sensitive animal species. Screening for thalidomide analogs devoid of teratogenicity/toxicity-attributable to drug metabolism and disposition, but having immunomodulatory properties-is a strategic pathway towards development of new anticancer drugs. Plasma concentrations of lenalidomide were investigated in immunodeficient control and humanized-liver mice following oral administration of lenalidomide (50 mg/kg). Plasma concentrations of lenalidomide (1-2 hr after administration) were slightly but significantly higher in humanized-liver mice than in control mice (p < 0.05). Human albumin mRNA, a liver-specific toxicity marker, was found in the blood of humanized-liver mice 24 hr after lenalidomide administration. Simulations of human plasma concentrations of lenalidomide were achieved with simplified physiologically-based pharmacokinetic models in control and humanized-liver mice or by the direct fitting analysis of reported human data, in accordance with reported lenalidomide concentrations after low dose administration in humans. The results indicate that pharmacokinetic profiles of lenalidomide, a compound resulting from introducing one aromatic amino group into thalidomide and removing one keto group, resulted in less species variation in in vivo pharmacokinetics in control and humanized-liver mice and that immunodeficient humanized-liver mice can serve as experimental model animals for human liver injury in drug development at high doses, with human albumin RNA analysis in plasma.

Cytochromes P450 (P450) largely remain to be characterized in great apes. Comparative immunochemical detection of drug metabolizing forms of P450s 1A, 2A, 2B, 2C, 2D, 2E, 2J, 3A, 4A, and 4F in liver microsomes from chimpanzees, gorillas, orangutans, gibbons, cynomolgus and rhesus macaques, and common marmosets were carried out., (© 2017 John Wiley & Sons A/S. Published by John Wiley & Sons Ltd.)

A β -blocker, metoprolol, is one of the in vivo probes for human cytochrome P450 (P450) 2D6. Investigation of nonhuman primate P450 enzymes helps to improve the accuracy of the extrapolation of pharmacokinetic data from animals into humans. Common marmosets ( Callithrix jacchus ) are a potential primate model for preclinical research, but the detailed roles of marmoset P450 enzymes in metoprolol oxidation remain unknown. In this study, regio- and stereo-selectivity of metoprolol oxidations by a variety of P450 enzymes in marmoset and human livers were investigated in vitro. Although liver microsomes from cynomolgus monkeys and rats preferentially mediated S -metoprolol O -demethylation and R -metoprolol α -hydroxylation, respectively, those from humans, marmosets, minipigs, and dogs preferentially mediated R -metoprolol O -demethylation, in contrast to the slow rates of R - and S -metoprolol oxidation in mouse liver microsomes. R - and S -metoprolol O -demethylation activities in marmoset livers were strongly inhibited by quinidine and ketoconazole, and were significantly correlated with bufuralol 1'-hydroxylation and midazolam 1'-hydroxylation activities and also with P450 2D and 3A4 contents, which is different from the case in human livers that did not have any correlations with P450 3A-mediated midazolam 1'-hydroxylation. Recombinant human P450 2D6 enzyme and marmoset P450 2D6/3A4 enzymes effectively catalyzed R -metoprolol O -demethylation, comparable to the activities of human and marmoset liver microsomes, respectively. These results indicated that the major roles of P450 2D enzymes for the regio- and stereo-selectivity of metoprolol oxidation were similar between human and marmoset livers, but the minor roles of P450 3A enzymes were unique to marmosets., (Copyright © 2017 by The American Society for Pharmacology and Experimental Therapeutics.)

1. A cDNA encoding novel cytochrome P450 (P450) 4A enzyme was cloned from marmoset livers by reverse transcription (RT)-polymerase chain reaction (PCR) based on the marmoset genome sequences. The amino acid sequence deduced from P450 4A11 cDNA contained consensus sequences of six substrate recognition sites and one heme-binding domain. 2. Marmoset P450 4A11, highly identical (85-88%) to cynomolgus monkey and human P450 4A enzymes, was grouped into the same cluster as cynomolgus monkey and human P450 4A enzymes by phylogenetic analysis. 3. The tissue distribution analyses by real-time RT PCR and immunoblotting demonstrated that marmoset P450 4A11 mRNA and proteins were expressed in kidneys and livers. Marmoset P450 4A11 enzyme heterologously expressed in Escherichia coli preferentially catalyzed the ω-hydroxylation of arachidonic acid and lauric acid, similar to cynomolgus monkey and human P450 4A11 enzymes. However, lauric acid ω-hydroxylation activity of marmoset P450 4A11 was low compared with those of marmoset liver microsomes. 4. These results indicated that novel marmoset P450 4A11 was also a fatty acid ω-hydroxylase expressed in kidneys and livers, with the same regioselectivity (at ω-position of fatty acid) as cynomolgus monkey and human P450 4A enzymes.

Common marmosets ( Callithrix jacchus ) are potentially primate models for preclinical drug metabolism studies because there are similarities in the molecular characteristics of cytochrome P450 enzymes between this species and humans. However, characterization of non-cytochrome P450 enzymes has not been clarified in marmosets. Here, we report characterization of flavin-containing monooxygenases FMO1-FMO5 identified in marmoset tissues. Marmoset FMO forms shared high amino acid sequence identities (93%-95%) and phylogenetic closeness with human homologous FMO forms. FMO1 and FMO3 mRNA were abundantly expressed in the liver and kidneys among five marmoset tissues examined, where FMO3 protein was detected by immunoblotting. FMO inhibition assays using preheated tissue microsomes indicated that benzydamine N -oxygenation and sulindac sulfide S -oxygenation in the marmoset liver was mainly catalyzed by FMO3, the major hepatic FMO. Marmoset FMO3 protein heterologously expressed in Escherichia coli effectively catalyzed benzydamine N -oxygenation and sulindac sulfide S -oxygenation comparable to marmoset liver microsomes. These results indicate that the FMO3 enzyme expressed in marmoset livers mainly metabolizes benzydamine and sulindac sulfide (typical human FMO substrates), suggesting its importance for FMO-dependent drug metabolism in marmosets., (Copyright © 2017 by The American Society for Pharmacology and Experimental Therapeutics.)

Common marmosets ( Callithrix jacchus ), small New World primates, are increasingly attracting attention as potentially useful animal models for drug development. However, characterization of cytochrome P450 (P450) 3A enzymes involved in the metabolism of a wide variety of drugs has not investigated in marmosets. In this study, sequence homology, tissue distribution, and enzymatic properties of marmoset P450 3A4 ortholog, 3A5 ortholog, and 3A90 were investigated. Marmoset P450 3A forms exhibited high amino acid sequence identities (88-90%) to the human and cynomolgus monkey P450 3A orthologs and evolutionary closeness to human and cynomolgus monkey P450 3A orthologs compared with other P450 3A enzymes. Among the five marmoset tissues examined, P450 3A4 ortholog mRNA was abundant in livers and small intestines where P450 3A4 ortholog proteins were immunologically detected. Three marmoset P450 3A proteins heterologously expressed in Escherichia coli membranes catalyzed midazolam 1'- and 4-hydroxylation, alprazolam 4-hydroxylation, nifedipine oxidation, and testosterone 6 β -hydroxylation, similar to cynomolgus monkey and human P450 3A enzymes. Among the marmoset P450 3A enzymes, P450 3A4 ortholog effectively catalyzed midazolam 1'-hydroxylation, comparable to microsomes from marmoset livers and small intestines. Correlation analyses with 23 individual marmoset liver microsomes suggested contributions of P450 3A enzymes to 1'-hydroxylation of both midazolam (human P450 3A probe) and bufuralol (human P450 2D6 probe), similar to cynomolgus monkey P450 3A enzymes. These results indicated that marmoset P450 3A forms had functional characteristics roughly similar to cynomolgus monkeys and humans in terms of tissue expression patterns and catalytic activities, suggesting marmosets as suitable animal models for P450 3A-dependent drug metabolism., (Copyright © 2017 by The American Society for Pharmacology and Experimental Therapeutics.)

Sex and age differences in hepatic expression of drug-metabolizing enzyme genes could cause variations in drug metabolism, but has not been fully elucidated, especially in Asian population. In this study, the global expression of human hepatic genes was analyzed by microarrays in 40 Japanese subjects (27 males and 13 females). Thirty-five sex-biased genes were identified (P < 0.005). Whereas, 60 age-biased genes in two age groups, <60 years and ≥70 years (P < 0.001), were identified in males. By Gene Ontology analysis, the sex-biased genes were related to protein catabolism and modification, while the age-biased genes were related to transcription regulation and cell death. Quantitative polymerase chain reaction confirmed the female-biased expression of drug-metabolizing enzyme genes BChE, CYP4X1, and SULT1E1 (≥1.5-fold, P < 0.05). Further analysis of drug-metabolizing enzyme genes indicated that expression of CYP2A6 and CYP3A4 in females in the ≥70 age group was less than in the <60 age group (≥1.5-fold, P < 0.05), and this trend was also observed for PXR expression in males (≥1.5-fold, P < 0.05). The results presented provide important insights into hepatic physiology and function, especially drug metabolism, with respect to sex and age., (Copyright © 2016 The Japanese Society for the Study of Xenobiotics. Published by Elsevier Ltd. All rights reserved.)

Species differences exist in terms of drug oxidation activities, which are mediated mainly by cytochrome P450 (P450) enzymes. To overcome the problem of species extrapolation, transchromosomic mice containing a human P450 3A cluster or chimeric mice transplanted with human hepatocytes have been introduced into the human toxicology research area. In this review, drug metabolism and disposition mediated by humanized livers in chimeric mice are summarized in terms of biliary/urinary excretions of phthalate and bisphenol A and plasma clearances of the human cocktail probe drugs caffeine, warfarin, omeprazole, metoprolol, and midazolam. Simulation of human plasma concentrations of the teratogen thalidomide and its human metabolites is possible with a simplified physiologically based pharmacokinetic model based on data obtained in chimeric mice, in accordance with reported clinical thalidomide concentrations. In addition, in vivo nonspecific hepatic protein binding parameters of metabolically activated 14 C-drug candidate and hepatotoxic medicines in humanized liver mice can be analyzed by accelerator mass spectrometry and are useful for predictions in humans.

Marmoset cytochrome P450 2C19, highly homologous to human P450 2C9 and 2C19, has been identified in common marmosets (Callithrix jacchus), a nonhuman primate species used in drug metabolism studies. Although genetic variants in human and macaque P450 2C genes account for the interindividual variability in drug metabolism, genetic variants have not been investigated in the marmoset P450 2C19 In this study, sequencing of P450 2C19 in 24 marmosets identified three variants p.[(Phe7Leu; Ser254Leu; Ile469Thr)], which showed substantially reduced metabolic capacity of S-warfarin compared with the wild-type group in vivo and in vitro. Although mean plasma concentrations of R-warfarin in marmosets determined after chiral separation were similar between the homozygous mutant and wild-type groups up to 24 hours after the intravenous and oral administrations of racemic warfarin, S-warfarin depletion from plasma was significantly faster in the three wild-type marmosets compared with the three homozygous mutant marmosets. These variants, cosegregating in the marmosets analyzed, influenced metabolic activities in 18 marmoset liver microsomes because the homozygotes and heterozygotes showed significantly reduced catalytic activities in liver microsomes toward S-warfarin 7-hydroxylation compared with the wild-type group. Kinetic analysis for S-warfarin 7-hydroxylation indicated that the recombinant P450 2C19 Ser254Leu variant would change the metabolic capacity. These results indicated that the interindividual variability of P450 2C-dependent drug metabolism such as S-warfarin clearance is at least partly accounted for by P450 2C19 variants in marmosets, suggesting that polymorphic P450 2C-dependent catalytic functions are relatively similar between marmosets and humans., (Copyright © 2016 by The American Society for Pharmacology and Experimental Therapeutics.)

Common marmosets (Callithrix jacchus) are attracting attention as animal models in preclinical studies for drug development. However, cytochrome P450s (P450s), major drug-metabolizing enzymes, have not been fully identified and characterized in marmosets. In this study, based on the four novel P450 4F genes found on the marmoset genome, we successfully isolated P450 4F2, 4F3B, 4F11, and 4F12 cDNAs in marmoset livers. Deduced amino acid sequences of the four marmoset P450 4F forms exhibited high sequence identities (87%-93%) to the human and cynomolgus monkey P450 4F homologs. Marmoset P450 4F3B and 4F11 mRNAs were predominantly expressed in livers, whereas marmoset P450 4F2 and 4F12 mRNAs were highly expressed in small intestines and livers. Four marmoset P450 4F proteins heterologously expressed in Escherichia coli catalyzed the ω-hydroxylation of leukotriene B4 In addition, marmoset P450 4F12 effectively catalyzed the hydroxylation of antiallergy drug ebastine, a human P450 2J/4F probe substrate. Ebastine hydroxylation activities by small intestine and liver microsomes from marmosets and cynomolgus monkeys showed greatly higher values than those of humans. Ebastine hydroxylation activities by marmoset and cynomolgus monkey small intestine microsomes were inhibited (approximately 60%) by anti-P450 4F antibodies, unlike human small intestine microsomes, suggesting that contribution of P450 4F enzymes for ebastine hydroxylation in the small intestine might be different between marmosets/cynomolgus monkeys and humans. These results indicated that marmoset P450 4F2, 4F3B, 4F11, and 4F12 were expressed in livers and/or small intestines and were functional in the metabolism of endogenous and exogenous compounds, similar to those of cynomolgus monkeys and humans., (Copyright © 2016 by The American Society for Pharmacology and Experimental Therapeutics.)

The microminipig is one of the smallest minipigs that has emerged as a possible experimental animal model, because it shares many anatomical and/or physiological similarities with humans, including the coronary artery distribution in the heart, the digestive physiology, the kidney size and its structure, and so on. However, information on gene expression profiles, including those on drug-metabolizing phase I and II enzymes, in the microminipig is limited. Therefore, the aim of the present study was to identify transcripts in microminipig livers and to determine gene expression profiles. De novo assembly and expression analyses of microminipig transcripts were conducted with liver samples from three male and three female microminipigs using parallel long-read and short-read sequencing technologies. After unique sequences had been automatically aligned by assembling software, the mean contig length of 50843 transcripts was 707 bp. The expression profiles of cytochrome P450 (P450) 1A2, 2C, 2E1 and 3A genes in livers in microminipigs were similar to those in humans. Liver carboxylesterase (CES) precursor, liver CES-like, UDP-glucuronosyltransferase (UGT) 2C1-like, amine sulfotransferase (SULT)-like, N-acetyltransferases (NAT8) and glutathione S-transferase (GST) A2 genes, which are relatively unknown genes in pigs and/or humans, were expressed strongly. Furthermore, no significant gender differences were observed in the gene expression profiles of phase I enzymes, whereas UGT2B17, SULT1E1, SULT2A1, amine SULT-like, NAT8 and GSTT4 genes were different between males and females among phase II enzyme genes under the present sample conditions. These results provide a foundation for mechanistic studies and the use of microminipigs as model animals for drug development in the future. Copyright © 2016 John Wiley & Sons, Ltd., (Copyright © 2016 John Wiley & Sons, Ltd.)

Macaques, including cynomolgus and rhesus macaques, are important animal species used in drug metabolism studies. CYP2D44 is expressed in cynomolgus macaque liver and encodes a functional drug metabolizing enzyme, metabolizing typical human CYP2D substrates such as bufuralol and dextromethorphan. CYP2D44 is highly homologous to human CYP2D6 that is known to be polymorphic with a large inter-individual variation in metabolic activities, however, genetic polymorphisms have not been investigated in macaque CYP2D44. In the present study, screening of 78 cynomolgus and 40 rhesus macaques found a total of 67 variants, including 64 non-synonymous variants, 1 nonsense mutation, and 2 frameshift mutations, and 1 gene conversion, of which 14, 19, and 15 variants were unique to Indochinese cynomolgus macaques, Indonesian cynomolgus macaques, and Chinese rhesus macaques, respectively. Eleven of the 64 non-synonymous variants were located in substrate recognition sites, the regions important for protein function. By site-directed mutagenesis and metabolic assays, S175N, V185L, A235G, R242G, R245K, and N337D showed substantially decreased activity in bufuralol 1'-hydroxylation as compared with wild-type proteins. Moreover, two null alleles (c.128T>del and c.664G>T) were found in Indonesian cynomolgus macaques, but not in Indochinese cynomolgus macaques or Chinese rhesus macaques. These results suggest that genetic polymorphisms might account for the variability of CYP2D44-dependent metabolism in macaques.

The aim of this study was to extrapolate to humans the pharmacokinetics of estrogen analog bisphenol A determined in chimeric mice transplanted with human hepatocytes. Higher plasma concentrations and urinary excretions of bisphenol A glucuronide (a primary metabolite of bisphenol A) were observed in chimeric mice than in control mice after oral administrations, presumably because of enterohepatic circulation of bisphenol A glucuronide in control mice. Bisphenol A glucuronidation was faster in mouse liver microsomes than in human liver microsomes. These findings suggest a predominantly urinary excretion route of bisphenol A glucuronide in chimeric mice with humanized liver. Reported human plasma and urine data for bisphenol A glucuronide after single oral administration of 0.1mg/kg bisphenol A were reasonably estimated using the current semi-physiological pharmacokinetic model extrapolated from humanized mice data using algometric scaling. The reported geometric mean urinary bisphenol A concentration in the U.S. population of 2.64μg/L underwent reverse dosimetry modeling with the current human semi-physiological pharmacokinetic model. This yielded an estimated exposure of 0.024μg/kg/day, which was less than the daily tolerable intake of bisphenol A (50μg/kg/day), implying little risk to humans. Semi-physiological pharmacokinetic modeling will likely prove useful for determining the species-dependent toxicological risk of bisphenol A., (Copyright © 2015 Elsevier Inc. All rights reserved.)