Machine Learning Prediction of the Three Main Input Parameters of a Simplified Physiologically Based Pharmacokinetic Model Subsequently Used to Generate Time-Dependent Plasma Concentration Data in Humans after Oral Doses of 212 Disparate Chemicals.

Physiologically based pharmacokinetic (PBPK) modeling has the potential to play significant roles in estimating internal chemical exposures. The three major PBPK model input parameters (i.e., absorption rate constants, volumes of the systemic circulation, and hepatic intrinsic clearances) were generated in silico for 212 chemicals using machine learning algorithms. These input parameters were calculated based on sets of between 17 and 65 chemical properties that were generated by in silico prediction tools before being processed by machine learning algorithms. The resulting simplified PBPK models were used to estimate plasma concentrations after virtual oral administrations in humans. The estimated absorption rate constants, volumes of the systemic circulation, and hepatic intrinsic clearance values for the 212 test compounds determined traditionally (i.e., based on fitting to measured concentration profiles) and newly estimated had correlation coefficients of 0.65, 0.68, and 0.77 (p < 0.01, n = 212), respectively. When human plasma concentrations were modeled using traditionally determined input parameters and again using in silico estimated input parameters, the two sets of maximum plasma concentrations (r = 0.85, p < 0.01, n = 212) and areas under the curve (r = 0.80, p < 0.01, n = 212) were correlated. Virtual chemical exposure levels in liver and kidney were also estimated using these simplified PBPK models along with human plasma levels. These results indicate that the PBPK model input parameters for humans of a diverse set of compounds can be reliability estimated using chemical descriptors calculated using in silico tools.

In Silico Prediction of Input Parameters for Simplified Physiologically Based Pharmacokinetic Models for Estimating Plasma, Liver, and Kidney Exposures in Rats after Oral Doses of 246 Disparate Chemicals.

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.

An Updated In Silico Prediction Method for Volumes of Systemic Circulation of 323 Disparate Chemicals for Use in Physiologically Based Pharmacokinetic Models to Estimate Plasma and Tissue Concentrations after Oral Doses in Rats.

Updated algorithms for predicting the volumes of systemic circulation (V1), along with absorption rate constants and hepatic intrinsic clearances, as input parameters for physiologically based pharmacokinetic (PBPK) models were established to improve the accuracy of estimated plasma and tissue concentrations of 323 chemicals after virtual oral administrations in rats. Using ridge regression with an enlarged set of chemical descriptors (up to 99), the estimated input V1 values resulted in an improved correlation coefficient (from 246 compounds) with the traditionally determined values. The PBPK model input parameters for rats of diverse compounds can be precisely estimated by increasing the number of descriptors.

Metabolic Profiles of Tetrabromobisphenol A in Humans Extrapolated from Humanized-Liver Mouse Data Using a Simplified Physiologically Based Pharmacokinetic Model.

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.

Differences in Pharmacokinetics and Haematotoxicities of Aniline and Its Dimethyl Derivatives Orally Administered in Rats.

Aniline and its dimethyl derivatives reportedly become haematotoxic after metabolic N-hydroxylation of their amino groups. The plasma concentrations of aniline and its dimethyl derivatives after single oral doses of 25 mg/kg in rats were quantitatively measured and semi-quantitatively estimated using LC-tandem mass spectrometry. The quantitatively determined elimination rates of aniline; 2,4-dimethylaniline; and 3,5-dimethylaniline based on rat plasma versus time curves were generally rapid compared with those of 2,3-; 2,5-; 2,6-; and N,2-dimethylaniline. The primary acetylated metabolites of aniline; 2,4-dimethylaniline; and 3,5-dimethylaniline, as semi-quantitatively estimated based on their peak areas in LC analyses, were more extensively formed than those of 2,3-; 2,5-; 2,6-; and N,2-dimethylaniline. The areas under the curve of unmetabolized (remaining) aniline and its dimethyl derivatives estimated using simplified physiologically based pharmacokinetic models (that were set up using the experimental plasma concentrations) showed an apparently positive correlation with the reported lowest-observed-effect levels for haematotoxicity of these chemicals. In the case of 2,4-dimethylaniline, a methyl group at another C4-positon would be one of the determinant factors for rapid metabolic elimination to form aminotoluic acid. These results suggest that rapid and extensive metabolic activation of aniline and its dimethyl derivatives occurred in rats and that the presence of a methyl group at the C2-positon may generally suppress fast metabolic rates of dimethyl aniline derivatives that promote metabolic activation reactions at NH2 moieties.

Hepatotoxicological potential of P-toluic acid in humanised-liver mice investigated using simplified physiologically based pharmacokinetic models.

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.

Plasma and hepatic concentrations of acetaminophen and its primary conjugates after oral administrations determined in experimental animals and humans and extrapolated by pharmacokinetic modeling.

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.

Different Roles of Human Cytochrome P450 2C9 and 3A Enzymes in Diclofenac 4'- and 5-Hydroxylations Mediated by Metabolically Inactivated Human Hepatocytes in Previously Transplanted Chimeric Mice.

To investigate the respective roles of cytochromes P450 2C9 and 3A in drug oxidation in human livers, the in vivo pharmacokinetics of S-warfarin and diclofenac were analyzed after intravenous administrations in chimeric mice that had been transplanted with human hepatocytes. P450 2C9 was metabolically inactivated in the humanized mice by orally pretreating them with tienilic acid. After intravenous administration of S-warfarin, a significant difference in the concentration-time profiles of the primary metabolite 7-hydroxywarfarin between untreated mice and mice treated with tienilic acid was observed. In contrast, there were no apparent differences in the profiles for S-warfarin between the treated and untreated groups. The mean values of the maximum concentrations (Cmax) and the areas under the plasma concentration versus time curves (AUCinfinity) for 7-hydroxywarfarin were significantly lower (22 and 16% of the untreated values, respectively) in the treated group. This presumably resulted from suppressed P450 2C9 activity in the primary oxidative metabolism in vivo in the treated group. After diclofenac administration, plasma levels of diclofenac, 5-hydroxydiclofenac, and diclofenac acylglucuronide were roughly similar in pretreated and untreated mice. However, the mean Cmax and AUCinfinity values for 4'-hydroxydiclofenac were significantly lower (38 and 53% of the untreated group, respectively) in the treated group. The reported value of ∼0.8 for the fraction of S-warfarin metabolized to 7-hydroxywarfarin mediated by P450 2C9 in in vitro systems was similar to the value implied by the present humanized-liver mouse model pretreated with tienilic acid in which the AUC of 7-hydroxywarfarin was reduced by 84%. In contrast, the fractions of diclofenac metabolized to 4'-hydroxydiclofenac in in vitro and in vivo experiments were inconsistent. These results suggested that humanized-liver mice orally treated with tienilic acid might constitute an in vivo model for metabolically inactivated P450 2C9 in human hepatocytes transplanted into chimeric mice. Moreover, diclofenac, a typical in vitro P450 2C9 probe substrate, was cleared differently in vitro and in humanized-liver mice in vivo.

Different Hepatic Concentrations of Bromobenzene, 1,2-Dibromobenzene, and 1,4-Dibromobenzene in Humanized-Liver Mice Predicted Using Simplified Physiologically Based Pharmacokinetic Models as Putative Markers of Toxicological Potential.

Bromobenzene is an industrial solvent that elicits toxicity predominantly in the liver. In this study, the hepatic concentrations of bromobenzene and its related compounds 1,2-dibromobenzene and 1,4-dibromobenzene in humanized-liver mice were predicted after single oral administrations by simplified physiologically based pharmacokinetic (PBPK) models that had been set up on experimental plasma concentrations after single oral doses of 100 mg/kg to rats and 100-250 mg/kg to control mice and humanized-liver mice. The output values by simplified PBPK models were consistent with measured blood substrate concentrations in rats, control mice, and humanized-liver mice with suitable input parameter values derived from in silico prediction and the literature or estimated by fitting the measured plasma substrate concentrations. The predicted time-dependent hepatic concentrations after virtual administrations in humanized-liver mice were partly confirmed with single measured hepatic concentrations of bromobenzene and 1,4-dibromobenzene 2 h after oral doses of 150-250 mg/kg to humanized-liver mice. Moreover, leaked human albumin mRNA, a marker of the extent of human hepatic injuries, in humanized-liver mouse plasma was detected after oral administration of bromobenzene, 1,2-dibromobenzene, and 1,4-dibromobenzene. These results suggest that dosimetry approaches for determining tissue and/or blood exposures of hepatic toxicants bromobenzene, 1,2-dibromobenzene, and 1,4-dibromobenzene in humanized-liver mice were useful after virtual oral doses using simplified PBPK models. Using simplified PBPK models and plasma data from humanized-liver mice has potential to predict and evaluate the hepatic toxicity of bromobenzenes and related compounds in humanized-liver mice and in humans.

Physiologically Based Pharmacokinetic Models Predicting Renal and Hepatic Concentrations of Industrial Chemicals after Virtual Oral Doses in Rats.

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.

Steady-State Human Pharmacokinetics of Monobutyl Phthalate Predicted by Physiologically Based Pharmacokinetic Modeling Using Single-Dose Data from Humanized-Liver Mice Orally Administered with Dibutyl Phthalate.

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.

Plasma and Hepatic Concentrations of Chemicals after Virtual Oral Administrations Extrapolated Using Rat Plasma Data and Simple Physiologically Based Pharmacokinetic Models.

Only a small fraction of chemicals possesses adequate in vivo toxicokinetic data for assessing potential hazards. The aim of the present study was to model the plasma and hepatic pharmacokinetics of more than 50 disparate types of chemicals and drugs after virtual oral administrations in rats. The models were based on reported pharmacokinetics determined after oral administration to rats. An inverse relationship was observed between no-observed-effect levels after oral administration and chemical absorbance rates evaluated for cell permeability ( r = -0.98, p < 0.001, n = 17). For a varied selection of more than 30 chemicals, the plasma concentration curves and the maximum concentrations obtained using a simple one-compartment model (recently recommended as a high-throughput toxicokinetic model) and a simple physiologically based pharmacokinetic (PBPK) model (consisting of chemical receptor, metabolizing, and central compartments) were highly consistent. The hepatic and plasma concentrations and the hepatic and plasma areas under the concentration-time curves of more than 50 chemicals were roughly correlated; however, differences were evident between the PBPK-modeled values in livers and empirically obtained values in plasma. Of the compounds selected for analysis, only seven had the lowest observed effect level (LOEL) values for hepatoxicity listed in the Hazard Evaluation Support System Integrated Platform in Japan. For these seven compounds, the LOEL values and the areas under the hepatic concentration-time curves estimated using PBPK modeling were inversely correlated ( r = -0.78, p < 0.05, n = 7). This study provides important information to help simulate the high hepatic levels of potent hepatotoxic compounds. Using suitable PBPK parameters, the present models could estimate the plasma/hepatic concentrations of chemicals and drugs after oral doses using both PBPK forward and reverse dosimetry, thereby indicating the potential value of this modeling approach in predicting hepatic toxicity as a part of risk assessments of chemicals absorbed in the human body.

Human urinary concentrations of monoisononyl phthalate estimated using physiologically based pharmacokinetic modeling and experimental pharmacokinetics in humanized-liver mice orally administered with diisononyl phthalate.

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-14C]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.

R-warfarin clearances from plasma associated with polymorphic cytochrome P450 2C19 and simulated by individual physiologically based pharmacokinetic models for 11 cynomolgus monkeys.

1. Cynomolgus monkey cytochrome P450 2C19 (formerly known as P450 2C75), homologous to human P450 2C19, has been identified as R-warfarin 7-hydroxylase. In this study, simulations of R-warfarin clearance in individual cynomolgus monkeys genotyped for P450 2C19 p.[(Phe100Asn; Ala103Val; Ile112Leu)] were performed using individual simplified physiologically based pharmacokinetic (PBPK) modeling. 2. Pharmacokinetic parameters and absorption rate constants, volumes of the systemic circulation, and hepatic intrinsic clearances for individual PBPK models were estimated for eleven cynomolgus monkeys. 3. One-way ANOVA revealed significant effects of the genotype (p < 0.01) on the observed elimination half-lives and areas under the curves of R-warfarin among the homozygous mutant, heterozygous mutant, and wild-type groups. R-Warfarin clearances in individual cynomolgus monkeys genotyped for P450 2C19 were simulated by simplified PBPK modeling. The modeled hepatic intrinsic clearances were significantly associated with the P450 2C19 genotypes. The liver microsomal elimination rates of R-warfarin for individual animals after in vivo administration showed significant reductions associated with the genotype (p < 0.01). 4. This study provides important information to help simulate clearances of R-warfarin and related medicines associated with polymorphic P450 2C19 in individual cynomolgus monkeys, thereby facilitating calculation of the fraction of hepatic clearance.

Efavirenz clearances in vitro and in vivo in six cynomolgus monkeys associated with polymorphic cytochrome P450 2C9 and simulated by individual physiologically based pharmacokinetic models.

Cynomolgus monkey cytochrome P450 2C9 (formerly known as P450 2C43) variation was reportedly associated with metabolic clearance of the antiretroviral drug efavirenz in vivo (in three wild-type, one heterozygote and two homozygote animals), being unlikely in the case of human P450 2B6-dependent efavirenz clearance. In this study, the liver microsomal elimination rates of efavirenz for the same individual animals previously treated with intravenous/oral administrations of efavirenz showed significant reductions associated with the P450 2C9 p.[(I112L)] genotype (p < 0.05). Simulations of efavirenz clearance after oral administrations in individual cynomolgus monkeys were performed using individual simplified physiologically based pharmacokinetic (PBPK) modeling consisting of gut, liver and central compartments. The modeled hepatic intrinsic clearances were also significantly associated with the P450 2C9 genotypes, however, absorption rate constants or volumes of the systemic circulation were not likely determining factors for the individual efavirenz clearance variations in the six cynomolgus monkeys. This study provides important information to help simulate the clearances of efavirenz and related medicines associated with polymorphic P450 2C9 in individual cynomolgus monkeys, thereby facilitating the calculation of the fraction of liver microsomal clearance for estimating in vivo drug clearance with simplified PBPK modeling.

Forward and reverse dosimetry for aniline and 2,6-dimethylaniline in humans extrapolated from humanized-liver mouse data using simplified physiologically based pharmacokinetic models.

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.

Plasma Concentration Profiles for Hepatotoxic Pyrrolizidine Alkaloid Senkirkine in Humans Extrapolated from Rat Data Sets Using a Simplified Physiologically Based Pharmacokinetic Model.

AIM: The main aim of the current study was to obtain forward dosimetry assessments of pyrrolizidine alkaloid senkirkine plasma and liver concentrations by setting up a human physiologically based pharmacokinetic (PBPK) model based on the limited information available. BACKGROUND: The risks associated with plant-derived pyrrolizidine alkaloids as natural toxins have been assessed. OBJECTIVE: The pyrrolizidine alkaloid senkirkine was investigated because it was analyzed in a European transcriptomics study of natural hepatotoxins and in a study of the alkaloidal constituents of traditional Japanese food plants Petasites japonicus. The in silico human plasma and liver concentrations of senkirkine were modeled using doses reported for acute-term toxicity in humans. METHODS: Using a simplified PBPK model established using rat pharmacokinetic data, forward dosimetry was conducted. Since in vitro rat and human intrinsic hepatic clearances were similar; an allometric scaling approach was applied to rat parameters to create a human PBPK model. RESULTS: After oral administration of 1.0 mg/kg in rats in vivo, water-soluble senkirkine was absorbed and cleared from plasma to two orders of magnitude below the maximum concentration in 8 h. Human in silico senkirkine plasma concentration curves were generated after virtual daily oral administrations of 3.0 mg/kg senkirkine (the dose involved in an acute fatal hepatotoxicity case). A high concentration of senkirkine in the culture medium caused in vitro hepatotoxicity as evidenced by lactate dehydrogenase leakage from human hepatocyte-like HepaRG cells. CONCLUSION: Higher virtual concentrations of senkirkine in human liver and plasma than those in rat plasma were estimated using the current rat and human PBPK models. Current simulations suggest that if P. japonicus (a water-soluble pyrrolizidine alkaloid-producing plant) is ingested daily as food, hepatotoxic senkirkine could be continuously present in human plasma and liver.

Pharmacokinetics of primary oxidative metabolites of thalidomide in rats and in chimeric mice humanized with different human hepatocytes.

The approved drug thalidomide is teratogenic in humans, nonhuman primates, and rabbits but not in rodents. The extensive biotransformation of 5'-hydroxythalidomide after oral administration of thalidomide (250 mg/kg) in rats was investigated in detail using liquid chromatography-tandem mass spectrometry. Probable metabolites 5'-hydroxythalidomide sulfate and glucuronide were extensively formed, with approximately tenfold and onefold peak areas, respectively, to the primary 5'-hydroxythalidomide measured using authentic standards. As a minor metabolite, 5-hydroxythalidomide was also detected. The output of simplified physiologically based pharmacokinetic rat models was consistent with the observed in vivo data under a metabolic ratio of 0.05 for the hepatic intrinsic clearance of thalidomide to unconjugated 5'-hydroxythalidomide. The aggregate of unconjugated and sulfate/glucuronide conjugated 5'-hydroxythalidomide forms appear to be the predominant metabolites in rats. Two hours after oral administration of thalidomide (100 mg/kg) to chimeric mice humanized with four different batches of genotyped human hepatocytes, the plasma concentration ratios of 5-hydroxythalidomide to 5'-hydroxythalidomide were correlated with replacement indexes of human liver cells previously transplanted in immunodeficient mice. These results indicate that rodent livers mediate thalidomide primary oxidation, leading to extensive deactivation in vivo to unconjugated/conjugated 5'-hydroxythalidomide and suggest that thalidomide activation might be dependent on the humanized livers in mice transplanted with human hepatocytes.