One-compartment Model Research Articles

#2422 Iohexol plasma clearance measurement protocol standardization—a consensus of the European Kidney Function Consortium

Abstract Background and Aims International consensus supports the development of standardized protocols for measured glomerular filtration rate (mGFR) to facilitate the adoption of mGFR determination in both clinical and research settings. Method The European Kidney Function Consortium (EKFC) convened an international group of mGFR experts. The working group performed an extensive literature search to inform the development of recommendations for mGFR determination using one-compartment plasma clearance models and iohexol as the exogenous filtration marker. Iohexol was selected as it is non-radio labeled, inexpensive, widely available, and can be assayed at a central laboratory. Other commonly used non-radio labeled tracers have been (inulin) or are soon to be (iothalamate) discontinued. A plasma clearance model was selected as it requires no timed urine collections. A one-compartment model was preferred to a two compartment model as it requires fewer samples. The resulting recommendations were based on published evidence complemented by expert opinion. Results Recommendations include practical advice for patients and health professionals, preparation, administration and safety aspects of iohexol, laboratory analysis, blood sample collection and sample timing for both multiple and single sample protocols, description of the mGFR mathematical calculations as well as implementation strategies (Fig.). EKFC will provide materials such as patient information sheets, standard operating procedures, a study protocol template, and support for mGFR calculation. Conclusion Members of the EKFC have evaluated the evidence and formulated recommendations that include all aspects of iohexol plasma clearance measurement using a one compartment model. *Author list on behalf of the European Kidney Function Consortium (EKFC).

Identifiability of spatiotemporal tissue perfusion models

Objective. Standard models for perfusion quantification in DCE-MRI produce a bias by treating voxels as isolated systems. Spatiotemporal models can remove this bias, but it is unknown whether they are fundamentally identifiable. The aim of this study is to investigate this question in silico using one-dimensional toy systems with a one-compartment blood flow model and a two-compartment perfusion model. Approach. For each of the two models, identifiability is explored theoretically and in-silico for three systems. Concentrations over space and time are simulated by forward propagation. Different levels of noise and temporal undersampling are added to investigate sensitivity to measurement error. Model parameters are fitted using a standard gradient descent algorithm, applied iteratively with a stepwise increasing time window. Model fitting is repeated with different initial values to probe uniqueness of the solution. Reconstruction accuracy is quantified for each parameter by comparison to the ground truth. Main results. Theoretical analysis shows that flows and volume fractions are only identifiable up to a constant, and that this degeneracy can be removed by proper choice of parameters. Simulations show that in all cases, the tissue concentrations can be reconstructed accurately. The one-compartment model shows accurate reconstruction of blood velocities and arterial input functions, independent of the initial values and robust to measurement error. The two-compartmental perfusion model was not fully identifiable, showing good reconstruction of arterial velocities and input functions, but multiple valid solutions for the perfusion parameters and venous velocities, and a strong sensitivity to measurement error in these parameters. Significance. These results support the use of one-compartment spatiotemporal flow models, but two-compartment perfusion models were not sufficiently identifiable. Future studies should investigate whether this degeneracy is resolved in more realistic 2D and 3D systems, by adding physically justified constraints, or by optimizing experimental parameters such as injection duration or temporal resolution.

Population pharmacokinetics and dose rationale for aciclovir in term and pre-term neonates with herpes.

Aciclovir is considered the first-line treatment against Herpes simplex virus (HSV) infections in new-borns and infants. As renal excretion is the major route of elimination, in renally-impaired patients, aciclovir doses are adjusted according to the degree of impairment. However, limited attention has been given to the implications of immature renal function or dysfunction due to the viral disease itself. The aim of this investigation was to characterize the pharmacokinetics of aciclovir taking into account maturation and disease processes in the neonatal population. Pharmacokinetic data obtained from 2 previously published clinical trials (n = 28) were analyzed using a nonlinear mixed effects modeling approach. Post-menstrual age (PMA) and creatinine clearance (CLCR) were assessed as descriptors of maturation and renal function. Simulation scenarios were also implemented to illustrate the use of pharmacokinetic data to extrapolate efficacy from adults. Aciclovir pharmacokinetics was described by a one-compartment model with first-order elimination. Body weight and diagnosis (systemic infection) were statistically significant covariates on the volume of distribution, whereas body weight, CLCR and PMA had a significant effect on clearance. Median clearance varied from 0.2 to 1.0 L/h in subjects with PMA <34 or ≥34 weeks, respectively. Population estimate for volume of distribution was 1.93 L with systemic infection increasing this value by almost 3-fold (2.67 times higher). A suitable model parameterization was identified, which discriminates the effects of developmental growth, maturation, and organ function. Exposure to aciclovir was found to increase with decreasing PMA and renal function (CLCR), suggesting different dosing requirement for pre-term neonates.

Development of a Web Application for Simulating Plasma Drug Concentrations in Patients with Zolpidem Intoxication.

Zolpidem is a widely prescribed hypnotic Z-drug used to treat short-term insomnia. However, a growing number of individuals intentionally overdose on these drugs. This study aimed to develop a predictive tool for physicians to assess patients with zolpidem overdose. A population pharmacokinetic (PK) model was established using digitized data obtained from twenty-three healthy volunteers after a single oral administration of zolpidem. Based on the final PK model, a web application was developed using open-source R packages such as rxode2, nonmem2rx, and shiny. The final model was a one-compartment model with first-order absorption and elimination with PK parameters, including clearance (CL, 16.9 L/h), absorption rate constant (Ka, 5.41 h-1), volume of distribution (Vd, 61.7 L), and lag time (ALAG, 0.394 h). Using the established population PK model in the current study, we developed a web application that enables users to simulate plasma zolpidem concentrations and visualize their profiles. This user-friendly web application may provide essential clinical information to physicians, ultimately helping in the management of patients with zolpidem intoxication.

Integrating Clopidogrel's First-Pass Effect in a Joint Semi-Physiological Population Pharmacokinetic Model of the Drug and Its Inactive Carboxylic Acid Metabolite.

Clopidogrel (CLO), a pro-drug for preventing thrombotic events, undergoes rapid absorption and extensive metabolism, with approximately 85-90% converted to an inactive carboxylic acid metabolite (CLO-CA) and the remaining to an active thiol (CLO-TH). Few pharmacokinetic models for the drug and its metabolites exist, with most focusing on CLO-TH. Although CLO-CA is inactive, its predominant (compared to its parent drug and metabolites) presence in plasma underscores the importance of characterizing its formation and pharmacokinetic profile. This study aimed to characterize the process of the absorption of CLO and its conversion to CLO-CA via developing a population pharmacokinetic model. Individual participants' data from two bioequivalence studies were utilized. Extensive blood samples were collected at predetermined intervals, including 841 concentrations of CLO and 1149 of CLO-CA. A nonlinear, mixed-effects modelling approach using NONMEM® software (v 7.5) was applied. A one-compartment model was chosen for CLO, while a two-compartment proved optimal for CLO-CA. Absorption from the depot compartment was modeled via two transit compartments, incorporating transit rate constants (Ktr). A semi-physiological model explained the first-pass effect of CLO, integrating a liver compartment. The estimated mean transit times (MTTs) for the studies were 0.470 and 0.410 h, respectively. The relative bioavailability for each study's generic medicine compared to the reference were 1.08 and 0.960, respectively. Based on the estimated parameters, the fractions metabolized to inactive metabolites (FiaM_st1 and FiaM_st2) were determined to be 87.27% and 86.87% for the two studies, respectively. The appropriateness of the final model was confirmed. Our model offers a robust framework for elucidating the pharmacokinetic profiles of CLO and CLO-CA.

Modelling the In Vivo and Ex Vivo DNA Damage Response after Internal Irradiation of Blood from Patients with Thyroid Cancer.

This work reports on a model that describes patient-specific absorbed dose-dependent DNA damage response in peripheral blood mononuclear cells of thyroid cancer patients during radioiodine therapy and compares the results with the ex vivo DNA damage response in these patients. Blood samples of 18 patients (nine time points up to 168 h post-administration) were analyzed for radiation-induced γ-H2AX + 53BP1 DNA double-strand break foci (RIF). A linear one-compartment model described the absorbed dose-dependent time course of RIF (Parameters: c characterizes DSB damage induction; k1 and k2 are rate constants describing fast and slow repair). The rate constants were compared to ex vivo repair rates. A total of 14 patient datasets could be analyzed; c ranged from 0.012 to 0.109 mGy-1, k2 from 0 to 0.04 h-1. On average, 96% of the damage is repaired quickly with k1 (range: 0.19-3.03 h-1). Two patient subgroups were distinguished by k1-values (n = 6, k1 > 1.1 h-1; n = 8, k1 < 0.6 h-1). A weak correlation with patient age was observed. While induction of RIF was similar among ex vivo and in vivo, the respective repair rates failed to correlate. The lack of correlation between in vivo and ex vivo repair rates and the applicability of the model to other therapies will be addressed in further studies.

Assessment of body mass-related covariates for rifampicin pharmacokinetics in healthy Caucasian volunteers

PurposeCurrently, body weight-based dosing of rifampicin is recommended. But lately, fat-free mass (FFM) was reported to be superior to body weight (BW). The present evaluation aimed to assess the influence of body mass-related covariates on rifampicin’s pharmacokinetics (PK) parameters in more detail using non-linear mixed effects modeling (NLMEM).MethodsTwenty-four healthy Caucasian volunteers were enrolled in a bioequivalence study, each receiving a test and a reference tablet of 600 mg of rifampicin separated by a wash-out period of at least 9 days. Monolix version 2023R1 was used for NLMEM. Monte Carlo simulations (MCS) were performed to visualize the relationship of body size descriptors to the exposure to rifampicin.ResultsA one-compartment model with nonlinear (Michaelis–Menten) elimination and zero-order absorption kinetics with a lag time best described the data. The covariate model including fat-free mass (FFM) on volume of distribution (V/F) and on maximum elimination rate (Vmax/F) lowered the objective function value (OFV) by 56.4. The second-best covariate model of sex on V/F and Vmax/F and BW on V/F reduced the OFV by 51.2. The decrease in unexplained inter-individual variability on Vmax/F in both covariate models was similar. For a given dose, MCS showed lower exposure to rifampicin with higher FFM and accordingly in males compared to females with the same BW and body height.ConclusionOur results indicate that beyond BW, body composition as reflected by FFM could also be relevant for optimized dosing of rifampicin. This assumption needs to be studied further in patients treated with rifampicin.

Population Pharmacokinetics of Sertraline in Psychiatric and Substance Use Disorders.

This study aimed to characterize the population pharmacokinetics of sertraline in Mexican patients with psychiatric and substance use disorders. Fifty-nine patients (13 to 76 years old) treated with doses of sertraline between 12.5 and 100 mg/day were included. Plasma sertraline concentrations were determined in blood samples and five of the main substances of abuse were determined by rapid tests in urine samples. Demographic, clinical, and pharmacogenetic factors were also evaluated. Population pharmacokinetic analysis was performed using NONMEM software with first-order conditional estimation method. A one-compartment model with proportional residual error adequately described the sertraline concentrations versus time. CYP2D6*2 polymorphism and CYP2C19 phenotypes significantly influenced sertraline clearance, which had a population mean value of 66 L/h in the final model. The absorption constant and volume of distribution were fixed at 0.855 1/h and 20.2 L/kg, respectively. The model explained 11.3% of the interindividual variability in sertraline clearance. The presence of the CYP2D6*2 polymorphism caused a 23.1% decrease in sertraline clearance, whereas patients with intermediate and poor phenotype of CYP2C19 showed 19.06% and 48.26% decreases in sertraline clearance, respectively. The model was internally validated by bootstrap and visual predictive check. Finally, stochastic simulations were performed to propose dosing regimens to achieve therapeutic levels that contribute to improving treatment response.

A switching lung mechanics model for detection of expiratory flow limitation

Abstract Expiratory flow limitation (EFL) is an often unrecognized clinical condition with a multitude of negative implications. A mathematical EFL model is proposed to detect flow limitations automatically. The EFL model is a switching one-compartment lung mechanics model with a volume-dependent airway resistance to simulate the dynamic behavior during expiration. The EFL detection is based on a breath-by-breath model parameter identification and validated on clinical data of mechanically ventilated patients. In the severe flow limitation group 93.9 % ± 5 % and in the no limitation group 10.2 % ± 13.7 % of the breaths are detected as EFL. Based on the high detection rate of EFL, these results support the usefulness of the EFL detection. It is a first step toward an automated detection of EFL in clinical applications and may help to reduce underdiagnosis of EFL.

Towards process-based modelling and parameterisation of bioaccumulation in humans across PFAS congeners

Per- and polyfluoroalkyl substances (PFAS) are a large class of stable toxic chemicals which have ended up in the environment and in organisms in significant concentrations. Toxicokinetic models are needed to facilitate extrapolation of bioaccumulation data across PFAS congeners and species. For the present study, we carried out an inventory of accumulation processes specific for PFAS, deviating from traditional Persistent Organic Pollutants (POPs). In addition, we reviewed toxicokinetic models on PFAS reported in literature, classifying them according to the number of compartments distinguished as a one-compartment model (1-CM), two-compartment model (2- CM) or a multi-compartment model, (multi-CM) as well as the accumulation processes included and the parameters used. As the inventory showed that simple 1-CMs were lacking, we developed a generic 1-CM of ourselves to include PFAS specific processes and validated the model for legacy perfluoroalkyl acids. Predicted summed elimination constants were accurate for long carbon chains (>C6), indicating that the model properly represented toxicokinetic processes for most congeners. Results for urinary elimination rate constants were mixed, which might be caused by the exclusion of reabsorption processes (renal reabsorption, enterohepatic circulation). The 1-CM needs to be improved further in order to better predict individual elimination pathways. Besides that, more data on PFAS-transporter specific processes are needed to extrapolate across PFAS congeners and species.

Pharmacokinetics of oral clonazepam in growing commercial pigs (Sus scrofa domestica).

Clonazepam causes sedation and psychomotor impairment in people. Due to similarities between people and swine in response to benzodiazepines, clonazepam may represent a viable option to produce mild-to-moderate tranquillization in pigs. The objective of this study was to determine the pharmacokinetic profile of a single oral dose (0.5 mg/kg) of clonazepam in eight healthy, growing commercial cross pigs. Serial plasma samples were collected at baseline and up to 96 h after administration. Plasma concentrations were quantified using reverse-phase high-performance liquid chromatography, and compartment models were fit to time-concentration data. A one-compartment first-order model best fits the data. Maximum plasma concentration was 99.5 ng/mL, and time to maximum concentration was 3.4 h. Elimination half-life was 7.3 h, mean residence time 7.4 h, and apparent volume of distribution 5.7 L/kg. Achieved plasma concentrations exceeded those associated with psychomotor impairment in people although pharmacodynamic effects have not been investigated in pigs. A simulated oral regimen consisting of 0.35 mg/kg administered every 8 h to pigs would achieve plasma concentrations above 32 ng/mL which are shown to produce psychomotor impairment in people. Further studies to test the clinical efficacy of these dosages in commercial and miniature pigs are warranted.

Population pharmacokinetics of nirmatrelvir in Chinese patients with COVID-19: Therapeutic drug monitoring and dosing regimen selection in clinical practice

ObjectivesTo establish a population pharmacokinetics (PopPK) model of nirmatrelvir in Chinese COVID-19 patients and provide reference for refining the dosing strategy of nirmatrelvir in patients confirmed to be infected with SARS-CoV-2. MethodsA total of 80 blood samples were obtained from 35 mild to moderate COVID-19 patients who were orally administered nirmatrelvir/ritonavir tablets. The PopPK model of nirmatrelvir was developed using a nonlinear mixed effects modelling approach. The stability and prediction of the final model were assessed through a combination of goodness-of-fit and bootstrap method. The exposure of nirmatrelvir across various clinical scenarios was simulated using Monte Carlo simulations. ResultsThe pharmacokinetics of nirmatrelvir was well characterised by a one-compartment model with first-order absorption, and with creatinine clearance (Ccr) as the significant covariate. Typical population parameter estimates of apparent clearance and distribution volume for a patient with a Ccr of 95.5 mL·min−1were 3.45 L·h−1 and 48.71 L, respectively. The bootstrap and visual predictive check procedures demonstrated satisfactory predictive performance and robustness of the final model. ConclusionThe final model was capable of offering an early prediction of drug concentration ranges for different nirmatrelvir dosing regimens and optimise the dose regimen of nirmatrelvir in individuals with confirmed SARS-CoV-2 infection.

Per- and polyfluoroalkyl substances (PFAS) and fetal growth: A nation-wide register-based study on PFAS in drinking water

BackgroundThere is inconclusive evidence for an association between per- and polyfluoroalkyl substances (PFAS) and fetal growth. ObjectivesWe conducted a nation-wide register-based cohort study to assess the associations of the estimated maternal exposure to the sum (PFAS4) of perfluorooctane sulfonic acid (PFOS), perfluorooctanoic acid (PFOA), perfluorononanoic acid (PFNA) and perfluorohexane sulfonic acid (PFHxS) with birthweight as well as risk of small- (SGA) and large-for-gestational-age (LGA). Materials and methodsWe included all births in Sweden during 2012–2018 of mothers residing ≥ four years prior to partus in localities served by municipal drinking water where PFAS were measured in raw and drinking water. Using a one-compartment toxicokinetic model we estimated cumulative maternal blood levels of PFAS4 during pregnancy by linking residential history, municipal PFAS water concentration and year-specific background serum PFAS concentrations in Sweden. Individual birth outcomes and covariates were obtained via register linkage. Mean values and 95 % confidence intervals (CI) of β coefficients and odds ratios (OR) were estimated by linear and logistic regressions, respectively. Quantile g-computation regression was conducted to assess the impact of PFAS4 mixture. ResultsAmong the 248,804 singleton newborns included, no overall association was observed for PFAS4 and birthweight or SGA. However, an association was seen for LGA, multivariable-adjusted OR 1.08 (95% CI: 1.01–1.16) when comparing the highest PFAS4 quartile to the lowest. These associations remained for mixture effect approach where all PFAS, except for PFOA, contributed with a positive weight. DiscussionsWe observed an association of the sum of PFAS4 – especially PFOS – with increased risk of LGA, but not with SGA or birthweight. The limitations linked to the exposure assessment still require caution in the interpretation.

Assessing Pharmacogenomic loci Associated with the Pharmacokinetics of Vamorolone in Boys with Duchenne Muscular Dystrophy.

Human genetic variation (polymorphisms) in genes coding proteins involved in the absorption, distribution, metabolism, and elimination (ADME) of drugs can have a strong effect on drug exposure and downstream efficacy and safety outcomes. Vamorolone, a dissociative steroidal anti-inflammatory drug for treating Duchenne muscular dystrophy (DMD), primarily undergoes oxidation by CYP3A4 and CYP3A5 and glucuronidation by UDP-glucuronosyltransferases. This work assesses the pharmacokinetics (PKs) of vamorolone and sources of interindividual variability (IIV) in 81 steroid-naïve boys with DMD aged 4 to <7 years old considering the genetic polymorphisms of CYPS3A4 (CYP3A4*22, CYP3A4*1B), CYP3A5 (CYP3A5*3), and UGT1A1 (UGT1A1*60) utilizing population PK modeling. A one-compartment model with zero-order absorption (Tk0, duration of absorption), linear clearance (CL/F), and volume (V/F) describes the plasma PK data for boys with DMD receiving a wide range of vamorolone doses (0.25-6 mg/kg/day). The typical CL/F and V/F values of vamorolone were 35.8 L/h and 119 L, with modest IIV. The population Tk0 was 3.14 h yielding an average zero-order absorption rate (k0) of 1.16 mg/kg/h with similar absorption kinetics across subjects at the same vamorolone dose (i.e., no IIV on Tk0). The covariate analysis showed that none of the genetic covariates had any significant impact on the PKs of vamorolone in boys with DMD. Thus, the PKs of vamorolone is very consistent in these young boys with DMD.

Genetic and clinical predictors of rifapentine and isoniazid pharmacokinetics in paediatrics with tuberculosis infection.

Twelve weekly doses of rifapentine and isoniazid (3HP regimen) are recommended for TB preventive therapy in children with TB infection. However, they present with variability in the pharmacokinetic profiles. The current study aimed to develop a pharmacokinetic model of rifapentine and isoniazid in 12 children with TB infection using NONMEM. Ninety plasma and 41 urine samples were collected at Week 4 of treatment. Drug concentrations were measured using a validated HPLC-UV method. MassARRAY® SNP genotyping was used to investigate genetic factors, including P-glycoprotein (ABCB1), solute carrier organic anion transporter B1 (SLCO1B1), arylacetamide deacetylase (AADAC) and N-acetyl transferase (NAT2). Clinically relevant covariates were also analysed. A two-compartment model for isoniazid and a one-compartment model for rifapentine with transit compartment absorption and first-order elimination were the best models for describing plasma and urine data. The estimated (relative standard error, RSE) of isoniazid non-renal clearance was 3.52 L·h-1 (23.1%), 2.91 L·h-1 (19.6%), and 2.58 L·h-1 (20.0%) in NAT2 rapid, intermediate and slow acetylators. A significant proportion of the unchanged isoniazid was cleared renally (2.7 L·h-1; 8.0%), while the unchanged rifapentine was cleared primarily through non-renal routes (0.681 L·h-1; 3.6%). Participants with the ABCB1 mutant allele had lower bioavailability of rifapentine, while food prolonged the mean transit time of isoniazid. ABCB1 mutant allele carriers may require higher rifapentine doses; however, this must be confirmed in larger trials. Food did not affect overall exposure to isoniazid and only delayed absorption time.

Development of a Weight-Band Dosing Approach for Vosoritide in Children with Achondroplasia Using a Population Pharmacokinetic Model.

Vosoritide is a recently approved therapy for achondroplasia, the most common form of disproportionate short stature, that has been shown to be well tolerated and effective in increasing linear growth. This study aimed to develop a population pharmacokinetic (PPK) model to characterize pharmacokinetics (PK) of vosoritide and establish a weight-band dosing regimen. A PPK model was developed using data from five clinical trials in children with achondroplasia (aged 0.95-15 years) who received daily per-kg doses of vosoritide. The model was used to simulate expected exposures in children with a refined weight-band dosing regimen. Simulated exposure was compared with the observed exposure from the pivotal clinical trial to evaluate appropriateness of the weight-band dosing regimen. A one-compartment model with a change-point first-order absorption and first-order elimination accurately described PK of vosoritide in children with achondroplasia. Body weight was found to be a predictor of vosoritide's clearance and volume of distribution. Additionally, it was observed that dosing solution concentration and duration of treatment influenced bioavailability. The weight-band dosing regimen resulted in simulated exposures that were within the range demonstrated to be well tolerated and effective in the pivotal clinical trial and showed improved consistency in drug exposure across the achondroplasia population. The weight-band dosing regimen reduced the number of recommended dose levels by body weight and is expected to simplify dosing for children with achondroplasia and their caregivers. NCT02055157, NCT02724228, NCT03197766, NCT03424018, and NCT03583697.

Rapid identification of potential nonsteroidal anti-inflammatory drug overdose-induced liver toxicity and prediction of follow-up exposure: Integrating bioanalytical and population pharmacokinetic assay.

Nonsteroidal anti-inflammatory drugs (NSAIDs) are among the most frequently used drugs that can cause liver toxicity. The aim of this study was to integrate bioanalytical and population pharmacokinetic (PopPK) assay to rapidly screen and quantify the concentrations of NSAIDs in plasma and monitor clinical safety. A liquid chromatography-tandem mass spectrometry (LC-MS/MS) method was developed for the simultaneous quantification of acetaminophen (APAP), flurbiprofen (FLB), aspirin (ASP), and ibuprofen (IBP), four commonly used NSAIDs. The PopPK model of the signature toxicant was analyzed based on the published literature. The LC-MS/MS method was successfully validated and applied to determine NSAID concentrations in patient plasma samples. APAP, ASP, and IBP data were best fitted using a one-compartment model, and FLB data were best fitted using a two-compartment model. Bootstrapping and visual predictive checks suggested that a robust and reliable pharmacokinetic model was developed. A fast, simple, and sensitive LC-MS/MS method was developed and validated for determining APAP, FLB, ASP, and IBP in human plasma. Combined with the PopPK model, this method was applied to rapidly analyze the concentrations of NSAIDs in clinical samples from patients presenting to the emergency department with acute liver dysfunction and monitored NSAIDs clinical safety.

Time-Varying Clearance in Milrinone Pharmacokinetics from Premature Neonates to Adolescents.

Milrinone is an inotrope and vasodilator used for prophylaxis or treatment of low cardiac output syndrome after weaning from cardiopulmonary bypass (CPB). It is renally eliminated and has an acceptable therapeutic range of 100-300μg/L, but weight-based dosing alone is associated with poor target attainment. We aimed to develop a population pharmacokinetic model for milrinone from premature neonates to adolescents, and to evaluate how age, renal function and recovery from CPB may impact dose selection. Fifty paediatric patients (aged 4 days to 16 years) were studied after undergoing cardiac surgery supported by CPB. Data from 29 premature neonates (23-28 weeks' postmenstrual age) treated for prophylaxis of low systemic blood flow were available for a pooled pharmacokinetic analysis. Population parameters were estimated using non-linear mixed effects modelling (NONMEM 7.5.1). There were 369 milrinone measurements available for analysis. A one-compartment model with zero-order input and first-order elimination was used to describe milrinone disposition. Population parameters were clearance 17.8L/70kg [95% CI 15.8-19.9] and volume 20.4L/h/70kg [95% CI 17.8-22.1]. Covariates included size, postmenstrual age and renal function for clearance, and size and postnatal age for volume. Milrinone clearance is reduced by 39.5% [95% CI 24.0-53.7] immediately after bypass, and recovers to baseline clearance with a half-time of 12.0h [95% CI 9.7-15.2]. Milrinone volume was 2.07 [95% CI 1.87-2.27] times greater at birth than the population standard and decreased over the first days of life with a half-time of 0.977days [95% CI 0.833-1.12]. Milrinone is predominately renally eliminated and so renal function is an important covariate describing variability in clearance. Increasing clearance over time likely reflects increasing cardiac output and renal perfusion due to milrinone and return to baseline following CPB.

Impact of covariate model building methods on their clinical relevance evaluation in population pharmacokinetic analyses: comparison of the full model, stepwise covariate model (SCM) and SCM+ approaches.

Covariate analysis in population pharmacokinetics is key for adjusting doses for patients. The main objective of this work was to compare the adequacy of various modeling approaches on covariate clinical relevance decision-making. The full model, stepwise covariate model (SCM) and SCM+ PsN algorithms were compared in a clinical trial simulation of a 383-patient population pharmacokinetic study mixing rich and sparse designs. A one-compartment model with first-order absorption was used. A base model including a body weight effect on CL/F and V/F and a covariate model including 4 additional covariates-parameters relationships were simulated. As for forest plots, ratios between covariates at a specific value and that of a typical individual were calculated with their 90% confidence interval (CI90) using standard errors. Covariates on CL, V and KA were considered relevant if their CI90 fell completely outside the reference area [0.8-1.2]. All approaches provided unbiased covariate ratio estimates. For covariates with a simulated effect, the 3 approaches correctly identify their clinical relevance. However, significant covariates were missed in up to 15% of cases with SCM/SCM+. For covariate with no simulated effects, the full model mainly identified them as non-relevant or with insufficient information while SCM/SCM+ mainly did not select them. SCM/SCM+ assume that non-selected covariates are non-relevant when it could be due to insufficient information, whereas the full model does not make this assumption and is faster. This study must be extended to other methods and completed by a more complex high-dimensional simulation framework.

Model-based comparison of subcutaneous versus sublingual apomorphine administration in the treatment of motor fluctuations in Parkinson’s disease

The objective of this study was to compare the effectiveness of subcutaneous (SC) and sublingual (SL) formulations of apomorphine for the treatment of motor fluctuations in Parkinson’s disease using a pharmacokinetics (PK)/pharmacodynamics (PD) modeling approach. The PK of SC and SL apomorphine are best described by a one-compartment model with first-order absorption and a two-compartment model with delayed absorption, respectively. The PK/PD model relating apomorphine plasma concentrations to the Unified Parkinson’s Disease Rating Scale (UPDRS) motor scores was described by a sigmoidal Emax model assuming effective concentration = drug concentration in an effect compartment. Apomorphine concentrations and UPDRS motor scores were simulated from the PK/PD models using 500 hypothetical subjects. UPDRS motor score change from baseline was evaluated using time to clinically relevant response, response duration, area under the curve, maximal response, and time to maximal response. Higher doses of each apomorphine formulation were associated with shorter time to response, longer response duration, and greater maximal response. Although the mean maximal responses to SC and SL apomorphine were comparable, the time to response was four times shorter (7 vs. 31 min) and time to maximal response was two times shorter (27 vs. 61 min) for 4 mg SC vs. 50 mg SL. Thus, faster onset of action was observed for the SC formulation compared to SL. These data may be useful for physicians when selecting “on demand” therapy for patients with Parkinson’s disease experiencing motor fluctuations.