Population-based analysis of methadone distribution and metabolism using an age-dependent physiologically based pharmacokinetic model

Mathematical Modeling of the Gastrointestinal System for Preliminary Drug Absorption Assessment

The objective of this study is to demonstrate the potential of a multicompartmental mathematical model to simulate the activity of the gastrointestinal system after the intake of drugs, with a limited number of parameters. The gastrointestinal system is divided into five compartments, modeled as both continuous systems with discrete events (stomach and duodenum) and systems with delay (jejunum, ileum, and colon). The dissolution of the drug tablet occurs in the stomach and is described through the Noyes-Whitney equation, with pH dependence expressed through the Henderson-Hasselbach relationship. The boluses resulting from duodenal activity enter the jejunum, ileum, and colon compartments, where drug absorption takes place as blood flows countercurrent. The model includes only three parameters with assigned physiological meanings. It was tested and validated using data from in vivo experiments. Specifically, the model was tested with the concentration profiles of nine different drugs and validated using data from two drugs with varying initial concentrations. Overall, the outputs of the model are in good agreement with experimental data, particularly with regard to the time of peak concentration. The primary sources of discrepancy were identified in the concentration decay. The model's main strength is its relatively low computational cost, making it a potentially excellent tool for in silico assessment and prediction of drug adsorption in the intestine.

Short-acting versus long-acting opioids for pediatric postoperative pain management

INTRODUCTION

Opioids are potent analgesics commonly used to manage children's moderate to severe perioperative pain in children. A wide range of short and long-acting opioids are used to treat surgical pain and will be reviewed in this article.

AREAS COVERED

Both short- and long-acting opioids contain unique therapeutic benefits and adverse effects; however, due to the side effect profile and safety concerns, lack of familiarity, and evidence with long-acting opioids to treat surgical pain, shorter-acting opioids have traditionally been used in children. Almost all opioids work by binding to the mu receptor. Methadone, a long-acting opioid, is an exception because it also has beneficial N-methyl-D-aspartate antagonist properties. Clinically methadone's properties could translate to improved analgesic outcomes, reduced risk of adverse events, less risk for acute hyperalgesia, tolerance and abuse potential, faster recovery, and reduced risk for chronic persistent surgical pain. This review article summarizes and compares the evidence of commonly used short and long-acting opioids for perioperative pain control in the pediatric population.

EXPERT OPINION

Individualized methadone therapy using pharmacogenomics has the potential to transform opioid use in pain management by improving patient safety and analgesic outcomes, thereby addressing the gaps in current standardized ERAS protocols.

An Overview of Physiologically-Based Pharmacokinetic Models for Forensic Science

A physiologically-based pharmacokinetic (PBPK) model represents the structural components of the body with physiologically relevant compartments connected via blood flow rates described by mathematical equations to determine drug disposition. PBPK models are used in the pharmaceutical sector for drug development, precision medicine, and the chemical industry to predict safe levels of exposure during the registration of chemical substances. However, one area of application where PBPK models have been scarcely used is forensic science. In this review, we give an overview of PBPK models successfully developed for several illicit drugs and environmental chemicals that could be applied for forensic interpretation, highlighting the gaps, uncertainties, and limitations.

A pediatric quantitative systems pharmacology model of neurofilament trafficking in spinal muscular atrophy treated with the antisense oligonucleotide nusinersen

Phosphorylated neurofilament heavy subunit (pNfH) has been recently identified as a promising biomarker of disease onset and treatment efficacy in spinal muscular atrophy (SMA). This study introduces a quantitative systems pharmacology model representing the SMA pediatric scenario in the age range of 0-20 years with and without treatment with the antisense oligonucleotide nusinersen. Physiological changes typical of the pediatric age and the contribution of SMA and its treatment to the peripheral pNfH levels were included in the model by extending the equations of a previously developed mathematical model describing the neurofilament trafficking in healthy adults. All model parameters were estimated by fitting data from clinical trials that enrolled SMA patients treated with nusinersen. The data from the control group of the study was employed to build an in silico population of untreated subjects, and the parameters related to the treatment were estimated by fitting individual pNfH time series of SMA patients followed during the treatment. The final model reproduces well the pNfH levels in the presence of SMA in both the treated and untreated conditions. The results were validated by comparing model predictions with the data obtained from an additional cohort of SMA patients. The reported good predictive model performance makes it a valuable tool for investigating pNfH as a biomarker of disease progression and treatment response in SMA and for the in silico evaluation of novel treatment protocols.

Rational use of antibiotics in children with diabetic ketoacidosis needs attention

Diabetic ketoacidosis (DKA) in children may lead to acute kidney injury (AKI). Among 45 children with DKA in our center, eight cases had AKI on admission, and in one child, his kidney function did not recover until 3 mo after discharge. This child was treated with antibiotics (cephalosporin), and we cannot rule out delayed AKI recovery due to the combined effects of the drug and the disease. Pediatricians should be concerned about the impact of nephrotoxic drug and disease interactions on children's kidney function, and need to follow up children with DKA and AKI to determine the development of AKI.

Pharmacokinetic modeling of R and S-Methadone and their metabolites to study the effects of various covariates in post-operative children

Methadone is a synthetic opioid used as an analgesic and for the treatment of opioid abuse disorder. The analgesic dose in the pediatric population is not well-defined. The pharmacokinetics (PKs) of methadone is highly variable due to the variability in alpha-1 acid glycoprotein (AAG) and genotypic differences in drug-metabolizing enzymes. Additionally, the R and S enantiomers of methadone have unique PK and pharmacodynamic properties. This study aims to describe the PKs of R and S methadone and its metabolite 2-ethylidene-1,5-dimethyl-3,3-diphenylpyrrolidine (EDDP) in pediatric surgical patients and to identify sources of inter- and intra-individual variability. Children aged 8-17.9 years undergoing orthopedic surgeries received intravenous methadone 0.1 mg/kg intra-operatively followed by oral methadone 0.1 mg/kg postoperatively every 12 h. Pharmacokinetics of R and S methadone and EDDP were determined using liquid chromatography tandem mass spectrometry assays and the data were modeled using nonlinear mixed-effects modeling in NONMEM. R and S methadone PKs were well-described by two-compartment disposition models with first-order absorption and elimination. EDDP metabolites were described by one compartment disposition models with first order elimination. Clearance of both R and S methadone were allometrically scaled by bodyweight. CYP2B6 phenotype was a determinant of the clearance of both the enantiomers in an additive gene model. The intronic CYP3A4 single-nucleotide polymorphism (SNP) rs2246709 was associated with decreased clearance of R and S methadone. Concentrations of AAG and the SNP of AAG rs17650 independently increased the volume of distribution of both the enantiomers. The knowledge of these important covariates will aid in the optimal dosing of methadone in children.

Opioid Treatment for Neonatal Opioid Withdrawal Syndrome: Current Challenges and Future Approaches

Chronic intrauterine exposure to psychoactive drugs often results in neonatal opioid withdrawal syndrome (NOWS). When nonpharmacologic measures are insufficient in controlling NOWS, morphine, methadone, and buprenorphine are first-line medications commonly used to treat infants with NOWS because of in utero exposure to opioids. Research suggests that buprenorphine may be the leading drug therapy used to treat NOWS when compared with morphine and methadone. Currently, there are no consensus or standardized treatment guidelines for medications prescribed for NOWS. Opioids used to treat NOWS exhibit large interpatient variability in pharmacokinetics (PK) and pharmacodynamic (PD) response in neonates. Organ systems undergo rapid maturation after birth that may alter drug disposition and exposure for any given dose during development. Data regarding the PK and PD of opioids in neonates are sparse. Pharmacometric methods such as physiologically based pharmacokinetic and population pharmacokinetic modeling can be used to explore factors predictive of some of the variability associated with the PK/PD of opioids in newborns. This review discusses the utility of pharmacometric techniques for enhancing precision dosing in infants requiring opioid treatment for NOWS. Applying these approaches may contribute to optimizing the outcome by reducing cumulative drug exposure, mitigating adverse drug effects, and reducing the burden of NOWS in neonates.

Integrating in vitro data and physiologically based kinetic modeling-facilitated reverse dosimetry to predict human cardiotoxicity of methadone

Development of novel testing strategies to detect adverse human health effects is of interest to replace in vivo-based drug and chemical safety testing. The aim of the present study was to investigate whether physiologically based kinetic (PBK) modeling-facilitated conversion of in vitro toxicity data is an adequate approach to predict in vivo cardiotoxicity in humans. To enable evaluation of predictions made, methadone was selected as the model compound, being a compound for which data on both kinetics and cardiotoxicity in humans are available. A PBK model for methadone in humans was developed and evaluated against available kinetic data presenting an adequate match. Use of the developed PBK model to convert concentration–response curves for the effect of methadone on human-induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CM) in the so-called multi electrode array (MEA) assay resulted in predictions for in vivo dose–response curves for methadone-induced cardiotoxicity that matched the available in vivo data. The results also revealed differences in protein plasma binding of methadone to be a potential factor underlying variation between individuals with respect to sensitivity towards the cardiotoxic effects of methadone. The present study provides a proof-of-principle of using PBK modeling-based reverse dosimetry of in vitro data for the prediction of cardiotoxicity in humans, providing a novel testing strategy in cardiac safety studies.

A Sensitive LC-MS/MS Assay for the Quantification of Methadone and its Metabolites in Dried Blood Spots: Comparison With Plasma

INTRODUCTION

Methadone, a synthetic narcotic, is widely used both in adults and children for pain control and as a replacement drug in opioid use disorder to prevent craving and withdrawal. To support clinical pharmacokinetic trials in neonates, infants, and children, the authors developed and validated a novel, automated, highly sensitive liquid chromatography-electrospray-tandem mass spectrometry ionization (LC-ESI-MS/MS) method for the quantification of methadone and its metabolites, 2-ethylidene-1,5-dimethyl-3,3-diphenylpyrrolidine (EDDP) and 2-ethyl-5-methyl-3,3-diphenylpyraline (EMDP), in samples collected as dried blood spots.

METHODS

Blood was spiked with different concentrations of methadone, EDDP, and EMDP, and blood drops were applied to filter paper cards. Punches of 6.4 mm were removed from the cards, and 600 µL of protein precipitation solution (methanol/0.2M ZnSO4, 7:3, vol/vol) containing the internal standards (methadone-d9 and EDDP-d5) at a concentration of 1 mcg/L was added. The extracts were analyzed using LC-ESI-MS/MS in combination with online extraction. The mass spectrometer was run in the positive multiple reaction monitoring mode, and the total run time was 3.2 minutes.

RESULTS

For the dried blood spots, the assay has a lower limit of quantification of 0.1 mcg/L for methadone, EDDP, and EMDP. The range of reliable response for methadone for the ion transition m/z = 310.2→265.1 was 0.1-100 mcg/L and for the ion transition m/z = 310.2→223.1 5-1000 mcg/L. For EDDP, on the range of reliable response for the ion transition, m/z = 278.2→234.3 was 0.1-100 mcg/L and for the ion transition m/z = 278.2→186.1 5-1000 mcg/L. The calibration range for EMDP was 0.1-100 mcg/L. Accuracy (85%-115%) and imprecision (<15%) met predefined acceptance criteria.

DISCUSSION

This assay allows for the measurement of small volume blood samples without the need for an intravenous blood draw, and thus, it is suitable for pharmacokinetics studies and therapeutic drug monitoring in pediatric patients.

Pediatric Age Groups and Approach to Studies

Pediatric clinical trials are often requested according to specific age ranges. In the past and still today, these ages may correspond to developmental stages, such as newborn, infancy, childhood, and adolescence. Selection of ages for pediatric participation in medication studies should correspond to ages of rapid changes in pharmacokinetics and pharmacodynamics. Age-related changes in several enzymes involved in drug metabolism and glomerular filtration are described as examples of optimal ages for study of specific drugs according to their pathways of disposition.

Methadone dosing strategies in preterm neonates can be simplified

AIMS

A dramatic increase in newborn infants with neonatal abstinence syndrome has been observed and these neonates are frequently treated with complex methadone dosing schemes to control their withdrawal symptoms. Despite its abundant use, hardly any data on the pharmacokinetics (PK) of methadone is available in preterm neonates. Therefore we investigated developmental PK of methadone and evaluated current dosing strategies and possible simplification in this vulnerable population.

METHODS

A single-centre open-label prospective study was performed to collect PK data after a single oral dose of methadone in preterm neonates. A population PK model was built to characterize developmental PK of (R)- and (S)-methadone. Model-based simulations were performed to identify a simplified dosing strategy to reach and maintain target methadone exposure.

RESULTS

A total of 121 methadone concentrations were collected from 31 preterm neonates. A one-compartment model with first order absorption and elimination kinetics best described PK data for (R)- and (S)-methadone. Clearance increases with advancing gestational age and differs between R- and S-enantiomer, being slightly higher for the former (0.244 vs 0.167 L/h). Preterm neonates reached target exposure after 48 hours with currently used dosing schedules. Output from simulations revealed that target exposures can be achieved with a simplified dosing strategy during the first 4 days of treatment.

CONCLUSION

Methadone clearance in preterm neonates increases with advancing gestational age and its disposition is influenced by its chirality. Simulations that account for developmental PK changes indicate a shorter methadone dosing strategy can maintain target exposure to control withdrawal symptoms.

The Impact of Scaling Factor Variability on Risk-Relevant Pharmacokinetic Outcomes in Children: A Case Study Using Bromodichloromethane (BDCM)

Biotransformation rates extrapolated from in vitro data are used increasingly in human physiologically based pharmacokinetic (PBPK) models. This practice requires use of scaling factors, including microsomal content (mg of microsomal protein/g liver, MPPGL), enzyme specific content, and liver mass as a fraction of body weight (FVL). Previous analyses indicated that scaling factor variability impacts pharmacokinetic (PK) outcomes used in adult population dose-response studies. This analysis was extended to pediatric populations because large inter-individual differences in enzyme ontogeny likely would further contribute to scaling factor variability. An adult bromodichloromethane (BDCM) model (Kenyon, E. M., Eklund, C., Leavens, T. L., and Pegram, R. A. (2016a). Development and application of a human PBPK model for bromodichloromethane (BDCM) to investigate impacts of multi-route exposure. J. Appl. Toxicol. 36, 1095-1111) was re-parameterized for neonates, infants, and toddlers. Monte Carlo analysis was used to assess the impact of pediatric scaling factor variation on model-derived PK outcomes compared with adult findings. BDCM dose metrics were estimated following a single 0.05-liter drink of water or a 20-min bath, under typical (5 µg/l) and plausible higher (20 µg/l) BDCM concentrations. MPPGL, CYP2E1, and FVL values reflected the distribution of reported pediatric population values. The impact of scaling factor variability on PK outcome variation was different for each exposure scenario, but similar for each BDCM water concentration. The higher CYP2E1 expression variability during early childhood was reflected in greater variability in predicted PK outcomes in younger age groups, particularly for the oral exposure route. Sensitivity analysis confirmed the most influential parameter for this variability was CYP2E1, particularly in neonates. These findings demonstrate the importance of age-dependent scaling factor variation used for in vitro to in vivo extrapolation of biotransformation rates.

A Spatio-Temporal Exposure-Hazard Model for Assessing Biological Risk and Impact

We developed a simulation model for quantifying the spatio-temporal distribution of contaminants (e.g., xenobiotics) and assessing the risk of exposed populations at the landscape level. The model is a spatio-temporal exposure-hazard model based on (i) tools of stochastic geometry (marked polygon and point processes) for structuring the landscape and describing the exposed individuals, (ii) a dispersal kernel describing the dissemination of contaminants from polygon sources, and (iii) an (eco)toxicological equation describing the toxicokinetics and dynamics of contaminants in affected individuals. The model was implemented in the briskaR package (biological risk assessment with R) of the R software. This article presents the model background, the use of the package in an illustrative example, namely, the effect of genetically modified maize pollen on nontarget Lepidoptera, and typical comparisons of landscape configurations that can be carried out with our model (different configurations lead to different mortality rates in the treated example). In real case studies, parameters and parametric functions encountered in the model will have to be precisely specified to obtain realistic measures of risk and impact and accurate comparisons of landscape configurations. Our modeling framework could be applied to study other risks related to agriculture, for instance, pathogen spread in crops or livestock, and could be adapted to cope with other hazards such as toxic emissions from industrial areas having health effects on surrounding populations. Moreover, the R package has the potential to help risk managers in running quantitative risk assessments and testing management strategies.

Physiologically based pharmacokinetic (PBPK) modeling and simulation in neonatal drug development: how clinicians can contribute

Introduction: Legal initiatives to stimulate neonatal drug development should be accompanied by development of valid research tools. Physiologically based (PB)-pharmacokinetic (PK) modeling and simulation are established tools, accepted by regulatory authorities. Consequently, PBPK holds promise to be a strong research tool to support neonatal drug development. Area covered: The currently available PBPK models still have poor predictive performance in neonates. Using an illustrative approach on distinct PK processes of absorption, distribution, metabolism, excretion, and real-world data in neonates, we provide evidence on the need to further refine available PBPK system parameters through generation and integration of new knowledge. This necessitates cross talk between clinicians and modelers to integrate knowledge (PK datasets, system knowledge, maturational physiology) or test and refine PBPK models. Expert opinion: Besides refining these models for 'small molecules', PBPK model development should also be more widely applied for therapeutic proteins and to determine exposure through breastfeeding. Researchers should also be aware that PBPK modeling in combination with clinical observations can also be used to elucidate age-related changes that are almost impossible to study based on in vivo or in vitro data. This approach has been explored for hepatic biliary excretion, renal tubular activity, and central nervous system exposure.

Hybrid Statistical and Mechanistic Mathematical Model Guides Mobile Health Intervention for Chronic Pain

Nearly a quarter of visits to the emergency department are for conditions that could have been managed via outpatient treatment; improvements that allow patients to quickly recognize and receive appropriate treatment are crucial. The growing popularity of mobile technology creates new opportunities for real-time adaptive medical intervention, and the simultaneous growth of "big data" sources allows for preparation of personalized recommendations. Here we focus on the reduction of chronic suffering in the sickle cell disease (SCD) community. SCD is a chronic blood disorder in which pain is the most frequent complication. There currently is no standard algorithm or analytical method for real-time adaptive treatment recommendations for pain. Furthermore, current state-of-the-art methods have difficulty in handling continuous-time decision optimization using big data. Facing these challenges, in this study, we aim to develop new mathematical tools for incorporating mobile technology into personalized treatment plans for pain. We present a new hybrid model for the dynamics of subjective pain that consists of a dynamical systems approach using differential equations to predict future pain levels, as well as a statistical approach tying system parameters to patient data (both personal characteristics and medication response history). Pilot testing of our approach suggests that it has significant potential to well predict pain dynamics, given patients reported pain levels and medication usages. With more abundant data, our hybrid approach should allow physicians to make personalized, data-driven recommendations for treating chronic pain.

Conversion from prolonged intravenous fentanyl infusion to enteral methadone in critically ill children

AIM

To describe our institutional experience with conversion from intravenous (IV) fentanyl infusion directly to enteral methadone and occurrence of withdrawal in critically ill mechanically ventilated children exposed to prolonged sedation and analgesia.

METHODS

With Institutional Review Board approval, we retrospectively studied consecutively admitted invasively mechanically ventilated children (0-18 years) sedated with IV fentanyl infusion > 5 d and subsequently converted directly to enteral methadone. Data were obtained on subject demographics, illness severity, daily IV fentanyl and enteral methadone dosing, time to complete conversion, withdrawal scores (WAT-1), pain scores, and need for rescue opioids. Patients were classified as rapid conversion group (RCG) if completely converted ≤ 48 h and slow conversion group (SCG) if completely converted in > 48 h. Primary outcome was difference in WAT-1 scores at 7 d. Secondary outcomes included differences in overall pain scores, and differences in daily rescue opioids.

RESULTS

Compared to SCG (n = 21), RCG (n = 21) had lower median WAT-1 scores at 7 d (2.5 vs 5, P = 0.027). Additionally, RCG had lower overall median pain scores (3 vs 6, P = 0.007), and required less median daily rescue opioids (3 vs 12, P = 0.003) than SCG. The starting daily median methadone dose was 2.3 times the daily median fentanyl dose in the RCG, compared to 1.1 times in the SCG (P = 0.049).

CONCLUSION

We observed wide variation in conversion from IV fentanyl infusion directly to enteral methadone and variability in withdrawal in critically ill mechanically ventilated children exposed to prolonged sedation. In those children who converted successfully from IV fentanyl infusion to enteral methadone within a period of 48 h, a methadone:fentanyl dose conversion ratio of approximately 2.5:1 was associated with less withdrawal and reduced need for rescue opioids.

Physiologically-based pharmacokinetic models: approaches for enabling personalized medicine

Personalized medicine strives to deliver the 'right drug at the right dose' by considering inter-person variability, one of the causes for therapeutic failure in specialized populations of patients. Physiologically-based pharmacokinetic (PBPK) modeling is a key tool in the advancement of personalized medicine to evaluate complex clinical scenarios, making use of physiological information as well as physicochemical data to simulate various physiological states to predict the distribution of pharmacokinetic responses. The increased dependency on PBPK models to address regulatory questions is aligned with the ability of PBPK models to minimize ethical and technical difficulties associated with pharmacokinetic and toxicology experiments for special patient populations. Subpopulation modeling can be achieved through an iterative and integrative approach using an adopt, adapt, develop, assess, amend, and deliver methodology. PBPK modeling has two valuable applications in personalized medicine: (1) determining the importance of certain subpopulations within a distribution of pharmacokinetic responses for a given drug formulation and (2) establishing the formulation design space needed to attain a targeted drug plasma concentration profile. This review article focuses on model development for physiological differences associated with sex (male vs. female), age (pediatric vs. young adults vs. elderly), disease state (healthy vs. unhealthy), and temporal variation (influence of biological rhythms), connecting them to drug product formulation development within the quality by design framework. Although PBPK modeling has come a long way, there is still a lengthy road before it can be fully accepted by pharmacologists, clinicians, and the broader industry.

Neonatal abstinence syndrome: Pharmacologic strategies for the mother and infant

Opioid use in pregnancy has increased dramatically over the past decade. Since prenatal opioid use is associated with numerous obstetrical and neonatal complications, this now has become a major public health problem. In particular, in utero opioid exposure can result in neonatal abstinence syndrome (NAS) which is a serious condition characterized by central nervous system hyperirritability and autonomic nervous system dysfunction. The present review seeks to define current practices regarding the approach to the pregnant mother and neonate with prenatal opiate exposure. Although the cornerstone of prenatal management of opioid dependence is opioid maintenance therapy, the ideal agent has yet to be definitively established. Pharmacologic management of NAS is also highly variable and may include an opioid, barbiturate, and/or α-agonist. Genetic factors appear to be associated with the incidence and severity of NAS. Establishing pharmacogenetic risk factors for the development of NAS has the potential for creating opportunities for "personalized genomic medicine" and novel, individualized therapeutic interventions.

Pharmacokinetics of Oral Methadone in the Treatment of Neonatal Abstinence Syndrome: A Pilot Study

OBJECTIVE

To characterize the population pharmacokinetics of oral methadone in neonates requiring pharmacologic treatment of neonatal abstinence syndrome and to develop a pharmacokinetic (PK) model toward an evidence-based treatment protocol.

STUDY DESIGN

Based on a methadone dosing protocol, serum concentrations of methadone and its metabolites were assessed by high performance liquid chromatography-tandem mass spectrometry from dried blood spots. Population PK analysis was performed to determine the volume of distribution and clearance of oral methadone. Methadone plasma concentration-time profiles were simulated from the deduced PK model to optimize the dosing regimen.

RESULTS

There was substantial interindividual variability in methadone concentrations. Blood concentrations of methadone were best described by a 1-compartment model with first-order absorption. The population mean estimates (coefficient of variation percentage) for oral clearance and volume of distribution were 8.94 (103%) L/h/70 kg and 177 (133%) L/70 kg, respectively. Optimized dosing strategies were developed based on the simulated PK profiles. We suggest a starting dose of 0.1 mg/kg per dose every 6 hours for most patients requiring pharmacologic treatment of neonatal abstinence syndrome followed by an expedited weaning phase.

CONCLUSIONS

The proposed dosing regimen may reduce the cumulative dose of opioid and shorten the length of hospitalization. Future studies should aim to validate the simulated dosing schemes with clinical data and expand our understanding of the between-patient PK variability.

TRIAL REGISTRATION

ClinicalTrials.gov: NCT01754324.

Effect of Kidney Function on Drug Kinetics and Dosing in Neonates, Infants, and Children

Neonates, infants, and children differ from adults in many aspects, not just in age, weight, and body composition. Growth, maturation and environmental factors affect drug kinetics, response and dosing in pediatric patients. Almost 80% of drugs have not been studied in children, and dosing of these drugs is derived from adult doses by adjusting for body weight/size. As developmental and maturational changes are complex processes, such simplified methods may result in subtherapeutic effects or adverse events. Kidney function is impaired during the first 2 years of life as a result of normal growth and development. Reduced kidney function during childhood has an impact not only on renal clearance but also on absorption, distribution, metabolism and nonrenal clearance of drugs. 'Omics'-based technologies, such as proteomics and metabolomics, can be leveraged to uncover novel markers for kidney function during normal development, acute kidney injury, and chronic diseases. Pharmacometric modeling and simulation can be applied to simplify the design of pediatric investigations, characterize the effects of kidney function on drug exposure and response, and fine-tune dosing in pediatric patients, especially in those with impaired kidney function. One case study of amikacin dosing in neonates with reduced kidney function is presented. Collaborative efforts between clinicians and scientists in academia, industry, and regulatory agencies are required to evaluate new renal biomarkers, collect and share prospective pharmacokinetic, genetic and clinical data, build integrated pharmacometric models for key drugs, optimize and standardize dosing strategies, develop bedside decision tools, and enhance labels of drugs utilized in neonates, infants, and children.

Bayesian Population Physiologically-Based Pharmacokinetic (PBPK) Approach for a Physiologically Realistic Characterization of Interindividual Variability in Clinically Relevant Populations

Interindividual variability in anatomical and physiological properties results in significant differences in drug pharmacokinetics. The consideration of such pharmacokinetic variability supports optimal drug efficacy and safety for each single individual, e.g. by identification of individual-specific dosings. One clear objective in clinical drug development is therefore a thorough characterization of the physiological sources of interindividual variability. In this work, we present a Bayesian population physiologically-based pharmacokinetic (PBPK) approach for the mechanistically and physiologically realistic identification of interindividual variability. The consideration of a generic and highly detailed mechanistic PBPK model structure enables the integration of large amounts of prior physiological knowledge, which is then updated with new experimental data in a Bayesian framework. A covariate model integrates known relationships of physiological parameters to age, gender and body height. We further provide a framework for estimation of the a posteriori parameter dependency structure at the population level. The approach is demonstrated considering a cohort of healthy individuals and theophylline as an application example. The variability and co-variability of physiological parameters are specified within the population; respectively. Significant correlations are identified between population parameters and are applied for individual- and population-specific visual predictive checks of the pharmacokinetic behavior, which leads to improved results compared to present population approaches. In the future, the integration of a generic PBPK model into an hierarchical approach allows for extrapolations to other populations or drugs, while the Bayesian paradigm allows for an iterative application of the approach and thereby a continuous updating of physiological knowledge with new data. This will facilitate decision making e.g. from preclinical to clinical development or extrapolation of PK behavior from healthy to clinically significant populations.

Modeling and simulation in pediatric drug therapy: Application of pharmacometrics to define the right dose for children

During the past decades significant progress has been made in our understanding of the importance of age-appropriate development of new drug therapies in children. Importantly, several regulatory initiatives in Europe and the US have provided a framework for a rationale. In the US, most notably the enactment of the Best Pharmaceuticals for Children Act (BPCA) and Product Research and Equity Act (PREA) has facilitated the studying of on-patent and off-patent drugs in children. The biggest challenge in pediatric studies is defining a safe and effective dose or dose range in a patient population that can span from premature neonates to adolescents. From a mechanism-based perspective, advances in the science of quantitative pharmacology and pharmacometrics have resulted in the development of model-based approaches to better describe and understand important age-related factors influencing drug disposition and response in pediatric patients. The application of modeling and simulation has been shown to result in better estimates of pediatric doses as evidenced by several studies, although the optimal approach is still being debated. The extrapolation of efficacy findings from adults to the pediatric population has streamlined the development process especially for studies in older children. However, a focus on developmental changes in neonates and infants as well as further developing a paradigm for conducting pharmacodynamic studies in neonates, infants, and children remain important unmet needs. In this overview we will review current approaches for age-appropriate dose selection and highlight ongoing efforts to define exposure-response and clinical outcome relationships across the pediatric age spectrum.

Variation in treatment of neonatal abstinence syndrome in US children's hospitals, 2004-2011

OBJECTIVE

Neonatal abstinence syndrome (NAS) is a drug withdrawal syndrome experienced by opioid-exposed infants. There is no standard treatment for NAS and surveys suggest wide variation in pharmacotherapy for NAS. Our objective was to determine whether different pharmacotherapies for NAS are associated with differences in outcomes and to determine whether pharmacotherapy and outcome vary by hospital.

STUDY DESIGN

We used the Pediatric Health Information System Database from 2004 to 2011 to identify a cohort of infants with NAS requiring pharmacotherapy. Mixed effects hierarchical negative binomial models evaluated the association between pharmacotherapy and hospital with length of stay (LOS), length of treatment (LOT) and hospital charges, after adjusting for socioeconomic variables and comorbid clinical conditions.

RESULT

Our cohort included 1424 infants with NAS from 14 children's hospitals. Among hospitals in our sample, six used morphine, six used methadone and two used phenobarbital as primary initial treatment for NAS. In multivariate analysis, when compared with NAS patients initially treated with morphine, infants treated with methadone had shorter LOT (incidence rate ratio (IRR) = 0.55; P < 0.0001) and LOS (IRR = 0.60; P < 0.0001). Phenobarbital as a second-line agent was associated with increased LOT (IRR = 2.09; P<0.0001), LOS (IRR = 1.78; P < 0.0001) and higher hospital charges (IRR = 1.84; P < 0.0001). After controlling for case-mix, hospitals varied in LOT, LOS and hospital charges.

CONCLUSION

We found variation in hospital in treatment for NAS among major US children's hospitals. In analyses controlling for possible confounders, methadone as initial treatment was associated with reduced LOT and hospital stay.

Developmental pharmacokinetics in pediatric populations

Information on drug absorption and disposition in infants and children has increased considerably over the past 2 decades. However, the impact of specific age-related effects on pharmacokinetics, pharmacodynamics, and dose requirements remains poorly understood. Absorption can be affected by the differences in gastric pH and stomach emptying time that have been observed in the pediatric population. Low plasma protein concentrations and a higher body water composition can change drug distribution. Metabolic processes are often immature at birth, which can lead to a reduced clearance and a prolonged half-life for those drugs for which metabolism is a significant mechanism for elimination. Renal excretion is also reduced in neonates due to immature glomerular filtration, tubular secretion, and reabsorption. Limited data are available on the pharmacodynamic behavior of drugs in the pediatric population. Understanding these age effects provide a mechanistic way to identify initial doses for the pediatric population. The various factors that impact pharmacokinetics and pharmacodynamics mature towards adult values at different rates, thus requiring continual modification of drug dose regimens in neonates, infants, and children. In this paper, the age-related changes in drug absorption, distribution, metabolism, and elimination in infants and children are reviewed, and the age-related dosing regimens for this population are discussed.

The effectiveness of opioid substitution treatments for patients with opioid dependence: a systematic review and multiple treatment comparison protocol

BACKGROUND

Opioids are psychoactive analgesic drugs prescribed for pain relief and palliative care. Due to their addictive potential, effort and vigilance in controlling prescriptions is needed to avoid misuse and dependence. Despite the effort, the prevalence of opioid use disorder continues to rise. Opioid substitution therapies are commonly used to treat opioid dependence; however, there is minimal consensus as to which therapy is most effective. Available treatments include methadone, heroin, buprenorphine, as well as naltrexone. This systematic review aims to assess and compare the effect of all available opioid substitution therapies on the treatment of opioid dependence.

METHODS/DESIGN

The authors will search Medline, EMBASE, PubMed, PsycINFO, Web of Science, Cochrane Library, Cochrane Clinical Trials Registry, World Health Organization International Clinical Trials Registry Platform Search Portal, and the National Institutes for Health Clinical Trials Registry. The title, abstract, and full-text screening will be completed in duplicate. When appropriate, multiple treatment comparison Bayesian meta-analytic methods will be performed to deduce summary statistics estimating the effectiveness of all opioid substitution therapies in terms of retention and response to treatment (as measured through continued opioid abuse).

DISCUSSION

Using evidence gained from this systematic review, we anticipate disseminating an objective review of the current available literature on the effectiveness of all opioid substitution therapies for the treatment of opioid use disorder. The results of this systematic review are imperative to the further enhancement of clinical practice in addiction medicine.

SYSTEMATIC REVIEW REGISTRATION

PROSPERO CRD42013006507.

The pharmacokinetics of methadone and its metabolites in neonates, infants, and children

BACKGROUND

The lack of methadone pharmacokinetic data in children and neonates restrains dosing to achieve the target concentration in these populations. A minimum effective analgesic concentration of methadone in opioid naïve adults is 0.058 mg·l(-1) , while no withdrawal symptoms were observed in neonates suffering opioid withdrawal if plasma concentrations of methadone were above 0.06 mg·l(-1) . The racemate of methadone which is commonly used in pediatric and anesthetic care is metabolized to 2-ethylidine-1,5-dimethyl-3,3-diphenylpyrrolidine (EDDP) and 2-ethyl-5-methyl-3,3-diphenylpyrroline (EMDP).

METHODS

Data from four studies (age 33-week PMA-15 years) were pooled (n = 56) for compartment analysis using nonlinear mixed effects modeling. Parameter estimates were standardized to a 70-kg person using an allometric model approach. Investigation was made of the racemate and metabolite (EDDP and EMDP) dispositions. In addition, neonatal data (n = 7) allowed further study of R- and S-enantiomer pharmacokinetics.

RESULTS

A three-compartment linear disposition model best described the observed time-concentration profiles with additional compartments for metabolites. Population parameter estimates (between-subject variability) were central volume (V1) 21.5 (29%) l.70 kg(-1) , peripheral volumes of distribution V2 75.1 (23%) l.70 kg(-1) and V3 484 (8%) l.70 kg(-1) , clearance (CL) 9.45 (11%) l·h(-1) .70 kg(-1) , and intercompartment clearances Q2 325 (21%) l·h(-1) .70 kg(-1) and Q3 136 (14%) l·h(-1) .70 kg(-1) . EDDP formation clearance was 9.1 (11%) l·h(-1) .70 kg(-1) , formation clearance of EMDP from EDDP 7.4 (63%) l·h(-1) .70 kg(-1) , elimination clearance of EDDP was 40.9 (26%) l·h(-1) .70 kg(-1) and the rate constant for intermediate compartments 2.17 (43%) h(-1) .

CONCLUSIONS

Current pharmacokinetic parameter estimates in children and neonates are similar to those reported in adults. There was no clearance maturation with age. Neonatal enantiomer clearances were similar to those described in adults. A regimen of 0.2 mg·kg(-1) per 8 h in neonates achieves a target concentration of 0.06 mg·l(-1) within 36 h. Infusion, rather than intermittent dosing, should be considered if this target is to be achieved in older children after cardiac surgery.

Evaluation of a pharmacist-managed methadone taper*

OBJECTIVES

To evaluate the efficacy of a pharmacist-managed methadone taper as compared to previous prescribing practices.

DESIGN

Retrospective chart review with comparison to historical controls.

SETTING

Ninety-six-bed pediatric quaternary care facility with the majority of the patients in a 24-bed ICU.

PATIENTS

Thirty-two patients, 1 month to 16 years old, weaned off opioids using methadone prior to protocol initiation were compared with 20 patients, 1 month to 15 years old, weaned per the pharmacist-managed methadone taper protocol.

INTERVENTIONS

Implementation of a pharmacist-managed methadone tapering protocol.

MEASUREMENTS AND MAIN RESULTS

Patient age, gender, weight, and risk of withdrawal were similar between the groups (p = not significant). The average taper length before pharmacist intervention was 24.7 days; after implementation it was 15 days (p = 0.0026). There was no statistical difference in the number of additional doses of opioid required or withdrawal scores between the populations. Opioid infusions were stopped, on average, 1.54 days sooner in the intervention group (p = 0.0039). A decrease in hospital length of stay was also seen in the intervention group (p = 0.023).

CONCLUSIONS

A pharmacist-managed protocolized methadone taper facilitates discontinuing methadone sooner, discontinuing additional opioid infusions sooner, and may impact on the length of hospital stay. The protocol developed at the University of Minnesota Amplatz Children's Hospital has demonstrated that it is an effective method to wean pediatric patients from opioids.

Optimal sampling times for a drug and its metabolite using SIMCYP(®) simulations as prior information

BACKGROUND

Since 2007, it is mandatory for the pharmaceutical companies to submit a Paediatric Investigation Plan to the Paediatric Committee at the European Medicines Agency for any drug in development in adults, and it often leads to the need to conduct a pharmacokinetic study in children. Pharmacokinetic studies in children raise ethical and methodological issues. Because of limitation of sampling times, appropriate methods, such as the population approach, are necessary for analysis of the pharmacokinetic data. The choice of the pharmacokinetic sampling design has an important impact on the precision of population parameter estimates. Approaches for design evaluation and optimization based on the evaluation of the Fisher information matrix (M(F)) have been proposed and are now implemented in several software packages, such as PFIM in R.

OBJECTIVES

The objectives of this work were to (1) develop a joint population pharmacokinetic model to describe the pharmacokinetic characteristics of a drug S and its active metabolite in children after intravenous drug administration from simulated plasma concentration-time data produced using physiologically based pharmacokinetic (PBPK) predictions; (2) optimize the pharmacokinetic sampling times for an upcoming clinical study using a multi-response design approach, considering clinical constraints; and (3) evaluate the resulting design taking data below the lower limit of quantification (BLQ) into account.

METHODS

Plasma concentration-time profiles were simulated in children using a PBPK model previously developed with the software SIMCYP(®) for the parent drug and its active metabolite. Data were analysed using non-linear mixed-effect models with the software NONMEM(®), using a joint model for the parent drug and its metabolite. The population pharmacokinetic design, for the future study in 82 children from 2 to 18 years old, each receiving a single dose of the drug, was then optimized using PFIM, assuming identical times for parent and metabolite concentration measurements and considering clinical constraints. Design evaluation was based on the relative standard errors (RSEs) of the parameters of interest. In the final evaluation of the proposed design, an approach was used to assess the possible effect of BLQ concentrations on the design efficiency. This approach consists of rescaling the M(F), using, at each sampling time, the probability of observing a concentration BLQ computed from Monte-Carlo simulations.

RESULTS

A joint pharmacokinetic model with three compartments for the parent drug and one for its active metabolite, with random effects on four parameters, was used to fit the simulated PBPK concentration-time data. A combined error model best described the residual variability. Parameters and dose were expressed per kilogram of bodyweight. Reaching a compromise between PFIM results and clinical constraints, the optimal design was composed of four samples at 0.1, 1.8, 5 and 10 h after drug injection. This design predicted RSE lower than 30 % for the four parameters of interest. For this design, rescaling M(F) for BLQ data had very little influence on predicted RSE.

CONCLUSION

PFIM was a useful tool to find an optimal sampling design in children, considering clinical constraints. Even if it was not forecasted initially by the investigators, this approach showed that it was really necessary to include a late sampling time for all children. Moreover, we described an approach to evaluate designs assuming expected proportions of BLQ data are omitted.

Attitudes and practices of pediatric oncologists regarding methadone use in the treatment of cancer-related pain: results of a North American Survey

Methadone is effective in the treatment of cancer-related pain in adults. Pediatric oncologists may be reluctant to use methadone, given the paucity of existing research and a lack of familiarity with its use. This study's purpose was to assess pediatric oncologists' experience, comfort and practice of methadone prescription, and determine interest in and appropriate venues for education on methadone. A 22-item survey was sent by electronic mail to 1912 practicing pediatric oncologists. Six hundred thirty-one pediatric oncologists (33%) responded to the survey. Seventy-two percent of respondents reported they prescribe methadone to <10% of their patients receiving opioids. Physicians practicing ≥10 years (84% vs. 76%, P=0.01), at centers that see ≥100 new patients per year (86% vs. 76%, P=0.003), or who have received prior education on methadone (89% vs. 54%, P<0.001) were more likely to have prescribed methadone. The primary reasons respondents did not utilize methadone included a lack of knowledge of methadone's pharmacodynamics (39%), effectiveness (39%), and/or dosing equivalence (34%). Perceived competence with dose equivalence, belief that methadone is effective, and working in a division where >20 patients per year died were all independently associated with having prescribed methadone to >10% of patients on opioids. Eighty-five percent of respondents would like additional education on methadone. Many pediatric oncologists lack experience and education in the use of methadone. Formal education initiatives should be implemented to enhance pediatric oncologists' comfort and expertise in methadone use.

Strategic biomarkers for drug development in treating rare diseases and diseases in neonates and infants

There are similar challenges in developing a product designed to treat patients with a rare disease and drugs to treat critically ill neonates and infants. Part of the challenge in developing such products as well as identifying the optimal dosing regimen for the treatment of young children arises from the complex interrelationship between developmental changes and changes in biomarkers responsive to drug therapy. These difficulties are further compounded by our lack of understanding of the key physiological factors that cause the differences in clinical responses between adults and neonates and infants. Regulatory efforts have succeeded in overcoming these challenges in many areas of pediatric and orphan drug development. Strategic applications of biomarkers and surrogate endpoints for the development and approval of a product used to treat an orphan disease will be highlighted with examples of approved products. Continued efforts are still needed to fill in our knowledge gap and to strategically link biomarkers and surrogate endpoints to clinical responses for rare diseases and diseases affecting neonates and infants.

The pharmacokinetics of methadone in adolescents undergoing posterior spinal fusion

BACKGROUND

The optimal methadone dosing regimen for children undergoing spinal surgery is uncertain because of sparse pediatric pharmacokinetic data and a paucity of analgesic effect data. The minimum effective analgesic concentration of methadone in opioid naïve adults is 58 mcg · L(-1).

METHODS

Adolescents aged 12-19 years undergoing idiopathic scoliosis correction were administered 0.25 mg · kg(-1) racemic methadone IV prior to surgical incision. Arterial blood samples for methadone assay were obtained at 0 min, 5 min, 10 min, 15 min, 20 min, 40 min, 1 h, 2 h, 4 h, 5 h, 6 h, 8 h, 10 h, 12 h, 24 h, and 48 h. Compartment analysis was undertaken using nonlinear mixed effects models. Parameter estimates were standardized to a 70-kg person using allometric models.

RESULTS

A three-compartment linear disposition model best described observed time-concentration profiles. Population parameter estimates (between-subjects variability) were central volume (V1) 19.1 (126%) L 70 kg(-1), peripheral volumes of distribution V2 65.5 (60%) L 70 kg(-1), V3 485 (23%) L 70 kg(-1), clearance (CL) 9.3 (11%) L · h(-1) · 70 kg(-1), and inter-compartment clearances Q2 282 (95%) L · h(-1) 70 kg(-1), Q3 139 (42%) L · h(-1) 70 kg(-1). The terminal elimination half-life was 44.4 h. The mean observed methadone concentration was <58 mcg · L(-1) by the first hour after administration.

CONCLUSIONS

Current pharmacokinetic parameter estimates in adolescents are similar to those reported in adults. Methadone undergoes rapid redistribution after bolus administration. This may result in plasma concentrations that provide inadequate analgesia postoperatively. We would suggest following the bolus (0.25 mg.kg(-1)) with an infusion (0.1-0.15 mg · kg(-1) · h(-1) for 4 h) during spinal surgery to ensure adequate plasma concentrations for 24 h.

Pharmacologic management of the opioid neonatal abstinence syndrome

Opioid use in pregnant women has increased over the last decade. Following birth, infants with in utero exposure demonstrate signs and symptoms of withdrawal known as the neonatal abstinence syndrome (NAS). Infants express a spectrum of disease, with most requiring the administration of pharmacologic therapy to ensure proper growth and development. Treatment often involves prolonged hospitalization. There is a general lack of high-quality clinical trial data to guide optimal therapy, and significant heterogeneity in treatment approaches. Emerging trends in the treatment of infants with NAS include the use of sublingual buprenorphine, transition to outpatient therapy, and pharmacogenetic risk stratification.

Modeling interindividual variability in physiologically based pharmacokinetics and its link to mechanistic covariate modeling

Covariate modeling is a key step in the analysis of clinical data and is essential for establishing dosing recommendations for specific populations, e.g., in obese individuals and children. So far, no systematic approach exists to leverage the knowledge inherent in physiologically based pharmacokinetic (PBPK) models in this context. We introduce (i) a novel approach to model interindividual variability in PBPK models based on lean body weight (LBW); and (ii) a systematic approach to translate interindividual variability into the design of mechanistic covariate models. We derive a new covariate relation for the volume of distribution at steady state (V_ss) that seamlessly integrates body weight and LBW as covariates, with a weighting factor depending on the physicochemical properties of the drug. We further show that for children, PBPK-based extrapolation and allometric scaling result in very similar predictions for V_ss and blood clearance.

Perioperative pharmacokinetics of methadone in adolescents

BACKGROUND

Methadone is frequently administered to adults experiencing anesthesia and receiving pain treatment. Methadone pharmacokinetics in adults are well characterized, including the perioperative period. Methadone is also used in children. There is, however, no information on methadone pharmacokinetics in children of any age. The purpose of this investigation was to determine the pharmacokinetics of intravenous methadone in children undergoing surgery. Perioperative opioid-sparing effects were also assessed.

METHODS

Eligible subjects were children 5-18 yr undergoing general anesthesia and surgery, with an anticipated postoperative inpatient stay exceeding 3 days. Three groups of 10 to 11 patients each received intravenous methadone hydrochloride after anesthetic induction in ascending dose groups of 0.1, 0.2, and 0.3 mg/kg (up to 20 mg). Anesthetic care was not otherwise changed. Venous blood was obtained for 4 days, for stereoselective determination of methadone and metabolites. Pain assessments were made each morning. Daily and total opioid consumption was determined. Perioperative opioid consumption and pain was determined in a second cohort, which was matched to age, sex, race, ethnicity, surgical procedure, and length of stay, but not receiving methadone.

RESULTS

The final methadone study cohort was 31 adolescents (14 ± 2 yr, range 10-18) undergoing major spine surgery for a diagnosis of scoliosis. Methadone pharmacokinetics were linear over the dose range 0.1-0.3 mg/kg. Disposition was stereoselective. Methadone administration did not dose-dependently affect postoperative pain scores, and did not dose-dependently decrease daily or total postoperative opioid consumption in spinal fusion patients.

CONCLUSIONS

Methadone enantiomer disposition in adolescents undergoing surgery was similar to that in healthy adults.

Methadone use in children and young adults at a cancer center: a retrospective study

OBJECTIVE

To augment the literature on methadone applications in pediatric oncology, the authors reviewed the use of methadone at a pediatric cancer center over a 5-year period.

DESIGN AND SETTING

Forty-one patients received methadone for inpatient or outpatient pain management. The authors retrospectively reviewed their demographic characteristics, diagnoses, type of pain (nociceptive, neuropathic, or mixed) and causes of pain, and the indications, dose regimens, adverse effects, and outcomes of methadone treatment.

RESULTS

There were four types of clinical uses for methadone in 41 patients (10 patients had two): nociceptive pain unresponsive to other opioids (17 patients, 33.3 percent), neuropathic pain (20 patients, 39.2 percent), facilitation of weaning from opioids (11 patients, 21.6 percent), and end-of-life pain management (3 patients, 5.9 percent). The mean age of the 24 males (58.5 percent) and 17 females (41.5 percent) at the start of treatment was 15.7 years (range, 0.6-23 years). The most common diagnoses were leukemia (n = 10, 24.4 percent), osteosarcoma (n = 7, 17.0 percent), and rhabdomyosarcoma (n = 5, 12.2 percent). The causes of pain were bone marrow transplant (n = 13, 31.7 percent), amputation (n = 6, 14.6 percent), chemotherapy (n = 5, 12.2 percent), tumor (n = 5, 12.2 percent), limb-sparing surgery (n = 4, 9.8 percent), and other (n = 8, 19.5 percent). Efficacy was assessed at the end (or after 6 months) of methadone treatment. For many patients (43.1 percent), methadone showed efficacy in achieving the purpose for which it was prescribed, including reduction of nociceptive or neuropathic pain and prevention of opioid withdrawal. Sedation was the most common side effect (24.4 percent).

CONCLUSIONS

Methadone was effective for pediatric patients with neuropathic pain or nociceptive pain unresponsive to other opioids, and it effectively prevented opioid withdrawal.

Use of methadone for prevention of opioid withdrawal in critically ill children

BACKGROUND

Opioids are commonly administered to critically ill children for analgesia and sedation, but many patients experience opioid withdrawal upon discontinuation. The authors' institution developed a protocol for using methadone to prevent opioid withdrawal in children who have received morphine by continuous IV infusion for 5 days or longer in the pediatric intensive care unit (PICU).

OBJECTIVES

The primary objectives were to determine if opioids were tapered according to the protocol and to determine the conversion ratio for IV morphine to oral methadone that was used. Secondary objectives were to describe the methadone dosage used and the clinical outcomes, to evaluate adjustments to methadone dosing, and to report the incidence of adverse effects.

METHODS

A retrospective analysis of charts was conducted for pediatric patients who had received morphine by continuous IV infusion for 5 days or longer followed by methadone in the PICU between May 2008 and August 2009. Validated scoring systems (the Withdrawal Assessment Tool and the State Behavioral Scale) were used to assess symptoms of withdrawal and degree of sedation, respectively.

RESULTS

Forty-three patients were included in the study, with median age of 8 months (range 0.25-201 months). For 31 patients (72%), the protocol was not used, and there were no patients for whom the protocol was followed to completion. The median duration of weaning was 10 days (range 0-91 days). The conversion ratio for IV morphine to oral methadone was 1:0.78 for anticipated 5-day weaning and 1:0.98 for anticipated 10-day weaning. During the first 10 days of weaning, 18 patients (42%) experienced withdrawal symptoms. The methadone dose was increased for 11 (26%) of the 43 patients. Patients were sedated for a median of 1 day (range 0-9 days), were comfortable for a median of 6.5 days (range 1-64 days), and were agitated for a median of 2.5 days (range 0-23 days). Naloxone was required for 2 patients.

CONCLUSIONS

The institution's methadone protocol was not followed consistently during the study period, and practices for transitioning from morphine by continuous IV infusion to methadone with tapering were also inconsistent. Further studies are needed to determine the optimal conversion ratio for morphine to methadone and the optimal tapering regimen to minimize withdrawal symptoms and adverse events.

Physiologically based pharmacokinetic tissue compartment model selection in drug development and risk assessment

A well-stirred tank (WST) has been the predominant flow-limited tissue compartment model in physiologically based pharmacokinetic (PBPK) modeling. Recently, we developed a two-region asymptotically reduced (TAR) PBPK tissue compartment model through an asymptotic approximation to a two-region vascular-extravascular system to incorporate more biophysical detail than the WST model. To determine the relevance of a flow-limited TAR (F-TAR) approach, 75 structurally diverse drugs were evaluated herein using a priori predicted tissue:plasma partition coefficients along with hybrid and whole-body PBPK of eight rat tissues to determine the impact of model selection on simulation and optimization. Simulations showed that the F-TAR model significantly improved the ability to predict drug exposure, with hybrid and whole-body WST model error approaching 50% for tissues with larger vascular volumes. When optimization was used to fit F-TAR and WST models to pseudo data, WST-optimized drug partition coefficients more appropriately represented curve-fitting parameters rather than biophysically meaningful partition coefficients. Median F-TAR-optimized error ranged from -0.4% to +0.3%, whereas WST-optimized median error ranged from -22.2% to +1.8%. These studies demonstrated that the use of F-TAR represents a more accurate, biophysically realistic PBPK tissue model for predicting tissue exposure to drug and that it should be considered for use in drug development and regulatory review.

Physiologically based pharmacokinetic modeling: methodology, applications, and limitations with a focus on its role in pediatric drug development

The concept of physiologically based pharmacokinetic (PBPK) modeling was introduced years ago, but it has not been practiced significantly. However, interest in and implementation of this modeling technique have grown, as evidenced by the increased number of publications in this field. This paper demonstrates briefly the methodology, applications, and limitations of PBPK modeling with special attention given to discuss the use of PBPK models in pediatric drug development and some examples described in detail. Although PBPK models do have some limitations, the potential benefit from PBPK modeling technique is huge. PBPK models can be applied to investigate drug pharmacokinetics under different physiological and pathological conditions or in different age groups, to support decision-making during drug discovery, to provide, perhaps most important, data that can save time and resources, especially in early drug development phases and in pediatric clinical trials, and potentially to help clinical trials become more “confirmatory” rather than “exploratory”.

Development of appropriate equations for physiologically based pharmacokinetic modeling of permeability-limited and flow-limited transport

Although the implementation of a flow-limited, well-stirred tank (WST) single-compartment tissue model in pharmacokinetics and toxicokinetics is widespread, its use is not always justified biophysically or physiologically. The WST model introduces a loss of biophysical detail, specifically the vascular space, which is present in the standard permeability-limited two-subcompartment (PLT) tissue model. To address this loss of detail when evaluating the in vivo kinetics of drugs, toxins, nutrients, and endogenous metabolites, a novel set of physiologically based pharmacokinetic tissue compartment equations is developed through application of an asymptotic approximation to a two-region vascular-extravascular system to arrive at a permeability-limited two-region asymptotically reduced (P-TAR) model and a flow-limited (F-TAR) model. Development of the TAR modeling approach illustrates the importance of relative timescales in PBPK tissue compartment model selection and the conditions under which improved biophysical realism is advantageous. In the permeability-limited regime, the TAR model formulations enable drug or toxicant concentration to be modeled in the vascular and extravascular spaces equivalent to the PLT tissue model while invoking only one state variable to represent the vascular and extravascular spaces. In the flow-limited regime, the F-TAR model is more biophysically realistic than the WST model because it maintains the anatomical distinction between the vascular and extravascular spaces, and hence offers greater pharmacological and physiological insight than the WST model, without introducing additional computational complexity.