Assessing potential drug-drug interactions between clofazimine and other frequently used agents to treat drug-resistant tuberculosis

Clofazimine is included in drug regimens to treat rifampicin/drug-resistant tuberculosis (DR-TB), but there is little information about its interaction with other drugs in DR-TB regimens. We evaluated the pharmacokinetic interaction between clofazimine and isoniazid, linezolid, levofloxacin, and cycloserine, dosed as terizidone. Newly diagnosed adults with DR-TB at Klerksdorp/Tshepong Hospital, South Africa, were started on the then-standard treatment with clofazimine temporarily excluded for the initial 2 weeks. Pharmacokinetic sampling was done immediately before and 3 weeks after starting clofazimine, and drug concentrations were determined using validated liquid chromatography-tandem mass spectrometry assays. The data were interpreted with population pharmacokinetics in NONMEM v7.5.1 to explore the impact of clofazimine co-administration and other relevant covariates on the pharmacokinetics of isoniazid, linezolid, levofloxacin, and cycloserine. Clofazimine, isoniazid, linezolid, levofloxacin, and cycloserine data were available for 16, 27, 21, 21, and 6 participants, respectively. The median age and weight for the full cohort were 39 years and 52 kg, respectively. Clofazimine exposures were in the expected range, and its addition to the regimen did not significantly affect the pharmacokinetics of the other drugs except levofloxacin, for which it caused a 15% reduction in clearance. A posteriori power size calculations predicted that our sample sizes had 97%, 90%, and 87% power at P < 0.05 to detect a 30% change in clearance of isoniazid, linezolid, and cycloserine, respectively. Although clofazimine increased the area under the curve of levofloxacin by 19%, this is unlikely to be of great clinical significance, and the lack of interaction with other drugs tested is reassuring.

The development and evaluation of dose-prediction tools for allopurinol therapy (Easy-Allo tools)

AIMS

Dose escalation at the initiation of allopurinol therapy can be protracted and resource intensive. Tools to predict the allopurinol doses required to achieve target serum urate concentrations would facilitate the implementation of more efficient dose-escalation strategies. The aim of this research was to develop and externally evaluate allopurinol dosing tools, one for use when the pre-urate-lowering therapy serum urate is known (Easy-Allo1) and one for when it is not known (Easy-Allo2).

METHODS

A revised population pharmacokinetic-pharmacodynamic model was developed using data from 653 people with gout. Maintenance doses to achieve the serum urate target of <0.36 mmol L-1 in >80% of individuals were simulated and evaluated against external data. The predicted and observed allopurinol doses were compared using the mean prediction error (MPE) and root mean square error (RMSE). The proportion of Easy-Allo predicted doses within 100 mg of the observed was quantified.

RESULTS

Allopurinol doses were predicted by total body weight, baseline urate, ethnicity and creatinine clearance. Easy-Allo1 produced unbiased and suitably precise dose predictions (MPE 2 mg day-1 95% confidence interval [CI] -13-17, RMSE 91%, 90% within 100 mg of the observed dose). Easy-Allo2 was positively biased by about 70 mg day-1 and slightly less precise (MPE 70 mg day-1 95% CI 52-88, RMSE 131%, 71% within 100 mg of the observed dose).

CONCLUSIONS

The Easy-Allo tools provide a guide to the allopurinol maintenance dose requirement to achieve the serum urate target of <0.36 mmol L-1 and will aid in the development of novel dose-escalation strategies for allopurinol therapy.

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

PURPOSE

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.

METHODS

The study included 50 patients of both sexes (aged 1-18 years) who received 3 or 5 g/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.

RESULTS

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.

CONCLUSION

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.

Evaluation of ChatGPT and Gemini large language models for pharmacometrics with NONMEM

To assess ChatGPT 4.0 (ChatGPT) and Gemini Ultra 1.0 (Gemini) large language models on NONMEM coding tasks relevant to pharmacometrics and clinical pharmacology. ChatGPT and Gemini were assessed on tasks mimicking real-world applications of NONMEM. The tasks ranged from providing a curriculum for learning NONMEM, an overview of NONMEM code structure to generating code. Prompts in lay language to elicit NONMEM code for a linear pharmacokinetic (PK) model with oral administration and a more complex model with two parallel first-order absorption mechanisms were investigated. Reproducibility and the impact of "temperature" hyperparameter settings were assessed. The code was reviewed by two NONMEM experts. ChatGPT and Gemini provided NONMEM curriculum structures combining foundational knowledge with advanced concepts (e.g., covariate modeling and Bayesian approaches) and practical skills including NONMEM code structure and syntax. ChatGPT provided an informative summary of the NONMEM control stream structure and outlined the key NONMEM Translator (NM-TRAN) records needed. ChatGPT and Gemini were able to generate code blocks for the NONMEM control stream from the lay language prompts for the two coding tasks. The control streams contained focal structural and syntax errors that required revision before they could be executed without errors and warnings. The code output from ChatGPT and Gemini was not reproducible, and varying the temperature hyperparameter did not reduce the errors and omissions substantively. Large language models may be useful in pharmacometrics for efficiently generating an initial coding template for modeling projects. However, the output can contain errors and omissions that require correction.

Quantification of biochemical PSA dynamics after radioligand therapy with [177Lu]Lu-PSMA-I&T using a population pharmacokinetic/pharmacodynamic model

BACKGROUND

There is an unmet need for prediction of treatment outcome or patient selection for [177Lu]Lu-PSMA therapy in patients with metastatic castration-resistant prostate cancer (mCRPC). Quantification of the tumor exposure-response relationship is pivotal for further treatment optimization. Therefore, a population pharmacokinetic (PK) model was developed for [177Lu]Lu-PSMA-I&T using SPECT/CT data and, subsequently, related to prostate-specific antigen (PSA) dynamics after therapy in patients with mCRPC using a pharmacokinetic/pharmacodynamic (PKPD) modelling approach.

METHODS

A population PK model was developed using quantitative SPECT/CT data (406 scans) of 76 patients who received multiple cycles [177Lu]Lu-PSMA-I&T (± 7.4 GBq with either two- or six-week interval). The PK model consisted of five compartments; central, salivary glands, kidneys, tumors and combined remaining tissues. Covariates (tumor volume, renal function and cycle number) were tested to explain inter-individual variability on uptake into organs and tumors. The final PK model was expanded with a PD compartment (sequential fitting approach) representing PSA dynamics during and after treatment. To explore the presence of a exposure-response relationship, individually estimated [177Lu]Lu-PSMA-I&T tumor concentrations were related to PSA changes over time.

RESULTS

The population PK model adequately described observed data in all compartments (based on visual inspection of goodness-of-fit plots) with adequate precision of parameters estimates (< 36.1% relative standard error (RSE)). A significant declining uptake in tumors (k14) during later cycles was identified (uptake decreased to 73%, 50% and 44% in cycle 2, 3 and 4-7, respectively, compared to cycle 1). Tumor growth (defined by PSA increase) was described with an exponential growth rate (0.000408 h-1 (14.2% RSE)). Therapy-induced PSA decrease was related to estimated tumor concentrations (MBq/L) using both a direct and delayed drug effect. The final model adequately captured individual PSA concentrations after treatment (based on goodness-of-fit plots). Simulation based on the final PKPD model showed no evident differences in response for the two different dosing regimens currently used.

CONCLUSIONS

Our population PK model accurately described observed [177Lu]Lu-PSMA-I&T uptake in salivary glands, kidneys and tumors and revealed a clear declining tumor uptake over treatment cycles. The PKPD model adequately captured individual PSA observations and identified population response rates for the two dosing regimens. Hence, a PKPD modelling approach can guide prediction of treatment response and thus identify patients in whom radioligand therapy is likely to fail.

Model-based precision dosing and remedial dosing recommendations for delayed or missed doses of isoniazid in Chinese patients with tuberculosis

AIMS

Isoniazid (INH) has been used as a first-line drug to treat tuberculosis (TB) for more than 50 years. However, large interindividual variability was found in its pharmacokinetics, and effects of nonadherence to INH treatment and corresponding remedy regime remain unclear. This study aimed to develop a population pharmacokinetic (PPK) model of INH in Chinese patients with TB to provide model-informed precision dosing and explore appropriate remedial dosing regimens for nonadherent patients.

METHODS

In total, 1012 INH observations from 736 TB patients were included. A nonlinear mixed-effects modelling was used to analyse the PPK of INH. Using Monte Carlo simulations to determine optimal dosage regimens and design remedial dosing regimens.

RESULTS

A 2-compartmental model, including first-order absorption and elimination with allometric scaling, was found to best describe the PK characteristics of INH. A mixture model was used to characterize dual rates of INH elimination. Estimates of apparent clearance in fast and slow eliminators were 28.0 and 11.2 L/h, respectively. The proportion of fast eliminators in the population was estimated to be 40.5%. Monte Carlo simulations determined optimal dosage regimens for slow and fast eliminators with different body weight. For remedial dosing regimens, the missed dose should be taken as soon as possible when the delay does not exceed 12 h, and an additional dose is not needed. delay for an INH dose exceeds 12 h, the patient only needs to take the next single dose normally.

CONCLUSION

PPK modelling and simulation provide valid evidence on the precision dosing and remedial dosing regimen of INH.

Population pharmacokinetics of a novel PCSK9 antisense oligonucleotide

AIMS

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.5 S-constrained ethyl motif (cET) chemistry and is conjugated to a triantennary N-acetylgalactosamine (GalNAc3) ligand for targeted hepatocyte uptake.

METHODS

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).

RESULTS

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.9 mg/h [12% relative standard error]; concentration yielding half-maximal elimination rate, 4.8 mg/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.

CONCLUSIONS

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.

A population pharmacokinetic model of methotrexate in Korean patients with haematologic malignancy

AIMS

This study was conducted to develop a population pharmacokinetic (PK) model of methotrexate in Korean patients with haematologic malignancy, identify factors affecting methotrexate PK, and propose an optimal dosage regimen for the Korean population.

METHODS

Data were retrospectively collected from 188 patients with acute leukaemia or non-Hodgkin's lymphoma who were admitted to Severance Hospital during the period from November 2005 to January 2016. Using demographic factors and laboratory results as potential covariates for PK parameters, model development was performed using NONMEM and optimal dosing regimens were developed using the final PK model.

RESULTS

A two-compartment model incorporating body weight via allometry best described the data, yielding typical parameter values of 25.09 L for central volume of distribution ( V 1 ), 17.65 L for peripheral volume of distribution ( V 2 ), 12.89 L/h for clearance (CL) and 0.655 L/h for inter-compartmental clearance in a 50 kg patient. Covariate analyses showed that, at the weight of 50 kg, CL decreased by 0.11 L/h for each 1-year increase in age above 14 years old and decreased 0.8-fold when serum creatinine level doubled, indicating the importance of age-specific dose individualization in methotrexate treatment. Volume of distribution at steady state derived from PK parameters (=  V 1  +  V 2 ) was 0.85 L/kg, which was similar to those in the Western or Chinese populations. Optimal doses simulated from the final model successfully produced the PK measures close to the target chosen.

CONCLUSIONS

The population PK model and optimal dosage regimens developed in this study can be used as a basis to achieve precision dosing in Korean patients with haematologic malignancy.

The cycle effect quantified: reduced tumour uptake in subsequent cycles of [177Lu]Lu-HA-DOTATATE during peptide receptor radionuclide therapy

BACKGROUND

Clear evidence regarding the effect of reduced tumour accumulation in later peptide receptor radionuclide therapy (PRRT) cycles is lacking. Therefore, we aimed to quantify potential cycle effects for patients treated with [177Lu]Lu-HA-DOTATATE using a population pharmacokinetic (PK) modelling approach.

METHODS

A population PK model was developed using imaging data from 48 patients who received multiple cycles of [177Lu]Lu-HA-DOTATATE. The five-compartment model included a central, kidney, spleen, tumour and lumped rest compartment. Tumour volume and continued use of long-acting somatostatin analogues (SSAs) were tested as covariates in the model. In addition, the presence of a cycle effect was evaluated by relating the uptake rate in a specific cycle as a fraction of the (tumour or organ) uptake rate in the first cycle.

RESULTS

The final PK model adequately captured observed radioactivity accumulation in kidney, spleen and tumour. A higher tumour volume was identified to increase the tumour uptake rate, where a twofold increase in tumour volume resulted in a 2.3-fold higher uptake rate. Also, continued use of long-acting SSAs significantly reduced the spleen uptake rate (68.4% uptake compared to SSA withdrawal (10.5% RSE)). Lastly, a cycle effect was significantly identified, where tumour uptake rate decreased to 86.9% (5.3% RSE) in the second cycle and even further to 79.7% (5.6% RSE) and 77.6% (6.2% RSE) in the third and fourth cycle, respectively, compared to cycle one.

CONCLUSIONS

Using a population PK modelling approach, the cycle effect of reduced tumour uptake in subsequent PRRT cycles was quantified. Our findings implied that downregulation of target receptors is probably not the major cause of the cycle effect, due to a plateau in the decrease of tumour uptake in the fourth cycle.

Population pharmacokinetics of antibacterial agents in the older population: a literature review

INTRODUCTION

Older individuals face an elevated risk of developing bacterial infections. The optimal use of antibacterial agents in this population is challenging because of age-related physiological alterations, changes in pharmacokinetics (PK) and pharmacodynamics (PD), and the presence of multiple underlying diseases. Therefore, population pharmacokinetics (PPK) studies are of great importance for optimizing individual treatments and prompt identification of potential risk factors.

AREA COVERED

Our search involved keywords such as 'elderly,' 'old people,' and 'geriatric,' combined with 'population pharmacokinetics' and 'antibacterial agents.' This comprehensive search yielded 11 categories encompassing 28 antibacterial drugs, including vancomycin, ceftriaxone, meropenem, and linezolid. Out of 127 studies identified, 26 (20.5%) were associated with vancomycin, 14 (11%) with meropenem, and 14 (11%) with piperacillin. Other antibacterial agents were administered less frequently.

EXPERT OPINION

PPK studies are invaluable for elucidating the characteristics and relevant factors affecting the PK of antibacterial agents in the older population. Further research is warranted to develop and validate PPK models for antibacterial agents in this vulnerable population.

External evaluation of the predictive performance of published population pharmacokinetic models of linezolid in adult patients

OBJECTIVES

Several linezolid population pharmacokinetic (popPK) models have been established to facilitate optimal therapy; however, their extrapolated predictive performance to other clinical sites is unknown. This study aimed to externally evaluate the predictive performance of published pharmacokinetic models of linezolid in adult patients.

METHODS

For the evaluation dataset, 150 samples were collected from 70 adult patients (72.9% of which were critically ill) treated with linezolid at our center. Twenty-five published popPK models were identified from PubMed and Embase. Model predictability was evaluated using prediction-based, simulation-based, and Bayesian forecasting-based approaches to assess model predictability.

RESULTS

Prediction-based diagnostics found that the prediction error within ±30% (F30) was less than 40% in all models, indicating unsatisfactory predictability. The simulation-based prediction- and variability-corrected visual predictive check and normalized prediction distribution error test indicated large discrepancies between the observations and simulations in most of the models. Bayesian forecasting with one or two prior observations significantly improved the models' predictive performance.

CONCLUSION

The published linezolid popPK models showed insufficient predictive ability. Therefore, their sole use is not recommended, and incorporating therapeutic drug monitoring of linezolid in clinical applications is necessary.

A semimechanistic pharmacokinetic/pharmacodynamic model for alanine aminotransferase-based hepatotoxicity of methotrexate in paediatric patients with acute lymphoid leukaemia

AIMS

Methotrexate (MTX) is recognized for its potential to induce hepatotoxicity, commonly manifested by elevated alanine aminotransferase (ALT) levels. However, the quantitative relationship between the pharmacokinetics (PK) of MTX and ALT-based hepatotoxicity remains unclear. This study aimed to develop a semimechanistic PK/pharmacodynamic (PD) model to characterize the MTX-induced hepatotoxicity based on ALT in paediatric patients with acute lymphoid leukaemia.

METHODS

A retrospective study was conducted on paediatric patients who received high-dose (3-5 g/m2 ) MTX treatment. MTX concentrations were assessed at 24-h intervals until the concentration dropped below 0.1 μmol/L. ALT concentrations were measured both before and after MTX administration. A population PK model was initially developed, which was later connected to a semimechanistic hepatotoxicity model.

RESULTS

The PK model was developed using 354 MTX concentrations obtained from 51 patients, while the PD model was constructed using 379 ALT concentrations collected from 48 patients. The optimal PK model for MTX consisted of a 2-compartment structure, where body surface area served as a covariate for clearance and central volume of distribution. An indirect response model coupled to a liver injury signal transduction model was developed to describe the dynamics of ALT after MTX administration. The drug effect was adequately described by a linear model, exhibiting considerable interoccasion variability for each treatment session. No significant covariates were identified to have an impact on the PD parameters.

CONCLUSION

A semimechanistic model was developed to describe ALT-based hepatotoxicity of MTX, and it has the potential to serve as a valuable tool for characterizing drug-induced hepatotoxicity.

Population pharmacokinetics and model-based dosing optimization of teicoplanin in elderly critically ill patients with pneumonia

PURPOSE

To evaluate the population pharmacokinetics and pharmacodynamics of teicoplanin in elderly critically ill patients with pneumonia for optimal dosages.

METHODS

Fifteen critically ill patients (9 men) ≥ 60 years received teicoplanin 6 mg/kg for three doses followed by standard maintenance doses (6 mg/kg q24h) with renal dosing adjustment. Serial plasma samples from all patients were analyzed simultaneously by population pharmacokinetic modeling using NONMEM. Probability of target attainment (PTA) was calculated through Monte Carlo simulations for various dosing regimens to achieve adequate systemic exposures.

RESULTS

The median (interquartile range, IQR) age, body mass index, and creatinine clearance (CrCl) was 75 (64-78) years, 22.5 (20.8-25.4) kg/m2, and 64 (47-106) mL/min, respectively. The median (IQR) peak and trough concentration was 46.5 (42.7-51.0) and 8.7 (7.2-9.5) mg/L. The population pharmacokinetic model showed slower clearance (CL) and larger peripheral volume of distribution (V2) in patients with reduced CrCl: CL (L/h) = 0.629 × (CrCl/64)0.656, V2 (L) = 55.7 × (CrCl/64)-0.665. Model-based simulations showed PTAs ≥85% only for higher-dose regimens (12 mg/kg) up to an MIC of 0.5 mg/L.

CONCLUSIONS

Standard teicoplanin dosages for pneumonia may provide inadequate systemic exposures in elderly critically ill patients. High-dose regimens should be considered as empiric therapy or for less susceptible pathogens.

Abacavir Drug Exposures in African Children Under 14 kg Using Pediatric Solid Fixed Dose Combinations According to World Health Organization Weight Bands

BACKGROUND

The pharmacokinetics of abacavir (ABC) in African children living with HIV (CLHIV) weighing <14 kg and receiving pediatric fixed dose combinations (FDC) according to WHO weight bands dosing are limited. An ABC population pharmacokinetic model was developed to evaluate ABC exposure across different World Health Organization (WHO) weight bands.

METHODS

Children enrolled in the LIVING study in Kenya and Uganda receiving ABC/lamivudine (3TC) dispersible tablets (60/30 mg) according to WHO weight bands. A population approach was used to determine the pharmacokinetic parameters. Monte Carlo simulations were conducted using an in silico population with demographic characteristics associated with African CLHIV. ABC exposures (AUC0-24) of 6.4-50.4 mg h/L were used as targets.

RESULTS

Plasma samples were obtained from 387 children. A 1-compartment model with allometric scaling of clearance (CL/F) and volume of distribution (V/F) according to body weight best characterized the pharmacokinetic data of ABC. The maturation of ABC CL/F was characterized using a sigmoidal Emax model dependent on postnatal age (50% of adult CL/F reached by 0.48 years of age). Exposures to ABC were within the target range for children weighing 6.0-24.9 kg, but children weighing 3-5.9 kg were predicted to be overexposed.

CONCLUSIONS

Lowering the ABC dosage to 30 mg twice daily or 60 mg once daily for children weighing 3-5.9 kg increased the proportion of children within the target and provided comparable exposures. Further clinical study is required to investigate clinical implications and safety of the proposed alternative ABC doses.

Dolutegravir pharmacokinetics in Ugandan patients with TB and HIV receiving standard- versus high-dose rifampicin

Higher rifampicin doses may improve tuberculosis treatment outcomes. This could however exacerbate the existing drug interaction with dolutegravir. Moreover, the metabolism of dolutegravir may also be affected by polymorphism of UGT1A1, a gene that codes for uridine diphosphate glucuronosyltransferase. We used population pharmacokinetic modeling to compare the pharmacokinetics of dolutegravir when coadministered with standard- versus high-dose rifampicin in adults with tuberculosis and HIV, and investigated the effect of genetic polymorphisms. Data from the SAEFRIF trial, where participants were randomized to receive first-line tuberculosis treatment with either standard- 10 mg/kg or high-dose 35 mg/kg rifampicin alongside antiretroviral therapy, were used. The dolutegravir model was developed with 211 plasma concentrations from 44 participants. The median (interquartile range) rifampicin area under the curve (AUC) in the standard- and high-dose arms were 32.3 (28.7-36.7) and 153 (138-175) mg·h/L, respectively. A one-compartment model with first-order elimination and absorption through transit compartments best described dolutegravir pharmacokinetics. For a typical 56 kg participant, we estimated a clearance, absorption rate constant, and volume of distribution of 1.87 L/h, 1.42 h-1, and 12.4 L, respectively. Each 10 mg·h/L increase in the AUC of coadministered rifampicin from 32.3 mg·h/L led to a 2.3 (3.1-1.4) % decrease in dolutegravir bioavailability. Genetic polymorphism of UGT1A1 did not significantly affect dolutegravir pharmacokinetics. Simulations of trough dolutegravir concentrations show that the 50 mg twice-daily regimen attains both the primary and secondary therapeutic targets of 0.064 and 0.3 mg/L, respectively, regardless of the dose of coadministered rifampicin, unlike the once-daily regimen.

Population Pharmacodynamic Modelling of the CD19+ Suppression Effects of Rituximab in Paediatric Patients with Neurological and Autoimmune Diseases

BACKGROUND

Limited pharmacotherapy and the failure of conventional treatments in complex pathologies in children lead to increased off-label use of rituximab. We aimed to characterize the time course of CD19+ B lymphocytes (CD19+) under treatment with intravenous rituximab in children with neurologic and autoimmune diseases and to evaluate the impact of covariates (i.e., demographics, diagnosis and substitution between innovator and biosimilar product) on rituximab pharmacodynamics and disease activity.

METHODS

Pre- and post-drug infusion CD19+ in peripheral blood were prospectively registered. A population pharmacodynamic model describing the time course of CD19+ was developed with NONMEM v7.4. Simulations of three different rituximab regimens were performed to assess the impact on CD19+. Logistic regression analysis was performed to identify predictors of clinical response recorded through disease activity scores.

RESULTS

281 measurements of CD19+ lymphocyte counts obtained from 63 children with neurologic (n = 36) and autoimmune (n = 27) diseases were available. The time course of CD19+ was described with a turn-over model in which the balance between synthesis and degradation rates is disrupted by rituximab, increasing the latter process. The model predicts half-lives (percent coefficient of variation, CV(%)) of rituximab and CD19+ of 11.6 days (17%) and 173.3 days (22%), respectively. No statistically significant effect was found between any of the studied covariates and model parameters (p > 0.05). Simulations of different regimens showed no clinically significant differences in terms of CD19+ repopulation times. A trend towards a lack of clinical response was observed in patients with lower CD19+ repopulation times and higher areas under the CD19+ versus time curve.

CONCLUSIONS

Rituximab pharmacodynamics was described in a real-world setting in children suffering from autoimmune and neurologic diseases. Diagnosis, substitution between innovator rituximab and its biosimilars or type of regimen did not affect rituximab-induced depletion of CD19+ nor the clinical response in this cohort of patients. According to this study, rituximab frequency and dosage may be chosen based on clinical convenience or safety reasons without affecting CD19+ repopulation times. Further studies in larger populations are required to confirm these results.

Model-informed precision dosing in vancomycin treatment

Introduction: While vancomycin remains a widely prescribed antibiotic, it can cause ototoxicity and nephrotoxicity, both of which are concentration-associated. Overtreatment can occur when the treatment lasts for an unnecessarily long time. Using a model-informed precision dosing scheme, this study aims to develop a population pharmacokinetic (PK) and pharmacodynamic (PD) model for vancomycin to determine the optimal dosage regimen and treatment duration in order to avoid drug-induced toxicity. Methods: The data were obtained from electronic medical records of 542 patients, including 40 children, and were analyzed using NONMEM software. For PK, vancomycin concentrations were described with a two-compartment model incorporating allometry scaling. Results and discussion: This revealed that systemic clearance decreased with creatinine and blood urea nitrogen levels, history of diabetes and renal diseases, and further decreased in women. On the other hand, the central volume of distribution increased with age. For PD, C-reactive protein (CRP) plasma concentrations were described by transit compartments and were found to decrease with the presence of pneumonia. Simulations demonstrated that, given the model informed optimal doses, peak and trough concentrations as well as the area under the concentration-time curve remained within the therapeutic range, even at doses smaller than routine doses, for most patients. Additionally, CRP levels decreased more rapidly with the higher dose starting from 10 days after treatment initiation. The developed R Shiny application efficiently visualized the time courses of vancomycin and CRP concentrations, indicating its applicability in designing optimal treatment schemes simply based on visual inspection.

Population pharmacokinetics, pharmacodynamics and pharmacogenetics modelling of oxypurinol in Hmong adults with gout and/or hyperuricemia

AIMS

The aim of this study was to quantify identifiable sources of variability, including key pharmacogenetic variants in oxypurinol pharmacokinetics and their pharmacodynamic effect on serum urate (SU).

METHODS

Hmong participants (n = 34) received 100 mg allopurinol twice daily for 7 days followed by 150 mg allopurinol twice daily for 7 days. A sequential population pharmacokinetic pharmacodynamics (PKPD) analysis with non-linear mixed effects modelling was performed. Allopurinol maintenance dose to achieve target SU was simulated based on the final PKPD model.

RESULTS

A one-compartment model with first-order absorption and elimination best described the oxypurinol concentration-time data. Inhibition of SU by oxypurinol was described with a direct inhibitory Emax model using steady-state oxypurinol concentrations. Fat-free body mass, estimated creatinine clearance and SLC22A12 rs505802 genotype (0.32 per T allele, 95% CI 0.13, 0.55) were found to predict differences in oxypurinol clearance. Oxypurinol concentration required to inhibit 50% of xanthine dehydrogenase activity was affected by PDZK1 rs12129861 genotype (-0.27 per A allele, 95% CI -0.38, -0.13). Most individuals with both PDZK1 rs12129861 AA and SLC22A12 rs505802 CC genotypes achieve target SU (with at least 75% success rate) with allopurinol below the maximum dose, regardless of renal function and body mass. In contrast, individuals with both PDZK1 rs12129861 GG and SLC22A12 rs505802 TT genotypes would require more than the maximum dose, thus requiring selection of alternative medications.

CONCLUSIONS

The proposed allopurinol dosing guide uses individuals' fat-free mass, renal function and SLC22A12 rs505802 and PDZK1 rs12129861 genotypes to achieve target SU.

Population pharmacokinetic analysis of lopinavir in HIV negative individuals exposed to SARS-CoV-2: a COPEP (COronavirus Post-Exposure Prophylaxis) sub-study

BACKGROUND

Lopinavir/ritonavir (LPV/r) is a drug traditionally used for the treatment of HIV that has been repurposed as a potential post-exposure prophylaxis agent against COVID-19 in the COronavirus Post-Exposure Prophylaxis (COPEP) study. The present analysis aims to evaluate LPV levels in individuals exposed to SARS-CoV-2 versus people living with HIV (PLWH) by developing a population pharmacokinetic (popPK) model, while characterizing external and patient-related factors that might affect LPV exposure along with dose-response association.

METHODS

We built a popPK model on 105 LPV concentrations measured in 105 HIV-negative COPEP individuals exposed to SARS-CoV-2, complemented with 170 LPV concentrations from 119 PLWH followed in our routine therapeutic drug-monitoring programme. Published LPV popPK models developed in PLWH and in COVID-19 patients were retrieved and validated in our study population by mean prediction error (MPE) and root mean square error (RMSE). The association between LPV model-predicted residual concentrations (Cmin) and the appearance of the COVID-19 infection in the COPEP participants was investigated.

RESULTS

A one-compartment model with linear absorption and elimination best described LPV concentrations in both our analysis and in the majority of the identified studies. Globally, similar PK parameters were found in all PK models, and provided close MPEs (from -19.4% to 8.0%, with a RMSE of 3.4% to 49.5%). No statistically significant association between Cmin and the occurrence of a COVID-19 infection could be detected.

CONCLUSION

Our analysis indicated that LPV circulating concentrations were similar between COPEP participants and PLWH, and that published popPK models described our data in a comparable way.

Population pharmacokinetics of Amisulpride in Chinese patients with schizophrenia with external validation: the impact of renal function

Introduction: Amisulpride is primarily eliminated via the kidneys. Given the clear influence of renal clearance on plasma concentration, we aimed to explicitly examine the impact of renal function on amisulpride pharmacokinetics (PK) via population PK modelling and Monte Carlo simulations. Method: Plasma concentrations from 921 patients (776 in development and 145 in validation) were utilized. Results: Amisulpride PK could be described by a one-compartment model with linear elimination where estimated glomerular filtration rate, eGFR, had a significant influence on clearance. All PK parameters (estimate, RSE%) were precisely estimated: apparent volume of distribution (645 L, 18%), apparent clearance (60.5 L/h, 2%), absorption rate constant (0.106 h-1, 12%) and coefficient of renal function on clearance (0.817, 10%). No other significant covariate was found. The predictive performance of the model was externally validated. Covariate analysis showed an inverse relationship between eGFR and exposure, where subjects with eGFR= 30 mL/min/1.73 m2 had more than 2-fold increase in AUC, trough and peak concentration. Simulation results further illustrated that, given a dose of 800 mg, plasma concentrations of all patients with renal impairment would exceed 640 ng/mL. Discussion: Our work demonstrated the importance of renal function in amisulpride dose adjustment and provided a quantitative framework to guide individualized dosing for Chinese patients with schizophrenia.

Population Pharmacokinetic Analysis of Perampanel in Portuguese Patients Diagnosed with Refractory Epilepsy

Perampanel is a promising antiepileptic drug (AED) for refractory epilepsy treatment due to its innovative mechanism of action. This study aimed to develop a population pharmacokinetic (PopPK) model to be further used in initial dose optimization of perampanel in patients diagnosed with refractory epilepsy. A total of seventy-two plasma concentrations of perampanel obtained from forty-four patients were analyzed through a population pharmacokinetic approach by means of nonlinear mixed effects modeling (NONMEM). A one-compartment model with first-order elimination best described the pharmacokinetic profiles of perampanel. Interpatient variability (IPV) was entered on clearance (CL), while the residual error (RE) was modeled as proportional. The presence of enzyme-inducing AEDs (EIAEDs) and body mass index (BMI) were found as significant covariates for CL and volume of distribution (V), respectively. The mean (relative standard error) estimates for CL and V of the final model were 0.419 L/h (5.56%) and 29.50 (6.41%), respectively. IPV was 30.84% and the proportional RE was 6.44%. Internal validation demonstrated an acceptable predictive performance of the final model. A reliable population pharmacokinetic model was successfully developed, and it is the first enrolling real-life adults diagnosed with refractory epilepsy.

External validation of population pharmacokinetic models of gentamicin in paediatric population from preterm newborns to adolescents

The aim of this study was to externally validate the predictive performance of published population pharmacokinetic models of gentamicin in all paediatric age groups, from preterm newborns to adolescents. We first selected published population pharmacokinetic models of gentamicin developed in the paediatric population with a wide age range. The parameters of the literature models were then re-estimated using the PRIOR subroutine in NONMEM^®. The predictive ability of the literature and the tweaked models was evaluated. Retrospectively collected data from a routine clinical practice (512 concentrations from 308 patients) were used for validation. The models with covariates characterising developmental changes in clearance and volume of distribution had better predictive performance, which improved further after re-estimation. The tweaked model by Wang 2019 performed best, with suitable accuracy and precision across the complete paediatric population. For patients treated in the intensive care unit, a lower proportion of patients would be expected to reach the target trough concentration at standard dosing. The selected model could be used for model-informed precision dosing in clinical settings where the entire paediatric population is treated. However, for use in clinical practice, the next step should include additional analysis of the impact of intensive care treatment on gentamicin pharmacokinetics, followed by prospective validation.

Donepezil Brain and Blood Pharmacokinetic Modeling after Nasal Film and Oral Solution Administration in Mice

Intranasal delivery is a non-invasive mode of administration, gaining popularity due to its potential for targeted delivery to the brain. The anatomic connection of the nasal cavity with the central nervous system (CNS) is based on two nerves: olfactory and trigeminal. Moreover, the high vasculature of the respiratory area enables systemic absorption avoiding possible hepatic metabolism. Due to these physiological peculiarities of the nasal cavity, compartmental modeling for nasal formulation is considered a demanding process. For this purpose, intravenous models have been proposed, based on the fast absorption from the olfactory nerve. However, most of the sophisticated approaches are required to describe the different absorption events occurring in the nasal cavity. Donepezil was recently formulated in the form of nasal film ensuring drug delivery in both bloodstream and the brain. In this work, a three-compartment model was first developed to describe donepezil oral brain and blood pharmacokinetics. Subsequently, using parameters estimated by this model, an intranasal model was developed dividing the administered dose into three fractions, corresponding to absorption directly to the bloodstream and brain, as well as indirectly to the brain expressed through transit compartments. Hence, the models of this study aim to describe the drug flow on both occasions and quantify the direct nose-to-brain and systemic distribution.

A model-based approach for a practical dosing strategy for the short, intensive treatment regimen for paediatric tuberculous meningitis

Following infection with Mycobacterium tuberculosis, young children are at high risk of developing severe forms of tuberculosis (TB) disease, including TB meningitis (TBM), which is associated with significant morbidity and mortality. In 2022, the World Health Organization (WHO) conditionally recommended that a 6-month treatment regimen composed of higher doses of isoniazid (H) and rifampicin (R), with pyrazinamide (Z) and ethionamide (Eto) (6HRZEto), be used as an alternative to the standard 12-month regimen (2HRZ-Ethambutol/10HR) in children and adolescents with bacteriologically confirmed or clinically diagnosed TBM. This regimen has been used in South Africa since 1985, in a complex dosing scheme across weight bands using fixed-dose combinations (FDC) available locally at the time. This paper describes the methodology used to develop a new dosing strategy to facilitate implementation of the short TBM regimen based on newer globally available drug formulations. Several dosing options were simulated in a virtual representative population of children using population PK modelling. The exposure target was in line with the TBM regimen implemented in South Africa. The results were presented to a WHO convened expert meeting. Given the difficulty to achieve simple dosing using the globally available RH 75/50 mg FDC, the panel expressed the preference to target a slightly higher rifampicin exposure while keeping isoniazid exposures in line with those used in South Africa. This work informed the WHO operational handbook on the management of TB in children and adolescents, in which dosing strategies for children with TBM using the short TBM treatment regimen are provided.

Population pharmacokinetics of anlotinib, a multiple receptor tyrosine kinase inhibitor, in chinese patients with malignant tumors

The objective of this study was to investigate the pharmacokinetic behavior of anlotinib in Chinese patients with malignant tumors using the population approach. A total of 407 anlotinib plasma concentrations from 16 patients were analyzed in this study. Anlotinib was administered orally 12 or 16 mg in the single-dose phase and 12 mg once daily in the multiple-dose phase. A population pharmacokinetic (PopPK) model was established using nonlinear mixed-effects model (NONMEM) method. The potential influence of demographic and pathophysiological factors on oral anlotinib pharmacokinetics was investigated in a covariate analysis. The final model was evaluated using goodness-of-fit plots, visual predictive check, and bootstrap methods. The pharmacokinetic profile of anlotinib was best described by a one-compartment model with first-order absorption and first-order elimination. The population estimates of the apparent total clearance (CL/F), apparent volume of distribution (V/F) and absorption rate constant (Ka) were 8.91 L/h, 1950 L and 0.745 h^-1 , respectively. Body weight was identified as a significant covariate on V/F. Patients with low body weight tended to show higher exposure to anlotinib than those with high body weight. However, these differences were not clinically significant based on the simulations of the individual body weight effects. Taken together, this popPK model adequately described the pharmacokinetics of anlotinib in patients with malignant tumors and supports the same starting dose among them. This article is protected by copyright. All rights reserved.

Evaluation of treosulfan cumulative exposure in paediatric patients through population pharmacokinetics and dosing simulations

AIM

To investigate the pharmacokinetics (PK) of intravenous treosulfan in paediatric patients undergoing haematopoietic stem cell transplantation (HSCT) for a broad range of diseases and to explore the impact of different dosing regimens on treosulfan exposure (area under the concentration-time curve, AUC_0→∞ ) through dosing simulations.

METHODS

A prospective multicentre PK study was conducted using treosulfan concentration data (n = 423) collected from 53 children (median age 3.5, range 0.2-17.0 years) receiving three daily age-guided doses (10-14 g/m² ). Population PK modelling was performed using NONMEM software, utilising a stepwise forward selection backward elimination method and likelihood-ratio test for screening covariates to describe PK variability. Monte Carlo simulation was used to generate patient PK data for 10 000 virtual paediatric patients and cumulative AUC_0→∞ values were evaluated using age, body surface area (BSA) and model-based dosing regimens, targeting 4800 mg*h/L.

RESULTS

Treosulfan concentration data were described using a one-compartment PK model with first-order elimination. Population mean (95% CI) estimates for clearance (CL) and volume of distribution (V) were 16.3 (14.9-18.1) L/h and 41.9 (38.8-45.1) L, respectively. Allometrically scaled body weight was the best covariate descriptor for CL and V, and maturational age further explained variability in CL. Dosing simulations indicated that in young patient groups (<2 years), a model-based dosing regimen more accurately achieved the target AUC_0→∞ (58.3%) over the age (42.6%) and BSA-based (51.3%) regimens.

CONCLUSION

Treosulfan disposition was described through allometric body weight and maturational age descriptors. Model-informed dosing is recommended for patients under 2 years. Treosulfan PK parameters and AUC_0→∞ were not influenced by patient disease.

Population Pharmacokinetic Study of Benzylpenicillin in Critically Unwell Adults

Pharmacokinetics are highly variable in critical illness, and suboptimal antibiotic exposure is associated with treatment failure. Benzylpenicillin is a commonly used beta-lactam antibiotic, and pharmacokinetic data of its use in critically ill adults are lacking. We performed a pharmacokinetic study of critically unwell patients receiving benzylpenicillin, using data from the ABDose study. Population pharmacokinetic modelling was undertaken using NONMEM version 7.5, and simulations using the final model were undertaken to optimize the pharmacokinetic profile. We included 77 samples from 12 participants. A two-compartment structural model provided the best fit, with allometric weight scaling for all parameters and a creatinine covariate effect on clearance. Simulations (n = 10,000) demonstrated that 25% of simulated patients receiving 2.4 g 4-hourly failed to achieve a conservative target of 50% of the dosing interval with free drug above the clinical breakpoint MIC (2 mg/L). Simulations demonstrated that target attainment was improved with continuous or extended dosing. To our knowledge, this study represents the first full population PK analysis of benzylpenicillin in critically ill adults.

Population pharmacokinetic/pharmacodynamic evaluations of amikacin dosing in critically ill patients undergoing continuous venovenous hemodiafiltration

OBJECTIVES

The pharmacokinetics/pharmacodynamics (PK/PD) of amikacin in critically ill patients undergoing continuous venovenous hemodiafiltration (CVVHDF) are poorly described, and appropriate dosing is unclear in this patient population. This study aimed to develop a population PK model of amikacin and to provide systemic PK/PD evaluations for different dosing regimens in CVVHDF patients.

METHODS

One hundred and sixty-one amikacin concentration observations from thirty-three CVVHDF patients were pooled to develop the population PK model. Monte Carlo simulations were performed to assess the PK/PD index-based efficacy (Cmax/minimal inhibitory concentration (MIC) > 8 and AUC/MIC > 58.3), nonrisk of drug resistance (T>MIC > 60%) and risk of toxicity (trough concentration > 5 mg/l) for different dosing regimens.

KEY FINDINGS

A two-compartment model adequately described the concentration data of amikacin. A loading dose of at least 25 mg/kg amikacin is needed to reach the efficacy targets in CVVHDF patients for an MIC of 4 mg/l, and the studied doses could not provide adequate drug exposure and T>MIC > 60% for an MIC ≥ 8 mg/l. The risk of toxicity for amikacin was unacceptably high for the patient population with low clearance.

CONCLUSIONS

Our study demonstrated that a loading dose of 25-30 mg/kg amikacin is needed to provide adequate PK/PD target attainment in CVVHDF patients for an MIC ≤ 4 mg/l.

Population pharmacokinetics and pharmacodynamics of cefazolin using total and unbound serum concentrations in patients with high body weight

The objective of this study was to evaluate the steady state pharmacokinetics and pharmacodynamics of cefazolin in patients with a high body weight. Cefazolin was administered by 0.5-h infusions to 11 patients with total body weight (TBW) ≥120 kg receiving 3 g q8h, and 12 patients with TBW <120 kg receiving 2 g q8h. Total and unbound serum concentration-time data obtained from serial blood samples were analysed simultaneously by population pharmacokinetic modelling using NONMEM. Probability of target attainment (PTA) was calculated for various dosing regimens through Monte Carlo simulations based on the cumulative percentage of the dosing interval that the unbound concentration exceeds the minimum inhibitory concentration (MIC) value for the pathogen at steady state (fT_MIC) ≥40%, ≥60% and 100%. A two-compartment model with non-linear protein binding and allometric scaling of the central volume of distribution using TBW best characterized both total and unbound concentration-time data. Unbound clearance was significantly associated with creatinine clearance, and maximum protein binding constant was significantly associated with serum albumin concentration and body mass index (P <0.05). Based on unbound concentration-time profiles, all simulated regimens achieved PTA >90% at MIC values ≤2 mg/L using fT_MIC ≥40%, at MIC values ≤1 mg/L using fT_MIC ≥60%, and at MIC values ≤0.5 mg/L using fT_MIC of 100%. At fT_MIC ≥60%, 0.5-h infusion of cefazolin 1 g q8h achieved PTA <90% at MIC values ≥2 mg/L in patients with TBW≥120 kg; however, prolonged-infusion and higher-dose regimens improved PTA to >90%. Overall, cefazolin pharmacokinetics are altered considerably in obese patients. Higher-dose and/or prolonged-infusion cefazolin regimens should be considered in patients with TBW ≥120 kg, particularly those with less-susceptible Gram-negative infections.

Prediction of the pharmacokinetics of pemetrexed with a low test dose: A proof-of-concept study

PURPOSE

Pemetrexed is a cytotoxic drug used for the treatment of lung cancer and mesothelioma. The use of a low test dosing of cytotoxic drugs may aid in dose individualization without causing harm. The aim of this proof-of-concept study was to assess if the pharmacokinetics (PKs) of a test dose could predict the PKs of a therapeutic pemetrexed dose.

METHODS

Ten patients received both a low test dose (100 μg) and a therapeutic dose of pemetrexed after which plasma concentrations pemetrexed were measured. PK analysis was performed by means of nonlinear mixed-effects modelling. The predictive performances of test dose clearance and renal function towards a therapeutic dose were assessed.

RESULTS

The PKs of a pemetrexed test dose were best described by a one-compartment model with linear elimination. A high variability in the administered dose was observed for the test dose, but not for the therapeutic dose. A statistically significant correlation between test dose clearance and therapeutic dose clearance was observed (Spearman's rho: 0.758, P = 0.02). The predictive performance of test dose clearance was worse than renal function: mean predictive error (+95% confidence interval [CI]) 53.9% (50.1-57.6%) vs 19.4% (12.4-26.4%) and normalized root-mean square error (+95% CI) 57.8% (30.5-85.1%) vs 25.7% (20.3-31.0%).

CONCLUSION

We show that test dosing of pemetrexed is feasible, but there seems no added value for a low test dosing in the dose individualization of pemetrexed.

Pharmacokinetic comparability between two populations using nonlinear mixed effect models: a Monte Carlo study

'Are two populations the same or are they different' is a question that is often faced in clinical pharmacology trials e.g., a pharmacokinetic trial studying a particular drug in racially different groups. To address this question, concentration-time data were simulated from a reference and test population, where in the latter the clearance, sample size, and sampling design were systematically varied. It was of interest to determine whether the estimates of clearance from the two groups were the same or different. Two approaches were used to estimate the empirical Bayes estimates (EBEs) for clearance. One approach developed a population pharmacokinetic model for the reference population and the EBEs for the reference population were estimated from this model. The parameters of the reference population were fixed to their maximum likelihood estimates. The model was then applied to the test population dataset to estimate the EBEs of the test population using the MAXEVAL = 0 option in NONMEM. A second approach, the theta approach, combined the reference and test datasets into a single dataset and used population as a covariate in the model; the EBEs were estimated from this combined model. The power and type I error rate of each approach were calculated for each treatment combination using a variety of statistical tests to determine whether there was a difference in the distribution of the EBEs in the reference population compared to the test population. Our results suggest that either MAXEVAL or theta approaches can be used with informative sampling designs. In addition to reasonable power and type I error, both approaches gave almost identical results under a dense sampling design. To statistically compare the distribution of EBEs of pharmacokinetic parameters from a reference group to that of a test group, a T-test and DTS eCDF test are equally useful.

Latent variable indirect response modeling of clinical efficacy endpoints with combination therapy: application to guselkumab and golimumab in patients with ulcerative colitis

Accurate characterization of longitudinal exposure-response of clinical trial endpoints is important in optimizing dose and dosing regimens in drug development. Clinical endpoints are often categorical, for which much progress has been made recently in latent variable indirect response (IDR) modeling with single drugs. However, such applications have not yet been used for trials employing multiple drugs administered concurrently. This study aims to demonstrate that the latent variable IDR approach provides a convenient longitudinal exposure-response modeling framework to assess potential interaction effects of combination therapies. This is illustrated by an application to the exposure-response modeling of guselkumab, a monoclonal antibody in clinical development that blocks the interleukin-23p19 subunit, and golimumab, a monoclonal antibody that binds with high affinity to tumor necrosis factor-alpha. A Phase 2a study was conducted in 214 patients with moderate-to severe active ulcerative colitis for which longitudinal assessments of disease severity based on patient-reported measures of rectal bleeding, stool frequency, and symptomatic remission were evaluated as categorical endpoints, and fecal calprotectin as a continuous endpoint. The modeling results suggested independent pharmacodynamic guselkumab and golimumab effects on fecal calprotectin as a continuous endpoint, as well as interaction effects on the categorical endpoints that may be explained by an additional pathway of competitive interaction.

Influence of Renal Function on Phosphoramide Mustard Exposure: A Nonlinear Mixed-Effects Analysis

Phosphoramide mustard (PM) is the final cytotoxic metabolite formed from the parent compound cyclophosphamide through a complex metabolic pathway, primarily through hepatic metabolism. Little is known about the effect of renal elimination on the disposition of PM. We evaluated the effect of renal function on PM exposure after single doses of cyclophosphamide in 85 patients undergoing allogeneic hematopoietic cell transplantation using nonlinear mixed-effects modeling. Mixed linear and nonlinear elimination pathways were required to adequately describe the disposition of PM. Creatinine clearance (CrCL) was incorporated as a covariate associated with first-order elimination, representing renal clearance (Cl_R ) of PM. For a 70-kg patient, Cl_R was 14.9 L/h, Volume of distribution was 525 L, maximum rate was 81.2 mg/h, and the concentration to achieve 50% of maximum rate was 0.51 mg/L. We conducted simulations to explore the impact of CrCL as a measure of renal function and observed that when CrCL decreases from 120 to 40 mL/min, PM area under the plasma concentration-time curve (AUC) from time 0 to 8 hours and AUC increases by 9.2% and 80.9% on average after a single dose, respectively. Our data suggest that renal function has limited influence on PM exposure during the first 8 hours after dosing but has a large impact on the total exposure. Dose adjustment of cyclophosphamide may not be necessary in hematopoietic cell transplant recipients with moderate to severe kidney dysfunction to attain targeted exposures based on AUC from time 0 to 8 hours. However, dose reduction may be necessary if demonstrated at some future time that total AUC is a better surrogate for safety or toxicity.

Pharmacokinetic and Pharmacodynamic Modeling and Simulation Analysis of Prasugrel in Healthy Male Volunteers

This study evaluated the pharmacokinetics and pharmacodynamics of the antiplatelet agent prasugrel, and explored its optimal dose regimens via modeling and simulation using NONMEM. We measured platelet aggregation and the serial plasma concentrations of the inactive (R-95913) and active metabolites (R-138727) of prasugrel after a single oral dose of 10-60 mg in 20 healthy adult male volunteers. A pharmacokinetic model for R-95913 and R-138727, and a pharmacodynamic model between the concentration of R-138727 and maximal platelet aggregation measured by light transmittance were constructed. The predictability of the model for platelet aggregation was evaluated by comparing the model prediction values with the observed ones not used in the construction of the model. Pharmacokinetic data were best described by a 3-compartment models for R-95913, a 1-compartment model for R-138727 with transit compartment model for absorption delay, and first-pass metabolic conversion of R-95913 into R-138727 during absorption. The association-dissociation model between R-138727 and its receptor in platelets was applied for the inhibitory effect of prasugrel on platelet aggregation. Prasugrel rapidly inhibited platelet aggregation after oral administration, with a prolonged duration of action; however, the concentration of the active metabolite decreased rapidly, which may have been due to the slow dissociation rate of R-138727 from its target receptor in platelets. The external validation suggests that our model could be used to individualize prasugrel treatment in various clinical situations. Simulation showed rapid onset of inhibitory effect with great magnitude and consistent inhibition after therapeutic dose of prasugrel.

Population pharmacokinetics of nalbuphine in patients undergoing general anesthesia surgery

Purpose: The aim of this study was to build a population pharmacokinetics (PopPK) model of nalbuphine and to estimate the suitability of bodyweight or fixed dosage regimen. Method: Adult patients who were undergoing general anesthetic surgery using nalbuphine for induction of anesthesia were included. Plasma concentrations and covariates information were analyzed by non-linear mixed-effects modeling approach. Goodness-of-fit (GOF), non-parametric bootstrap, visual predictive check (VPC) and external evaluation were applied for the final PopPK model evaluation. Monte Carlo simulation was conducted to assess impact of covariates and dosage regimens on the plasma concentration to nalbuphine. Results: 47 patients aged 21-78 years with a body weight of 48-86 kg were included in the study. Among them, liver resection accounted for 14.8%, cholecystectomy for 12.8%, pancreatic resection for 36.2% and other surgeries for 36.2%. 353 samples from 27 patients were enrolled in model building group; 100 samples from 20 patients were enrolled in external validation group. The results of model evaluation showed that the pharmacokinetics of nalbuphine was adequately described by a two-compartment model. The hourly net fluid volume infused (HNF) was identified as a significant covariate about the intercompartmental clearance (Q) of nalbuphine with objective function value (OFV) decreasing by 9.643 (p < 0.005, df = 1). Simulation results demonstrated no need to adjust dosage based on HNF, and the biases of two dosage methods were less than 6%. The fixed dosage regimen had lower PK variability than the bodyweight regimen. Conclusion: A two-compartment PopPK model adequately described the concentration profile of nalbuphine intravenous injection for anesthesia induction. While HNF can affect the Q of nalbuphine, the magnitude of the effect was limited. Dosage adjustment based on HNF was not recommended. Furthermore, fixed dosage regimen might be better than body weight dosage regimen.

Predictive model of recurrent ischemic stroke: model development from real-world data

Background

There are established correlations between risk factors and ischemic stroke (IS) recurrence; however, does the hazard of recurrent IS change over time? What is the predicted baseline hazard of recurrent IS if there is no influence of variable predictors? This study aimed to quantify the hazard of recurrent IS when the variable predictors were set to zero and quantify the secondary prevention influence on the hazard of recurrent ischemic stroke.

Methods

In the population cohort involved in this study, data were extracted from 7,697 patients with a history of first IS attack registered with the National Neurology Registry of Malaysia from 2009 to 2016. A time-to-recurrent IS model was developed using NONMEM version 7.5. Three baseline hazard models were fitted into the data. The best model was selected using maximum likelihood estimation, clinical plausibility, and visual predictive checks.

Results

Within the maximum 7.37 years of follow-up, 333 (4.32%) patients had at least one incident of recurrent IS. The data were well described by the Gompertz hazard model. Within the first 6 months after the index IS, the hazard of recurrent IS was predicted to be 0.238, and 6 months after the index attack, it reduced to 0.001. The presence of typical risk factors such as hyperlipidemia [HR, 2.22 (95%CI: 1.81-2.72)], hypertension [HR, 2.03 (95%CI: 1.52-2.71)], and ischemic heart disease [HR, 2.10 (95%CI: 1.64-2.69)] accelerated the hazard of recurrent IS, but receiving antiplatelets (APLTs) upon stroke decreased this hazard [HR, 0.59 (95%CI: 0.79-0.44)].

Conclusion

The hazard of recurrent IS magnitude differs during different time intervals based on the concomitant risk factors and secondary prevention.

Predicting treatment response to vancomycin using bacterial DNA load as a pharmacodynamic marker in premature and very low birth weight neonates: A population PKPD study

Background: While positive blood cultures are the gold standard for late-onset sepsis (LOS) diagnosis in premature and very low birth weight (VLBW) newborns, these results can take days, and early markers of possible treatment efficacy are lacking. The objective of the present study was to investigate whether the response to vancomycin could be quantified using bacterial DNA loads (BDLs) determined by real-time quantitative polymerase chain reaction (RT-qPCR). Methods: VLBW and premature neonates with suspected LOS were included in a prospective observational study. Serial blood samples were collected to measure BDL and vancomycin concentrations. BDLs were measured with RT-qPCR, whereas vancomycin concentrations were measured by LC-MS/MS. Population pharmacokinetic-pharmacodynamic modeling was performed with NONMEM. Results: Twenty-eight patients with LOS treated with vancomycin were included. A one-compartment model with post-menstrual age (PMA) and weight as covariates was used to describe the time PK profile of vancomycin concentrations. In 16 of these patients, time profiles of BDL could be described with a pharmacodynamic turnover model. The relationship between vancomycin concentration and first-order BDL elimination was described with a linear-effect model. Slope S increased with increasing PMA. In 12 patients, no decrease in BDL over time was observed, which corresponded with clinical non-response. Discussion: BDLs determined through RT-qPCR were adequately described with the developed population PKPD model, and treatment response to vancomycin using BDL in LOS can be assessed as early as 8 h after treatment initiation.

Optimal dosing of enoxaparin in overweight and obese children

AIM

Current enoxaparin dosing guidelines in children are based on total body weight. This is potentially inappropriate in obese children as it may overestimate the drug clearance. Current evidence suggests that obese children may require lower initial doses of enoxaparin, therefore the aim of this work was to characterise the pharmacokinetics of enoxaparin in obese children and to propose a more appropriate dosing regimen.

METHODS

Data from 196 unique encounters of 160 children who received enoxaparin treatment doses were analysed. Enoxaparin concentration was quantified using the chromogenic anti factor Xa (anti-Xa) assay. Patients provided a total of 552 anti-Xa samples. Existing published pharmacokinetic (PK) models were fitted and evaluated against our dataset using prediction-corrected visual predictive check plots (pcVPCs). A PK model was fitted using a nonlinear mixed-effects modelling approach. The fitted model was used to evaluate the current standard dosing and identify an optimal dosing regimen for obese children.

RESULTS

Published models of enoxaparin pharmacokinetics in children did not capture the pharmacokinetics of enoxaparin in obese children as shown by pcVPCs. A one-compartment model with linear elimination best described the pharmacokinetics of enoxaparin. Allometrically scaled fat-free mass with an estimated exponent of 0.712 (CI 0.66-0.76) was the most influential covariate on clearance while linear fat-free mass was selected as the covariate on volume. Simulations from the model showed that fat-free mass-based dosing could achieve the target anti-Xa activity at steady state in 77.5% and 78.2% of obese and normal-weight children, respectively, compared to 65.2% and 75.5% for standard total body weight-based dosing.

CONCLUSIONS

A population PK model that describes the time course of anti-Xa activity of enoxaparin was developed in a paediatric population. Based on this model, a unified dosing regimen was proposed that will potentially improve the success rate of target attainment in overweight/obese patients without the need for patient body size categorisation. Therefore, prospective validation of the proposed approach is warranted.

Pharmacogenetics-based population pharmacokinetic analysis and dose optimization of valproic acid in Chinese southern children with epilepsy: Effect of ABCB1 gene polymorphism

Objective: The aim of this study was to establish a population pharmacokinetic (PPK) model of valproic acid (VPA) in pediatric patients with epilepsy in southern China, and provide guidance for individualized medication of VPA therapy. Methods: A total of 376 VPA steady-state trough concentrations were collected from 103 epileptic pediatric patients. The PPK parameter values for VPA were calculated by using the nonlinear mixed-effects modeling (NONMEM) method, and a one-compartment model with first-order absorption and elimination processes was applied. Covariates included demographic information, concomitant medications and selected gene polymorphisms. Goodness-of-fit (GOF), bootstrap analysis, and visual predictive check (VPC) were used for model evaluation. In addition, we used Monte Carlo simulations to propose dose recommendations for different subgroup patients. Results: A significant effect of the patient age and ABCB1 genotypes was observed on the VPA oral clearance (CL/F) in the final PPK model. Compared with patients with the ABCB1 rs3789243 AA genotype, CL/F in patients with GG and AG genotypes was increased by 8% and reduced by 4.7%, respectively. The GOF plots indicated the satisfactory predictive performance of the final model, and the evaluation by bootstrap and VPC showed that a stable model had been developed. A table of individualized dosing regimens involving age and ABCB1 genotype was constructed based on the final PPK model. Conclusion: This study quantitatively investigated the effects of patient age and ABCB1 rs3789243 variants on the pharmacokinetic variability of VPA. The PPK models could be beneficial to individual dose optimization in epileptic children on VPA therapy.

Population Pharmacokinetic Meta-Analysis and Dosing Recommendation for Meropenem in Critically Ill Patients Receiving Continuous Renal Replacement Therapy

The optimal dosing regimen for meropenem in critically ill patients undergoing continuous renal replacement therapy (CRRT) remains undefined due to small studied sample sizes and uninformative pharmacokinetic (PK)/pharmacodynamic (PD) analyses in reported studies. The present study aimed to perform a population PK/PD meta-analysis of meropenem using available literature data to suggest the optimal treatment regimen. A total of 501 meropenem concentration measurements from 78 adult CRRT patients pooled from nine published studies were used to develop the population PK model for meropenem. PK/PD target (40% and 100% of the time with the unbound drug plasma concentration above the MIC) marker-based efficacy and risk of toxicity (trough concentrations of >45 mg/L) for short-term (30 min), prolonged (3 h), and continuous (24 h) infusion dosing strategies for meropenem were investigated. The impact of CRRT dose and identified covariates on the PD probability of target attainment (PTA) and predicted toxicity was also examined. Meropenem concentration data were adequately described by a two-compartment model with linear elimination. Trauma was identified as a pronounced modifier for endogenous clearance of meropenem. Simulations demonstrated that adequate PK/PD targets and low risk of toxicity could be achieved in non-trauma CRRT patients receiving meropenem regimens of 1 g every 6 h infused over 30 min, 1 g every 8 h infused over 3 h, and 2 to 4 g every 24 h infused over 24 h. The impact of CRRT dose (25 to 50 mL/kg/h) on PTA was clinically irrelevant, and continuous infusion of 3 to 4 g every 24 h was suitable for trauma CRRT patients (MICs of ≤0.5 mg/L). A population PK model was developed for meropenem in CRRT patients, and different dosing regimens were proposed for non-trauma and trauma CRRT patients.

Population Pharmacokinetics of Z-Endoxifen in Patients With Advanced Solid Tumors

The purpose of this study was to develop and validate a population pharmacokinetic model for Z-endoxifen in patients with advanced solid tumors and to identify clinical variables that influence pharmacokinetic parameters. Z-endoxifen-HCl was administered orally once a day on a 28-day cycle (±3 days) over 11 dose levels ranging from 20 to 360 mg. A total of 1256 Z-endoxifen plasma concentration samples from 80 patients were analyzed using nonlinear mixed-effects modeling to develop a population pharmacokinetic model for Z-endoxifen. A 2-compartment model with oral depot and linear elimination adequately described the data. The estimated apparent total clearance, apparent central volume of distribution, and apparent peripheral volume of distribution were 4.89 L/h, 323 L, and 39.7 L, respectively, with weight-effect exponents of 0.75, 1, and 1, respectively. This model was used to explore the effects of clinical and demographic variables on Z-endoxifen pharmacokinetics. Weight, race on clearance, and aspartate aminotransferase on the absorption rate constant were identified as significant covariates in the final model. This novel population pharmacokinetic model provides insight regarding factors that may affect the pharmacokinetics of Z-endoxifen and may assist in the design of future clinical trials.

Population Pharmacodynamic Analyses of Human Anti-Rabies Virus Monoclonal Antibody (Ormutivimab) in Healthy Adult Subjects

Ormutivimab is the first recombinant human anti-rabies monoclonal antibody (rhRIG) approved for clinical application in China. In this study, a population pharmacodynamic (PPD) model was established to compare the neutralizing antibody activities of Ormutivimab and human rabies immunoglobulin (HRIG), alone or combined with human rabies vaccine (Vero), in a phase II clinical trial, and to recommend a target dose for the phase III trial. The model was verified to fit the PPD data well. The stability of the model was verified by the bootstrap method. The level of neutralizing antibodies in vivo increased rapidly after administration of Ormutivimab or HRIG. Neutralizing antibodies with a strong activity were produced at 7 days (Ormutivimab + vaccine) or 10 days (HRIG + vaccine) after induction by the vaccine in vivo. Compared to that induced by HRIG + vaccine, the level of the neutralizing antibodies induced by Ormutivimab + vaccine peaked higher and faster. The levels of neutralizing antibodies induced by Ormutivimab + vaccine and HRIG + vaccine were similar within 21 days after administration. According to these results and the safety data, 20 IU·kg^-1 was recommended as the target dose in the confirmatory study of Ormutivimab. Registration: ClinicalTrials.gov #NCT02559921.

Population pharmacokinetic modelling of imatinib in healthy subjects receiving a single dose of 400 mg

Purpose

Imatinib is indicated for treatment of CML, GIST, etc. The population pharmacokinetics (popPK) of imatinib in patients under long-term treatment are reported in literature. Data obtained from bioequivalence trials for healthy subjects were used to evaluate the influence of demographic and pharmacogenetic factors on imatinib pharmacokinetics (PK) in a collective without concurrent drugs, organ dysfunction, inflammation etc. In addition, the differences in PK between the healthy subjects and a patient cohort was examined to identify possible disease effects.

Methods

26 volunteers were administered orally with single dose of 400 mg imatinib. 16–19 plasma samples per volunteer were collected from 0.5 up to 72 h post-dose. The popPK was built and post hoc estimates were compared with previously published PK parameters evaluated by non-compartmental analysis in the same cohort. The predictivity of the model for data collected from 40 patients with gastrointestinal stromal tumors at steady state was evaluated.

Results

The popPK was best described by a two-compartment transit model with first-order elimination. No significant covariates were identified, probably due to the small cohort and the narrow range of demographic covariates; CYP3A5 phenotypes appeared to have some influence on the clearance of imatinib. Good agreement between non-compartment and popPK analyses was observed with the differences of the geometric means/ median of PK estimates below 10%. The model indicated lower clearance for patients compared to healthy volunteers (p value < 0.01).

Conclusion

The two-compartment transit model adequately describes the absorption and distribution of imatinib in healthy volunteers. For patients, a lower clearance of imatinib compared to healthy volunteer was estimated by the model. The model can be applied for dose individualization based on trough concentrations assuming no significant differences in absorption between patients and healthy volunteers.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00280-022-04454-y.

Improving categorical endpoint longitudinal exposure-response modeling through the joint modeling with a related endpoint

Exposure-response modeling is important to optimize dose and dosing regimens in clinical drug development. While primary clinical trial endpoints often have few categories and thus provide only limited information, sometimes there may be additional, more informative endpoints. Benefits of fully incorporating relevant information in longitudinal exposure-response modeling through joint modeling have recently been shown. This manuscript aims to further investigate the benefit of joint modeling of an ordered categorical primary endpoint with a related near-continuous endpoint, through the sharing of model parameters in the latent variable indirect response (IDR) modeling framework. This is illustrated by analyzing the data collected through up to 116 weeks from a phase 3b response-adaptive trial of ustekinumab in patients with psoriasis. The primary endpoint was based on the 6-point physician's global assessment (PGA) score. The Psoriasis area and severity Index (PASI) data, ranging from 0 to 72 with 0.1 increments, were also available. Separate and joint latent variable Type I IDR models of PGA and PASI scores were developed and compared. The results showed that the separate PGA model had a substantial structural bias, which was corrected by the joint modeling of PGA and PASI scores.

Population pharmacokinetics and IVIVC for mesalazine enteric-coated tablets

Although in-vivo bioequivalence (BE) study serves as a golden standard for establishing interchangeability of oral dosage forms, it remains challenging for products with high inter-subject variability such as mesalazine enteric-coated tablet to fulfil the BE criteria set by regulatory authorities. Mesalazine, as a BCS class IV drug, targets to be delivered to distal ileum or colon with a pH-sensitive polymer coating for the remission of ulcerative colitis. Through population pharmacokinetic (PK) analysis and in-vitro in-vivo correlation (IVIVC) modeling on the dissolution and BE data of a generic enteric-coated product (EM) and its reference Salofalk® 250 mg tablet (SM), we for the first time revealed the underlying mechanism of the high inter-subject variability for such delayed-release formulation. It was also noted that the in-vivo start time of absorption (T_s) for EM and SM was positively correlated with their in-vitro lag time (T_lag) under the USP three-stage dissolution condition and reversely correlated with their in-vivo bioavailability. The varied oral bioavailability of mesalazine enteric-coated tablet was mainly due to the varied N-acetyltransferase activities along GI tract. Although such extensive intestinal first-pass metabolism with large individual differences led to a significant variation of mesalazine C_max (coefficient of variation: 60-63.5%) and AUC_0-t (coefficient of variation: 37.5-46.9%), the corresponding variations in the total absorbed mesalazine (mesalazine and its metabolite N-acetyl mesalazine) were significantly reduced by 12 to 45%. Since the BE purpose for mesalazine enteric-coated tablet focused on their comparable safety profiles, total absorbed mesalazine was recommended to be adopted for the development of the IVIVC model and BE evaluation for EM. All in all, our model-based approach has not only successfully identified the key factors that affect the BE of EM to guide its further formulation optimization, but also demonstrated the indispensable role of modeling in the development of generic pharmaceutical product at its early stages.

Model-Based Anticancer Effect of Botulinum Neurotoxin Type A1 on Syngeneic Melanoma Mice

In recent, Botulinum Neurotoxin A1 (BoNT/A1) has been suggested as a potential anticancer agent due to neuronal innervation in tumor cells. Although potential BoNT/A1's mechanism of action for the tumor suppression has been gradually revealed so far, there were no reports to figure out the exposure-response relationships because of the difficulty of its quantitation in the biological matrix. The main objectives of this study were to measure the anticancer effect of BoNT/A1 using a syngeneic mouse model transplanted with melanoma cells (B16-F10) and developed a kinetic-pharmacodynamic (K-PD) model for quantitative exposure-response evaluation. To overcome the lack of exposure information, the K-PD model was implemented by the virtual pharmacokinetic compartment link to the pharmacodynamic compartment of Simeoni's tumor growth inhibition model and evaluated using curve-fitting for the tumor growth-time profile after intratumoral injection of BoNT/A1. The final K-PD model was adequately explained for a pattern of tumor growth depending on represented exposure parameters and simulation studies were conducted to determine the optimal dose under various scenarios considering dose strength and frequency. The optimal dose range and regimen of ≥13.8 units kg^-1 once a week or once every 3 days was predicted using the final model in B16-F10 syngeneic model and it was demonstrated with an extra in-vivo experiment. In conclusion, the K-PD model of BoNT/A1 was well developed to optimize the dosing regimen for evaluation of anticancer effect and this approach could be expandable to figure out quantitative interpretation of BoNT/A1's efficacy in various xenograft and/or syngeneic models.

Establishment and application of population pharmacokinetics model of vancomycin in infants with meningitis

BACKGROUND

To establish a population pharmacokinetics (PPK) model of vancomycin (VCM) for dose individualization in Chinese infants with meningitis.

METHODS

We collected the data of 82 children with meningitis in hospital from July 2014 to June 2016. The initial vancomycin dosage regimen for children was 10 or 15 mg/kg for q12 h, q8 h or q6 h. Serum concentrations were determined by Viva-E Analyzer before and after the fifth administration. The PPK model was developed by nonlinear mixed-effect model software, assessed by the bootstrap method and then tested in 20 infant patients.

RESULTS

The VCM clearance (CL) was increased by body weight (WT) and decreased by blood urea nitrogen (BUN). Pharmacokinetic parameters of VCM were not influenced by co-administered drugs. The trough concentrations of VCM were accurately predicted by the PPK model, with the prediction errors less than 32%.

CONCLUSION

A new individual strategy for VCM regimens was proposed and validated by the PPK model.

A longitudinal model for the Mayo Clinical Score and its sub-components in patients with ulcerative colitis

Clinical trials in patients with ulcerative colitis (UC) face the challenge of high and variable placebo response rates. The Mayo Clinical Score (MCS) is used widely as the primary endpoint in clinical trials to describe the clinical status of patients with UC. The MCS is comprised of four subscores, each scored 0, 1, 2 and 3: rectal bleeding (RB), stool frequency (SF), physician's global assessment (PGA), and endoscopy (ENDO) subscore. Excluding the PGA subscore gives the modified MCS. Quantitative insight on the placebo response, and its impact on the components of the MCS over time, can better inform clinical trial design and interpretation. Longitudinal modeling of the MCS, and the modified MCS, can be challenging due to complex clinical trial design, population heterogeneity, and limited assessments for the ENDO subscore. The current study pooled patient-level placebo/standard of care (SoC) arm data from five clinical trials in the TransCelerate database to develop a longitudinal placebo response model that describes the MCS over time in patients with UC. MCS subscores were modeled using proportional odds models, and the removal of patients from the placebo/SoC arm, or "dropout", was modeled using logistic regression models. The subscore and dropout models were linked to allow for the prediction of the MCS and the modified MCS. Stepwise covariate modeling identified prior exposure to TNF-α antagonists as a statistically significant predictor on the RB + SF subscore. Patients with prior exposure to TNF-α antagonists had higher post-baseline RB + SF subscores than naive patients.

Predictive Value of C-Reactive Protein and Albumin for Temporal Within-Individual Pharmacokinetic Variability of Voriconazole in Pediatric Hematopoietic Cell Transplant Patients

Voriconazole is a widely used antifungal agent in immunocompromised patients, but its utility is limited by its variable exposure and narrow therapeutic index. Population pharmacokinetic (PK) models have been used to characterize voriconazole PK and derive individualized dosing regimens. However, determinants of temporal within-patient variability of voriconazole PK were not well-established. We aimed to characterize temporal variability of voriconazole PK within individuals and identify predictive clinical factors. This study was conducted as a part of a single-institution, phase I study of intravenous voriconazole in children undergoing HCT (NCT02227797). We analyzed voriconazole PK study data collected at week 1 and again at week 2 after the start of voriconazole therapy in 59 pediatric HCT patients (age <21 years). Population PK analysis using nonlinear mixed effect modeling was performed to analyze temporal within-individual variability of voriconazole PK by incorporating a between-occasion variability term in the model. A two-compartment linear elimination model incorporating body weight and CYP2C19 phenotype described the data. Ratio of individual voriconazole clearance between weeks 1 to 2 ranged from 0.11 to 3.3 (-9.1 to +3.3-fold change). Incorporation of covariate effects by serum C-reactive protein (CRP) and albumin levels decreased between-occasion variability of clearance (coefficient of variation: from 59.5% to 41.2%) and improved the model fit (p<0.05). As significant covariates on voriconazole PK, CRP and albumin concentrations may potentially serve as useful biomarkers as part of therapeutic drug monitoring. This article is protected by copyright. All rights reserved.

Evaluation of sex differences in the pharmacokinetics of oral sumatriptan in healthy Korean subjects using population pharmacokinetic modeling

Sumatriptan was introduced in 1983, as the first of the triptans, selective 5‐hydroxytryptamine (5‐HT_1B/1D) receptor agonists, to treat moderate to severe migraine. Migraine predominates in females. Although there have been reports of sex differences in migraine‐associated features and pharmacokinetics (PKs) of some triptans, sex differences in the PKs of oral sumatriptan have never been evaluated in Korean. We conducted this study of oral sumatriptan to assess the sex differences in Korean population. Thirty‐eight healthy Korean subjects who participated in two separate clinical studies receiving a single oral dose of 50 mg sumatriptan with the same protocols were included in this analysis. A total of 532 sumatriptan concentration observations were used for a population PK modeling. Validation of final population PK model of sumatriptan was performed using bootstrap and visual predictive check. The PK profile of oral sumatriptan was adequately described by a one‐compartmental model with combined transit compartment model and a first‐order absorption. The covariate analysis showed that the clearance of oral sumatriptan was significantly higher in males than in females (male: 444 L/h, female: 281 L/h). Our results showed that there were sex differences in the clearance of oral sumatriptan. These results encourage further studies to establish the sumatriptan pharmacokinetic–pharmacodynamic model considering sex‐related PK differences, which may help to determine optimal dosing regimens for effective treatment of migraine in males and females. Clinical trial registration: CRIS Registration No. KCT0001784.