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

Physiologically Based Pharmacokinetic Modeling in Neonates: Current Status and Future Perspectives

Rational drug use in special populations is a clinical problem that doctors and pharma-cists must consider seriously. Neonates are the most physiologically immature and vulnerable to drug dosing. There is a pronounced difference in the anatomical and physiological profiles be-tween neonates and older people, affecting the absorption, distribution, metabolism, and excretion of drugs in vivo, ultimately leading to changes in drug concentration. Thus, dose adjustments in neonates are necessary to achieve adequate therapeutic concentrations and avoid drug toxicity. Over the past few decades, modeling and simulation techniques, especially physiologically based pharmacokinetic (PBPK) modeling, have been increasingly used in pediatric drug development and clinical therapy. This rigorously designed and verified model can effectively compensate for the deficiencies of clinical trials in neonates, provide a valuable reference for clinical research design, and even replace some clinical trials to predict drug plasma concentrations in newborns. This review introduces previous findings regarding age-dependent physiological changes and pathological factors affecting neonatal pharmacokinetics, along with their research means. The application of PBPK modeling in neonatal pharmacokinetic studies of various medications is also reviewed. Based on this, we propose future perspectives on neonatal PBPK modeling and hope for its broader application.

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.

Drug-drug interactions that alter the exposure of glucuronidated drugs: Scope, UDP-glucuronosyltransferase (UGT) enzyme selectivity, mechanisms (inhibition and induction), and clinical significance

Drug-drug interactions (DDIs) arising from the perturbation of drug metabolising enzyme activities represent both a clinical problem and a potential economic loss for the pharmaceutical industry. DDIs involving glucuronidated drugs have historically attracted little attention and there is a perception that interactions are of minor clinical relevance. This review critically examines the scope and aetiology of DDIs that result in altered exposure of glucuronidated drugs. Interaction mechanisms, namely inhibition and induction of UDP-glucuronosyltransferase (UGT) enzymes and the potential interplay with drug transporters, are reviewed in detail, as is the clinical significance of known DDIs. Altered victim drug exposure arising from modulation of UGT enzyme activities is relatively common and, notably, the incidence and importance of UGT induction as a DDI mechanism is greater than generally believed. Numerous DDIs are clinically relevant, resulting in either loss of efficacy or an increased risk of adverse effects, necessitating dose individualisation. Several generalisations relating to the likelihood of DDIs can be drawn from the known substrate and inhibitor selectivities of UGT enzymes, highlighting the importance of comprehensive reaction phenotyping studies at an early stage of drug development. Further, rigorous assessment of the DDI liability of new chemical entities that undergo glucuronidation to a significant extent has been recommended recently by regulatory guidance. Although evidence-based approaches exist for the in vitro characterisation of UGT enzyme inhibition and induction, the availability of drugs considered appropriate for use as 'probe' substrates in clinical DDI studies is limited and this should be research priority.

Infant Exposure to Antituberculosis Drugs via Breast Milk and Assessment of Potential Adverse Effects in Breastfed Infants: Critical Review of Data

Infants of mothers treated for tuberculosis might be exposed to drugs via breast milk. The existing information on the exposure of breastfed infants lacks a critical review of the published data. We aimed to evaluate the quality of the existing data on antituberculosis (anti-TB) drug concentrations in the plasma and milk as a methodologically sound basis for the potential risk of breastfeeding under therapy. We performed a systematic search in PubMed for bedaquiline, clofazimine, cycloserine/terizidone, levofloxacin, linezolid, pretomanid/pa824, pyrazinamide, streptomycin, ethambutol, rifampicin and isoniazid, supplemented with update references found in LactMed^®. We calculated the external infant exposure (EID) for each drug and compared it with the recommended WHO dose for infants (relative external infant dose) and assessed their potential to elicit adverse effects in the breastfed infant. Breast milk concentration data were mainly not satisfactory to properly estimate the EID. Most of the studies suffer from limitations in the sample collection, quantity, timing and study design. Infant plasma concentrations are extremely scarce and very little data exist documenting the clinical outcome in exposed infants. Concerns for potential adverse effects in breastfed infants could be ruled out for bedaquiline, cycloserine/terizidone, linezolid and pyrazinamide. Adequate studies should be performed covering the scenario in treated mothers, breast milk and infants.

Contribution of UGT Enzymes to Human Drug Metabolism Stereoselectivity: A Case Study of Medetomidine, RO5263397, Propranolol, and Testosterone

The enantiomeric forms of chiral compounds have identical physical properties but may vary greatly in their metabolism by individual enzymes. Enantioselectivity in UDP-glucuronosyl transferase (UGT) metabolism has been reported for a number of compounds and with different UGT isoforms involved. However, the impact of such individual enzyme results on overall clearance stereoselectivity is often not clear. The enantiomers of medetomidine, RO5263397, and propranolol and the epimers testosterone and epitestosterone exhibit more than a 10-fold difference in glucuronidation rates by individual UGT enzymes. In this study, we examined the translation of human UGT stereoselectivity to hepatic drug clearance considering the combination of multiple UGTs to overall glucuronidation, the contribution of other metabolic enzymes such as cytochrome P450s (P450s), and the potential for differences in protein binding and blood/plasma partitioning. For medetomidine and RO5263397, the high individual enzyme (UGT2B10) enantioselectivity translated into ∼3- to >10-fold differences in predicted human hepatic in vivo clearance. For propranolol, the UGT enantioselectivity was irrelevant in the context of high P450 metabolism. For testosterone, a complex picture emerged due to differential epimeric selectivity of various contributing enzymes and potential for extrahepatic metabolism. Quite different patterns of P450- and UGT-mediated metabolism were observed across species, as well as differences in stereoselectivity, indicating that extrapolation from human enzyme and tissue data are essential when predicting human clearance enantioselectivity. SIGNIFICANCE STATEMENT: Individual enzyme stereoselectivity illustrates the importance of three-dimensional drug-metabolizing enzyme-substrate interactions and is essential when considering the clearance of racemic drugs. However, translation from in vitro to in vivo can be challenging as contributions from multiple enzymes and enzyme classes must be combined with protein binding and blood/plasma partitioning data to estimate the net intrinsic clearance for each enantiomer. Preclinical species may be misleading as enzyme involvement and metabolism stereoselectivity can differ substantially.

In silico modeling for the hepatic circulation and transport: From the liver organ to lobules

The function of the liver depends critically on its blood supply. Numerous in silico models have been developed to study various aspects of the hepatic circulation, including not only the macro-hemodynamics at the organ level, but also the microcirculation at the lobular level. In addition, computational models of blood flow and bile flow have been used to study the transport, metabolism, and clearance of drugs in pharmacokinetic studies. These in silico models aim to provide insights into the liver organ function under both healthy and diseased states, and to assist quantitative analysis for surgical planning and postsurgery treatment. The purpose of this review is to provide an update on state-of-the-art in silico models of the hepatic circulation and transport processes. We introduce the numerical methods and the physiological background of these models. We also discuss multiscale frameworks that have been proposed for the liver, and their linkage with the large context of systems biology, systems pharmacology, and the Physiome project. This article is categorized under: Metabolic Diseases > Computational Models Metabolic Diseases > Biomedical Engineering Cardiovascular Diseases > Computational Models.

Population pharmacokinetics of mycophenolate mofetil in pediatric patients early after liver transplantation

Objective: To investigate the factors influencing the pharmacokinetics of mycophenolate mofetil (MMF) in pediatric patients after liver transplantation, and to establish a population pharmacokinetics model, which can provide a reference for clinical dosage adjustment.

Methods: A prospective study in a single center was performed on pediatric patients who were administrated with mycophenolate mofetil dispersible tablets (MMFdt) for at least 4 days after liver transplantation continuously. Blood samples were collected in ethylene diamine tetraacetic acid anticoagulant tubes before dosing and 0.5, 1, 2, 4, 8, and 12 h after the morning intake of MMFdt. The concentrations of mycophenolic acid (MPA) in plasma were assayed with a validated reverse-phase high-performance liquid chromatography method. UGT1A8 518C > G, UGT1A9 -275T > A, UGT1A9 -2152C > T, UGT2B7 211G > T, SLC O 1B1 521T > C polymorphism were determined by Sanger sequencing. Nonlinear mixed effects modeling was used to establish the population pharmacokinetics (PPK) model. The predictability and stability of the model were internally evaluated by the goodness of fit plots, visual prediction check, normalized prediction errors, and bootstraps.

Results: A two-compartment model with first-order absorption and first-order elimination was established with 115 MPA concentrations from 20 pediatric patients. The final model were: CL/F (L/h) = 14.8×(WT/7.5)^0.75×(DOSE/11.16)^0.452×е^0.06, Ka (h⁻¹) = 2.02×(WT/7.5)^−0.25, Vc/F (L) = 6.01×(WT/7.5), Vp/F (L) = 269 (fixed), Q/F (L/h) = 15.4×(WT/7.5)^0.75×е^1.39. Where CL/F was the apparent clearance rate, Ka was the absorption rate constant, Vc/F was the apparent distribution volume of the central compartment, Vp/F was the apparent distribution volume of the peripheral compartment, Q/F was the atrioventricular clearance rate, WT was the body weight of the subject, and DOSE was the MMFdt administered dose. The model indicated there was large inter-individual variability in CL/F and Q/F after multiple dosing of MMFdt. Internal evaluation results showed that the final model had good stability and prediction performance.

Conclusion: A stable and predictive population pharmacokinetic model of MMFdt in pediatric patients after the early stage of liver transplantation was established. The pediatric patient’s weight and the dose of MMFdt can be a reference to adjust the MMFdt dose.

Applications, Challenges, and Outlook for PBPK Modeling and Simulation: A Regulatory, Industrial and Academic Perspective

Several regulatory guidances on the use of physiologically based pharmacokinetic (PBPK) analyses and physiologically based biopharmaceutics model(s) (PBBM(s)) have been issued. Workshops are routinely held, demonstrating substantial interest in applying these modeling approaches to address scientific questions in drug development. PBPK models and PBBMs have remarkably contributed to model-informed drug development (MIDD) such as anticipating clinical PK outcomes affected by extrinsic and intrinsic factors in general and specific populations. In this review, we proposed practical considerations for a "base" PBPK model construction and development, summarized current status, challenges including model validation and gaps in system models, and future perspectives in PBPK evaluation to assess a) drug metabolizing enzyme(s)- or drug transporter(s)- mediated drug-drug interactions b) dosing regimen prediction, sampling timepoint selection and dose validation in pediatric patients from newborns to adolescents, c) drug exposure in patients with renal and/or and hepatic organ impairment, d) maternal-fetal drug disposition during pregnancy, and e) pH-mediated drug-drug interactions in patients treated with proton pump inhibitors/acid-reducing agents (PPIs/ARAs) intended for gastric protection. Since PBPK can simulate outcomes in clinical studies with enrollment challenges or ethical issues, the impact of PBPK models on waivers and how to strengthen study waiver is discussed.

A potential implication of UDP-glucuronosyltransferase 2B10 in the detoxification of drugs used in pediatric hematopoietic stem cell transplantation setting: an in silico investigation

Background

Sinusoidal occlusion syndrome (SOS) is a potentially severe complication following hematopoietic stem cell transplantation (HSCT) in pediatric patients. Treatment related risk factors such as intensity of conditioning, hepatotoxic co-medication and patient related factors such as genetic variants predispose individuals to develop SOS. The variant allele for SNP rs17146905 in UDP-glucuronosyl transferase 2B10 (UGT2B10) gene was correlated with the occurrence of SOS in an exome-wide association study. UGT2B10 is a phase II drug metabolizing enzyme involved in the N-glucuronidation of tertiary amine containing drugs.

Methods

To shed light on the functionality of UGT2B10 enzyme in the metabolism of drugs used in pediatric HSCT setting, we performed in silico screening against custom based library of putative ligands. First, a list of potential substrates for in silico analysis was prepared using a systematic consensus-based strategy. The list comprised of drugs and their metabolites used in pediatric HSCT setting. The three-dimensional structure of UGT2B10 was not available from the Research Collaboratory Structural Bioinformatics - Protein Data Bank (RCSB - PDB) repository and thus we predicted the first human UGT2B10 3D model by using multiple template homology modeling with MODELLER Version 9.2 and molecular docking calculations with AutoDock Vina Version 1.2 were implemented to quantify the estimated binding affinity between selected putative substrates or ligands and UGT2B10. Finally, we performed molecular dynamics simulations using GROMACS Version 5.1.4 to confirm the potential UGT2B10 ligands prioritized after molecular docking (exhibiting negative free binding energy).

Results

Four potential ligands for UGT2B10 namely acetaminophen, lorazepam, mycophenolic acid and voriconazole n-oxide intermediate were identified. Other metabolites of voriconazole satisfied the criteria of being possible ligands of UGT2B10. Except for bilirubin and 4-Hydroxy Voriconazole, all the ligands (particularly voriconazole and hydroxy voriconazole) are oriented in substrate binding site close to the co-factor UDP (mean ± SD; 0.72 ± 0.33 nm). Further in vitro screening of the putative ligands prioritized by in silico pipeline is warranted to understand the nature of the ligands either as inhibitors or substrates of UGT2B10.

Conclusions

These results may indicate the clinical and pharmacological relevance UGT2B10 in pediatric HSCT setting. With this systematic computational methodology, we provide a rational-, time-, and cost-effective way to identify and prioritize the interesting putative substrates or inhibitors of UGT2B10 for further testing in in vitro experiments.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12860-021-00402-5.

Feasibility of a Pragmatic PBPK Modeling Approach: Towards Model-Informed Dosing in Pediatric Clinical Care

BACKGROUND AND OBJECTIVE

More than half of all drugs are still prescribed off-label to children. Pharmacokinetic (PK) data are needed to support off-label dosing, however for many drugs such data are either sparse or not representative. Physiologically-based pharmacokinetic (PBPK) models are increasingly used to study PK and guide dosing decisions. Building compound models to study PK requires expertise and is time-consuming. Therefore, in this paper, we studied the feasibility of predicting pediatric exposure by pragmatically combining existing compound models, developed e.g. for studies in adults, with a pediatric and preterm physiology model.

METHODS

Seven drugs, with various PK characteristics, were selected (meropenem, ceftazidime, azithromycin, propofol, midazolam, lorazepam, and caffeine) as a proof of concept. Simcyp^® v20 was used to predict exposure in adults, children, and (pre)term neonates, by combining an existing compound model with relevant virtual physiology models. Predictive performance was evaluated by calculating the ratios of predicted-to-observed PK parameter values (0.5- to 2-fold acceptance range) and by visual predictive checks with prediction error values.

RESULTS

Overall, model predicted PK in infants, children and adolescents capture clinical data. Confidence in PBPK model performance was therefore considered high. Predictive performance tends to decrease when predicting PK in the (pre)term neonatal population.

CONCLUSION

Pragmatic PBPK modeling in pediatrics, based on compound models verified with adult data, is feasible. A thorough understanding of the model assumptions and limitations is required, before model-informed doses can be recommended for clinical use.

Physiologically Based Pharmacokinetic Models Are Effective Support for Pediatric Drug Development

Pediatric drug development faces many difficulties. Traditionally, pediatric drug doses are simply calculated linearly based on the body weight, age, and body surface area of adults. Due to the ontogeny of children, this simple linear scaling may lead to drug overdose in pediatric patients. The physiologically based pharmacokinetic (PBPK) model, as a mathematical model, contributes to the research and development of pediatric drugs. An example of a PBPK model guiding drug dose selection in pediatrics has emerged and has been approved by the relevant regulatory agencies. In this review, we discuss the principle of the PBPK model, emphasize the necessity of establishing a pediatric PBPK model, introduce the absorption, distribution, metabolism, and excretion of the pediatric PBPK model, and understand the various applications and related prospects of the pediatric PBPK model.

Characterization of Hepatic UDP-Glucuronosyltransferase Enzyme Abundance-Activity Correlations and Population Variability Using a Proteomics Approach and Comparison with Cytochrome P450 Enzymes

The expression of ten major drug-metabolizing UDP-glucuronosyltransferase (UGT) enzymes in a panel of 130 human hepatic microsomal samples was measured using a liquid chromatography-tandem mass spectrometry-based approach. Simultaneously, ten cytochromes P450 and P450 reductase were also measured, and activity-expression relationships were assessed for comparison. The resulting data sets demonstrated that, with the exception of UGT2B17, 10th to 90th percentiles of UGT expression spanned 3- to 8-fold ranges. These ranges were small relative to ranges of reported mean UGT enzyme expression across different laboratories. We tested correlation of UGT expression with enzymatic activities using selective probe substrates. A high degree of abundance-activity correlation (Spearman's rank correlation coefficient > 0.6) was observed for UGT1As (1A1, 3, 4, 6) and cytochromes P450. In contrast, protein abundance and activity did not correlate strongly for UGT1A9 and UGT2B enzymes (2B4, 7, 10, 15, and 17). Protein abundance was strongly correlated for UGTs 2B7, 2B10, and 2B15. We suggest a number of factors may contribute to these differences including incomplete selectivity of probe substrates, correlated expression of these UGT2B isoforms, and the impact of splice and polymorphic variants on the peptides used in proteomics analysis, and exemplify this in the case of UGT2B10. Extensive correlation analyses identified important criteria for validating the fidelity of proteomics and enzymatic activity approaches for assessing UGT variability, population differences, and ontogenetic changes. SIGNIFICANCE STATEMENT: Protein expression data allow detailed assessment of interindividual variability and enzyme ontogeny. This study has observed that expression and enzyme activity are well correlated for hepatic UGT1A enzymes and cytochromes P450. However, for the UGT2B family, caution is advised when assuming correlation of expression and activity as is often done in physiologically based pharmacokinetic modeling. This can be due to incomplete probe substrate specificities, but may also be related to presence of inactive UGT protein materials and the effect of splicing variations.

Evaluation of Physiologically Based Pharmacokinetic Models to Predict the Absorption of BCS Class I Drugs in Different Pediatric Age Groups

Age-related changes in many parameters affecting drug absorption remain poorly characterized. The objective of this study was to apply physiologically based pharmacokinetic (PBPK) models in pediatric patients to investigate the absorption and pharmacokinetics of 4 drugs belonging to the Biopharmaceutics Classification System (BCS) class I administered as oral liquid formulations. Pediatric PBPK models built with PK-Sim/MoBi were used to predict the pharmacokinetics of acetaminophen, emtricitabine, theophylline, and zolpidem in different pediatric populations. The model performance for predicting drug absorption and pharmacokinetics was assessed by comparing the predicted absorption profile with the deconvoluted dose fraction absorbed over time and predicted with observed plasma concentration-time profiles. Sensitivity analyses were performed to analyze the effects of changes in relevant input parameters on the model output. Overall, most pharmacokinetic parameters were predicted within a 2-fold error range. The absorption profiles were generally reasonably predicted, but relatively large differences were observed for acetaminophen. Sensitivity analyses showed that the predicted absorption profile was most sensitive to changes in the gastric emptying time (GET) and the specific intestinal permeability. The drug's solubility played only a minor role. These findings confirm that gastric emptying time, more than intestinal permeability or solubility, is a key factor affecting BCS class I drug absorption in children. As gastric emptying time is prolonged in the fed state, a better understanding of the interplay between food intake and gastric emptying time in children is needed, especially in the very young in whom the (semi)fed condition is the prevailing prandial state, and hence prolonged gastric emptying time seems more plausible than the fasting state.

Ontogeny of Hepatic Transporters and Drug-Metabolizing Enzymes in Humans and in Nonclinical Species

The liver represents a major eliminating and detoxifying organ, determining exposure to endogenous compounds, drugs, and other xenobiotics. Drug transporters (DTs) and drug-metabolizing enzymes (DMEs) are key determinants of disposition, efficacy, and toxicity of drugs. Changes in their mRNA and protein expression levels and associated functional activity between the perinatal period until adulthood impact drug disposition. However, high-resolution ontogeny profiles for hepatic DTs and DMEs in nonclinical species and humans are lacking. Meanwhile, increasing use of physiologically based pharmacokinetic (PBPK) models necessitates availability of underlying ontogeny profiles to reliably predict drug exposure in children. In addition, understanding of species similarities and differences in DT/DME ontogeny is crucial for selecting the most appropriate animal species when studying the impact of development on pharmacokinetics. Cross-species ontogeny mapping is also required for adequate translation of drug disposition data in developing nonclinical species to humans. This review presents a quantitative cross-species compilation of the ontogeny of DTs and DMEs relevant to hepatic drug disposition. A comprehensive literature search was conducted on PubMed Central: Tables and graphs (often after digitization) in original manuscripts were used to extract ontogeny data. Data from independent studies were standardized and normalized before being compiled in graphs and tables for further interpretation. New insights gained from these high-resolution ontogeny profiles will be indispensable to understand cross-species differences in maturation of hepatic DTs and DMEs. Integration of these ontogeny data into PBPK models will support improved predictions of pediatric hepatic drug disposition processes. SIGNIFICANCE STATEMENT: Hepatic drug transporters (DTs) and drug-metabolizing enzymes (DMEs) play pivotal roles in hepatic drug disposition. Developmental changes in expression levels and activities of these proteins drive age-dependent pharmacokinetics. This review compiles the currently available ontogeny profiles of DTs and DMEs expressed in livers of humans and nonclinical species, enabling robust interpretation of age-related changes in drug disposition and ultimately optimization of pediatric drug therapy.

Model-Informed Pediatric Dose Selection for Dapagliflozin by Incorporating Developmental Changes

This analysis reports a quantitative modeling and simulation approach for oral dapagliflozin, a primarily uridine diphosphate-glucuronosyltransferase (UGT)-metabolized human sodium-glucose cotransporter 2 selective inhibitor. A mechanistic dapagliflozin physiologically based pharmacokinetic (PBPK) model was developed using in vitro metabolism and clinical pharmacokinetic (PK) data and verified for context of use (e.g., exposure predictions in pediatric subjects aged 1 month to 18 years). Dapagliflozin exposure is challenging to predict in pediatric populations owing to differences in UGT1A9 ontogeny maturation and paucity of clinical PK data in younger age groups. Based on the exposure-response relationship of dapagliflozin, twofold acceptance criteria were applied between model-predicted and observed drug exposures and PK parameters (area under the curve and maximum drug concentration) in various scenarios, including monotherapy in healthy adults (single/multiple dose), monotherapy in hepatically or renally impaired patients, and drug-drug interactions with UGT1A9 modulators, such as mefenamic acid and rifampin. The PBPK model captured the observed exposure within twofold of the observed monotherapy data in adults and adolescents and in special population. As a guide to determining dosing regimens in pediatric studies, the verified PBPK model, along with UGT enzyme ontogeny maturation understanding, was used for predictions of dapagliflozin monotherapy exposures in pediatric subjects aged 1 month to 18 years that best matched exposure in adult patients with a 10-mg single dose of dapagliflozin.

Enzyme Kinetics of Uridine Diphosphate Glucuronosyltransferases (UGTs)

Glucuronidation, catalyzed by uridine diphosphate glucuronosyltransferases (UGTs), is an important process for the metabolism and clearance of many lipophilic chemicals, including drugs, environmental chemicals, and endogenous compounds. Glucuronidation is a bisubstrate reaction that requires the aglycone and the cofactor, UDP-GlcUA. Accumulating evidence suggests that the bisubstrate reaction follows a compulsory-order ternary mechanism. To simplify the kinetic modeling of glucuronidation reactions in vitro, UDP-GlcUA is usually added to incubations in large excess. Many factors have been shown to influence UGT activity and kinetics in vitro, and these must be accounted for during experimental design and data interpretation. While the assessment of drug-drug interactions resulting from UGT inhibition has been challenging in the past, the increasing availability of UGT enzyme-selective substrate and inhibitor "probes" provides the prospect for more reliable reaction phenotyping and assessment of drug-drug interaction potential. Although extrapolation of the in vitro intrinsic clearance of a glucuronidated drug often underpredicts in vivo clearance, careful selection of in vitro experimental conditions and inclusion of extrahepatic glucuronidation may improve the predictivity of in vitro-in vivo extrapolation. Physiologically based pharmacokinetic (PBPK) modeling has also shown to be of value for predicting PK of drugs eliminated by glucuronidation.

Infant exposure to dolutegravir through placental and breastmilk transfer: a population pharmacokinetic analysis of DolPHIN-1

BACKGROUND

Rapid reduction of HIV viral load is paramount to prevent peripartum transmission in women diagnosed late in pregnancy. We investigated dolutegravir population pharmacokinetics in maternal plasma, cord, breastmilk and infant plasma of DolPHIN-1 participants (NCT02245022) presenting with untreated HIV late in pregnancy (28-36 weeks gestation).

METHODS

Pregnant women from Uganda and South Africa were randomised (1:1) to daily dolutegravir (50 mg) or efavirenz-based therapy. Dolutegravir pharmacokinetic sampling (0-24 hours) was undertaken 14 days after treatment initiation and within 1-3 weeks of delivery, with matched maternal and cord samples at delivery. Mothers switched to efavirenz and maternal and infant plasma and breastmilk samples taken 24, 48 or 72 hours post-switch. Nonlinear mixed effects (NONMEM v. 7.4) was used to describe dolutegravir in all matrices and to evaluate covariates.

RESULTS

Twenty-eight women and 22 infants were included. Maternal dolutegravir was described by a two-compartment model linked to a fetal and breastmilk compartment. Cord and breastmilk to maternal plasma ratios were 1.279 (1.209-1.281) and 0.033 (0.021-0.050), respectively. Infant dolutegravir was described by breastmilk-to-infant and infant elimination rate constants. No covariate effects were observed. Predicted infant dolutegravir half-life and time to protein adjusted-IC90 (0.064 mg/L) for those above this threshold were 37.9 hours (22.1-63.5) and 108.9 hours [(18.6-129.6); 4.5 days (0.8-5.4); n=13].

CONCLUSIONS

Breastfeeding contributed relatively little to infant plasma exposures but a median of 4.5 days additional prophylaxis to some of the breastfed infants was observed following maternal dolutegravir cessation (3-15 days postpartum), which waned with time postpartum as transplacental dolutegravir cleared.

Population pharmacokinetics of propofol in neonates and infants: Gestational and postnatal age to determine clearance maturation

AIMS

Develop a population pharmacokinetic model describing propofol pharmacokinetics in (pre)term neonates and infants, that can be used for precision dosing (e.g. during target-controlled infusion) of propofol in this population.

METHODS

A nonlinear mixed effects pharmacokinetic analysis (Monolix 2018R2) was performed, based on a pooled study population in 107 (pre)term neonates and infants.

RESULTS

In total, 836 blood samples were collected from 66 (pre)term neonates and 41 infants originating from 3 studies. Body weight (BW) of the pooled study population was 3.050 (0.580-11.440) kg, postmenstrual age (PMA) was 36.56 (27.00-43.00) weeks and postnatal age (PNA) was 1.14 (0-104.00) weeks (median and min-max range). A 3-compartment structural model was identified and the effect of BW was modelled using fixed allometric exponents. Elimination clearance maturation was modelled accounting for the maturational effect on elimination clearance until birth (by gestational age [GA]) and postpartum (by PNA and GA). The extrapolated adult (70 kg) population propofol elimination clearance (1.64 L min^-1 , estimated relative standard error = 6.02%) is in line with estimates from previous population pharmacokinetic studies. Empirical scaling of BW on the central distribution volume in function of PNA improved the model fit.

CONCLUSIONS

It is recommended to describe elimination clearance maturation by GA and PNA instead of PMA on top of size effects when analyzing propofol pharmacokinetics in populations including preterm neonates. Changes in body composition in addition to weight changes or other physio-anatomical changes may explain the changes in central distribution volume. The developed model may serve as a prior for propofol dose finding and target-controlled infusion in (preterm) neonates.

Evidence-based strategies for the characterisation of human drug and chemical glucuronidation in vitro and UDP-glucuronosyltransferase reaction phenotyping

Enzymes of the UDP-glucuronosyltransferase (UGT) superfamily contribute to the elimination of drugs from almost all therapeutic classes. Awareness of the importance of glucuronidation as a drug clearance mechanism along with increased knowledge of the enzymology of drug and chemical metabolism has stimulated interest in the development and application of approaches for the characterisation of human drug glucuronidation in vitro, in particular reaction phenotyping (the fractional contribution of the individual UGT enzymes responsible for the glucuronidation of a given drug), assessment of metabolic stability, and UGT enzyme inhibition by drugs and other xenobiotics. In turn, this has permitted the implementation of in vitro - in vivo extrapolation approaches for the prediction of drug metabolic clearance, intestinal availability, and drug-drug interaction liability, all of which are of considerable importance in pre-clinical drug development. Indeed, regulatory agencies (FDA and EMA) require UGT reaction phenotyping for new chemical entities if glucuronidation accounts for ≥25% of total metabolism. In vitro studies are most commonly performed with recombinant UGT enzymes and human liver microsomes (HLM) as the enzyme sources. Despite the widespread use of in vitro approaches for the characterisation of drug and chemical glucuronidation by HLM and recombinant enzymes, evidence-based guidelines relating to experimental approaches are lacking. Here we present evidence-based strategies for the characterisation of drug and chemical glucuronidation in vitro, and for UGT reaction phenotyping. We anticipate that the strategies will inform practice, encourage development of standardised experimental procedures where feasible, and guide ongoing research in the field.

Characterization of the Ontogeny of Hepatic UDP-Glucuronosyltransferase Enzymes Based on Glucuronidation Activity Measured in Human Liver Microsomes

An understanding of the postnatal development of hepatic UDP-glucuronosyltransferase (UGT) enzymes is required for accurate prediction of the age-dependent changes in pharmacokinetics of many drugs used in children. However, the maturation rate of hepatic UGT isoforms remains a major knowledge gap. This study aimed to establish the age-associated changes in glucuronidation activity of 10 major hepatic UGT isoforms in humans, namely, UGT1A1, UGT1A3, UGT1A4, UGT1A6, UGT1A9, UGT2B4, UGT2B7, UGT2B10, UGT2B15, and UGT2B17. Human liver microsomes from pediatric and adult donors were incubated under optimized incubation conditions to assess the activity rates of hepatic UGT isoforms using a panel of 19 in vitro UGT probe substrates and clinically used drugs. Statistically strong correlations of glucuronidation activities allowed the ontogeny of UGT1A1, UGT1A4, UGT2B7, UGT2B10, and UGT2B15 to be established using multiple selective UGT substrates and matched human liver microsome samples. The postnatal development of hepatic UGTs is isoform-dependent using either individual or cross-correlated selective isoform substrates. Maximal adult activity was reached at different times ranging from within a month (UGT1A1, UGT2B4, UGT2B7, UGT2B10, and UGT2B15), during infancy (UGT1A3, UGT1A4, and UGT1A9), to adolescence (UGT1A6 and UGT2B17). This study provides an extensive characterization of the postnatal ontogeny profiles of hepatic UGT enzymes that are instrumental for predicting drug disposition via in vitro-in vivo extrapolation algorithms and verifying pharmacokinetic predictions against in vivo observations via pediatric physiologically based pharmacokinetic modeling in pediatric patients.

Predictive Pediatric Modeling and Simulation Using Ontogeny Information

Food and Drug Administration submissions of physiologically based pharmacokinetic (PBPK) modeling and simulation of small-molecule drugs document the relevance of pediatric drug development and, in particular, information on dosing strategies in children. The most relevant prerequisite for reliable PBPK-based translation of adult pharmacokinetics of a small molecule to children is knowledge of the drug-specific absorption, distribution, metabolism, and elimination (ADME) processes in adults together with existing information about ontogeny of ADME processes relevant for the drug. All mechanisms driving a drug's clearance are of specific importance. For other drug modalities, our knowledge of ADME processes and ontogeny is still limited. More research is required, for example, to understand why some therapeutic proteins show complex differences in pharmacokinetics between adults and children, whereas other proteins seem to follow simple allometric scaling rules. Ontogeny information originates from various sources, such as (semi)quantitative mRNA expression, in vitro activity data, and deconvolution of in vivo pharmacokinetic data. The workflow for pediatric predictions is well described in several articles documenting successful translation from adults to children. The technical hurdles for PBPK modeling are low. State-of-the-art PBPK modeling software tools provide integrated pediatric translation workflows. For example, PK-Sim and MoBi are freely available as fully transparent open-source software via Open Systems Pharmacology (OSP). With the latest 2019 software release, version 8.0, OSP even provides a fully integrated technical framework for the qualification (and requalification) of any specific intended PBPK use in line with Food and Drug Administration and European Medicines Agency PBPK guidance. Qualification packages for pediatric translation are available on the OSP platform.

Use of computational toxicology (CompTox) tools to predict in vivo toxicity for risk assessment

Computational Toxicology tools were used to predict toxicity for three pesticides: propyzamide (PZ), carbaryl (CB) and chlorpyrifos (CPF). The tools used included: a) ToxCast/Tox21 assays (AC₅₀ s μM: concentration 50% maximum activity); b) in vitro-to-in vivo extrapolation (IVIVE) using ToxCast/Tox21 AC₅₀s to predict administered equivalent doses (AED: mg/kg/d) to compare to known in vivo Lowest-Observed-Effect-Level (LOEL)/Benchmark Dose (BMD); c) high throughput toxicokinetics population based (HTTK-Pop) using AC₅₀s for endpoints associated with the mode of action (MOA) to predict age-adjusted AED for comparison with in vivo LOEL/BMDs. ToxCast/Tox21 active-hit-calls for each chemical were predictive of targets associated with each MOA, however, assays directly relevant to the MOAs for each chemical were limited. IVIVE AEDs were predictive of in vivo LOEL/BMD₁₀s for all three pesticides. HTTK-Pop was predictive of in vivo LOEL/BMD₁₀s for PZ and CPF but not for CB after human age adjustments 11-15 (PZ) and 6-10 (CB) or 6-10 and 11-20 (CPF) corresponding to treated rat ages (in vivo endpoints). The predictions of computational tools are useful for risk assessment to identify targets in chemical MOAs and to support in vivo endpoints. Data can also aid is decisions about the need for further studies.

Physiologically-Based Pharmacokinetic Predictions of the Effect of Curcumin on Metabolism of Imatinib and Bosutinib: In Vitro and In Vivo Disconnect

PURPOSE

This study aimed to investigate the potential pharmacokinetic interactions between curcumin, imatinib and bosutinib, combining In Vitro and in silico methods.

METHODS

In Vitro metabolism of imatinib and bosutinib were investigated in pooled human liver microsomes and recombinant CYP3A4 enzyme in the presence and absence of curcumin and curcumin glucuronide using an LC-MS/MS assay for N-desmethyl metabolites. A physiologically-based pharmacokinetic (PBPK) model for curcumin formulated as solid lipid nanoparticles (SLN) was constructed using In Vitro glucuronidation kinetics and published clinical pharmacokinetic data. The potential effects of curcumin coadministration on systemic exposures of imatinib and bosutinib were predicted in silico using PBPK simulations.

RESULTS

Curcumin demonstrated potent reversible inhibition of cytochrome P450 (CYP)3A4-mediated N-demethylation of imatinib and bosutinib and CYP2C8-mediated metabolism of imatinib with inhibitory constants (k_i,u) of ≤1.5 μmol. L^-1. A confirmatory In Vitro study with paclitaxel, the 6α-hydroxylation of which is exclusively mediated by CYP2C8, was consistent with a potent inhibition of this enzyme by curcumin. Curcumin glucuronide also inhibited both CYP enzymes In Vitro, albeit to a lesser extent than that of curcumin. PBPK model simulations predicted that at recommended dosing regimens of SLN curcumin, coadministration would result in an increase in systemic exposures of imatinib and bosutinib of up to only 10%.

CONCLUSION

A PBPK model for curcumin in a SLN formulation was successfully developed. Although curcumin possesses a strong In Vitro inhibitory activity towards CYP3A4 and CYP2C8 enzymes, its interactions with imatinib and bosutinib were unlikely to be of clinical importance due to curcumin's poor bioavailability.

Extrapolation for a pharmacokinetic model for acetaminophen from adults to neonates: A Latin Hypercube Sampling analysis

Physiological and drug-specific parameters need to be adjusted when extrapolating a pharmacokinetic (PK) model from adults to neonates, so as to reproduce the time profiles of the studied drug(s) consistent with clinical, in vivo data or in vitro cell line measurements. In this paper we present a parameter analysis method, i.e. the Latin Hypercube Sampling (LHS) method for an acetaminophen (APAP) PK model. The original model consists of two compartments (the blood and the urine) with Michaelis-Menten kinetic parameters determined for APAP and its metabolites. The physiological parameters are scaled through allometric laws from adults to neonates, and APAP-specific parameters are adjusted for enzymatic maturational changes. The LHS method is used to statistically investigate the interplay between these parameters. The results for the extrapolated APAP model are consistent with published APAP PK data in neonates. We found the sulphation clearance parameter played a crucial role in the neonatal PK model, but its influence was weakened if the volume of distribution parameters were included. We suggest that this kind of in silico experiment could be valuable as the first step in PK model extrapolation between different ages.

Optimization of Experimental Conditions of Automated Glucuronidation Assays in Human Liver Microsomes Using a Cocktail Approach and Ultra-High Performance Liquid Chromatography-Tandem Mass Spectrometry

UDP-glucuronosyltransferase (UGT)-mediated metabolism is possibly the most important conjugation reaction for marketed drugs. However, there are currently no generally accepted standard incubation conditions for UGT microsomal assays, and substantial differences in experimental design and methodology between laboratories hinder cross-study comparison of in vitro activities. This study aimed to define optimal experimental conditions to determine glucuronidation activity of multiple UGT isoforms simultaneously using human liver microsomes. Hepatic glucuronidation activities of UGT1A1, UGT1A3, UGT1A4, UGT1A6, UGT1A9, UGT2B4, UGT2B7, UGT2B10, UGT2B15, and UGT2B17 were determined using cocktail incubations of 10 UGT probe substrates. Buffer components and cosubstrates were assessed over a range of concentrations including magnesium chloride (MgCl₂; 0-10 mM) and uridine 5'-diphosphoglucuronic acid (UDPGA; 1-25 mM) with either Tris-HCl or potassium phosphate buffer (100 mM, pH 7.4). Greater microsomal glucuronidation activity by different hepatic UGT isoforms was obtained using 10 mM MgCl₂ and 5 mM UDPGA with 100 mM Tris-HCl buffer. The influence of bovine serum albumin (BSA; 0.1%-2% w/v) on glucuronidation activity was also assessed. Enzyme- and substrate-dependent effects of BSA were observed, resulting in decreased total activity of UGT1A1, UGT1A3, and UGT2B17 and increased total UGT1A9 and UGT2B7 activity. The inclusion of BSA did not significantly reduce the between-subject variability of UGT activity. Future in vitro UGT profiling studies under the proposed optimized experimental conditions would allow high-quality positive control data to be generated across laboratories, with effective control of a high degree of between-donor variability for UGT activity and for chemical optimization toward lower-clearance drug molecules in a pharmaceutical drug discovery setting.