Physiologically based pharmacokinetic modeling to predict the impact of liver cirrhosis on glucuronidation via UGT1A4 & UGT2B7/2B4 - a case study with midazolam

Hepatic impairment, due to liver cirrhosis, decreases the activity of cytochrome P450 enzymes (CYPs). The use of physiologically-based pharmacokinetic (PBPK) modeling to predict this effect for CYP substrates has been well-established, but the effect of cirrhosis on uridine-glucuronosyltransferase (UGT) activities is less studied and few PBPK models have been reported. UGT enzymes are involved in primary N-glucuronidation of midazolam and glucuronidation of 1'-OH-midazolam following CYP3A hydroxylation. In this study Simcyp® was used to establish PBPK models for midazolam, its primary metabolites midazolam-N-glucuronide (UGT1A4) and 1'-OH midazolam (CYP3A4/3A5) and the secondary metabolite 1'-OH-midazolam-O-glucuronide (UGT2B7/2B4), allowing to simulate the impact of liver cirrhosis on the primary and secondary glucuronidation of midazolam. The model was verified in non-cirrhotic subjects before extrapolation to cirrhotic patients of Child-Pugh (CP) classes A, B, and C. Our model successfully predicted the exposures of midazolam and its metabolites in non-cirrhotic and cirrhotic patients, with 86% of observed plasma concentrations within 5th-95th percentiles of predictions and observed geometrical mean of AUCinf and Cmax within 0.7-1.43-fold of predictions. The simulated metabolic ratio (AUCglucuronide/AUCparent, MR), was calculated for midazolam-N-glucuronide to midazolam (indicative of UGT1A4 activity) and decreased by 40% (CP A), 48% (CP B), and 75% (CP C). For 1'-OH-midazolam-O-glucuronide to 1'-OH-midazolam the MR (indicative of UGT2B7/2B4 activity) dropped by 35% (CP A), 51% (CP B), and 64% (CP C). These predicted MRs were corroborated by the observed data. This work thus increases confidence in Simcyp® predictions of the effect of liver cirrhosis on the pharmacokinetics of UGT1A4 and UGT2B7/UGT2B4 substrates. Significance Statement This paper presents a PBPK model for midazolam and its metabolites and verifies the accurate simulation of pharmacokinetic profiles when using the Simcyp® hepatic impairment population models. Exposure changes of midazolam-N-glucuronide and 1'-OH-midazolam-O-glucuronide reflect the impact of decreases in UGT1A4 and UGT2B7/2B4 glucuronidation activity in cirrhotic patients. The approach used in this study may be extended to verify the modeling of other UGT enzymes affected by liver cirrhosis.

Determining the exposure of maternal medicines through breastfeeding: the UmbrelLACT study protocol-a contribution from the ConcePTION project

INTRODUCTION

Breastfeeding is beneficial for the health of the mother and child. However, at least 50% of postpartum women need pharmacotherapy, and this number is rising due to the increasing prevalence of chronic diseases and pregnancies at a later age. Making informed decisions on medicine use while breastfeeding is often challenging, considering the extensive information gap on medicine exposure and safety during lactation. This can result in the unnecessary cessation of breastfeeding, the avoidance of pharmacotherapy or the off-label use of medicines. The UmbrelLACT study aims to collect data on human milk transfer of maternal medicines, child exposure and general health outcomes. Additionally, the predictive performance of lactation and paediatric physiologically based pharmacokinetic (PBPK) models, a promising tool to predict medicine exposure in special populations, will be evaluated.

METHODS AND ANALYSIS

Each year, we expect to recruit 5-15 breastfeeding mothers using pharmacotherapy via the University Hospitals Leuven, the BELpREG project (pregnancy registry in Belgium) or external health facilities. Each request and compound will be evaluated on relevance (ie, added value to available scientific evidence) and feasibility (including access to analytical assays). Participants will be requested to complete at least one questionnaire on maternal and child's general health and collect human milk samples over 24 hours. Optionally, two maternal and one child's blood samples can be collected. The maternal medicine concentration in human milk will be determined along with the estimation of the medicine intake (eg, daily infant dose and relative infant dose) and systemic exposure of the breastfed child. The predictive performance of PBPK models will be assessed by comparing the observed concentrations in human milk and plasma to the PBPK predictions.

ETHICS AND DISSEMINATION

This study has been approved by the Ethics Committee Research UZ/KU Leuven (internal study number S67204). Results will be published in peer-reviewed journals and presented at (inter)national scientific meetings.

TRIAL REGISTRATION NUMBER

NCT06042803.

Second generation Al18F-labeled D-amino acid peptide for CXCR4 targeted molecular imaging

BACKGROUND

The C-X-C chemokine receptor type 4 (CXCR4) is overexpressed in many cancers, e.g. multiple myeloma and acute leukemia, yet solely [68Ga]PentixaFor is used for clinical PET imaging. The aim of this study was to develop and assess a second generation Al18F-labeled D-amino acid peptide based on the viral macrophage inflammatory protein II for CXCR4 targeted molecular imaging.

METHODS

We designed a library of monomer and multimer constructs and evaluated their binding affinity for human and mouse CXCR4. Based on these results, we selected the best vector molecule for development of an Al18F-labeled ligand, [18F]AlF-NOTA-2xDV1(c11sc12s), which was further evaluated in a cell-based binding assay to assess its binding properties and specificity for CXCR4. Next, pharmacokinetics and tumor uptake of [18F]AlF-NOTA-2xDV1(c11sc12s) were evaluated in naïve mice and mice with xenografts derived from U87.CXCR4 cells. Finally, we performed an imaging study in a non-human primate to assess the in vivo distribution of this novel radioligand in a species closely related to humans.

RESULTS

The lead ligand AlF-NOTA-2xDV1(c11sc12s) showed six-fold higher affinity for human CXCR4 compared to Ga-Pentixafor. The corresponding radiotracer was obtained in a good radiochemical yield of 40.1 ± 13.5 % (n = 4) and apparent molar activity of 20.4 ± 3.3 MBq/nmol (n = 4) after optimization. In U87.CD4.CXCR4 cell binding assays, the total bound fraction of [18F]AlF-NOTA-(2×)DV1(c11sc12s) was 32.4 ± 1.8 %. This fraction could be reduced by 82.5 % in the presence of 75 μM AMD3100. In naïve mice, [18F]AlF-NOTA-2xDV1(c11sc12s) accumulated in organs expressing mouse CXCR4, e.g. the liver (SUVmean (mean standardized uptake value) 75 min p.i. 11.7 ± 0.6), which was blockable by co-injecting AMD3100 (5 mg/kg). In U87.CXCR4 xenografted tumor mice, the tumor uptake of [18F]AlF-NOTA-2xDV1(c11sc12s) remained low (SUVmean 0.5 ± 0.1), but was reduced by co-administration of AMD3100. Surprisingly, [18F]AlF-NOTA-2xDV1(c11sc12s) exhibited a similar biodistribution in a non-human primate as in mice indicating off-target binding of [18F]AlF-NOTA-2xDV1(c11sc12s) in liver tissue. We confirmed that [18F]AlF-NOTA-2xDV1(c11sc12s) is taken up by hepatocytes using in vitro studies and that the uptake can be blocked with AMD3100 and rifampicin, a potent organic anion-transporting-polypeptide (OATP)1B1 and OATP1B3 inhibitor.

CONCLUSION

The second generation D-peptide AlF-NOTA-2xDV1(c11sc12s) showed high affinity for human CXCR4 and the corresponding radiotracer was produced in good radiochemical yields. However, [18F]AlF-NOTA-2xDV1(c11sc12s) is not specific for CXCR4 and is also a substrate for OATP1B1 and/or OATP1B3, known to mediate hepatic uptake. Therefore, D-amino acid peptides, based on the viral macrophage inflammatory protein II, are not the prefered vector molecule for the development of CXCR4 targeting molecular imaging tools.

Interdisciplinary Collaboration on Real World Data to Close the Knowledge Gap: A Reflection on "De Sutter et al. Predicting Volume of Distribution in Neonates: Performance of Physiologically Based Pharmacokinetic Modelling"

This commentary further reflects on the paper of De Sutter et al. on predicting volume of distribution in neonates, and the performance of physiologically based pharmacokinetic models We hereby stressed the add on value to collaborate on real world data to further close this knowledge gap. We illustrated this by weight distribution characteristics in breastfed (physiology) and in asphyxiated (pathophysiology), with additional reflection on their kidney and liver function.

A Population Model of Time-Dependent Changes in Serum Creatinine in (Near)term Neonates with Hypoxic-Ischemic Encephalopathy During and After Therapeutic Hypothermia

The objective was to apply a population model to describe the time course and variability of serum creatinine (sCr) in (near)term neonates with moderate to severe encephalopathy during and after therapeutic hypothermia (TH). The data consisted of sCr observations up to 10 days of postnatal age in neonates who underwent TH during the first 3 days after birth. Available covariates were birth weight (BWT), gestational age (GA), survival, and acute kidney injury (AKI). A previously published population model of sCr kinetics in neonates served as the base model. This model predicted not only sCr but also the glomerular filtration rate normalized by its value at birth (GFR/GFR0). The model was used to compare the TH neonates with a reference full term non-asphyxiated population of neonates. The estimates of the model parameters had good precision and showed high between subject variability. AKI influenced most of the estimated parameters denoting a strong impact on sCr kinetics and GFR. BWT and GA were not significant covariates. TH transiently increased [Formula: see text] in TH neonates over the first days compared to the reference group. Asphyxia impacted not only GFR, but also the [Formula: see text] synthesis rate. We also observed that AKI neonates exhibit a delayed onset of postnatal GFR increase and have a higher [Formula: see text] synthesis rate compared to no-AKI patients. Our findings show that the use of [Formula: see text] as marker of renal function in asphyxiated neonates treated with TH to guide dose selection for renally cleared drugs is challenging, while we captured the postnatal sCr patterns in this specific population.

Challenges Related to Acquisition of Physiological Data for Physiologically Based Pharmacokinetic (PBPK) Models in Postpartum, Lactating Women and Breastfed Infants-A Contribution from the ConcePTION Project

Physiologically based pharmacokinetic (PBPK) modelling is a bottom-up approach to predict pharmacokinetics in specific populations based on population-specific and medicine-specific data. Using an illustrative approach, this review aims to highlight the challenges of incorporating physiological data to develop postpartum, lactating women and breastfed infant PBPK models. For instance, most women retain pregnancy weight during the postpartum period, especially after excessive gestational weight gain, while breastfeeding might be associated with lower postpartum weight retention and long-term weight control. Based on a structured search, an equation for human milk intake reported the maximum intake of 153 mL/kg/day in exclusively breastfed infants at 20 days, which correlates with a high risk for medicine reactions at 2-4 weeks in breastfed infants. Furthermore, the changing composition of human milk and its enzymatic activities could affect pharmacokinetics in breastfed infants. Growth in breastfed infants is slower and gastric emptying faster than in formula-fed infants, while a slower maturation of specific metabolizing enzymes in breastfed infants has been described. The currently available PBPK models for these populations lack structured systematic acquisition of population-specific data. Future directions include systematic searches to fully identify physiological data. Following data integration as mathematical equations, this holds the promise to improve postpartum, lactation and infant PBPK models.

Automated Generation of hiPSC-Derived Hepatic Progeny by Cost-Efficient Compounds

Human pluripotent stem cell (hPSC)-derived hepatocyte-like cells (HLCs) hold great promise for liver disease modeling, drug discovery, and drug toxicity screens. Yet, several hurdles still need to be overcome, including among others decrease in the cost of goods to generate HLCs and automation of the differentiation process. We here describe that the use of an automated liquid handling system results in highly reproducible HLC differentiation from hPSCs. This enabled us to screen 92 chemicals to replace expensive growth factors at each step of the differentiation protocol to reduce the cost of goods of the differentiation protocol by approximately 79%. In addition, we also evaluated several recombinant extracellular matrices to replace Matrigel. We demonstrated that differentiation of hPSCs on Laminin-521 using an optimized small molecule combination resulted in HLCs that were transcriptionally identical to HLCs generated using the growth factor combinations. In addition, the HLCs created using the optimized small molecule combination secreted similar amounts of albumin and urea, and relatively low concentrations of alfa-fetoprotein, displayed similar CYP3A4 functionality, and a similar drug toxicity susceptibility as HLCs generated with growth factor cocktails. The broad applicability of the new differentiation protocol was demonstrated for 4 different hPSC lines. This allowed the creation of a scalable, xeno-free, and cost-efficient hPSC-derived HLC culture, suitable for high throughput disease modeling and drug screenings, or even for the creation of HLCs for regenerative therapies.

PBPK Modelling for Drugs Cleared by Non-CYP Enzymes: State-of-the-Art and Future Perspectives

Physiologically-based pharmacokinetic (PBPK) modeling has become the established method for predicting human pharmacokinetics (PK) and drug-drug interactions (DDI). The number of drugs cleared by non-CYP enzyme metabolism has increased steadily and to date, there is no consolidated overview of PBPK modeling for drugs cleared by non-CYP enzymes. This review aims to describe the state-of-the-art for PBPK modeling for drugs cleared via non-CYP enzymes, to identify successful strategies, to describe gaps and to provide suggestion to overcome them. To this end, we conducted a detailed literature search and found 58 articles published before the 1st of January 2023 containing 95 examples of clinical PBPK models for 62 non-CYP enzyme substrates. Reviewed articles covered the drug clearance by uridine 5'-diphospho-glucuronosyltransferases (UGTs), aldehyde oxidase (AO), flavin-containing monooxygenases (FMOs), sulfotransferases (SULTs) and carboxylesterases (CES), with UGT2B7, UGT1A9, CES1, FMO3 and AO being the enzymes most frequently involved. In vitro-in vivo extrapolation (IVIVE) of intrinsic clearance and the bottom-up PBPK modeling involving non-CYP enzymes remains challenging. We observed that the middle-out modeling approach was applied in 80% of the cases, with metabolism parameters optimized in 73% of the models. Our review could not identify a standardized approach used for model optimization based on clinical data, with manual optimization employed most frequently. Successful development of models for UGT2B7, UGT1A9, CES1, and FMO3 substrates provides a foundation for other drugs metabolized by these enzymes and guides the way forward in creating PBPK models for other enzymes in these families. Significance Statement Our review charts the rise of PBPK modeling for drugs cleared by non-CYP enzymes. Analyzing 58 articles and 62 non-CYP enzyme substrates, we found that UGTs, AO, FMOs, SULTs, and CES were the main enzyme families involved and that UGT2B7, UGT1A9, CES1, FMO3 and AO are the individual enzymes with the strongest PBPK modeling precedents. Approaches established for these enzymes can now be extended to additional substrates and to drugs metabolized by enzymes that are similarly well characterized.

Physiologically based pharmacokinetic modelling to inform combination dosing regimens of ceftaroline and daptomycin in special populations

AIMS

The combination of daptomycin and ceftaroline used as salvage therapy is associated with higher survival and decreased clinical failure in complicated methicillin-resistant Staphylococcus aureus (MRSA) infections that are resistant to standard MRSA treatment. This study aimed to evaluate dosing regimens for coadministration of daptomycin and ceftaroline in special populations including paediatrics, renally impaired (RI), obese and geriatrics that generate sufficient coverage against daptomycin-resistant MRSA.

METHODS

Physiologically based pharmacokinetic models were developed from pharmacokinetic studies of healthy adults, geriatric, paediatric, obese and RI patients. The predicted profiles were used to evaluate joint probability of target attainment (PTA), as well as tissue-to-plasma ratios.

RESULTS

The adult dosing regimens of 6 mg/kg every (q)24h or q48h daptomycin and 300-600 mg q12h ceftaroline fosamil by RI categories achieved ≥90% joint PTA when the minimum inhibitory concentrations in the combination are at or below 1 and 4 μg/mL against MRSA. In paediatrics, wherein there is no recommended daptomycin dosing regimen for S. aureus bacteraemia, ≥90% joint PTA is achieved when the minimum inhibitory concentrations in the combination are up to 0.5 and 2 μg/mL for standard paediatric dosing regimens of 7 mg/kg q24h daptomycin and 12 mg/kg q8h ceftaroline fosamil. Model predicted tissue-to-plasma ratios of 0.3 and 0.7 in the skin and lung, respectively, for ceftaroline and 0.8 in the skin for daptomycin.

CONCLUSION

Our work illustrates how physiologically based pharmacokinetic modelling can inform appropriate dosing of adult and paediatric patients and thereby enable prediction of target attainment in the patients during multitherapies.

Sensitive and rapid method for the quantitation of amoxicillin in minipig plasma and milk by LC-MS/MS: A contribution from the IMI ConcePTION project

The IMI project ConcePTION was launched to fill the knowledge gap of using medicines during pregnancy and lactation. To achieve this goal, several studies are being conducted, including the bioanalysis of amoxicillin in minipig plasma and milk. A high-throughput, robust and reliable liquid chromatography tandem mass spectrometry method was developed and validated according to FDA and EMA guidelines to determine the concentrations of amoxicillin in a large number of minipig plasma and milk samples. Chromatographic separation was achieved on a Luna® Omega Polar C18, 1.6 μm, 100 × 2.1 mm column, with a mobile phase consisting of 0.1% formic acid in water and acetonitrile. Mass spectrometry used in a positive ionization mode and the transitions m/z 366.1 → 349.2 was selected to monitor amoxicillin, while m/z 370.1 → 114.15 was selected for the stable isotope labelled internal standard. This method features a linear quantification range of 10 ng/mL - 10 μg/mL, recovery of not less than 94.1%, a single sample extraction method for both plasma and milk matrices, and an analysis runtime of 5 min.

Pharmacokinetics during therapeutic hypothermia in neonates: from pathophysiology to translational knowledge and physiologically-based pharmacokinetic (PBPK) modeling

INTRODUCTION

Perinatal asphyxia (PA) still causes significant morbidity and mortality. Therapeutic hypothermia (TH) is the only effective therapy for neonates with moderate to severe hypoxic-ischemic encephalopathy after PA. These neonates need additional pharmacotherapy, and both PA and TH may impact physiology and, consequently, pharmacokinetics (PK) and pharmacodynamics (PD).

AREAS COVERED

This review provides an overview of the available knowledge in PubMed (until November 2022) on the pathophysiology of neonates with PA/TH. In vivo pig models for this setting enable distinguishing the effect of PA versus TH on PK and translating this effect to human neonates. Available asphyxia pig models and methodological considerations are described. A summary of human neonatal PK of supportive pharmacotherapy to improve neurodevelopmental outcomes is provided.

EXPERT OPINION

To support drug development for this population, knowledge from clinical observations (PK data, real-world data on physiology), preclinical (in vitro and in vivo (minipig)) data, and molecular and cellular biology insights can be integrated into a predictive physiologically-based PK (PBPK) framework, as illustrated by the I-PREDICT project (Innovative physiology-based pharmacokinetic model to predict drug exposure in neonates undergoing cooling therapy). Current knowledge, challenges, and expert opinion on the future directions of this research topic are provided.

Unravelling the Hepatic Elimination Mechanisms of Colistin

PURPOSE

Colistin is an antibiotic which is increasingly used as a last-resort therapy in critically-ill patients with multidrug resistant Gram-negative infections. The purpose of this study was to evaluate the mechanisms underlying colistin's pharmacokinetic (PK) behavior and to characterize its hepatic metabolism.

METHODS

In vitro incubations were performed using colistin sulfate with rat liver microsomes (RLM) and with rat and human hepatocytes (RH and HH) in suspension. The uptake of colistin in RH/HH and thefraction of unbound colistin in HH (fu,hep) was determined. In vitro to in vivo extrapolation (IVIVE) was employed to predict the hepatic clearance (CLh) of colistin.

RESULTS

Slow metabolism was detected in RH/HH, with intrinsic clearance (CLint) values of 9.34± 0.50 and 3.25 ± 0.27 mL/min/kg, respectively. Assuming the well-stirred model for hepatic drug elimination, the predicted rat CLh was 3.64± 0.22 mL/min/kg which could explain almost 70% of the reported non-renal in vivo clearance. The predicted human CLh was 91.5 ± 8.83 mL/min, which was within two-fold of the reported plasma clearance in healthy volunteers. When colistin was incubated together with the multidrug resistance-associated protein (MRP/Mrp) inhibitor benzbromarone, the intracellular accumulation of colistin in RH/HH increased significantly.

CONCLUSION

These findings indicate the major role of hepatic metabolism in the non-renal clearance of colistin, while MRP/Mrp-mediated efflux is involved in the hepatic disposition of colistin. Our data provide detailed quantitative insights into the hereto unknown mechanisms responsible for non-renal elimination of colistin.

Safety of medicines during breastfeeding - from case report to modeling : A contribution from the ConcePTION project

INTRODUCTION

Despite many research efforts, current data on the safety of medicines during breastfeeding are either fragmented or lacking, resulting in restrictive labeling of most medicines. In the absence of pharmacoepidemiologic safety studies, risk estimation for breastfed infants is mainly derived from pharmacokinetic (PK) information on the medicine. This manuscript provides a description and a comparison of the different methodological approaches that can yield reliable information on medicine transfer into human milk and the resulting infant exposure.

AREA COVERED

Currently, most information on medicine transfer in human milk relies on case reports or traditional PK studies, which generate data that can hardly be generalized to the population. Some methodological approaches, such as population PK (popPK) and physiologically-based PK (PBPK) modeling, can be used to provide a more complete characterization of infant medicine exposure through human milk and simulate the most extreme situations, while decreasing the burden of sampling in breastfeeding women.

EXPERT OPINION

PBPK and popPK modeling are promising approaches to fill the gap of knowledge in medicine safety in breastfeeding, as illustrated with our escitalopram example.

Generic Workflow to Predict Medicine Concentrations in Human Milk Using Physiologically-Based Pharmacokinetic (PBPK) Modelling-A Contribution from the ConcePTION Project

Women commonly take medication during lactation. Currently, there is little information about the exposure-related safety of maternal medicines for breastfed infants. The aim was to explore the performance of a generic physiologically-based pharmacokinetic (PBPK) model to predict concentrations in human milk for ten physiochemically diverse medicines. First, PBPK models were developed for "non-lactating" adult individuals in PK-Sim/MoBi v9.1 (Open Systems Pharmacology). The PBPK models predicted the area-under-the-curve (AUC) and maximum concentrations (C_max) in plasma within a two-fold error. Next, the PBPK models were extended to include lactation physiology. Plasma and human milk concentrations were simulated for a three-months postpartum population, and the corresponding AUC-based milk-to-plasma (M/P) ratios and relative infant doses were calculated. The lactation PBPK models resulted in reasonable predictions for eight medicines, while an overprediction of human milk concentrations and M/P ratios (>2-fold) was observed for two medicines. From a safety perspective, none of the models resulted in underpredictions of observed human milk concentrations. The present effort resulted in a generic workflow to predict medicine concentrations in human milk. This generic PBPK model represents an important step towards an evidence-based safety assessment of maternal medication during lactation, applicable in an early drug development stage.

Development of a neonatal Göttingen Minipig model for dose precision in perinatal asphyxia: technical opportunities, challenges, and potential further steps

Animal models provide useful information on mechanisms in human disease conditions, but also on exploring (patho)physiological factors affecting pharmacokinetics, safety, and efficacy of drugs in development. Also, in pediatric patients, nonclinical data can be critical for better understanding the disease conditions and developing new drug therapies in this age category. For perinatal asphyxia (PA), a condition defined by oxygen deprivation in the perinatal period and possibly resulting in hypoxic ischemic encephalopathy (HIE) or even death, therapeutic hypothermia (TH) together with symptomatic drug therapy, is the standard approach to reduce death and permanent brain damage in these patients. The impact of the systemic hypoxia during PA and/or TH on drug disposition is largely unknown and an animal model can provide useful information on these covariates that cannot be assessed separately in patients. The conventional pig is proven to be a good translational model for PA, but pharmaceutical companies do not use it to develop new drug therapies. As the Göttingen Minipig is the commonly used pig strain in nonclinical drug development, the aim of this project was to develop this animal model for dose precision in PA. This experiment consisted of the instrumentation of 24 healthy male Göttingen Minipigs, within 24 h of partus, weighing approximately 600 g, to allow the mechanical ventilation and the multiple vascular catheters inserted for maintenance infusion, drug administration and blood sampling. After premedication and induction of anesthesia, an experimental protocol of hypoxia was performed, by decreasing the inspiratory oxygen fraction (FiO₂) at 15%, using nitrogen gas. Blood gas analysis was used as an essential tool to evaluate oxygenation and to determine the duration of the systemic hypoxic insult to approximately 1 h. The human clinical situation was mimicked for the first 24 h after birth in case of PA, by administering four compounds (midazolam, phenobarbital, topiramate and fentanyl), frequently used in a neonatal intensive care unit (NICU). This project aimed to develop the first neonatal Göttingen Minipig model for dose precision in PA, allowing to separately study the effect of systemic hypoxia versus TH on drug disposition. Furthermore, this study showed that several techniques that were thought to be challenging or even impossible in these very small animals, such as endotracheal intubation and catheterization of several veins, are feasible by trained personnel. This is relevant information for laboratories using the neonatal Göttingen Minipig for other disease conditions or drug safety testing.

Dose adjustment in patients with liver cirrhosis - comparison of two different modeling approaches

Failure to perform adequate dose adjustment in patients with liver cirrhosis may be associated with increased toxicity. We compared the prediction of AUC and clearance for the six compounds of the Basel phenotyping cocktail (caffeine, efavirenz, flurbiprofen, omeprazole, metoprolol, midazolam) using a well-known physiology-based pharmacokinetic approach ("PBPK approach", Simcyp®) and a novel top-down method based on the systemic clearance in healthy volunteers adjusted for markers of liver and renal dysfunction ("top-down approach"). With few exceptions, plasma concentration-time curves were accurately predicted by the PBPK approach. In comparison to the measured AUC and clearance of these drugs in patients with liver cirrhosis and healthy controls, except for efavirenz, the estimates of both approaches were within 2 standard deviations of the mean for total and free drug concentrations. For both approaches, a correction factor for dose adjustment in patients with liver cirrhosis could be calculated for the drugs administered. AUCs calculated using the adjusted doses were comparable to the AUCs measured in control subjects, with slightly more accurate predictions generated by the PBPK approach. For drugs with a free fraction <50%, predictions using free drug concentrations were more accurate than with total drug concentrations. In conclusion, both methods provided good qualitative predictions of the changes by liver cirrhosis in the pharmacokinetics of the six compounds investigated. The top-down approach is easier to implement but the PBPK approach predicted changes in drug exposure more accurately than the top-down approach and provided reliable estimates for plasma concentrations.

PBPK-based dose finding for sildenafil in pregnant women for antenatal treatment of congenital diaphragmatic hernia

Sildenafil is a potent vasodilator and phosphodiesterase type five inhibitor, commercially known as Revatio® and approved for the treatment of pulmonary arterial hypertension. Maternal administration of sildenafil during pregnancy is being evaluated for antenatal treatment of several conditions, including the prevention of pulmonary hypertension in fetuses with congenital diaphragmatic hernia. However, determination of a safe and effective maternal dose to achieve adequate fetal exposure to sildenafil remains challenging, as pregnancy almost always is an exclusion criterion in clinical studies. Physiologically-based pharmacokinetic (PBPK) modelling offers an attractive approach for dose finding in this specific population. The aim of this study is to exploit physiologically-based pharmacokinetic modelling to predict the required maternal dose to achieve therapeutic fetal exposure for the treatment congenital diaphragmatic hernia. A full-PBPK model was developed for sildenafil and N-desmethyl-sildenafil using the Simcyp simulator V21 platform, and verified in adult reference individuals, as well as in pregnant women, taking into account maternal and fetal physiology, along with factors known to determine hepatic disposition of sildenafil. Clinical pharmacokinetic data in mother and fetus were previously obtained in the RIDSTRESS study and were used for model verification purposes. Subsequent simulations were performed relying either on measured values for fetal fraction unbound (fu = 0.108) or on values predicted by the simulator (fu = 0.044). Adequate doses were predicted according to the efficacy target of 15 ng/mL (or 38 ng/mL) and safety target of 166 ng/mL (or 409 ng/mL), assuming measured (or predicted) fu values, respectively. Considering simulated median profiles for average steady state sildenafil concentrations, dosing regimens of 130 mg/day or 150 mg/day (administered as t.i.d.), were within the therapeutic window, assuming either measured or predicted fu values, respectively. For safety reasons, dosing should be initiated at 130 mg/day, under therapeutic drug monitoring. Additional experimental measurements should be performed to confirm accurate fetal (and maternal) values for fu. Additional characterization of pharmacodynamics in this specific population is required and may lead to further optimization of the dosing regimen.

Clinicopathological Features and Survival of Patients with Hepatocellular Carcinoma in Ethiopia: A Multicenter Study

(1) Background: Hepatocellular carcinoma (HCC) is one of the deadliest cancers globally, killing over 700,000 people each year. Despite the rising incidence and mortality rates of HCC in Ethiopia, only few single-centered studies have been conducted; therefore, we aimed to explore the clinicopathological characteristics and survival of patients with HCC in multicenter settings. (2) Methods: We conducted a retrospective analysis of 369 patients with confirmed HCC diagnosed between 2016 and 2021. The survival of patients weas determined using the Kaplan−Meier method, and hazard ratios of the prognostic factors were estimated in Cox proportional hazard models. (3) Results: Majority patients were male (67%) and had a mean age of 52.0 ± 15.6 years. The majority of patients (87%) had a large tumor size (>5 cm) at diagnosis and presented with an advanced-stage condition. Cirrhosis (58%) and viral hepatitis (46.5%) were the main risk factors associated with HCC. The median overall survival was 141 days (95% CI: 117−165). Patients who took antivirals for HBV had a higher survival benefit compared to the untreated group (469 vs. 104 days; p < 0.001). The risk of death was 12 times higher in patients with Barcelona Clinic Liver Cancer-D (BCLC-D) terminal stage HCC compared to patients with an early stage (BCLC-A) HCC. The stage of HCC and treatment against HBV are the most significant survival predictors. (4) Conclusions: The overall survival of HCC patients in Ethiopia is poor. Cirrhosis and viral hepatitis are the primary risk factors linked with HCC. Patients who received antiviral therapy for HBV had a better survival outcome.

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.

Ceftriaxone dosing based on the predicted probability of augmented renal clearance in critically ill patients with pneumonia

OBJECTIVES

PTA of protein-unbound ceftriaxone may be compromised in critically ill patients with community-acquired pneumonia (CAP) with augmented renal clearance (ARC). We aimed to determine an optimized ceftriaxone dosage regimen based on the probability of developing ARC on the next day (PARC,d+1; www.arcpredictor.com).

PATIENTS AND METHODS

Thirty-three patients enrolled in a prospective cohort study were admitted to the ICU with severe CAP and treated with ceftriaxone 2 g once daily. Patients contributed 259 total ceftriaxone concentrations, collected during 1 or 2 days (±7 samples/day). Unbound fractions of ceftriaxone were determined in all peak and trough samples (n = 76). Population pharmacokinetic modelling and simulation were performed using NONMEM7.4. Target attainment was defined as an unbound ceftriaxone concentration >4 mg/L throughout the dosing interval.

RESULTS

A two-compartment population pharmacokinetic model described the data well. The maximal protein-bound ceftriaxone concentration decreased with lower serum albumin. Ceftriaxone clearance increased with body weight and PARC,d+1 determined on the previous day. A high PARC,d+1 was identified as a clinically relevant predictor for underexposure on the next day (area under the receiver operating characteristics curve 0.77). Body weight had a weak predictive value and was therefore considered clinically irrelevant. Serum albumin had no predictive value. An optimal PARC,d+1 threshold of 5.7% was identified (sensitivity 73%, specificity 69%). Stratified once- or twice-daily 2 g dosing when below or above the 5.7% PARC,d+1 cut-off, respectively, was predicted to result in 81% PTA compared with 47% PTA under population-level once-daily 2 g dosing.

CONCLUSIONS

Critically ill patients with CAP with a high PARC,d+1 may benefit from twice-daily 2 g ceftriaxone dosing for achieving adequate exposure on the next day.

Current and future physiologically based pharmacokinetic (PBPK) modeling approaches to optimize pharmacotherapy in preterm neonates

INTRODUCTION

There is a need for structured approaches to inform on pharmacotherapy in preterm neonates. With their proven track record up to regulatory acceptance, physiologically based pharmacokinetic (PBPK) modeling and simulation provide such a structured approach, and hold the promise to support drug development in preterm neonates.

AREAS COVERED

Compared to the general and pediatric use of PBPK modeling, its use to inform pharmacotherapy in preterms is limited. Using a systematic search (PBPK + preterm), we retained 25 records (20 research papers, 2 letters, 3 abstracts). We subsequently collated the published information on PBPK software packages (PK-Sim®, Simcyp®), and their applications and optimization efforts in preterm neonates. It is encouraging that these applications cover a broad range of scenarios (pharmacokinetic-dynamic analyses, drug-drug interactions, developmental pharmacogenetics, lactation related exposure) and compounds (small molecules, proteins). Furthermore, specific compartments (cerebrospinal fluid, tissue) or (patho)physiologic processes (cardiac output, biliary excretion, first pass metabolism) are considered.

EXPERT OPINION

Knowledge gaps exist, giving rise to various levels of model uncertainty in PBPK applications in preterm neonates. To improve this setting, we need cross talk between clinicians and modelers to generate and integrate knowledge (PK datasets, system knowledge, maturational physiology and pathophysiology) to further refine PBPK models.

Case Report: Bosentan and Sildenafil Exposure in Human Milk - A Contribution From the ConcePTION Project

Introduction: Quantitative information on disposition of maternal medicines in human milk remains a major knowledge gap. This case report presents the clinical and pharmacokinetic data of a single mother-infant pair exposed to bosentan and sildenafil for the treatment of pulmonary arterial hypertension (PAH) during lactation.

Case presentation: A 43-year old mother was treated with sildenafil (20 mg, 3x/day) and bosentan (125 mg, 2x/day) for PAH. Her 21-months old infant received breastfeeding in combination with adequate complementary foods. Milk samples were collected over 24 h, at day 637 and 651 after delivery. The observed average steady-state concentrations of sildenafil (2.84 μg/L) and bosentan (49.0 μg/L) in human milk were low. The Daily Infant Dosage ingested by the nursing infant through human milk was 0.02 μg/kg/day for sildenafil and 0.29 μg/kg/day for bosentan at day 637, and 0.03 μg/kg/day and 0.60 μg/kg/day at day 651. The Relative Infant Dose calculated for an exclusively breastfed infant with an estimated milk intake of 150 ml/kg/day, was 0.06% for sildenafil and 0.24% for bosentan. General health outcome of the infant, reported by the mother, was uneventful until the sampling days.

Conclusion: Low medicine concentrations were found in human milk expressed 21 months after delivery after maternal intake of 20 mg sildenafil three times daily and 125 mg bosentan twice daily. General health of the nursing infant until sampling was reported as optimal by the mother.

Effects of Drugs Formerly Proposed for COVID-19 Treatment on Connexin43 Hemichannels

Connexin43 (Cx43) hemichannels form a pathway for cellular communication between the cell and its extracellular environment. Under pathological conditions, Cx43 hemichannels release adenosine triphosphate (ATP), which triggers inflammation. Over the past two years, azithromycin, chloroquine, dexamethasone, favipiravir, hydroxychloroquine, lopinavir, remdesivir, ribavirin, and ritonavir have been proposed as drugs for the treatment of the coronavirus disease 2019 (COVID-19), which is associated with prominent systemic inflammation. The current study aimed to investigate if Cx43 hemichannels, being key players in inflammation, could be affected by these drugs which were formerly designated as COVID-19 drugs. For this purpose, Cx43-transduced cells were exposed to these drugs. The effects on Cx43 hemichannel activity were assessed by measuring extracellular ATP release, while the effects at the transcriptional and translational levels were monitored by means of real-time quantitative reverse transcriptase polymerase chain reaction analysis and immunoblot analysis, respectively. Exposure to lopinavir and ritonavir combined (4:1 ratio), as well as to remdesivir, reduced Cx43 mRNA levels. None of the tested drugs affected Cx43 protein expression.

Predicting model‐informed precision dosing: A test‐case in tacrolimus dose adaptation for kidney transplant recipients

Before investing resources into the development of a precision dosing (model‐informed precision dosing [MIPD]) tool for tacrolimus, the performance of the tool was evaluated in silico. A retrospective dataset of 315 de novo kidney transplant recipients was first used to identify a one‐compartment pharmacokinetic (PK) model with time‐dependent clearance. MIPD performance was subsequently evaluated by calculating errors to predict future concentrations, which is directly related to dosing precision and probability of target attainment (PTA). Based on the identified model residual error, the theoretical upper limit was 45% PTA for a target of 13.5 ng/ml and an acceptable range of 12–15 ng/ml. Using empirical Bayesian estimation, this limit was reached on day 5 post‐transplant and beyond. By incorporating correlated within‐patient variability when predicting future individual concentrations, PTA improved beyond the theoretical upper limit. This yielded a Bayesian feedback dosing algorithm accurately predicting future trough concentrations and adapting each dose to reach a target concentration. Simulated concentration‐time profiles were then used to quantify MIPD‐based improvement on three end points: average PTA increased from 28% to 39%, median time to three concentrations in target decreased from 10 to 7 days, and mean log‐squared distance to target decreased from 0.080 to 0.055. A study of 200 patients was predicted to have sufficient power to demonstrate these nuanced PK end points reliably. These simulations supported our decision to develop a precision dosing tool for tacrolimus and test it in a prospective trial.

Meropenem Target Attainment and Population Pharmacokinetics in Critically Ill Septic Patients with Preserved or Increased Renal Function

Purpose

Critically ill patients with preserved or increased renal function have been shown to be at risk of underexposure to meropenem. Although many meropenem population pharmacokinetic (PK) models have been published, there is no large prospective population PK study with rich sampling focusing on patients most at risk of suboptimal pharmacokinetic/pharmacodynamic (PK/PD) target attainment. Therefore, the aim of the present study was to evaluate PK/PD target attainment and to perform a thorough covariate screening using population PK modelling of meropenem in septic patients with preserved or increased renal function.

Patients and Methods

A single-centre prospective observational PK study was performed in the intensive care unit (ICU) of the University Hospitals Leuven. Patients with severe sepsis or septic shock and treated with meropenem in the ICU were screened for inclusion. Patients were excluded if they received renal replacement therapy or had an estimated glomerular filtration rate according to the Chronic Kidney Disease Epidemiology collaboration equation <70 mL/min/1.73m² on the day of PK sampling. Successful PK/PD target attainment was defined as an unbound meropenem trough concentration above 2 mg/L or 8 mg/L. Population PK modelling was performed with NONMEM7.4.

Results

In total, 58 patients were included, contributing 345 plasma samples over 70 dosing intervals. The 2 mg/L and 8 mg/L targets were successfully attained in 46% and 11% of all dosing intervals, respectively. A two-compartment population PK model with linear elimination and interindividual variability on clearance best described meropenem PK. The estimated creatinine clearance according to the Cockcroft-Gault equation was the only covariate retained during population PK analysis.

Conclusion

This study provided detailed insight into meropenem PK in critically ill patients with preserved or increased renal function. We observed poor PK/PD target attainment, for which renal function was the only significant covariate.

Trial Registration

This study is registered at ClinicalTrials.gov (NCT03560557).

Serum Creatinine Patterns in Neonates Treated with Therapeutic Hypothermia for Neonatal Encephalopathy

INTRODUCTION

There is large variability in kidney function and injury in neonates with neonatal encephalopathy (NE) treated with therapeutic hypothermia (TH). Acute kidney injury (AKI) definitions that apply categorical approaches may lose valuable information about kidney function in individual patients. Centile serum creatinine (SCr) over postnatal age (PNA) may provide more valuable information in TH neonates.

METHODS

Data from seven TH neonates and one non-TH-treated, non-NE control cohorts were pooled in a retrospective study. SCr centiles over PNA, and AKI incidence (definition: SCr ↑≥0.3 mg/dL within 48 h, or ↑ ≥1.5 fold vs. the lowest prior SCr within 7 days) and mortality were calculated. Repeated measurement linear models were applied to SCr trends, modeling SCr on PNA, birth weight or gestational age (GA), using heterogeneous autoregressive residual covariance structure and maximum likelihood methods. Findings were compared to patterns in the control cohort.

RESULTS

Among 1,136 TH neonates, representing 4,724 SCr observations, SCr (10th-25th-50th-75th-90th-95th) PNA centiles (day 1-10) were generated. In TH neonates, the AKI incidence was 132/1,136 (11.6%), mortality 193/1,136 (17%). AKI neonates had a higher mortality (37.2-14.3%, p < 0.001). Median SCr patterns over PNA were significantly higher in nonsurvivors (p < 0.01) or AKI neonates (p < 0.001). In TH-treated neonates, PNA and GA or birth weight explained SCr variability. Patterns over PNA were significantly higher in TH neonates to controls (801 neonates, 2,779 SCr).

CONCLUSIONS

SCr patterns in TH-treated NE neonates are specific. Knowing PNA-related patterns enable clinicians to better assess kidney function and tailor pharmacotherapy, fluids, or kidney supportive therapies.

Antihistamine use during breastfeeding with focus on breast milk transfer and safety in humans: A systematic literature review

Current data on use of antihistamines during breastfeeding and risks to the breastfed infant are insufficient. The aim of this systematic review was to provide an overview of studies measuring the levels of antihistamines in human breast milk, estimating the exposure for breastfed infants and/or reporting possible adverse effects on the breastfed infant. An additional aim was to review the antihistamine product labels available in the European Union (EU) and the United States. We searched seven online databases and identified seven human lactation studies that included 25 mother-infant pairs covering cetirizine, clemastine, ebastine, epinastine, loratadine, terfenadine and triprolidine. In addition, one study investigated the impact of chlorpheniramine or promethazine on prolactin levels among 17 women, and one study investigated possible adverse drug reactions in 85 breastfed infants exposed to various antihistamines. The relative infant dose was below 5% for all antihistamines, ranging from 0.3% for terfenadine to 4.5% for clemastine. Most product labels of the 10 antihistamines with available information in both the EU and the United States reported lack of evidence and recommended to avoid use during breastfeeding. The knowledge gap on antihistamines and lactation is extensive, and further human studies are warranted to ensure optimal treatment of breastfeeding women with allergy.

Pharmacokinetics and pharmacodynamics of sildenafil in fetal lambs on extracorporeal support

BACKGROUND

Maternal transplacental administration of sildenafil is being considered for a variety of fetal conditions. Clinical translation also requires evaluation of fetal safety in a higher species, such as the fetal lamb. Experiments with the pregnant ewe are curtailed by minimal transplacental transfer as well as limited access to the fetus. The EXTra-uterine Environment for Neonatal Development (EXTEND) model renders the isolated fetal lamb readily accessible and allows for direct fetal administration of sildenafil.

METHODS

Five fetal lambs were placed on extracorporeal support in the EXTEND device and received continuous intravenous (IV) sildenafil (0.3-0.5-0.7 mg/kg/24hr) for a duration of one to seven days. Plasma sildenafil concentrations were sampled at regular intervals to establish the pharmacokinetic profile using population pharmacokinetic modeling. Serial Doppler ultrasound examination, continuous non-invasive hemodynamic monitoring and blood gas analysis were done to evaluate the pharmacodynamic effects and fetal response.

FINDINGS

The target concentration range (47-500 ng/mL) was attained with all doses. Sildenafil induced an immediate and temporary reduction of pulmonary vascular resistance, mean arterial pressure and circuit flow, without change in fetal lactate levels and acid-base status. The duration of the systemic effects increased with the dose.

INTERPRETATION

Immediate temporary pulmonary vascular and systemic hemodynamic changes induced by sildenafil were biochemically well tolerated by fetal lambs on extracorporeal support, with the 0.5 mg/kg/24 h dose balancing rapid attainment of target concentrations with short-lived systemic effects.

RESEARCH IN CONTEXT

None.

SEARCH STRATEGY BEFORE UNDERTAKING THE STUDY

A literature review was conducted searching online databases (Medline, Embase and Cochrane), using search terms: fetal OR prenatal OR antenatal AND sildenafil, without time-limit and excluding human studies. Where relevant, investigators were contacted in order to avoid duplication of work.

EVIDENCE BEFORE THIS STUDY

Prenatal therapy with sildenafil, a phosphodiesterase-5 inhibitor with vasodilatory and anti-remodeling effects on vascular smooth muscle cells, has been considered for a variety of fetal conditions. One multicenter clinical trial investigating the benefit of sildenafil in severe intrauterine growth restriction (the STRIDER-trial) was halted early due to excess mortality in the sildenafil-exposed arm at one treatment site. Such findings demonstrate the importance of extensive preclinical safety assessment in relevant animal models. Transplacentally administered sildenafil leads to decreased pulmonary arterial muscularization, preventing or reducing the occurrence of pulmonary hypertension in rat and rabbit fetuses with diaphragmatic hernia (DH). Validation of these results in a higher and relevant animal model, e.g. fetal lambs, is the next step to advance clinical translation. We recently demonstrated that, in contrast to humans, transplacental transfer of sildenafil in sheep is minimal, precluding the in vivo study of fetal effects at target concentrations using the conventional pregnant ewe model.

ADDED VALUE OF THIS STUDY

We therefore used the extracorporeal support model for fetal lambs, referred to as the EXTra-uterine Environment for Neonatal Development (EXTEND) system, bypassing placental and maternal metabolism, to investigate at what dose the target concentrations are reached, and what the fetal hemodynamic impact and response are. Fetal hemodynamic and metabolic tolerance to sildenafil are a crucial missing element on the road to clinical translation. This is therefore the first study investigating the pharmacokinetics, hemodynamic and biochemical effects of clinical-range concentrations of sildenafil in fetal lambs, free from placental and maternal interference.

IMPLICATIONS OF ALL THE AVAILABLE EVIDENCE

We demonstrated self-limiting pulmonary vasodilation, a decrease of both systemic arterial pressures and circuit flows, induced by clinical range concentrations of sildenafil, without the development of fetal acidosis. This paves the way for further investigation of prenatal sildenafil in fetal lambs on extracorporeal support. A dose of 0.5 mg/kg/24 h offered the best trade-off between rapid achievement of target concentrations and shortest duration of systemic effects. This is also the first study using the EXTEND as a model for pharmacotherapy during pregnancy.

Pharmacokinetics of 2 oral paracetamol formulations in hospitalized octogenarians

AIMS

It is currently unclear how paracetamol should be dosed in order to increase its efficacy while warranting safety in very old adults. The objective was to evaluate the pharmacokinetics of 2 oral paracetamol formulations and its metabolites in hospitalized octogenarians.

METHODS

Geriatric inpatients aged 80 years and older received a 1000-mg paracetamol tablet or granulate at 08.00, 14.00 and 20.00. After at least 4 consecutive gifts, plasma samples were collected around the 08.00 dose (trough, +0.5, +1, +2, +4, +5 and +6 h). Plasma concentrations of paracetamol and its metabolites were determined and individual pharmacokinetic parameters were derived. The Edmonton Frail Scale was used to assess frailty. An analgesic plasma target was defined as an average plasma concentration (C_avg ) of 10 mg/L.

RESULTS

The mean (±standard deviation) age was 86.78 (±4.20) years. The majority (n = 26/36, 72%) received the tablet, 10 (28%) the granulate. Thirty patients (85%) were classified with moderate to severe frailty. Seven (21%) patients had a C_avg above 10 mg/L. The median [interquartile range] time to reach the peak concentration was 50.5 [31.50-92.50] and 42.50 [33.75-106.75] min for the tablet and granulate, respectively. The coefficient of variation was 95% for time to reach the peak concentration and 30% for C_avg of paracetamol. A correlation of C_avg of paracetamol was observed with female sex and total serum bilirubin.

CONCLUSION

Large interindividual differences were found for pharmacokinetic parameters of oral paracetamol in frail inpatients after multiple dosing. Female sex and higher total serum bilirubin concentrations were associated with paracetamol exposure. No significant differences were observed between the tablet and granulate.

Bosentan Alters Endo- and Exogenous Bile Salt Disposition in Sandwich-Cultured Human Hepatocytes

Bosentan, a well-known cholestatic agent, was not identified as cholestatic at concentrations up to 200 µM based on the drug-induced cholestasis (DIC) index value, determined in a sandwich-cultured human hepatocyte (SCHH)-based DIC assay. To obtain further quantitative insights into the effects of bosentan on cellular bile salt handling by human hepatocytes, the present study determined the effect of 2.5-25 µM bosentan on endogenous bile salt levels and on the disposition of 10 µM chenodeoxycholic acid (CDCA) added to the medium in SCHHs. Bosentan reduced intracellular as well as extracellular concentrations of both endogenous glycochenodeoxycholic acid (GCDCA) and glycocholic acid in a concentration-dependent manner. When exposed to 10 µM CDCA, bosentan caused a shift from canalicular efflux to sinusoidal efflux of GCDCA. CDCA levels were not affected. Our mechanistic model confirmed the inhibitory effect of bosentan on canalicular GCDCA clearance. Moreover, our results in SCHHs also indicated reduced GCDCA formation. We confirmed the direct inhibitory effect of bosentan on CDCA conjugation with glycine in incubations with liver S9 fraction. SIGNIFICANCE STATEMENT: Bosentan was evaluated at therapeutically relevant concentrations (2.5-25 µM) in sandwich-cultured human hepatocytes. It altered bile salt disposition and inhibited canalicular secretion of glycochenodeoxycholic acid (GCDCA). Within 24 hours, bosentan caused a shift from canalicular to sinusoidal efflux of GCDCA. These results also indicated reduced GCDCA formation. This study confirmed a direct effect of bosentan on chenodeoxycholic acid conjugation with glycine in liver S9 fraction.

Determination of tacrolimus, three mono-demethylated metabolites and a M1 tautomer in human whole blood by liquid chromatography - tandem mass spectrometry

The immunosuppressant tacrolimus is the primary drug used in kidney transplantation to prevent organ rejection. A sensitive ultra-performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS) method was developed to measure tacrolimus and its three known mono-demethylated metabolites 13-O-desmethyl tacrolimus (M1), 31-O-desmethyl tacrolimus (M2), 15-O-desmethyl tacrolimus (M3). By generating the metabolites to use as standards after incubation of tacrolimus with rat liver microsomes, we discovered multiple M1 peaks which we identified as two tautomers of M1. The M1 tautomer II was also successfully validated in this method. The separation and purification of the metabolites and tautomers were performed by semi-preparative liquid chromatography with UV-detection, while confirmation was done by UPLC-MS/MS and Nuclear Magnetic Resonance. For quantification an easy sample preparation was performed with zinc sulfate and acetonitrile as cell lyses and precipitation. Detection was performed in positive electrospray ionization. By better characterization of the metabolites and the tautomers, we could possibly explain insight into the clinical condition and thus adjust the immunosuppressant therapy individually per patient. Calibration curves were linear for all compounds. Precision was assessed according to the NCCLS EP5-T guideline, being below 15 % and mean recoveries were between 93 and 110 % for tacrolimus, its three metabolites and the M1 tautomer II. The validated method was successfully applied in a cohort of 20 patients after kidney transplantation.

A Population Pharmacokinetic Modeling and Simulation Study of Posaconazole Oral Suspension in Immunocompromised Pediatric Patients: A Short Communication

BACKGROUND

Posaconazole oral suspension emerged as a promising candidate for prophylaxis of invasive fungal infections in immunocompromised children. Its pharmacodynamic advantages include a broad-spectrum activity and a favorable safety profile; however, they are overshadowed by its large pharmacokinetic (PK) variability, which might cause subtherapeutic exposure. The aim of this study was to develop a population (pop) PK model based on rich sampling data to better understand the PK of posaconazole oral suspension in pediatric patients.

METHODS

Data were obtained from a prospective interventional study involving hospitalized pediatric patients with a hematologic malignancy and prophylactically treated with posaconazole oral suspension. After constructing the popPK model, the probability of target attainment (PTA; 100% T ≥ 0.7 mg/L) for prophylaxis under fixed, body weight-based, and body surface area-based dosing was evaluated using Monte Carlo simulation.

RESULTS

Fourteen patients contributed 112 posaconazole plasma concentrations. The PK of posaconazole was adequately described by a 1-compartment model with lag time 2.71 hours [13%]; nonlinear bioavailability ED50 99.1 mg/m2 (fixed); first-order absorption rate constant 0.325 hour-1 [27%]; apparent volume of distribution 1150 L [34%]; and apparent clearance 15.4 L/h [24%] (∼70-kg individual). The bioavailability decreased in the presence of diarrhea and co-treatment with a proton pump inhibitor (PPI). The unexplained interindividual variability in posaconazole PK remained large. The PTA was <85%, irrespective of the simulated dosing strategy. Patients without diarrhea and not administered a PPI had the highest PTA (85% under the fixed 300-mg dosing 4 times per day).

CONCLUSIONS

Therapeutic drug monitoring is recommended during prophylactic posaconazole therapy in immunocompromised pediatric patients. Large-scale comparative studies are needed to characterize the PK variability between different posaconazole formulations in this cohort.

Quantification and Explanation of the Variability of First-Dose Amikacin Concentrations in Critically Ill Patients Admitted to the Emergency Department: A Population Pharmacokinetic Analysis

BACKGROUND

There may be a difference between the determinants of amikacin exposure in emergency department (ED) versus intensive care (ICU) patients, and the peak amikacin concentration varies widely between patients. Moreover, when the first dose of antimicrobials is administered to septic patients admitted to the ED, fluid resuscitation and vasopressors have just been initiated. Nevertheless, population pharmacokinetic modelling data for amikacin in ED patients are unavailable.

OBJECTIVE

The aim of this study was to quantify the interindividual variability (IIV) in the pharmacokinetics of amikacin in patients admitted to the ED and to identify the patient characteristics that explain this IIV.

METHODS

Patients presenting at the ED with severe sepsis or septic shock were randomly assigned to receive amikacin 25 mg/kg or 15 mg/kg intravenously. Blood samples were collected at 1, 6 and 24 h after the onset of the first amikacin infusion. Data were analysed using nonlinear mixed-effects modelling.

RESULTS

A two-compartment population pharmacokinetic model was developed based on 279 amikacin concentrations from 97 patients. The IIV in clearance (CL) and central distribution volume (V₁) were 71% and 26%, respectively. Body mass index (BMI), serum total protein level, serum sodium level, and fluid balance 24 h after amikacin administration explained 30% of the IIV in V₁, leaving 18% of the IIV unexplained. BMI and creatinine clearance according to the Chronic Kidney Disease Epidemiology Collaboration (CKD-EPI) equation 24 h after amikacin administration explained 46% of the IIV in CL, and 39% remained unexplained.

CONCLUSION

The IIV of amikacin pharmacokinetics in ED patients is large. Higher doses may be considered in patients with low serum sodium levels, low total protein levels, or a high fluid balance.

TRIAL REGISTRATION

ClinicalTrials.gov ID: NCT02365272.

Pharmacokinetics in Zebrafish Embryos (ZFE) Following Immersion and Intrayolk Administration: A Fluorescence-Based Analysis

Zebrafish embryos (ZFE) have increasingly gained in popularity as a model to perform safety screenings of compounds. Although immersion of ZFE is the main route of exposure used, evidence shows that not all small molecules are equally absorbed, possibly resulting in false-negative readouts and incorrect conclusions. In this study, we compared the pharmacokinetics of seven fluorescent compounds with known physicochemical properties that were administered to two-cell stage embryos by immersion or by IY microinjection. Absorption and distribution of the dyes were followed at various timepoints up to 120 hpf by spatiotemporal fluorescence imaging. The concentration (10 µM) and dose (2 mg/kg) used were selected as quantities typically applied in preclinical experiments and zebrafish studies. The data show that in the case of a lipophilic compound (log D: 1.73) the immersion procedure resulted in an intrabody exposure which is similar or higher than that seen after the IY microinjection. In contrast, zero to low intrabody exposure was reached after immersion of the embryos with less lipophilic compounds. In the latter case IY microinjection, a technical procedure that can be easily automated, is highly recommended.

Meropenem Pharmacokinetics and Target Attainment in Critically Ill Patients Are Not Affected by Extracorporeal Membrane Oxygenation: A Matched Cohort Analysis

Existing evidence is inconclusive whether meropenem dosing should be adjusted in patients receiving extracorporeal membrane oxygenation (ECMO). Therefore, the aim of this observational matched cohort study was to evaluate the effect of ECMO on pharmacokinetic (PK) variability and target attainment (TA) of meropenem. Patients admitted to the intensive care unit (ICU) simultaneously treated with meropenem and ECMO were eligible. Patients were matched 1:1, based on renal function and body weight, with non-ECMO ICU patients. Meropenem blood sampling was performed over one or two dosing intervals. Population PK modelling was performed using NONMEM7.5. TA was defined as free meropenem concentrations >2 or 8 mg/L (i.e., 1 or 4× minimal inhibitory concentration, respectively) throughout the whole dosing interval. In total, 25 patients were included, contributing 27 dosing intervals. The overall TA was 56% and 26% for the 2 mg/L and 8 mg/L target, respectively. Population PK modelling identified estimated glomerular filtration rate according to the Chronic Kidney Disease Epidemiology equation and body weight, but not ECMO, as significant predictors. In conclusion, TA of meropenem was confirmed to be poor under standard dosing in critically ill patients but was not found to be influenced by ECMO. Future studies should focus on applying dose optimisation strategies for meropenem based on renal function, regardless of ECMO.

Spatiotemporal imaging and pharmacokinetics of fluorescent compounds in zebrafish eleuthero-embryos after different routes of administration

Zebrafish (Danio rerio) is increasingly used to assess the pharmacological activity and toxicity of compounds. The spatiotemporal distribution of seven fluorescent alkyne compounds was examined during 48 h after immersion (10 µM) or microinjection (2 mg/kg) in the pericardial cavity (PC), intraperitoneally (IP) and yolk sac (IY) of 3 dpf zebrafish eleuthero-embryos. By modelling the fluorescence of whole-body contours present in fluorescence images, the main pharmacokinetic (PK) parameter values of the compounds were determined. It was demonstrated that especially in case of short incubations (1-3 h) immersion can result in limited intrabody exposure to compounds. In this case, PC and IP microinjections represent excellent alternatives. Significantly, IY microinjections did not result in a suitable intrabody distribution of the compounds. Performing a QSPkR (quantitative structure-pharmacokinetic relationship) analysis, LogD was identified as the only molecular descriptor that explains the final uptake of the selected compounds. It was also shown that combined administration of compounds (immersion and microinjection) provides a more stable intrabody exposure, at least in case of a prolonged immersion and compounds with LogD value > 1. These results will help reduce the risk of false negative results and can offer an invaluable input for future translational research and safety assessment applications.

Creatinine Trends and Patterns in Neonates Undergoing Whole Body Hypothermia: A Systematic Review

Many neonates undergoing whole body hypothermia (WBH) following moderate to severe perinatal asphyxia may also suffer from renal impairment. While recent data suggest WBH-related reno-protection, differences in serum creatinine (Scr) patterns to reference patterns were not yet reported. We therefore aimed to document Scr trends and patterns in asphyxiated neonates undergoing WBH and compared these to centiles from a reference Scr data set of non-asphyxiated (near)term neonates. Using a systematic review strategy, reports on Scr trends (mean ± SD, median or interquartile range) were collected (day 1-7) in WBH cohorts and compared to centiles of an earlier reported reference cohort of non-asphyxia cases. Based on 13 papers on asphyxia + WBH cases, a pattern of postnatal Scr trends in asphyxia + WBH cases was constructed. Compared to the reference 50th centile Scr values, mean or median Scr values at birth and up to 48 h were higher in asphyxia + WBH cases with a subsequent uncertain declining trend towards, at best, high or high-normal creatinine values afterwards. Such patterns are valuable for anticipating average changes in renal drug clearance but do not yet cover the relevant inter-patient variability observed in WBH cases, as this needs pooling of individual Screa profiles, preferably beyond the first week of life.

Pharmacokinetic/Pharmacodynamic Target Attainment Based on Measured versus Predicted Unbound Ceftriaxone Concentrations in Critically Ill Patients with Pneumonia: An Observational Cohort Study

The impact of ceftriaxone pharmacokinetic alterations on protein binding and PK/PD target attainment still remains unclear. We evaluated pharmacokinetic/pharmacodynamic (PK/PD) target attainment of unbound ceftriaxone in critically ill patients with severe community-acquired pneumonia (CAP). Besides, we evaluated the accuracy of predicted vs. measured unbound ceftriaxone concentrations, and its impact on PK/PD target attainment. A prospective observational cohort study was carried out in adult patients admitted to the intensive care unit with severe CAP. Ceftriaxone 2 g q24h intermittent infusion was administered to all patients. Successful PK/PD target attainment was defined as unbound trough concentrations above 1 or 4 mg/L throughout the whole dosing interval. Acceptable overall PK/PD target attainment was defined as successful target attainment in ≥90% of all dosing intervals. Measured unbound ceftriaxone concentrations (CEFu) were compared to unbound concentrations predicted from various protein binding models. Thirty-one patients were included. The 1 mg/L and 4 mg/L targets were reached in 26/32 (81%) and 15/32 (47%) trough samples, respectively. Increased renal function was associated with the failure to attain both PK/PD targets. Unbound ceftriaxone concentrations predicted by the protein binding model developed in the present study showed acceptable bias and precision and had no major impact on PK/PD target attainment. We showed suboptimal (i.e., <90%) unbound ceftriaxone PK/PD target attainment when using a standard 2 g q24h dosing regimen in critically ill patients with severe CAP. Renal function was the major driver for the failure to attain the predefined targets, in accordance with results found in general and septic ICU patients. Interestingly, CEFu was reliably predicted from CEFt without major impact on clinical decisions regarding PK/PD target attainment. This suggests that, when carefully selecting a protein binding model, CEFu does not need to be measured. As a result, the turn-around time and cost for ceftriaxone quantification can be substantially reduced.

Identification of novel inhibitors of rat Mrp3

Multidrug resistance-associated protein (MRP; ABCC gene family) mediated efflux transport plays an important role in the systemic and tissue exposure profiles of many drugs and their metabolites, and also of endogenous compounds like bile acids and bilirubin conjugates. However, potent and isoform-selective inhibitors of the MRP subfamily are currently lacking. Therefore, the purpose of the present work was to identify novel rat Mrp3 inhibitors. Using 5(6)-carboxy-2',7'-dichlorofluorescein diacetate (CDFDA) as a model-(pro)substrate for Mrp3 in an oil-spin assay with primary rat hepatocytes, the extent of inhibition of CDF efflux was determined for 1584 compounds, yielding 59 hits (excluding the reference inhibitor) that were identified as new Mrp3 inhibitors. A naive Bayesian prediction model was constructed in Pipeline Pilot to elucidate physicochemical and structural features of compounds causing Mrp3 inhibition. The final Bayesian model generated common physicochemical properties of Mrp3 inhibitors. For instance, more than half of the hits contain a phenolic structure. The identified compounds have an AlogP between 2 and 4.5, between 5 to 8 hydrogen bond acceptor atoms, a molecular weight between 260 and 400, and 2 or more aromatic rings. Compared to the depleted dataset (i.e. 90% remaining compounds), the Mrp3 hit rate in the enriched set was 7.5-fold higher (i.e. 17.2% versus 2.3%). Several hits from this first screening approach were confirmed in an additional study using Mrp3 transfected inside-out membrane vesicles. In conclusion, several new and potent inhibitors of Mrp3 mediated efflux were identified in an optimized in vitro rat hepatocyte assay and confirmed using Mrp3 transfected inside-out membrane vesicles. A final naive Bayesian model was developed in an iterative way to reveal common physicochemical and structural features for Mrp3 inhibitors. The final Bayesian model will enable in silico screening of larger libraries and in vitro identification of more potent Mrp3 inhibitors.

Biomarkers of cholestasis

Cholestasis is a major pathological manifestation, often resulting in detrimental liver conditions, which occurs in a variety of indications collectively termed cholestatic liver diseases. The frequent asymptomatic character and complexity of cholestasis, together with the lack of a straightforward biomarker, hampers early detection and treatment of the condition. The 'omics' era, however, has resulted in a plethora of cholestatic indicators, yet a single clinically applicable biomarker for a given cholestatic disease remains missing. The criteria to fulfil as an ideal biomarker as well as the challenging molecular pathways in cholestatic liver diseases advocate for a scenario in which multiple biomarkers, originating from different domains, will be assessed concomitantly. This review gives an overview of classical clinical and novel molecular biomarkers in cholestasis, focusing on their benefits and drawbacks.

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.

A comprehensive review on non-clinical methods to study transfer of medication into breast milk - A contribution from the ConcePTION project

Breastfeeding plays a major role in the health and wellbeing of mother and infant. However, information on the safety of maternal medication during breastfeeding is lacking for most medications. This leads to discontinuation of either breastfeeding or maternal therapy, although many medications are likely to be safe. Since human lactation studies are costly and challenging, validated non-clinical methods would offer an attractive alternative. This review gives an extensive overview of the non-clinical methods (in vitro, in vivo and in silico) to study the transfer of maternal medication into the human breast milk, and subsequent neonatal systemic exposure. Several in vitro models are available, but model characterization, including quantitative medication transport data across the in vitro blood-milk barrier, remains rather limited. Furthermore, animal in vivo models have been used successfully in the past. However, these models don't always mimic human physiology due to species-specific differences. Several efforts have been made to predict medication transfer into the milk based on physicochemical characteristics. However, the role of transporter proteins and several physiological factors (e.g., variable milk lipid content) are not accounted for by these methods. Physiologically-based pharmacokinetic (PBPK) modelling offers a mechanism-oriented strategy with bio-relevance. Recently, lactation PBPK models have been reported for some medications, showing at least the feasibility and value of PBPK modelling to predict transfer of medication into the human milk. However, reliable data as input for PBPK models is often missing. The iterative development of in vitro, animal in vivo and PBPK modelling methods seems to be a promising approach. Human in vitro models will deliver essential data on the transepithelial transport of medication, whereas the combination of animal in vitro and in vivo methods will deliver information to establish accurate in vitro/in vivo extrapolation (IVIVE) algorithms and mechanistic insights. Such a non-clinical platform will be developed and thoroughly evaluated by the Innovative Medicines Initiative ConcePTION.

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.

The Neonatal and Juvenile Pig in Pediatric Drug Discovery and Development

Pharmacotherapy in pediatric patients is challenging in view of the maturation of organ systems and processes that affect pharmacokinetics and pharmacodynamics. Especially for the youngest age groups and for pediatric-only indications, neonatal and juvenile animal models can be useful to assess drug safety and to better understand the mechanisms of diseases or conditions. In this respect, the use of neonatal and juvenile pigs in the field of pediatric drug discovery and development is promising, although still limited at this point. This review summarizes the comparative postnatal development of pigs and humans and discusses the advantages of the juvenile pig in view of developmental pharmacology, pediatric diseases, drug discovery and drug safety testing. Furthermore, limitations and unexplored aspects of this large animal model are covered. At this point in time, the potential of the neonatal and juvenile pig as nonclinical safety models for pediatric drug development is underexplored.