skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on March 1, 2026

Title: Intranasal diamorphine population pharmacokinetics modeling and simulation in pediatric breakthrough pain
Abstract Intranasal diamorphine (IND), approved for managing breakthrough pain in the UK, has been identified as an acceptable alternative offering effective, expedient, and less traumatic analgesia for children. However, the current dose regimen in pediatric populations relies on clinical expertise while the pharmacokinetics properties are poorly understood. This study aimed to develop diamorphine population pharmacokinetics (pop‐PK) models and simulate the IND dosing in virtual pediatric subjects. An integrated four‐compartment pop‐PK model with first‐order absorption and elimination provided an appropriate fit and characterized publicly available 385 concentration measurements of diamorphine, 6‐monoacetylmorphine, and morphine collected from adults. Body weight allometry and renal function maturation (age) were incorporated into the final model, serving as two covariates. The estimated IND relative bioavailability was around 52% compared with intramuscularly injected diamorphine. Using this final model, the morphine plasma concentrations, as the active metabolite for pain relief, were simulated in virtual subjects. The utility of model extrapolation was supported by external verification with acceptable average fold errors of 1.06 ± 0.30 and 0.83 ± 0.07 for morphine maximum concentration and exposures. Meanwhile, the simulated morphine concentration–time profiles could recover the PK profiles observed in children after a single dose of IND. The model‐based dosing simulations were therefore assessed in four children age groups to match the therapeutic window of morphine concentrations in steady state (10–20 μg/L). Our study demonstrates that the dose regimen of 0.3 mg/kg loading dose plus 0.1 mg/kg hourly maintenance dose is generally appropriate for multiple pediatric populations with breakthrough pain, in the view of PK.  more » « less
Award ID(s):
1955260
PAR ID:
10627177
Author(s) / Creator(s):
; ; ; ; ; ;
Publisher / Repository:
ASCPT
Date Published:
Journal Name:
CPT: Pharmacometrics & Systems Pharmacology
Volume:
14
Issue:
3
ISSN:
2163-8306
Page Range / eLocation ID:
435 to 447
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Malaria is a severe parasite infectious disease with high fatality. As one of the approved treatments of this disease, hydroxychloroquine (HCQ) lacks clinical administration guidelines for patients with special health conditions and co-morbidities. This may result in improper dosing for different populations and lead them to suffer from severe side effects. One of the most important toxicities of HCQ overdose is cardiotoxicity. In this study, we built and validated a physiologically based pharmacokinetic modeling (PBPK) model for HCQ. With the full-PBPK model, we predicted the pharmacokinetic (PK) profile for malaria patients without other co-morbidities under the HCQ dosing regimen suggested by Food and Drug Administration (FDA) guidance. The PK profiles for different special populations were also predicted and compared to the normal population. Moreover, we proposed a series of adjusted dosing regimens for different populations with special health conditions and predicted the concentration-time (C-T) curve of the drug plasma concentration in these populations which include the pregnant population, elderly population, RA patients, and renal impairment populations. The recommended special population-dependent dosage regimens can maintain the similar drug levels observed in the virtual healthy population under the original dosing regimen provided by FDA. Last, we developed mathematic formulas for predicting dosage based on a patient’s body measurements and two indexes of renal function (glomerular filtration rate and serum creatine level) for the pediatric and morbidly obese populations. Those formulas can facilitate personalized treatment of this disease. We hope to provide some advice to clinical practice when taking HCQ as a treatment for malaria patients with special health conditions or co-morbidities so that they will not suffer from severe side effects due to higher drug plasma concentration, especially cardiotoxicity. 
    more » « less
  2. Unique characteristics of the naked mole-rat (NMR) have made it increasingly popular as a laboratory animal model. These rodents are used to study many fields of research including longevity and aging, cancer, circadian rhythm, pain, and metabolism. Currently, the analgesic dosing regimens used in the NMR mirror those used in other rodent species. However, there is no pharmacokinetic (PK) data supporting the use of injectable analgesics in the NMR. Therefore, we conducted two independent PK studies to evaluate two commonly used analgesics in the NMR; meloxicam (2 mg/kg SC) and buprenorphine (0.1 mg/kg SC). In each study, blood was collected at 8 time points after subcutaneous injection of meloxicam or buprenorphine (0 (pre-dose), 0.25, 0.5, 1, 2, 4, 8, and 24 hrs). Three NMRs were used per time point for a total of 24 animals per PK study. Plasma concentrations of meloxicam were highest between 0.5 hrs and 1 hr post-injection. Levels remained above the extrapolated dog and cat therapeutic threshold levels (390-911 ng/mL) for at least 24 hrs. Plasma concentrations of buprenorphine were highest between 0.25 and 0.5 hrs post-injection. Levels remained above the human therapeutic threshold (1 ng/mL) for up to 21 hrs. No skin reactions were seen in association with injection of either drug. In summary, this data supports dosing meloxicam (2 mg/kg SC) once every 24 hrs and buprenorphine (0.1 mg/kg SC) once every 8-12 hrs in the NMR. Further studies should be performed to evaluate the clinical efficacy of these drugs by correlating plasma concentrations with post-operative pain assessments. 
    more » « less
  3. Suckow, Mark (Ed.)
    Abstract Unique characteristics of the naked mole-rat (NMR) have made it increasingly popular as a laboratory animal model. These rodents are used to study many fields of research including longevity and aging, cancer, circadian rhythm, pain, and metabolism. Currently, the analgesic dosing regimens used in the NMR mirror those used in other rodent species. However, there is no pharmacokinetic (PK) data supporting the use of injectable analgesics in the NMR. Therefore, we conducted 2 independent PK studies to evaluate 2 commonly used analgesics in the NMR: meloxicam (2 mg/kg SC) and buprenorphine (0.1 mg/kg SC). In each study, blood was collected at 8 time points after subcutaneous injection of meloxicam or buprenorphine (0 [predose], 0.25, 0.5, 1, 2, 4, 8, and 24 h). Three NMRs were used per time point for a total of 24 animals per PK study. Plasma concentrations of meloxicam were highest between 0.5 and 1 h postinjection. Levels remained above the extrapolated dog and cat therapeutic threshold levels (390 to 911 ng/mL) for at least 24 h. Plasma concentrations of buprenorphine were highest between 0.25 and 0.5 h postinjection. Levels remained above the human therapeutic threshold (1 ng/mL) for up to 21 h. No skin reactions were seen in association with injection of either drug. In summary, these data support dosing meloxicam (2 mg/kg SC) once every 24 h and buprenorphine (0.1 mg/kg SC) once every 8 to 12 h in the NMR. Further studies should be performed to evaluate the clinical efficacy of these drugs by correlating plasma concentrations with postoperative pain assessments. 
    more » « less
  4. Population pharmacokinetic (PK) modeling has become a cornerstone of drug development and optimal patient dosing. This approach offers great benefits for datasets with sparse sampling, such as in pediatric patients, and can describe between-patient variability. While most current algorithms assume normal or log-normal distributions for PK parameters, we present a mathematically consistent nonparametric maximum likelihood (NPML) method for estimating multivariate mixing distributions without any assumption about the shape of the distribution. This approach can handle distributions with any shape for all PK parameters. It is shown in convexity theory that the NPML estimator is discrete, meaning that it has finite number of points with nonzero probability. In fact, there are at most N points where N is the number of observed subjects. The original infinite NPML problem then becomes the finite dimensional problem of finding the location and probability of the support points. In the simplest case, each point essentially represents the set of PK parameters for one patient. The probability of the points is found by a primal-dual interior-point method; the location of the support points is found by an adaptive grid method. Our method is able to handle high-dimensional and complex multivariate mixture models. An important application is discussed for the problem of population pharmacokinetics and a nontrivial example is treated. Our algorithm has been successfully applied in hundreds of published pharmacometric studies. In addition to population pharmacokinetics, this research also applies to empirical Bayes estimation and many other areas of applied mathematics. Thereby, this approach presents an important addition to the pharmacometric toolbox for drug development and optimal patient dosing. 
    more » « less
  5. Jenner, Adrianne (Ed.)
    With the recent approval by the FDA of the first disease-modifying drug for Alzheimer’s Disease (AD), personalized medicine will be increasingly important for appropriate management and counseling of patients with AD and those at risk. The growing availability of clinical biomarker data and data-driven computational modeling techniques provide an opportunity for new approaches to individualized AD therapeutic planning. In this paper, we develop a new mathematical model, based on AD cognitive, cerebrospinal fluid (CSF) and MRI biomarkers, to provide a personalized optimal treatment plan for individuals. This model is parameterized by biomarker data from the AD Neuroimaging Initiative (ADNI) cohort, a large multi-institutional database monitoring the natural history of subjects with AD and mild cognitive impairment (MCI). Optimal control theory is used to incorporate time-varying treatment controls and side-effects into the model, based on recent clinical trial data, to provide a personalized treatment regimen with anti-amyloid-beta therapy. In-silico treatment studies were conducted on the approved treatment, aducanumab, as well as on another promising anti-amyloid-beta therapy under evaluation, donanemab. Clinical trial simulations were conducted over both short-term (78 weeks) and long-term (10 years) periods with low-dose (6 mg/kg) and high-dose (10 mg/kg) regimens for aducanumab, and a single-dose regimen (1400 mg) for donanemab. Results confirm those of actual clinical trials showing a large and sustained effect of both aducanumab and donanemab on amyloid beta clearance. The effect on slowing cognitive decline was modest for both treatments, but greater for donanemab. This optimal treatment computational modeling framework can be applied to other single and combination treatments for both prediction and optimization, as well as incorporate new clinical trial data as it becomes available. 
    more » « less