A quantitative modeling framework to understand the physiology of the hypothalamic-pituitary-adrenal axis and interaction with cortisol replacement therapy.

Congenital adrenal hyperplasia (CAH) is characterized by impaired adrenal cortisol production. Hydrocortisone (synthetic cortisol) is the drug-of-choice for cortisol replacement therapy, aiming to mimic physiological cortisol circadian rhythm. The hypothalamic-pituitary-adrenal (HPA) axis controls cortisol production through the pituitary adrenocorticotropic hormone (ACTH) and feedback mechanisms. The aim of this study was to quantify key mechanisms involved in the HPA axis activity regulation and their interaction with hydrocortisone therapy. Data from 30 healthy volunteers was leveraged: Endogenous ACTH and cortisol concentrations without any intervention as well as cortisol concentrations measured after dexamethasone suppression and single dose administration of (i) 0.5-10 mg hydrocortisone as granules, (ii) 20 mg hydrocortisone as granules and intravenous bolus. A stepwise model development workflow was used: A newly developed model for endogenous ACTH and cortisol was merged with a refined hydrocortisone pharmacokinetic model. The joint model was used to simulate ACTH and cortisol trajectories in CAH patients with varying degrees of enzyme deficiency, with or without hydrocortisone administration, and healthy individuals. Time-dependent ACTH-driven endogenous cortisol production and cortisol-mediated feedback inhibition of ACTH secretion processes were quantified and implemented in the model. Comparison of simulated ACTH and cortisol trajectories between CAH patients and healthy individuals showed the importance of administering hydrocortisone before morning ACTH secretion peak time to suppress ACTH overproduction observed in untreated CAH patients. The developed framework allowed to gain insights on the physiological mechanisms of the HPA axis regulation, its perturbations in CAH and interaction with hydrocortisone administration, paving the way towards cortisol replacement therapy optimization.

Population pharmacokinetics of TLD-1, a novel liposomal doxorubicin, in a phase I trial.

STUDY OBJECTIVES: TLD-1 is a novel pegylated liposomal doxorubicin (PLD) formulation aiming to optimise the PLD efficacy-toxicity ratio. We aimed to characterise TLD-1's population pharmacokinetics using non-compartmental analysis and nonlinear mixed-effects modelling. METHODS: The PK of TLD-1 was analysed by performing a non-compartmental analysis of longitudinal doxorubicin plasma concentration measurements obtained from a clinical trial in 30 patients with advanced solid tumours across a 4.5-fold dose range. Furthermore, a joint parent-metabolite PK model of doxorubicinentrapped, doxorubicinfree, and metabolite doxorubicinol was developed. Interindividual and interoccasion variability around the typical PK parameters and potential covariates to explain parts of this variability were explored. RESULTS: Medians  ± standard deviations of dose-normalised doxorubicinentrapped+free Cmax and AUC0-∞ were 0.342 ± 0.134 mg/L and 40.1 ± 18.9 mg·h/L, respectively. The median half-life (95 h) was 23.5 h longer than the half-life of currently marketed PLD. The novel joint parent-metabolite model comprised a one-compartment model with linear release (doxorubicinentrapped), a two-compartment model with linear elimination (doxorubicinfree), and a one-compartment model with linear elimination for doxorubicinol. Body surface area on the volumes of distribution for free doxorubicin was the only significant covariate. CONCLUSION: The population PK of TLD-1, including its release and main metabolite, were successfully characterised using non-compartmental and compartmental analyses. Based on its long half-life, TLD-1 presents a promising candidate for further clinical development. The PK characteristics form the basis to investigate TLD-1 exposure-response (i.e., clinical efficacy) and exposure-toxicity relationships in the future. Once such relationships have been established, the developed population PK model can be further used in model-informed precision dosing strategies. CLINICAL TRIAL REGISTRATION: ClinicalTrials.gov-NCT03387917-January 2, 2018.

Personalized Dosing of Medicines for Children: A Primer on Pediatric Pharmacometrics for Clinicians.

The widespread use of drugs for unapproved purposes remains common in children, primarily attributable to practical, ethical, and financial constraints associated with pediatric drug research. Pharmacometrics, the scientific discipline that involves the application of mathematical models to understand and quantify drug effects, holds promise in advancing pediatric pharmacotherapy by expediting drug development, extending applications, and personalizing dosing. In this review, we delineate the principles of pharmacometrics, and explore its clinical applications and prospects. The fundamental aspect of any pharmacometric analysis lies in the selection of appropriate methods for quantifying pharmacokinetics and pharmacodynamics. Population pharmacokinetic modeling is a data-driven method ('top-down' approach) to approximate population-level pharmacokinetic parameters, while identifying factors contributing to inter-individual variability. Model-informed precision dosing is increasingly used to leverage population pharmacokinetic models and patient data, to formulate individualized dosing recommendations. Physiologically based pharmacokinetic models integrate physicochemical drug properties with biological parameters ('bottom-up approach'), and is particularly valuable in situations with limited clinical data, such as early drug development, assessing drug-drug interactions, or adapting dosing for patients with specific comorbidities. The effective implementation of these complex models hinges on strong collaboration between clinicians and pharmacometricians, given the pivotal role of data availability. Promising advancements aimed at improving data availability encompass innovative techniques such as opportunistic sampling, minimally invasive sampling approaches, microdialysis, and in vitro investigations. Additionally, ongoing research efforts to enhance measurement instruments for evaluating pharmacodynamics responses, including biomarkers and clinical scoring systems, are expected to significantly bolster our capacity to understand drug effects in children.

Nonlinear Mixed-Effects Model of Z-Endoxifen Concentrations in Tamoxifen-Treated Patients from the CEPAM Cohort.

Tamoxifen is widely used in patients with hormone receptor-positive breast cancer. The polymorphic enzyme CYP2D6 is primarily responsible for metabolic activation of tamoxifen, resulting in substantial interindividual variability of plasma concentrations of its most important metabolite, Z-endoxifen. The Z-endoxifen concentration thresholds below which tamoxifen treatment is less efficacious have been proposed but not validated, and prospective trials of individualized tamoxifen treatment to achieve Z-endoxifen concentration thresholds are considered infeasible. Therefore, we aim to validate the association between Z-endoxifen concentration and tamoxifen treatment outcomes, and identify a Z-endoxifen concentration threshold of tamoxifen efficacy, using pharmacometric modeling and simulation. As a first step, the CYP2D6 Endoxifen Percentage Activity Model (CEPAM) cohort was created by pooling data from 28 clinical studies (> 7,000 patients) with measured endoxifen plasma concentrations. After cleaning, data from 6,083 patients were used to develop a nonlinear mixed-effect (NLME) model for tamoxifen and Z-endoxifen pharmacokinetics that includes a conversion factor to allow inclusion of studies that measured total endoxifen but not Z-endoxifen. The final parent-metabolite NLME model confirmed the primary role of CYP2D6, and contributions from body weight, CYP2C9 phenotype, and co-medication with CYP2D6 inhibitors, on Z-endoxifen pharmacokinetics. Future work will use the model to simulate Z-endoxifen concentrations in patients receiving single agent tamoxifen treatment within large prospective clinical trials with long-term survival to identify the Z-endoxifen concentration threshold below which tamoxifen is less efficacious. Identification of this concentration threshold would allow personalized tamoxifen treatment to improve outcomes in patients with hormone receptor-positive breast cancer.

Bayesian framework for multi-source data integration-Application to human extrapolation from preclinical studies.

In preclinical investigations, for example, in in vitro, in vivo, and in silico studies, the pharmacokinetic, pharmacodynamic, and toxicological characteristics of a drug are evaluated before advancing to first-in-man trial. Usually, each study is analyzed independently and the human dose range does not leverage the knowledge gained from all studies. Taking into account all preclinical data through inferential procedures can be particularly interesting in obtaining a more precise and reliable starting dose and dose range. Our objective is to propose a Bayesian framework for multi-source data integration, customizable, and tailored to the specific research question. We focused on preclinical results extrapolated to humans, which allowed us to predict the quantities of interest (e.g. maximum tolerated dose, etc.) in humans. We build an approach, divided into four steps, based on a sequential parameter estimation for each study, extrapolation to human, commensurability checking between posterior distributions and final information merging to increase the precision of estimation. The new framework is evaluated via an extensive simulation study, based on a real-life example in oncology. Our approach allows us to better use all the information compared to a standard framework, reducing uncertainty in the predictions and potentially leading to a more efficient dose selection.

C-Reactive Protein as an Early Predictor of Efficacy in Advanced Non-Small-Cell Lung Cancer Patients: A Tumor Dynamics-Biomarker Modeling Framework.

In oncology, longitudinal biomarkers reflecting the patient's status and disease evolution can offer reliable predictions of the patient's response to treatment and prognosis. By leveraging clinical data in patients with advanced non-small-cell lung cancer receiving first-line chemotherapy, we aimed to develop a framework combining anticancer drug exposure, tumor dynamics (RECIST criteria), and C-reactive protein (CRP) concentrations, using nonlinear mixed-effects models, to evaluate and quantify by means of parametric time-to-event models the significance of early longitudinal predictors of progression-free survival (PFS) and overall survival (OS). Tumor dynamics was characterized by a tumor size (TS) model accounting for anticancer drug exposure and development of drug resistance. CRP concentrations over time were characterized by a turnover model. An x-fold change in TS from baseline linearly affected CRP production. CRP concentration at treatment cycle 3 (day 42) and the difference between CRP concentration at treatment cycles 3 and 2 were the strongest predictors of PFS and OS. Measuring longitudinal CRP allows for the monitoring of inflammatory levels and, along with its reduction across treatment cycles, presents a promising prognostic marker. This framework could be applied to other treatment modalities such as immunotherapies or targeted therapies allowing the timely identification of patients at risk of early progression and/or short survival to spare them unnecessary toxicities and provide alternative treatment decisions.

Towards Model-Informed Precision Dosing of Voriconazole: Challenging Published Voriconazole Nonlinear Mixed-Effects Models with Real-World Clinical Data.

BACKGROUND AND OBJECTIVES: Model-informed precision dosing (MIPD) frequently uses nonlinear mixed-effects (NLME) models to predict and optimize therapy outcomes based on patient characteristics and therapeutic drug monitoring data. MIPD is indicated for compounds with narrow therapeutic range and complex pharmacokinetics (PK), such as voriconazole, a broad-spectrum antifungal drug for prevention and treatment of invasive fungal infections. To provide guidance and recommendations for evidence-based application of MIPD for voriconazole, this work aimed to (i) externally evaluate and compare the predictive performance of a published so-called 'hybrid' model for MIPD (an aggregate model comprising features and prior information from six previously published NLME models) versus two 'standard' NLME models of voriconazole, and (ii) investigate strategies and illustrate the clinical impact of Bayesian forecasting for voriconazole. METHODS: A workflow for external evaluation and application of MIPD for voriconazole was implemented. Published voriconazole NLME models were externally evaluated using a comprehensive in-house clinical database comprising nine voriconazole studies and prediction-/simulation-based diagnostics. The NLME models were applied using different Bayesian forecasting strategies to assess the influence of prior observations on model predictivity. RESULTS: The overall best predictive performance was obtained using the aggregate model. However, all NLME models showed only modest predictive performance, suggesting that (i) important PK processes were not sufficiently implemented in the structural submodels, (ii) sources of interindividual variability were not entirely captured, and (iii) interoccasion variability was not adequately accounted for. Predictive performance substantially improved by including the most recent voriconazole observations in MIPD. CONCLUSION: Our results highlight the potential clinical impact of MIPD for voriconazole and indicate the need for a comprehensive (pre-)clinical database as basis for model development and careful external model evaluation for compounds with complex PK before their successful use in MIPD.

Model-Informed Target Morning 17α-Hydroxyprogesterone Concentrations in Dried Blood Spots for Pediatric Congenital Adrenal Hyperplasia Patients.

Monitoring cortisol replacement therapy in congenital adrenal hyperplasia (CAH) patients is vital to avoid serious adverse events such as adrenal crises due to cortisol underexposure or metabolic consequences due to cortisol overexposure. The less invasive dried blood spot (DBS) sampling is an advantageous alternative to traditional plasma sampling, especially in pediatric patients. However, target concentrations for important disease biomarkers such as 17α-hydroxyprogesterone (17-OHP) are unknown using DBS. Therefore, a modeling and simulation framework, including a pharmacokinetic/pharmacodynamic model linking plasma cortisol concentrations to DBS 17-OHP concentrations, was used to derive a target morning DBS 17-OHP concentration range of 2-8 nmol/L in pediatric CAH patients. Since either capillary or venous DBS sampling is becoming more common in the clinics, the clinical applicability of this work was shown by demonstrating the comparability of capillary and venous cortisol and 17-OHP concentrations collected by DBS sampling, using a Bland-Altman and Passing-Bablok analysis. The derived target morning DBS 17-OHP concentration range is a first step towards providing improved therapy monitoring using DBS sampling and adjusting hydrocortisone (synthetic cortisol) dosing in children with CAH. In the future, this framework can be used to assess further research questions, e.g., target replacement ranges for the entire day.

Ulcerative Colitis and Acute Severe Ulcerative Colitis Patients Are Overlooked in Infliximab Population Pharmacokinetic Models: Results from a Comprehensive Review.

Ulcerative colitis (UC) is part of the inflammatory bowels diseases, and moderate to severe UC patients can be treated with anti-tumour necrosis α monoclonal antibodies, including infliximab (IFX). Even though treatment of UC patients by IFX has been in place for over a decade, many gaps in modelling of IFX PK in this population remain. This is even more true for acute severe UC (ASUC) patients for which early prediction of IFX pharmacokinetic (PK) could highly improve treatment outcome. Thus, this review aims to compile and analyse published population PK models of IFX in UC and ASUC patients, and to assess the current knowledge on disease activity impact on IFX PK. For this, a semi-systematic literature search was conducted, from which 26 publications including a population PK model analysis of UC patients receiving IFX therapy were selected. Amongst those, only four developed a model specifically for UC patients, and only three populations included severe UC patients. Investigations of disease activity impact on PK were reported in only 4 of the 14 models selected. In addition, the lack of reported model codes and assessment of predictive performance make the use of published models in a clinical setting challenging. Thus, more comprehensive investigation of PK in UC and ASUC is needed as well as more adequate reports on developed models and their evaluation in order to apply them in a clinical setting.

Quantification of the Time Course of CYP3A Inhibition, Activation, and Induction Using a Population Pharmacokinetic Model of Microdosed Midazolam Continuous Infusion.

BACKGROUND: Cytochrome P450 (CYP) 3A contributes to the metabolism of many approved drugs. CYP3A perpetrator drugs can profoundly alter the exposure of CYP3A substrates. However, effects of such drug-drug interactions are usually reported as maximum effects rather than studied as time-dependent processes. Identification of the time course of CYP3A modulation can provide insight into when significant changes to CYP3A activity occurs, help better design drug-drug interaction studies, and manage drug-drug interactions in clinical practice. OBJECTIVE: We aimed to quantify the time course and extent of the in vivo modulation of different CYP3A perpetrator drugs on hepatic CYP3A activity and distinguish different modulatory mechanisms by their time of onset, using pharmacologically inactive intravenous microgram doses of the CYP3A-specific substrate midazolam, as a marker of CYP3A activity. METHODS: Twenty-four healthy individuals received an intravenous midazolam bolus followed by a continuous infusion for 10 or 36 h. Individuals were randomized into four arms: within each arm, two individuals served as a placebo control and, 2 h after start of the midazolam infusion, four individuals received the CYP3A perpetrator drug: voriconazole (inhibitor, orally or intravenously), rifampicin (inducer, orally), or efavirenz (activator, orally). After midazolam bolus administration, blood samples were taken every hour (rifampicin arm) or every 15 min (remaining study arms) until the end of midazolam infusion. A total of 1858 concentrations were equally divided between midazolam and its metabolite, 1'-hydroxymidazolam. A nonlinear mixed-effects population pharmacokinetic model of both compounds was developed using NONMEM®. CYP3A activity modulation was quantified over time, as the relative change of midazolam clearance encountered by the perpetrator drug, compared to the corresponding clearance value in the placebo arm. RESULTS: Time course of CYP3A modulation and magnitude of maximum effect were identified for each perpetrator drug. While efavirenz CYP3A activation was relatively fast and short, reaching a maximum after approximately 2-3 h, the induction effect of rifampicin could only be observed after 22 h, with a maximum after approximately 28-30 h followed by a steep drop to almost baseline within 1-2 h. In contrast, the inhibitory impact of both oral and intravenous voriconazole was prolonged with a steady inhibition of CYP3A activity followed by a gradual increase in the inhibitory effect until the end of sampling at 8 h. Relative maximum clearance changes were +59.1%, +46.7%, -70.6%, and -61.1% for efavirenz, rifampicin, oral voriconazole, and intravenous voriconazole, respectively. CONCLUSIONS: We could distinguish between different mechanisms of CYP3A modulation by the time of onset. Identification of the time at which clearance significantly changes, per perpetrator drug, can guide the design of an optimal sampling schedule for future drug-drug interaction studies. The impact of a short-term combination of different perpetrator drugs on the paradigm CYP3A substrate midazolam was characterized and can define combination intervals in which no relevant interaction is to be expected. CLINICAL TRIAL REGISTRATION: The trial was registered at the European Union Drug Regulating Authorities for Clinical Trials (EudraCT-No. 2013-004869-14).

Microdialysis of Voriconazole and its N-Oxide Metabolite: Amalgamating Knowledge of Distribution and Metabolism Processes in Humans.

PURPOSE: Voriconazole is an essential antifungal drug whose complex pharmacokinetics with high interindividual variability impedes effective and safe therapy. By application of the minimally-invasive sampling technique microdialysis, interstitial space fluid (ISF) concentrations of VRC and its potentially toxic N-oxide metabolite (NO) were assessed to evaluate target-site exposure for further elucidating VRC pharmacokinetics. METHODS: Plasma and ISF samples of a clinical trial with an approved VRC dosing regimen were analyzed for VRC and NO concentrations. Concentration-time profiles, exposure assessed as area-under-the-curve (AUC) and metabolic ratios of four healthy adults in plasma and ISF were evaluated regarding the impact of multiple dosing and CYP2C19 genotype. RESULTS: VRC and NO revealed distribution into ISF with AUC values being ≤2.82- and 17.7-fold lower compared to plasma, respectively. Intraindividual variability of metabolic ratios was largest after the first VRC dose administration while interindividual variability increased with multiple dosing. The CYP2C19 genotype influenced interindividual differences with a maximum 6- and 24-fold larger AUCNO/AUCVRC ratio between the intermediate and rapid metabolizer in plasma and ISF, respectively. VRC metabolism was saturated/auto-inhibited indicated by substantially decreasing metabolic concentration ratios with increasing VRC concentrations and after multiple dosing. CONCLUSION: The feasibility of the simultaneous microdialysis of VRC and NO in vivo was demonstrated and provided new quantitative insights by leveraging distribution and metabolism processes of VRC in humans. The exploratory analysis suggested substantial dissimilarities of VRC and NO pharmacokinetics in plasma and ISF. Ultimately, a thorough understanding of target-site pharmacokinetics might contribute to the optimization of personalized VRC dosing regimens.

The Impact of Low Cardiac Output on Propofol Pharmacokinetics across Age Groups-An Investigation Using Physiologically Based Pharmacokinetic Modelling.

BACKGROUND: pathophysiological changes such as low cardiac output (LCO) impact pharmacokinetics, but its extent may be different throughout pediatrics compared to adults. Physiologically based pharmacokinetic (PBPK) modelling enables further exploration. METHODS: A validated propofol model was used to simulate the impact of LCO on propofol clearance across age groups using the PBPK platform, Simcyp® (version 19). The hepatic and renal extraction ratio of propofol was then determined in all age groups. Subsequently, manual infusion dose explorations were conducted under LCO conditions, targeting a 3 µg/mL (80-125%) propofol concentration range. RESULTS: Both hepatic and renal extraction ratios increased from neonates, infants, children to adolescents and adults. The relative change in clearance following CO reductions increased with age, with the least impact of LCO in neonates. The predicted concentration remained within the 3 µg/mL (80-125%) range under normal CO and LCO (up to 30%) conditions in all age groups. When CO was reduced by 40-50%, a dose reduction of 15% is warranted in neonates, infants and children, and 25% in adolescents and adults. CONCLUSIONS: PBPK-driven, the impact of reduced CO on propofol clearance is predicted to be age-dependent, and proportionally greater in adults. Consequently, age group-specific dose reductions for propofol are required in LCO conditions.

Microdialysis of Drug and Drug Metabolite: a Comprehensive In Vitro Analysis for Voriconazole and Voriconazole N-oxide.

PURPOSE: Voriconazole is a therapeutically challenging antifungal drug associated with high interindividual pharmacokinetic variability. As a prerequisite to performing clinical trials using the minimally-invasive sampling technique microdialysis, a comprehensive in vitro microdialysis characterization of voriconazole (VRC) and its potentially toxic N-oxide metabolite (NO) was performed. METHODS: The feasibility of simultaneous microdialysis of VRC and NO was explored in vitro by investigating the relative recovery (RR) of both compounds in the absence and presence of the other. The dependency of RR on compound combination, concentration, microdialysis catheter and study day was evaluated and quantified by linear mixed-effects modeling. RESULTS: Median RR of VRC and NO during individual microdialysis were high (87.6% and 91.1%). During simultaneous microdialysis of VRC and NO, median RR did not change (87.9% and 91.1%). The linear mixed-effects model confirmed the absence of significant differences between RR of VRC and NO during individual and simultaneous microdialysis as well as between the two compounds (p > 0.05). No concentration dependency of RR was found (p = 0.284). The study day was the main source of variability (46.3%) while the microdialysis catheter only had a minor effect (4.33%). VRC retrodialysis proved feasible as catheter calibration for both compounds. CONCLUSION: These in vitro microdialysis results encourage the application of microdialysis in clinical trials to assess target-site concentrations of VRC and NO. This can support the generation of a coherent understanding of VRC pharmacokinetics and its sources of variability. Ultimately, a better understanding of human VRC pharmacokinetics might contribute to the development of personalized dosing strategies.

Evaluation of a Meropenem and Piperacillin Monitoring Program in Intensive Care Unit Patients Calls for the Regular Assessment of Empirical Targets and Easy-to-Use Dosing Decision Tools.

The drug concentrations targeted in meropenem and piperacillin/tazobactam therapy also depend on the susceptibility of the pathogen. Yet, the pathogen is often unknown, and antibiotic therapy is guided by empirical targets. To reliably achieve the targeted concentrations, dosing needs to be adjusted for renal function. We aimed to evaluate a meropenem and piperacillin/tazobactam monitoring program in intensive care unit (ICU) patients by assessing (i) the adequacy of locally selected empirical targets, (ii) if dosing is adequately adjusted for renal function and individual target, and (iii) if dosing is adjusted in target attainment (TA) failure. In a prospective, observational clinical trial of drug concentrations, relevant patient characteristics and microbiological data (pathogen, minimum inhibitory concentration (MIC)) for patients receiving meropenem or piperacillin/tazobactam treatment were collected. If the MIC value was available, a target range of 1-5 × MIC was selected for minimum drug concentrations of both drugs. If the MIC value was not available, 8-40 mg/L and 16-80 mg/L were selected as empirical target ranges for meropenem and piperacillin, respectively. A total of 356 meropenem and 216 piperacillin samples were collected from 108 and 96 ICU patients, respectively. The vast majority of observed MIC values was lower than the empirical target (meropenem: 90.0%, piperacillin: 93.9%), suggesting empirical target value reductions. TA was found to be low (meropenem: 35.7%, piperacillin 50.5%) with the lowest TA for severely impaired renal function (meropenem: 13.9%, piperacillin: 29.2%), and observed drug concentrations did not significantly differ between patients with different targets, indicating dosing was not adequately adjusted for renal function or target. Dosing adjustments were rare for both drugs (meropenem: 6.13%, piperacillin: 4.78%) and for meropenem irrespective of TA, revealing that concentration monitoring alone was insufficient to guide dosing adjustment. Empirical targets should regularly be assessed and adjusted based on local susceptibility data. To improve TA, scientific knowledge should be translated into easy-to-use dosing strategies guiding antibiotic dosing.

Exploring Dried Blood Spot Cortisol Concentrations as an Alternative for Monitoring Pediatric Adrenal Insufficiency Patients: A Model-Based Analysis.

Congenital adrenal hyperplasia (CAH) is the most common form of adrenal insufficiency in childhood; it requires cortisol replacement therapy with hydrocortisone (HC, synthetic cortisol) from birth and therapy monitoring for successful treatment. In children, the less invasive dried blood spot (DBS) sampling with whole blood including red blood cells (RBCs) provides an advantageous alternative to plasma sampling. Potential differences in binding/association processes between plasma and DBS however need to be considered to correctly interpret DBS measurements for therapy monitoring. While capillary DBS samples would be used in clinical practice, venous cortisol DBS samples from children with adrenal insufficiency were analyzed due to data availability and to directly compare and thus understand potential differences between venous DBS and plasma. A previously published HC plasma pharmacokinetic (PK) model was extended by leveraging these DBS concentrations. In addition to previously characterized binding of cortisol to albumin (linear process) and corticosteroid-binding globulin (CBG; saturable process), DBS data enabled the characterization of a linear cortisol association with RBCs, and thereby providing a quantitative link between DBS and plasma cortisol concentrations. The ratio between the observed cortisol plasma and DBS concentrations varies highly from 2 to 8. Deterministic simulations of the different cortisol binding/association fractions demonstrated that with higher blood cortisol concentrations, saturation of cortisol binding to CBG was observed, leading to an increase in all other cortisol binding fractions. In conclusion, a mathematical PK model was developed which links DBS measurements to plasma exposure and thus allows for quantitative interpretation of measurements of DBS samples.

Towards the Elucidation of the Pharmacokinetics of Voriconazole: A Quantitative Characterization of Its Metabolism.

The small-molecule drug voriconazole (VRC) shows a complex and not yet fully understood metabolism. Consequently, its in vivo pharmacokinetics are challenging to predict, leading to therapy failures or adverse events. Thus, a quantitative in vitro characterization of the metabolism and inhibition properties of VRC for human CYP enzymes was aimed for. The Michaelis-Menten kinetics of voriconazole N-oxide (NO) formation, the major circulating metabolite, by CYP2C19, CYP2C9 and CYP3A4, was determined in incubations of human recombinant CYP enzymes and liver and intestine microsomes. The contribution of the individual enzymes to NO formation was 63.1% CYP2C19, 13.4% CYP2C9 and 29.5% CYP3A4 as determined by specific CYP inhibition in microsomes and intersystem extrapolation factors. The type of inhibition and inhibitory potential of VRC, NO and hydroxyvoriconazole (OH-VRC), emerging to be formed independently of CYP enzymes, were evaluated by their effects on CYP marker reactions. Time-independent inhibition by VRC, NO and OH-VRC was observed on all three enzymes with NO being the weakest and VRC and OH-VRC being comparably strong inhibitors of CYP2C9 and CYP3A4. CYP2C19 was significantly inhibited by VRC only. Overall, the quantitative in vitro evaluations of the metabolism contributed to the elucidation of the pharmacokinetics of VRC and provided a basis for physiologically-based pharmacokinetic modeling and thus VRC treatment optimization.

Plasma concentration guided dosing of drugs used for the treatment of childhood leukaemias: protocol for a systematic review.

INTRODUCTION: Childhood leukaemia is the most common type of cancer in children and represents among 25% of the diagnoses in children <15 years old. Childhood survival rates have significantly improved within the last 40 years due to a rapid advancement in therapeutic interventions. However, in high-risk groups, survival rates remain poor. Pharmacokinetic (PK) data of cancer medications in children are limited and thus current dosing regimens are based on studies with small sample sizes. In adults, large variability in PK is observed and dose individualisation (plasma concentration guided dosing) has been associated with improved clinical outcomes; whether this is true for children is still unknown. This provides an opportunity to explore this strategy in children to potentially reduce toxicities and ensure optimal dosing. This paper will provide a protocol to systematically review studies that have used dose individualisation of drugs used in the treatment of childhood leukaemias. METHODS AND ANALYSIS: Systematic review methodology will be applied to identify, select and extract data from published plasma guided dosing studies conducted in a paediatric leukaemia cohort. Databases (eg, Ovid Embase, Ovid MEDLINE, Ovid Cochrane) and clinical trial registries (CENTRAL, ClinicalTrials.gov and ISRCTN) will be used to perform the systematic literature search (up until February 2021). Only full empirical studies will be included, with primary clinical outcomes (progression-free survival, toxicities, minimal residual disease status, complete cytogenetic response, partial cytogenetic response and major molecular response) being used to decide whether the study will be included. The quality of included studies will be undertaken, with a subgroup analysis where appropriate. ETHICS AND DISSEMINATION: This systematic review will not require ethics approval as there will not be collection of primary data. Findings of this review will be made available through publications in peer-reviewed journals and conference presentations. Gaps will be identified in current literature to inform future-related research. PROSPERO REGISTRATION NUMBER: CRD42021225045.

A versatile high-performance LC-MS/MS assay for the quantification of voriconazole and its N-oxide metabolite in small sample volumes of multiple human matrices for biomedical applications.

Voriconazole (VRC) pharmacokinetics, in particular its complex metabolism, is still not fully understood which challenges its optimal therapeutic use. To increase knowledge on the pharmacokinetics of this antifungal drug, it is essential to broaden the perspective and expand in vitro and clinical in vivo investigations in particular to aspects such as unbound plasma, target-site and metabolite concentrations. Innovative sampling approaches such as microdialysis, a minimally-invasive technique for the analysis of compound concentrations in target-site human tissue fluids, are associated with bioanalytical challenges, i.e. small sample volumes and low concentrations. Thus, a bioanalytical LC-MS/MS assay for the simultaneous quantification of VRC and its main N-oxide (NO) metabolite in human plasma, ultrafiltrate and microdialysate was developed and validated according to the European Medicines Agency guideline. Quantification was rapid, simple and feasible for clinically relevant concentrations from 5 to 5000 ng/mL in plasma and ultrafiltrate as well as from 4 to 4000 ng/mL in microdialysate. Due to the high sensitivity of the assay, only 20 µL of plasma or ultrafiltrate and 5 µL of microdialysate were required. For VRC and NO in all matrices, between-run accuracy was high with a maximum mean deviation of 7.0% from the nominal value and between-run precision was demonstrated by ≤ 11.8% coefficient of variation. Both compounds proved stable under various conditions. The assay suitability was demonstrated by the application to a clinical study quantifying simultaneously VRC and NO concentrations in plasma, ultrafiltrate and microdialysate. Additionally, the assay was successfully adapted for pharmacokinetic analyses in human tissue-derived in vitro experiments. Overall, by reducing the required sample volume, the bioanalytical method allows for an increased number of plasma samples in vulnerable populations, e.g. infants, and enables the generation of concentration-time profiles with a higher temporal resolution in microdialysis studies. Consequently, the developed assay is apt to elucidate the complex pharmacokinetics of VRC in clinical settings as prerequisite for therapy optimisation.

Insights in the maturational processes influencing hydrocortisone pharmacokinetics in congenital adrenal hyperplasia patients using a middle-out approach.

Introduction: Hydrocortisone is the standard of care in cortisol replacement therapy for congenital adrenal hyperplasia patients. Challenges in mimicking cortisol circadian rhythm and dosing individualization can be overcome by the support of mathematical modelling. Previously, a non-linear mixed-effects (NLME) model was developed based on clinical hydrocortisone pharmacokinetic (PK) pediatric and adult data. Additionally, a physiologically-based pharmacokinetic (PBPK) model was developed for adults and a pediatric model was obtained using maturation functions for relevant processes. In this work, a middle-out approach was applied. The aim was to investigate whether PBPK-derived maturation functions could provide a better description of hydrocortisone PK inter-individual variability when implemented in the NLME framework, with the goal of providing better individual predictions towards precision dosing at the patient level. Methods: Hydrocortisone PK data from 24 adrenal insufficiency pediatric patients and 30 adult healthy volunteers were used for NLME model development, while the PBPK model and maturation functions of clearance and cortisol binding globulin (CBG) were developed based on previous studies published in the literature. Results: Clearance (CL) estimates from both approaches were similar for children older than 1 year (CL/F increasing from around 150 L/h to 500 L/h), while CBG concentrations differed across the whole age range (CBGNLME stable around 0.5 µM vs. steady increase from 0.35 to 0.8 µM for CBG PBPK). PBPK-derived maturation functions were subsequently included in the NLME model. After inclusion of the maturation functions, none, a part of, or all parameters were re-estimated. However, the inclusion of CL and/or CBG maturation functions in the NLME model did not result in improved model performance for the CL maturation function (ΔOFV > -15.36) and the re-estimation of parameters using the CBG maturation function most often led to unstable models or individual CL prediction bias. Discussion: Three explanations for the observed discrepancies could be postulated, i) non-considered maturation of processes such as absorption or first-pass effect, ii) lack of patients between 1 and 12 months, iii) lack of correction of PBPK CL maturation functions derived from urinary concentration ratio data for the renal function relative to adults. These should be investigated in the future to determine how NLME and PBPK methods can work towards deriving insights into pediatric hydrocortisone PK.