Tissue-specific mRNA expression profiles of human phase I metabolizing enzymes except for cytochrome P450 and phase II metabolizing enzymes

Physiologically-based pharmacokinetic models versus allometric scaling for prediction of tyrosine-kinase inhibitor exposure from adults to children

PURPOSE

Model-based methods can predict pediatric exposure and support initial dose selection. The aim of this study was to evaluate the performance of allometric scaling of population pharmacokinetic (popPK) versus physiologically based pharmacokinetic (PBPK) models in predicting the exposure of tyrosine kinase inhibitors (TKIs) for pediatric patients (≥ 2 years), based on adult data. The drugs imatinib, sunitinib and pazopanib were selected as case studies due to their complex PK profiles including high inter-patient variability, active metabolites, time-varying clearances and non-linear absorption.

METHODS

Pediatric concentration measurements and adult popPK models were derived from the literature. Adult PBPK models were generated in PK-Sim® using available physicochemical properties, calibrated to adult data when needed. PBPK and popPK models for the pediatric populations were translated from the models for adults and were used to simulate concentration-time profiles that were compared to the observed values.

RESULTS

Ten pediatric datasets were collected from the literature. While both types of models captured the concentration-time profiles of imatinib, its active metabolite, sunitinib and pazopanib, the PBPK models underestimated sunitinib metabolite concentrations. In contrast, allometrically scaled popPK simulations accurately predicted all concentration-time profiles. Trough concentration (Ctrough) predictions from the popPK model fell within a 2-fold range for all compounds, while 3 out of 5 PBPK predictions exceeded this range for the imatinib and sunitinib metabolite concentrations.

CONCLUSION

Based on the identified case studies it appears that allometric scaling of popPK models is better suited to predict exposure of TKIs in pediatric patients ≥ 2 years. This advantage may be attributed to the stable enzyme expression patterns from 2 years old onwards, which can be easily related to adult levels through allometric scaling. In some instances, both methods performed comparably. Understanding where discrepancies between the model methods arise, can further inform model development and ultimately support pediatric dose selection.

PBPK modeling to predict the pharmacokinetics of venlafaxine and its active metabolite in different CYP2D6 genotypes and drug-drug interactions with clarithromycin and paroxetine

Venlafaxine, a serotonin-norepinephrine reuptake inhibitor (SNRI), is indicated for the treatment of major depressive disorder, social anxiety disorder, generalized anxiety disorder, and panic disorder. Venlafaxine is metabolized to the active metabolite desvenlafaxine mainly by CYP2D6. Genetic polymorphism of CYP2D6 and coadministration with other medications can significantly affect the pharmacokinetics and/or pharmacodynamics of venlafaxine and its active metabolite. This study aimed to establish the PBPK models of venlafaxine and its active metabolite related to CYP2D6 genetic polymorphism and to predict drug-drug interactions (DDIs) with clarithromycin and paroxetine in different CYP2D6 genotypes. Clinical pharmacogenomic data for venlafaxine and desvenlafaxine were collected to build the PBPK model. Physicochemical and absorption, distribution, metabolism, and excretion (ADME) characteristics of respective compounds were obtained from previously reported data, predicted by the PK-Sim® software, or optimized to capture the plasma concentration-time profiles. Model evaluation was performed by comparing the predicted pharmacokinetic parameters and plasma concentration-time profiles to the observed data. Predicted plasma concentration-time profiles of venlafaxine and its active metabolite were visually similar to the observed profiles and all predicted AUC and Cmax values for respective compounds were included in the twofold error range of observed values in non-genotyped populations and different CYP2D6 genotypes. When clarithromycin or clarithromycin plus paroxetine was concomitantly administered, predicted plasma concentration-time profiles of venlafaxine properly captured the observed profiles in two different CYP2D6 genotypes and all predicted DDI ratios for AUC and Cmax were included within the acceptance range. Consequently, the present model successfully captured the pharmacokinetic alterations of venlafaxine and its active metabolite according to CYP2D6 genetic polymorphism as well as the DDIs between venlafaxine and two CYP inhibitors. The present model can be used to predict the pharmacokinetics of venlafaxine and its active metabolite considering different races, ages, coadministered drugs, and CYP2D6 activity of individuals and it can contribute to individualized pharmacotherapy of venlafaxine.

Maximum likelihood estimation of renal transporter ontogeny profiles for pediatric PBPK modeling

Optimal treatment of infants with many renally cleared drugs must account for maturational differences in renal transporter (RT) activity. Pediatric physiologically-based pharmacokinetic (PBPK) models may incorporate RT activity, but this requires ontogeny profiles for RT activity in children, especially neonates, to predict drug disposition. Therefore, RT expression measurements from human kidney postmortem cortical tissue samples were normalized to represent a fraction of mature RT activity. Using these data, maximum likelihood estimated the distributions of RT activity across the pediatric age spectrum, including preterm and term neonates. PBPK models of four RT substrates (acyclovir, ciprofloxacin, furosemide, and meropenem) were evaluated with and without ontogeny profiles using average fold error (AFE), absolute average fold error (AAFE), and proportion of observations within the 5-95% prediction interval. Novel maximum likelihood profiles estimated ontogeny distributions for the following RT: OAT1, OAT3, OCT2, P-gp, URAT1, BCRP, MATE1, MRP2, MRP4, and MATE-2 K. Profiles for OAT3, P-gp, and MATE1 improved infant furosemide and neonate meropenem PBPK model AFE from 0.08 to 0.70 and 0.53 to 1.34 and model AAFE from 12.08 to 1.44 and 2.09 to 1.36, respectively, and improved the percent of data within the 5-95% prediction interval from 48% to 98% for neonatal ciprofloxacin simulations, respectively. Even after accounting for other critical population-specific maturational differences, novel RT ontogeny profiles substantially improved neonatal PBPK model performance, providing validated estimates of maturational differences in RT activity for optimal dosing in children.

Maternal Ezetimibe Concentrations Measured in Breast Milk and Its Use in Breastfeeding Infant Exposure Predictions

BACKGROUND

Lactating mothers taking ezetimibe, an antihyperlipidemic agent, may be hesitant to breastfeed despite the known benefit of breastfeeding to both mother and infant. Currently, no data exist on the presence or concentration of ezetimibe and its main active metabolite, ezetimibe-glucuronide (EZE-glucuronide), in human breast milk.

METHODS

Voluntary breast milk samples containing ezetimibe and EZE-glucuronide were attained from lactating mothers taking ezetimibe as part of their treatment. An assay was developed and validated to measure ezetimibe and EZE-glucuronide concentrations in breast milk. A workflow that utilized a developed and evaluated pediatric physiologically based pharmacokinetic (PBPK) model, the measured concentrations in milk, and weight-normalized breast milk intake volumes was applied to predict infant exposures and determine the upper area under the curve ratio (UAR).

RESULTS

Fifteen breast milk samples from two maternal-infant pairs were collected. The developed liquid chromatography-tandem mass spectrometry (LC-MS/MS) assay showed an analytical range of 0.039-5.0 ng/mL and 0.39-50.0 ng/mL for ezetimibe and EZE-glucuronide, respectively. The measured concentrations in the breast milk samples were 0.17-1.02 ng/mL and 0.42-2.65 ng/mL of ezetimibe and EZE-glucuronide, respectively. The evaluated pediatric PBPK model demonstrated minimal exposure overlap in adult therapeutic dose and breastfed infant simulated area under the concentration-time curve from time zero to 24 h (AUC24). Calculated UAR across infant age groups ranged from 0.0015 to 0.0026.

CONCLUSIONS

PBPK model-predicted ezetimibe and EZE-glucuronide exposures and UAR suggest that breastfeeding infants would receive non-therapeutic exposures. Future work should involve a 'mother-infant pair study' to ascertain breastfed infant plasma ezetimibe and EZE-glucuronide concentrations to confirm the findings of this work.

Physiologically Based Pharmacokinetic (PBPK) Modeling to Predict CYP3A-Mediated Drug Interaction between Saxagliptin and Nicardipine: Bridging Rat-to-Human Extrapolation

The aim of this study was to predict the cytochrome P450 3A (CYP3A)-mediated drug-drug interactions (DDIs) between saxagliptin and nicardipine using a physiologically based pharmacokinetic (PBPK) model. Initially, in silico and in vitro parameters were gathered from experiments or the literature to construct PBPK models for each drug in rats. These models were integrated to predict the DDIs between saxagliptin, metabolized via CYP3A2, and nicardipine, exhibiting CYP3A inhibitory activity. The rat DDI PBPK model was completed by optimizing parameters using experimental rat plasma concentrations after co-administration of both drugs. Following co-administration in Sprague-Dawley rats, saxagliptin plasma concentration significantly increased, resulting in a 2.60-fold rise in AUC, accurately predicted by the rat PBPK model. Subsequently, the workflow of the rat PBPK model was applied to humans, creating a model capable of predicting DDIs between the two drugs in humans. Simulation from the human PBPK model indicated that nicardipine co-administration in humans resulted in a nearly unchanged AUC of saxagliptin, with an approximate 1.05-fold change, indicating no clinically significant changes and revealing a lack of direct translation of animal interaction results to humans. The animal-to-human PBPK model extrapolation used in this study could enhance the reliability of predicting drug interactions in clinical settings where DDI studies are challenging.

Physiologically based pharmacokinetic (PBPK) modeling of pitavastatin in relation to SLCO1B1 genetic polymorphism

Pitavastatin, a potent 3-hydroxymethylglutaryl coenzyme A reductase inhibitor, is indicated for the treatment of hypercholesterolemia and mixed dyslipidemia. Hepatic uptake of pitavastatin is predominantly occupied by the organic anion transporting polypeptide 1B1 (OATP1B1) and solute carrier organic anion transporter family member 1B1 (SLCO1B1) gene, which is a polymorphic gene that encodes OATP1B1. SLCO1B1 genetic polymorphism significantly alters the pharmacokinetics of pitavastatin. This study aimed to establish the physiologically based pharmacokinetic (PBPK) model to predict pitavastatin pharmacokinetics according to SLCO1B1 genetic polymorphism. PK-Sim® version 10.0 was used to establish the whole-body PBPK model of pitavastatin. Our pharmacogenomic data and a total of 27 clinical pharmacokinetic data with different dose administration and demographic properties were used to develop and validate the model, respectively. Physicochemical properties and disposition characteristics of pitavastatin were acquired from previously reported data or optimized to capture the plasma concentration-time profiles in different SLCO1B1 diplotypes. Model evaluation was performed by comparing the predicted pharmacokinetic parameters and profiles to the observed data. Predicted plasma concentration-time profiles were visually similar to the observed profiles in the non-genotyped populations and different SLCO1B1 diplotypes. All fold error values for AUC and Cmax were included in the two fold range of observed values. Thus, the PBPK model of pitavastatin in different SLCO1B1 diplotypes was properly established. The present study can be useful to individualize the dose administration strategy of pitavastatin in individuals with various ages, races, and SLCO1B1 diplotypes.

Applied physiologically-based pharmacokinetic modeling to assess uridine diphosphate-glucuronosyltransferase-mediated drug-drug interactions for Vericiguat

Vericiguat (Verquvo; US: Merck, other countries: Bayer) is a novel drug for the treatment of chronic heart failure. Preclinical studies have demonstrated that the primary route of metabolism for vericiguat is glucuronidation, mainly catalyzed by uridine diphosphate-glucuronosyltransferase (UGT)1A9 and to a lesser extent UGT1A1. Whereas a drug-drug interaction (DDI) study of the UGT1A9 inhibitor mefenamic acid showed a 20% exposure increase, the effect of UGT1A1 inhibitors has not been assessed clinically. This modeling study describes a physiologically-based pharmacokinetic (PBPK) approach to complement the clinical DDI liability assessment and support prescription labeling. A PBPK model of vericiguat was developed based on in vitro and clinical data, verified against data from the mefenamic acid DDI study, and applied to assess the UGT1A1 DDI liability by running an in silico DDI study with the UGT1A1 inhibitor atazanavir. A minor effect with an area under the plasma concentration-time curve (AUC) ratio of 1.12 and a peak plasma concentration ratio of 1.04 was predicted, which indicates that there is no clinically relevant DDI interaction anticipated. Additionally, the effect of potential genetic polymorphisms of UGT1A1 and UGT1A9 was evaluated, which showed that an average modest increase of up to 1.7-fold in AUC may be expected in the case of concomitantly reduced UGT1A1 and UGT1A9 activity for subpopulations expressing non-wild-type variants for both isoforms. This study is a first cornerstone to qualify the PK-Sim platform for use of UGT-mediated DDI predictions, including PBPK models of perpetrators, such as mefenamic acid and atazanavir, and sensitive UGT substrates, such as dapagliflozin and raltegravir.

Physiologically based pharmacokinetic (PBPK) modeling to predict the pharmacokinetics of irbesartan in different CYP2C9 genotypes

Irbesartan, a potent and selective angiotensin II type-1 (AT1) receptor blocker (ARB), is one of the representative medications for the treatment of hypertension. Cytochrome P450 (CYP) 2C9 is primarily involved in the oxidation of irbesartan. CYP2C9 is highly polymorphic, and genetic polymorphism of this enzyme is the leading cause of significant alterations in the pharmacokinetics of irbesartan. This study aimed to establish the physiologically based pharmacokinetic (PBPK) model to predict the pharmacokinetics of irbesartan in different CYP2C9 genotypes. The irbesartan PBPK model was established using the PK-Sim® software. Our previously reported pharmacogenomic data for irbesartan was leveraged in the development of the PBPK model and collected clinical pharmacokinetic data for irbesartan was used for the validation of the model. Physicochemical and ADME properties of irbesartan were obtained from previously reported data, predicted by the modeling software, or optimized to fit the observed plasma concentration-time profiles. Model evaluation was performed by comparing the predicted plasma concentration-time profiles and pharmacokinetic parameters to the observed results. Predicted plasma concentration-time profiles were visually similar to observed profiles. Predicted AUCinf in CYP2C9*1/*3 and CYP2C9*1/*13 genotypes were increased by 1.54- and 1.62-fold compared to CYP2C9*1/*1 genotype, respectively. All fold error values for AUC and Cmax in non-genotyped and CYP2C9 genotyped models were within the two-fold error criterion. We properly established the PBPK model of irbesartan in different CYP2C9 genotypes. It can be used to predict the pharmacokinetics of irbesartan for personalized pharmacotherapy in individuals of various races, ages, and CYP2C9 genotypes.

Aldehyde oxidase 1 activity and protein expression in human, rabbit, and pig ocular tissues

Aldehyde oxidase (AOX) is a cytosolic drug-metabolizing enzyme which has attracted increasing attention in drug development due to its high hepatic expression, broad substrate profile and species differences. In contrast, there is limited information on the presence and activity of AOX in extrahepatic tissues including ocular tissues. Because several ocular drugs are potential substrates for AOX, we performed a comprehensive analysis of the AOX1 expression and activity profile in seven ocular tissues from humans, rabbits, and pigs. AOX activities were determined using optimized assays for the established human AOX1 probe substrates 4-dimethylamino-cinnamaldehyde (DMAC) and phthalazine. Inhibition studies were undertaken in conjunctival and retinal homogenates using well-established human AOX1 inhibitors menadione and chlorpromazine. AOX1 protein contents were quantitated with targeted proteomics and confirmed by immunoblotting. Overall, DMAC oxidation rates varied over 10-fold between species (human ˃˃ rabbit ˃ pig) and showed 2- to 6-fold differences between tissues from the same species. Menadione seemed a more potent inhibitor of DMAC oxidation across species than chlorpromazine. Human AOX1 protein levels were highest in the conjunctiva, followed by most posterior tissues, whereas anterior tissues showed low levels. The rabbit AOX1 expression was high in the conjunctiva, retinal pigment epithelial (RPE), and choroid while lower in the anterior tissues. Quantification of pig AOX1 was not successful but immunoblotting confirmed the presence of AOX1 in all species. DMAC oxidation rates and AOX1 contents correlated quite well in humans and rabbits. This study provides, for the first time, insights into the ocular expression and activity of AOX1 among multiple species.

Challenges and Opportunities for In Vitro-In Vivo Extrapolation of Aldehyde Oxidase-Mediated Clearance: Toward a Roadmap for Quantitative Translation

Underestimation of aldehyde oxidase (AO)-mediated clearance by current in vitro assays leads to uncertainty in human dose projections, thereby reducing the likelihood of success in drug development. In the present study we first evaluated the current drug development practices for AO substrates. Next, the overall predictive performance of in vitro-in vivo extrapolation of unbound hepatic intrinsic clearance (CLint,u) and unbound hepatic intrinsic clearance by AO (CLint,u,AO) was assessed using a comprehensive literature database of in vitro (human cytosol/S9/hepatocytes) and in vivo (intravenous/oral) data collated for 22 AO substrates (total of 100 datapoints from multiple studies). Correction for unbound fraction in the incubation was done by experimental data or in silico predictions. The fraction metabolized by AO (fmAO) determined via in vitro/in vivo approaches was found to be highly variable. The geometric mean fold errors (gmfe) for scaled CLint,u (mL/min/kg) were 10.4 for human hepatocytes, 5.6 for human liver cytosols, and 5.0 for human liver S9, respectively. Application of these gmfe's as empirical scaling factors improved predictions (45%-57% within twofold of observed) compared with no correction (11%-27% within twofold), with the scaling factors qualified by leave-one-out cross-validation. A road map for quantitative translation was then proposed following a critical evaluation on the in vitro and clinical methodology to estimate in vivo fmAO In conclusion, the study provides the most robust system-specific empirical scaling factors to date as a pragmatic approach for the prediction of in vivo CLint,u,AO in the early stages of drug development. SIGNIFICANCE STATEMENT: Confidence remains low when predicting in vivo clearance of AO substrates using in vitro systems, leading to de-prioritization of AO substrates from the drug development pipeline to mitigate risk of unexpected and costly in vivo impact. The current study establishes a set of empirical scaling factors as a pragmatic tool to improve predictability of in vivo AO clearance. Developing clinical pharmacology strategies for AO substrates by utilizing mass balance/clinical drug-drug interaction data will help build confidence in fmAO.

Verifying in vitro-determined enzyme contributions to cannabidiol clearance for exposure predictions in human through physiologically-based pharmacokinetic modeling

Cannabidiol (CBD) is approved for treatment of seizures associated with two forms of epilepsy that become apparent in infancy or early childhood. To consider an adult physiologically-based pharmacokinetic (PBPK) model for pediatric scaling, we assessed in vitro-derived cytochrome P450 (CYP) and uridine 5'-diphospho-glucuronosyltransferase (UGT) enzyme contributions to CBD clearance in human. An i.v. PBPK model was constructed using CBD physicochemical properties and knowledge of disposition. The i.v. datasets were used for model building and evaluation. Oral PBPK models for CBD administered in fasted and fed states were developed using single dose oral datasets and parameters optimized from the i.v. model and evaluated with multiple dose datasets. Relative contributions of CBD metabolizing enzymes were partitioned according to in vitro studies. Clinical drug-drug interaction (DDI) studies were simulated using CBD fed state, itraconazole, fluconazole, and rifampicin PBPK models. Linear mixed effect modeling was used to estimate area under the concentration-time curve from zero to infinity (AUC0-∞ ) perpetrator + CBD versus CBD alone. The i.v. and oral datasets used in model evaluation produced acceptable average fold error (AFE) of 1.28 and absolute AFE of 1.65. Relative contributions of drug-metabolizing enzymes to CBD clearance were proposed from in vitro data: UGT1A7 4%, UGT1A9 16%, UGT2B7 10%, CYP3A4 38%, CYP2C19 21%, and CYP2C9 11%. The simulated DDI studies using the in vitro-derived values produced AUC0-∞ treatment ratios comparable to observed: itraconazole 1.24 versus 1.07, fluconazole 1.45 versus 1.22, and rifampicin 0.49 versus 0.69. The constructed CBD PBPK models can predict adult exposures and have potential for use in pediatrics where exposure estimates are limited.

Pharmacokinetics, mass balance, and metabolism of [14C]TPN171, a novel PDE5 inhibitor, in humans for the treatment of pulmonary arterial hypertension

TPN171 is a novel phosphodiesterase-5 (PDE5) inhibitor used to treat pulmonary arterial hypertension (PAH) and erectile dysfunction (ED), which currently is undergoing phase II clinical trials in China. In this single-center, single-dose, nonrandomized, and open design study, radiolabeled [¹⁴C]TPN171 was used to investigate the metabolic mechanism, pharmacokinetic characteristics, and clearance pathways of TPN171 in 6 healthy Chinese male volunteers. Each volunteer was administered a single oral suspension of 10 mg (100 μCi) of [¹⁴C]TPN171. We found that TPN171 was absorbed rapidly in humans with a peak time (T_max) of 0.667 h and a half-life (t_1/2) of approximately 9.89 h in plasma. Excretion of radiopharmaceutical-related components was collected 216 h after administration, accounting for 95.21% of the dose (46.61% in urine and 48.60% in feces). TPN171 underwent extensive metabolism in humans. Twenty-two metabolites were detected in human plasma, urine, and feces using a radioactive detector combined with a high-resolution mass spectrometer. According to radiochromatograms, a glucuronide metabolite of O-dealkylated TPN171 exceeded 10% of the total drug-related components in human plasma. However, according to the Food and Drug Administration (FDA) guidelines, no further tests are needed to evaluate the safety of this metabolite because it is a phase II metabolite, but the compound is still worthy of attention. The main metabolic biotransformation of TPN171 was mono-oxidation (hydroxylation and N-oxidation), dehydrogenation, N-dealkylation, O-dealkylation, amide hydrolysis, glucuronidation, and acetylation. Cytochrome P450 3A4 (CYP3A4) mainly catalyzed the formation of metabolites, and CYP2E1 and CYP2D6 were involved in the oxidative metabolism of TPN171 to a lesser extent. According to the incubation data, M1 was mainly metabolized to M1G by UDP-glucuronosyltransferase 1A9 (UGT1A9), followed by UGT1A7 and UGT1A10.

The Somatic Mutation Landscape of UDP-Glycosyltransferase (UGT) Genes in Human Cancers

The human UDP-glycosyltransferase (UGTs) superfamily has a critical role in the metabolism of anticancer drugs and numerous pro/anti-cancer molecules (e.g., steroids, lipids, fatty acids, bile acids and carcinogens). Recent studies have shown wide and abundant expression of UGT genes in human cancers. However, the extent to which UGT genes acquire somatic mutations within tumors remains to be systematically investigated. In the present study, our comprehensive analysis of the somatic mutation profiles of 10,069 tumors from 33 different TCGA cancer types identified 3427 somatic mutations in UGT genes. Overall, nearly 18% (1802/10,069) of the assessed tumors had mutations in UGT genes with huge variations in mutation frequency across different cancer types, ranging from over 25% in five cancers (COAD, LUAD, LUSC, SKCM and UCSC) to less than 5% in eight cancers (LAML, MESO, PCPG, PAAD, PRAD, TGCT, THYM and UVM). All 22 UGT genes showed somatic mutations in tumors, with UGT2B4, UGT3A1 and UGT3A2 showing the largest number of mutations (289, 307 and 255 mutations, respectively). Nearly 65% (2260/3427) of the mutations were missense, frame-shift and nonsense mutations that have been predicted to code for variant UGT proteins. Furthermore, about 10% (362/3427) of the mutations occurred in non-coding regions (5' UTR, 3' UTR and splice sites) that may be able to alter the efficiency of translation initiation, miRNA regulation or the splicing of UGT transcripts. In conclusion, our data show widespread somatic mutations of UGT genes in human cancers that may affect the capacity of cancer cells to metabolize anticancer drugs and endobiotics that control pro/anti-cancer signaling pathways. This highlights their potential utility as biomarkers for predicting therapeutic efficacy and clinical outcomes.

Sulfation of 12-hydroxy-nevirapine by human SULTs and the effects of genetic polymorphisms of SULT1A1 and SULT2A1

Nevirapine (NVP) is an effective drug for the treatment of HIV infections, but its use is limited by a high incidence of severe skin rash and liver injury. 12-Hydroxynevirapine (12-OH-NVP) is the major metabolite of nevirapine. There is strong evidence that the sulfate of 12-OH-NVP is responsible for the skin rash. While several cytosolic sulfotransferases (SULTs) have been shown to be capable of sulfating 12-OH-NVP, the exact mechanism of sulfation in vivo is unclear. The current study aimed to clarify human SULT(s) and human organs that are capable of sulfating 12-OH-NVP and investigate the metabolic sulfation of 12-OH-NVP using cultured HepG2 human hepatoma cells. Enzymatic assays revealed that of the thirteen human SULTs, SULT1A1 and SULT2A1 displayed strong 12-OH-NVP-sulfating activity. 1-Phenyl-1-hexanol (PHHX), which applied topically prevents the skin rash in rats, inhibited 12-OH-NVP sulfation by SULT1A1 and SULT2A1, implying the involvement of these two enzymes in the sulfation of 12-OH-NVP in vivo. Among five human organ cytosols analyzed, liver cytosol displayed the strongest 12-OH-NVP-sulfating activity, while a low but significant activity was detected with skin cytosol. Cultured HepG2 cells were shown to be capable of sulfating 12-OH-NVP. The effects of genetic polymorphisms of SULT1A1 and SULT2A1 genes on the sulfation of 12-OH-NVP by SULT1A1 and SULT2A1 allozymes were investigated. Two SULT1A1 allozymes, Arg37Asp and Met223Val, showed no detectable 12-OH-NVP-sulfating activity, while a SULT2A1 allozyme, Met57Thr, displayed significantly higher 12-OH-NVP-sulfating activity compared with the wild-type enzyme. Collectively, these results contribute to a better understanding of the involvement of sulfation in NVP-induced skin rash and provide clues to the possible role of SULT genetic polymorphisms in the risk of this adverse reaction.

Physiologically based pharmacokinetic (PBPK) modeling of flurbiprofen in different CYP2C9 genotypes

The aim of this study was to establish the physiologically based pharmacokinetic (PBPK) model of flurbiprofen related to CYP2C9 genetic polymorphism and describe the pharmacokinetics of flurbiprofen in different CYP2C9 genotypes. PK-Sim® software was used for the model development and validation. A total of 16 clinical pharmacokinetic data for flurbiprofen in different CYP2C9 genotypes, dose regimens, and age groups were used for the PBPK modeling. Turnover number (k_cat) of CYP2C9 values were optimized to capture the observed profiles in different CYP2C9 genotypes. In the simulation, predicted fraction metabolized by CYP2C9, fraction excreted to urine, bioavailability, and volume of distribution were similar to previously reported values. Predicted plasma concentration-time profiles in different CYP2C9 genotypes were visually similar to the observed profiles. Predicted AUC_inf in CYP2C9*1/*2, CYP2C9*1/*3, and CYP2C9*3/*3 genotypes were 1.44-, 2.05-, and 3.67-fold higher than the CYP2C9*1/*1 genotype. The ranges of fold errors for AUC_inf, C_max, and t_1/2 were 0.84-1.00, 0.61-1.22, and 0.74-0.94 in development and 0.59-0.98, 0.52-0.97, and 0.61-1.52 in validation, respectively, which were within the acceptance criterion. Thus, the PBPK model was successfully established and described the pharmacokinetics of flurbiprofen in different CYP2C9 genotypes, dose regimens, and age groups. The present model could guide the decision-making of tailored drug administration strategy by predicting the pharmacokinetics of flurbiprofen in various clinical scenarios.

Physiologically based pharmacokinetic (PBPK) modeling of piroxicam with regard to CYP2C9 genetic polymorphism

Piroxicam is a non-steroidal anti-inflammatory drug used to alleviate symptoms of osteoarthritis and rheumatoid arthritis. CYP2C9 genetic polymorphism significantly influences the pharmacokinetics of piroxicam. The objective of this study was to develop and validate the piroxicam physiologically based pharmacokinetic (PBPK) model related to CYP2C9 genetic polymorphism. PK-Sim^® version 10.0 was used for the PBPK modeling. The PBPK model was evaluated by predicted and observed plasma concentration-time profiles, fold errors of predicted to observed pharmacokinetic parameters, and a goodness-of-fit plot. The turnover number (k_cat) of CYP2C9 was adjusted to capture the pharmacokinetics of piroxicam in different CYP2C9 genotypes. The population PBPK model overall accurately described and predicted the plasma concentration-time profiles in different CYP2C9 genotypes. In our simulations, predicted AUC_inf in CYP2C9*1/*2, CYP2C9*1/*3, and CYP2C9*3/*3 genotypes were 1.83-, 2.07-, and 6.43-fold higher than CYP2C9*1/*1 genotype, respectively. All fold error values for AUC, C_max, and t_1/2 were included in the acceptance criterion with the ranges of 0.57-1.59, 0.63-1.39, and 0.65-1.51, respectively. The range of fold error values for predicted versus observed plasma concentrations was 0.11-3.13. 93.9% of fold error values were within the two-fold range. Average fold error, absolute average fold error, and root mean square error were 0.93, 1.27, and 0.72, respectively. Our model accurately captured the pharmacokinetic alterations of piroxicam according to CYP2C9 genetic polymorphism.

Odorant Metabolism in Humans

Odorants are relatively small molecules which are easily taken up and distributed in the human body. Despite their relevance in everyday life, however, only a limited amount of evidence about their metabolism, pathways, and bioactivities in the human body exists. With this Review, we aim to encourage future interdisciplinary research on the function and mechanisms of the biotransformation of odorants, involving different disciplines such as nutrition, medicine, biochemistry, chemistry, and sensory sciences. Starting with a general overview of the different ways of odorant uptake and enzymes involved in the metabolism of odorants, a more precise description of biotransformation processes and their function in the oral cavity, the nose, the lower respiratory tract (LRT), and the gastrointestinal tract (GIT) is given together with an overview of the different routes of odorant excretion. Finally, perspectives for future research are discussed.

Physiologically Based Pharmacokinetic (PBPK) Modeling of Clopidogrel and Its Four Relevant Metabolites for CYP2B6, CYP2C8, CYP2C19, and CYP3A4 Drug–Drug–Gene Interaction Predictions

The antiplatelet agent clopidogrel is listed by the FDA as a strong clinical index inhibitor of cytochrome P450 (CYP) 2C8 and weak clinical inhibitor of CYP2B6. Moreover, clopidogrel is a substrate of—among others—CYP2C19 and CYP3A4. This work presents the development of a whole-body physiologically based pharmacokinetic (PBPK) model of clopidogrel including the relevant metabolites, clopidogrel carboxylic acid, clopidogrel acyl glucuronide, 2-oxo-clopidogrel, and the active thiol metabolite, with subsequent application for drug–gene interaction (DGI) and drug–drug interaction (DDI) predictions. Model building was performed in PK-Sim^® using 66 plasma concentration-time profiles of clopidogrel and its metabolites. The comprehensive parent-metabolite model covers biotransformation via carboxylesterase (CES) 1, CES2, CYP2C19, CYP3A4, and uridine 5′-diphospho-glucuronosyltransferase 2B7. Moreover, CYP2C19 was incorporated for normal, intermediate, and poor metabolizer phenotypes. Good predictive performance of the model was demonstrated for the DGI involving CYP2C19, with 17/19 predicted DGI AUC_last and 19/19 predicted DGI C_max ratios within 2-fold of their observed values. Furthermore, DDIs involving bupropion, omeprazole, montelukast, pioglitazone, repaglinide, and rifampicin showed 13/13 predicted DDI AUC_last and 13/13 predicted DDI C_max ratios within 2-fold of their observed ratios. After publication, the model will be made publicly accessible in the Open Systems Pharmacology repository.

Investigation of Janus Kinase (JAK) Inhibitors for Lung Delivery and the Importance of Aldehyde Oxidase Metabolism

The Janus family of tyrosine kinases (JAK1, JAK2, JAK3, and TYK2) play an essential role in the receptor signaling of cytokines that have been implicated in the pathogenesis of severe asthma, and there is emerging interest in the development of small-molecule-inhaled JAK inhibitors as treatments. Here, we describe the optimization of a quinazoline series of JAK inhibitors and the results of mouse lung pharmacokinetic (PK) studies where only low concentrations of parent compound were observed. Subsequent investigations revealed that the low exposure was due to metabolism by aldehyde oxidase (AO), so we sought to identify quinazolines that were not metabolized by AO. We found that specific substituents at the quinazoline 2-position prevented AO metabolism and this was rationalized through computational docking studies in the AO binding site, but they compromised kinome selectivity. Results presented here highlight that AO metabolism is a potential issue in the lung.

Physiologically based pharmacokinetic (PBPK) modeling of meloxicam in different CYP2C9 genotypes

Meloxicam, a non-steroidal anti-inflammatory drug, is used for the treatment of rheumatoid arthritis and osteoarthritis. Cytochrome P450 (CYP) 2C9 and CYP3A4 are major and minor enzymes involved in the metabolism of meloxicam. Impaired enzyme activity of CYP2C9 variants increases the plasma exposures of meloxicam and the risk of adverse events. The objective of our study is to develop and validate the physiologically based pharmacokinetic (PBPK) model of meloxicam related to CYP2C9 genetic polymorphism using the PK-Sim^® software. In vitro k_cat of CYP2C9 was optimized in different CYP2C9 genotypes. The demographic and pharmacokinetic dataset for the development of the PBPK model was extracted from two previous clinical pharmacokinetic studies. Thirty-one clinical datasets, representing different dose regimens and demographic characteristics, were utilized to validate the PBPK model. The shapes of simulated plasma concentration-time profiles in each CYP2C9 genotype were visually similar to observed profiles. The predicted exposures (AUC_inf) of meloxicam in CYP2C9*1/*3, CYP2C9*1/*13, and CYP2C9*3/*3 genotypes were increased by 1.77-, 2.91-, and 8.35-fold compared to CYP2C9*1/*1 genotype, respectively. In all datasets for the development and validations, fold errors between predicted and observed pharmacokinetic parameters were within the two-fold error criteria. As a result, the PBPK model was appropriately established and properly described the pharmacokinetics of meloxicam in different CYP2C9 genotypes. This study is expected to contribute to reducing the risk of adverse events of meloxicam through optimization of meloxicam dosing in different CYP2C9 genotypes.

Whole Body PBPK Modeling of Remdesivir and its Metabolites to Aid in Estimating Active Metabolite Exposure in the Lung and Liver in Patients with Organ Dysfunction

Remdesivir (RDV) is the first drug approved by the US Food and Drug Administration (FDA) for the treatment of coronavirus disease 2019 (COVID-19) in certain patients requiring hospitalization. As a nucleoside analogue prodrug, RDV undergoes intracellular multistep activation to form its pharmacologically active species, GS-443902, which is not detectable in the plasma. A question arises that whether the observed plasma exposure of RDV and its metabolites would correlate with or be informative about the exposure of GS-443902 in tissues. A whole body physiologically-based pharmacokinetic (PBPK) modeling and simulation approach was utilized to elucidate the disposition mechanism of RDV and its metabolites in the lungs and liver and explore the relationship between plasma and tissue pharmacokinetics (PK) of RDV and its metabolites in healthy subjects. In addition, the potential alteration of plasma and tissue PK of RDV and its metabolites in patients with organ dysfunction was explored. Our simulation results indicated that intracellular exposure of GS-443902 was decreased in the liver and increased in the lungs in subjects with hepatic impairment relative to the subjects with normal liver function. In subjects with severe renal impairment, the exposure of GS-443902 in the liver was slightly increased, whereas the lung exposure of GS-443902 was not impacted. These predictions along with the organ impairment study results may be used to support decision making regarding the RDV dosage adjustment in these patient subgroups. The modeling exercise illustrated the potential of whole body PBPK modeling to aid in decision making for nucleotide analogue prodrugs, particularly when the active metabolite exposure in the target tissues is not available.

UGT1A Gene Family Members Serve as Potential Targets and Prognostic Biomarkers for Pancreatic Cancer

Background

Pancreatic cancer (PC) is one of the most common cancers worldwide, with high mortality. The UGT1A gene family plays important roles in pharmacology and toxicology, contributing to interindividual differences in drug disposition. However, mRNA expression and prognostic value of the UGT1A gene family in PC have not been identified.

Methods

Oncomine, GEPIA2, DAVID 6.8, Metascape, Kaplan-Meier plotter, cBioPortal, GeneMANIA, TRRUST v2, TIMER, and R software were used in our study.

Results

The transcriptional levels of UGT1A1/3/6/8/9/10 in PC tissues were significantly higher than those in normal tissues. These results were further validated using five pairs of PC tumor tissues and adjacent nontumor tissues. A significant correlation was found between the expression of UGT1A1/6/10 and the pathological stage of PC. PC patients with lower transcriptional levels of UGT1A1/4/5/6/10 were associated with a better prognosis. The differentially expressed UGT1A gene family functions were primarily related to the glucuronidation pathway, cytokine-cytokine receptor interactions, and the ILK signaling pathway. Our data suggest that HNF1A, AHR, and CDX2 are key transcription factors for the UGT1A gene family. Furthermore, the expression levels of UGT1A1/3/8/9/10 were positively correlated with the activities of tumor-infiltrating immune cells, especially B cells. The expression levels of UGT1A6/9 were negatively correlated with macrophage infiltration levels.

Conclusions

These results suggest that the UGT1A gene family could serve as a potential prognostic biomarker and target for PC. However, future studies are required to validate our findings and promote the clinical utility of the UGT1A gene family in PC.

Non-P450 Drug-Metabolizing Enzymes: Contribution to Drug Disposition, Toxicity, and Development

Most clinically used drugs are metabolized in the body via oxidation, reduction, or hydrolysis reactions, which are considered phase I reactions. Cytochrome P450 (P450) enzymes, which primarily catalyze oxidation reactions, contribute to the metabolism of over 50% of clinically used drugs. In the last few decades, the function and regulation of P450s have been extensively studied, whereas the characterization of non-P450 phase I enzymes is still incomplete. Recent studies suggest that approximately 30% of drug metabolism is carried out by non-P450 enzymes. This review summarizes current knowledge of non-P450 phase I enzymes, focusing on their roles in controlling drug efficacy and adverse reactions as an important aspect of drug development. Expected final online publication date for the Annual Review of Pharmacology and Toxicology, Volume 62 is January 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

The Expression Profiles and Deregulation of UDP-Glycosyltransferase (UGT) Genes in Human Cancers and Their Association with Clinical Outcomes

The human UDP-glycosyltransferase (UGTs) superfamily has 22 functional enzymes that play a critical role in the metabolism of small lipophilic compounds, including carcinogens, drugs, steroids, lipids, fatty acids, and bile acids. The expression profiles of UGT genes in human cancers and their impact on cancer patient survival remains to be systematically investigated. In the present study, a comprehensive analysis of the RNAseq and clinical datasets of 9514 patients from 33 different TCGA (the Genome Cancer Atlas) cancers demonstrated cancer-specific UGT expression profiles with high interindividual variability among and within individual cancers. Notably, cancers derived from drug metabolizing tissues (liver, kidney, gut, pancreas) expressed the largest number of UGT genes (COAD, KIRC, KIRP, LIHC, PAAD); six UGT genes (1A6, 1A9, 1A10, 2A3, 2B7, UGT8) showed high expression in five or more different cancers. Kaplan-Meier plots and logrank tests revealed that six UGT genes were significantly associated with increased overall survival (OS) rates [UGT1A1 (LUSC), UGT1A6 (ACC), UGT1A7 (ACC), UGT2A3 (KIRC), UGT2B15 (BLCA, SKCM)] or decreased OS rates [UGT2B15 (LGG), UGT8 (UVM)] in specific cancers. Finally, differential expression analysis of 611 patients from 12 TCGA cancers identified 16 UGT genes (1A1, 1A3, 1A6, 1A7, 1A8, 1A9, 1A10, 2A1, 2A3, 2B4, 2B7, 2B11, 2B15, 3A1, 3A2, UGT8) that were up/downregulated in at least one cancer relative to normal tissues. In conclusion, our data show widespread expression of UGT genes in cancers, highlighting the capacity for intratumoural drug metabolism through the UGT conjugation pathway. The data also suggests the potentials for specific UGT genes to serve as prognostic biomarkers or therapeutic targets in cancers.

Physiologically based pharmacokinetic (PBPK) modeling for prediction of celecoxib pharmacokinetics according to CYP2C9 genetic polymorphism

Celecoxib is a non-steroidal anti-inflammatory drug (NSAID) and a representative selective cyclooxygenase (COX)-2 inhibitor, which is commonly prescribed for osteoarthritis, rheumatoid arthritis, ankylosing spondylitis, acute pain, and primary dysmenorrhea. It is mainly metabolized by CYP2C9 and partly by CYP3A4 after oral administration. Many studies reported that CYP2C9 genetic polymorphism has significant effects on the pharmacokinetics of celecoxib and the occurrence of adverse drug reactions. The aim of this study was to develop a physiologically based pharmacokinetic (PBPK) model of celecoxib according to CYP2C9 genetic polymorphism for personalized pharmacotherapy. Initially, a clinical pharmacokinetic study was conducted where a single dose (200 mg) of celecoxib was administered to 39 healthy Korean subjects with CYP2C9*1/*1 or CYP2C9*1/*3 genotypes to obtain data for PBPK development. Based on the conducted pharmacokinetic study and a previous pharmacokinetic study involving subjects with CYP2C9*1/*13 and CYP2C9*3/*3 genotype, PBPK model for celecoxib was developed. A PBPK model for CYP2C9*1/*1 genotype group was developed and then scaled to other genotype groups (CYP2C9*1/*3, CYP2C9*1/*13 and CYP2C9*3/*3). After model development, model validation was performed with comparison of five pharmacokinetic studies. As a result, the developed PBPK model of celecoxib successfully described the pharmacokinetics of each CYP2C9 genotype group and its predicted values were within the acceptance criterion. Additionally, all the predicted values were within two-fold error range in comparison to the previous pharmacokinetic studies. This study demonstrates the possibility of determining the appropriate dosage of celecoxib for each individual through the PBPK modeling with CYP2C9 genomic information. This approach could contribute to the reduction of adverse drug reactions of celecoxib and enable precision medicine.

Body Composition and Dose-limiting Toxicity in Colorectal Cancer Chemotherapy Treatment; a Systematic Review of the Literature. Could Muscle Mass be the New Body Surface Area in Chemotherapy Dosing?

Chemotherapy dosing is traditionally based on body surface area calculations; however, these calculations ignore separate tissue compartments, such as the lean body mass (LBM), which is considered a big pool of drug distribution. In our era, colorectal cancer patients undergo a plethora of computed tomography scans as part of their diagnosis, staging and monitoring, which could easily be used for body composition analysis and LBM calculation, allowing for personalised chemotherapy dosing. This systematic review aims to evaluate the effect of muscle mass on dose-limiting toxicity (DLT), among different chemotherapy regimens used in colorectal cancer patients. This review was carried out according to the PRISMA guidelines. MEDLINE and EMBASE databases were searched from 1946 to August 2019. The primary search terms were 'sarcopenia', 'myopenia', 'chemotherapy toxicity', 'chemotherapy dosing', 'dose limiting toxicity', 'colorectal cancer', 'primary colorectal cancer' and 'metastatic colorectal cancer'. Outcomes of interest were - DLT and chemotoxicity related to body composition, and chemotherapy dosing on LBM. In total, 363 studies were identified, with 10 studies fulfilling the selection criteria. Seven studies were retrospective and three were prospective. Most studies used the same body composition analysis software but the chemotherapy regimens used varied. Due to marked study heterogeneity, quantitative data synthesis was not possible. Two studies described a toxicity cut-off value for 5-fluorouracil and one for oxaliplatin based on LBM. The rest of the studies showed an association between different body composition metrics and DLTs. Prospective studies are required with a larger colorectal cancer cohort, longitudinal monitoring of body composition changes during treatment, similar body composition analysis techniques, agreed cut-off values and standardised chemotherapy regimens. Incorporation of body composition analysis in the clinical setting will allow early identification of sarcopenic patients, personalised dosing based on their LBM and early optimisation of these patients undergoing chemotherapy.

A novel metabolic disorder in the degradation pathway of endogenous methanol due to a mutation in the gene of alcohol dehydrogenase

BACKGROUND

A small amount of methanol is produced endogenously in the human body but it is efficiently metabolized by alcohol dehydrogenase (ADH) and other enzymes, and the products eliminated without harm. In this study, we present a new entity of inborn error of methanol metabolism due to a mutation in the ADH1C gene coding for the γ subunit that is part of several ADH isoenzymes.

RESULTS

This disorder was discovered in an 11.58-year-old boy. During one 9-month hospital admission, he had periods of 1-4 days during which he was comatose, and between these periods he was sometimes verbose and euphoric, and had ataxia, dysarthria. Following hemodialysis treatments, he became conscious and appeared healthy. Organ evaluations and his laboratory tests were normal. Toxicological evaluation of his blood showed a high methanol level [12.2 mg/dL (3.8 mmol/L), normal range up to 3.5 mg/dL (1.09 mmol/L) while the formaldehyde level was undetectable. The finding of liver function tests that were within normal limits, coupled with a normal eye examination and size of the liver, elevated blood methanol levels and an undetectable formaldehyde level, suggested ADH insufficiency. Adding zinc to the drug regimen 15 mg/daily dramatically reduced the patient's methanol level and alleviated the abnormal symptoms. When zinc supplementation was discontinued, the patient relapsed into a coma and hemodialysis was once again required. A homozygous mutation in ADH1C gene located at exon 3 was found, and both parents were heterozygous for this mutation.

CONCLUSION

Accumulation of methanol due to mutation in ADH1C gene may result in drunkenness and ataxia, and leads to coma. This condition can be successfully treated with zinc supplementation as the cofactor of ADH.

Human superoxide dismutase 1 attenuates quinoneimine metabolite formation from mefenamic acid

Mefenamic acid (MFA), one of the nonsteroidal anti-inflammatory drugs (NSAIDs), sometimes causes liver injury. Quinoneimines formed by cytochrome P450 (CYP)-mediated oxidation of MFA are considered to be causal metabolites of the toxicity and are detoxified by glutathione conjugation. A previous study reported that NAD(P)H:quinone oxidoreductase 1 (NQO1) can reduce the quinoneimines, but NQO1 is scarcely expressed in the human liver. The purpose is to identify enzyme(s) responsible for the decrease in MFA-quinoneimine formation in the human liver. The formation of MFA-quinoneimine by recombinant CYP1A2 and CYP2C9 was significantly decreased by the addition of human liver cytosol, and the extent of the decrease in the metabolite formed by CYP1A2 was larger than that by CYP2C9. By column chromatography, superoxide dismutase 1 (SOD1) was identified from the human liver cytosol as an enzyme decreasing MFA-quinoneimine formation. Addition of recombinant SOD1 into the reaction mixture decreased the formation of MFA-quinoneimine from MFA by recombinant CYP1A2. By a structure-activity relationship study, we found that SOD1 decreased the formation of quinoneimines from flufenamic acid and tolfenamic acid, but did not affect those produced from acetaminophen, amodiaquine, diclofenac, and lapatinib. Thus, SOD1 may selectively decrease the quinoneimine formation from fenamate-class NSAIDs. To examine whether SOD1 can attenuate cytotoxicity caused by MFA, siRNA for SOD1 was transfected into CYP1A2-overexpressed HepG2 cells. The leakage of lactate dehydrogenase caused by MFA treatment was significantly increased by knockdown of SOD1. In conclusion, we found that SOD1 can serve as a detoxification enzyme for quinoneimines to protect from drug-induced toxicity.

Altered Metabolism of Polycyclic Aromatic Hydrocarbons by UDP-Glycosyltransferase 3A2 Missense Variants

The UDP-glycosyltransferase (UGT) family of enzymes are important in the metabolism of a variety of exogenous substances including polycyclic aromatic hydrocarbons (PAHs), a potent class of environmental carcinogens. As compared to the majority of UGT enzymes, which utilize UDP-glucuronic acid as a cosubstrate, UGT3A2 utilizes alternative cosubstrates (UDP-glucose and UDP-xylose). UGT3A2 is expressed in aerodigestive tract tissues and was highly active against multiple PAHs with both cosubstrates. The goal of the present study was to assess the functional effects of UGT3A2 missense variants (MAF ≥ 0.005) on PAH metabolism and the utilization of cosubstrates. The glycosylation activity (V_max/K_m) of all variants against simple PAHs using both cosubstrates was significantly (P < 0.05) decreased by 42-100% when compared to wild-type UGT3A2. When utilizing UDP-glucose, the variant isoforms exhibited up to a 362-fold decrease in V_max/K_m when compared to wild-type UGT3A2, with a 3.1- to 14-fold decrease for D140N, A344T, and S435Y, a 24- and 43-fold decrease for A436T and R445C, respectively, and a 147- and 362-fold decrease for Y474C and Y74N, respectively. When utilizing UDP-xylose, the variants exhibited up to a 4.0-fold decrease in V_max/K_m when compared to wild-type UGT3A2; Y74N did not exhibit activity, and Y474C did not reach saturation (K_m > 4000 μM). Additionally, both wild-type and variant UGT3A2 exhibited a significant (P < 0.05) difference in their utilization of UDP-glucose vs UDP-xylose as cosubstrates using 1-OH-pyrene as substrate. These data suggest that UGT3A2 missense variants decrease the detoxification of PAHs, potentially resulting in altered individual risk for PAH-related cancers.

The Ontogeny of UDP-glucuronosyltransferase Enzymes, Recommendations for Future Profiling Studies and Application Through Physiologically Based Pharmacokinetic Modelling

Limited understanding of drug pharmacokinetics in children is one of the major challenges in paediatric drug development. This is most critical in neonates and infants owing to rapid changes in physiological functions, especially in the activity of drug-metabolising enzymes. Paediatric physiologically based pharmacokinetic models that integrate ontogeny functions for cytochrome P450 enzymes have aided our understanding of drug exposure in children, including those under the age of 2 years. Paediatric physiologically based pharmacokinetic models have consequently been recognised by the European Medicines Agency and the US Food and Drug Administration as innovative tools in paediatric drug development and regulatory decision making. However, little is currently known about age-related changes in UDP-glucuronosyltransferase-mediated metabolism, which represents the most important conjugation reaction for xenobiotics. Therefore, the objective of the review was to conduct a thorough literature survey to summarise our current understanding of age-related changes in UDP-glucuronosyltransferases as well as associated clinical and experimental sources of variance. Our findings indicate that there are distinct differences in UDP-glucuronosyltransferase expression and activity between isoforms for different age groups. In addition, there is substantial variability between individuals and laboratories reported for human liver microsomes, which results in part from a lack of standardised experimental conditions. Therefore, we provide a number of best practice recommendations for experimental conditions, which ultimately may help improve the quality of data used for quantitative clinical pharmacology approaches, and thus for safe and effective pharmacotherapy in children.

Sources of all-trans retinal oxidation independent of the aldehyde dehydrogenase 1A isozymes exist in the postnatal testis†

Despite the essential role of the active metabolite of vitamin A, all-trans retinoic acid (atRA) in spermatogenesis, the enzymes, and cellular populations responsible for its synthesis in the postnatal testis remain largely unknown. The aldehyde dehydrogenase 1A (ALDH1A) family of enzymes residing within Sertoli cells is responsible for the synthesis of atRA, driving the first round of spermatogenesis. Those studies also revealed that the atRA required to drive subsequent rounds of spermatogenesis is possibly derived from the ALDH1A enzymes residing within the meiotic and post-meiotic germ cells. Three ALDH1A isozymes (ALDH1A1, ALDH1A2, and ALDH1A3) are present in the testis. Although, ALDH1A1 is expressed in adult Sertoli cells and is suggested to contribute to the atRA required for the pre-meiotic transitions, ALDH1A2 is proposed to be the essential isomer involved in testicular atRA biosynthesis. In this report, we first examine the requirement for ALDH1A2 via the generation and analysis of a conditional Aldh1a2 germ cell knockout and a tamoxifen-induced Aldh1a2 knockout model. We then utilized the pan-ALDH1A inhibitor (WIN 18446) to test the collective contribution of the ALDH1A enzymes to atRA biosynthesis following the first round of spermatogenesis. Collectively, our data provide the first in vivo evidence demonstrating that animals severely deficient in ALDH1A2 postnatally proceed normally through spermatogenesis. Our studies with a pan-ALDH1A inhibitor (WIN 18446) also suggest that an alternative source of atRA biosynthesis independent of the ALDH1A enzymes becomes available to maintain atRA levels for several spermatogenic cycles following an initial atRA injection.

The role of xenobiotic-metabolizing enzymes in the placenta: a growing research field

Introduction: The placenta is a temporary and unique organ that allows for the physical connection between a mother and fetus; this organ regulates the transport of gases and nutrients mediating the elimination of waste products contained in the fetal circulation. The placenta performs metabolic and excretion functions, on the basis of multiple enzymatic systems responsible for the oxidation, reduction, hydrolysis, and conjugation of xenobiotics. These mechanisms give the placenta a protective role that limits the fetal exposure to harmful compounds. During pregnancy, some diseases require uninterrupted treatment even if it is detrimental to the fetus. Drugs and other xenobiotics alter gene expression in the placenta with repercussions for the fetus and mother's well-being.Areas covered: This review provides a brief description of the human placental structure and function, the main drug and xenobiotic transporters and metabolizing enzymes, placenta-metabolized substrates, and alterations in gene expression that the exposure to xenobiotics may cause.Expert opinion: Research should be focused on the identification and validation of biological markers for the assessment of the harmful effects of some drugs in pregnancy, including the evaluation of polymorphisms and methylation patterns in chorionic villous samples and/or amniotic fluid.

UDP-Glycosyltransferase 3A Metabolism of Polycyclic Aromatic Hydrocarbons: Potential Importance in Aerodigestive Tract Tissues

Polycyclic aromatic hydrocarbons (PAHs) are potent carcinogens and are a primary risk factor for the development of lung and other aerodigestive tract cancers in smokers. The detoxification of PAHs by glucuronidation is well-characterized for the UDP-glycosyltransferase (UGT) 1A, 2A, and 2B subfamilies; however, the role of the UGT3A subfamily in PAH metabolism remains poorly understood. UGT3A enzymes are functionally distinct from other UGT subfamilies (which use UDP-glucuronic acid as a cosubstrate) due to their utilization of alternative cosubstrates (UDP-N-acetylglucosamine for UGT3A1, and UDP-glucose and UDP-xylose for UGT3A2). The goal of the present study was to characterize UGT3A glycosylation activity against PAHs and examine their expression in human aerodigestive tract tissues. In vitro metabolism assays using UGT3A2-overexpressing cell microsomes indicated that UGT3A2 exhibits glycosylation activity against all of the simple and complex PAHs tested. The V max/K m ratios for UGT3A2 activity with UDP-xylose versus UDP-glucose as the cosubstrate ranged from 0.65 to 4.4 for all PAHs tested, demonstrating that PAH glycosylation may be occurring at rates up to 4.4-fold higher with UDP-xylose than with UDP-glucose. Limited glycosylation activity was observed against PAHs with UGT3A1-overexpressing cell microsomes. While UGT3A2 exhibited low levels of hepatic expression, it was shown by western blot analysis to be widely expressed in aerodigestive tract tissues. Conversely, UGT3A1 exhibited the highest expression in liver with lower expression in aerodigestive tract tissues. These data suggest that UGT3A2 plays an important role in the detoxification of PAHs in aerodigestive tract tissues, and that there may be cosubstrate-dependent differences in the detoxification of PAHs by UGT3A2. SIGNIFICANCE STATEMENT: UGT3A2 is highly active against PAHs with either UDP-glucose or UDP-xylose as a cosubstrate. UGT3A1 exhibited low levels of activity against PAHs. UGT3A1 is highly expressed in liver while UGT3A2 is well expressed in extrahepatic tissues. UGT3A2 may be an important detoxifier of PAHs in humans.

Transcriptomic modifications of the thyroid gland upon exposure to phytosanitary-grade fipronil: Evidence for the activation of compensatory pathways

Fipronil is a phenylpyrazole insecticide used for the control of a variety of pest for domestic, veterinary and agricultural uses. Fipronil exposure is associated to thyroid disruption in the rat. It increases thyroid hormone (TH) hepatic clearance. The effect on thyroxine (T4) clearance is about four fold higher than the effect on T4 plasma concentrations suggesting that the thyroid gland might develop compensatory mechanisms. The aim of this study was to document the potential effects of fipronil treatment on the thyroid transcriptome together with its effects on TSH and TH blood levels under well characterized internal exposure to fipronil and its main metabolite fipronil sulfone. Fipronil (3 mg/kg/d by gavage for 14 days) clearance increased while its half-life decreased (about 10 fold) throughout treatment. Fipronil treatment in adult female rats significantly decreased total T4 and free triiodothyronine (T3) concentrations. Key genes related to thyroid hormone synthesis and/or cellular dynamic were modulated by fipronil exposure. RT-PCR confirmed that thyroglobulin gene expression was upregulated. A trend toward higher Na/I symporter expression was also noted, while sulfotransferase 1a1 gene expression was down-regulated. The expression of genes potentially involved in thyroid cell dynamic were upregulated (e.g. prostaglandin synthase 1, amphiregulin and Rhoa). Our results indicate that both pathways of TH synthesis and thyroid cell dynamics are transcriptional targets of fipronil and/or its main sulfone metabolite. The underlying mechanisms remain to be elucidated.

Detection of ALDH3B2 in Human Placenta

Aldehyde dehydrogenase 3B2 (ALDH3B2) gene contains a premature termination codon, which can be skipped or suppressed resulting in full-length protein expression. Alternatively, the longest putative open reading frame starting with the second in-frame start codon would encode short isoform. No unequivocal evidence of ALDH3B2 expression in healthy human tissues is available. The aim of this study was to confirm its expression in human placenta characterized by the highest ALDH3B2 mRNA abundance. ALDH3B2 DNA and mRNA were sequenced. The expression was investigated using western blot. The identity of the protein was confirmed using mass spectrometry (MS). The predicted tertiary and quaternary structures, subcellular localization, and phosphorylation sites were assessed using bioinformatic analyses. All DNA and mRNA isolates contained the premature stop codon. In western blot analyses, bands corresponding to the mass of full-length protein were detected. MS analysis led to the identification of two unique peptides, one of which is encoded by the nucleotide sequence located upstream the second start codon. Bioinformatic analyses suggest cytoplasmic localization and several phosphorylation sites. Despite premature stop codon in DNA and mRNA sequences, full-length ALDH3B2 was found. It can be formed as a result of premature stop codon readthrough, complex phenomenon enabling stop codon circumvention.

Xenobiotica-metabolizing enzymes in the lung of experimental animals, man and in human lung models

The xenobiotic metabolism in the lung, an organ of first entry of xenobiotics into the organism, is crucial for inhaled compounds entering this organ intentionally (e.g. drugs) and unintentionally (e.g. work place and environmental compounds). Additionally, local metabolism by enzymes preferentially or exclusively occurring in the lung is important for favorable or toxic effects of xenobiotics entering the organism also by routes other than by inhalation. The data collected in this review show that generally activities of cytochromes P450 are low in the lung of all investigated species and in vitro models. Other oxidoreductases may turn out to be more important, but are largely not investigated. Phase II enzymes are generally much higher with the exception of UGT glucuronosyltransferases which are generally very low. Insofar as data are available the xenobiotic metabolism in the lung of monkeys comes closed to that in the human lung; however, very few data are available for this comparison. Second best rate the mouse and rat lung, followed by the rabbit. Of the human in vitro model primary cells in culture, such as alveolar macrophages and alveolar type II cells as well as the A549 cell line appear quite acceptable. However, (1) this generalization represents a temporary oversimplification born from the lack of more comparable data; (2) the relative suitability of individual species/models is different for different enzymes; (3) when more data become available, the conclusions derived from these comparisons quite possibly may change.

The UGTome: The expanding diversity of UDP glycosyltransferases and its impact on small molecule metabolism

The UDP glycosyltransferase (UGT) superfamily of enzymes is responsible for the metabolism and clearance of thousands of lipophilic chemicals including drugs, toxins and endogenous signaling molecules. They provide a protective interface between the organism and its chemical-rich environment, as well as controlling critical signaling pathways to maintain healthy tissue function. UGTs are associated with drug responses and interactions, as well as a wide range of diseases including cancer. The human genome contains 22 UGT genes; however as befitting their exceptionally diverse substrate ranges and biological activities, the output of these UGT genes is functionally diversified by multiple processes including alternative splicing, post-translational modification, homo- and hetero-oligomerization, and interactions with other proteins. All UGT genes are subject to extensive alternative splicing generating variant/truncated UGT proteins with altered functions including the capacity to dominantly modulate/inhibit cognate full-length forms. Heterotypic oligomerization of different UGTs can alter kinetic properties relative to monotypic complexes, and potentially produce novel substrate specificities. Moreover, the recently profiled interactions of UGTs with non-UGT proteins may facilitate coordination between different metabolic processes, as well as providing opportunities for UGTs to engage in novel 'moonlighting' functions. Herein we provide a detailed and comprehensive review of all known modes of UGT functional diversification and propose a UGTome model to describe the resulting expansion of metabolic capacity and its potential to modulate drug/xenobiotic responses and cell behaviours in normal and disease contexts.

Aldehyde oxidase and its role as a drug metabolizing enzyme

Aldehyde oxidase (AO) is a cytosolic enzyme that belongs to the family of structurally related molybdoflavoproteins like xanthine oxidase (XO). The enzyme is characterized by broad substrate specificity and marked species differences. It catalyzes the oxidation of aromatic and aliphatic aldehydes and various heteroaromatic rings as well as reduction of several functional groups. The references to AO and its role in metabolism date back to the 1950s, but the importance of this enzyme in the metabolism of drugs has emerged in the past fifteen years. Several reviews on the role of AO in drug metabolism have been published in the past decade indicative of the growing interest in the enzyme and its influence in drug metabolism. Here, we present a comprehensive monograph of AO as a drug metabolizing enzyme with emphasis on marketed drugs as well as other xenobiotics, as substrates and inhibitors. Although the number of drugs that are primarily metabolized by AO are few, the impact of AO on drug development has been extensive. We also discuss the effect of AO on the systemic exposure and clearance these clinical candidates. The review provides a comprehensive analysis of drug discovery compounds involving AO with the focus on developmental candidates that were reported in the past five years with regards to pharmacokinetics and toxicity. While there is only one known report of AO-mediated clinically relevant drug-drug interaction (DDI), a detailed description of inhibitors and inducers of AO known to date has been presented here and the potential risks associated with DDI. The increasing recognition of the importance of AO has led to significant progress in predicting the site of AO-mediated metabolism using computational methods. Additionally, marked species difference in expression of AO makes it is difficult to predict human clearance with high confidence. The progress made towards developing in vivo, in vitro and in silico approaches for predicting AO metabolism and estimating human clearance of compounds that are metabolized by AO have also been discussed.

The trichloroethylene metabolite S-(1,2-dichlorovinyl)-L-cysteine induces progressive mitochondrial dysfunction in HTR-8/SVneo trophoblasts

Trichloroethylene is an industrial solvent and common environmental pollutant. Despite efforts to ban trichloroethylene, its availability and usage persist globally, constituting a hazard to human health. Recent studies reported associations between maternal trichloroethylene exposure and increased risk for low birth weight. Despite these associations, the toxicological mechanism underlying trichloroethylene adverse effects on pregnancy remains largely unknown. The trichloroethylene metabolite S-(1,2-dichlorovinyl)-L-cysteine (DCVC) induces mitochondrial-mediated apoptosis in a trophoblast cell line. To gain further understanding of mitochondrial-mediated DCVC placental toxicity, this study investigated the effects of DCVC exposure on mitochondrial function using non-cytolethal concentrations in placental cells. Human trophoblasts, HTR-8/SVneo, were exposed in vitro to a maximum of 20 μM DCVC for up to 12 h. Cell-based oxygen consumption and extracellular acidification assays were used to evaluate key aspects of mitochondrial function. Following 6 h of exposure to 20 μM DCVC, elevated oxygen consumption, mitochondrial proton leak and sustained energy coupling deficiency were observed. Similarly, 12 h of exposure to 20 μM DCVC decreased mitochondrial-dependent basal, ATP-linked and maximum oxygen consumption rates. Using the fluorochrome TMRE, dissipation of mitochondrial membrane potential was detected after a 12-h exposure to 20 μM DCVC, and (±)-α-tocopherol, a known suppressor of lipid peroxidation, attenuated DCVC-stimulated mitochondrial membrane depolarization but failed to rescue oxygen consumption perturbations. Together, these results suggest that DCVC caused progressive mitochondrial dysfunction, resulting in lipid peroxidation-associated mitochondrial membrane depolarization. Our findings contribute to the biological plausibility of DCVC-induced placental impairment and provide new insights into the role of the mitochondria in DCVC-induced toxicity.

The UDP-Glycosyltransferase (UGT) Superfamily: New Members, New Functions, and Novel Paradigms

UDP-glycosyltransferases (UGTs) catalyze the covalent addition of sugars to a broad range of lipophilic molecules. This biotransformation plays a critical role in elimination of a broad range of exogenous chemicals and by-products of endogenous metabolism, and also controls the levels and distribution of many endogenous signaling molecules. In mammals, the superfamily comprises four families: UGT1, UGT2, UGT3, and UGT8. UGT1 and UGT2 enzymes have important roles in pharmacology and toxicology including contributing to interindividual differences in drug disposition as well as to cancer risk. These UGTs are highly expressed in organs of detoxification (e.g., liver, kidney, intestine) and can be induced by pathways that sense demand for detoxification and for modulation of endobiotic signaling molecules. The functions of the UGT3 and UGT8 family enzymes have only been characterized relatively recently; these enzymes show different UDP-sugar preferences to that of UGT1 and UGT2 enzymes, and to date, their contributions to drug metabolism appear to be relatively minor. This review summarizes and provides critical analysis of the current state of research into all four families of UGT enzymes. Key areas discussed include the roles of UGTs in drug metabolism, cancer risk, and regulation of signaling, as well as the transcriptional and posttranscriptional control of UGT expression and function. The latter part of this review provides an in-depth analysis of the known and predicted functions of UGT3 and UGT8 enzymes, focused on their likely roles in modulation of levels of endogenous signaling pathways.

Selectivity of Dietary Phenolics for Inhibition of Human Monoamine Oxidases A and B

Monoamine oxidases (MAOs) regulate local levels of neurotransmitters such as dopamine, norepinephrine, and serotonin and thus have been targeted by drugs for the treatment of certain CNS disorders. However, recent studies have shown that these enzymes are upregulated with age in nervous and cardiac tissues and may be involved in degeneration of these tissues, since their metabolic mechanism releases hydrogen peroxide leading to oxidative stress. Thus, targeting these enzymes may be a potential anti-aging strategy. The purpose of this study was to compare the MAO inhibition and selectivity of selected dietary phenolic compounds, using a previously validated assay that would avoid interference from the compounds. Kynuramine metabolism by human recombinant MAO-A and MAO-B leads to formation of 4-hydroxyquinoline, with Vmax values of 10.2±0.2 and 7.35±0.69 nmol/mg/min, respectively, and Km values of 23.1±0.8 μM and 18.0±2.3 μM, respectively. For oral dosing and interactions with the gastrointestinal tract, curcumin, guaiacol, isoeugenol, pterostilbene, resveratrol, and zingerone were tested at their highest expected luminal concentrations from an oral dose. Each of these significantly inhibited both enzymes except for zingerone, which only inhibited MAO-A. The IC50 values were determined, and selectivity indices (MAO-A/MAO-B IC50 ratios) were calculated. Resveratrol and isoeugenol were selective for MAO-A, with IC50 values of 0.313±0.008 and 3.72±0.20 μM and selectivity indices of 50.5 and 27.4, respectively. Pterostilbene was selective for MAO-B, with IC50 of 0.138±0.013 μM and selectivity index of 0.0103. The inhibition of resveratrol (MAO-A) and pterostilbene (MAO-B) was consistent with competitive time-independent mechanisms. Resveratrol 4'-glucoside was the only compound which inhibited MAO-A, but itself, resveratrol 3-glucoside, and pterostilbene 4'-glucoside failed to inhibit MAO-B. Additional studies are needed to establish the effects of these compounds on MAO-A and/or MAO-B in humans.

Comparison of Rat and Human Pulmonary Metabolism Using Precision-cut Lung Slices (PCLS)

BACKGROUND

Although the liver is the primary organ of drug metabolism, the lungs also contain drug-metabolizing enzymes and may, therefore, contribute to the elimination of drugs. In this investigation, the Precision-cut Lung Slice (PCLS) technique was standardized with the aims of characterizing and comparing rat and human pulmonary drug metabolizing activity.

METHOD

Due to the limited availability of human lung tissue, standardization of the PCLS method was performed with rat lung tissue. Pulmonary enzymatic activity was found to vary significantly with rat age and rat strain. The Dynamic Organ Culture (DOC) system was superior to well-plates for tissue incubations, while oxygen supply appeared to have a limited impact within the 4h incubation period used here.

RESULTS

The metabolism of a range of phase I and phase II probe substrates was assessed in rat and human lung preparations. Cytochrome P450 (CYP) activity was relatively low in both species, whereas phase II activity appeared to be more significant.

CONCLUSION

PCLS is a promising tool for the investigation of pulmonary drug metabolism. The data indicates that pulmonary CYP activity is relatively low and that there are significant differences in enzyme activity between rat and human lung.

Biomarkers for In Vivo Assessment of Transporter Function

Drug-drug interactions are a major concern not only during clinical practice, but also in drug development. Due to limitations of in vitro-in vivo predictions of transporter-mediated drug-drug interactions, multiple clinical Phase I drug-drug interaction studies may become necessary for a new molecular entity to assess potential drug interaction liabilities. This is a resource-intensive process and exposes study participants, who frequently are healthy volunteers without benefit from study treatment, to the potential risks of a new drug in development. Therefore, there is currently a major interest in new approaches for better prediction of transporter-mediated drug-drug interactions. In particular, researchers in the field attempt to identify endogenous compounds as biomarkers for transporter function, such as hexadecanedioate, tetradecanedioate, coproporphyrins I and III, or glycochenodeoxycholate sulfate for hepatic uptake via organic anion transporting polypeptide 1B or N¹-methylnicotinamide for multidrug and toxin extrusion protein-mediated renal secretion. We summarize in this review the currently proposed biomarkers and potential limitations of the substances identified to date. Moreover, we suggest criteria based on current experiences, which may be used to assess the suitability of a biomarker for transporter function. Finally, further alternatives and supplemental approaches to classic drug-drug interaction studies are discussed.

Expression of clinically relevant drug-metabolizing enzymes along the human intestine and their correlation to drug transporters and nuclear receptors: An intra-subject analysis

The oral bioavailability of many drugs is highly influenced not only by hepatic but also by intestinal biotransformation. To estimate the impact of intestinal phase I and II metabolism on oral drug absorption, knowledge on the expression levels of the respective enzymes is an essential prerequisite. In addition, the potential interplay of metabolism and transport contributes to drug disposition. Both mechanisms may be subjected to coordinative regulation by nuclear receptors, leading to unwanted drug-drug interactions due to induction of intestinal metabolism and transport. Thus, it was the aim of this study to comprehensively analyse the regional expression of clinically relevant phase I and II enzymes along the entire human intestine and to correlate these data to expression data of drug transporters and nuclear receptors of pharmacokinetic relevance. Gene expression of 11 drug-metabolizing enzymes (CYP2B6, 2C8, 2C9, 2C19, 2D6, 3A4, 3A5, SULT1A, UGT1A, UGT2B7, UGT2B15) was studied in duodenum, jejunum, ileum and colon from six organ donors by real-time RT-PCR. Enzyme expression was correlated with expression data of the nuclear receptors PXR, CAR and FXR as well as drug transporters observed in the same cohort. Intestinal expression of all studied metabolizing enzymes was significantly higher in the small intestine compared to colonic tissue. CYP2B6, CYP2C9, CYP2C19, CYP2D6, CYP3A4/5, SULT1A, UGT1A and UGT2B7 expression increased from the duodenum to jejunum but was markedly lower in the ileum. In the small intestine, that is, the predominant site of drug absorption, the highest expression has been observed for CYP3A4, CYP2C9, SULT1A and UGT1A. In addition, significant correlations were found between several enzymes and PXR as well as ABC transporters in the small intestine. In conclusion, the observed substantial site-dependent intestinal expression of several enzymes may explain regional differences in intestinal drug absorption. The detected correlations between intestinal enzymes, transporters and nuclear receptors provide indirect evidence for their coordinative expression, regulation and function in the human small intestine.