Determination of Human Hepatic CYP2C8 and CYP1A2 Age-Dependent Expression to Support Human Health Risk Assessment for Early Ages

Advancing understanding of human variability through toxicokinetic modeling, in vitro-in vivo extrapolation, and new approach methodologies

The merging of physiology and toxicokinetics, or pharmacokinetics, with computational modeling to characterize dosimetry has led to major advances for both the chemical and pharmaceutical research arenas. Driven by the mutual need to estimate internal exposures where in vivo data generation was simply not possible, the application of toxicokinetic modeling has grown exponentially in the past 30 years. In toxicology the need has been the derivation of quantitative estimates of toxicokinetic and toxicodynamic variability to evaluate the suitability of the tenfold uncertainty factor employed in risk assessment decision-making. Consideration of a host of physiologic, ontogenetic, genetic, and exposure factors are all required for comprehensive characterization. Fortunately, the underlying framework of physiologically based toxicokinetic models can accommodate these inputs, in addition to being amenable to capturing time-varying dynamics. Meanwhile, international interest in advancing new approach methodologies has fueled the generation of in vitro toxicity and toxicokinetic data that can be applied in in vitro-in vivo extrapolation approaches to provide human-specific risk-based information for historically data-poor chemicals. This review will provide a brief introduction to the structure and evolution of toxicokinetic and physiologically based toxicokinetic models as they advanced to incorporate variability and a wide range of complex exposure scenarios. This will be followed by a state of the science update describing current and emerging experimental and modeling strategies for population and life-stage variability, including the increasing application of in vitro-in vivo extrapolation with physiologically based toxicokinetic models in pharmaceutical and chemical safety research. The review will conclude with case study examples demonstrating novel applications of physiologically based toxicokinetic modeling and an update on its applications for regulatory decision-making. Physiologically based toxicokinetic modeling provides a sound framework for variability evaluation in chemical risk assessment.

Caffeine: The Story beyond Oxygen-Induced Lung and Brain Injury in Neonatal Animal Models-A Narrative Review

Caffeine is one of the most commonly used drugs in intensive care to stimulate the respiratory control mechanisms of very preterm infants. Respiratory instability, due to the degree of immaturity at birth, results in apnea of prematurity (AOP), hyperoxic, hypoxic, and intermittent hypoxic episodes. Oxidative stress cannot be avoided as a direct reaction and leads to neurological developmental deficits and even a higher prevalence of respiratory diseases in the further development of premature infants. Due to the proven antioxidant effect of caffeine in early use, largely protective effects on clinical outcomes can be observed. This is also impressively observed in experimental studies of caffeine application in oxidative stress-adapted rodent models of damage to the developing brain and lungs. However, caffeine shows undesirable effects outside these oxygen toxicity injury models. This review shows the effects of caffeine in hyperoxic, hypoxic/hypoxic-ischemic, and intermittent hypoxic rodent injury models, but also the negative effects on the rodent organism when caffeine is administered without exogenous oxidative stress. The narrative analysis of caffeine benefits in cerebral and pulmonary preterm infant models supports protective caffeine use but should be given critical consideration when considering caffeine treatment beyond the recommended corrected gestational age.

A physiologically-based pharmacokinetic modeling approach for dosing amiodarone in children on ECMO

Extracorporeal membrane oxygenation (ECMO) is a cardiopulmonary bypass device commonly used to treat cardiac arrest in children. The American Heart Association guidelines for cardiopulmonary resuscitation (CPR) and emergency cardiovascular care recommend using amiodarone as a first-line agent to treat ventricular arrhythmias in children with cardiac arrest. However, there are no dosing recommendations for amiodarone to treat ventricular arrhythmias in pediatric patients on ECMO. Amiodarone has a high propensity for adsorption to the ECMO components due to its physicochemical properties leading to altered pharmacokinetics (PK) in ECMO patients. The change in amiodarone PK due to interaction with ECMO components may result in a difference in optimal dosing in patients on ECMO when compared with non-ECMO patients. To address this clinical knowledge gap, a physiologically-based pharmacokinetic model of amiodarone was developed in adults and scaled to children, followed by the addition of an ECMO compartment. The pediatric model included ontogeny functions of cytochrome P450 (CYP450) enzyme maturation across various age groups. The ECMO compartment was parameterized using the adsorption data of amiodarone obtained from ex vivo studies. Model predictions captured observed concentrations of amiodarone in pediatric patients with ECMO well with an average fold error between 0.5 and 2. Model simulations support an amiodarone intravenous (i.v) bolus dose of 22 mg/kg (neonates), 13 mg/kg (infants), 8 mg/kg (children), and 6 mg/kg (adolescents). This PBPK modeling approach can be applied to explore the dosing of other drugs used in children on ECMO.

Activity-Based Protein Profiling to Probe Relationships between Cytochrome P450 Enzymes and Early-Age Metabolism of Two Polycyclic Aromatic Hydrocarbons (PAHs): Phenanthrene and Retene

A growing body of literature has linked early-life exposures to polycyclic aromatic hydrocarbons (PAH) with adverse neurodevelopmental effects. Once in the body, metabolism serves as a powerful mediator of PAH toxicity by bioactivating and detoxifying PAH metabolites. Since enzyme expression and activity vary considerably throughout human development, we evaluated infant metabolism of PAHs as a potential contributing factor to PAH susceptibility. We measured and compared rates of phenanthrene and retene (two primary PAH constituents of woodsmoke) metabolism in human hepatic microsomes from individuals ≤21 months of age to a pooled sample (n = 200) consisting primarily of adults. We used activity-based protein profiling (ABPP) to characterize cytochrome P450 enzymes (CYPs) in the same hepatic microsome samples. Once incubated in microsomes, phenanthrene demonstrated rapid depletion. Best-fit models for phenanthrene metabolism demonstrated either 1 or 2 phases, depending on the sample, indicating that multiple enzymes could metabolize phenanthrene. We observed no statistically significant differences in phenanthrene metabolism as a function of age, although samples from the youngest individuals had the slowest phenanthrene metabolism rates. We observed slower rates of retene metabolism compared with phenanthrene also in multiple phases. Rates of retene metabolism increased in an age-dependent manner until adult (pooled) metabolism rates were achieved at ∼12 months. ABPP identified 28 unique CYPs among all samples, and we observed lower amounts of active CYPs in individuals ≤21 months of age compared to the pooled sample. Phenanthrene metabolism correlated to CYPs 1A1, 1A2, 2C8, 4A22, 3A4, and 3A43 and retene metabolism correlated to CYPs 1A1, 1A2, and 2C8 measured by ABPP and vendor-supplied substrate marker activities. These results will aid efforts to determine human health risk and susceptibility to PAHs exposure during early life.

Single-cell metabolic profiling reveals subgroups of primary human hepatocytes with heterogeneous responses to drug challenge

BACKGROUND

Xenobiotics are primarily metabolized by hepatocytes in the liver, and primary human hepatocytes are the gold standard model for the assessment of drug efficacy, safety, and toxicity in the early phases of drug development. Recent advances in single-cell genomics demonstrate liver zonation and ploidy as main drivers of cellular heterogeneity. However, little is known about the impact of hepatocyte specialization on liver function upon metabolic challenge, including hepatic metabolism, detoxification, and protein synthesis.

RESULTS

Here, we investigate the metabolic capacity of individual human hepatocytes in vitro. We assess how chronic accumulation of lipids enhances cellular heterogeneity and impairs the metabolisms of drugs. Using a phenotyping five-probe cocktail, we identify four functional subgroups of hepatocytes responding differently to drug challenge and fatty acid accumulation. These four subgroups display differential gene expression profiles upon cocktail treatment and xenobiotic metabolism-related specialization. Notably, intracellular fat accumulation leads to increased transcriptional variability and diminishes the drug-related metabolic capacity of hepatocytes.

CONCLUSIONS

Our results demonstrate that, upon a metabolic challenge such as exposure to drugs or intracellular fat accumulation, hepatocyte subgroups display different and heterogeneous transcriptional responses.

Integration of Biorelevant Pediatric Dissolution Methodology into PBPK Modeling to Predict In Vivo Performance and Bioequivalence of Generic Drugs in Pediatric Populations: a Carbamazepine Case Study

This study investigated the impact of gastro-intestinal fluid volume and bile salt (BS) concentration on the dissolution of carbamazepine (CBZ) immediate release (IR) 100 mg tablets and to integrate these in vitro biorelevant dissolution profiles into physiologically based pharmacokinetic modelling (PBPK) in pediatric and adult populations to determine the biopredictive dissolution profile. Dissolution profiles of CBZ IR tablets (100 mg) were generated in 50-900 mL biorelevant adult fasted state simulated gastric and intestinal fluid (Ad-FaSSGF and Ad-FaSSIF), also in three alternative compositions of biorelevant pediatric FaSSGF and FaSSIF medias at 200 mL. This study found that CBZ dissolution was poorly sensitive to changes in the composition of the biorelevant media, where dissimilar dissolution (F2 = 46.2) was only observed when the BS concentration was changed from 3000 to 89 μM (Ad-FaSSIF vs Ped-FaSSIF 50% 14 BS). PBPK modeling demonstrated the most predictive dissolution volume and media composition to forecast the PK was 500 mL of Ad-FaSSGF/Ad-FaSSIF media for adults and 200 mL Ped-FaSSGF/FaSSIF media for pediatrics. A virtual bioequivalence simulation was conducted by using Ad-FaSSGF and/or Ad-FaSSIF 500 mL or Ped-FaSSGF and/or Ped-FaSSIF 200 mL dissolution data for CBZ 100 mg (reference and generic test) IR product. The CBZ PBPK models showed bioequivalence of the product. This study demonstrates that the integration of biorelevant dissolution data can predict the PK profile of a poorly soluble drug in both populations. Further work using more pediatric drug products is needed to verify biorelevant dissolution data to predict the in vivo performance in pediatrics.

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.

Therapeutic drug monitoring of caffeine and its primary metabolites in plasma using LC-ESI-MS/MS for apnea of prematurity treatment: Evaluation of ultrapure water as a surrogate matrix

The growing evidence has endorsed the view that therapeutic drug monitoring (TDM) of caffeine for apnea of prematurity is helpful for dose tailoring when the therapeutic response is lacking or toxicity is suspected. However, the plasma without caffeine is difficult to obtain. Therefore, a method was developed and validated to measure caffeine and its three primary metabolites (paraxanthine, theobromine, and theophylline) using LC-ESI-MS/MS in human plasma and several surrogate matrices. The chromatographic separation of analytes was finally achieved on a Waters Symmetry C18 (4.6 × 75 mm, 3.5 μm) column. Several strategies were successfully applied to overcome the matrix effects: 1) appropriate dilution for sample cleanup; 2) a starting lower proportion of organic phase; 3) multiple individual stable-labeled isotopic internal standards. The parallelism between the authentic matrix and surrogate matrices was convincing. The recovery of the analytes in both human plasma and rat plasma was acceptable over the linear range (0.500 to 50.0 μg/mL for caffeine, and 0.0100 to 1.00 μg/mL for three metabolites). The method was successfully applied in 118 samples from 74 preterm infants with apnea of prematurity. The rat plasma or ultrapure water as surrogate matrix is worthy of being recommended for routine TDM of caffeine.

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

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

[Recent research on gene polymorphisms related to caffeine therapy in preterm infants with apnea of prematurity]

Apnea of prematurity (AOP) is one of the common diseases in preterm infants. The main cause of AOP is immature development of the respiratory control center. If AOP is not treated timely and effectively, it will lead to respiratory failure, hypoxic brain injury, and even death in severe cases. Caffeine is the first choice for the treatment of AOP, but its effectiveness varies in preterm infants. With the deepening of AOP research, more and more genetic factors have been confirmed to play important roles in the pathogenesis and treatment of AOP; in particular, the influence of single nucleotide polymorphism on the efficacy of caffeine has become a research hotspot in recent years. This article reviews the gene polymorphisms that affect the efficacy of caffeine, in order to provide a reference for individualized caffeine therapy. Citation.

Physiologically-Based Pharmacokinetic Modeling to Investigate the Effect of Maturation on Buprenorphine Pharmacokinetics in Newborns with Neonatal Opioid Withdrawal Syndrome

Neonatal opioid withdrawal syndrome (NOWS) is a major public health concern whose incidence has paralleled the opioid epidemic in the United States. Sublingual buprenorphine is an emerging treatment for NOWS, but given concerns about long-term adverse effects of perinatal opioid exposure, precision dosing of buprenorphine is needed. Buprenorphine pharmacokinetics (PK) in newborns, however, is highly variable. To evaluate underlying sources of PK variability, a neonatal physiologically-based pharmacokinetic (PBPK) model of sublingual buprenorphine was developed using Simcyp (version 19.1). The PBPK model included metabolism by cytochrome P450 (CYP) 3A4, CYP2C8, UDP-glucuronosyltransferase (UGT) 1A1, UGT1A3, UGT2B7, and UGT2B17, with additional biliary excretion. Maturation of metabolizing enzymes was incorporated, and default CYP2C8 and UGT2B7 ontogeny profiles were updated according to recent literature. A biliary clearance developmental profile was outlined using clinical data from neonates receiving sublingual buprenorphine as NOWS treatment. Extensive PBPK model validation in adults demonstrated good predictability, with geometric mean (95% confidence interval (CI)) predicted/observed ratios (P/O ratios) of area under the curve from zero to infinity (AUC0-∞ ), peak concentration (Cmax ), and time to reach peak concentration (Tmax ) equaling 1.00 (0.74-1.33), 1.04 (0.84-1.29), and 0.95 (0.72-1.26), respectively. In neonates, the geometric mean (95% CI) P/O ratio of whole blood concentrations was 0.75 (95% CI 0.64-0.87). PBPK modeling and simulation demonstrated that variability in biliary clearance, sublingual absorption, and CYP3A4 abundance are likely important drivers of buprenorphine PK variability in neonates. The PBPK model could be used to guide development of improved buprenorphine starting dose regimens for the treatment of NOWS.

Plasma Drug Concentration of Propranolol and Genetic Study in Chinese Han Patients With Infantile Haemangioma

BACKGROUND AND PURPOSE

This study was conducted to explore the plasma drug concentration of propranolol in Chinese Han patients with infantile haemangioma (IH) and the influencing factors, as well as the relationship among plasma drug concentrations of propranolol, β1-AR mutation and CYP2D6 188C>T, efficacy, and safety.

EXPERIMENTAL APPROACH

From January 2018 to April 2019, 140 patients with IH who were admitted to the hospital for oral propranolol and agreed to have their plasma concentration of propranolol tested, including 112 patients with β1-AR and CYP2D6 gene tested.

KEY RESULTS AND CONCLUSIONS AND IMPLICATIONS

The mean peak blood levels of propranolol, 4-hydroxypropranolol (4-OH-P), and N-deisopropylpropranolol (NDP) were 60.35 ± 37.90, 1.90 ± 2.37, and 0.24 ± 0.18 ng/ml, respectively. The mean trough blood levels of propranolol, 4-OH-P, and NDP were 24.98 ± 17.68, 0.45 ± 0.52, and 0.05±0.05 ng/ml, respectively. The higher the dose of propranolol, the higher the plasma concentration of propranolol (p = 0.031). The plasma concentration of propranolol was not related to the treatment efficacy.

Pharmacokinetic Evaluation of Eltrombopag in ITP Pediatric Patients

Background: Eltrombopag (EPAG) is an oral thrombopoietin receptor agonist, approved for refractory primary immune thrombocytopenia (ITP) in pediatric patients. In two pediatric RCTs, EPAG led to an improvement of platelet counts and a reduction in bleeding severity. However, a significant number of pediatric patients did not achieve the primary endpoints. We performed a pharmacokinetic evaluation of EPAG in pediatric patients with refractory ITP. Methods: Outpatients aged from 1 to 17 y, affected by refractory ITP to first-line treatment, were enrolled for a pharmacokinetic assessment. The analysis of drug plasma concentration was performed by the LC-MS/MS platform. Non-compartmental and statistical subgroup analyses were carried out using the R package ncappc. Results: Among 36 patients eligible for PK analysis, the median dose of EPAG given once daily was 50 mg. The EPAG peak occurs between 2 and 4 h with a population Cmax and AUC 0-24 geo-mean of 23, 38 μg/ml, and 275, 4 µg*h/mL, respectively. The pharmacokinetic profile of EPAG did not show a dose proportionality. Female patients showed a statistically significant increase of dose-normalized exposure parameters, increasing by 110 and 123% for Cmax and AUC 0-24, respectively, when compared to male patients. Patients aged 1-5 y showed values increased by more than 100% considering both exposure parameters, compared to older children. Furthermore, patients presenting complete response (83%), showed augmented EPAG exposure parameters compared to subjects with partial or no response. Conclusion: These data highlight the need to further explore the variability of EPAG exposure and its pharmacokinetic/pharmacodynamic profile in pediatric patients also in a real-life setting.

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.

Caffeine for the Pharmacological Treatment of Apnea of Prematurity in the NICU: Dose Selection Conundrum, Therapeutic Drug Monitoring and Genetic Factors

Caffeine citrate is the drug of choice for the pharmacological treatment of apnea of prematurity. Factors such as maturity and genetic variation contribute to the interindividual variability in the clinical response to caffeine therapy in preterm infants, making the optimal dose administered controversial. Moreover, the necessity for therapeutic drug monitoring (TDM) of caffeine is still worth discussing due to the need to achieve the desired target concentrations as well as concerns about the safety of higher doses. Therefore, we reviewed the pharmacokinetic profile of caffeine in preterm infants, evidence of the safety and efficacy of different doses of caffeine, therapeutic concentration ranges of caffeine and impact of genetic variability on caffeine therapy. Whereas the safety and efficacy of standard-dose caffeine have been demonstrated, evidence for the safety of higher administered doses is insufficient. Thus, preterm infants who lack clinical response to standard-dose caffeine therapy are of interest for TDM when dose optimization is performed. Polymorphisms in pharmacodynamics-related genes, but not in pharmacokinetics-related genes, have a significant impact on the interindividual variability in clinical response to caffeine therapy. For preterm infants lacking clinical response, how to develop individualized medication regimens for caffeine remains to be explored.

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

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

Development of a physiologically based pharmacokinetic (PBPK) population model for Chinese elderly subjects

Aims

This study aims to develop and verify a physiologically based pharmacokinetic (PBPK) population model for the Chinese geriatric population in Simcyp.

Methods

Firstly, physiological information for the Chinese geriatric population was collected and later employed to develop the Chinese geriatric population model by recalibration of corresponding physiological parameters in the Chinese adult population model available in Simcyp (i.e., Chinese healthy volunteer model). Secondly, drug‐dependent parameters were collected for six drugs with different elimination pathways (i.e., metabolized by CYP1A2, CYP3A4 or renal excretion). The drug models were then developed and verified by clinical data from Chinese adults, Caucasian adults and Caucasian elderly subjects to ensure that drug‐dependent parameters are correctly inputted. Finally, the tested drug models in combination with the newly developed Chinese geriatric population model were applied to simulate drug concentration in Chinese elderly subjects. The predicted results were then compared with the observations to evaluate model prediction performance.

Results

Ninety‐eight per cent of predicted AUC, 95% of predicted C_max, and 100% of predicted CL values were within two‐fold of the observed values, indicating all drug models were properly developed. The drug models, combined with the newly developed population model, were then used to predict pharmacokinetics in Chinese elderly subjects aged 60–93. The predicted AUC, C_max, and CL values were all within two‐fold of the observed values.

Conclusion

The population model for the Chinese elderly subjects appears to adequately predict the concentration of the drug that was metabolized by CYP1A2, CYP3A4 or eliminated by renal clearance.

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.

Does Prenatal Exposure to CNS Stimulants Increase the Risk of Cardiovascular Disease in Adult Offspring?

Prenatal exposure to an adverse uterine environment can have long lasting effects on adult offspring through DNA methylation, histone acetylation, and other epigenetic effects that alter gene expression and physiology. It is well-known that consumption of CNS stimulants such as caffeine, nicotine, amphetamines, and cocaine during pregnancy can adversely impact the offspring. However, most work in this area has focused on neurological and behavioral outcomes and has been limited to assessments in young offspring. The impact of prenatal exposure to these agents on the adult cardiovascular system has received relatively little attention. Evidence from both animal and human studies indicate that exposure to CNS stimulants during the gestational period can negatively impact the adult heart and vasculature, potentially leading to cardiovascular diseases later in life. This review discusses our current understanding of the impact of prenatal exposure to cocaine, methamphetamine, nicotine, and caffeine on the adult cardiovascular system.

Pharmacokinetic profile of amodiaquine and its active metabolite desethylamodiaquine in Ghanaian patients with uncomplicated falciparum malaria

BACKGROUND

Accurate measurement of anti-malarial drug concentrations in therapeutic efficacy studies is essential to distinguish between inadequate drug exposure and anti-malarial drug resistance, and to inform optimal anti-malarial dosing in key target population groups.

METHODS

A sensitive and selective LC-MS/MS method was developed and validated for the simultaneous determination of amodiaquine and its active metabolite, desethylamodiaquine, and used to describe their pharmacokinetic parameters in Ghanaian patients with uncomplicated falciparum malaria treated with the fixed-dose combination, artesunate-amodiaquine.

RESULTS

The day-28 genotype-adjusted adequate clinical and parasitological response rate in 308 patients studied was > 97% by both intention-to-treat and per-protocol analysis. After excluding 64 patients with quantifiable amodiaquine concentrations pre-treatment and 17 with too few quantifiable concentrations, the pharmacokinetic analysis included 227 patients (9 infants, 127 aged 1-4 years, 91 aged ≥ 5 years). Increased median day-3 amodiaquine concentrations were associated with a lower risk of treatment failure [HR 0.87 (95% CI 0.78-0.98), p = 0.021]. Amodiaquine exposure (median AUC0-∞) was significantly higher in infants (4201 ng h/mL) and children aged 1-5 years (1994 ng h/mL) compared to older children and adults (875 ng h/mL, p = 0.001), even though infants received a lower mg/kg amodiaquine dose (median 25.3 versus 33.8 mg/kg in older patients). Desethylamodiaquine AUC0-∞ was not significantly associated with age. No significant safety concerns were identified.

CONCLUSIONS

Efficacy of artesunate-amodiaquine at currently recommended dosage regimens was high across all age groups. Reassuringly, amodiaquine and desethylamodiaquine exposure was not reduced in underweight-for-age young children or those with high parasitaemia, two of the most vulnerable target populations. A larger pharmacokinetic study with close monitoring of safety, including full blood counts and liver function tests, is needed to confirm the higher amodiaquine exposure in infants, understand any safety implications and assess whether dose optimization in this vulnerable, understudied population is needed.

Therapy for Apnoea of Prematurity: A Retrospective Study on Effects of Standard Dose and Genetic Variability on Clinical Response to Caffeine Citrate in Chinese Preterm Infants

INTRODUCTION

Apnoea of prematurity (AOP) is among the most common diagnoses in the neonatal intensive care unit. Caffeine treatment is a preferred treatment choice. However, neonatal caffeine therapy results in significant intersubject variability. This study aimed to determine the effects of plasma caffeine levels based on standard dose and genetic variability on clinical response to caffeine citrate in Chinese preterm infants.

METHODS

This single-center and retrospective study examined data from 112 preterm infants (< 35 weeks gestational age) between July 2017 and July 2018. Subjects were divided into apnoea-free (n = 48) and apnoeic (n = 64) groups, and their clinical outcomes were summarized. Liquid chromatography-tandem mass spectrometry was used to measure levels of caffeine and its primary metabolites. Eighty-eight single-nucleotide polymorphisms were chosen for genotyping by a MassARRAY system.

RESULTS

Preterm infants in the apnoea-free group were associated with a reduction in the incidence of bronchopulmonary dysplasia and a reduced requirement for patent ductus arteriosus ligation. No significant association was observed between plasma-trough-concentration-to-dose (C₀/D) ratio and birth weight, gestational age, or postnatal age in either group. Polymorphisms in CYP1A2 and aryl hydrocarbon receptor (AHR) genes did not affect plasma caffeine levels. Polymorphisms in adenosine receptor genes ADORA1 (rs10920568 and rs12744240), ADORA2A (rs34923252 and rs5996696), and ADORA3 (rs10776727 and rs2298191), especially in AHR (rs4410790) and adenosine deaminase (rs521704), play critical roles in the interindividual response to caffeine therapy.

CONCLUSIONS

Genetic polymorphisms in caffeine's target receptors, but not the exposure levels based on the standard dosing, were associated with variable responses to caffeine therapy in preterm neonates. Future studies are needed to uncover how these genetic variants affect responses to caffeine therapy in this patient population.

Pharmacotherapy of neonatal opioid withdrawal syndrome: a review of pharmacokinetics and pharmacodynamics

INTRODUCTION

Neonatal opioid withdrawal syndrome (NOWS) often arises in infants born to mothers who used opioids during pregnancy. Morphine, methadone, and buprenorphine are the most common first-line treatments, whereas clonidine and phenobarbital are generally reserved for adjunctive therapy. These drugs exhibit substantial pharmacokinetic (PK) and pharmacodynamic (PD) variability. Current pharmacological treatments for NOWS are based on institutional protocols and largely rely on empirical treatment of patient symptoms.

AREAS COVERED

This article reviews the PK/PD of NOWS pharmacotherapies with a focus on the implication of physiological development and maturation. Body size-standardized clearance is consistently low in neonates, except for methadone. This can be ascribed to underdeveloped metabolic and elimination pathways. The effects of pharmacogenetics have been clarified especially for morphine. The PK/PD relationship of medications used in the treatment of NOWS is generally understudied.

EXPERT OPINION

Providing an appropriate opioid dose in neonates is challenging. Advancements in quantitative pharmacology and PK/PD modeling approaches facilitate identification of key factors driving PK/PD variability and characterization of exposure-response relationships. PK/PD model-informed simulations have been widely employed to define age-appropriate pediatric dosing regimens. The model-informed approach holds promise to aid more rational use of medications in the treatment of NOWS.

Treprostinil Attains Clinically Therapeutic Concentrations in Neonates with Pulmonary Hypertension on Extracorporeal Membrane Oxygenation Support

BACKGROUND

TREPROSTINIL IS A PROSTACYCLIN ANALOG USED FOR TREATMENT OF PULMONARY HYPERTENSION (PH) IN ADULTS AND CHILDREN, CURRENTLY AWAITING CLINICAL ASSESSMENT FOR USE IN NEONATES.: OBJECTIVES: WE AIMED TO INVESTIGATE THE USE OF TREPROSTINIL IN NEONATES WITH PH ON EXTRACORPOREAL MEMBRANE OXYGENATION (ECMO) SUPPORT AND MEASURE PLASMA CONCENTRATIONS OF THE DRUG.: METHODS: THIS IS A RETROSPECTIVE CASE-SERIES WITH PROSPECTIVELY COLLECTED BLOOD SAMPLES, CONDUCTED IN A QUATERNARY CARE NEONATAL INTENSIVE CARE UNIT. BRAIN NATRIURETIC PEPTIDE, CARDIAC FUNCTION ON DOPPLER ECHOCARDIOGRAPHY, AND THE OCCURRENCE OF ADVERSE EFFECTS WAS MONITORED. PLASMA CONCENTRATIONS WERE MEASURED USING HIGH-PERFORMANCE LIQUID CHROMATOGRAPHY AND MASS SPECTROMETRY.: RESULTS: FOUR PATIENTS WITH PH REQUIRING ECMO THERAPY WERE STUDIED. TREPROSTINIL DOSES OF 20-58 NG/KG/MIN REACHED CONCENTRATIONS OF 0.99-4.39 NG/ML AND INDUCED CLINICAL IMPROVEMENT. INFUSION OF TREPROSTINIL WAS ASSOCIATED WITH IMPROVED RIGHT VENTRICULAR FUNCTION, REVERSED RIGHT-TO-LEFT SHUNTING THROUGH THE DUCTUS ARTERIOSUS, AND STABLE OR DECREASING NEED FOR VASOPRESSOR SUPPORT.

NO ADVERSE EFFECTS WERE OBSERVED

CONCLUSIONS: THIS IS THE FIRST STUDY TO REPORT CLINICALLY THERAPEUTIC TREPROSTINIL CONCENTRATIONS IN CIRCULATING PLASMA AFTER TREPROSTINIL ADMINISTRATION IN NEONATES ON ECMO, WITH ASSOCIATED CLINICAL IMPROVEMENT OF PH AND NO SIGNS OF HEMODYNAMIC INSTABILITY.

Use of a physiologically based pharmacokinetic-pharmacodynamic model for initial dose prediction and escalation during a paediatric clinical trial

AIMS

To build and verify a physiologically based pharmacokinetic (PBPK) model for radiprodil in adults and link this to a pharmacodynamic (PD) receptor occupancy (RO) model derived from in vitro data. Adapt this model to the paediatric population and predict starting and escalating doses in infants based on RO. Use the model to guide individualized dosing in a clinical trial in 2- to 14-month-old children with infantile spasms.

METHODS

A PBPK model for radiprodil was developed to investigate the systemic exposure of the drug after oral administration in fasted and fed adults; this was then linked to RO via a PD model. The model was then expanded to include developmental physiology and ontogeny to predict escalating doses in infants that would result in a specific RO of 20, 40 and 60% based on average unbound concentration following a twice daily (b.i.d.) dosing regimen. Dose progression in the clinical trial was based on observed concentration-time data against PBPK predictions.

RESULTS

For paediatric predictions, the elimination of radiprodil, based on experimental evidence, had no ontogeny. Predicted b.i.d. doses ranged from 0.04 mg/kg for 20% RO, 0.1 mg/kg for 40% RO to 0.21 mg/kg for 60% RO. For all infants recruited in the study, observed concentration-time data following the 0.04 mg/kg and subsequent doses were within the PBPK model predicted 5^th and 95^th percentiles.

CONCLUSION

To our knowledge, this is the first time a PBPK model linked to RO has been used to guide dose selection and escalation in the live phase of a paediatric clinical trial.

Toxicokinetics of cis- and trans-Permethrin: Influence of Isomer, Maturation, and Sex

Permethrin exposure of children and adults is widespread in many populations, but knowledge of its relative toxicokinetics (TK) and health risks in immature age groups is lacking. Studies were conducted in rats to determine the influence of immaturity and sex (on plasma and target organ dosimetry of each of the insecticide’s 2 isomers, cis- and trans-permethrin [CIS and TRANS]). Postnatal day 15, 21, and 90 (adult), Sprague Dawley rats were orally administered a graduated series of doses of CIS and TRANS in corn oil. Serial sacrifices were conducted over 24 h to obtain plasma, brain, liver, skeletal muscle, and fat profiles of CIS and TRANS. Levels of TRANS decreased relatively rapidly, despite administration of relatively high doses. Concentrations of each isomer in plasma, brain, and other tissues monitored were inversely proportional to the animals’ age. The youngest pups exhibited 4-fold higher plasma and brain area under the curves than did adults. Little difference was observed in the TK of CIS or TRANS between adult male and female rats, other than higher initial plasma and liver CIS levels in females. Elevated exposure of the immature brain appears to be instrumental in increased susceptibility to the acute neurotoxicity of high-dose permethrin (Cantalamessa [1993]), but it remains to be established whether age-dependent TK is relevant to long-term, low-level risks.

A Time-Embedding Network Models the Ontogeny of 23 Hepatic Drug Metabolizing Enzymes

Pediatric patients are at elevated risk of adverse drug reactions, and there is insufficient information on drug safety in children. Complicating risk assessment in children, there are numerous age-dependent changes in the absorption, distribution, metabolism, and elimination of drugs. A key contributor to age-dependent drug toxicity risk is the ontogeny of drug metabolism enzymes, the changes in both abundance and type throughout development from the fetal period through adulthood. Critically, these changes affect not only the overall clearance of drugs but also exposure to individual metabolites. In this study, we introduce time-embedding neural networks in order to model population-level variation in metabolism enzyme expression as a function of age. We use a time-embedding network to model the ontogeny of 23 drug metabolism enzymes. The time-embedding network recapitulates known demographic factors impacting 3A5 expression. The time-embedding network also effectively models the nonlinear dynamics of 2D6 expression, enabling a better fit to clinical data than prior work. In contrast, a standard neural network fails to model these features of 3A5 and 2D6 expression. Finally, we combine the time-embedding model of ontogeny with additional information to estimate age-dependent changes in reactive metabolite exposure. This simple approach identifies age-dependent changes in exposure to valproic acid and dextromethorphan metabolites and suggests potential mechanisms of valproic acid toxicity. This approach may help researchers evaluate the risk of drug toxicity in pediatric populations.

Toxicokinetics of Deltamethrin: Dosage Dependency, Vehicle Effects, and Low-Dose Age-Equivalent Dosimetry in Rats

There is increasing concern that infants and children may be at increased risk of neurological effects of pyrethroids, the most widely used class of insecticide. The objectives of this investigation were to (1) characterize the dose-dependent toxicokinetics (TK) of deltamethrin (DLM) for exposures ranging from environmentally relevant to acutely toxic; (2) determine the influence of an aqueous versus oil vehicle on oral absorption and bioavailability; and (3) determine whether DLM exhibits low-dose, age-equivalent internal dosimetry. Serial arterial plasma samples were obtained for 72 h from adult, male Sprague Dawley rats given 0.05-5.0 mg DLM/kg as an oral bolus in corn oil (CO). DLM exhibited linear, absorption rate-limited TK. Increases in maximum plasma concentration (Cmax) and AUC∘∞ were directly proportional to the dose. Oral bioavailability was quite limited. The vehicle and its volume had modest effect on the rate and extent of systemic absorption in adult rats. Postnatal day (PND) 15, 21, and 90 (adult) rats received 0.10, 0.25, or 0.50 mg DLM/kg orally in CO and were sacrificed periodically for plasma, brain, and liver collection. Age-dependent differences between PND 15 and 90 plasma Cmax and AUC∘24 values progressively diminished as the dose decreased, but there was a lack of low dose age equivalence in these brain and liver dosimeters. Other maturational factors may account for the lack of the low-dose age equivalence in brain and liver. This investigation provides support for the premise that the relatively low metabolic capacity of immature subjects may be adequate to effectively eliminate trace amounts of DLM and other pyrethroids from the plasma.

The Impact of Scaling Factor Variability on Risk-Relevant Pharmacokinetic Outcomes in Children: A Case Study Using Bromodichloromethane (BDCM)

Biotransformation rates extrapolated from in vitro data are used increasingly in human physiologically based pharmacokinetic (PBPK) models. This practice requires use of scaling factors, including microsomal content (mg of microsomal protein/g liver, MPPGL), enzyme specific content, and liver mass as a fraction of body weight (FVL). Previous analyses indicated that scaling factor variability impacts pharmacokinetic (PK) outcomes used in adult population dose-response studies. This analysis was extended to pediatric populations because large inter-individual differences in enzyme ontogeny likely would further contribute to scaling factor variability. An adult bromodichloromethane (BDCM) model (Kenyon, E. M., Eklund, C., Leavens, T. L., and Pegram, R. A. (2016a). Development and application of a human PBPK model for bromodichloromethane (BDCM) to investigate impacts of multi-route exposure. J. Appl. Toxicol. 36, 1095-1111) was re-parameterized for neonates, infants, and toddlers. Monte Carlo analysis was used to assess the impact of pediatric scaling factor variation on model-derived PK outcomes compared with adult findings. BDCM dose metrics were estimated following a single 0.05-liter drink of water or a 20-min bath, under typical (5 µg/l) and plausible higher (20 µg/l) BDCM concentrations. MPPGL, CYP2E1, and FVL values reflected the distribution of reported pediatric population values. The impact of scaling factor variability on PK outcome variation was different for each exposure scenario, but similar for each BDCM water concentration. The higher CYP2E1 expression variability during early childhood was reflected in greater variability in predicted PK outcomes in younger age groups, particularly for the oral exposure route. Sensitivity analysis confirmed the most influential parameter for this variability was CYP2E1, particularly in neonates. These findings demonstrate the importance of age-dependent scaling factor variation used for in vitro to in vivo extrapolation of biotransformation rates.

Developmental regulation of kidney and liver solute carrier and ATP-binding cassette drug transporters and drug metabolizing enzymes: the role of remote organ communication

INTRODUCTION

The ontogeny of drug transport and metabolism is generally studied independently in tissues, yet in the immediate postnatal period the developmental regulation of SLC and ABC transporters and metabolizing enzymes must be coordinated. Using the Remote Sensing and Signaling Hypothesis as a framework, we describe how a systems physiology view helps to make sense of how inter-organ communication via hepatic, renal, and intestinal transporters and drug metabolizing enzymes (DMEs) is regulated from the immediate postnatal period through adulthood. Areas covered: This review examines patterns of developmental expression and function of transporters and DMEs with a focus on how cross-talk between these proteins in the kidney, liver and other organs (e.g., intestine) may be coordinated postnatally to optimize levels of metabolites and endogenous signaling molecules as well as gut-microbiome products. Expert opinion/commentary: Developmental expression is considered in terms of the Remote Sensing and Signaling Hypothesis, which addresses how transporters and DMEs participate in inter-organ and inter-organism small molecule communication in health, development, and disease. This hypothesis, for which there is growing support, is particularly relevant to the 'birth transition' and post-natal developmental physiology when organs must deal with critical physiological tasks distinct from the fetal period and where remote inter-organ and possibly inter-organismal (e.g. infant-gut microbiome) communication is likely to be critical to maintain homeostasis.