Ontogeny of Human Liver Aldehyde Oxidase: Developmental Changes and Implications for Drug Metabolism.

Despite the increasing importance of aldehyde oxidase (AO) in the drug metabolism of clinical candidates, ontogeny data for AO are limited. The objective of our study was to characterize the age-dependent AO content and activity in the human liver cytosolic fraction (HLC) and human hepatocytes (HH). HLC (n = 121 donors) and HH (n = 50 donors) were analyzed for (1) AO protein content by quantitative proteomics and (2) enzyme activity using carbazeran as a probe substrate. AO activity showed high technical variability and poor correlation with the content in HLC samples, whereas hepatocyte samples showed a strong correlation between the content and activity. Similarly, AO content and activity showed no significant age-dependent differences in HLC samples, whereas the average AO content and activity in hepatocytes increased significantly (∼20-40-fold) from the neonatal levels (0-28 days). Based on the hepatocyte data, the age at which 50% of the adult AO content is reached (age50) was 3.15 years (0.32-13.97 years, 95% CI). Metabolite profiling of carbazeran revealed age-dependent metabolic switching and the role of non-AO mechanisms (glucuronidation and desmethylation) in carbazeran elimination. The content-activity correlation in hepatocytes improved significantly (R2 = 0.95; p < 0.0001) in samples showing <10% contribution of glucuronidation toward the overall metabolism, confirming that AO-mediated oxidation and glucuronidation are the key routes of carbazeran metabolism. Considering the confounding effect of glucuronidation on AO activity, AO content-based ontogeny data are a more direct reflection of developmental changes in protein expression. The comprehensive ontogeny data of AO in HH samples are more reliable than HLC data, which are important for developing robust physiologically based pharmacokinetic models for predicting AO-mediated metabolism in children.

Ontogeny of Scaling Factors for Pediatric Physiologically Based Pharmacokinetic Modeling and Simulation: Cytosolic Protein Per Gram of Liver.

Scaling factors are necessary for translating in vitro drug biotransformation data to in vivo clearance values in physiologically-based pharmacokinetic modeling and simulation. Values for microsomal protein per gram of liver are available from several sources for use as a scaling factor to estimate hepatic clearance from microsomal drug biotransformation data. However, data regarding the distribution of cytosolic protein per gram of liver (CPPGL) values across the lifespan are limited, and sparse pediatric data have been published to date. Thus, CPPGL was determined in 160 liver samples from pediatric (n = 129) and adult (n = 31) donors obtained from multiple sources: the University of Maryland Brain and Tissue Bank, tissue retrieval services at the University of Minnesota and University of Pittsburgh, and Sekisui-XenoTech. Tissues were homogenized and subjected to differential centrifugation to isolate cytosolic fractions. Cytosolic protein content was determined by BCA assay. CPPGL varied from two- to sixfold within each age group/developmental stage. Tissue source and sex did not contribute substantially to variability in protein content. Regression analyses revealed minimal change in CPPGL over the first two decades of life (logCPPGL increases 0.1 mg/g per decade). A mean ± S.D. CPPGL value of 44.4 ± 17.4 mg/g or median 41.0 mg/g is representative of values observed between birth and early adulthood (0-18 years, n = 129). SIGNIFICANCE STATEMENT: Cytosolic protein per gram of liver (CPPGL) is a scaling factor required for physiologically based pharmacokinetic modeling and simulation of drug biotransformation by cytosolic enzymes, but pediatric data are limited. Although CPPGL varies from two- to sixfold within developmental stages, a value of 44.4 ± 17.4 mg/g (mean ± S.D.) is representative of the pediatric period (0-18 years, n = 129).

Isolation and Identification of 3,4-Seco-Solanidine-3,4-dioic acid (SSDA) as a Urinary Biomarker of Cytochrome P450 2D6 (CYP2D6) Activity.

Cytochrome P450 2D6 (CYP2D6), is responsible for the metabolism and elimination of approximately 25% of clinically used drugs, including antidepressants and antipsychotics, and its activity varies considerably on a population basis primary due to genetic variation. CYP2D6 phenotype can be assessed in vivo following administration of an exogenous probe compound, such as dextromethorphan or debrisoquine, but use of a biomarker that does not require administration of an exogenous compound (i.e., drug) has considerable appeal for assessing CYP2D6 activity in vulnerable populations, such as children. The goal of this study was to isolate, purify and identify an "endogenous" urinary biomarker (M1; m/z 444.3102) of CYP2D6 activity reported previously. Several chromatographic separation techniques (reverse phase HPLC, cation exchange and analytical reverse phase UPLC) were used to isolate and purify 96 µg of M1 from 40 L of urine. Subsequently, 1D and 2D NMR, and functional group modification reactions were used to elucidate its structure. Analysis of mass spectrometry and NMR data revealed M1 to have similar spectroscopic features to the nitrogen-containing steroidal alkaloid, solanidine. 2D NMR characterization by HMBC, COSY, TOCSY, and HSQC-TOCSY proved to be invaluable in the structural elucidation of M1; derivatization of M1 revealed the presence of two carboxylic acid moieties. M1 was determined to be a steroidal alkaloid with a solanidine backbone that had undergone C-C bond scission to yield 3,4-seco-solanidine-3,4-dioic acid (SSDA). SSDA may have value as a dietary biomarker of CYP2D6 activity in populations where potato consumption is common. Significance Statement Endogenous biomarkers of processes involved in drug disposition and response may allow improved individualization of drug treatment, especially in vulnerable populations, such as children. Given that several CYP2D6 substrates are commonly used in pediatrics and the ubiquitous nature of potato consumption in western diets, SSDA has considerable appeal as non-invasive biomarker of CYP2D6 activity to guide treatment with CYP2D6 substrates in children and adults.

Ontogeny of Scaling Factors for Pediatric Physiology-Based Pharmacokinetic Modeling and Simulation: Microsomal Protein Per Gram of Liver.

Microsomal protein per gram of liver (MPPGL) is an important scaling factor for bottom-up physiology-based pharmacokinetic modeling and simulation, but data in pediatrics are limited. Therefore, MPPGL was determined in 160 liver samples from pediatric (n = 129) and adult (n = 31) donors obtained from four sources: the University of Maryland Brain and Tissue Bank (UMBTB), tissue retrieval services at the University of Minnesota and University of Pittsburgh, and Sekisui-Xenotech. Tissues were homogenized and subjected to differential centrifugation to prepare microsomes, and cytochrome c reductase activities in tissue homogenates and microsomes were used to estimate cytochrome P450 reductase (POR) activity as a marker of microsomal recovery; microsomal POR content was also assessed by quantitative proteomics. MPPGL values varied 5- to 10-fold within various age groups/developmental stages, and tissue source was identified as a contributing factor. Using a "trimmed" dataset comprised of samples ranging from 3 to 18 years of age common to the four sources, POR protein abundance and activity in microsomes and POR activity in homogenates was lower in UMBTB samples (autopsy) compared with other sources (perfused/flash-frozen). Regression analyses revealed that the UMBTB samples were driving an apparent age effect as no effect of age on log-transformed MPPGL values was observed when the UMBTB samples were excluded. We conclude that a mean±SD MPPGL value of 30.4±1.7 mg/g is representative between one month postnatal age and early adulthood. Potential source effects should be considered for studies involving tissue samples from multiple sources with different procurement and processing procedures. SIGNIFICANCE STATEMENT: Microsomal protein per gram of liver (MPPGL) is an important scaling factor for bottom up PBPK modeling and simulation, but data in pediatrics are limited. Although MPPGL varies 5- to 10-fold at a given developmental stage, a value of 30.4 ± 1.7 mg/g (mean ± SD) is representative between one month postnatal age and early adulthood. However, when tissue samples are obtained from multiple sources, different procurement and processing procedures may influence the results and should be taken into consideration.

The Impact of Age and Genetics on Naltrexone Biotransformation.

Naltrexone, an opioid antagonist primarily metabolized by aldo-keto reductase 1C4 (AKR1C4), treats pediatric conditions involving compulsiveness (e.g., autism spectrum, Prader-Willi, eating disorders, non-suicidal self-injury). Pharmacokinetic variability is apparent in adults, yet no data are available for children. This study aimed to examine the impact of age and genetic variation on naltrexone biotransformation. Human liver cytosol (HLC) samples (n = 158) isolated from children and adult organ donors were incubated with therapeutically relevant concentrations of naltrexone (0.1, 1 µM). Naltrexone biotransformation was determined by ultraperformance mass spectrometry quantification of the primary metabolite, 6-beta-naltrexol (6ßN), and 6ßN formation rates (pmol/mg protein/min) were calculated. HLCs from organ donors, age range 0-79 y (mean 16.0 ± 18.2 y), 37% (n = 60) female, 20% (n = 33) heterozygous and 1.2% (n = 2) homozygous for co-occurring AKR1C4 variants (S145C/L311V) showed >200-fold range in 6ßN formation (0.37-76.5 pmol/mg protein/min). Source of donor samples was found to be a substantial contributor to variability. Model estimates for a trimmed data set of source-adjusted pediatric samples (aged 0-18 y) suggested that AKR1C4 genetic variation, age, and sex explained 36% of the variability in 6ßN formation. Although activity increased steadily from birth and peaked in middle childhood (2-5 years), genetic variation (S145C/L311V) demonstrated a greater effect on activity than did age. Naltrexone biotransformation is highly variable in pediatric and adult livers and can be partly accounted for by individual factors feasible to obtain (e.g., genetic variability, age, sex). These data may inform a precision therapeutics approach (e.g., exposure optimization) to further study Naltrexone responsiveness in children and adults. SIGNIFICANCE STATEMENT: Biotransformation of the commonly used opioid antagonist naltrexone is highly variable and may contribute to reduced therapeutic response. Age, sex, and genetic variation in the drug-metabolizing enzyme, AKR1C4, are potential factors contributing to this variability. In pediatric samples, genetic variation (S145C/L311V) demonstrates a greater impact on activity than age. Additionally, the source of donor samples was identified as an important contributor and must be accounted for to confidently elucidate the biological variables most impactful to drug biotransformation.

Interindividual Variability and Differential Tissue Abundance of Mitochondrial Amidoxime Reducing Component Enzymes in Humans.

Mitochondrial amidoxime-reducing component (mARC) enzymes are molybdenum-containing proteins that metabolize a number of endobiotics and xenobiotics. The interindividual variability and differential tissue abundance of mARC1 and mARC2 were quantified using targeted proteomics in three types of tissue fractions: 1) pediatric liver tissue homogenates, 2) total membrane fraction of the paired liver and kidney samples from pediatric and adult donors, and 3) pooled S9 fractions of the liver, intestine, kidney, lung, and heart. The absolute levels of mARC1 and mARC2 in the pediatric liver homogenate were 40.08 ± 4.26 and 24.58 ± 4.02 pmol/mg homogenate protein, respectively, and were independent of age and sex. In the total membrane fraction of the paired liver and kidney samples, the abundance of hepatic mARC1 and mARC2 was comparable, whereas mARC2 abundance in the kidney was approximately 9-fold higher in comparison with mARC1. The analysis of the third set of samples (i.e., S9 fraction) revealed that mARC1 abundance in the kidney, intestine, and lung was 5- to 13-fold lower than the liver S9 abundance, whereas mARC2 abundance was approximately 3- and 16-fold lower in the intestine and lung than the liver S9, respectively. In contrast, the kidney mARC2 abundance in the S9 fraction was approximately 2.5-fold higher as compared with the hepatic mARC2 abundance. The abundance of mARC enzymes in the heart was below the limit of quantification (∼0.6 pmol/mg protein). The mARC enzyme abundance data presented here can be used to develop physiologically based pharmacokinetic models for the prediction of in vivo pharmacokinetics of mARC substrates. SIGNIFICANCE STATEMENT: A precise targeted quantitative proteomics method was developed and applied to quantify newly discovered drug-metabolizing enzymes, mARC1 and mARC2, in pediatric and adult tissue samples. The data suggest that mARC enzymes are ubiquitously expressed in an isoform-specific manner in the human liver, kidney, intestine, and lung, and the enzyme abundance is not associated with age and sex. These data are important for developing physiologically based pharmacokinetic models for the prediction of in vivo pharmacokinetics of mARC substrates.

Recent advances in the ontogeny of drug disposition.

Developmental changes that occur throughout childhood have long been known to impact drug disposition. However, pharmacokinetic studies in the paediatric population have historically been limited due to ethical concerns arising from incorporating children into clinical trials. As such, much of the early work in the field of developmental pharmacology was reliant on difficult-to-interpret in vitro and in vivo animal studies. Over the last 2 decades, our understanding of the mechanistic processes underlying age-related changes in drug disposition has advanced considerably. Progress has largely been driven by technological advances in mass spectrometry-based methods for quantifying proteins implicated in drug disposition, and in silico tools that leverage these data to predict age-related changes in pharmacokinetics. This review summarizes our current understanding of the impact of childhood development on drug disposition, particularly focusing on research of the past 20 years, but also highlighting select examples of earlier foundational research. Equally important to the studies reviewed herein are the areas that we cannot currently describe due to the lack of research evidence; these gaps provide a map of drug disposition pathways for which developmental trends still need to be characterized.

Functional Consequences of Pravastatin Isomerization on OATP1B1-Mediated Transport.

Pravastatin acid (PVA) can be isomerized to its inactive metabolite 3'α-iso-pravastatin acid (3αPVA) under acidic pH conditions. Previous studies reported interindividual differences in circulating concentrations of PVA and 3αPVA. This study investigated the functional consequences of PVA isomerization on OATP1B1-mediated transport. We characterized 3αPVA inhibition of OATP1B1-mediated PVA uptake into human embryonic kidney 293 cells expressing the four different OATP1B1 proteins (*1a, *1b, *5, and *15). 3αPVA inhibited OATP1B1-mediated PVA uptake in all four OATP1B1 gene products but with lower IC50/Ki values for OATP1B1*5 and *15 than for the reference proteins (*1a and *1b). PVA and 3αPVA were transported by all four OATP1B1 proteins. Kinetic experiments revealed that maximal transport rates (Vmax values) for OATP1B1 variants *5 and *15 were lower than for *1a and *1b for both substrates. Apparent affinities for 3αPVA transport were similar for all four variants. However, the apparent affinity of OATP1B1*5 for 3αPVA was higher (lower Km value) than for PVA. These data confirm that PVA conversion to 3αPVA can have functional consequences on PVA uptake and impacts OATP1B1 variants more than the reference protein, thus highlighting another source variation that must be taken into consideration when optimizing the PVA dose-exposure relationship for patients. SIGNIFICANCE STATEMENT: 3'α-iso-pravastatin acid inhibits pravastatin uptake for all OATP1B1 protein types; however, the IC50 values were significantly lower in OATP1B1*5 and *15 transfected cells. This suggests that a lower concentration of 3'α-iso-pravastatin is needed to disrupt OATP1B1-mediated pravastatin uptake, secondary to decreased cell surface expression of functional OATP1B1 in variant-expressing cells. These data will refine previous pharmacokinetic models that are utilized to characterize pravastatin interindividual variability with an ultimate goal of maximizing efficacy at the lowest possible risk for toxicity.

The Respective Roles of CYP3A4 and CYP2D6 in the Metabolism of Pimozide to Established and Novel Metabolites.

Pimozide is a dopamine receptor antagonist indicated for the treatment of Tourette syndrome. Prior in vitro studies characterized N-dealkylation of pimozide to 1,3-dihydro-1-(4-piperidinyl)-2H-benzimidazol-2-one (DHPBI) via CYP3A4 and, to a lesser extent, CYP1A2 as the only notable routes of pimozide biotransformation. However, drug-drug interactions between pimozide and CYP2D6 inhibitors and CYP2D6 genotype-dependent effects have since been observed. To reconcile these incongruities between the prior in vitro and in vivo studies, we characterized two novel pimozide metabolites: 5-hydroxypimozide and 6-hydroxypimozide. Notably, 5-hydroxypimozide was the major metabolite produced by recombinant CYP2D6 (Km ∼82 nM, V max ∼0.78 pmol/min per picomoles), and DHPBI was the major metabolite produced by recombinant CYP3A4 (apparent Km ∼1300 nM, V max ∼2.6 pmol/min per picomoles). Kinetics in pooled human liver microsomes (HLMs) for the 5-hydroxylation (Km ∼2200 nM, V max ∼59 pmol/min per milligram) and N-dealkylation (Km ∼3900 nM, V max ∼600 pmol/min per milligram) reactions were also determined. Collectively, formation of DHPBI, 5-hydroxypimozide, and 6-hydroxypimozide accounted for 90% of pimozide depleted in incubations of NADPH-supplemented pooled HLMs. Studies conducted in HLMs isolated from individual donors with specific cytochrome P450 isoform protein abundances determined via mass spectrometry revealed that 5-hydroxypimozide (r 2 = 0.94) and 6-hydroxypimozide (r 2 = 0.86) formation rates were correlated with CYP2D6 abundance, whereas the DHPBI formation rate (r 2 = 0.98) was correlated with CYP3A4 abundance. Furthermore, the HLMs differed with respect to their capacity to form 5-hydroxypimozide relative to DHPBI. Collectively, these data confirm a role for CYP2D6 in pimozide clearance via 5-hydroxylation and provide an explanation for a lack of involvement when only DHPBI formation was monitored in prior in vitro studies. SIGNIFICANCE STATEMENT: Current CYP2D6 genotype-guided dosing information in the pimozide label is discordant with available knowledge regarding the primary biotransformation pathways. Herein, we characterize the CYP2D6-dependent biotransformation of pimozide to previously unidentified metabolites. In human liver microsomes, formation rates for the novel metabolites and a previously identified metabolite were determined to be a function of CYP2D6 and CYP3A4 content, respectively. These findings provide a mechanistic basis for observations of CYP2D6 genotype-dependent pimozide clearance in vivo.

Novel Protective Role of Nicotinamide Phosphoribosyltransferase in Acetaminophen-Induced Acute Liver Injury in Mice.

Acetaminophen overdose is the most common cause of acute liver injury (ALI) or acute liver failure in the United States. Its pathogenetic mechanisms are incompletely understood. Additional studies are warranted to identify new genetic risk factors for more mechanistic insights and new therapeutic target discoveries. The objective of this study was to explore the role and mechanisms of nicotinamide phosphoribosyltransferase (NAMPT) in acetaminophen-induced ALI. C57BL/6 Nampt gene wild-type (Nampt+/+), heterozygous knockout (Nampt+/-), and overexpression (NamptOE) mice were treated with overdose of acetaminophen, followed by histologic, biochemical, and transcriptomic evaluation of liver injury. The mechanism of Nampt in acetaminophen-induced hepatocytic toxicity was also explored in cultured primary hepatocytes. Three lines of evidence have convergently demonstrated that acetaminophen overdose triggers the most severe oxidative stress and necrosis, and the highest expression of key necrosis driving genes in Nampt+/- mice, whereas the effects in NamptOE mice were least severe relative to Nampt+/+ mice. Treatment of P7C3-A20, a small chemical molecule up-regulator of Nampt, ameliorated acetaminophen-induced mouse ALI over the reagent control. These findings support the fact that NAMPT protects against acetaminophen-induced ALI.

Ontogeny of Hepatic Sulfotransferases and Prediction of Age-Dependent Fractional Contribution of Sulfation in Acetaminophen Metabolism.

Cytosolic sulfotransferases (SULTs), including SULT1A, SULT1B, SULT1E, and SULT2A isoforms, play noteworthy roles in xenobiotic and endobiotic metabolism. We quantified the protein abundances of SULT1A1, SULT1A3, SULT1B1, and SULT2A1 in human liver cytosol samples (n = 194) by liquid chromatography-tandem mass spectrometry proteomics. The data were analyzed for their associations by age, sex, genotype, and ethnicity of the donors. SULT1A1, SULT1B1, and SULT2A1 showed significant age-dependent protein abundance, whereas SULT1A3 was invariable across 0-70 years. The respective mean abundances of SULT1A1, SULT1B1, and SULT2A1 in neonatal samples was 24%, 19%, and 38% of the adult levels. Interestingly, unlike UDP-glucuronosyltransferases and cytochrome P450 enzymes, SULT1A1 and SULT2A1 showed the highest abundance during early childhood (1 to <6 years), which gradually decreased by approx. 40% in adolescents and adults. SULT1A3 and SULT1B1 abundances were significantly lower in African Americans compared with Caucasians. Multiple linear regression analysis further confirmed the association of SULT abundances by age, ethnicity, and genotype. To demonstrate clinical application of the characteristic SULT ontogeny profiles, we developed and validated a proteomics-informed physiologically based pharmacokinetic model of acetaminophen. The latter confirmed the higher fractional contribution of sulfation over glucuronidation in the metabolism of acetaminophen in children. The study thus highlights that the ontogeny-based age-dependent fractional contribution (fm) of individual drug-metabolizing enzymes has better potential in prediction of drug-drug interactions and the effect of genetic polymorphisms in the pediatric population.

Evaluating metronidazole as a novel, safe CYP2A6 phenotyping probe in healthy adults.

AIMS: CYP2A6 is a genetically polymorphic enzyme resulting in differential substrate metabolism and health behaviours. Current phenotyping probes for CYP2A6 exhibit limitations related to procurement (deuterated cotinine), toxicity (coumarin), specificity (caffeine) and age-appropriate administration (nicotine, NIC). In vitro, CYP2A6 selectively forms 2-hydroxymetronidazole (2HM) from metronidazole (MTZ). The purpose of this study was to evaluate MTZ as a CYP2A6 phenotyping probe drug in healthy adults against the well-established method of measuring trans-3-hydroxycotinine (3HC)/cotinine (COT). METHODS: A randomized, cross-over, pharmacokinetic study was completed in 16 healthy, nonsmoking adults. Separated by a washout period of at least 2 weeks, MTZ 500 mg and NIC gum 2 mg were administered and plasma was sampled over 48 hours and 8 hours, respectively. Correlations of plasma metabolite/parent ratios (2HM/MTZ; 3HC/COT) were assessed by Pearson coefficient. CYP2A6 genotyping was conducted and incorporated as a variable of plasma ratio response. RESULTS: Correlations between the plasma ratio 2HM/MTZ and 3HC/COT were ≥ 0.9 at multiple time points (P < 0.001), demonstrating a wide window during which 2HM/MTZ can be queried post-MTZ dose. CYP2A6 genotype had significant impacts on both MTZ and NIC phenotyping ratios with decreased activity predicted phenotypes demonstrating 2HM/MTZ ratios ≤58% and 3HC/COT ratios ≤56% compared with extensive activity predicted phenotypes at all time points examined in the study (P < 0.05). No adverse events were reported in the MTZ arm while 38% (n = 6) of participants reported mild adverse events in the NIC arm. CONCLUSIONS: Metronidazole via 2HM/MTZ performed well as a novel, safe phenotyping probe for CYP2A6 in healthy adults.

A population pharmacokinetic model for simvastatin and its metabolites in children and adolescents.

PURPOSE: Poor adherence to dietary/behaviour modifications as interventions for hypercholesterolemia in paediatric patients often necessitates the initiation of statin therapy. The aim of this study was to develop a joint population pharmacokinetic model for simvastatin and four metabolites in children and adolescents to investigate sources of variability in simvastatin acid exposure in this patient population, in addition to SLCO1B1 genotype status. METHODS: Plasma concentrations of simvastatin and its four metabolites, demographic and polymorphism data for OATP1B1 and CYP3A5 were analysed utilising a population pharmacokinetic modelling approach from an existing single oral dose (10 mg < 17 years and 20 mg ≥ 18 years) pharmacokinetic dataset of 32 children and adolescents. RESULTS: The population PK model included a one compartment disposition model for simvastatin with irregular oral absorption described by two parallel absorption processes each consisting of sequential zero and first-order processes. The data for each metabolite were described by a one-compartment disposition model with the formation and elimination apparent parameters estimated. The model confirmed the statistically significant effect of c.521T>C (rs4149056) on the pharmacokinetics of the active metabolite simvastatin acid in children/adolescents, consistent with adult data. In addition, age was identified as a covariate affecting elimination clearances of 6-hydroxymethyl simvastatin acid and 3, 5 dihydrodiol simvastatin metabolites. CONCLUSION: The model developed describes the pharmacokinetics of simvastatin and its metabolites in children/adolescents capturing the effects of both c.521T>C and age on variability in exposure in this patient population. This joint simvastatin metabolite model is envisaged to facilitate optimisation of simvastatin dosing in children/adolescents.

Ontogeny-related pharmacogene changes in the pediatric liver transcriptome.

OBJECTIVES: The majority of drug dosing studies are based on adult populations, with modification of the dosing for children based on size and weight. This rudimentary approach for drug dosing children is limited, as biologically a child can differ from an adult in far more aspects than just size and weight. Specifically, understanding the ontogeny of childhood liver development is critical in dosing drugs that are metabolized through the liver, as the rate of metabolism determines the duration and intensity of a drug's pharmacologic action. Therefore, we set out to determine pharmacogenes that change over childhood development, followed by a secondary agnostic analysis, assessing changes transcriptome wide. MATERIALS AND METHODS: A total of 47 human liver tissue samples, with between 10 and 13 samples in four age groups spanning childhood development, underwent pair-end sequencing. Kruskal-Wallis and Spearman's rank correlation tests were used to determine the association of gene expression levels with age. Gene set analysis based on the pathways in KEGG utilized the gamma method. Correction for multiple testing was completed using q-values. RESULTS: We found evidence for increased expression of 'very important pharmacogenes', for example, coagulation factor V (F5) (P=6.7×10(-7)), angiotensin I converting enzyme (ACE) (P=6.4×10(-3)), and solute carrier family 22 member 1 (SLC22A1) (P=7.0×10(-5)) over childhood development. In contrast, we observed a significant decrease in expression of two alternative CYP3A7 transcripts (P=1.5×10(-5) and 3.0×10(-5)) over development. The analysis of genome-wide changes detected transcripts in the following genes with significant changes in mRNA expression (P<1×10(-9) with false discovery rate<5×0(-5)): ADCY1, PTPRD, CNDP1, DCAF12L1 and HIP1. Gene set analysis determined ontogeny-related transcriptomic changes in the renin-angiotensin pathway (P<0.002), with lower expression of the pathway, in general, observed in liver samples from younger participants. CONCLUSION: Considering that the renin-angiotensin pathway plays a central role in blood pressure and plasma sodium concentration, and our observation that ACE and PTPRD expression increased over the spectrum of childhood development, this finding could potentially impact the dosing of an entire class of drugs known as ACE-inhibitors in pediatric patients.

Nicotinamide Phosphoribosyltransferase Deficiency Potentiates the Antiproliferative Activity of Methotrexate through Enhanced Depletion of Intracellular ATP.

Lower plasma nicotinamide phosphoribosyltransferase (NAMPT) levels are associated with improved response to methotrexate (MTX) in patients with juvenile idiopathic arthritis. Cell-based studies confirmed that reduced cellular NAMPT activity potentiates the pharmacologic activity of MTX; however, the mechanism of this interaction has yet to be defined. Therefore, in this study, we investigate the mechanism of enhanced pharmacologic activity of MTX in NAMPT-deficient A549 cells. Small interfering RNA-based silencing of NAMPT expression resulted in a greater than 3-fold increase in sensitivity to MTX (P < 0.005) that was completely reversed by supplementation with folinic acid. Despite a 68% reduction in cellular NAD levels in NAMPT-deficient cells, no change in expression or activity of dihydrofolate reductase was observed and uptake of MTX was not significantly altered. MTX did not potentiate the depletion of cellular NAD levels, but NAMPT-deficient cells had significant elevations in levels of intermediates of de novo purine biosynthesis and were 4-fold more sensitive to depletion of ATP by MTX (P < 0.005). Supplementation with hypoxanthine and thymidine completely reversed the antiproliferative activity of MTX in NAMPT-deficient cells and corresponded to repletion of the cellular ATP pool without any effect on NAD levels. Together, these findings demonstrate that increased MTX activity with decreased NAMPT expression is dependent on the antifolate activity of MTX and is driven by enhanced sensitivity to the ATP-depleting effects of MTX. For the first time, these findings provide mechanistic details to explain the increase in pharmacological activity of MTX under conditions of reduced NAMPT activity.

Hepatic Abundance and Activity of Androgen- and Drug-Metabolizing Enzyme UGT2B17 Are Associated with Genotype, Age, and Sex.

The major objective of this study was to investigate the association of genetic and nongenetic factors with variability in protein abundance and in vitro activity of the androgen-metabolizing enzyme UGT2B17 in human liver microsomes (n = 455). UGT2B17 abundance was quantified by liquid chromatography-tandem mass spectrometry proteomics, and enzyme activity was determined by using testosterone and dihydrotestosterone as in vitro probe substrates. Genotyping or gene resequencing and mRNA expression were also evaluated. Multivariate analysis was used to test the association of UGT2B17 copy number variation, single nucleotide polymorphisms (SNPs), age, and sex with its mRNA expression, abundance, and activity. UGT2B17 gene copy number and SNPs (rs7436962, rs9996186, rs28374627, and rs4860305) were associated with gene expression, protein levels, and androgen glucuronidation rates in a gene dose-dependent manner. UGT2B17 protein (mean ± S.D. picomoles per milligram of microsomal protein) is sparsely expressed in children younger than 9 years (0.12 ± 0.24 years) but profoundly increases from age 9 years to adults (∼10-fold) with ∼2.6-fold greater abundance in males than in females (1.2 vs. 0.47). Association of androgen glucuronidation with UGT2B15 abundance was observed only in the low UGT2B17 expressers. These data can be used to predict variability in the metabolism of UGT2B17 substrates. Drug companies should include UGT2B17 in early phenotyping assays during drug discovery to avoid late clinical failures.

Genetic and Nongenetic Factors Associated with Protein Abundance of Flavin-Containing Monooxygenase 3 in Human Liver.

Hepatic flavin-containing mono-oxygenase 3 (FMO3) metabolizes a broad array of nucleophilic heteroatom (e.g., N or S)-containing xenobiotics (e.g., amphetamine, sulindac, benzydamine, ranitidine, tamoxifen, nicotine, and ethionamide), as well as endogenous compounds (e.g., catecholamine and trimethylamine). To predict the effect of genetic and nongenetic factors on the hepatic metabolism of FMO3 substrates, we quantified FMO3 protein abundance in human liver microsomes (HLMs; n = 445) by liquid chromatography-tandem mass chromatography proteomics. Genotyping/gene resequencing, mRNA expression, and functional activity (with benzydamine as probe substrate) of FMO3 were also evaluated. FMO3 abundance increased 2.2-fold (13.0 ± 11.4 pmol/mg protein vs. 28.0 ± 11.8 pmol/mg protein) from neonates to adults. After 6 years of age, no significant difference in FMO3 abundance was found between children and adults. Female donors exhibited modestly higher mRNA fragments per kilobase per million reads values (139.9 ± 76.9 vs. 105.1 ± 73.1; P < 0.001) and protein FMO3 abundance (26.7 ± 12.0 pmol/mg protein vs. 24.1 ± 12.1 pmol/mg protein; P < 0.05) compared with males. Six single nucleotide polymorphisms (SNPs), including rs2064074, rs28363536, rs2266782 (E158K), rs909530 (N285N), rs2266780 (E308G), and rs909531, were associated with significantly decreased protein abundance. FMO3 abundance in individuals homozygous and heterozygous for haplotype 3 (H3), representing variant alleles for all these SNPs (except rs2066534), were 50.8% (P < 0.001) and 79.5% (P < 0.01), respectively, of those with the reference homozygous haplotype (H1, representing wild-type). In summary, FMO3 protein abundance is significantly associated with age, gender, and genotype. These data are important in predicting FMO3-mediated heteroatom-oxidation of xenobiotics and endogenous biomolecules in the human liver.

Prediction of CYP2D6 phenotype from genotype across world populations.

PURPOSE: Owing to its highly polymorphic nature and major contribution to the metabolism and bioactivation of numerous clinically used drugs, CYP2D6 is one of the most extensively studied drug-metabolizing enzymes and pharmacogenes. CYP2D6 alleles confer no, decreased, normal, or increased activity and cause a wide range of activity among individuals and between populations. However, there is no standard approach to translate diplotypes into predicted phenotype. METHODS: We exploited CYP2D6 allele-frequency data that have been compiled for Clinical Pharmacogenetics Implementation Consortium (CPIC) guidelines (>60,000 subjects, 173 reports) in order to estimate genotype-predicted phenotype status across major world populations based on activity score (AS) assignments. RESULTS: Allele frequencies vary considerably across the major ethnic groups predicting poor metabolizer status (AS = 0) between 0.4 and 5.4% across world populations. The prevalence of genotypic intermediate (AS = 0.5) and normal (AS = 1, 1.5, or 2) metabolizers ranges between 0.4 and 11% and between 67 and 90%, respectively. Finally, 1 to 21% of subjects (AS >2) are predicted to have ultrarapid metabolizer status. CONCLUSIONS: This comprehensive study summarizes allele frequencies, diplotypes, and predicted phenotype across major populations, providing a rich data resource for clinicians and researchers. Challenges of phenotype prediction from genotype data are highlighted and discussed.Genet Med 19 1, 69-76.

Age-dependent Protein Abundance of Cytosolic Alcohol and Aldehyde Dehydrogenases in Human Liver.

Hepatic cytosolic alcohol and aldehyde dehydrogenases (ADHs and ALDHs) catalyze the biotransformation of xenobiotics (e.g., cyclophosphamide and ethanol) and vitamin A. Because age-dependent hepatic abundance of these proteins is unknown, we quantified protein expression of ADHs and ALDH1A1 in a large cohort of pediatric and adult human livers by liquid chromatography coupled with tandem mass spectrometry proteomics. Purified proteins were used as calibrators. Two to three surrogate peptides per protein were quantified in trypsin digests of liver cytosolic samples and calibrator proteins under optimal conditions of reproducibility. Neonatal levels of ADH1A, ADH1B, ADH1C, and ALDH1A1 were 3-, 8-, 146-, and 3-fold lower than the adult levels, respectively. For all proteins, the abundance steeply increased during the first year of life, which mostly reached adult levels during early childhood (age between 1 and 6 years). Only for ADH1A protein abundance in adults (age > 18 year) was ∼40% lower relative to the early childhood group. Abundances of ADHs and ALDH1A1 were not associated with sex in samples with age > 1 year compared with males. Known single nucleotide polymorphisms had no effect on the protein levels of these proteins. Quantification of ADHs and ALDH1A1 protein levels could be useful in predicting disposition and response of substrates of these enzymes in younger children.

Age-Dependent Absolute Abundance of Hepatic Carboxylesterases (CES1 and CES2) by LC-MS/MS Proteomics: Application to PBPK Modeling of Oseltamivir In Vivo Pharmacokinetics in Infants.

The age-dependent absolute protein abundance of carboxylesterase (CES) 1 and CES2 in human liver was investigated and applied to predict infant pharmacokinetics (PK) of oseltamivir. The CES absolute protein abundance was determined by liquid chromatography-tandem mass spectrometry proteomics in human liver microsomal and cytosolic fractions prepared from tissue samples obtained from 136 pediatric donors and 35 adult donors. Two surrogate peptides per protein were selected for the quantification of CES1 and CES2 protein abundance. Purified CES1 and CES2 protein standards were used as calibrators, and the heavy labeled peptides were used as the internal standards. In hepatic microsomes, CES1 and CES2 abundance (in picomoles per milligram total protein) increased approximately 5-fold (315.2 vs. 1664.4) and approximately 3-fold (59.8 vs. 174.1) from neonates to adults, respectively. CES1 protein abundance in liver cytosol also showed age-dependent maturation. Oseltamivir carboxylase activity was correlated with protein abundance in pediatric and adult liver microsomes. The protein abundance data were then used to model in vivo PK of oseltamivir in infants using pediatric physiologically based PK modeling and incorporating the protein abundance-based ontogeny function into the existing pediatric Simcyp model. The predicted pediatric area under the curve, maximal plasma concentration, and time for maximal plasma concentration values were below 2.1-fold of the clinically observed values, respectively.