Ontogeny of human hepatic and intestinal transporter gene expression during childhood: age matters

Expression of intestinal drug transporter proteins and metabolic enzymes in neonatal and pediatric patients

The development of pediatric oral drugs is hampered by a lack of predictive simulation tools. These tools, in turn, require data on the physiological variables that influence oral drug absorption, including the expression of drug transporter proteins (DTPs) and drug-metabolizing enzymes (DMEs) in the intestinal tract. The expression of hepatic DTPs and DMEs shows age-related changes, but there are few data on protein levels in the intestine of children. In this study, tissue was collected from different regions of the small and large intestine from neonates (i.e., surgically removed tissue) and from pediatric patients (i.e., gastroscopic duodenal biopsies). The protein expression of clinically relevant DTPs and DMEs was determined using a targeted mass spectrometry approach. The regional distribution of DTPs and DMEs was similar to adults. Most DTPs, with the exception of MRP3, MCT1, and OCT3, and all DMEs showed the highest protein expression in the proximal small intestine. Several proteins (i.e., P-gp, ASBT, CYP3A4, CYP3A5, CYP2C9, CYP2C19, and UGT1A1) showed an increase with age. Such increase appeared to be even more pronounced for DMEs. This exploratory study highlights the developmental changes in DTPs and DMEs in the intestinal tract of the pediatric population. Additional evaluation of protein function in this population would elucidate the implications of the presented changes in protein expression on absorption of orally administered drugs in neonates and pediatric patients.

Have We Neglected to Study Target-Site Drug Exposure in Children? A Systematic Review of the Literature

BACKGROUND AND OBJECTIVE

Drug dosing should ideally be based on the drug concentrations at the target site, which, for most drugs, corresponds to the tissue. The exact influence of growth and development on drug tissue distribution is unclear. This systematic review compiles the current knowledge on the tissue distribution of systemically applied drugs in children, with the aim to identify priorities in tissue pharmacokinetic (PK) research in this population.

METHODS

A systematic literature search was performed in the MEDLINE and Embase databases.

RESULTS

Forty-two relevant articles were identified, of which 71% investigated antibiotics, while drug classes from the other studies were anticancer drugs, antifungals, anthelmintics, sedatives, thyreostatics, immunomodulators, antiarrhythmics, and exon skipping therapy. The majority of studies (83%) applied tissue biopsy as the sampling technique. Tonsil and/or adenoid tissue was most frequently examined (70% of all included patients). The majority of studies had a small sample size (median 9, range 1-93), did not include the youngest age categories (neonates and infants), and were of low reporting quality. Due to the heterogeneous data from different study compounds, dosing schedules, populations, and target tissues, the possibility for comparison of PK data between studies was limited.

CONCLUSION

The influence of growth and development on drug tissue distribution continues to be a knowledge gap, due to the paucity of tissue PK data in children, especially in the younger age categories. Future research in this field should be encouraged as techniques to safely investigate drug tissue disposition in children are available.

Challenges of pediatric pharmacotherapy: A narrative review of pharmacokinetics, pharmacodynamics, and pharmacogenetics

PURPOSE

Personalized pharmacotherapy, including for the pediatric population, provides optimal treatment and has emerged as a major trend owing to advanced drug therapeutics and diversified drug selection. However, it is essential to understand the growth and developmental characteristics of this population to provide appropriate drug therapy. In recent years, clinical pharmacogenetics has accumulated knowledge in pediatric pharmacotherapy, and guidelines from professional organizations, such as the Clinical Pharmacogenetics Implementation Consortium, can be consulted to determine the efficacy of specific drugs and the risk of adverse effects. However, the existence of a large knowledge gap hinders the use of these findings in clinical practice.

METHODS

We provide a narrative review of the knowledge gaps in pharmacokinetics (PK) and pharmacodynamics (PD) in the pediatric population, focusing on the differences from the perspective of growth and developmental characteristics. In addition, we explored PK/PD in relation to pediatric clinical pharmacogenetics.

RESULTS

The lack of direct and indirect biomarkers for more accurate assessment of the effects of drug administration limits the current knowledge of PD. In addition, incorporating pharmacogenetic insights as pivotal covariates is indispensable in this comprehensive synthesis for precision therapy; therefore, we have provided recommendations regarding the current status and challenges of personalized pediatric pharmacotherapy. The integration of clinical pharmacogenetics with the health care system and institution of educational programs for health care providers is necessary for its safe and effective implementation. A comprehensive understanding of the physiological and genetic complexities of the pediatric population will facilitate the development of effective and personalized pharmacotherapeutic strategies.

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.

Quantification of drug metabolising enzymes and transporter proteins in the paediatric duodenum via LC-MS/MS proteomics using a QconCAT technique

Characterising the small intestine absorptive membrane is essential to enable prediction of the systemic exposure of oral formulations. In particular, the ontogeny of key intestinal Drug Metabolising Enzymes and Transporter (DMET) proteins involved in drug disposition needs to be elucidated to allow for accurate prediction of the PK profile of drugs in the paediatric cohort. Using pinch biopsies from the paediatric duodenum (n = 36; aged 11 months to 15 years), the abundance of 21 DMET proteins and two enterocyte markers were quantified via LC-MS/MS. An established LCMS nanoflow method was translated to enable analysis on a microflow LC system, and a new stable-isotope-labelled QconCAT standard developed to enable quantification of these proteins. Villin-1 was used to standardise abundancy values. The observed abundancies and ontogeny profiles, agreed with adult LC-MS/MS-based data, and historic paediatric data obtained via western blotting. A linear trend with age was observed for duodenal CYP3A4 and CES2 only. As this work quantified peptides on a pinch biopsy coupled with a microflow method, future studies using a wider population range are very feasible. Furthermore, this DMET ontogeny data can be used to inform paediatric PBPK modelling and to enhance the understanding of oral drug absorption and gut bioavailability in paediatric populations.

Development and Verification of a Japanese Pediatric Physiologically Based Pharmacokinetic Model with Emphasis on Drugs Eliminated by Cytochrome P450 or Renal Excretion

Physiologically based pharmacokinetic (PBPK) models are useful in bridging drug exposure in different ethnic groups, and there is increasing regulatory application of this approach in adults. Reported pediatric PBPK models tend to focus on the North European population, with few examples in other ethnic groups. This study describes the development and verification of a Japanese pediatric PBPK population. The development of the model was based on the existing North European pediatric population. Japanese systems and clinical data were collated from public databases and the literature, and the underlying demographics and equations were optimized so that physiological outputs represented the Japanese pediatric population. The model was tested using 14 different small molecule drugs, eliminated by a variety of pathways, including cytochrome P450 3A4 (CYP3A4) metabolism and renal excretion. Given the limitations of the clinical data, the overall performance of the model was good, with 44/62 predictions for PK parameters (area under the plasma drug concentration-time curve, AUC; maximum serum concentration, Cmax ; clearance, CL) being within 0.8- to 1.25-fold, 56/62 within 0.67- to 1.5-fold, and 61/62 within 0.5- to 2.0-fold of the observed values. Specific results for the 5 CYP3A4 substrates showed 20/31 cases were predicted within 0.8- to 1.25-fold, 27/31 within 0.67- to 1.5-fold, and all were within 0.5- to 2.0-fold of the observed values. Given the increased regulatory use of pediatric PBPK in drug development, expanding these models to other ethnic groups are important. Considering qualifying these models based on the context of use, there is a need to expand on the current research to include a larger range of drugs with different elimination pathways. Collaboration among academic, industry, model providers, and regulators will facilitate further development.

Significant Effects of Renal Function on Mycophenolic Acid Total Clearance in Pediatric Kidney Transplant Recipients with Population Pharmacokinetic Modeling

BACKGROUND AND OBJECTIVES

Mycophenolic acid (MPA) is an immunosuppressant commonly prescribed in pediatric kidney transplantation to prevent graft rejection. Large variabilities in MPA plasma exposures have been observed in this population, which could result in severe adverse effects. The majority of the MPA pharmacokinetic data have been reported in adult populations, whereas information in pediatric patients is still very limited. The objective of this study was to establish a novel, nonlinear mixed-effects model for MPA and investigate the clinical variables affecting MPA population pharmacokinetics in pediatric kidney transplant recipients.

METHODS

Data were collected retrospectively from pediatric kidney transplant patients (≤ 18 years when MPA concentrations were initially collected; on oral administration of mycophenolate mofetil) in Calgary, Alberta, Canada. Nonlinear mixed-effect modeling was conducted using stochastic approximation expectation-maximization in Monolix 2021R2 (Lixoft SAS, France) to determine population pharmacokinetic estimates, interindividual variabilities, and interoccasional variabilities. Covariate models were constructed using the Model Proposal function in Monolix in conjunction with a systematic stepwise inclusion/elimination protocol. The best model was selected based on objective function values, relative standard errors, goodness-of-fit plots, prediction-corrected visual predictive checks, and numerical predictive checks.

RESULTS

A total of 50 pediatric kidney transplant patients (25 female) with 219 MPA plasma concentration-time profiles were included. The average age (± standard deviation) and posttransplant time for the sample population were 12.8 ± 4.8 years and 762 ± 1160 days, respectively. The majority of study subjects (i.e., > 85% based on all occasions) were co-administered tacrolimus. A two-compartment, first-order absorption with lag time and linear elimination structural model with lognormal distributed proportional residual errors best described the MPA concentration-time data. The absorption rate constant (2.52 h-1 or 0.042 min-1), lag time (0.166 h or 9.96 min), volumes of distributions of the central (22.8 L) and peripheral (216 L) compartments, and intercompartment clearance (17.6 L h-1 or 0.293 L min-1) were consistent with literature values; whereas total MPA clearance (0.72 L h-1 or 0.012 L min-1) was relatively reduced, likely due to the general lack of cyclosporine interactions and the stabilized graft functions from significantly longer posttransplant time in our sample population. Of the clinical variables tested, only estimated glomerular filtration rate (eGFR) was identified a significant covariate affecting total MPA clearance with a positive, exponential relationship. The final population pharmacokinetic model was successfully evaluated/validated using a variety of complementary methods.

CONCLUSION

We have successfully constructed and validated a novel population pharmacokinetic model of MPA in pediatric kidney transplant patients. A positive, nonlinear relationship between eGFR and total MPA clearance identified in our model is likely attributed to multiple concurrent mechanisms, which warrant further systematic investigations.

Developmental Pharmacokinetics of Antibiotics Used in Neonatal ICU: Focus on Preterm Infants

Neonatal Infections are among the most common reasons for admission to the intensive care unit. Neonatal sepsis (NS) significantly contributes to mortality rates. Empiric antibiotic therapy of NS recommended by current international guidelines includes benzylpenicillin, ampicillin/amoxicillin, and aminoglycosides (gentamicin). The rise of antibacterial resistance precipitates the growth of the use of antibiotics of the Watch (second, third, and fourth generations of cephalosporines, carbapenems, macrolides, glycopeptides, rifamycins, fluoroquinolones) and Reserve groups (fifth generation of cephalosporines, oxazolidinones, lipoglycopeptides, fosfomycin), which are associated with a less clinical experience and higher risks of toxic reactions. A proper dosing regimen is essential for effective and safe antibiotic therapy, but its choice in neonates is complicated with high variability in the maturation of organ systems affecting drug absorption, distribution, metabolism, and excretion. Changes in antibiotic pharmacokinetic parameters result in altered efficacy and safety. Population pharmacokinetics can help to prognosis outcomes of antibiotic therapy, but it should be considered that the neonatal population is heterogeneous, and this heterogeneity is mainly determined by gestational and postnatal age. Preterm neonates are common in clinical practice, and due to the different physiology compared to the full terms, constitute a specific neonatal subpopulation. The objective of this review is to summarize the evidence about the developmental changes (specific for preterm and full-term infants, separately) of pharmacokinetic parameters of antibiotics used in neonatal intensive care units.

Cell-type dependent regulation of pluripotency and chromatin remodeling genes by hydralazine

BACKGROUND

The generation of induced pluripotent stem cells has opened the field of study for stem cell research, disease modeling and drug development. However, the epigenetic signatures present in somatic cells make cell reprogramming still an inefficient process. This epigenetic memory constitutes an obstacle in cellular reprogramming. Here, we report the effect of hydralazine (HYD) and valproic acid (VPA), two small molecules with proven epigenetic activity, on the expression of pluripotency genes in adult (aHF) and neonatal (nbHF) human fibroblasts.

METHODS

aHF and nbHF were treated with HYD and/or VPA, and viability and gene expression assays for OCT4, NANOG, c-MYC, KLF4, DNMT1, TET3, ARID1A and ARID2 by quantitative PCR were performed. aHF and nbHF were transfected with episomal plasmid bearing Yamanaka factors (OCT4, SOX2, KLF4 and c-MYC) and exposed to HYD and VPA to determine the reprogramming efficiency. Methylation sensitive restriction enzyme (MSRE) qPCR assays were performed on OCT4 and NANOG promoter regions. Immunofluorescence assays were carried out for pluripotency genes on iPSC derived from aHF and nbHF.

RESULTS

HYD upregulated the expression of OCT4 (2.5-fold) and NANOG (fourfold) genes but not c-Myc or KLF4 in aHF and had no significant effect on the expression of all these genes in nbHF. VPA upregulated the expression of NANOG (twofold) in aHF and c-MYC in nbHF, while it downregulated the expression of NANOG in nbHF. The combination of HYD and VPA canceled the OCT4 and NANOG overexpression induced by HYD in aHF, while it reinforced the effects of VPA on c-Myc expression in nbHF. The HYD-induced overexpression of OCT4 and NANOG in aHDF was not dependent on demethylation of gene promoters, and no changes in the reprogramming efficiency were observed in both cell populations despite the downregulation of epigenetic genes DNMT1, ARID1A, and ARID2 in nbHF.

CONCLUSIONS

Our data provide evidence that HYD regulates the expression of OCT4 and NANOG pluripotency genes as well as ARID1A and ARID2 genes, two members of the SWI/SNF chromatin remodeling complex family, in normal human dermal fibroblasts.

Current knowledge, challenges and innovations in developmental pharmacology: A combined conect4children Expert Group and European Society for Developmental, Perinatal and Paediatric Pharmacology White Paper

Developmental pharmacology describes the impact of maturation on drug disposition (pharmacokinetics, PK) and drug effects (pharmacodynamics, PD) throughout the paediatric age range. This paper, written by a multidisciplinary group of experts, summarizes current knowledge, and provides suggestions to pharmaceutical companies, regulatory agencies and academicians on how to incorporate the latest knowledge regarding developmental pharmacology and innovative techniques into neonatal and paediatric drug development. Biological aspects of drug absorption, distribution, metabolism and excretion throughout development are summarized. Although this area made enormous progress during the last two decades, remaining knowledge gaps were identified. Minimal risk and burden designs allow for optimally informative but minimally invasive PK sampling, while concomitant profiling of drug metabolites may provide additional insight in the unique PK behaviour in children. Furthermore, developmental PD needs to be considered during drug development, which is illustrated by disease- and/or target organ-specific examples. Identifying and testing PD targets and effects in special populations, and application of age- and/or population-specific assessment tools are discussed. Drug development plans also need to incorporate innovative techniques such as preclinical models to study therapeutic strategies, and shift from sequential enrolment of subgroups, to more rational designs. To stimulate appropriate research plans, illustrations of specific PK/PD-related as well as drug safety-related challenges during drug development are provided. The suggestions made in this joint paper of the Innovative Medicines Initiative conect4children Expert group on Developmental Pharmacology and the European Society for Developmental, Perinatal and Paediatric Pharmacology, should facilitate all those involved in drug development.

Attenuated Ochratoxin A Transporter Expression in a Mouse Model of Nonalcoholic Steatohepatitis Protects against Proximal Convoluted Tubule Toxicity

Ochratoxin A (OTA) is an abundant mycotoxin, yet the toxicological impact of its disposition is not well studied. OTA is an organic anion transporter (OAT) substrate primarily excreted in urine despite a long half-life and extensive protein binding. Altered renal transporter expression during disease, including nonalcoholic steatohepatitis (NASH), may influence response to OTA exposure, but the impact of NASH on OTA toxicokinetics, tissue distribution, and associated nephrotoxicity is unknown. By inducing NASH in fast food-dieted/thioacetamide-exposed mice, we evaluated the effect of NASH on a bolus OTA exposure (12.5 mg/kg by mouth) after 3 days. NASH mice presented with less gross toxicity (44% less body weight loss), and kidney and liver weights of NASH mice were 11% and 24% higher, respectively, than healthy mice. Organ and body weight changes coincided with reduced renal proximal tubule cells vacuolation, degeneration, and necrosis, though no OTA-induced hepatic lesions were found. OTA systemic exposure in NASH mice increased modestly from 5.65 ± 1.10 to 7.95 ± 0.61 mg*h/ml per kg BW, and renal excretion increased robustly from 5.55% ± 0.37% to 13.11% ± 3.10%, relative to healthy mice. Total urinary excretion of OTA increased from 24.41 ± 1.74 to 40.07 ± 9.19 µg in NASH mice, and kidney-bound OTA decreased by ∼30%. Renal OAT isoform expression (OAT1-5) in NASH mice decreased by ∼50% with reduced OTA uptake by proximal convoluted cells. These data suggest that NASH-induced OAT transporter reductions attenuate renal secretion and reabsorption of OTA, increasing OTA urinary excretion and reducing renal exposure, thereby reducing nephrotoxicity in NASH. SIGNIFICANCE STATEMENT: These data suggest a disease-mediated transporter mechanism of altered tissue-specific toxicity after mycotoxin exposure, despite minimal systemic changes to ochratoxin A (OTA) concentrations. Further studies are warranted to evaluate the clinical relevance of this functional model and the potential effect of human nonalcoholic steatohepatitis on OTA and other organic anion substrate toxicity.

A proof of concept using the Ussing chamber methodology to study pediatric intestinal drug transport and age-dependent differences in absorption

Little is known about the impact of age on the processes governing human intestinal drug absorption. The Ussing chamber is a system to study drug transport across tissue barriers, but it has not been used to study drug absorption processes in children. This study aimed to explore the feasibility of the Ussing chamber methodology to assess pediatric intestinal drug absorption. Furthermore, differences between intestinal drug transport processes of children and adults were explored as well as the possible impact of age. Fresh terminal ileal leftover tissues from both children and adults were collected during surgery and prepared for Ussing chamber experiments. Paracellular (enalaprilat), transcellular (propranolol), and carrier-mediated drug transport by MDR1 (talinolol) and BCRP (rosuvastatin) were determined with the Ussing chamber methodology. We calculated apparent permeability coefficients and efflux ratios and explored their relationship with postnatal age. The success rate for the Ussing chamber experiments, as determined by electrophysiological measurements, was similar between children (58%, N = 15, median age: 44 weeks; range 8 weeks to 17 years) and adults (67%, N = 13). Mean serosal to mucosal transport of talinolol by MDR1 and rosuvastatin by BCRP was higher in adult than in pediatric tissues (p = 0.0005 and p = 0.0091). In contrast, within our pediatric cohort, there was no clear correlation for efflux transport across different ages. In conclusion, the Ussing chamber is a suitable model to explore pediatric intestinal drug absorption and can be used to further elucidate ontogeny of individual intestinal pharmacokinetic processes like drug metabolism and transport.

Role of Drug Transporters in Elucidating Inter-Individual Variability in Pediatric Chemotherapy-Related Toxicities and Response

Pediatric cancer treatment has evolved significantly in recent decades. The implementation of risk stratification strategies and the selection of evidence-based chemotherapy combinations have improved survival outcomes. However, there is large interindividual variability in terms of chemotherapy-related toxicities and, sometimes, the response among this population. This variability is partly attributed to the functional variability of drug-metabolizing enzymes (DME) and drug transporters (DTS) involved in the process of absorption, distribution, metabolism and excretion (ADME). The DTS, being ubiquitous, affects drug disposition across membranes and has relevance in determining chemotherapy response in pediatric cancer patients. Among the factors affecting DTS function, ontogeny or maturation is important in the pediatric population. In this narrative review, we describe the role of drug uptake/efflux transporters in defining pediatric chemotherapy-treatment-related toxicities and responses. Developmental differences in DTS and the consequent implications are also briefly discussed for the most commonly used chemotherapeutic drugs in the pediatric population.

Important roles of transporters in the pharmacokinetics of anti-viral nucleoside/nucleotide analogs

INTRODUCTION

Nucleoside analogs are an important class of antiviral agents. Due to the high hydrophilicity and limited membrane permeability of antiviral nucleoside/nucleotide analogs (AVNAs), transporters play critical roles in AVNA pharmacokinetics. Understanding the properties of these transporters is important to accelerate translational research for AVNAs.

AREAS COVERED

The roles of key transporters in the pharmacokinetics of 25 approved AVNAs were reviewed. Clinically relevant information that can be explained by the modulation of transporter functions is also highlighted.

EXPERT OPINION

Although the roles of transporters in the intestinal absorption and renal excretion of AVNAs have been well identified, more research is warranted to understand their roles in the distribution of AVNAs, especially to immune privileged compartments where treatment of viral infection is challenging. P-gp, MRP4, BCRP, and nucleoside transporters have shown extensive impacts in the disposition of AVNAs. It is highly recommended that the role of transporters should be investigated during the development of novel AVNAs. Clinically, co-administered inhibitors and genetic polymorphism of transporters are the two most frequently reported factors altering AVNA pharmacokinetics. Physiopathology conditions also regulate transporter activities, while their effects on pharmacokinetics need further exploration. Pharmacokinetic models could be useful for elucidating these complicated factors in clinical settings.

Pharmacokinetics of free and total mycophenolic acid in paediatric and adult renal transplant recipients: Exploratory analysis of the effects of clinical factors and gene variants

Clinical and genetic influencing factors on free fraction of mycophenolic acid (MPA) have rarely been discussed. The present study investigated whether the clinical and genetic factors could explain the variability in the pharmacokinetics of free MPA (fMPA) and total MPA (tMPA) in Chinese paediatric and adult renal transplant recipients. Twenty-eight paediatric and 31 adult patients were enrolled, and the concentrations of tMPA and fMPA were determined at 0 h (predose) and 0.5, 1, 1.5, 2, 4, 5, 8, 9, 10 and 12 h after mycophenolate mofetil administration. Genetic polymorphisms of UGTs (rs671448, rs1042597, rs2741049, rs62298861, rs7439366, rs12233719) and ABCC2 (rs717620) were simultaneously determined. The clinical and genetic data were analysed and reported. tMPA and fMPA concentrations adjusted for dose per body weight were consistently higher in adults than in paediatric patients. In the paediatric group, only albumin and time after transplantation correlated significantly with the MPA-free fraction variation, which could explain 32.4% of the variability. Besides, ABCC2 polymorphism, albumin and time after transplantation correlated significantly with the MPA-free fraction variation in adults, which could explain 56.9% of the variability. The influencing factors in the paediatric group are different from those in adults, which may be due to age-related transporter expression.

Impact of pharmacogenetics on variability in exposure to oral vinorelbine among pediatric patients: a model-based population pharmacokinetic analysis

PURPOSE

Better understanding of pharmacokinetics of oral vinorelbine (VNR) in children would help predicting drug exposure and, beyond, clinical outcome. Here, we have characterized the population pharmacokinetics of oral VNR and studied the factors likely to explain the variability observed in VNR exposure among young patients.

DESIGN/METHODS

We collected blood samples from 36 patients (mean age 11.6 years) of the OVIMA multicentric phase II study in children with recurrent/progressive low-grade glioma. Patients received 60 mg/m² of oral VNR on days 1, 8, and 15 during the first 28-day treatment cycle and 80 mg/m², unless contraindicated, from cycle 2-12. Population pharmacokinetic analysis was performed using nonlinear mixed-effects modeling within the Monolix^® software. Fifty SNPs of pharmacokinetic-related genes were genotyped. The influence of demographic, biological, and pharmacogenetic covariates on pharmacokinetic parameters was investigated using a stepwise multivariate procedure.

RESULTS

A three-compartment model, with a delayed double zero-order absorption and a first-order elimination, best described VNR pharmacokinetics in children. Typical population estimates for the apparent central volume of distribution (V_c/F) and elimination rate constant were 803 L and 0.60 h^-1, respectively. Following covariate analysis, BSA, leukocytes count, and drug transport ABCB1-rs2032582 SNP showed a dramatic impact on V_c/F. Conversely, age and sex had no significant effect on VNR pharmacokinetics.

CONCLUSION

Beyond canonical BSA and leukocytes, ABCB1-rs2032582 polymorphism showed a meaningful impact on VNR systemic exposure. Simulations showed that the identified covariates could have an impact on both efficacy and toxicity outcomes. Thus, a personalized dosing strategy, using those covariates, could help to optimize the efficacy/toxicity balance of VNR in children.

P-glycoprotein: new insights into structure, physiological function, regulation and alterations in disease

The multidrug resistance phenomenon presents a major threat to the pharmaceutical industry. This resistance is a common occurrence in several diseases and is mediated by multidrug transporters that actively pump substances out of the cell and away from their target regions. The most well-known multidrug transporter is the P-glycoprotein transporter. The binding sites within P-glycoprotein can accommodate a variety of compounds with diverse structures. Hence, numerous drugs are P-glycoprotein substrates, with new ones being identified every day. For many years, the mechanisms of action of P-glycoprotein have been shrouded in mystery, and scientists have only recently been able to elucidate certain structural and functional aspects of this protein. Although P-glycoprotein is highly implicated in multidrug resistant diseases, this transporter also performs various physiological roles in the human body and is expressed in several tissues, including the brain, kidneys, liver, gastrointestinal tract, testis, and placenta. The expression levels of P-glycoprotein are regulated by different enzymes, inflammatory mediators and transcription factors; alterations in which can result in the generation of a disease phenotype. This review details the discovery, the recently proposed structure and the regulatory functions of P-glycoprotein, as well as the crucial role it plays in health and disease.

Neonatal Dubin-Johnson syndrome: biochemical parameters, characteristics, and genetic variants study

BACKGROUND

The clinical characteristics and gene mutation characteristics of children with Dubin-Johnson syndrome (DJS) need in-depth study.

METHODS

The clinical and genomic data of neonatal Dubin-Johnson syndrome (NDJS) and 155 cases with idiopathic cholestasis (IC) were analyzed from June 2016 to August 2020 RESULTS: ABCC2 gene variants were identified in eight patients, including one patient with homozygous variants and seven patients with compound heterozygous variants. A total of 13 different ABCC variants were detected in the NDJS patients, including three nonsense variants, six missense variants, three frameshift variants, and a splice site variant. The variant c.2443C > T (p.R815X), c.4237_4238insCT (p.H1414Lfs*17), c.960_961insGT (p.L322Cfs*3), c.4250delC (p.S1417Ffs*14), c.2224G > A (p.D742N), c.4020G > C (p.K1340N), and c.2439 + 5G > A were not reported in the Human Gene Variant Database. There was no significance in the sex, birth weight, and onset age between the NDJS and IC groups. Compared with the IC group, the NDJS group had significantly higher levels of total bilirubin (TB), but a significantly lower level of alanine transaminase and a ratio of direct bilirubin (DB) to TB. There is no significance in total bile acid, gamma-glutamyl-transpeptidase, albumin, or international normalized ratio between the two groups.

CONCLUSIONS

NDJS should be considered in prolonged neonatal intrahepatic cholestasis, especially in infants with normal or slightly elevated transaminase levels.

IMPACT

Explore the biochemical parameters, characteristics, and genetic profile of NDJS. By summarizing the characteristics of biochemical indicators, seven new mutation types of the ABCC2 gene were detected, which expanded the mutation spectrum of the ABCC2 gene. NDJS should be considered in prolonged neonatal intrahepatic cholestasis, especially in infants with normal or slightly elevated transaminase levels.

Review of paediatric gastrointestinal physiology relevant to the absorption of orally administered medicines

Despite much progress in regulations to improve paediatric drug development, there remains a significant need to develop better medications for children. For the design of oral dosage forms, a detailed understanding of the specific gastrointestinal (GI) conditions in children of different age categories and how they differ from GI conditions in adults is essential. Several review articles have been published addressing the ontogeny of GI characteristics, including luminal conditions in the GI tract of children. However, the data reported in most of these reviews are of limited quality because (1) information was cited from very old publications and sometimes low quality sources, (2) data gaps in the original data were filled with textbook knowledge, (3) data obtained on healthy and sick children were mixed, (4) average data obtained on groups of patients were mixed with data obtained on individual patients, and (5) results obtained using investigative techniques that may have altered the outcome of the respective studies were considered. Consequently, many of these reviews draw conclusions that may be incorrect. The aim of the present review was to provide a comprehensive and updated overview of the available original data on the ontogeny of GI luminal conditions relevant to oral drug absorption in the paediatric population. To this end, the PubMed and Web of Science metadatabases were searched for appropriate studies that examined age-related conditions in the oral cavity, esophagus, stomach, small intestine, and colon. Maturation was observed for several GI parameters, and corresponding data sets were identified for each paediatric age group. However, it also became clear that the ontogeny of several GI traits in the paediatric population is not yet known. The review article provides a robust and valuable data set for the development of paediatric in vitro and in silico biopharmaceutical tools to support the development of age-appropriate dosage forms. In addition, it provides important information on existing data gaps and should provide impetus for further systematic and well-designed in vivo studies on GI physiology in children of specific age groups in order to close existing knowledge gaps and to sustainably improve oral drug therapy in children.

Recent lessons learned from population pharmacokinetic studies of mycophenolic acid: physiological, genomic, and drug interactions leading to the prediction of drug effects

INTRODUCTION

Mycophenolic acid (MPA) is a widely used immunosuppressant in transplantation and autoimmune disease. Highly variable pharmacokinetics have been observed with MPA, but the exact mechanisms remain largely unknown.

AREAS COVERED

The current review provided a critical, comprehensive update of recently published population pharmacokinetic/dynamic models of MPA (n=16 papers identified from PubMed and Embase, inclusive from January 2017 to August 2021), with specific emphases on the intrinsic and extrinsic factors influencing the pharmacology of MPA. The significance of the identified covariates, potential mechanisms, and comparisons to historical literature have been provided.

EXPERT OPINION

While select covariates affecting the population pharmacokinetics of MPA are consistently observed and mechanistically supported, some variables have not been regularly reported and/or lacked mechanistic explanation. Very few pharmacodynamic models were available, pointing to the need to extrapolate pharmacokinetic findings. Ideal models of MPA should consist of: i) utilizing optimal sampling points to allow the characterizations of absorption, re-absorption, and elimination phases; ii) characterizing unbound/total MPA, MPA metabolites, plasma/urinary concentrations, and genetic polymorphisms to facilitate mechanistic interpretations; and iii) incorporating actual outcomes and pharmacodynamic data to establish clinical relevance. We anticipate the field will continue to expand in the next 5 to 10 years.

Application of proteomic data in the translation of in vitro observations to associated clinical outcomes

Translation of information on drug exposure and effect is facilitated by in silico models that enable extrapolation of in vitro measurements to in vivo clinical outcomes. These models integrate drug-specific data with information describing physiological processes and pathological changes, including alterations to proteins involved in drug absorption, distribution and elimination. Over the past 15 years, quantitative proteomics has contributed a wealth of protein expression data, which are currently used for a variety of systems pharmacology applications, as a complement or a surrogate for activity of the corresponding proteins. In this review, we explore current and emerging applications of targeted and global (untargeted) proteomics in translational pharmacology as well as strategies for improved integration into model-based drug development.

Ontogeny of Small Intestinal Drug Transporters and Metabolizing Enzymes Based on Targeted Quantitative Proteomics

Most drugs are administered to children orally. An information gap remains on the protein abundance of small intestinal drug-metabolizing enzymes (DMEs) and drug transporters (DTs) across the pediatric age range, which hinders precision dosing in children. To explore age-related differences in DMEs and DTs, surgical leftover intestinal tissues from pediatric and adult jejunum and ileum were collected and analyzed by targeted quantitative proteomics for apical sodium-bile acid transporter, breast cancer resistance protein (BCRP), monocarboxylate transporter 1 (MCT1), multidrug resistance protein 1 (MDR1), multidrug resistance-associated protein (MRP) 2, MRP3, organic anion-transporting polypeptide 2B1, organic cation transporter 1, peptide transporter 1 (PEPT1), CYP2C19, CYP3A4, CYP3A5, UDP glucuronosyltransferase (UGT) 1A1, UGT1A10, and UGT2B7. Samples from 58 children (48 ileums, 10 jejunums, age range: 8 weeks to 17 years) and 16 adults (8 ileums, 8 jejunums) were analyzed. When comparing age groups, BCRP, MDR1, PEPT1, and UGT1A1 abundance was significantly higher in adult ileum as compared with the pediatric ileum. Jejunal BCRP, MRP2, UGT1A1, and CYP3A4 abundance was higher in the adults compared with children 0-2 years of age. Examining the data on a continuous age scale showed that PEPT1 and UGT1A1 abundance was significantly higher, whereas MCT1 and UGT2B7 abundance was lower in adult ileum as compared with the pediatric ileum. Our data contribute to the deeper understanding of the ontogeny of small intestinal drug-metabolizing enzymes and drug transporters and shows DME-, DT-, and intestinal location-specific, age-related changes. SIGNIFICANCE STATEMENT: This is the first study that describes the ontogeny of small intestinal DTs and DMEs in human using liquid chromatography with tandem mass spectrometry-based targeted quantitative proteomics. The current analysis provides a detailed picture about the maturation of DT and DME abundances in the human jejunum and ileum. The presented results supply age-related DT and DME abundance data for building more accurate PBPK models that serve to support safer and more efficient drug dosing regimens for the pediatric population.

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.

Bile formation and secretion: An update

Bile formation is a fundamental physiological process that is vital to the survival of all vertebrates. However, little was known about the mechanisms of this secretion until after World War II. Initial studies involved classic physiologic studies in animal models and humans, which progressed to include studies in isolated cells and membrane vesicles. The advent of molecular biology then led to the identification of specific transport systems that are the determinants of this secretion. Progress in this field was reviewed in the American Physiologic Society's series on "Comprehensive Physiology" in 2013. Herein, we provide an in-depth update of progress since that time.

Impact of gastrointestinal tract variability on oral drug absorption and pharmacokinetics: An UNGAP review

The absorption of oral drugs is frequently plagued by significant variability with potentially serious therapeutic consequences. The source of variability can be traced back to interindividual variability in physiology, differences in special populations (age- and disease-dependent), drug and formulation properties, or food-drug interactions. Clinical evidence for the impact of some of these factors on drug pharmacokinetic variability is mounting: e.g. gastric pH and emptying time, small intestinal fluid properties, differences in pediatrics and the elderly, and surgical changes in gastrointestinal anatomy. However, the link of colonic factors variability (transit time, fluid composition, microbiome), sex differences (male vs. female) and gut-related diseases (chronic constipation, anorexia and cachexia) to drug absorption variability has not been firmly established yet. At the same time, a way to decrease oral drug pharmacokinetic variability is provided by the pharmaceutical industry: clinical evidence suggests that formulation approaches employed during drug development can decrease the variability in oral exposure. This review outlines the main drivers of oral drug exposure variability and potential approaches to overcome them, while highlighting existing knowledge gaps and guiding future studies in this area.

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.

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.

Population pharmacokinetics of mycophenolic acid in paediatric patients

Mycophenolic acid (MPA) is widely used in paediatric kidney transplant patients and sometimes prescribed for additional indications. Population pharmacokinetic or pharmacodynamic modelling has been frequently used to characterize the fixed, random and covariate effects of MPA in adult patients. However, MPA population pharmacokinetic data in the paediatric population have not been systematically summarized. The objective of this narrative review was to provide an up-to-date critique of currently available paediatric MPA population pharmacokinetic models, with emphases on modelling techniques, pharmacological findings and clinical relevance. PubMed and EMBASE were searched from inception of database to May 2020, where a total of 11 studies have been identified representing kidney transplant (n = 4), liver transplant (n = 1), haematopoietic stem cell transplant (n = 1), idiopathic nephrotic syndrome (n = 2), systemic lupus erythematosus (n = 2), and a combined population consisted of kidney, liver and haematopoietic stem cell transplant patients (n = 1). Critical analyses were provided in the context of MPA absorption, distribution, metabolism, excretion and bioavailability in this paediatric database. Comparisons to adult patients were also provided. With respect to clinical utility, Bayesian estimation models (n = 6) with acceptable accuracy and precision for MPA exposure determination have also been identified and systematically evaluated. Overall, our analyses have identified unique features of MPA clinical pharmacology in the paediatric population, while recognizing several gaps that still warrant further investigations. This review can be used by pharmacologists and clinicians for improving MPA pharmacokinetic-pharmacodynamic modelling and patient care.

Opioid Treatment for Neonatal Opioid Withdrawal Syndrome: Current Challenges and Future Approaches

Chronic intrauterine exposure to psychoactive drugs often results in neonatal opioid withdrawal syndrome (NOWS). When nonpharmacologic measures are insufficient in controlling NOWS, morphine, methadone, and buprenorphine are first-line medications commonly used to treat infants with NOWS because of in utero exposure to opioids. Research suggests that buprenorphine may be the leading drug therapy used to treat NOWS when compared with morphine and methadone. Currently, there are no consensus or standardized treatment guidelines for medications prescribed for NOWS. Opioids used to treat NOWS exhibit large interpatient variability in pharmacokinetics (PK) and pharmacodynamic (PD) response in neonates. Organ systems undergo rapid maturation after birth that may alter drug disposition and exposure for any given dose during development. Data regarding the PK and PD of opioids in neonates are sparse. Pharmacometric methods such as physiologically based pharmacokinetic and population pharmacokinetic modeling can be used to explore factors predictive of some of the variability associated with the PK/PD of opioids in newborns. This review discusses the utility of pharmacometric techniques for enhancing precision dosing in infants requiring opioid treatment for NOWS. Applying these approaches may contribute to optimizing the outcome by reducing cumulative drug exposure, mitigating adverse drug effects, and reducing the burden of NOWS in neonates.

Mechanism of bilirubin elimination in urine: insights and prospects for neonatal jaundice

Despite a century of research, bilirubin metabolism and the transport mechanisms responsible for homeostasis of bilirubin in serum remain controversial. Emerging evidence on the hepatic membrane transporters and inherited disorders of bilirubin metabolism have contributed to a greater understanding of the various steps involved in bilirubin homeostasis and its associated excretory pathways. We discuss these recent research findings on hepatic membrane transporters and evaluate their significance on the newborn bilirubin metabolism and excretion. New insights gained speculate that a proportion of conjugated bilirubin is excreted via the renal system, as an alternative to the intestinal excretion, even in normal physiological jaundice with no associated pathological concerns. Finally, this paper discusses the clinical relevance of targeting the altered renal excretory pathway, as bilirubin in urine may hold diagnostic importance in screening for neonatal jaundice.

Pediatric Pharmacokinetics and Dose Predictions: A Report of a Satellite Meeting to the 10th Juvenile Toxicity Symposium

On April 24, 2019, a symposium on Pediatric Pharmacokinetics and Dose Predictions was held as a satellite meeting to the 10th Juvenile Toxicity Symposium. This symposium brought together scientists from academia, industry, and clinical research organizations with the aim to update each other on the current knowledge on pediatric drug development. Through more knowledge on specific ontogeny profiles of drug metabolism and transporter proteins, integrated into physiologically-based pharmacokinetic (PBPK) models, we have gained a more integrated understanding of age-related differences in pharmacokinetics (PKs), Relevant examples were presented during the meeting. PBPK may be considered the gold standard for pediatric PK prediction, but still it is important to know that simpler methods, such as allometry, allometry combined with maturation function, functions based on the elimination pathway, or linear models, also perform well, depending on the age range or the mechanisms involved. Knowledge from different methods and information sources should be combined (e.g., microdosing can reveal early read-out of age-related differences in exposure), and such results can be a value to verify models. To further establish best practices for dose setting in pediatrics, more in vitro and in vivo research is needed on aspects such as age-related changes in the exposure-response relationship and the impact of disease on PK. New information coupled with the refining of model-based drug development approaches will allow faster targeting of intended age groups and allow more efficient design of pediatric clinical trials.

No population left behind: Improving paediatric drug safety using informatics and systems biology

Adverse drugs effects (ADEs) in children are common and may result in disability and death. The current paediatric drug safety landscape, including clinical trials, is limited as it rarely includes children and relies on extrapolation from adults. Children are not small adults but go through an evolutionarily conserved and physiologically dynamic process of growth and maturation. Novel quantitative approaches, integrating observations from clinical trials and drug safety databases with dynamic mechanisms, can be used to systematically identify ADEs unique to childhood. In this perspective, we discuss three critical research directions using systems biology methodologies and novel informatics to improve paediatric drug safety, namely child versus adult drug safety profiles, age-dependent drug toxicities and genetic susceptibility of ADEs across childhood. We argue that a data-driven framework that leverages observational data, biomedical knowledge and systems biology modelling will reveal previously unknown mechanisms of pediatric adverse drug events and lead to improved paediatric drug safety.

Mass spectrometry-based abundance atlas of ABC transporters in human liver, gut, kidney, brain and skin

ABC transporters (ATP-binding cassette transporter) traffic drugs and their metabolites across membranes, making ABC transporter expression levels a key factor regulating local drug concentrations in different tissues and individuals. Yet, quantification of ABC transporters remains challenging because they are large and low-abundance transmembrane proteins. Here, we analysed 200 samples of crude and membrane-enriched fractions from human liver, kidney, intestine, brain microvessels and skin, by label-free quantitative mass spectrometry. We identified 32 (out of 48) ABC transporters: ABCD3 was the most abundant in liver, whereas ABCA8, ABCB2/TAP1 and ABCE1 were detected in all tissues. Interestingly, this atlas unveiled that ABCB2/TAP1 may have TAP2-independent functions in the brain and that biliary atresia (BA) and control livers have quite different ABC transporter profiles. We propose that meaningful biological information can be derived from a direct comparison of these data sets.

Proof of Concept: First Pediatric [14 C]microtracer Study to Create Metabolite Profiles of Midazolam

Growth and development affect drug-metabolizing enzyme activity thus could alter the metabolic profile of a drug. Traditional studies to create metabolite profiles and study the routes of excretion are unethical in children due to the high radioactive burden. To overcome this challenge, we aimed to show the feasibility of an absorption, distribution, metabolism, and excretion (ADME) study using a [14 C]midazolam microtracer as proof of concept in children. Twelve stable, critically ill children received an oral [14 C]midazolam microtracer (20 ng/kg; 60 Bq/kg) while receiving intravenous therapeutic midazolam. Blood was sampled up to 24 hours after dosing. A time-averaged plasma pool per patient was prepared reflecting the mean area under the curve plasma level, and subsequently one pool for each age group (0-1 month, 1-6 months, 0.5-2 years, and 2-6 years). For each pool [14 C]levels were quantified by accelerator mass spectrometry, and metabolites identified by high resolution mass spectrometry. Urine and feces (n = 4) were collected up to 72 hours. The approach resulted in sufficient sensitivity to quantify individual metabolites in chromatograms. [14 C]1-OH-midazolam-glucuronide was most abundant in all but one age group, followed by unchanged [14 C]midazolam and [14 C]1-OH-midazolam. The small proportion of unspecified metabolites most probably includes [14 C]midazolam-glucuronide and [14 C]4-OH-midazolam. Excretion was mainly in urine; the total recovery in urine and feces was 77-94%. This first pediatric pilot study makes clear that using a [14 C]midazolam microtracer is feasible and safe to generate metabolite profiles and study recovery in children. This approach is promising for first-in-child studies to delineate age-related variation in drug metabolite profiles.

Incorporating Ontogeny in Physiologically Based Pharmacokinetic Modeling to Improve Pediatric Drug Development: What We Know About Developmental Changes in Membrane Transporters

Developmental changes in the biological processes involved in the disposition of drugs, such as membrane transporter expression and activity, may alter the drug exposure and clearance in pediatric patients. Physiologically based pharmacokinetic (PBPK) models take these age-dependent changes into account and may be used to predict drug exposure in children. As a result, this mechanistic-based tool has increasingly been applied to improve pediatric drug development. Under the Prescription Drug User Fee Act VI, the US Food and Drug Administration has committed to facilitate the advancement of PBPK modeling in the drug application review process. Yet, significant knowledge gaps on developmental biology still exist, which must be addressed to increase the confidence of prediction. Recently, more data on ontogeny of transporters have emerged and supplied a missing piece of the puzzle. This article highlights the recent findings on the ontogeny of transporters specifically in the intestine, liver, and kidney. It also provides a case study that illustrates the utility of incorporating this information in predicting drug exposure in children using a PBPK approach. Collaborative work has greatly improved the understanding of the interplay between developmental physiology and drug disposition. Such efforts will continue to be needed to address the remaining knowledge gaps to enhance the application of PBPK modeling in drug development for children.

Pharmacokinetics in children with chronic kidney disease

In children, the main causes of chronic kidney disease (CKD) are congenital diseases and glomerular disorders. CKD is associated with multiple physiological changes and may therefore influence various pharmacokinetic (PK) parameters. A well-known consequence of CKD on pharmacokinetics is a reduction in renal clearance due to a decrease in the glomerular filtration rate. The impact of renal impairment on pharmacokinetics is, however, not limited to a decreased elimination of drugs excreted by the kidney. In fact, renal dysfunction may lead to modifications in absorption, distribution, transport, and metabolism as well. Currently, insufficient evidence is available to guide dosing decisions on many commonly used drugs. Moreover, the impact of maturation on drug disposition and action should be taken into account when selecting and dosing drugs in the pediatric population. Clinicians should take PK changes into consideration when selecting and dosing drugs in pediatric CKD patients in order to avoid toxicity and increase efficiency of drugs in this population. The aim of this review is to summarize known PK changes in relation to CKD and to extrapolate available knowledge to the pediatric CKD population to provide guidance for clinical practice.

Ceftriaxone-associated biliary pseudolithiasis in children: do we know enough?

Ceftriaxone is an antibiotic agent frequently used in paediatric hospital practice for the treatment of severe bacterial infections. The use of this agent can result in cholelithiasis and/or biliary sludge, more commonly in children than in adults. This systematic review was aimed at analysing available literature concerning ceftriaxone-associated biliary pseudolithiasis in paediatric patients, with a special emphasis on the clinical aspects. A literature analysis was performed using Medline and Embase electronic databases (articles published in English up to December 2019), with the search terms and combinations as follows:'ceftriaxone', 'cholelithiasis', 'biliary sludge' 'gallstones' 'neonates' 'children' 'clinical aspects' 'management'. Several case reports, case series and prospective/retrospective studies have documented a relationship between ceftriaxone treatment and biliary pseudolithiasis in the paediatric population, even though literature data regarding neonates and infants are scarce. Ceftriaxone-associated biliary pseudolithiasis is dose-dependent and usually asymptomatic but, sometimes, it may present with abdominal pain, nausea and emesis. Abdominal ultrasonography should be performed when this complication is suspected. Generally, ceftriaxone-associated cholelithiasis resolves over a variable period of time (days to months) after cessation of therapy. Therefore, a conservative approach to this condition is advocated, but a prolonged follow-up may be necessary. A personalized assessment of factors predisposing to ceftriaxone-associated biliary pseudolithiasis before prescribing the drug can allow to minimize the risk of developing it, with significant advantages in terms of human and economic costs.

Developmental patterns in human blood-brain barrier and blood-cerebrospinal fluid barrier ABC drug transporter expression

When drugs exert their effects in the brain, linear extrapolation of doses from adults could be harmful for children as the blood-brain barrier (BBB) and blood-CSF barrier (BCSFB) function is still immature. More specifically, age-related variation in membrane transporters may impact brain disposition. As human data on brain transporter expression is scarce, age dependent [gestational age (GA), postnatal age (PNA), and postmenstrual age (PMA)] variation in immunohistochemical localization and staining intensity of the ABC transporters P-glycoprotein (Pgp), breast cancer resistance protein (BCRP), and multidrug resistance-associated proteins 1, 2, 4, and 5 (MRP1/2/4/5) was investigated. Post mortem brain cortical and ventricular tissue was derived from 23 fetuses (GA range 12.9-39 weeks), 17 neonates (GA range 24.6-41.3 weeks, PNA range 0.004-3.5 weeks), 8 children (PNA range 0.1-3 years), and 4 adults who died from a wide variety of underlying conditions. In brain cortical BBB, immunostaining increased with age for Pgp and BCRP, while in contrast, MRP1 and MRP2 staining intensity appeared higher in fetuses, neonates, and children, as compared to adults. BCSFB was positively stained for Pgp, MRP1, and MRP2 and appeared stable across age, while BCRP was not detected. MRP4 and MRP5 were not detected in BBB or BCSFB. In conclusion, human BBB and BCSFB ABC membrane transporters show brain location and transporter-specific maturation.

Alternative Splicing of the SLCO1B1 Gene: An Exploratory Analysis of Isoform Diversity in Pediatric Liver

The hepatic influx transporter OATP1B1 (SLCO1B1) plays an important role in the disposition of endogenous substrates and drugs prescribed to children. Alternative splicing increases the diversity of protein products from > 90% of human genes and may be triggered by developmental signals. As concentrations of several endogenous OATP1B1 substrates change during growth and development, with this exploratory study we investigated age-dependent alternative splicing of SLCO1B1 mRNA in 97 postmortem livers (fetus-adolescents). Twenty-seven splice variants were detected; 10 were confirmed by additional bioinformatic analyses and verified by quantitative polymerase chain reaction, and selected for detailed analysis based on relative abundance, association with age, and overlap with an adjacent gene. Two splice variants code for reference OATP1B1 protein, and eight code for truncated proteins. The expression of eight isoforms was associated with age. We conclude that alternative splicing of SLCO1B1 occurs frequently in children; although the functional consequences remain unknown, the data raise the possibility of a regulatory role for alternative splicing in mediating developmental changes in drug disposition.

Distinguishing Kawasaki Disease from Febrile Infectious Disease Using Gene Pair Signatures

Kawasaki disease (KD) is an acute systemic vasculitis of childhood with prolonged fever, and the diagnosis of KD is mainly based on clinical criteria, which is prone to misdiagnosis with other febrile infectious (FI) diseases. Currently, there remain no effective molecular markers for KD diagnosis. In this study, we aimed to use a relative-expression-based method k-TSP and resampling framework to identify robust gene pair signatures to distinguish KD from bacterial and virus febrile infectious diseases. Our study pool consisted of 808 childhood patients from several studies and assigned to three groups, namely, the discovery set (n = 224), validation set-1 (n = 197), and validation set-2 (n = 387). We had identified 60 biologically relevant gene pairs and developed a top-ranked gene pair classifier (TRGP) using the first seven signatures, with the area under the receiver-operating characteristic curves (AUROC) of 0.947 (95% CI, 0.918-0.976), a sensitivity of 0.936 (95% CI, 0.872-0.987), and a specificity of 0.774 (95% CI, 0.705-0.836) in the discovery set. In the validation set-1, the TRGP classifier distinguished KD from FI with AUROC of 0.955 (95% CI, 0.919-0.991), a sensitivity of 0.959 (95% CI, 0.925-0.986), and a specificity of 0.863 (95% CI, 0.764-0.961). In the validation set-2, the predictive performance of classification was with an AUROC of 0.796 (95% CI, 0.747-0.845), a sensitivity of 0.797 (95% CI, 0.720-0.864), and a specificity of 0.661 (95% CI, 0.606-0.717). Our study reveals that gene pair signatures are robust across diverse studies and can be utilized as objective biomarkers to distinguish KD from FI, helping to develop a fast, simple, and effective molecular approach to improve the diagnosis of KD.

Model-Informed Drug Development for Everolimus Dosing Selection in Pediatric Infant Patients

Everolimus is currently approved in Europe as an adjunctive therapy for patients aged ≥ 2 years with tuberous sclerosis complex (TSC)–associated treatment‐refractory partial‐onset seizures, based on the EXIST‐3 study (NCT01713946) results. As TSC‐associated seizures can also affect children aged between 6 months and 2 years, a modeling and simulation (M&S) approach was undertaken to extrapolate exposure (trough plasma concentration (C_min)) after a dose of 6 mg/m² and reduction in seizure frequency (RSF). A physiologically based pharmacokinetic model using Simcyp was developed to predict C_min in adult and pediatric patients, which was then used by a population pharmacodynamic model and a linear mixed effect model to predict short‐term and long‐term efficacy in adults (for validation) and in children, respectively. Based on the results of the M&S study, everolimus at the dose of 6 mg/m² is anticipated to be an efficacious treatment in children 6 months to 2 years of age (up to 77.8% RSF) with concentrations within the recommended target range.

Implementation of a Physiologically Based Pharmacokinetic Modeling Approach to Guide Optimal Dosing Regimens for Imatinib and Potential Drug Interactions in Paediatrics

Long-term use of imatinib is effective and well-tolerated in children with chronic myeloid leukaemia (CML) yet defining an optimal dosing regimen for imatinib in younger patients is a challenge. The potential interactions between imatinib and coadministered drugs in this “special” population also remains largely unexplored. This study implements a physiologically based pharmacokinetic (PBPK) modeling approach to investigate optimal dosing regimens and potential drug interactions with imatinib in the paediatric population. A PBPK model for imatinib was developed in the Simcyp Simulator (version 17) utilizing in silico, in vitro drug metabolism, and in vivo pharmacokinetic data and verified using an independent set of published clinical pharmacokinetic data. The model was then extrapolated to children and adolescents (aged 2–18 years) by incorporating developmental changes in organ size and maturation of drug-metabolising enzymes and plasma protein responsible for imatinib disposition. The PBPK model described imatinib pharmacokinetics in adult and paediatric populations and predicted drug interaction with carbamazepine, a cytochrome P450 (CYP)3A4 and 2C8 inducer, with a good accuracy (evaluated by visual inspections of the simulation results and predicted pharmacokinetic parameters that were within 1.25-fold of the clinically observed values). The PBPK simulation suggests that the optimal dosing regimen range for imatinib is 230–340 mg/m²/d in paediatrics, which is supported by the recommended initial dose for treatment of childhood CML. The simulations also highlighted that children and adults being treated with imatinib have similar vulnerability to CYP modulations. A PBPK model for imatinib was successfully developed with an excellent performance in predicting imatinib pharmacokinetics across age groups. This PBPK model is beneficial to guide optimal dosing regimens for imatinib and predict drug interactions with CYP modulators in the paediatric population.

WB-PBPK approach in predicting zidovudine pharmacokinetics in preterm neonates

Antiretroviral therapy has been the mainstay of treatment for neonates born to HIV infected mothers. Neonates born prematurely to HIV positive mothers are underdeveloped not only in anatomical terms but also in their physiological systems. Zidovudine, the first antiretroviral drug in clinical therapy for the treatment of HIV has been approved for use in preterm neonates both prophylactically and therapeutically. The present work describes the whole body physiologically based pharmacokinetic (WB-PBPK) model development for zidovudine in preterm neonates of varying gestational ages, to observe the pharmacokinetic behavior of the drug in this vulnerable group of the population. Along with the height, weight, post-natal, and gestational ages of the preterm neonates, metabolic enzymes CYP2A6, CYP2C8, etc. were incorporated for each neonate. The composition of the different organs in terms of water and lipid components, blood flow rates, etc. were specified during simulations according to the gestational ages of these neonates. The following PK parameters were estimated for preterm neonates using simulated plasma profiles: AUC 2686.41 ± 123.49 μmol min/L, C_max 6.46 ± 0.74 μmol/L, half-life 8.98 ± 2.36 hr, mean residence time 12.23 ± 3.43 hr, and total plasma clearance 1.48 ± 0.19 ml/min/kg in comparison with the observed PK parameters of a clinical study by Mirochknic et al. in preterm neonates with AUC 2020.04 μmol/min/L, C_max 6.10 μmol/L, and total plasma clearance 1.62 ml/min/kg. PBPK simulations provide an opportunity to visualize the possible impact of physiological maturity levels at varying gestational ages on the pharmacokinetic behavior of zidovudine in preterm neonates.

Clinical and molecular characterization of four patients with NTCP deficiency from two unrelated families harboring the novel SLC10A1 variant c.595A>C (p.Ser199Arg)

Sodium taurocholate cotransporting polypeptide (NTCP), a carrier protein encoded by solute carrier family 10 member 1 (SLC10A1), is expressed in the basolateral membrane of hepatocytes, where it is responsible for the uptake of bile acids from plasma into hepatocytes. The first patient with NTCP deficiency was described in 2015. A limited number of such patients have been reported in the literature and their genotypic and phenotypic features require further investigation. The current study investigated 4 patients with NTCP deficiency from two unrelated families. The patients were subjected to SLC10A1 genetic analysis and it was revealed that all patients were compound heterozygous for the c.800C>T (p.Ser267Phe) and c.595A>C (p.Ser199Arg) SLC10A1 variants. To the best of the authors' knowledge, the latter variant had not been previously reported. Further analysis in 50 healthy individuals did not identify carriers. The c.595A>C (p.Ser199Arg) variant exhibited co-segregation with hypercholanemia and exhibited a relatively conserved amino acid when compared with homologous peptides. Moreover, SWISS-MODEL prediction revealed that the mutation affected the conformation of the NTCP molecule. The 4 patients demonstrated varying degrees of hypercholanemia while a downward trend in the plasma levels of total bile acids (TBA) in 2 pediatric patients and occasionally normal TBA level in an adult case were observed. The results indicated an autosomal recessive trait for NTCP deficiency, supported the primary role of NTCP in the uptake of bile acids from plasma and suggested that hepatic uptake of bile acids may occur by means other than NTCP uptake. Moreover, the novel missense variant c.595A>C(p.Ser199Arg) enriched the SLC10A1 mutation spectrum and may serve as a new genetic marker for the molecular diagnosis and genetic counseling of NTCP deficiency.

Population Pharmacokinetic Study of Prophylactic Fluconazole in Preterm Infants for Prevention of Invasive Candidiasis

Fluconazole is an antifungal agent with reported evidence for its prophylactic effect against systemic fungal infection in preterm infants. The aim of this study was to build a population pharmacokinetic model to evaluate the pharmacokinetic characteristics of intravenous and oral fluconazole in preterm infants with the current prophylactic fluconazole dosing regimen. A pharmacokinetic model was developed using 301 fluconazole concentrations from 75 preterm infants with a baseline body weight (WT) ranging from 0.5 to 1.5 kg and an estimated glomerular filtration rate (eGFR) ranging from 12.9 to 58.5 ml/min/1.73 m² Eligible infants received an intravenous or oral dose of 3 mg/kg of body weight of fluconazole, twice weekly with a ≥72-h dose interval, for 4 weeks. The model was qualified with basic goodness-of-fit diagnostics, visual predictive checks, and bootstrapping. The fluconazole pharmacokinetics was well described with a one-compartment linear model with a proportional residual error. The population clearance (CL) and volume of distribution (V) were derived as 0.0197 × (WT/1.00)^0.746 × (eGFR/25.0)^0.463 × exp(η) and 1.04 × WT × exp(η), respectively. Such covariate analyses augment the awareness of the need for personalized dosing in preterm infants. (This study has been registered at ClinicalTrials.gov under identifier NCT01683760).