Nonlinear Mixed-effect Modeling Approach Research Articles

Five multivariate Duchenne muscular dystrophy progression models bridging six-minute walk distance and MRI relaxometry of leg muscles.

The study aimed to provide quantitative information on the utilization of MRI transverse relaxation time constant (MRI-T2) of leg muscles in DMD clinical trials by developing multivariate disease progression models of Duchenne muscular dystrophy (DMD) using 6-min walk distance (6MWD) and MRI-T2. Clinical data were collected from the prospective and longitudinal ImagingNMD study. Disease progression models were developed by a nonlinear mixed-effect modeling approach. Univariate models of 6MWD and MRI-T2 of five muscles were developed separately. Age at assessment was the time metric. Multivariate models were developed by estimating the correlation of 6MWD and MRI-T2 modelvariables. Full model estimation approach for covariate analysis and five-fold cross validation were conducted. Simulations were performed to compare the models and predict the covariate effects on the trajectories of 6MWD and MRI-T2. Sigmoid Imax and Emax models best captured the profiles of 6MWD and MRI-T2 over age. Steroid use, baseline 6MWD, and baseline MRI-T2 were significant covariates. The median age at which 6MWD is half of its maximum decrease in the five models was similar, while the median age at which MRI-T2 is half of its maximum increase varied depending on the type of muscle. The models connecting 6MWD and MRI-T2 successfully quantified how individual characteristics alter disease trajectories. The models demonstrate a plausible correlation between 6MWD and MRI-T2, supporting the use of MRI-T2. The developed models will guide drug developers in using the MRI-T2 to most efficient use in DMD clinical trials.

Population pharmacokinetics of a novel PCSK9 antisense oligonucleotide.

The aim of this study was to characterize the population pharmacokinetics of AZD8233, an antisense oligonucleotide (ASO) that targets the PCSK9 transcript to reduce hepatocyte PCSK9 protein production and plasma levels. AZD8233 utilizes generation 2.5S-constrained ethyl motif (cET) chemistry and is conjugated to a triantennary N-acetylgalactosamine (GalNAc3) ligand for targeted hepatocyte uptake. A non-linear mixed-effect modelling approach utilizing NONMEM software was applied to AZD8233 concentration-time data from 3416 samples in 219 participants from four phase 1-2 studies, one in healthy volunteers (NCT03593785) and three in patients with dyslipidaemia (NCT04155645, NCT04641299 and NCT04823611). The final model described the AZD8233 plasma concentration-time profile from four phase 1-2 studies in healthy volunteers or participants with dyslipidaemia, covering a dose range of 4 to 120 mg. The pharmacokinetics of AZD8233 were adequately described by a two-compartment model with first-order absorption. The supra-proportional increase in maximum plasma concentration (Cmax ) across the observed dose range was described by non-linear Michaelis-Menten elimination (maximum elimination rate, 9.9mg/h [12% relative standard error]; concentration yielding half-maximal elimination rate, 4.8mg/L [18% relative standard error]). Body weight, sex, estimated glomerular filtration rate and disease status (healthy participant vs. patient with dyslipidaemia) were identified as factors affecting exposure to AZD8233. Covariate analysis showed body weight to be the main factor affecting exposure to AZD8233, which largely explained the higher Cmax observed in the Asian population relative to non-Asians.

PK Modeling of L-4-Boronophenylalanine and Development of Bayesian Predictive Platform for L-4-Boronophenylalanine PKs for Boron Neutron Capture Therapy.

L-4-[(10B)]Boronophenylalanine (BPA) is an amino acid analogue with a boron-10 moiety. It is most widely used as a boron carrier in boron neutron capture therapy. In this study, a Bayesian predictive platform of blood boron concentration based on a BPA pharmacokinetic (PK) model was developed. This platform is user-friendly and can predict the individual boron PK and optimal time window for boron neutron capture therapy in a simple way. The present study aimed to establish a PK model of L-4-boronophenylalanine and develop a Bayesian predictive platform for blood boron PKs for user-friendly estimation of boron concentration during neutron irradiation of neutron capture therapy. Whole blood boron concentrations from seven previous reports were graphically extracted and analyzed using the nonlinear mixed-effects modeling (NONMEM) approach. Model robustness was assessed using nonparametric bootstrap and visual predictive check approaches. The visual predictive check indicated that the final PK model is able to adequately predict observed concentrations. The Shiny package was used to input real-time blood boron concentration data, and during the following irradiation session blood boron was estimated with an acceptably short calculation time for the determination of irradiation time. Finally, a user-friendly Bayesian estimation platform for BPA PKs was developed to optimize individualized therapy for patients undergoing BNCT.

Understanding hemoglobin contribution to high-dose methotrexate disposition-population pharmacokinetics in pediatric patients with hematological malignancies.

The aim of the present study was to develop a population pharmacokinetic model for methotrexate (MTX) during high-dose treatment (HDMTX) in pediatric patients with acute lymphoblastic leukemia (ALL) and non-Hodgkin's lymphoma (NHL) and to describe the influence of variability factors. The study included 50 patients of both sexes (aged 1-18years) who received 3 or 5g/m2 of HDMTX. A nonlinear mixed effect modeling approach was applied for data analysis. Parameter estimation was performed by first-order conditional estimation method with interaction (FOCEI), whereas stepwise covariate modeling was used to assess variability factors. The final model is a two-compartment model that incorporates the effect of body surface area and the influence of hemoglobin and serum creatinine on MTX clearance (CL). Population pharmacokinetic values for a typical subject were estimated at 5.75 L/h/m2 for clearance (CL), 21.3 L/m2 for volume of the central compartment (V1), 8.2 L/m2 for volume of the peripheral compartment (V2), and 0.087 L/h/m2 for intercompartmental clearance (Q). According to the final model, MTX CL decreases with increasing serum creatinine, whereas a positive effect was captured for hemoglobin. A difference of almost 32% in MTX CL was observed among patients' hemoglobin values reported in the study. The developed population pharmacokinetic model can contribute to the therapy optimization during HDMTX in pediatric patients with ALL and NHL. In addition to renal function and body weight, it describes the influence of hemoglobin on CL, allowing better understanding of its contribution to the disposition of HDMTX.

Radial increment dynamics of Maritime pine (Pinus pinaster Ait.) in pure and mixed stands with Scots pine (Pinus sylvestris L.) under changing environmental conditions

Forest dynamics are undergoing profound alteration due to the fact that climate change is increasing the frequency and severity of abiotic and biotic forest disturbances. Because of the unpredictable nature of the drought periods and the variation in their severity, Mediterranean forests are typically more vulnerable. Mediterranean Pinus pinaster forests are showing decay symptoms linked to climate change. There is clear evidence that promoting mixtures can serve as an effective forest adaptation strategy. In this regard, we sought to better understand the responses of P. pinaster radial growth dynamics to various factors, in both mixed and pure forest stands, and provide valuable insights into P. pinaster dynamics when mixed with P. sylvestris. In this study, we examined the variation of intra-annual cumulative radial increment patterns in response to the climate of P. pinaster between pure and mixed stands with P. sylvestris. Using data from band dendrometers collected over five consecutive climatically distinct years (2016–2020), a nonlinear mixed-effect model approach was used to analyze the differences in intra-annual cumulative radial increment patterns for P. pinaster between years in mixed and pure stands. The intra-annual radial increment pattern of P. pinaster showed significant year-to-year variation and varied with tree size, with greater increment in larger trees. Trees in mixed stands had a higher mean radial increment compared to corresponding ones in pure stands. Increased summer maximum temperatures negatively affected tree cumulative annual increment regardless of composition, but with a lower impact on trees in pure stands. Spring precipitation increased the length of the growing season, while higher spring maximum temperatures triggered an earlier inflection point. Our results highlight the high plasticity of P. pinaster in adapting to varying intra- and inter-annual environmental conditions and competition with other species and suggest that promoting mixtures with P. sylvestris may be an interesting management strategy for adaptation to climate change.

Population Pharmacokinetic Analysis of Atomoxetine and its Metabolites in Children and Adolescents with Attention-Deficit/Hyperactivity Disorder.

Atomoxetine (ATX) is a non-stimulant used to treat attention-deficit/hyperactivity disorder (ADHD) and systemic exposure is highly variable due to polymorphic cytochrome P450 2D6 (CYP2D6) activity. The objective of this study was to characterize the time course of ATX and metabolites (4-hydroxyatomoxetine (4-OH); N-desmethylatomoxetine (NDA); and 2-carboxymethylatomoxetine (2-COOH)) exposure following oral ATX dosing in children with ADHD to support individualized dosing. A nonlinear mixed-effect modeling approach was used to analyze ATX, 4-OH, and NDA plasma and urine, and 2-COOH urine profiles obtained over 24-72 hours from children with ADHD (n = 23) following a single oral ATX dose. Demographics and CYP2D6 activity score (AS) were evaluated as covariates. Simulations were performed to explore the ATX dosing in subjects with various CYP2D6 AS. A simultaneous pharmacokinetic modeling approach was used in which a model for ATX, 4-OH, and NDA in plasma and urine, and 2-COOH in urine was developed. Plasma ATX, 4-OH, and NDA were modeled using two-compartment models with first-order elimination. CYP2D6 AS was a significant determinant of ATX apparent oral clearance (CL/F), fraction metabolized to 4-OH, and systemic exposure of NDA. CL/F of ATX varied almost 7-fold across the CYP2D6 AS groups: AS 2: 20.02 L/hour; AS 1: 19.00 L/hour; AS 0.5: 7.47 L/hour; and AS 0: 3.10 L/hour. The developed model closely captures observed ATX, 4-OH, and NDA plasma and urine, and 2-COOH urine profiles. Application of the model shows the potential for AS-based dosing recommendations for improved individualized dosing.

Development and internal validation of a clinical prediction model for individualized dosing of BAY 81-8973, A full-length recombinant factor VIII, in pediatric patients with haemophilia A

BackgroundAs the most commonly used coagulation factor VIII (FVIII) concentrate in China, the individualized dosing prediction model of Kovaltry (BAY81–8973) is not fully investigated in pediatric patients. The prophylaxis tailored by population pharmacokinetic (PopPK) model can optimize dosing regimens. ObjectivesThis study aimed to develop PopPK models of BAY 81–8973 in pediatric patients, identify quantitative relationships of blood type (as a substitution for von Willebrand factor) on FVIII clearance and provide model-informed precision dosing (MIPD) procedures. MethodsPediatric patients with severe hemophilia A were enrolled and PK tests were conducted. The blood samples were collected at six time point. One-stage-based activated partial thromboplastin time was used for FVIII activity. Basic demographics and key covariates (blood type and von Willebrand factor antigen) were collected. A nonlinear mixed-effect modeling (NONMEM) approach was employed to establish PopPK model. Simulations were performed to evaluate current dosing regimens and present MIPD strategies. ResultsA total of 30 pediatric patients were included in the analysis. In the final model, Fat-free mass calculated from weight, age, and height was included as a size descriptor which affect FVIII apparent volume of distribution and clearance. Both von Willebrand factor antigen (VWF:Ag) and blood type accounted for the interindividual variability of FVIII clearance, but only one can retain in the final model. Therefore, two PopPK models based on VWF:Ag or blood type were developed. When the VWF:Ag value is doubled, the FVIII clearance is reduced by 35 %. Compared with blood type non-O, the clearance in pediatric patients with blood type O increased by 25.9 %. Weight-based dosing without regard to age and blood type resulted in large differences in FVIII trough activity. Patient demographics, dosing information, sparse blood samples and PopPK model together with Bayesian estimate constituted the MIPD workflow. Using it, the individual parameters and optimized dosing regimen could be achieved. ConclusionThis is the first predictive model designed to predict individualized dosing of BAY 81–8973 in pediatric patients with Hemophilia A. These results are useful in the PK-guided prophylaxis among pediatric patients and hold great potential to improve their long-term clinical outcomes.

Population pharmacokinetic modeling and dosing simulation of avalglucosidase alfa for selecting alternative dosing regimen in pediatric patients with late-onset pompe disease

Avalglucosidase alfa (AVAL) was approved in the United States (2021) for patients with late-onset Pompe disease (LOPD), aged ≥ 1 year. In the present study, pharmacokinetic (PK) simulations were conducted to propose alternative dosing regimens for pediatric LOPD patients based on a bodyweight cut-off. Population PK (PopPK) analysis was performed using nonlinear mixed effect modeling approach on pooled data from three clinical trials with LOPD patients, and a phase 2 study (NCT03019406) with infantile-onset Pompe disease (IOPD: 1–12 years) patients. A total of 2257 concentration-time points from 91 patients (LOPD, n = 75; IOPD, n = 16) were included in the analysis. The model was bodyweight dependent allometric scaling with time varying bodyweight included on clearance and distribution volume. Simulations were performed for two dosing regimens (20 mg/kg or 40 mg/kg) with different bodyweight cut-off (25, 30, 35 and 40 kg) by generating virtual pediatric (1–17 years) and adult patients. Corresponding simulated individual exposures (maximal concentration, Cmax and area under the curve in the 2-week dosing interval, AUC2W), and distributions were calculated. It was found that dosing of 40 mg/kg and 20 mg/kg in pediatric patients < 30 kg and ≥ 30 kg, respectively, achieved similar AVAL exposure (based on AUC2W) to adult patients receiving 20 mg/kg. PK simulations conducted on the basis of this model provided supporting data for the currently approved US labelling for dosing adapted bodyweight in LOPD patients ≥ 1 year by USFDA.

Vancomycin population pharmacokinetics analysis in Chinese paediatric patients with varying degrees of renal function and ages: development of new practical dosing recommendations.

To describe the pharmacokinetics of vancomycin in a large Chinese paediatric cohort with varying degrees of renal function and ages and to develop practical dosing guidelines. We conducted a retrospective population pharmacokinetic study using data from paediatric patients who received vancomycin between June 2013 and June 2022. A non-linear mixed-effect modelling approach with a one-compartment model structure was applied. Monte Carlo simulations were used to stimulate an optimal dosage regimen to achieve the target of AUC24/MIC between 400 and 650. We analysed a total of 673 paediatric patients and 1547 vancomycin serum concentrations. Covariate analysis revealed that physiological maturation, renal function, albumin and cardiothoracic surgery (CTS) significantly affected vancomycin pharmacokinetics. The typical clearance and volume of distribution, standardized to 70 kg, were 7.75 L/h (2.3% relative standard error, RSE) and 36.2 L (1.7% RSE), respectively. Based on the model, we proposed an optimal dosing regimen that considers the patient's age and estimate glomerular filtration rate (eGFR) to achieve a target AUC24/MIC for CTS and non-CTS patients. We also found that a loading dose of 20 mg/kg can help patients with an eGFR of <60 mL/min/1.73 m2 achieve the target AUC on the first day of treatment. We established vancomycin pharmacokinetic parameters in Chinese paediatric patients and proposed a dosing guideline integrating eGFR, age and CTS status, potentially improving clinical outcomes and reducing nephrotoxicity risk.

Population Pharmacokinetics of Total Rabbit Anti-thymocyte Globulin in Non-obese Adult Patients Undergoing Hematopoietic Cell Transplantation for Hematologic Malignancy.

Rabbit anti-thymocyte globulin (rATG), a therapeutic polyclonal antibody against human T cells, is commonly used in conditioning therapy prior to allogeneic hematopoietic cell transplantation (HCT). Previous studies successfully developed an individualized rATG dosing regimen based on "active" rATG population PK (popPK) analysis, while "total" rATG can be a more logistically favorable alternative for early HCT outcomes. We conducted a novel popPK analysis of total rATG. Total rATG concentration was measured in adult human-leukocyte antigen (HLA) mismatched HCT patients who received a low-dose rATG regimen (total 2.5-3mg/kg) within 3days prior to HCT. PopPK modeling and simulation was performed using nonlinear mixed effect modeling approach. A total of 504 rATG concentrations were available from 105 non-obese patients with hematologic malignancy (median age 47years) treated in Japan. The majority had acute leukemia or malignant lymphoma (94%). Total rATG PK was described by a two-compartment linear model. Influential covariate relations include ideal body weight [positively on both clearance (CL) and central volume of distribution], baseline serum albumin (negatively on CL), CD4+ T cell dose (positively on CL), and baseline serum IgG (positively on CL). Simulated covariate effects predicted that early total rATG exposures were affected by ideal body weight. This novel popPK model described the PK of total rATG in the adult HCT patients who received a low-dose rATG conditioning regimen. This model can be used for model-informed precision dosing in the settings with minimal baseline rATG targets (T cells), and early clinical outcomes are of interest.

A pharmacokinetic analysis of amisulpride in adult Chinese patients with schizophrenia: Impact of creatinine clearance.

To develop a stable population pharmacokinetic (PPK) model of amisulpride and to investigate the effects of covariates on the pharmacokinetic parameters in adult Chinese patients with schizophrenia. This retrospective study was carried out using 168 serum samples from 88 patients collected during routine clinical monitoring. Covariates recorded included demographic parameters (gender, age, weight), clinical parameters (serum creatinine, creatinine clearance), and intake of co-medications. The amisulpride PPK model was established using a nonlinear mixed effects modeling (NONMEM) approach. Goodness-of-fit (GOF) plots, bootstrap validation (1,000 runs), and normalized prediction distribution error (NPDE) were used in the evaluation of the final model. A one-compartment model with first-order absorption and elimination was developed. The population estimates for apparent clearance (CL/F) and apparent volume of distribution (V/F) were 32.6 L/h and 391 L, respectively. Estimated creatinine clearance (eCLcr) was a significant covariate for CL/F. The established model was: CL/F = 32.6 × (eCLcr/114.3)0.485 (L/h). The stability of the model was confirmed using GOF plots, bootstrap, and NPDE. Creatinine clearance is a major covariate which is positively correlated with CL/F. Therefore, additional dose adjustments of amisulpride may be required on the basis of eCLcr. An ethnic difference may exist in the pharmacokinetics of amisulpride, but further research is needed in order to confirm this possibility. The PPK model of amisulpride for adult Chinese schizophrenic patients established here using NONMEM, is potentially an important tool for individualizing drug dosage and therapeutic drug monitoring.

New Recommendations of a Height-Based Dosing Regimen of Tobramycin for Cystic Fibrosis in Adults: A Population Pharmacokinetic Analysis.

Acute pulmonary exacerbations (APEs) in patients with adult cystic fibrosis (CF) are treated with a beta-lactam and an aminoglycoside for activity against Pseudomonas aeruginosa (PA). Emerging drug resistance and changing pharmacokinetic profile in an aging population involve a reevaluation of tobramycin dosing recommendations. The objective of this study was to develop a population pharmacokinetic model and establish optimal dosing recommendations for tobramycin using Monte Carlo simulations. This retrospective clinical study and data collection were performed at the CF center of the McGill University Health Center (MUHC), Canada. Model development and simulations were performed using a nonlinear mixed-effect modeling approach (NONMEM, version 7.4.2). The ratios of maximal concentration (C max ) to the minimal inhibitory concentration (MIC) (C max /MIC ≥8 and ≥10) and area under the curve (AUC) to the MIC (AUC/MIC ≥70 and ≥100) were evaluated. Adult patients with CF (n = 51) treated with tobramycin were included in the study. Plasma concentrations of tobramycin were obtained for 699 samples from the MUHC database. The two-compartmental model best described the pharmacokinetics of tobramycin. The association of patient height with the central volume of distribution significantly improved this model. Height, rather than weight, induced the best reduction in objective function. According to simulations, doses between 3.4 mg/cm and 4.4 mg/cm were necessary to achieve C max /MIC values of ≥8 and ≥10, respectively. However, higher doses were required to achieve the AUC/MIC targets. This study demonstrated that height of the patients seems to be more suitable than their weight for dosing adjustments in adult patients with CF. According to this model, initial doses of tobramycin between 3.4 and 4.4 mg/cm should be recommended for patients with a median height of 164 cm and weight of 55 kg to achieve the target plasma concentrations.

Safety, pharmacokinetics and exploratory exposure-response analysis of CX3002, a novel inhibitor of Xa, in Chinese healthy subjects

Background and objectiveCX3002 is a structurally novel inhibitor of factor Xa, with promising prospects. This study aims to report the results of a first-in-human ascending-dose study of CX3002 in Chinese healthy subjects, and to establish an exploratory population pharmacokinetic/pharmacodynamic (PK/PD) model to investigate the exposure-response relationship of CX3002. MethodsThe randomized, double-blind, placebo-controlled study included six single-dose groups and three multiple-dose groups, with a dose range of 1–30 mg. Safety, tolerability, pharmacokinetics (PK) and pharmacodynamics (PD) of CX3002 were evaluated. PK of CX3002 was analyzed using both non-compartment method and population modeling. PK/PD model was developed using nonlinear mixed effect modeling approach and was evaluated by prediction-corrected visual predictive check and bootstrap methods. ResultsA total of 84 subjects were enrolled and all participants completed the study. CX3002 exhibited satisfactory safety and tolerability in healthy subjects. Cmax and AUC of CX3002 increased with dose from 1 to 30 mg, but less-than-proportional increases were observed. There was no obvious accumulation with multiple doses. Anti-Xa activity showed dose-related increases after administration of CX3002 but not placebo. The PK of CX3002 was well described by a two-compartment model with a modification of bioavailability according to dose, and anti-Xa activity was described by a Hill function. No covariate was identified significant based on the limited data in this study. ConclusionsCX3002 was well tolerated and resulted in dose-related anti-Xa activity across the dose range. The PK of CX3002 were predictable, and correlated with PD effects. Continued clinical investigation of CX3002 was supported. Chinadrugtrials.org.cn identifier: CTR20190153.

Investigating the contribution of residual unexplained variability components on bias and imprecision of parameter estimates in population pharmacokinetic mixed-effects modeling.

In a nonlinear mixed-effects modeling (NLMEM) approach of pharmacokinetic (PK) and pharmacodynamic (PD) data, two levels of random effects are generally modeled: between-subject variability (BSV) and residual unexplained variability (RUV). The goal of this simulation-estimation study was to investigate the extent to which PK and RUV model misspecification, errors in recording dosing and sampling times, and variability in drug content uniformity contribute to the estimated magnitude of RUV and PK parameter bias. A two-compartment model with first-order absorption and linear elimination was simulated as a true model. PK parameters were clearance 5.0 L/h; central volume of distribution 35 L; inter-compartmental clearance 50 L/h; peripheral volume of distribution 50 L. All parameters were assumed to have a 30% coefficient of variation (CV). One hundred in-silico subjects were administered a labeled dose of 120mg under 4 sample collection designs. PK and RUV model misspecifications were associated with relatively larger increases in the magnitude of RUV compared to other sources for all levels of sampling design. The contribution of dose and dosing time misspecifications have negligible effects on RUV but result in higher bias in PK parameter estimates. Inaccurate sampling time data results in biased RUV and increases with the magnitude of perturbations. Combined perturbation scenarios in the studied sources will propagate the variability and accumulate in RUV magnitude and bias in PK parameter estimates. This work provides insight into the potential contributions of many factors that comprise RUV and bias in PK parameters.

Population pharmacokinetics of ruxolitinib in children with hemophagocytic lymphohistiocytosis: focus on the drug-drug interactions.

The optimal dose regimen of ruxolitinib (RUX) in children with hemophagocytic lymphohistiocytosis (HLH) remains to be determined. The aim was to develop and verify a population pharmacokinetic (PPK) model, and then provide references for the optimization of dose regimen of RUX in children with HLH. A total of 189 RUX concentrations from 32 children were included. The PPK model was established using the nonlinear mixed-effects model approach. Predictive performance and stability of the final PPK model were evaluated. The exposure of RUX in different clinical scenarios was simulated through Monte Carlo simulations. A one-compartment model with first-order absorption and linear elimination was identified to describe the disposition of RUX. The absorption rate constant (Ka) in the final PPK model was 1.05h-1, and the apparent clearance (CL/F) and volume of distribution (V/F) were 9.80 L/h and 30.6 L, respectively. Coadministration with triazoles (TZS) and azithromycin (AZM) resulted in approximately 31.0% and 32.4% reductions in the CL/F of RUX, respectively. Multiple evaluation procedures showed satisfactory predictive performance and stability of the final model. Monte Carlo simulations showed that the exposure of RUX was significantly affected by the coadministration with TZS and/or AZM under different clinical scenarios. For the first time, a PPK model of RUX in children with HLH was developed and evaluated. The coadministration with TZS and/or AZM were found to reduce the clearance of RUX in children. These findings could provide new insights for the precise treatment of RUX in children with HLH.

Joint population pharmacokinetic modeling of venlafaxine and O-desmethyl venlafaxine in healthy volunteers and patients to evaluate the impact of morbidity and concomitant medication.

Introduction: Venlafaxine (VEN) is a widely used dual selective serotonin/noradrenaline reuptake inhibitor indicated for depression and anxiety. It undergoes first-pass metabolism to its active metabolite, O-desmethyl venlafaxine (ODV). The aim of the present study was to develop a joint population pharmacokinetic (PPK) model to characterize their pharmacokinetic characters simultaneously. Methods: Plasma concentrations with demographic and clinical data were derived from a bioequivalence study in 24 healthy subjects and a naturalistic TDM setting containing 127 psychiatric patients. A parent-metabolite PPK modeling was performed with NONMEM software using a non-linear mixed effect modeling approach. Goodness of fit plots and normalized prediction distribution error method were used for model validation. Results and conclusion: Concentrations of VEN and ODV were well described with a one-compartment model incorporating first-pass metabolism. The first-pass metabolism was modeled as a first-order conversion. The morbid state and concomitant amisulpride were identified as two significant covariates affecting the clearance of VEN and ODV, which may account for some of the variations in exposure. This model may contribute to the precision medication in clinical practice and may inspire other drugs with pre-system metabolism.

A Population Pharmacokinetic Analysis of Emicizumab in Patients with Acquired Hemophilia a Enrolled in a Phase III Study (AGEHA)

INTRODUCTION: Emicizumab is a bispecific antibody mimicking the cofactor function of activated factor VIII. A population pharmacokinetic (PK) analysis was performed to characterize the PK of emicizumab in patients with acquired hemophilia A (PwAHA) leveraging a model that had been developed in patients with congenital hemophilia A (PwCHA) [1]. METHODS: Plasma emicizumab concentration data from 12 PwAHA, who received emicizumab subcutaneously in a phase III study (AGEHA) and were included in its primary analysis [2], were analyzed by a nonlinear mixed-effect modeling approach. A linear one-compartment model with first-order absorption and elimination was employed for the structural model. Covariate investigation was limited to minimal given the possibly poor power for identifying right covariates and their right effects due to the much smaller sample size than in PwCHA (N = 381). A full covariate modeling approach was applied to include and infer the covariate effects that had been identified to be statistically significant in PwCHA, i.e., body weight (BW) on the apparent clearance (CL/F) and apparent volume of distribution (V/F), age on the bioavailability (parameterized as relative bioavailability [Frel] to a typical age), and albumin level (ALB) on the CL/F, except for race on the V/F because all patients in this analysis were Asian, without statistical testing for inclusion or exclusion. Because the age and ALB (50 to 92 years and 20 to 40 g/L, respectively) included some values out of the ranges that had been evaluated in PwCHA (1.22 to 77 years and 33 to 57 g/L, respectively), these covariate effects were separately estimated for ≤ 77 and > 77 years and for < 33 and ≥ 33 g/L, respectively, without statistical testing for separation or integration. Other covariate effects that had been identified not to be statistically significant in PwCHA were not repeatedly explored in this analysis. A stepwise covariate modeling approach was applied to explore the effects of a potential new covariate of sex, which had not been tested in PwCHA because all patients were male, on the CL/F, V/F, and absorption rate constant (KA) with statistical testing for inclusion and exclusion. Effects of renal and hepatic functions (normal to severely impaired and normal to moderately impaired, respectively; no patients had severe hepatic impairment) on the CL/F, V/F, and KA were explored graphically in a post hoc manner. Parameter estimates were corrected as appropriate when compared with PwCHA. RESULTS: In PwAHA, all prior identified covariate effects within the ranges previously evaluated in PwCHA exhibited consistent directions (positive or negative) in the point estimates with those in PwCHA, and their 95% confidence intervals (CIs) involved the point estimates in PwCHA. The effects of age > 77 years and ALB < 33 g/L had 95% CIs involving not only the point estimates for ≤ 77 years and ≥ 33 g/L, respectively, but also the null value. Effects of sex were not statistically significant. No consistent effects were observed with the severity of renal or hepatic impairment. The 95% CIs of the CL/F/Frel and V/F/Frel in PwAHA involved the point estimates in PwCHA regardless of the presence or absence of covariate effect correction, indicating no clear effects of being a PwAHA itself. The 95% CI of the KA in PwAHA fell below the point estimate in PwCHA, indicating a possibly slower absorption than in PwCHA; however, its impact on the steady-state exposure was predicted to be limited. Simulations validated the observed achievement of the efficacious mean steady-state trough concentration of > 30 μg/mL with the tested maintenance dose of 1.5 mg/kg once weekly in AGEHA as in PwCHA [2,3], supporting the appropriateness of applying the maintenance dose for PwAHA. CONCLUSIONS: No remarkable differences from PwCHA in the covariates or their effects on the PK of emicizumab were indicated in PwAHA. This finding would warrant application of the same covariate-based dose adjustment as in PwCHA (i.e., only BW is considered) for PwAHA. REFERENCES: 1. Clin Pharmacokinet 2020;59:1611-25. 2. ISTH 2022 Congress. Abstract OC 40.2. 3. Clin Pharmacokinet 2021;60:931-41.

Meropenem Population Pharmacokinetics and Dosing Regimen Optimization in Critically Ill Children Receiving Continuous Renal Replacement Therapy.

We aimed to develop a meropenem population pharmacokinetic model in critically ill children receiving continuous renal replacement therapy and simulate dosing regimens to optimize patient exposure. Meropenem plasma concentration was quantified by high-performance liquid chromatography. Meropenem pharmacokinetics was investigated using a non-linear mixed-effect modeling approach. Monte Carlo simulations were performed to compute the optimal scheme of administration, according to the target of a 100% inter-dose interval time in which concentration is one to four times above the minimum inhibitory concentration (100% fT>1-4×MIC). A total of 27 patients with a median age of 4 [interquartile range 0-11] years, a median body weight of 16 [range 7-35] kg receiving continuous renal replacement therapy were included. Concentration-time courses were best described by a one-compartment model with first-order elimination. Body weight (BW) produced significant effects on volume of distribution (V) and BW and continuous renal replacement therapy effluent flow rate (Qeff) produced significant effects on clearance (CL): [Formula: see text] and [Formula: see text], where Vpop and CLpop estimates were 32.5 L and 5.88 L/h, respectively, normalized to a 70-kg BW and median Qeff at 1200 mL/h. Using this final model and Monte Carlo simulations, for patients with Qeff over 1200 mL/h, meropenem continuous infusion was adequate in most cases to attain 100% fT>1-4xMIC. For bacterial infections with a low minimum inhibitory concentration (≤2 mg/L), meropenem intermitent administration was appropriate for patients weighing more than 20 kg with Qeff <500 mL/h and for patients weighing more than 10 kg with Qeff <100 mL/h. Meropenem exposure in critically ill children receiving continuous renal replacement therapy needs dosing adjustments to the minimum inhibitory concentration that take into account body weight and the continuous renal replacement therapy effluent flow rate.

Population Pharmacokinetics and Dose Optimization Based on Renal Function of Rivaroxaban in Thai Patients with Non-Valvular Atrial Fibrillation.

Low-dose rivaroxaban has been used in Asian patients with direct oral anticoagulants (DOACs) eligible for atrial fibrillation (AF). However, there are few pharmacokinetic (PK) data in Thai patients to support precise dosing. This study aimed to develop a population PK model and determine the optimal rivaroxaban doses in Thai patients. A total of 240 Anti-Xa levels of rivaroxaban from 60 Thai patients were analyzed. A population PK model was established using the nonlinear mixed-effect modeling approach. Monte Carlo simulations were used to predict drug exposures at a steady state for various dosages. Proportions of patients having rivaroxaban exposure within typical exposure ranges were determined. A one-compartment model with first-order absorption best described the data. Creatinine clearance (CrCl) and body weight significantly affected CL/F and V/F, respectively. Regardless of body weight, a higher proportion of patients with CrCl < 50 mL/min receiving the 10-mg once-daily dose had rivaroxaban exposures within the typical exposure ranges. In contrast, a higher proportion of patients with CrCl ≥ 50 mL/min receiving the 15-mg once-daily dose had rivaroxaban exposures within the typical exposure ranges. The study’s findings suggested that low-dose rivaroxaban would be better suited for Thai patients and suggested adjusting the medication’s dose in accordance with renal function.

Single-time-point estimation of absorbed doses in PRRT using a non-linear mixed-effects model

IntroductionEstimation of accurate time-integrated activity coefficients (TIACs) and radiation absorbed doses (ADs) is desirable for treatment planning in peptide-receptor radionuclide therapy (PRRT). This study aimed to investigate the accuracy of a simplified dosimetry using a physiologically-based pharmacokinetic (PBPK) model, a nonlinear mixed effect (NLME) model, and single-time-point imaging to calculate the TIACs and ADs of 90Y-DOTATATE in various organs of dosimetric interest and tumors. Materials & MethodsBiokinetic data of 111In-DOTATATE in tumors, kidneys, liver, spleen, and whole body were obtained from eight patients using planar scintigraphic imaging at T1 = (2.9 ± 0.6), T2 = (4.6 ± 0.4), T3 = (22.8 ± 1.6), T4 = (46.7 ± 1.7) and T5 = (70.9 ± 1.0) h post injection. Serum activity concentration was measured at 5 and 15 min; 0.5, 1, 2, and 4 h; and 1, 2, and 3 d p.i.. A published PBPK model for PRRT, NLME, and a single-time-point imaging datum at different time points were used to calculate TIACs in tumors, kidneys, liver, spleen, whole body, and serum. Relative deviations (RDs) (median [min, max]) between the calculated TIACs from single-time-point imaging were compared to the TIACs calculated from the all-time-points fit. The root mean square error (RMSE) of the difference between the computed ADs from the single-time-point imaging and reference ADs from the all-time point fittings were analyzed. A joint root mean square error RMSEjoint of the ADs was calculated with the RSME from both the tumor and kidneys to sort the time points concerning accurate results for the kidneys and tumor dosimetry. The calculations of TIACs and ADs from the single-time-point dosimetry were repeated using the sum of exponentials (SOE) approach introduced in the literature. The RDs and the RSME of the PBPK approach in our study were compared to the SOE approach. ResultsUsing the PBPK and NLME models and the biokinetic measurements resulted in a good fit based on visual inspection of the fitted curves and the coefficient of variation CV of the fitted parameters (<50%). T4 was identified being the time point with a relatively low median and range of TIACs RDs, i.e., 5 [1, 21]% and 2 [−15, 21]% for kidneys and tumors, respectively. T4 was found to be the time point with the lowest joint root mean square error RMSEjoint of the ADs. Based on the RD and RMSE, our results show a similar performance as the SOE and NLME model approach. SummaryIn this study, we introduced a simplified calculation of TIACs/ADs using a PBPK model, an NLME model, and a single-time-point measurement. Our results suggest a single measurement might be used to calculate TIACs/ADs in the kidneys and tumors during PRRT.