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).

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.

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.

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.

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.

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.

Characterization of Atomoxetine Biotransformation and Implications for Development of PBPK Models for Dose Individualization in Children.

Atomoxetine (ATX) is a second-line nonstimulant medication used to control symptoms of attention deficit hyperactivity disorder (ADHD). Inconsistent therapeutic efficacy has been reported with ATX, which may be related to variable CYP2D6-mediated drug clearance. We characterized ATX metabolism in a panel of human liver samples as a basis for a bottom-up PBPK model to aid in ATX exposure prediction and control. Km, Vmax, and Clint values in pooled human liver microsomes (HLMs) were 2.4 µM, 479 pmol/min/mg protein, and 202 µl/min/mg protein, respectively. Mean population values of kinetic parameters are not adequate to describe variability in a population, given that Km, Vmax, and Clint values from single-donor HLMs ranged from 0.93 to 79.2 µM, 20.0 to 1600 pmol/min/mg protein, and 0.3 to 936 µl/min/mg protein. All kinetic parameters were calculated from 4-hydroxyatomoxetine (4-OH-ATX) formation. CYP2E1 and CYP3A contributed to 4-OH-ATX formation in livers with CYP2D6 intermediate and poor metabolizer status. In HLMs with lower CYP2D6 activity levels, 2-hydroxymethylatomoxetine (2-CH2OH-ATX) formation became a more predominant pathway of metabolism, which appeared to be catalyzed by CYP2B6. ATX biotransformation at clinically relevant plasma concentrations was characterized in a panel of pediatric HLM (n = 116) samples by evaluating primary metabolites. Competing pathways of ATX metabolism [N-desmethylatomoxetine (NDM-ATX) and 2-CH2OH-ATX formation] had increasing importance in livers with lesser CYP2D6 activity, but, overall ATX clearance was still compromised. Modeling ATX exposure to individualize therapy would require comprehensive knowledge of factors that affect CYP2D6 activity as well as an understanding of competing pathways, particularly for individuals with lower CYP2D6 activity.

Developmental Expression of CYP2B6: A Comprehensive Analysis of mRNA Expression, Protein Content and Bupropion Hydroxylase Activity and the Impact of Genetic Variation.

Although CYP2B6 catalyzes the biotransformation of many drugs used clinically for children and adults, information regarding the effects of development on CYP2B6 expression and activity are scarce. Utilizing a large panel of human liver samples (201 donors: 24 fetal, 141 pediatric, and 36 adult), we quantified CYP2B6 mRNA and protein expression levels, characterized CYP2B6 (bupropion hydroxylase) activity in human liver microsomes (HLMs), and performed an extensive genotype analysis to differentiate CYP2B6 haplotypes such that the impact of genetic variation on these parameters could be assessed. Fetal livers contained extremely low levels of CYP2B6 mRNA relative to postnatal samples and fetal HLMs did not appear to catalyze bupropion hydroxylation; however, fetal CYP2B6 protein levels were not significantly different from postnatal levels. Considerable interindividual variation in CYP2B6 mRNA expression, protein levels, and activity was observed in postnatal HLMs (mRNA, ∼40,000-fold; protein, ∼300-fold; activity, ∼600-fold). The extremely wide range of interindividual variability in CYP2B6 expression and activity was significantly associated with age (P < 0.01) following log transformation of the data. Our data suggest that CYP2B6 activity appears as early as the first day of life, increases through infancy, and by 1 year of age, CYP2B6 levels and activity may approach those of adults. Surprisingly, CYP2B6 interindividual variability was not significantly associated with genetic variation in CYP2B6, nor was it associated with differences in gender or ethnicity, suggesting that factors other than these are largely responsible for the wide range of variability in CYP2B6 expression and activity observed among a large group of individuals/samples.

In Vitro Hepatic Oxidative Biotransformation of Trimethoprim.

Trimethoprim (TMP) has been widely used since the 1960s, both alone and in combination with sulfamethoxazole. Unfortunately, information regarding the role that cytochrome P450 enzymes (P450s) play in the formation of TMP primary metabolites is scarce. Hence, we undertook in vitro studies to identify and more fully characterize the P450s that catalyze formation of six TMP primary metabolites: TMP 1-N-oxide (1-NO-TMP) and 3-N-oxide (3-NO-TMP), 3'- and 4'-desmethyl-TMP, a benzylic alcohol (Cα-OH-TMP), and an N-acetyl cysteine (NAC) adduct of TMP (Cα-NAC-TMP). Formation kinetics for each TMP metabolite in human liver microsomes (HLMs) were consistent with single-enzyme Michaelis-Menten kinetics, and Km values were markedly above (≥10-fold) the therapeutic concentrations of TMP (50 µM). The combined results from correlation studies between rates of metabolite formation and marker P450 activities in a panel of HLMs along with inhibition studies utilizing selective P450 inhibitors incubated with pooled HLMs suggested that 1-NO-TMP, Cα-NAC-TMP, and Cα-OH-TMP were predominantly formed by CYP3A4. In contrast, 3-NO-TMP was formed predominantly by CYP1A2 in HLMs and inhibited by α-naphthoflavone. 4'-Desmethyl-TMP, which is believed to be a reactive TMP metabolite precursor, was formed by several P450s, including CYP3A4, correlated with multiple P450 activities, but was inhibited primarily by ketoconazole (up to 50%), suggesting that CYP3A4 makes a major contribution to TMP 4'-demethylation. TMP 3'-demethylation was catalyzed by multiple P450s, including CYP2C9, correlated with CYP2C9 activity, and was inhibited by sulfaphenazole (up to 40%). Overall, CYP2C9 and CYP3A4 appear to be the most significant contributors to TMP primary metabolism.

Variability in Expression of CYP3A5 in Human Fetal Liver.

Members of the cytochrome P450 3A (CYP3A) subfamily of drug metabolizing enzymes exhibit developmental changes in expression in human liver characterized by a transition between CYP3A7 and CYP3A4 over the first few years of life. In contrast, the developmental expression of CYP3A5 is less well understood due to polymorphic expression of the enzyme in human tissues as a result of the prevalence of the CYP3A5*3 allele, which leads to alternative splicing. We further explored the expression of CYP3A5 and the impact of alternative splicing on the variability of CYP3A5 functional activity in a large bank of human prenatal liver samples (7 to 32 weeks of age postconception). The expression of normally spliced CYP3A5 mRNA in all human fetal liver samples varied 235-fold whereas CYP3A5 SV1 mRNA was only detected in fetal liver samples with at least one CYP3A5*3 allele. Formation of 1'-OH midazolam (MDZ) varied 79-fold, and the ratio of 1'-OH MDZ to 4-OH MDZ varied 8-fold and depended on the presence or absence of the CYP3A5*3 allele. Formation of 4-OH MDZ was significantly associated with 1'-OH MDZ (r(2) = 0.76, P < 0.0001) but varied (36-fold) independently of CYP3A5 genotype or expression. The substantial interindividual variability that remains even after stratification for CYP3A5 genotype suggests that factors such as environmental exposure and epigenetic alterations act in addition to genetic variation to contribute to the variability of CYP3A5 expression in human prenatal liver.

Urinary biomarkers of trimethoprim bioactivation in vivo following therapeutic dosing in children.

The antimicrobial trimethoprim-sulfamethoxazole (TMP-SMX) is widely used for the treatment of skin and soft-tissue infections in the outpatient setting. Despite its therapeutic benefits, TMP-SMX has been associated with a number of adverse drug reactions, which have been primarily attributed to the formation of reactive metabolites from SMX. Recently, in vitro experiments have demonstrated that TMP may form reactive intermediates as well. However, evidence of TMP bioactivation in patients has not yet been demonstrated. In this study, we performed in vitro trapping experiments with N-acetyl-l-cysteine (NAC) to determine stable markers of reactive TMP intermediates, focusing on eight potential markers (NAC-TMP adducts), some of which were previously identified in vitro. We developed a specific and sensitive assay involving liquid chromatography followed by tandem mass spectrometry for measurement of these adducts in human liver microsomal samples and expanded the methodology toward the detection of these analytes in human urine. Urine samples from four patients receiving TMP-SMX treatment were analyzed, and all samples demonstrated the presence of six NAC-TMP adducts, which were also detected in vitro. These adducts are consistent with the formation of imino-quinone-methide and para-quinone-methide reactive intermediates in vivo. As a result, the TMP component of TMP-SMX should be considered as well when evaluating adverse drug reactions to TMP-SMX.

Age-Dependent Abundance of CYP450 Enzymes Involved in Metronidazole Metabolism: Application to Pediatric PBPK Modeling.

The expression of cytochrome P450 (CYP) enzymes is highly variable and associated with factors, such as age, genotype, sex, and disease states. In this study, quantification of metronidazole metabolizing CYP isoforms (CYP2A6, CYP2E1, CYP3A4, CYP3A5, and CYP3A7) in human liver microsomes from 115 children and 35 adults was performed using a quantitative proteomics method. The data confirmed age-dependent increase in CYP2A6, CYP2E1, and CYP3A4 abundance, whereas, as expected, CYP3A7 abundance showed postnatal decrease with age. In particular, the fold difference (neonatal to adulthood levels) in the protein abundance of CYP2A6, CYP2E1, and CYP3A4 was 14, 11, and 20, respectively. In contrast, protein abundance of CYP3A7 was > 125-fold higher in the liver microsomes of neonates than of adults. The abundance of CYP2A6 and CYP3A5 was associated with genotypes, rs4803381 and rs776746, respectively. A proteomics-informed physiologically based pharmacokinetic (PBPK) model was developed to describe the pharmacokinetics of metronidazole and its primary metabolite, 2-hydroxymethylmetronidazole. The model revealed an increase in the metabolite-to-parent ratio with age and showed a strong correlation between CYP2A6 abundance and metabolite formation (r2 = 0.75). Notably, the estimated contribution of CYP3A7 was ~ 75% in metronidazole clearance in neonates. These data suggest that variability in CYP2A6 and CYP3A7 in younger children poses the risk of variable pharmacokinetics of metronidazole and its active metabolite with a potential impact on drug efficacy and safety. No sex-dependent difference was observed in the protein abundance of the studied CYPs. The successful integration of hepatic CYP ontogeny data derived from a large liver bank into the pediatric PBPK model of metronidazole can be extended to other drugs metabolized by the studied CYPs.

Influence of novel CYP2C-haplotype on proton pump inhibitor pharmacokinetics in children.

In this brief report, we provide an analysis of the influence of a novel CYP2C haplotype (CYP2C:TG) on proton pump inhibitor (PPI) pharmacokinetics (PK) in children. The CYP2C:TG haplotype has been proposed to be associated with increased CYP2C19 activity. We sought to determine if this CYP2C:TG haplotype resulted in similar alterations in metabolism for proton pump inhibitors, which are primarily metabolized by CYP2C19. In a cohort of 41 children aged 6-21 participating in a PPI pharmacokinetic study, effects of the CYP2C:TG allele were assessed by fitting two linear regression models for each of the six PK outcomes assessed, the second of which accounted for the presence of the CYP2C:TG allele. The difference in R2 values between the two models was computed to quantify the variability in the outcome that could be accounted for by the CYP2C:TG allele after adjustment for the CYP2C19 genotype. We found the CYP2C:TG haplotype to have no measurable additive impact on CYP2C19-mediated metabolism of PPIs in vivo in older children and adolescents. The findings of this study do not support the clinical utility of routine testing for the CYP2C:TG haplotype to guide PPI dose adjustments in children.

The role of human cytochrome P450 enzymes in the formation of 2-hydroxymetronidazole: CYP2A6 is the high affinity (low Km) catalyst.

Despite metronidazole's widespread clinical use since the 1960s, the specific enzymes involved in its biotransformation have not been previously identified. Hence, in vitro studies were conducted to identify and characterize the cytochrome P450 enzymes involved in the formation of the major metabolite, 2-hydroxymetronidazole. Formation of 2-hydroxymetronidazole in human liver microsomes was consistent with biphasic, Michaelis-Menten kinetics. Although several cDNA-expressed P450 enzymes catalyzed 2-hydroxymetronidazole formation at a supratherapeutic concentration of metronidazole (2000 µM), at a "therapeutic concentration" of 100 µM only CYPs 2A6, 3A4, 3A5, and 3A7 catalyzed metronidazole 2-hydroxylation at rates substantially greater than control vector, and CYP2A6 catalyzed 2-hydroxymetronidazole formation at rates 6-fold higher than the next most active enzyme. Kinetic studies with these recombinant enzymes revealed that CYP2A6 has a Km = 289 µM which is comparable to the Km for the high-affinity (low-Km) enzyme in human liver microsomes, whereas the Km values for the CYP3A enzymes corresponded with the low-affinity (high-Km) component. The sample-to-sample variation in 2-hydroxymetronidazole formation correlated significantly with CYP2A6 activity (r ≥ 0.970, P < 0.001) at substrate concentrations of 100 and 300 µM. Selective chemical inhibitors of CYP2A6 inhibited metronidazole 2-hydroxylation in a concentration-dependent manner and inhibitory antibodies against CYP2A6 virtually eliminated metronidazole 2-hydroxylation (>99%). Chemical and antibody inhibitors of other P450 enzymes had little or no effect on metronidazole 2-hydroxylation. These results suggest that CYP2A6 is the primary catalyst responsible for the 2-hydroxylation of metronidazole, a reaction that may function as a marker of CYP2A6 activity both in vitro and in vivo.

The Impact of the CYP2D6 "Enhancer" Single Nucleotide Polymorphism on CYP2D6 Activity.

rs5758550 has been associated with enhanced transcription and suggested to be a useful marker of CYP2D6 activity. As there are limited and inconsistent data regarding the utility of this distant "enhancer" single nucleotide polymorphism (SNP), our goal was to further assess the impact of rs5758550 on CYP2D6 activity toward two probe substrates, atomoxetine (ATX) and dextromethorphan (DM), using in vivo urinary metabolite (DM; n = 188) and pharmacokinetic (ATX; n = 70) and in vitro metabolite formation (ATX and DM; n = 166) data. All subjects and tissues were extensively genotyped, the "enhancer" SNP phased with established CYP2D6 haplotypes either computationally or experimentally, and the impact on CYP2D6 activity investigated using several linear models of varying complexity to determine the proportion of variability in CYP2D6 activity captured by each model. For all datasets and models, the "enhancer" SNP had no or only a modest impact on CYP2D6 activity prediction. An increased effect, when present, was more pronounced for ATX than DM suggesting potential substate-dependency. In addition, CYP2D6*2 alleles with the "enhancer" SNP were associated with modestly higher metabolite formation rates in vitro, but not in vivo; no effect was detected for CYP2D6*1 alleles with "enhancer" SNP. In summary, it remains inconclusive whether the small effects detected in this investigation are indeed caused by the "enhancer" SNP or are rather due to its incomplete linkage with other variants within the gene. Taken together, there does not appear to be sufficient evidence to warrant the "enhancer" SNP be included in clinical CYP2D6 pharmacogenetic testing.

A longitudinal study of cytochrome P450 2D6 (CYP2D6) activity during adolescence.

CYP2D6 substrates are among the most highly prescribed medications in teenagers and also commonly associated with serious adverse events. To investigate the relative contributions of genetic variation, growth, and development on CYP2D6 activity during puberty, healthy children and adolescents 7-15 years of age at enrollment participated in a longitudinal phenotyping study involving administration of 0.3 mg/kg dextromethorphan (DM) and 4-h urine collection every 6 months for 3 years (7 total visits). At each visit, height, weight, and sexual maturity were recorded, and CYP2D6 activity was determined as the urinary molar ratio of DM to its metabolite dextrorphan (DX). A total of 188 participants completed at least one visit, and 102 completed all seven study visits. Following univariate analysis, only CYP2D6 activity score (p < 0.001), urinary pH (p < 0.001), weight (p = 0.018), and attention-deficit/hyperactivity disorder (ADHD) diagnosis (p < 0.001) were significantly correlated with log(DM/DX). Results of linear mixed model analysis with random intercept, random slope covariance structure revealed that CYP2D6 activity score had the strongest effect on log(DM/DX), with model-estimated average log(DM/DX) being 3.8 SDs higher for poor metabolizers than for patients with activity score 3. A moderate effect on log(DM/DX) was observed for sex, and smaller effects were observed for ADHD diagnosis and urinary pH. The log(DM/DX) did not change meaningfully with age or pubertal development. CYP2D6 genotype remains the single, largest determinant of variability in CYP2D6 activity during puberty. Incorporation of genotype-based dosing guidelines should be considered for CYP2D6 substrates given the prevalent use of these agents in this pediatric age group.

Utility of the 13 C-pantoprazole breath test as a CYP2C19 phenotyping probe for children.

The 13 C-pantoprazole breath test (PAN-BT) is a safe, noninvasive, in vivo CYP2C19 phenotyping probe for adults. Our objective was to evaluate PAN-BT performance in children, with a focus on discriminating individuals who, according to guidelines from the Clinical Pharmacology Implementation Consortium (CPIC), would benefit from starting dose escalation versus reduction for proton pump inhibitors (PPIs). Children (n = 65, 6-17 years) genotyped for CYP2C19 variants *2, *3, *4, and *17 received a single oral dose of 13 C-pantoprazole. Plasma concentrations of pantoprazole and its metabolites, and changes in exhaled 13 CO2 (termed delta-over-baseline or DOB), were measured 10 times over 8 h using high performance liquid chromatography with ultraviolet detection and spectrophotometry, respectively. Pharmacokinetic parameters of interest were generated and DOB features derived using feature engineering for the first 180 min postadministration. DOB features, age, sex, and obesity status were used to run bootstrap analysis at each timepoint (Ti ) independently. For each iteration, stratified samples were drawn based on genotype prevalence in the original cohort. A random forest was trained, and predictive performance of PAN-BT was evaluated. Strong discriminating ability for CYP2C19 intermediate versus normal/rapid metabolizer phenotype was noted at DOBT30 min (mean sensitivity: 0.522, specificity: 0.784), with consistent model outperformance over a random or a stratified classifier approach at each timepoint (p < 0.001). With additional refinement and investigation, the test could become a useful and convenient dosing tool in clinic to help identify children who would benefit most from PPI dose escalation versus dose reduction, in accordance with CPIC guidelines.