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

Deprescribing Antidepressants in Children and Adolescents: A Systematic Review of Discontinuation Approaches, Cross-Titration, and Withdrawal Symptoms

Background: Antidepressant medications, including selective serotonin reuptake inhibitors (SSRIs) and serotonin-norepinephrine reuptake inhibitors (SNRIs), are commonly used to treat depressive, anxiety, and obsessive-compulsive disorders in youth. Yet, data on discontinuing these medications, withdrawal symptoms, and strategies to switch between them are limited. Methods: We searched PubMed and ClinicalTrials.gov through June 1, 2024, to identify randomized controlled trials assessing antidepressant discontinuation in youth. We summarized pediatric pharmacokinetic data to inform tapering and cross-titration strategies for antidepressants and synthesized these data with reports of antidepressant withdrawal. Results: Our search identified 528 published articles, of which 28 were included. In addition, 19 records were obtained through other methods, with 14 included. The corpus of records included 13 randomized, double-blind, placebo-controlled trials (3026 patients), including SSRIs (K = 10), SNRIs (K = 4), and TCAs (K = 1), ranging from 4 to 35 weeks. Deprescribing antidepressants requires considering clinical status, treatment response, and, in cross-titration cases, the pharmacokinetics and pharmacodynamics of both medications. Antidepressant withdrawal symptoms are related to the pharmacokinetics of the medication, which vary across antidepressants and may include irritability, palpitations, anxiety, nausea, sweating, headaches, insomnia, paresthesia, and dizziness. These symptoms putatively involve changes in serotonin transporter expression and receptor sensitivity, impacting the serotonin, dopamine, and norepinephrine pathways. Conclusions: Although approaches to deprescribing antidepressants in pediatric patients are frequently empirically guided, accumulating data related to the course of relapse and withdrawal symptoms, as well as the pharmacokinetic and pharmacodynamic properties of medications, should inform these approaches. Recommendations within this review support data-informed discussions of deprescribing-including when and how-that are critically important in the clinician-family-patient relationship.

Sub- and supratherapeutic efavirenz plasma concentrations with risk for HIV therapy failure are mainly genetically explained in Ugandan children: The prospective GENEFA cohort study

AIMS

Interindividual variations in efavirenz (EFV) plasma concentrations are extensive, but paediatric data on its consequences for viral control are scarce. The aim of this study was to explore the role of genetic variation in achieving therapeutic efavirenz plasma concentrations in a cohort of Ugandan children and the linkage between genetic CYP2B6 variants, EFV plasma variability, viral resistance and viral outcome.

METHODS

Ninety-nine treatment-naïve children, aged 3-12 years and living with HIV, were followed for 24 weeks after ART initiation assessing mid-dose efavirenz plasma concentrations, HIV RNA, HIV drug resistance and adherence. Polymorphisms in genes coding for drug-metabolizing enzymes were genotyped. Efavirenz concentrations were determined by liquid chromatography coupled with high-resolution tandem mass spectrometry. Metabolizer phenotype was predicted from composite genotypes of CYP2B6 (c.516G>T and c.983 T>C). A mixed effects restricted maximum likelihood regression model was used to identify important factors for efavirenz exposure.

RESULTS

Efavirenz plasma concentrations were below the therapeutic interval (1000-4000 mg/mL) in 12-17% and above in 21-24% of measurements. Eight children had persisting subtherapeutic concentrations, five of which failed virologically and three acquired at least one new resistant mutation. Multivariate modelling explained 70% of interindividual variation in plasma concentration, with treatment duration, adherence, CYP2B6c.136A>G, and metabolizer phenotype as independent predictors of EFV concentration. In univariate analysis, metabolizer phenotype explained 50% of interindividual variation.

CONCLUSIONS

Metabolizer phenotype explained 50% of interindividual variation in efavirenz plasma concentration. Autoinduction was not confirmed and >33% of the concentrations were outside the therapeutic interval. Subtherapeutic concentrations worsened virological resistance and outcomes. Genotype-based dosing may help avert both sub- and supratherapeutic efavirenz plasma concentrations in Ugandan children.

Clinical Pharmacogenetics Implementation Consortium Guideline for CYP2B6 Genotype and Methadone Therapy

Methadone is a mu (μ) opioid receptor agonist used clinically in adults and children to manage opioid use disorder, neonatal abstinence syndrome, and acute and chronic pain. It is typically marketed as a racemic mixture of R- and S-enantiomers. R-methadone has 30-to 50-fold higher analgesic potency than S-methadone, and S-methadone has a greater adverse effect (prolongation) on the cardiac QTc interval. Methadone undergoes stereoselective metabolism. CYP2B6 is the primary enzyme responsible for catalyzing the metabolism of both enantiomers to the inactive metabolites, S- and R-2-ethylidene-1,5-dimethyl-3,3-diphenylpyrrolidine (S- and R-EDDP). Genetic variation in the CYP2B6 gene has been investigated in the context of implications for methadone pharmacokinetics, dose, and clinical outcomes. Most CYP2B6 variants result in diminished or loss of CYP2B6 enzyme activity, which can lead to higher plasma methadone concentrations (affecting S- more than R-methadone). However, the data do not consistently indicate that CYP2B6-based metabolic variability has a clinically significant effect on methadone dose, efficacy, or QTc prolongation. Expert analysis of the published literature does not support a change from standard methadone prescribing based on CYP2B6 genotype (updates at www.cpicpgx.org).

Physiologically Based Pharmacokinetics Modeling in the Neonatal Population-Current Advances, Challenges, and Opportunities

Physiologically based pharmacokinetic (PBPK) modeling is an approach to predicting drug pharmacokinetics, using knowledge of the human physiology involved and drug physiochemical properties. This approach is useful when predicting drug pharmacokinetics in under-studied populations, such as pediatrics. PBPK modeling is a particularly important tool for dose optimization for the neonatal population, given that clinical trials rarely include this patient population. However, important knowledge gaps exist for neonates, resulting in uncertainty with the model predictions. This review aims to outline the sources of variability that should be considered with developing a neonatal PBPK model, the data that are currently available for the neonatal ontogeny, and lastly to highlight the data gaps where further research would be needed.

Association between CYP2B6 genetic variability and cyclophosphamide therapy in pediatric patients with neuroblastoma

Cyclophosphamide, an oxazaphosphorine prodrug is frequently used in treatment of neuroblastoma, which is one of the most prevalent solid organ malignancies in infants and young children. Cytochrome P450 2B6 (CYP2B6) is the major catalyst and CYP2C19 is the minor enzyme in bioactivation and inactivation pathways of cyclophosphamide. CYP-mediated metabolism may contribute to the variable pharmacokinetics of cyclophosphamide and its toxic byproducts leading to insufficient response to the therapy and development of clinically significant side effects. The aim of the study was to reveal the contribution of pharmacogenetic variability in CYP2B6 and CYP2C19 to the treatment efficacy and cyclophosphamide-induced side effects in pediatric neuroblastoma patients under cyclophosphamide therapy (N = 50). Cyclophosphamide-induced hematologic toxicities were pivotal in all patients, whereas only moderate hepatorenal toxicity was developed. The patients' CYP2B6 metabolizer phenotypes were associated with the occurrence of lymphopenia, thrombocytopenia, and monocytopenia as well as of liver injury, but not with kidney or urinary bladder (hemorrhagic cystitis) toxicities. Furthermore, the patients' age (< 1.5 years, P = 0.03) and female gender (P ≤ 0.02), but not CYP2B6 or CYP2C19 metabolizer phenotypes appeared as significant prognostic factors in treatment outcomes. Our results may contribute to a better understanding of the impact of CYP2B6 variability on cyclophosphamide-induced side effects.

Simultaneous absolute protein quantification of seven cytochrome P450 isoforms in rat liver microsomes by LC-MS/MS-based isotope internal standard method

The cytochrome P450 (CYP) enzymes play a pivotal role in drug metabolism. LC-MS/MS-based targeting technology has been applied to the analysis of CYP enzymes, promoting drug development and drug-drug interaction studies. Rat is one of the most commonly used models for drug metabolism assessment, but LC-MS/MS assay quantifying the abundance of CYP enzymes in rats is rarely reported. Herein, an accurate and stable LC-MS/MS based method was developed and validated for the simultaneous quantification of seven major rat CYP isoforms (CYP1A2, 2B1, 2C6, 2C11, 2D1, 2E1, and 3A1) in liver microsomes. The careful optimization of trypsin digestion and chromatography combined with isotope-labeled peptide as internal standard improved the efficiency and accuracy of the analysis. Highly specific surrogate peptides were obtained by a procedure including trypsin digestion for six hours and separated on a Hypersil Gold C18 column (100 × 2.1 mm, 3 μm) using gradient elution for 15 min with a mobile phase of water containing 0.1% formic acid and acetonitrile. In the method validation, linearity, matrix effect, recovery, stability, accuracy, and precision all meet the requirements. Subsequently, this method was applied to detect seven enzymes in rat liver microsomes from four different sources, and the correlation between the abundance and activity of CYP enzymes was further analyzed. The high-throughput detection method provided in this study will provide support for pertinent pharmaceutical research based on rat models.

Cytochrome P450 enzymes in the pediatric population: Connecting knowledge on P450 expression with pediatric pharmacokinetics

Cytochrome P450 enzymes play an important role in the pharmacokinetics, efficacy, and toxicity of drugs. Age-dependent changes in P450 enzyme expression have been studied based on several detection systems, as well as by deconvolution of in vivo pharmacokinetic data observed in pediatric populations. The age-dependent changes in P450 enzyme expression can be important determinants of drug disposition in childhood, in addition to the changes in body size and the other physiological parameters, and effects of pharmacogenetics and disease on organ functions. As a tool incorporating drug-specific and body-specific factors, physiologically-based pharmacokinetic (PBPK) models have become increasingly used to characterize and explore mechanistic insights into drug disposition. Thus, PBPK models can be a bridge between findings from basic science and utilization in predictive science. Pediatric PBPK models incorporate additional system specific information on developmental physiology and ontogeny and have been used to predict pharmacokinetic parameters from preterm neonates onwards. These models have been advocated by regulatory authorities in order to support pediatric clinical trials. The purpose of this chapter is to highlight accumulated knowledge and findings from basic research focusing on P450 enzymes, as well as the current status and future challenges of expanding the utilization of pediatric PBPK modeling.

CYP2B6 allelic variants and non-genetic factors influence CYP2B6 enzyme function

Human CYP2B6 enzyme although constitutes relatively low proportion (1–4%) of hepatic cytochrome P450 content, it is the major catalyst of metabolism of several clinically important drugs (efavirenz, cyclophosphamide, bupropion, methadone). High interindividual variability in CYP2B6 function, contributing to impaired drug-response and/or adverse reactions, is partly elucidated by genetic polymorphisms, whereas non-genetic factors can significantly modify the CYP2B6 phenotype. The influence of genetic and phenoconverting non-genetic factors on CYP2B6-selective activity and CYP2B6 expression was investigated in liver tissues from Caucasian subjects (N = 119). Strong association was observed between hepatic S-mephenytoin N-demethylase activity and CYP2B6 mRNA expression (P < 0.0001). In less than one third of the tissue donors, the CYP2B6 phenotype characterized by S-mephenytoin N-demethylase activity and/or CYP2B6 expression was concordant with CYP2B6 genotype, whereas in more than 35% of the subjects, an altered CYP2B6 phenotype was attributed to phenoconverting non-genetic factors (to CYP2B6-specific inhibitors and inducers, non-specific amoxicillin + clavulanic acid treatment and chronic alcohol consumption, but not to the gender). Furthermore, CYP2B6 genotype–phenotype mismatch still existed in one third of tissue donors. In conclusion, identifying potential sources of CYP2B6 variability and considering both genetic variations and non-genetic factors is a pressing requirement for appropriate elucidation of CYP2B6 genotype–phenotype mismatch.

Ontogeny of Drug-Metabolizing Enzymes

Almost 50% of prescription drugs lack age-appropriate dosing guidelines and therefore are used "off-label." Only ~10% drugs prescribed to neonates and infants have been studied for safety or efficacy. Immaturity of drug metabolism in children is often associated with drug toxicity. This chapter summarizes data on the ontogeny of major human metabolizing enzymes involved in oxidation, reduction, hydrolysis, and conjugation of drugs. The ontogeny data of individual drug-metabolizing enzymes are important for accurate prediction of drug pharmacokinetics and toxicity in children. This information is critical for designing clinical studies to appropriately test pharmacological hypotheses and develop safer pediatric drugs, and to replace the long-standing practice of body weight- or surface area-normalized drug dosing. The application of ontogeny data in physiologically based pharmacokinetic model and regulatory submission are discussed.

CYP2B6 Functional Variability in Drug Metabolism and Exposure Across Populations-Implication for Drug Safety, Dosing, and Individualized Therapy

Adverse drug reactions (ADRs) are one of the major causes of morbidity and mortality worldwide. It is well-known that individual genetic make-up is one of the causative factors of ADRs. Approximately 14 million single nucleotide polymorphisms (SNPs) are distributed throughout the entire human genome and every patient has a distinct genetic make-up which influences their response to drug therapy. Cytochrome P450 2B6 (CYP2B6) is involved in the metabolism of antiretroviral, antimalarial, anticancer, and antidepressant drugs. These drug classes are commonly in use worldwide and face specific population variability in side effects and dosing. Parts of this variability may be caused by single nucleotide polymorphisms (SNPs) in the CYP2B6 gene that are associated with altered protein expression and catalytic function. Population variability in the CYP2B6 gene leads to changes in drug metabolism which may result in adverse drug reactions or therapeutic failure. So far more than 30 non-synonymous variants in CYP2B6 gene have been reported. The occurrence of these variants show intra and interpopulation variability, thus affecting drug efficacy at individual and population level. Differences in disease conditions and affordability of drug therapy further explain why some individuals or populations are more exposed to CYP2B6 pharmacogenomics associated ADRs than others. Variabilities in drug efficacy associated with the pharmacogenomics of CYP2B6 have been reported in various populations. The aim of this review is to highlight reports from various ethnicities that emphasize on the relationship between CYP2B6 pharmacogenomics variability and the occurrence of adverse drug reactions. In vitro and in vivo studies evaluating the catalytic activity of CYP2B6 variants using various substrates will also be discussed. While implementation of pharmacogenomic testing for personalized drug therapy has made big progress, less data on pharmacogenetics of drug safety has been gained in terms of CYP2B6 substrates. Therefore, reviewing the existing evidence on population variability in CYP2B6 and ADR risk profiles suggests that, in addition to other factors, the knowledge on pharmacogenomics of CYP2B6 in patient treatment may be useful for the development of personalized medicine with regards to genotype-based prescription.

Development and Application of a Life-Stage Physiologically Based Pharmacokinetic (PBPK) Model to the Assessment of Internal Dose of Pyrethroids in Humans

To address concerns around age-related sensitivity to pyrethroids, a life-stage physiologically based pharmacokinetic (PBPK) model, supported by in vitro to in vivo extrapolation (IVIVE) was developed. The model was used to predict age-dependent changes in target tissue exposure of 8 pyrethroids; deltamethrin (DLM), cis-permethrin (CPM), trans-permethrin, esfenvalerate, cyphenothrin, cyhalothrin, cyfluthrin, and bifenthrin. A single model structure was used based on previous work in the rat. Intrinsic clearance (CLint) of each individual cytochrome P450 or carboxylesterase (CES) enzyme that are active for a given pyrethroid were measured in vitro, then biologically scaled to obtain in vivo age-specific total hepatic CLint. These IVIVE results indicate that, except for bifenthrin, CES enzymes are largely responsible for human hepatic metabolism (>50% contribution). Given the high efficiency and rapid maturation of CESs, clearance of the pyrethroids is very efficient across ages, leading to a blood flow-limited metabolism. Together with age-specific physiological parameters, in particular liver blood flow, the efficient metabolic clearance of pyrethroids across ages results in comparable to or even lower internal exposure in the target tissue (brain) in children than that in adults in response to the same level of exposure to a given pyrethroid (Cmax ratio in brain between 1- and 25-year old = 0.69, 0.93, and 0.94 for DLM, bifenthrin, and CPM, respectively). Our study demonstrated that a life-stage PBPK modeling approach, coupled with IVIVE, provides a robust framework for evaluating age-related differences in pharmacokinetics and internal target tissue exposure in humans for the pyrethroid class of chemicals.

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.

PharmVar GeneFocus: CYP2B6

The Pharmacogene Variation Consortium (PharmVar) catalogs star (*) allele nomenclature for the polymorphic human CYP2B6 gene. Genetic variation within the CYP2B6 gene locus impacts the metabolism or bioactivation of clinically important drugs. Of particular importance are efficacy and safety concerns regarding: efavirenz, which is used for the treatment of HIV type-1 infection; methadone, a mainstay in the treatment of opioid use disorder and as an analgesic; ketamine, used as an antidepressant and analgesic; and bupropion, which is prescribed to treat depression and for smoking cessation. This GeneFocus provides a comprehensive overview and summary of CYP2B6 and describes how haplotype information catalogued by PharmVar is utilized by the Pharmacogenomics Knowledgebase (PharmGKB) and the Clinical Pharmacogenetics Implementation Consortium (CPIC).

Clinical pharmacology and dosing regimen optimization of neonatal opioid withdrawal syndrome treatments

In this paper, we review the management of neonatal opioid withdrawal syndrome (NOWS) and clinical pharmacology of primary treatment agents in NOWS, including morphine, methadone, buprenorphine, clonidine, and phenobarbital. Pharmacologic treatment strategies in NOWS have been mostly empirical, and heterogeneity in dosing regimens adds to the difficulty of extrapolating study results to broader patient populations. As population pharmacokinetics (PKs) of pharmacologic agents in NOWS become more well‐defined and knowledge of patient‐specific factors affecting treatment outcomes continue to accumulate, PK/pharmacodynamic modeling and simulation will be powerful tools to aid the design of optimal dosing regimens at the patient level. Although there is an increasing number of clinical trials on the comparative efficacy of treatment agents in NOWS, here, we also draw attention to the importance of optimizing the dosing regimen, which can be arguably equally important at identifying the optimal treatment agent.

Potential implications of DMET ontogeny on the disposition of commonly prescribed drugs in neonatal and pediatric intensive care units

Introduction: Pediatric patients, especially neonates and infants, are more susceptible to adverse drug events as compared to adults. In particular, immature small molecule drug metabolism and excretion can result in higher incidences of pediatric toxicity than adults if the pediatric dose is not adjusted.Area covered: We reviewed the top 29 small molecule drugs prescribed in neonatal and pediatric intensive care units and compiled the mechanisms of their metabolism and excretion. The ontogeny of Phase I and II drug metabolizing enzymes and transporters (DMETs), particularly relevant to these drugs, are summarized. The potential effects of DMET ontogeny on the metabolism and excretion of the top pediatric drugs were predicted. The current regulatory requirements and recommendations regarding safe and effective use of drugs in children are discussed. A few representative examples of the use of ontogeny-informed physiologically based pharmacokinetic (PBPK) models are highlighted.Expert opinion: Empirical prediction of pediatric drug dosing based on body weight or body-surface area from the adult parameters can be inaccurate because DMETs are not mature in children and the age-dependent maturation of these proteins is different. Ontogeny-informed-PBPK modeling provides a better alternative to predict the pharmacokinetics of drugs in children.

Sources of Interindividual Variability

The efficacy, safety, and tolerability of drugs are dependent on numerous factors that influence their disposition. A dose that is efficacious and safe for one individual may result in sub-therapeutic or toxic blood concentrations in others. A significant source of this variability in drug response is drug metabolism, where differences in presystemic and systemic biotransformation efficiency result in variable degrees of systemic exposure (e.g., AUC, C_max, and/or C_min) following administration of a fixed dose.Interindividual differences in drug biotransformation have been studied extensively. It is recognized that both intrinsic factors (e.g., genetics, age, sex, and disease states) and extrinsic factors (e.g., diet , chemical exposures from the environment, and the microbiome) play a significant role. For drug-metabolizing enzymes, genetic variation can result in the complete absence or enhanced expression of a functional enzyme. In addition, upregulation and downregulation of gene expression, in response to an altered cellular environment, can achieve the same range of metabolic function (phenotype), but often in a less predictable and time-dependent manner. Understanding the mechanistic basis for variability in drug disposition and response is essential if we are to move beyond the era of empirical, trial-and-error dose selection and into an age of personalized medicine that will improve outcomes in maintaining health and treating disease.

Ontogeny of Phase I Metabolism of Drugs

Capturing ontogeny of enzymes involved in phase I metabolism is crucial to improve prediction of dose-concentration and concentration-effect relationships throughout infancy and childhood. Once captured, these patterns can be integrated in semiphysiologically or physiology-based pharmacokinetic models to support predictions in specific pediatric settings or to support pediatric drug development. Although these translational efforts are crucial, isoenzyme-specific ontogeny-based models should also incorporate data on variability of maturational and nonmaturational covariates (eg, disease, treatment modalities, pharmacogenetics). Therefore, this review provides a summary of the ontogeny of phase I drug-metabolizing enzymes, indicating current knowledge gaps and recent progresses. Furthermore, we tried to illustrate that straightforward translation of isoenzyme-specific ontogeny to predictions does not allow full exploration of scenarios of potential variability related to maturational (non-age-related variability, other isoenzymes or transporters) or nonmaturational (disease, pharmacogenetics) covariates, and necessitates integration in a "systems" concept.

Impact of cytochrome P450 variation on meperidine N-demethylation to the neurotoxic metabolite normeperidine

1. Meperidine is an opioid analgesic that undergoes N-demethylation to form the neurotoxic metabolite normeperidine. Previous studies indicate that meperidine N-demethylation is catalyzed by cytochrome P450 2B6 (CYP2B6), CYP3A4, and CYP2C19.2. The purpose of this study was to examine the relative P450 contributions to meperidine N-demethylation and to evaluate the effect of CYP2C19 polymorphism on normeperidine generation. Experiments were performed using recombinant P450 enzymes, selective chemical inhibitors, enzyme kinetic assays, and correlation analysis with individual CYP2C19-genotyped human liver microsomes.3. The catalytic efficiency (k_cat/K_m) for meperidine N-demethylation was similar between recombinant CYP2B6 and CYP2C19, but markedly lower by CYP3A4.4. In CYP2C19-genotyped human liver microsomes, normeperidine formation was significantly correlated with CYP2C19 activity (S-mephenytoin 4´-hydroxylation).5. CYP2C19 inhibitor (+)-N-3-benzylnirvanol and CYP3A inhibitor ketoconazole significantly reduced microsomal normeperidine generation by an individual donor with high CYP2C19 activity, whereas donors with lower CYP2C19 activity were sensitive to inhibition by ketoconazole but not benzylnirvanol.6. These findings demonstrate that the relative CYP3A4, CYP2B6, and CYP2C19 involvement in meperidine N-demethylation depends on the enzyme activities in individual human liver microsomal samples. CYP2C19 is likely an important contributor to normeperidine generation in individuals with high CYP2C19 activity, but additional factors influence inter-individual metabolite accumulation.

Clinical Pharmacogenetics Implementation Consortium (CPIC) Guideline for CYP2B6 and Efavirenz-Containing Antiretroviral Therapy

The HIV type-1 nonnucleoside reverse transcriptase inhibitor, efavirenz, is widely used to treat HIV type-1 infection. Efavirenz is predominantly metabolized into inactive metabolites by cytochrome P450 (CYP)2B6, and patients with certain CYP2B6 genetic variants may be at increased risk for adverse effects, particularly central nervous system toxicity and treatment discontinuation. We summarize the evidence from the literature and provide therapeutic recommendations for efavirenz prescribing based on CYP2B6 genotypes.

Model qualification of the PK-Sim® pediatric module for pediatric exposure assessment of CYP450 metabolized compounds

Pediatric physiologically based pharmacokinetic (PBPK) models facilitate the estimation of pharmacokinetic (PK) parameters in children under specific exposure conditions. In human health risk assessment, PBPK modeling has been used to determine a chemical-specific human kinetic adjustment factor (HKAF). Due to increased demands in regulatory assessment, model evaluation and qualification have gained growing attention. The aim of this study was to undertake model qualification of pediatric PBPK models for compounds that are primarily metabolized by cytochrome P450 (CYP) enzymes. The objectives were to determine the appropriateness of the virtual individual creating algorithm in PK-Sim® in predicting PK parameters and their variability in children and identify critical system-specific inputs. PBPK models in adults were constructed for several pharmaceuticals (grouped by major clearance process such as CYP3A4). Several age groups of virtual individuals were created to represent children in pediatric clinical studies. The mean and variance of clearance (CL) from virtual populations were compared to observed values. Sensitivity analysis on area under the curve (AUC) was performed. System-specific parameters of virtual children that contribute to inter-individual PK properties were assessed. Eighty-one percent of the comparisons between simulated and observed clearance values were within twofold error. The mean fold errors were 1.1, 1, 0.7 and 1.8 in adolescents, children, infants and neonates, respectively. CL variability was reasonably predicted for 70% of the comparisons with comparable coefficients of variation between observed and predicted. The sensitivity analysis revealed that fraction unbound in plasma, parameters related to CYP enzyme-mediated metabolism and liver volumewere most important in the estimation of pediatric exposure. A comparison of variabilities in weight, height and liver volume in virtual children showed reliable agreement with observed data. The presented results of predictive performance and properties of virtual populations provide confidence in the use of PK-Sim for pediatric PBPK modeling in toxicological applications including PBPK-based-HKAF derivation.

Estimation of Drug Pharmacokinetics from Breast Feeding: A Simple Method Based on Meta-analysis

Background:

In resource-limited settings, breastfeeding is the healthiest source of nutrition for newborns. For economic/cultural reasons, breastfeeding is the preferred option for the majority of mothers, including HIV-positive mothers.

Objective:

The objective of this review is to document parameters characterizing antiretroviral therapy (ARV) diffusion into breast milk associated with the estimated ARV amount ingested by breastfed infant and clinical/biological abnormalities.

Data Source and Eligibility Criteria:

Twenty seven (27) published articles on the aspects of Pharmacokinetic parameters on ARV diffusion into breast milk have shown a large variability without clear interpretation on drugs diffusion. Using PubMed and Embase, we conducted a search to identify all published studies at 2015 that characterized antiretroviral drug diffusion from mother to infant via breast milk. We identified 27 published studies that characterized antiretroviral drug passage from mother to infant (drug concentrations in mother’s milk and breastfed plasma). Information was sufficiently complete for inclusion in the present analysis for only six antiretroviral drugs.

Results:

Finally, only data for nevirapine and efavirenz were exploitable because some of the studies found null or non-detectable levels, which were not suitable for simulations. Median (IQR) nevirapine CL/F were 0.022 (0.013–0.038) for newborns, 0.121 (0.116–0.125) for children and 0.056 (0.045–0.070) for mothers, all in L/h/kg. Efavirenz CL/F were 0.025 (0.016–0.039) for newborns, 0.273 (0.261–0.285) for children and 0.160 (0.153–0.167) for mothers, also in L/h/kg.

Conclusion:

Pharmacokinetics parameters of efavirenz and nevirapine are important to be determined in breastfed newborns.

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.

The importance of both CYP2C19 and CYP2B6 germline variations in cyclophosphamide pharmacokinetics and clinical outcomes

Cyclophosphamide is an alkylating agent used in the treatment of solid and haematological malignancies and as an immunosuppressive agent. As a prodrug, it is dependent on bioactivation to the active phosphoramide mustard metabolite to elicit its therapeutic effect. This focused review will highlight the evidence for the role of germline pharmacogenetic variation in both plasma pharmacokinetics and clinical outcomes. There is a substantial indication from 13 pharmacokinetic and 17 therapeutic outcome studies, in contexts as diverse as haematological malignancy, breast cancer, systemic lupus erythematosus and myeloablation, that pharmacogenetic variation in both CYP2C19 and CYP2B6 influence the bioactivation of cyclophosphamide. An additional role for pharmacogenetic variation in ALDH1A1 has also been reported. Future studies should comprehensively assess these 3 pharmacogenes and undertake appropriate statistical analysis of gene-gene interactions to confirm these findings and may allow personalised treatment regimens.

Methadone dosing strategies in preterm neonates can be simplified

AIMS

A dramatic increase in newborn infants with neonatal abstinence syndrome has been observed and these neonates are frequently treated with complex methadone dosing schemes to control their withdrawal symptoms. Despite its abundant use, hardly any data on the pharmacokinetics (PK) of methadone is available in preterm neonates. Therefore we investigated developmental PK of methadone and evaluated current dosing strategies and possible simplification in this vulnerable population.

METHODS

A single-centre open-label prospective study was performed to collect PK data after a single oral dose of methadone in preterm neonates. A population PK model was built to characterize developmental PK of (R)- and (S)-methadone. Model-based simulations were performed to identify a simplified dosing strategy to reach and maintain target methadone exposure.

RESULTS

A total of 121 methadone concentrations were collected from 31 preterm neonates. A one-compartment model with first order absorption and elimination kinetics best described PK data for (R)- and (S)-methadone. Clearance increases with advancing gestational age and differs between R- and S-enantiomer, being slightly higher for the former (0.244 vs 0.167 L/h). Preterm neonates reached target exposure after 48 hours with currently used dosing schedules. Output from simulations revealed that target exposures can be achieved with a simplified dosing strategy during the first 4 days of treatment.

CONCLUSION

Methadone clearance in preterm neonates increases with advancing gestational age and its disposition is influenced by its chirality. Simulations that account for developmental PK changes indicate a shorter methadone dosing strategy can maintain target exposure to control withdrawal symptoms.

Effect of Rifampin-Isoniazid-Containing Antituberculosis Therapy on Efavirenz Pharmacokinetics in HIV-Infected Children 3 to 14 Years Old

We compared efavirenz pharmacokinetic (PK) parameters in children with tuberculosis (TB)/human immunodeficiency virus (HIV) coinfection on and off first-line antituberculosis therapy to that in HIV-infected children. Children 3 to 14 years old with HIV infection, with and without TB, were treated with standard efavirenz-based antiretroviral therapy without any efavirenz dose adjustments. The new World Health Organization-recommended antituberculosis drug dosages were used in the coinfected participants. Steady-state efavirenz concentrations after 4 weeks of antiretroviral therapy were measured using validated liquid chromatography with tandem mass spectrometry (LC-MS/MS) assays. Pharmacokinetic parameters were calculated using noncompartmental analysis. Between groups, PK parameters were compared by Wilcoxon rank-sum test and within group by signed-rank test. Of the 105 participants, 43 (41.0%) had TB coinfection. Children with TB/HIV coinfection compared to those with HIV infection were younger, had lower median weight-for-age Z score, and received a higher median efavirenz weight-adjusted dose. Geometric mean (GM) efavirenz peak concentration (Cmax), concentration at 12 h (C12h), Cmin, and total area under the curve from time 0 to 24 h (AUC0-24h) values were similar in children with HIV infection and those with TB/HIV coinfection during anti-TB therapy. Geometric mean efavirenz C12h, Cmin, and AUC0-24h values were lower in TB/HIV-coinfected patients off anti-TB therapy than in the children with HIV infection or TB/HIV coinfection on anti-TB therapy. Efavirenz clearance was lower and AUC0-24h was higher on than in patients off anti-TB therapy. Reduced efavirenz clearance by first-line anti-TB therapy at the population level led to similar PK parameters in HIV-infected children with and without TB coinfection. Our findings do not support modification of efavirenz weight-band dosing guidelines based on TB coinfection status in children. (The study was registered with ClinicalTrials.gov under registration number NCT01704144.).

Age- and Genotype-Dependent Variability in the Protein Abundance and Activity of Six Major Uridine Diphosphate-Glucuronosyltransferases in Human Liver

The ontogeny of hepatic uridine diphosphate-glucuronosyltransferases (UGTs) was investigated by determining their protein abundance in human liver microsomes isolated from 136 pediatric (0-18 years) and 35 adult (age >18 years) donors using liquid chromatography / tandem mass spectrometry (LC-MS/MS) proteomics. Microsomal protein abundances of UGT1A1, UGT1A4, UGT1A6, UGT1A9, UGT2B7, and UGT2B15 increased by ∼8, 55, 35, 33, 8, and 3-fold from neonates to adults, respectively. The estimated age at which 50% of the adult protein abundance is observed for these UGT isoforms was between 2.6-10.3 years. Measured in vitro activity was generally consistent with the protein data. UGT1A1 protein abundance was associated with multiple single nucleotide polymorphisms exhibiting noticeable ontogeny-genotype interplay. UGT2B15 rs1902023 (*2) was associated with decreased protein activity without any change in protein abundance. Taken together, these data are invaluable to facilitate the prediction of drug disposition in children using physiologically based pharmacokinetic modeling as demonstrated here for zidovudine and morphine.

PBPK and its Virtual Populations: the Impact of Physiology on Pediatric Pharmacokinetic Predictions of Tramadol

In pediatric PBPK models, age-related changes in the body are known to occur. Given the sparsity of and the variability associated with relevant physiological parameters, different PBPK software providers may vary in their system's data. In this work, three commercially available PBPK software packages (PK-Sim®, Simcyp®, and Gastroplus®) were investigated regarding their differences in system-related information, possibly affecting clearance prediction. Three retrograde PBPK clearance models were set up to enable prediction of pediatric tramadol clearance. These models were qualified in terms of total, CYP2D6, and renal clearance in adults. Tramadol pediatric clearance predictions from PBPK were compared with a pooled popPK model covering clearance ranging from neonates to adults. Fold prediction errors were used to evaluate the results. Marked differences in liver clearance prediction between PBPK models were observed. In general, the prediction bias of total clearance was greatest at the youngest population and decreased with age. Regarding CYP2D6 and renal clearance, important differences exist between PBPK software tools. Interestingly, the PBPK model with the shortest CYP2D6 maturation half-life (PK-Sim) agreed best with the in vivo CYP2D6 maturation model. Marked differences in physiological data explain the observed differences in hepatic clearance prediction in early life between the various PBPK software providers tested. Consensus on the most suited pediatric data to use should harmonize and optimize pediatric clearance predictions. Moreover, the combination of bottom-up and top-down approaches, using a convenient probe substrate, has the potential to update system-related parameters in order to better represent pediatric physiology.

Efavirenz and Lopinavir Levels in HIV-Infected Women and Their Nursing Infants, in Mali

Limited data are currently available on antiretroviral pharmacokinetics in breast milk (BM) and in breastfed infants' blood. To explore these parameters in patients in Mali, we measured plasma antiretroviral levels in human immunodeficiency virus (HIV)-infected mothers and their breastfed infants over 6 months. We specifically analyzed the concentrations of efavirenz (EFV) and lopinavir (LPV) in the plasma of mothers living with HIV and their breastfed infants. Blood samples were collected at delivery and at month 1, 3, and 6 postpartum. EFV and LPV concentrations were measured by liquid chromatography-tandem mass spectrometry. HIV-1 RNA load was measured by Abbott M2000RT RealTime System at delivery and 6 months postpartum for mothers, and at 3 and 6 months postbirth for infants. The median duration of antiretroviral therapy at study inclusion was 57 months [interquartile range (IQR), 0-168 months]. The median EFV ratios of infant plasma/maternal plasma (MP) were 0.057 at month 1, 0.072 at month 3, and 0.048 at month 6. During the study period, the median BM/MP ratio of EFV was 1.16 (IQR, 0.96-20.62), which corresponds to a relative infant dose of 2.46% of the recommended weight-adjusted pediatric EFV dose at month 6. The apparent infant clearance of EFV was 0.146 l/h per kilogram at month 6. The LPV concentrations in the plasma of all infants were undetectable. No drug-related adverse reaction or toxicity was observed in any of the infants. The two women who presented a viral load of >50 copies/ml at month 6 had undetectable plasma drug concentrations at the same period. This study showed that breastfed infants received a low level of EFV but not LPV from their treated mothers.

Variants in the CYP2B6 3'UTR Alter In Vitro and In Vivo CYP2B6 Activity: Potential Role of MicroRNAs

CYP2B6*6 and CYP2B6*18 are the most clinically important variants causing reduced CYP2B6 protein expression and activity. However, these variants do not account for all variability in CYP2B6 activity. Emerging evidence has shown that genetic variants in the 3'UTR may explain variable drug response by altering microRNA regulation. Five 3'UTR variants were associated with significantly altered efavirenz AUC0-48 (8-OH-EFV/EFV) ratios in healthy human volunteers. The rs70950385 (AG>CA) variant, predicted to create a microRNA binding site for miR-1275, was associated with a 33% decreased CYP2B6 activity among normal metabolizers (AG/AG vs. CA/CA (P < 0.05)). In vitro luciferase assays were used to confirm that the CA on the variant allele created a microRNA binding site causing an 11.3% decrease in activity compared to the AG allele when treated with miR-1275 (P = 0.0035). Our results show that a 3'UTR variant contributes to variability in CYP2B6 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.

The Ontogeny of Cytochrome P450 Enzyme Activity and Protein Abundance in Conventional Pigs in Support of Preclinical Pediatric Drug Research

Since the implementation of several legislations to improve pediatric drug research, more pediatric clinical trials are being performed. In order to optimize these pediatric trials, adequate preclinical data are necessary, which are usually obtained by juvenile animal models. The growing piglet has been increasingly suggested as a potential animal model due to a high degree of anatomical and physiological similarities with humans. However, physiological data in pigs on the ontogeny of major organs involved in absorption, distribution, metabolism, and excretion of drugs are largely lacking. The aim of this study was to unravel the ontogeny of porcine hepatic drug metabolizing cytochrome P450 enzyme (CYP450) activities as well as protein abundances. Liver microsomes from 16 conventional pigs (8 males and 8 females) per age group: 2 days, 4 weeks, 8 weeks, and 6-7 months were prepared. Activity measurements were performed with substrates of major human CYP450 enzymes: midazolam (CYP3A), tolbutamide (CYP2C), and chlorzoxazone (CYP2E). Next, the hepatic scaling factor, microsomal protein per gram liver (MPPGL), was determined to correct for enzyme losses during the fractionation process. Finally, protein abundance was determined using proteomics and correlated with enzyme activity. No significant sex differences within each age category were observed in enzyme activity or MPPGL. The biotransformation rate of all three substrates increased with age, comparable with human maturation of CYP450 enzymes. The MPPGL decreased from birth till 8 weeks of age followed by an increase till 6-7 months of age. Significant sex differences in protein abundance were observed for CYP1A2, CYP2A19, CYP3A22, CYP4V2, CYP2C36, CYP2E_1, and CYP2E_2. Midazolam and tolbutamide are considered good substrates to evaluate porcine CYP3A/2C enzymes, respectively. However, chlorzoxazone is not advised to evaluate porcine CYP2E enzyme activity. The increase in biotransformation rate with age can be attributed to an increase in absolute amount of CYP450 proteins. Finally, developmental changes were observed regarding the involvement of specific CYP450 enzymes in the biotransformation of the different substrates.

The impact of CYP2B6 polymorphisms on the interactions of efavirenz with lumefantrine: Implications for paediatric antimalarial therapy

Lumefantrine is a widely used antimalarial in children in sub-Saharan Africa and is predominantly metabolised by CYP3A4. The concomitant use of lumefantrine with the antiretroviral efavirenz, which is metabolised by CYP2B6 and is an inducer of CYP3A4, increases the risk of lumefantrine failure and can result in an increased recrudescence rate in HIV-infected children. This is further confounded by CYP2B6 being highly polymorphic resulting in a 2-3 fold higher efavirenz plasma concentration in polymorphic subjects, which enhances the potential for an efavirenz-lumefantrine drug-drug interaction (DDI). This study developed a population-based PBPK model capable of predicting the impact of efavirenz-mediated DDIs on lumefantrine pharmacokinetics in African paediatric population groups, which also considered the polymorphic nature of CYP2B6. The validated model demonstrated a significant difference in lumefantrine target day 7 concentrations (C_d7) in the presence and absence of efavirenz and confirmed the capability of efavirenz to initiate this DDI. This was more apparent in the *6/*6 compared to *1/*1 population group and resulted in a significantly lower (P < 0.001) lumefantrine C_d7. A prospective change in dosing schedule from 3-days to 7-days resulted in a greater number of *6/*6 subjects (28-57%) attaining the target C_d7 across age bands (0.25-13 years), with the greatest increase evident in the 1-4 year old group (3-day: 1%; 7-day: 28%).

Challenges and Opportunities for Increasing the Knowledge Base Related to Drug Biotransformation and Pharmacokinetics during Growth and Development

It is generally acknowledged that there is a need and role for informative pharmacokinetic models to improve predictions and simulation as well as individualization of drug therapy in pediatric populations of different ages and developmental stages. This special issue contains more than 20 papers responding to the challenge of providing new information on scaling factors, ontogeny functions for drug metabolizing enzymes and transporters, the mechanisms underlying the observed developmental trajectories for these gene products, age-dependent changes in physiologic processes affecting drug disposition in children, as well as in vitro and in vivo studies describing the relative contribution of ontogeny and genetic factors as sources of variability in drug disposition in children. Considered together, these contributions serve to illustrate some of the current limitations regarding sample availability, number, and quality, but also provide a framework that allows for the potential value of the results of a given study to be interpreted within the context of these limitations. Among the challenges for the future are improving our understanding of the mechanisms regulating age-dependent changes in factors influencing drug disposition and response, thereby facilitating generalization to systems lacking detailed data, better integrating age-dependent changes in pharmacokinetics with age-dependent changes in pharmacodynamics, and allowing better predictability and individualization of drug disposition and response across the pediatric age spectrum.

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.