Morphine Glucuronidation and Elimination in Intensive Care Patients: A Comparison with Healthy Volunteers

The Drug Titration Paradox in the Presence of Intra-Individual Variation: Can we Estimate the True Concentration-Effect Relationship?

The drug titration paradox arises when higher drug concentrations are paradoxically associated with poorer efficacy outcomes, due to the titration of an individual's drug dose to achieve a desired effect. In cases with substantial intraindividual variability of the disease state, the drug titration paradox can also occur on the individual level (resulting in a higher dose when the individual has a worse disease state) and it has been suggested that it may not be possible to estimate the true exposure-response (ER) relationship in such situations. We simulated a titration study with strong intra-individual variability of disease state (causing the drug titration paradox at the individual level) and investigated the performance of four PKPD modelling methods in obtaining an unbiased estimate of the ER relationship. Strong bias in the estimated ER relationship was observed with two commonly used modelling methods: the model which only estimated inter-individual variability (IIV) and the model that included IIV and inter-occasion variability (IOV) on disease severity. In contrast, inclusion of stochastic differential equations (SDE) or accounting for the autocorrelation of the residual error between observations did yield successful estimation of the ER relationship without bias. The success of these methods can be understood from the principles of causal inference: confounding is avoided by controlling for the previous observations which drive the drug titration. Our results underline the importance of adequately characterizing intra-individual variability to avoid bias in PKPD modelling, especially for clinical areas where titration designs are common, such as analgesia.

Association of NR1I2 Polymorphism with Midazolam Clearance in Mechanically Ventilated ICU Patients: A Population Pharmacokinetic and Pharmacogenetic Study

Background

Significant variability in the metabolism of midazolam (MDZ) exists among mechanically ventilated (MV) patients in the intensive care unit (ICU) due to complex clinical conditions and genetic factors. The NR1I2 gene (PXR), which encodes a nuclear receptor that regulates drug-metabolizing enzymes like CYP3A4, plays a critical role in MDZ metabolism. Polymorphisms in NR1I2, along with variations in genes such as CYP3A4, CYP3A5, and ABCB1, may influence enzyme activity and MDZ pharmacokinetics (PK). Understanding these factors is essential for optimizing MDZ dosing in high-risk patient populations.

Methods

We studied 61 MV ICU patients receiving continuous MDZ infusion. A population pharmacokinetic (PopPK) model was used to assess MDZ PK, with genetic factors (NR1I2 rs2461817, CYP3A4, CYP3A5, ABCB1, and other PXR polymorphisms) and clinical covariates (body weight (BW), aspartate aminotransferase (AST) levels) evaluated for their impact on MDZ clearance (CL).

Results

The PK of MDZ and its metabolite, 1-hydroxymidazolam (1-OH-MDZ), were accurately described using a one-compartment model. The estimated population means for MDZ and 1-OH-MDZ CL were 22.6 L/h (inter-individual variability [IIV] 59.4%) and 67.1 L/h (IIV 57.7%), respectively. MDZ CL was significantly associated with the NR1I2 rs2461817 polymorphism and AST levels, accounting for 11.3% of the variability. MDZ CL decreased by 32.7% as AST increased from 22 IU/L to 60 IU/L, and by 40.7% in patients homozygous for the NR1I2 rs2461817 variant. BW also influenced the CL of 1-OH-MDZ, demonstrating a 34.2% increase as weight increased from 54 kg to 65 kg. Simulations confirmed the significant impact of NR1I2 rs2461817 on MDZ CL.

Conclusion

The PopPK model highlights the significant impact of NR1I2 rs2461817 polymorphism on MDZ CL in Chinese MV patients, emphasizing the need to consider genetic and clinical factors for optimizing MDZ dosing in ICU settings.

Pharmacokinetics of Opioid Infusions in the Adult Intensive Care Unit Setting-A Systematic Review

INTRODUCTION

Pharmacokinetics (PKs) of drugs are often altered in the intensive care unit (ICU). Opioids are often used in the ICU, particularly as continuous infusions, and their characteristics lead them to undergo PK alterations. We conducted a systematic review to assess the PK of opioid infusions in the ICU.

METHODS

Embase, MEDLINE, PubMed, CINAHL, and Evidence-Based Medicine Reviews (EBMR) were searched from inception to March 2024. Studies were included if they evaluated PKs of opioid infusions in adult patients in the ICU. Two reviewers independently selected and extracted data.

RESULTS

Out of the 1040 records screened, 17 studies were included. Five studies were conducted on fentanyl, seven on morphine, one on hydromorphone, two on remifentanil, two on alfentanil, and one on sufentanil. Most studies where observational studies or case series. The mean age was 56 years old. Duration of the infusion varied between 3 h and 20 days. PKs of fentanyl, sufentanil, and hydromorphone were significantly impaired, whereas the PKs of morphine, alfentanil, and remifentanil were impaired to a lesser degree. The PK parameter that was most affected by critical illness was the half-life (T½).

CONCLUSIONS

To counter these PK alterations, new therapeutic avenues must be further explored in the ICU to individualize opioid infusions.

Impact of CYP2D6*2, CYP2D6*35, rs5758550, and related haplotypes on risperidone clearance in vivo

PURPOSE

The CYP2D6 gene exhibits significant polymorphism, contributing to variability in responses to drugs metabolized by CYP2D6. While CYP2D6*2 and CYP2D6*35 are presently designated as alleles encoding normal metabolism, this classification is based on moderate level evidence. Additionally, the role of the formerly called "enhancer" single nucleotide polymorphism (SNP) rs5758550 is unclear. In this study, the impacts of CYP2D6*2, CYP2D6*35 and rs5758550 on CYP2D6 activity were investigated using risperidone clearance as CYP2D6 activity marker.

METHODS

A joint parent-metabolite population pharmacokinetic model was used to describe 1,565 serum concentration measurements of risperidone and 9-hydroxyrisperidone in 512 subjects. Risperidone population clearance was modeled as the sum of a CYP2D6-independent clearance term and the partial clearances contributed from each individually expressed CYP2D6 allele or haplotype. In addition to the well-characterized CYP2D6 alleles (*3-*6, *9, *10 and *41), *2, *35 and two haplotypes assigned as CYP2D6*2-rs5758550G and CYP2D6*2-rs5758550A were evaluated.

RESULTS

Each evaluated CYP2D6 allele was associated with significantly lower risperidone clearance than the reference normal function allele CYP2D6*1 (p < 0.001). Further, rs5758550 differentiated the effect of CYP2D6*2 (p = 0.005). The haplotype-specific clearances for CYP2D6*2-rs5758550A, CYP2D6*2-rs5758550G and CYP2D6*35 were estimated to 30%, 66% and 57%, respectively, relative to the clearance for CYP2D6*1. Notably, rs5758550 is in high linkage disequilibrium (R2 > 0.85) with at least 24 other SNPs and cannot be assigned as a functional SNP.

CONCLUSION

CYP2D6*2 and CYP2D6*35 encode reduced risperidone clearance, and the extent of reduction for CYP2D6*2 is differentiated by rs5758550. Genotyping of these haplotypes might improve the precision of genotype-guided prediction of CYP2D6-mediated clearance.

Pharmacokinetics and metabolism of lidocaine HCl 2% with epinephrine in horses following a palmar digital nerve block

BACKGROUND

Lidocaine is a local anesthetic that is sometimes administered in combination with epinephrine. The addition of epinephrine increases the time lidocaine remains at the site of administration, thus prolonging the duration of effect. Due to their potential to prevent the visual detection of lameness, the administration of local anesthetics is strictly regulated in performance and racehorses. Recent reports of positive regulatory findings for lidocaine in racehorses suggests a better understanding of the behavior of this drug is warranted. The objective of the current study was to describe serum and urine concentrations and the pharmacokinetics of lidocaine and its primary metabolites following administration in combination with epinephrine, as a palmar digital nerve block in horses. Twelve horses received a single administration of 1 mL of 2% lidocaine HCl (20 mg/horse) with epinephrine 1:100,000, over the palmar digital nerve. Blood samples were collected up to 30 h and urine samples up to 48 h post administration. Lidocaine and metabolite concentrations were determined by liquid chromatography- mass spectrometry and pharmacokinetic (non-compartmental and compartmental) analysis was performed.

RESULTS

Serum concentrations of lidocaine and 3-hydroxylidocaine were above the LOQ of the assay at 30 h post administration and monoethylglycinexylidide (MEGX) and glycinexylidide (GX) were below detectable levels by 24 and 48 h, respectively. In urine, lidocaine, MEGX and GX were all non-detectable by 48 h post administration while 3-hydroxylidocaine was above LOQ at 48 h post administration. The time of maximal concentration for lidocaine was 0.26 h (median) and the terminal half-life was 3.78 h (mean). The rate of absorption (Ka) was 1.92 1/h and the rate of elimination (Kel) was 2.21 1/h.

CONCLUSIONS

Compared to previous reports, the terminal half-life and subsequent detection time observed following administration of lidocaine in combination with epinephrine is prolonged. This is likely due to a decrease in systemic uptake of lidocaine because of epinephrine induced vasoconstriction. Results of the current study suggest it is prudent to use an extended withdrawal time when administering local anesthetics in combination with epinephrine to performance horses.

Population pharmacokinetics of buprenorphine and naloxone sublingual combination in Chinese healthy volunteers and patients with opioid use disorder: Model-based dose optimization

The sublingual combination of buprenorphine (BUP) and naloxone (NLX) is a new treatment option for opioid use disorder (OUD) and is effective in preventing drug abuse. This study aimed to explore rational dosing regimen for OUD patients in China via a model-based dose optimization approach. BUP, norbuprenorphine (norBUP), and NLX plasma concentrations of 34 healthy volunteers and 12 OUD subjects after single or repeated dosing were included. A parent-metabolite population pharmacokinetics (popPK) model with transit compartments for absorption was implemented to describe the pharmacokinetic profile of BUP-norBUP. In addition, NLX concentrations were well captured by a one-compartment popPK model. Covariate analysis showed that every additional swallow after the administration within the observed range (0-12) resulted in a 3.5% reduction in BUP bioavailability. This provides a possible reason for the less-than-dose proportionality of BUP. There were no differences in the pharmacokinetic characteristics between BUP or NLX in healthy volunteers and OUD subjects. Ethnic sensitivity analysis demonstrated that the dose-normalized peak concentration and area-under-the-curve of BUP in Chinese were about half of Puerto Ricans, which was consistent with a higher clearance observed in Chinese (166 L / h vs. 270 L / h ). Furthermore, Monte Carlo simulations showed that an 8 mg three-times daily dose was the optimized regimen for Chinese OUD subjects. This regimen ensured that opioid receptor occupancy remained at a maximum (70%) in more than 95% of subjects, at the same time, with NLX plasma concentrations below the withdrawal reaction threshold (4.6  n g / m L ).

Predicting Unacceptable Pain in Cardiac Surgery Patients Receiving Morphine Maintenance and Rescue Doses: A Model-Based Pharmacokinetic-Pharmacodynamic Analysis

BACKGROUND

Optimal analgesic treatment following cardiac surgery is crucial for both patient comfort and successful postoperative recovery. While knowledge of both the pharmacokinetics and pharmacodynamics of analgesics is required to predict optimal drug dosing, models quantifying the pharmacodynamics are scarce. Here, we quantify the pharmacodynamics of morphine by modeling the need for rescue morphine to treat unacceptable pain in 118 patients after cardiac surgery.

METHODS

The rescue morphine event data were analyzed with repeated time-to-event (RTTE) modeling using NONMEM. Postoperative pain titration protocol consisted of continuous morphine infusions (median duration 20.5 hours) with paracetamol 4 times daily and rescue morphine in case of unacceptable pain (numerical rating scale ≥4).

RESULTS

Patients had a median age of 73 years (interquartile range [IQR]: 63-77) and median bodyweight of 80 kg (IQR: 72-90 kg). Most patients (55%) required at least 1 rescue morphine dose. The hazard for rescue morphine following cardiac surgery was found to be significantly influenced by time after surgery, a day/night cycle with a peak at 23:00 (95% confidence interval [CI], 19:35-02:03) each day, and an effect of morphine concentration with 50% hazard reduction at 9.3 ng·mL-1 (95% CI, 6.7-16).

CONCLUSIONS

The pharmacodynamics of morphine after cardiac surgery was successfully quantified using RTTE modeling. Future studies can be used to expand the model to better predict morphine's pharmacodynamics on the individual level and to include the pharmacodynamics of other analgesics so that improved postoperative pain treatment protocols can be developed.

Enzyme Kinetics of Uridine Diphosphate Glucuronosyltransferases (UGTs)

Glucuronidation, catalyzed by uridine diphosphate glucuronosyltransferases (UGTs), is an important process for the metabolism and clearance of many lipophilic chemicals, including drugs, environmental chemicals, and endogenous compounds. Glucuronidation is a bisubstrate reaction that requires the aglycone and the cofactor, UDP-GlcUA. Accumulating evidence suggests that the bisubstrate reaction follows a compulsory-order ternary mechanism. To simplify the kinetic modeling of glucuronidation reactions in vitro, UDP-GlcUA is usually added to incubations in large excess. Many factors have been shown to influence UGT activity and kinetics in vitro, and these must be accounted for during experimental design and data interpretation. While the assessment of drug-drug interactions resulting from UGT inhibition has been challenging in the past, the increasing availability of UGT enzyme-selective substrate and inhibitor "probes" provides the prospect for more reliable reaction phenotyping and assessment of drug-drug interaction potential. Although extrapolation of the in vitro intrinsic clearance of a glucuronidated drug often underpredicts in vivo clearance, careful selection of in vitro experimental conditions and inclusion of extrahepatic glucuronidation may improve the predictivity of in vitro-in vivo extrapolation. Physiologically based pharmacokinetic (PBPK) modeling has also shown to be of value for predicting PK of drugs eliminated by glucuronidation.

Sedative and Analgesic Pharmacokinetics During Pediatric ECMO

Sedatives and analgesics are often administered to critically ill children supported by extracorporeal membrane oxygenation (ECMO) to facilitate comfort and to decrease risks of life-threatening complications. Optimization of sedative and analgesic dosing is necessary to achieve desired therapeutic benefits and must consider interactions between the circuit and patient that may affect drug metabolism, clearance, and impact on target organs. This paper reviews existing in vitro and pediatric in vivo literature concerning the effects of the ECMO circuit on sedative and analgesic disposition and offers dosing guidance for the management of critically ill children receiving these drugs.

Morphine Dose Optimization in Critically Ill Pediatric Patients With Acute Respiratory Failure: A Population Pharmacokinetic-Pharmacogenomic Study

OBJECTIVE

To develop a pharmacokinetic-pharmacogenomic population model of morphine in critically ill children with acute respiratory failure.

DESIGN

Prospective pharmacokinetic-pharmacogenomic observational study.

SETTING

Thirteen PICUs across the United States.

PATIENTS

Pediatric subjects (n = 66) mechanically ventilated for acute respiratory failure, weight greater than or equal to 7 kg, receiving morphine and/or midazolam continuous infusions.

INTERVENTIONS

Serial blood sampling for drug quantification and a single blood collection for genomic evaluation.

MEASUREMENTS AND MAIN RESULTS

Concentrations of morphine, the two main metabolites, morphine-3-glucuronide and morphine-6-glucuronide, were quantified by high-performance liquid chromatography tandem mass spectrometry/mass spectroscopy. Subjects were genotyped using the Illumina HumanOmniExpress genome-wide single nucleotide polymorphism chip. Nonlinear mixed-effects modeling was performed to develop the pharmacokinetic-pharmacogenomic model. A two-compartment model with linear elimination and two individual compartments for metabolites best describe morphine disposition in this population. Our analysis demonstrates that body weight and postmenstrual age are relevant predictors of pharmacokinetic parameters of morphine and its metabolites. Furthermore, our research shows that a duration of mechanical ventilation greater than or equal to 10 days reduces metabolite formation and elimination upwards of 30%. However, due to the small sample size and relative heterogeneity of the population, no heritable factors associated with uridine diphosphate glucuronyl transferase 2B7 metabolism of morphine were identified.

CONCLUSIONS

The results provide a better understanding of the disposition of morphine and its metabolites in critically ill children with acute respiratory failure requiring mechanical ventilation due to nonheritable factors. It also provides the groundwork for developing additional studies to investigate the role of heritable factors.

Midazolam Dose Optimization in Critically Ill Pediatric Patients With Acute Respiratory Failure: A Population Pharmacokinetic-Pharmacogenomic Study

OBJECTIVES

To develop a pharmacokinetic-pharmacogenomic population model of midazolam in critically ill children with primary respiratory failure.

DESIGN

Prospective pharmacokinetic-pharmacogenomic observational study.

SETTING

Thirteen PICUs across the United States.

PATIENTS

Pediatric subjects mechanically ventilated for acute respiratory failure, weight greater than or equal to 7 kg, receiving morphine and/or midazolam continuous infusions.

INTERVENTIONS

Serial blood sampling for drug quantification and a single blood collection for genomic evaluation.

MEASUREMENTS AND MAIN RESULTS

Concentrations of midazolam, the 1' (1`-hydroxymidazolam metabolite) and 4' (4`-hydroxymidazolam metabolite) hydroxyl, and the 1' and 4' glucuronide metabolites were measured. Subjects were genotyped using the Illumina HumanOmniExpress genome-wide single nucleotide polymorphism chip. Nonlinear mixed effects modeling was performed to develop the pharmacokinetic-pharmacogenomic model. Body weight, age, hepatic and renal functions, and the UGT2B7 rs62298861 polymorphism are relevant predictors of midazolam pharmacokinetic variables. The estimated midazolam clearance was 0.61 L/min/70kg. Time to reach 50% complete mature midazolam and 1`-hydroxymidazolam metabolite/4`-hydroxymidazolam metabolite clearances was 1.0 and 0.97 years postmenstrual age. The final model suggested a decrease in midazolam clearance with increase in alanine transaminase and a lower clearance of the glucuronide metabolites with a renal dysfunction. In the pharmacogenomic analysis, rs62298861 and rs28365062 in the UGT2B7 gene were in high linkage disequilibrium. Minor alleles were associated with a higher 1`-hydroxymidazolam metabolite clearance in Caucasians. In the pharmacokinetic-pharmacogenomic model, clearance was expected to increase by 10% in heterozygous and 20% in homozygous for the minor allele with respect to homozygous for the major allele.

CONCLUSIONS

This work leveraged available knowledge on nonheritable and heritable factors affecting midazolam pharmacokinetic in pediatric subjects with primary respiratory failure requiring mechanical ventilation, providing the basis for a future implementation of an individual-based approach to sedation.

Influence of Morbid Obesity on the Pharmacokinetics of Morphine, Morphine-3-Glucuronide, and Morphine-6-Glucuronide

Introduction

Obesity is associated with many pathophysiological changes that may result in altered drug metabolism. The aim of this study is to investigate the influence of obesity on the pharmacokinetics of morphine, morphine-3-glucuronide (M3G), and morphine-6-glucuronide (M6G) through a combined analysis in morbidly obese patients and non-obese healthy volunteers.

Methods

In this analysis, data from 20 morbidly obese patients [mean body mass index 49.9 kg/m² (range 37.6–78.6 kg/m²) and weight 151.3 kg (range 112–251.9 kg)] and 20 healthy volunteers [mean weight 70.6 kg (range 58–85 kg)] were included. Morbidly obese patients received 10 mg of intravenous (I.V.) morphine after gastric bypass surgery, with additional morphine I.V. doses as needed. Healthy volunteers received an I.V. bolus of morphine of 0.1 mg/kg followed by an infusion of 0.030 mg kg⁻¹ h⁻¹ for 1 h. Population pharmacokinetic modeling was performed using NONMEM 7.2.

Results

In morbidly obese patients, elimination clearance of M3G and M6G was decreased substantially compared with healthy volunteers (p < 0.001). Regarding glucuronidation, only a slight decrease in the formation of M6G and a delay in the formation of M3G was found (both p < 0.001). Obesity was also identified as a covariate for the peripheral volume of distribution of morphine (p < 0.001).

Conclusion

Metabolism of morphine is not altered in morbidly obese patients. However, decreased elimination of both M3G and M6G is evident, resulting in a substantial increase in exposure to these two metabolites. A rational explanation of this finding is that it results from alterations in membrane transporter function and/or expression in the liver.

ClinicalTrials.gov identifier: NCT01097148.

Electronic supplementary material

The online version of this article (doi:10.1007/s40262-017-0544-2) contains supplementary material, which is available to authorized users.

Endogenous Opiates and Behavior: 2016

This paper is the thirty-ninth consecutive installment of the annual review of research concerning the endogenous opioid system. It summarizes papers published during 2016 that studied the behavioral effects of molecular, pharmacological and genetic manipulation of opioid peptides, opioid receptors, opioid agonists and opioid antagonists. The particular topics that continue to be covered include the molecular-biochemical effects and neurochemical localization studies of endogenous opioids and their receptors related to behavior, and the roles of these opioid peptides and receptors in pain and analgesia, stress and social status, tolerance and dependence, learning and memory, eating and drinking, drug abuse and alcohol, sexual activity and hormones, pregnancy, development and endocrinology, mental illness and mood, seizures and neurologic disorders, electrical-related activity and neurophysiology, general activity and locomotion, gastrointestinal, renal and hepatic functions, cardiovascular responses, respiration and thermoregulation, and immunological responses.

Paracetamol and morphine for infant and neonatal pain; still a long way to go?

INTRODUCTION

Pharmacologic pain management in newborns and infants is often based on limited scientific data. To close the knowledge gap, drug-related research in this population is increasingly supported by the authorities, but remains very challenging. This review summarizes the challenges of analgesic studies in newborns and infants on morphine and paracetamol (acetaminophen). Areas covered: Aspects such as the definition and multimodal character of pain are reflected to newborn infants. Specific problems addressed include defining pharmacodynamic endpoints, performing clinical trials in this population and assessing developmental changes in both pharmacokinetics and pharmacodynamics. Expert commentary: Neonatal and infant pain management research faces two major challenges: lack of clear biomarkers and very heterogeneous pharmacokinetics and pharmacodynamics of analgesics. There is a clear call for integral research addressing the multimodality of pain in this population and further developing population pharmacokinetic models towards physiology-based models.

Pharmacodynamics and Pharmacokinetics of Morphine After Cardiac Surgery in Children With and Without Down Syndrome

OBJECTIVE

To compare the pharmacodynamics and pharmacokinetics of IV morphine after cardiac surgery in two groups of children-those with and without Down syndrome.

DESIGN

Prospective, single-center observational trial.

SETTING

PICU in a university-affiliated pediatric teaching hospital.

PATIENTS

Twenty-one children with Down syndrome and 17 without, 3-36 months old, scheduled for cardiac surgery with cardiopulmonary bypass.

INTERVENTIONS

A loading dose of morphine (100 μg/kg) was administered after coming off bypass; thereafter, morphine infusion was commenced at 40 μg/kg/hr. During intensive care, nurses regularly assessed pain and discomfort with validated observational instruments (COMFORT-Behavior scale and Numeric Rating Scale-for pain). These scores guided analgesic and sedative treatment. Plasma samples were obtained for pharmacokinetic analysis.

MEASUREMENTS AND MAIN RESULTS

Median COMFORT-Behavior and Numeric Rating Scale scores were not statistically significantly different between the two groups. The median morphine infusion rate during the first 24 hours after surgery was 31.3 μg/kg/hr (interquartile range, 23.4-36.4) in the Down syndrome group versus 31.7 μg/kg/hr (interquartile range, 25.1-36.1) in the control group (p = 1.00). Population pharmacokinetic analysis revealed no statistically significant differences in any of the pharmacokinetic variables of morphine between the children with and without Down syndrome.

CONCLUSIONS

This prospective trial showed that there are no differences in pharmacokinetics or pharmacodynamics between children with and without Down syndrome if pain and distress management is titrated to effect based on outcomes of validated assessment instruments. We have no evidence to adjust morphine dosing after cardiac surgery in children with Down syndrome.

Children in clinical trials: towards evidence-based pediatric pharmacotherapy using pharmacokinetic-pharmacodynamic modeling

INTRODUCTION

In pediatric pharmacotherapy, many drugs are still used off-label, and their efficacy and safety is not well characterized. Different efficacy and safety profiles in children of varying ages may be anticipated, due to developmental changes occurring across pediatric life.

AREAS COVERED

Beside pharmacokinetic (PK) studies, pharmacodynamic (PD) studies are urgently needed. Validated PKPD models can be used to derive optimal dosing regimens for children of different ages, which can be evaluated in a prospective study before implementation in clinical practice. Strategies should be developed to ensure that formularies update their drug dosing guidelines regularly according to the most recent advances in research, allowing for clinicians to integrate these guidelines in daily practice. Expert commentary: We anticipate a trend towards a systems-level approach in pediatric modeling to optimally use the information gained in pediatric trials. For this approach, properly designed clinical PKPD studies will remain the backbone of pediatric research.