Skip to main content
Advertisement

Main menu

  • Home
  • Articles
    • Current Issue
    • Fast Forward
    • Latest Articles
    • Archive
  • Information
    • Instructions to Authors
    • Submit a Manuscript
    • FAQs
    • For Subscribers
    • Terms & Conditions of Use
    • Permissions
  • Editorial Board
  • Alerts
    • Alerts
    • RSS Feeds
  • Virtual Issues
  • Feedback
  • Other Publications
    • Drug Metabolism and Disposition
    • Journal of Pharmacology and Experimental Therapeutics
    • Molecular Pharmacology
    • Pharmacological Reviews
    • Pharmacology Research & Perspectives
    • ASPET

User menu

  • My alerts
  • Log in
  • My Cart

Search

  • Advanced search
Drug Metabolism & Disposition
  • Other Publications
    • Drug Metabolism and Disposition
    • Journal of Pharmacology and Experimental Therapeutics
    • Molecular Pharmacology
    • Pharmacological Reviews
    • Pharmacology Research & Perspectives
    • ASPET
  • My alerts
  • Log in
  • My Cart
Drug Metabolism & Disposition

Advanced Search

  • Home
  • Articles
    • Current Issue
    • Fast Forward
    • Latest Articles
    • Archive
  • Information
    • Instructions to Authors
    • Submit a Manuscript
    • FAQs
    • For Subscribers
    • Terms & Conditions of Use
    • Permissions
  • Editorial Board
  • Alerts
    • Alerts
    • RSS Feeds
  • Virtual Issues
  • Feedback
  • Visit dmd on Facebook
  • Follow dmd on Twitter
  • Follow ASPET on LinkedIn
Research ArticleSpecial Section on Pediatric Drug Disposition and Pharmacokinetics—Minireview

Physiologically Based Pharmacokinetic Modeling in Pediatric Oncology Drug Development

Nathalie Rioux and Nigel J. Waters
Drug Metabolism and Disposition July 2016, 44 (7) 934-943; DOI: https://doi.org/10.1124/dmd.115.068031
Nathalie Rioux
Epizyme Inc., Cambridge, Massachusetts
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Nigel J. Waters
Epizyme Inc., Cambridge, Massachusetts
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

Abstract

Childhood cancer represents more than 100 rare and ultra-rare diseases, with an estimated 12,400 new cases diagnosed each year in the United States. As such, this much smaller patient population has led to pediatric oncology drug development lagging behind that for adult cancers. Developing drugs for pediatric malignancies also brings with it a number of unique trial design considerations, including flexible enrollment approaches, age-appropriate formulation, acceptable sampling schedules, and balancing the need for age-stratified dosing regimens, given the smaller patient populations. The regulatory landscape for pediatric pharmacotherapy has evolved with U.S. Food and Drug Administration (FDA) legislation such as the 2012 FDA Safety and Innovation Act. In parallel, regulatory authorities have recommended the application of physiologically based pharmacokinetic (PBPK) modeling, for example, in the recently issued FDA Strategic Plan for Accelerating the Development of Therapies for Pediatric Rare Diseases. PBPK modeling provides a quantitative and systems-based framework that allows the effects of intrinsic and extrinsic factors on drug exposure to be modeled in a mechanistic fashion. The application of PBPK modeling in drug development for pediatric cancers is relatively nascent, with several retrospective analyses of cytotoxic therapies, and latterly for targeted agents such as obatoclax and imatinib. More recently, we have employed PBPK modeling in a prospective manner to inform the first pediatric trials of pinometostat and tazemetostat in genetically defined populations (mixed lineage leukemia–rearranged and integrase interactor-1–deficient sarcomas, respectively). In this review, we evaluate the application of PBPK modeling in pediatric cancer drug development and discuss the important challenges that lie ahead in this field.

Introduction

Childhood cancer represents a collection of more than 100 rare and ultra-rare diseases, with an estimated 12,400 new cases diagnosed each year in the United States in patients aged <21 years (Ries et al., 1999). This much smaller patient population has, in large part, led to cancer drug development for pediatric patients being an afterthought, with approval for pediatric indications following the development and approval of the agent in adult cancers (Kearns et al., 2003). The rate of development has also been slow; from 1948 to 2003, there were 120 new cancer drug approvals, only 30 of which were used in pediatrics (Adamson et al., 2014). Developing drugs for pediatric malignancies also brings with it a number of unique challenges in terms of trial design, potential for disparity in genetics and pathophysiology relative to the adult disease, and what is considered acceptable clinical benefit. Trial design considerations include age-stratified dosing regimens, age-appropriate formulations (Breitkreutz and Boos, 2007), amenable sampling schemes (given the limitations on the number and volume of blood draws), as well as the necessity for flexible enrollment approaches in phase I, such as the rolling six design (Skolnik et al., 2008). The specific oncogenic pathways in adult carcinomas may not be active in the childhood malignancy, although common cellular pathways have emerged (Rossig et al., 2011). The benchmark for clinical benefit can also be quite different. Increasing overall survival by a few months may be acceptable in adults but the focus in pediatric oncology is cure, with cure rates in most cases of >70% and 5-year event free survival of approximately 80% (Adamson et al., 2014). However, novel drug approaches are needed for children with difficult-to-treat malignancies such as stage IV neuroblastoma, sarcomas, brain tumors, and relapsed leukemia (Horton and Berg, 2011). As a direct consequence of the unmet need in pediatric pharmacotherapy in general, the regulatory landscape has evolved with the introduction of U.S. Food and Drug Administration (FDA) legislation including the Best Pharmaceuticals for Children Act in 2002 and the Pediatric Research Equity Act in 2003 providing incentives for evaluating and developing drugs for children, which were both made permanent under the FDA Safety and Innovation Act in 2012.

In developing and optimizing pediatric pharmacotherapy, an appreciation of the comparative physiology and biochemistry between adults and children is paramount, being particularly acute in prenatal, infant, and toddler populations. As the long-standing adage states, “children are not small adults”; thus, simple dose reductions may not be sufficient and the nature of disease can also be divergent. As our knowledge of normal growth and development has increased, so too has our understanding of how these growth changes can greatly affect pharmacokinetics (PK) and ultimately response to therapy. Drugs in pediatric oncology are typically administered by either intravenous or oral routes; as such, there are a number of physical, chemical, and biologic barriers that must be overcome to reach the target site of action. The gastrointestinal absorption of drugs can be markedly affected by the developmental changes in absorptive surface area, intraluminal pH, splanchnic blood flow, gastric emptying time, and intestinal motility (Batchelor et al., 2014). The maturation profile in biliary function also leads to age-related differences in duodenal bile salt concentrations, which can affect drug dissolution and solubility at the site of absorption (Poley et al., 1964). The ontogeny in intestinal expression of drug-metabolizing enzymes and transporters can also be a contributing factor for oral bioavailability. This is illustrated with the alkylating antineoplastic agent, busulfan, in which biopsies of distal duodenum indicated that glutathione-S-transferase (GST) activity decreased from infancy through adolescence and manifested as a reduced oral clearance of busulfan, a GST substrate (Gibbs et al., 1999). Drug distribution can also show marked age-related changes as a result of physiologic and biochemical ontogeny. Infants and children have proportionally higher extracellular water and total body water than adults, in addition to the obvious differences in cardiac output, tissue perfusion, and organ size (Kearns et al., 2003). The major plasma proteins, albumin and α-1-acid glycoprotein (AAG), also show developmental changes, which can affect the distribution of highly protein bound drugs (Ehrnebo et al., 1971). The delayed maturation of drug-metabolizing enzymes (DMEs) and potential for toxicity was typified with the case of newborns treated with doses of the antimicrobial chloramphenicol, extrapolated from those shown to be safe and effective in adults, which led to cardiovascular collapse associated with Gray syndrome in neonates. The mechanistic basis for these drug-induced sequelae was later shown to be an immature uridine diphosphoglucuronosyl transferase (UGT) system resulting in impaired metabolic clearance. The maturation profile of DMEs, the most studied in this regard being the cytochrome P450 (P450) superfamily (Hines, 2009), can follow one of three possible trajectories: 1) the enzyme exhibits high expression during the first trimester and either remains high or decreases during gestation before a steep decline in expression level postnatally, within a few days or up to 2 years after birth (e.g., CYP3A7); 2) the enzyme shows relatively constant expression throughout gestation and postnatally (e.g., CYP2C19); and 3) the enzyme may be not expressed or present at low levels in the fetus, followed by a marked increase in expression level, depending on the isoform in question, between the second and third trimesters to the first years of life (e.g., CYP3A4). Similarly, isoforms of the other DME families can be categorized in this manner, such as the sulfotransferases, flavin monooxygenases, and UGTs. The prevalence of enzymes falling into the third class underpins the consistent observation in clinical studies of drugs metabolized by the liver showing an age-dependent increase in plasma clearance in children aged <10 years, such as dextromethorphan (Blake et al., 2007), caffeine (Pons et al., 1988), and midazolam (Altamimi et al., 2015). It should also be noted that excretory mechanisms also show an age-based relationship in terms of renal function, glomerular filtration rate, and renal blood flow that can all affect the elimination of renally excreted drugs. More detailed reviews of the physiologic and biochemical ontogeny influencing the absorption, distribution, metabolism, and excretion (ADME) of drugs were reported previously (Kearns et al., 2003; Hines, 2009).

It is clear that ontogenesis plays a key role in the metabolism and disposition of drugs in children; as a result, simplified dosing algorithms based on body size alone are usually insufficient. Traditional methods such as allometry have been used to scale pediatric doses, wherein absolute dose (or clearance) is related to the child-to-adult body weight ratio raised to the power of 3/4 (eq. 1) (Björkman, 2006; Mahmood, 2006). The basis for the 3/4 power law has been discussed elsewhere (West et al., 1999). Since this does not account for maturation processes and has a tendency to overpredict clearance in younger children, its use has been limited to children aged >1 to 2 years. Further refinement with the inclusion of terms for maturation or organ function has shown utility within predefined limits (eq. 2) but still lacks the holistic and mechanistic facets that are possible with physiologically based pharmacokinetic (PBPK) modeling.Embedded Image(1)Embedded Image(2)In eqs. 1 and 2, BW indicates body weight, CL is clearance, MF is maturation function, and OF is organ function.

PBPK modeling provides a quantitative and systems-based framework with each compartment representing a physiologic volume interconnected by flow rates that are anatomically representative of the circulatory system (Fig. 1). Mass balance equations describe the transfer of drug from the arterial blood into tissues and from tissues into venous blood. As such, simulations are dynamic in nature, with model fitting and prediction focused on time-concentration data. Importantly, PBPK models are more comprehensive than empirical PK models, incorporating both system-specific parameters as well as drug-specific parameters. Although it is not a new approach, PBPK modeling applied in pharmaceutical research and development has gained increasing attention in recent years with the availability of robust in vitro and in silico data, advancements in in vitro to in vivo extrapolation, and the substantial investment that has been made in the curated databases of system-dependent and probe drug-dependent parameters that are the foundation of the commercially available software packages. This technological progress has enhanced the number and breadth of robust applications in the area of clinical pharmacology, including first-in-human studies, special populations, drug interactions, and biopharmaceutics/formulations (Jones et al., 2015; Sager et al., 2015). Considerable experience has now accumulated with the development of PBPK models to describe drug concentration-time profiles in adults (Jones et al., 2009; Kostewicz et al., 2014), and, more recently, in cancer drug development (Block, 2015) and special populations such as pediatrics (Khalil and Läer, 2011; Barrett et al., 2012; Maharaj and Edginton, 2014). There have been multiple applications of PBPK in pediatric drug development, many of which are in support of the first pediatric trial, including starting dose selection, prediction of exposures across the age continuum, optimization of blood sampling strategy, prediction of target organ exposure [safety and pharmacodynamics (PD)], and evaluation of drug interaction potential.

Fig. 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 1.

(A) A typical workflow scheme for the development of a PBPK model in support of pediatric dose selection early in clinical development. (B) Representative drug-specific data inputs used in “bottom-up” PBPK model building. (C) Generalized scheme for a PBPK model. IV, intravenous; PO, oral.

The momentum in the PBPK modeling arena has evolved to a regulatory science. As mentioned earlier, the health authorities, FDA, European Medicines Agency, and the Ministry of Health, Labor and Welfare of Japan (European Medicines Agency, 2014; Zhao et al., 2012; Shepard et al., 2015; Wagner et al., 2015) have all highlighted PBPK modeling in guidances covering drug interactions and hepatic impairment. Furthermore, PBPK modeling and simulation has recently been used to directly support labeling statements (e.g., Imbruvica; related to the use of moderate CYP3A inhibitors with ibrutinib; Imbruvica product label; Janssen, Titusville, NJ). Between 2008 and 2012, the FDA received 33 Investigational New Drug/New Drug Application submissions containing PBPK modeling approaches (including six pediatric submissions, although none in oncology). In parallel with the increasing number of submissions, the agency has increasingly used de novo (i.e., FDA initiated) PBPK modeling in its reviews to help characterize PK in a variety of complex clinical scenarios (Huang et al., 2013). In 2013, the cumulative number of submissions containing PBPK modeling increased to 84, with 22% of them related to pediatric applications (Zhao, 2014). Furthermore, the FDA recently issued a “Strategic Plan for Accelerating the Development of Therapies for Pediatric Rare Diseases,” which specifically highlights the role of modeling and simulation approaches such as PBPK modeling to inform the design and conduct of PK/PD studies and other clinical trials for investigational drugs in pediatric rare disease populations (http://www.fda.gov). Given this regulatory focus, the PBPK community is redoubling efforts to formulate industry-wide best practices and standardize the level of rigor and verification that is needed in PBPK model submissions (Sager et al., 2015).

The construction of PBPK models has typically used the various in vitro and physicochemical data available in early preclinical drug development in what is often termed a “bottom-up” approach, in an attempt to recapitulate the concentration-time data (Fig. 1). Alternatively, there may be reasons to consider a “top-down” approach in which the underlying drug-specific parameterization of the PBPK model is accomplished by optimizing the fit of the time-concentration data in question. In practice, there is usually an approach that lands somewhere between the two, termed “middle out,” in which some of the drug-specific parameters based on in vitro or in silico data may not scale appropriately or there are missing data and incomplete information regarding some aspects of drug disposition, and an element of top-down model optimization is necessary. Furthermore, mechanistic understanding can be gained through parameter sensitivity analysis and asking “what if” questions of the optimized model.

Case Studies of PBPK Modeling in Pediatric Oncology

The application of PBPK modeling in pediatric oncology drug development is still quite nascent and has largely been focused on retrospective analyses in validating the approach. A summary of reported applications is shown in Table 1. The first published model in this area of drug research dates back to 1982, with this seminal PBPK model being built with a limited number of compartments from the Bischoff–Dedrick multiorgan model. The observed parent cisplatin and total platinum serum concentrations in 14 pediatric patients were adequately recapitulated in the model and helped to better understand changes in the renal clearance of total platinum (Evans et al., 1982). Similar reports based on a modified Bischoff–Dedrick model followed and were expanded to incorporate effusion spaces as a compartment of drug distribution and clearance (Li and Gwilt, 2002). Simulated methotrexate plasma concentrations in adults without effusions were adequately modeled prior to providing congruent simulated and measured plasma and effusion methotrexate concentrations in one pediatric patient with a malignant pleural effusion after intravenous administration. This early work has been expanded on with PBPK modeling applications reported for a number of drugs used in the pediatric oncology setting, including busulfan, docetaxel, etoposide, coadministered 6-mercaptopurine and methotrexate, and, more recently, imatinib, obatoclax, pinometostat (EPZ-5676; 9-[5-deoxy-5-[[cis-3-[2-[6-(1,1-dimethylethyl)-1H-benzimidazol-2-yl]ethyl]cyclobutyl](1-methylethyl)amino]-β-d-ribofuranosyl]-9H-purin-6-amine), and tazemetostat (EPZ-6438; N-[(1,2-dihydro-4,6-dimethyl-2-oxo-3-pyridinyl)methyl]-5-[ethyl(tetrahydro-2H-pyran-4-yl)amino]-4-methyl-4′-(4-morpholinylmethyl)-[1,1′-biphenyl]-3-carboxamide).

View this table:
  • View inline
  • View popup
TABLE 1

Summary of PBPK modeling applications in pediatric oncology drug development

Unlike most pediatric applications wherein adult PK data are leveraged in the PBPK model building process, the cisplatin model was initially developed with PK data in dogs and adjusted for human pediatric physiology (Evans et al., 1982). As summarized in Table 1, the more common approach to pediatric PBPK modeling is to scale a validated adult PBPK model to children by incorporating ontogeny and maturation processes, thereby negating the issue of potential species differences. Another exception to the adult-to-children model workflow was presented by Barrett et al. (2011). After establishing quantitative pharmacology relationships between in vitro activity data and a preclinical diseased mouse model, obatoclax exposures were recapitulated in the mouse prior to scaling to a human infant, with the goal of evaluating the potential to achieve targeted exposure in this population (Barrett et al., 2011). This approach allowed tissue distribution to be modeled and scaled, such that the PBPK model based on obatoclax physicochemical properties compared well with observed obatoclax levels in the plasma, spleen, liver, and kidney of leukemia-bearing mice for two single intravenous dose levels; by contrast, exposure values were overpredicted in the brain. After scaling to human infants, PBPK simulation of obatoclax administration to 1-year-old patients suggested that adequate drug exposure to target organs should be achievable at clinically relevant doses. To the best of our knowledge, this alternate approach for obatoclax has not been verified with clinical data from a pediatric population. By contrast, and in line with a recently proposed pediatric PBPK model development workflow (Maharaj and Edginton, 2014), Thai et al. (2015) first developed a model for docetaxel from in vitro ADME and intravenous PK data from more than 500 patients with cancer after either single or multiple dosing, to demonstrate how modeling can be applied to optimize dose and sampling times for a pediatric PK bridging study. Docetaxel is highly protein bound, has a high volume of distribution, and is highly metabolized. In the full-body PBPK model developed for adult patients, all 16 compartments were well stirred or perfusion limited, except for the liver, where organic anion transporting polypeptide 1B1 and1B3 activities were considered, and, for the muscle, where apparent passive diffusion was used. The tissue/plasma partition was determined using the volume of distribution as a primary predictor, in addition to physicochemical descriptors. Hepatic (CYP3A), biliary, and renal clearances were fit in a top-down approach. The model was customized for oncology patients by considering demographic and pathophysiological differences from healthy subjects. After scaling to a pediatric population, and as determined by visual inspection, the model adequately predicted docetaxel plasma concentration-time profiles in neonates to 18-year-old patients by accounting for age-dependent physiologic differences and CYP3A ontogeny, with predicted clearance and volume of distribution within 1.5-fold of observed data (Thai et al., 2015).

Although P450 ontogeny functions are implemented in commercial PBPK software, scaling of other metabolic pathways is less common. In that regard, the modeling exercise described by Diestelhorst et al. (2014) for the DNA-alkylating agent busulfan provides a proof of concept for GST substrates. In brief, an adult PBPK model was refined by implementing GST-A1 in 11 organs, using the PK-Sim integrated enzyme expression database, and adding irreversible DNA binding and plasma protein binding processes. Age-dependent enzyme activity and maturation factors were considered and the adult-to-child scaling indicated lower clearance values for children relative to adults. Intravenous administration of busulfan was simulated in pediatric patients, with a mean percentage error for all patients of 3.9%; 3 of 23 children demonstrated a mean percentage error of greater than ±30%, showing an adequate predictive performance of this retrospective model (Diestelhorst et al., 2014). In another study, both CYP3A4 and UGT1A1 ontogeny information was included to model etoposide exposure in children (Kersting et al., 2012). Etoposide is highly bound to plasma proteins, metabolized by CYP3A4 and UGT1A1, and eliminated in urine and bile. An adult PBPK model was developed using data from nine women with primary breast cancer receiving the drug intravenously, as part of a polychemotherapy regimen before stem cell transplantation. The physicochemical parameters of etoposide, along with the unbound plasma fraction at high and low doses in adults, kinetics for in vitro metabolic enzymes [CYP3A4 and UGT1A1], and active biliary and renal transporters [P-glycoprotein (P-gp), multidrug resistance-associated protein 2, and a hypothetical influx transporter] were incorporated into the adult model. The tissue/plasma partition coefficients were generated using the Rodgers and Rowland model. In adults, the simulated plasma concentration-time profiles of protein-bound and free etoposide were in good agreement with observed data, with mean relative deviation of 1.12 and 1.36 for low and high doses of etoposide, respectively. The model was scaled to children, incorporating ontogeny for both CYP3A4 and UGT1A1, glomerular filtration, and tubular excretion mediated by P-gp and adequately simulated the exposure seen in 18 children with various diagnoses and normal renal and liver function. Of interest, the effects of the coadministration of cyclosporine A on the metabolism and excretion of etoposide was determined in five patients, suggesting that the PBPK model could be useful for performing hypothesis testing on the effect of concomitant medications, a relatively poorly explored application of PBPK in pediatric drug development but extremely relevant given the polypharmacy in this patient population. A recent model developed for 6-mercaptopurine, a purine antimetabolite and prodrug used in combination with methotrexate as therapy for childhood acute lymphoblastic leukemia, is another example of drug–drug interaction (DDI) risk assessment by PBPK modeling in pediatric oncology (Ogungbenro et al., 2014a,b). 6-Mercaptopurine undergoes very extensive intestinal and hepatic metabolism after oral dosing, owing to the activity of xanthine oxidase leading to low and highly variable bioavailability. Dose adjustment during treatment is still based on toxicity rather than routine therapeutic drug monitoring and this work was an attempt to improve dose individualization and dosage regimen optimization through modeling and simulation, ultimately to achieve a better outcome in patients with childhood acute lymphoblastic leukemia. The PBPK model, based on the assumption of the same elimination pathways in adults and children with age-dependent parameters adjusted for pediatric scaling, adequately predicted plasma and red blood cell concentrations both in terms of population mean and variability versus observed data after intravenous or oral administration in children aged 3–18 years. Further model refinements incorporated parameters such as net secretion clearance, biliary transit time, and red blood cell distribution and binding and enabled the oral absorption of methotrexate to be well described with recapitulation of the nonlinear relationship between the fraction absorbed and the methotrexate dose. The inhibition of 6-mercaptopurine first-pass metabolism by methotrexate (via xanthine oxidase inhibition in the gut and liver) predicted an interaction but to a greater extent than that reported clinically. The predicted percentage increase in the area under the curve (AUC) and peak serum concentration (Cmax) was approximately 65% and 50% in 5- to 18-year-olds, respectively, whereas the observed increase in AUC and Cmax has been reported as 31% and 26%, respectively.

In the context of a Pediatric Investigation Plan, the European Medicines Agency (2013) requested PBPK simulations as part of the submission package for the use of imatinib (Gleevec; Novartis, Basel, Switzerland) in the treatment of pediatric patients with diagnosed Philadelphia chromosome–positive acute lymphoblastic leukemia combined with chemotherapy. PBPK modeling was used to predict area under the curve at steady state (AUCss) and concentration-time profiles in pediatric subjects, to evaluate factors influencing imatinib exposure in pediatric patients. An adult PBPK model was validated and then transformed using growth and maturation processes for pediatric scaling. The model evaluations compared the predicted versus observed imatinib AUCss and concentration-time profile from 67 pediatric patients. The simulations showed that 29 of 31 actual AUCss values normalized to 340 mg/m2 in pediatric patients fell within the 0.5 and 99.5 percentiles of the model-projected range scaled from adult measurements and the predicted plasma concentration-time profiles were generally in good agreement for the pediatric cohort (n = 67), except for subjects aged ≤2 years, for which the exposure appeared to be overpredicted. Nonetheless, the prediction was 1.5-fold of the adult value at 1 year of age. Of interest, the model was used to better understand the differences in prediction of children and adults, which seems to be the mixed effect of changing distribution volume and blood flow, in addition to clearance maturation with age. Overall, in conjunction with a population PK analysis, PBPK modeling suggested that the proposed posology of 340 mg/m2 in pediatric patients from 1 to 18 years of age was appropriate for imatinib.

Pediatric PBPK Modeling In-House at Epizyme

Building on many of the retrospective analyses described above and the recent regulatory precedent, we have applied PBPK modeling in a prospective manner to enable the first trials in pediatric patients for two first-in-class epigenetic therapies: pinometostat in mixed lineage leukemia–rearranged (MLL-r) and tazemetostat in integrase interactor (INI)-1–deficient solid tumors. The preemptive utilization of this approach early in development has informed some important aspects of phase I pediatric trial design, including starting dose regimen and blood sampling strategy. With model verification and optimization on availability of PK data in pediatric patients, the PBPK modeling workflow can be further applied to address additional clinical pharmacology questions such as drug interactions and PK/PD analysis, through the course of development of these investigational therapies in pediatric populations.

Relapsed/refractory MLL-r acute leukemia in children has a poor prognosis, with a reported 5-year overall survival for children aged <1 year with recurrent MLL-r leukemia being 4%–20%. The mixed lineage leukemia (MLL) gene, on chromosome 11q23, codes for a histone methyltransferase (HMT) that is responsible for methylation of H3K4, a modification associated with active transcription. Translocations of MLL result in the loss of the SET or catalytic domain of the protein with the most common translocation partners, AF4, AF9, and ENL, recruiting another HMT, the disruptor of telomeric silencing 1-like (DOT1L). The aberrant recruitment of DOT1L to MLL fusion target genes results in ectopic H3K79 methylation and increased expression of genes including HOXA9 and MEIS1, which are involved in leukemogenesis of MLL-r leukemias. Pinometostat (EPZ-5676) is a small molecule inhibitor of DOT1L with subnanomolar affinity and > 37,000-fold selectively against other HMTs. Preclinically, pinometostat selectively inhibits intracellular histone H3K79 methylation and downstream target gene expression and demonstrated complete tumor regressions in an MLL-r leukemia xenograft model (Daigle et al., 2013). Because of the unmet need of MLL-r in children, we employed a prospective PBPK modeling approach leveraging pinometostat preclinical data (Basavapathruni et al., 2014) and early clinical data from the first-in-human phase I open label study in adult patients with relapsed/refractory leukemia (Stein et al., 2014), to guide dose selection and trial design for a companion pediatric study (Waters et al., 2014; Shukla et al., 2015). A PBPK model describing the concentration-time profile of pinometostat after continuous intravenous administration in adult patients at dose levels of 24–90 mg/m2 per day was built using Simcyp (Simcyp Ltd., Sheffield, UK). The quantitative contribution of P450 enzymes to the metabolic clearance of pinometostat was based on intrinsic clearance (CLint) data derived from in vitro recombinant P450 systems using a relative activity factor approach. Pinometostat is a substrate for CYP3A4 and CYP2C19 in human and these enzymes are believed to account for 90% and 10% of the metabolic intrinsic clearance respectively. The metabolic intrinsic clearance was scaled using the well stirred model to generate a total hepatic in vivo CLint, which was then allocated to the P450 isoforms involved in the metabolism of pinometostat and thus allowed isoform-specific ontogeny functions to be used in the translation to the pediatric setting. The renal clearance in human was low (0.1 l/h) and was also incorporated into the model. Pinometostat was shown to preferentially bind to AAG in vitro (manuscript in preparation) and this allowed both the extent of binding and identity of plasma proteins involved to be included in the model, together with the maturation function for AAG. Early in model optimization, perfusion-limited kinetics were not able to recapitulate the plasma steady-state profile due to a marked overprediction of total clearance (several fold above the median observed clearance of 5.5 l/h). Using permeability-limited kinetics and a low steady-state volume of distribution (0.08 l/kg, derived from a separate compartmental analysis) was consistent with the short terminal half-life (t1/2) observed postinfusion and was able to fit the data well in terms of CL, steady-state plasma concentration, and the initial phase after the start of infusion. A sensitivity analysis of the passive diffusion clearance indicated that a value of 0.0075 mL/min/million hepatocytes gave the best model fit. This was plausible and consistent with the physicochemical properties and preclinical ADME data which showed low permeability in MDCK cell monolayers and a greater than 20-fold higher liver microsomal scaled CLint compared with hepatocyte scaled CLint (Basavapathruni et al., 2014). The adult model simulations and observed data are shown in Fig. 2A, with the Css predicted within 2-fold of the observed values and CL predicted within 20% of observed across the dose range 24–90 mg/m2 per day. Having built and qualified the PBPK model for pinometostat using adult PK data, the next step was to prospectively estimate exposures across the pediatric age range from 1 month to 18 years. The median clearance projected in pediatric patients ranged from 0.5 l/h in 1- to 3-month-olds to 4.2 l/h in 6- to 18-year-olds, and when normalized for body surface area (BSA) ranged from 1.8 l/h per square meter to 3.2 l/h per square meter, respectively. A modest 1.7-fold difference in BSA-normalized clearance between infant and adolescent suggested a dampening of the effect of P450 ontogeny on the clearance of pinometostat, as a direct result of age-independent, permeability-limited hepatic uptake incorporated into the model. This observation is in stark contrast to the clearance of prototypical CYP3A substrates in children which show an age-dependent relationship, concordant with the maturation of CYP3A expression during the first 2 years postnatal (Björkman, 2005; Hines, 2009). The model assumption with potential to affect predictive accuracy was clearly the permeability-limited metabolic clearance being independent of age. This was considered reasonable since the basis for the rate-limiting passive diffusion clearance of pinometostat was physicochemical in nature, as opposed to transporter-mediated active efflux, for example, which may show age dependent expression or activity. Simulations of the predicted steady-state systemic exposure of pinometostat in pediatric virtual populations were used to support starting dose selection. At a fixed BSA-normalized dose, adjustments of 55%–65% of the adult dose were projected in infants up to 2 years old to achieve equivalent exposures to adults; in children aged >6 years, the dose was predicted to be equivalent to the adult dose. For the practical purposes of trial conduct in this rare patient population, the starting dose level and age stratification was further simplified in a conservative manner to derive a pediatric starting dose of 80% and 50% of the highest adult dose (90 mg/m2 per day) in >12-month-olds and <12-month-olds, respectively. PK data recently obtained from this study (Shukla et al., 2015) have enabled model verification and this is summarized in Fig. 2, B and C. The observed Css at 70 mg/m2 per day (n = 6; age range, 1.25–15 years) was within 0.75- to 2.33-fold of the predicted Css across this age range, showing that the data were consistent with the postulated effect of dampened P450 ontogeny on the clearance of pinometostat. In this cohort, clearance ranged from 1.8 to 3.7 l/h per square meter, showing good concordance with the PBPK model CL projection of 2.4–3.2 l/h per square meter in subjects older than 1 year. The limited enrollment of patients less than 1 year of age in the relapsed/refractory setting precluded a more thorough assessment of P450 ontogeny on pinometostat clearance. Notwithstanding, the observed Css at 90 mg/m2/d (n = 5; age range, 1–15 years) showed a similar level of agreement with model predictions and was comparable with the plasma exposure observed in adult at this dose level (Fig. 2C; Shukla et al., 2015).

Fig. 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 2.

(A) Representative plot showing the model fit for interim adult PK data at 90 mg/m2 per day. Data at 24 mg/m2 per day (n = 5) and 36 mg/m2 per day (n = 3) showed similar fits. The black line represents the mean, the dotted line is the 95% confidence interval, and gray lines represent mean results from individual simulation trials (20 trials of four subjects). The majority of data fall within the 95% confidence interval of the model. (B) Representative plots of the simulated time-concentration profile for pinometostat administered at 70 mg/m2 per day in age groups ranging from 6 months to 18 years (solid gray lines), with observed plasma concentrations in pediatric patients receiving an infusion of 70 mg/m2 per day (open circles). Medians are shown as solid black symbols. (C) Box and whisker plot of the plasma Css of pinometostat administered at 45, 70, and 90 mg/m2 per day in pediatric patients (median age range, 0.33, 5.5, and 4 years, respectively). The dotted line corresponds to the observed median plasma Css in adults at 90 mg/m2 per day, and the shaded area represents the predicted Css range in pediatric patients aged >1 year at the dose level of 90 mg/m2 per day.

Tazemetostat (EPZ-6438) is a selective, orally active, small molecule inhibitor of the histone-lysine methyltransferase Enhancer of Zeste Homolog 2 (EZH2), which has been implicated in the pathogenesis of a variety of malignancies, including B-cell non-Hodgkin lymphoma and integrase interactor (INI)-1–deficient tumors (Knutson et al., 2013, 2014). INI-1–negative malignant rhabdoid tumors, for instance, occur mainly in young children. Other INI-1–negative tumors include epithelioid sarcomas, epithelioid malignant peripheral nerve sheath tumors, extraskeletal myxoid chondrosarcoma, renal medullary carcinoma, and myoepithelial carcinomas, characterized by their rarity and unmet medical need. INI-1–deficient tumors include several tumors for which there is no established standard of care, with a median survival of 0.3 years after recurrence in children. The primary objective of the PBPK modeling approach was to simulate tazemetostat exposures in adults prior to scaling to a pediatric population and to predict starting doses in children that would provide systemic exposures in the safe and efficacious range observed in adults. Tazemetostat PK data from patients enrolled in the dose escalation cohorts of the first-in-human phase I clinical study across the oral dose range of 100 mg (suspension) and 100, 200, 400, 800, and 1600 mg (tablet) twice daily, together with physicochemical and in vitro data including plasma protein binding, blood partitioning, metabolic stability, and P450 phenotyping, were used to simulate adult exposures by PBPK modeling using GastroPlus (version 8.5; Simulation Plus, Inc). In these adult patients, tazemetostat was rapidly absorbed with a maximum plasma concentration observed approximately 1 to 2 hours postdose. Plasma concentrations declined in a monoexponential manner with a mean t1/2 of approximately 3–6 hours. The tazemetostat AUC was slightly higher than dose proportional after a single dose, reasonably dose proportional at steady state, and comparable between the tablet and suspension formulations. After multiple dosing (day 15), the median time to reach maximum plasma concentrations (tmax) and t1/2 remained unchanged across the dose range, with a modest decrease in systemic exposure on day 15 and no further reduction onward (Ribrag et al., 2015). In the PBPK model, partition coefficients were set using the Lucakova (Rodgers–Rowland Single) equation for perfusion-limited tissues. Effective permeability was fitted to the clinical PK data in a top-down approach, refining the initial Cmax estimated from an in vitro permeability assay. Using permeability data from porcine kidney epithelial cells resulted in a 1.6-fold underestimation of Cmax at 800 mg. This difference likely reflects a concentration-dependent effect since the in vitro assay was performed at low micromolar concentrations, whereas the local gastrointestinal concentration could be as high as molar range (800 mg in 250 ml). Total clearance consisted of hepatic, intestinal, and renal components. Renal clearance was set as passive renal filtration (i.e., product of free fraction and glomerular filtration rate), which was in line with the preliminary estimate of renal clearance in the phase I study. Metabolic clearances in the gut and liver were based on system parameters for the expression levels of CYP3A4 in each gut compartment and the average expression of CYP3A4 in liver. Human liver microsome kinetics (Km and Vmax) were used as initial inputs to model tazemetostat hepatic clearance in humans, prior to refinement based on observed clinical data, and it was assumed that the fraction of the drug metabolized by CYP3A4 was close to unity. The refinement of Km and Vmax values was not substantial, being within 20% of the experimentally determined values in human liver microsomes. Although the tazemetostat accumulation ratio (AUCday15/AUCday1) was suggestive of a modest induction of metabolism, omitting to include CYP3A4 induction was not detrimental to the simulation of the steady-state (day 15) exposures (Fig. 3A). Because accumulation ratios remained modest throughout the dosing range (accumulation ratio of 0.8 to 0.5), and since the simulated concentration-time profile was very similar to the observed, CYP3A4 induction kinetics were not added to this model. This is a limitation of the model, which likely accounts for the nonoptimal capture of the time dependence of tazemetostat PK, particularly at higher doses (although predicted day 1 AUC and Cmax were within 2-fold of observed values at doses of ≤800 mg). Furthermore, the scaling of CYP3A autoinduction to a pediatric population (i.e., ontogeny of nuclear hormone receptors such as the pregnane X receptor) is not sufficiently understood to allow quantitative modeling. This would have made it difficult to account for autoinduction in the simulations of the pediatric population even if it had been included in the adult model. Similarly, because no transporter kinetic data were available at the time of this study and sufficient P-gp ontogeny data are lacking, the effect of potential efflux transport on tazemetostat PK was not considered in this model. Given these limitations, the potential effects of autoinduction should be considered in the interpretation of the model results. Based on visual inspection, the model fit of the adult exposures (n = 24) adequately described the time-concentration profiles of tazemetostat and resulted in prediction of mean AUCss and oral clearance (CL/F) within ±30% of the observed results across the dose range. In addition, mean Cmax,ss values were predicted within 2-fold of the observed values (Fig. 3B; Rioux et al., 2015). The resultant adult model was then used to simulate tazemetostat steady-state exposures in discrete pediatric age ranges (6 months to 1 year, >1–2 years, >2–6 years, >6–12 years, and >12–18 years) after twice-daily administration of an oral suspension. In addition to pediatric physiology, the simulations accounted for ontogeny in hematocrit, plasma protein levels, and P450 expression (liver and gut). Although CYP3A4 ontogeny was modeled, the propensity for CYP3A7-mediated metabolism of tazemetostat was unknown and so no interaction with this major fetal isoform was accounted for in the pediatric model. Nevertheless, the effect of potential metabolism by CYP3A7 is expected to be limited, since projection of pediatric doses was limited to patients aged ≥6 months (Hines, 2009). Using this exposure-based analysis, pediatric doses that afforded the target AUC (80% of adult AUCss at 800 mg or 390 mg/m2 twice daily) were identified. On a BSA-normalized basis, the projected doses were comparable across the age range (1–18 years), from 270 to 305 mg/m2 twice daily, with a slightly lower projected dose of 230 mg/m2 twice daily, for the group aged 6 months to 1 year. This lower dose was attributable to multiple factors, including CYP3A4 ontogeny (60%–80% of adult expression in the liver for the group aged 6 months to 1 year), change in free fraction, volume of distribution, and blood flow. By contrast, the tazemetostat fraction absorbed was not projected to change significantly from adults to children for doses of ≤300 mg/m2 twice daily. Because the projected doses by age were comparable, population simulations were performed to determine the corresponding exposures for each age range at two fixed doses (240 and 300 mg/m2 twice daily; Fig. 3B), since using an age-independent dose was expected to facilitate trial design in this rare patient population. In brief, at 300 mg/m2 twice daily, mean steady-state AUC0–τ values ranged from approximately 80% to 100% of that observed in adults at 800 mg twice daily. By contrast, the mean steady-state Cmax was projected to range from 110% to 190% of that observed at 800 mg twice daily in adults, but within the safe and efficacious exposure range defined in adults at doses up to 1600 mg twice daily. A lower dose of 240 mg/m2 twice daily was projected to maintain AUCss values ≤ 80% of the adult target exposure across the age range of 1–18 years.

Fig. 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 3.

(A) Representative plots showing the model fit for tazemetostat interim adult PK data at steady state (day 15, n = 24) showed a good fit for each dose group. Symbols represent observed data from individuals (± S.D., n = 3 or 6 per dose) and the solid line represents the mean profiles predicted from the PBPK model in GastroPlus. Dotted gray lines represent the 90% confidence interval. (B) Predicted mean total steady-state plasma concentration-time profiles of tazemetostat administered as a 300-mg/m2 twice-daily or a 240-mg/m2 twice-daily oral dose across the age ranges, with mean measured total steady-state plasma concentration-time profile of tazemetostat administered as a 800-mg (390 mg/m2) or 1600-mg (780 mg/m2) twice-daily oral dose in adults (n = 6 per dose). BID, twice daily.

These studies demonstrate the prospective application of PBPK early in clinical development to support clinical trial design, and they have provided dosing recommendations for the ongoing trials of pinometostat and tazemetostat in pediatric patients.

Current Challenges and Limitations in Pediatric PBPK Modeling

Although multiple ontogenic physiologic processes are fairly well characterized in PBPK models, knowledge gaps remain because some growth and maturation trajectories are still only partially or poorly characterized. Among these, the Pediatric Biopharmaceutics Classification System Working Group determined that additional research was required to fully understand age-based changes in gastrointestinal fluid composition, motility, absorptive surface area, and pH ranges encountered along the gastrointestinal tract, all of which are critical factors needed to better understand age-dependent drug absorption (Abdel-Rahman et al., 2012).

For clearance, incomplete knowledge of the developmental patterns for some P450 isoforms, UGTs, and other conjugative enzymes remains a challenge. Partial coverage of the pediatric age spectrum, sparse data in extrahepatic tissues, and an undetermined fraction metabolized by CYP3A7 (a major isoform in neonates) are some of the factors that can make scaling adult clearance to children a complex endeavor (Kearns et al., 2003; Abdel-Rahman et al., 2012). In addition, the parameterization of P450 ontogeny functions used in PBPK models may be based on data derived either in vitro (expression and/or activity) or in vivo (activity only), and recent studies have highlighted the potential for marked differences between the two methods (Salem et al., 2014; Upreti and Wahlstrom, 2016). In parallel, understanding of the pharmacokinetic effect of the transporter-specific expression level in discrete age strata is still limited and remains a challenge when performing PBPK modeling for pediatric applications. Screening of the mRNA expression of uptake and efflux transporters in the human liver reveals that only a small number of isoforms are detectable up to 1 month after birth; in general, although expression increased with age, the timing at which adult mRNA levels were reached did vary (Klaassen and Aleksunes, 2010; Abdel-Rahman et al., 2012; Mooij et al., 2014). Although more data are becoming available on transporter expression levels, studies on transporter activity in pediatric populations are still scarce. As data highlighting the importance of drug transporters in adult medicine continue to emerge, this critical knowledge gap in the pediatric population becomes even more evident, as raised by the National Institutes of Health Pediatric Transporter Working Group in their recent white paper (Brouwer et al., 2015).

Although PBPK modeling is increasingly used for prediction of DDI potential in adult populations, there are a number of gaps in the robust prediction of DDI potential in pediatrics. From the perspective of a victim drug interaction, the magnitude of a change in metabolic clearance depends on the fraction metabolized by the inhibited or induced pathway(s), and age-related changes owing to ontogeny may affect the predictive accuracy in the quantitative contribution of individual clearance mechanisms (Salem et al., 2013). Furthermore, nuclear hormone receptor mRNA expression (e.g., the pregnane X receptor and the constitutive androstane receptor, known to regulate expression of genes involved in drug disposition) demonstrate considerable interindividual variability in human fetal and pediatric livers (Vyhlidal et al., 2006), with the effect on P450 induction not yet fully understood. This confounds modeling of complex drug interactions such as simultaneous inhibition and induction of metabolic pathways and/or transporters by a drug (and/or its metabolites), presenting significant challenges to the prediction of clinical DDIs in adults (Varma et al., 2015) and even more so in pediatric populations.

Ultimately, the ability to model the effect of disease state on PK would further enhance the PBPK modeling approach by factoring in the physiologic sequelae of cancer in children relative to the current databases for healthy pediatric subjects; this is not currently a surmountable task, given the pleiotropic nature of malignant neoplastic disease. Furthermore, the extension of pediatric PBPK modeling to include a pharmacodynamics component remains an area of much-needed research.

The evaluation of the predictive performance of published PBPK models is an area that is in need of further scrutiny in moving toward a more standardized and systematic reporting of PBPK modeling and simulation efforts, and this issue has been raised by a number of groups from academia, industry, and regulatory authorities (Maharaj and Edginton, 2014; Zhao, 2014; Sager et al., 2015).

Summary and Conclusions

PBPK modeling and simulation in pediatric oncology is still quite nascent, limited to a few well characterized drugs and often performed as part of a retrospective analysis validating the methodology in this area of drug research. Nevertheless, these studies demonstrated the ability of PBPK modeling to simulate PK in pediatric patients, refine clinical trial design, inform on covariates influencing drug exposure levels in children, and evaluate the potential for clinically relevant DDIs. The case studies presented for pinometostat and tazemetostat provide further validation of their prospective utility a priori to the first-in-pediatrics clinical investigation, providing guidance on starting dose regimen and judicious selection of blood sampling times. Modeling and simulation endeavors in this area of drug research are particularly challenging since the patient population is rare. First-in-pediatric study populations tend to be relatively small, which impedes the statistical rigor usually associated with evaluating age as a covariate in later-stage trials. However, as detailed herein, PBPK modeling applied with a fit-for-purpose tenet can enable informed decisions on the starting dose in phase I pediatric oncology trials, in an effort to ensure that the age continuum is treated with a starting dose at or close to a safe and efficacious one. The future direction of PBPK modeling in the field of pediatric oncology drug development is likely to see expanded application in early clinical development with the assessment of DDI potential, influence of disease covariates and providing mechanistic insight into observed variability in the target patient population.

Acknowledgments

The authors thank Epizyme colleagues, clinical investigators, and their teams and, most importantly, the patients and families who participated in the studies reported herein.

Authorship Contributions

Wrote or contributed to the writing of the manuscript: Rioux, Waters

Footnotes

    • Received October 22, 2015.
    • Accepted March 1, 2016.
  • ↵1 Current affiliation: H3 Biomedicine, Cambridge, Massachusetts.

  • ↵2 Current affiliation: Syros Pharmaceuticals, Cambridge, Massachusetts.

  • dx.doi.org/10.1124/dmd.115.068031.

Abbreviations

AAG
α-1-acid glycoprotein
ADME
absorption, distribution, metabolism, and excretion
AUC
area under the curve
AUCss
area under the curve at steady state
BSA
body surface area
CL/F
oral clearance
Cmax
peak serum concentration
CLint
intrinsic clearance
DDI
drug–drug interaction
DME
drug-metabolizing enzyme
DOT1L
disruptor of telomeric silencing 1-like
EPZ-5676
9-[5-deoxy-5-[[cis-3-[2-[6-(1,1-dimethylethyl)-1H-benzimidazol-2-yl]ethyl]cyclobutyl](1-methylethyl)amino]-β-d-ribofuranosyl]-9H-purin-6-amine
EPZ-6438
N-[(1,2-dihydro-4,6-dimethyl-2-oxo-3-pyridinyl)methyl]-5-[ethyl(tetrahydro-2H-pyran-4-yl)amino]-4-methyl-4′-(4-morpholinylmethyl)-[1,1′-biphenyl]-3-carboxamide
FDA
Food and Drug Administration
GST
glutathione-S-transferase
HMT
histone methyltransferase
INI
integrase interactor
MLL
mixed lineage leukemia
MLL-r
mixed lineage leukemia–rearranged
P450
cytochrome P450
PBPK
physiologically based pharmacokinetics
P-gp
P-glycoprotein
PK
pharmacokinetics
t1/2
terminal half-life
UGT
uridine diphosphoglucuronosyl transferase
  • Copyright © 2016 by The American Society for Pharmacology and Experimental Therapeutics

References

  1. ↵
    1. Abdel-Rahman SM,
    2. Amidon GL,
    3. Kaul A,
    4. Lukacova V,
    5. Vinks AA,
    6. Knipp GT, and
    7. Members of the BCS Task Force
    (2012) Summary of the National Institute of Child Health and Human Development–Best Pharmaceuticals for Children Act Pediatric Formulation Initiatives Workshop–Pediatric Biopharmaceutics Classification System Working Group. Clin Ther 34:S11–S24.
    OpenUrlCrossRefPubMed
  2. ↵
    1. Adamson PC,
    2. Houghton PJ,
    3. Perilongo G, and
    4. Pritchard-Jones K
    (2014) Drug discovery in paediatric oncology: roadblocks to progress. Nat Rev Clin Oncol 11:732–739.
    OpenUrlCrossRef
  3. ↵
    1. Altamimi MI,
    2. Sammons H, and
    3. Choonara I
    (2015) Inter-individual variation in midazolam clearance in children. Arch Dis Child 100:95–100.
    OpenUrlAbstract/FREE Full Text
  4. ↵
    1. Barrett JS,
    2. Della Casa Alberighi O,
    3. Läer S, and
    4. Meibohm B
    (2012) Physiologically based pharmacokinetic (PBPK) modeling in children. Clin Pharmacol Ther 92:40–49.
    OpenUrlCrossRefPubMed
  5. ↵
    1. Barrett JS,
    2. Zhang AY,
    3. Urtishak KA,
    4. Danet-Desnoyers G,
    5. Carroll AJ,
    6. Hunger SP, and
    7. Felix CA
    (2011) Predicting clinical dose-exposure and exposure-response relationships of pan-antiapoptotic BCL-2 family inhibitor obatoclax in MLL rearranged infant leukemias from preclinical disease models and adult experience (Abstract). Blood 118:2580.
    OpenUrl
  6. ↵
    1. Basavapathruni A,
    2. Olhava EJ,
    3. Daigle SR,
    4. Therkelsen CA,
    5. Jin L,
    6. Boriack-Sjodin PA,
    7. Allain CJ,
    8. Klaus CR,
    9. Raimondi A,
    10. Scott MP,
    11. et al.
    (2014) Nonclinical pharmacokinetics and metabolism of EPZ-5676, a novel DOT1L histone methyltransferase inhibitor. Biopharm Drug Dispos 35:237–252.
    OpenUrlCrossRefPubMed
  7. ↵
    1. Batchelor HK,
    2. Fotaki N, and
    3. Klein S
    (2014) Paediatric oral biopharmaceutics: key considerations and current challenges. Adv Drug Deliv Rev 73:102–126.
    OpenUrlCrossRefPubMed
  8. ↵
    1. Björkman S
    (2005) Prediction of drug disposition in infants and children by means of physiologically based pharmacokinetic (PBPK) modelling: theophylline and midazolam as model drugs. Br J Clin Pharmacol 59:691–704.
    OpenUrlCrossRefPubMed
  9. ↵
    1. Björkman S
    (2006) Prediction of cytochrome p450-mediated hepatic drug clearance in neonates, infants and children: how accurate are available scaling methods? Clin Pharmacokinet 45:1–11.
    OpenUrlCrossRefPubMed
  10. ↵
    1. Blake MJ,
    2. Gaedigk A,
    3. Pearce RE,
    4. Bomgaars LR,
    5. Christensen ML,
    6. Stowe C,
    7. James LP,
    8. Wilson JT,
    9. Kearns GL, and
    10. Leeder JS
    (2007) Ontogeny of dextromethorphan O- and N-demethylation in the first year of life. Clin Pharmacol Ther 81:510–516.
    OpenUrlCrossRefPubMed
  11. ↵
    1. Block M
    (2015) Physiologically based pharmacokinetic and pharmacodynamic modeling in cancer drug development: status, potential and gaps. Expert Opin Drug Metab Toxicol 11:743–756.
    OpenUrlCrossRef
  12. ↵
    1. Breitkreutz J and
    2. Boos J
    (2007) Paediatric and geriatric drug delivery. Expert Opin Drug Deliv 4:37–45.
    OpenUrlCrossRefPubMed
  13. ↵
    1. Brouwer KLR,
    2. Aleksunes LM,
    3. Brandys B,
    4. Giacoia GP,
    5. Knipp G,
    6. Lukacova V,
    7. Meibohm B,
    8. Nigam SK,
    9. Rieder M,
    10. de Wildt SN, and
    11. Pediatric Transporter Working Group
    (2015) Human ontogeny of drug transporters: review and recommendations of the pediatric transporter working group. Clin Pharmacol Ther 98:266–287.
    OpenUrlCrossRefPubMed
  14. ↵
    1. Daigle SR,
    2. Olhava EJ,
    3. Therkelsen CA,
    4. Basavapathruni A,
    5. Jin L,
    6. Boriack-Sjodin PA,
    7. Allain CJ,
    8. Klaus CR,
    9. Raimondi A,
    10. Scott MP,
    11. et al.
    (2013) Potent inhibition of DOT1L as treatment of MLL-fusion leukemia. Blood 122:1017–1025.
    OpenUrlAbstract/FREE Full Text
  15. ↵
    1. Diestelhorst C,
    2. Boos J,
    3. McCune JS,
    4. Russell J,
    5. Kangarloo SB, and
    6. Hempel G
    (2014) Predictive performance of a physiologically based pharmacokinetic model of busulfan in children. Pediatr Hematol Oncol 31:731–742.
    OpenUrlCrossRef
  16. ↵
    1. Ehrnebo M,
    2. Agurell S,
    3. Jalling B, and
    4. Boréus LO
    (1971) Age differences in drug binding by plasma proteins: studies on human foetuses, neonates and adults. Eur J Clin Pharmacol 3:189–193.
    OpenUrlCrossRefPubMed
  17. ↵
    1. European Medicines Agency
    (2013) EMA/CHMP/161314/2013: Assessment Report for Glivec (Imatinib) (May 30, 2013), European Medicines Agency, London.
  18. ↵
    1. European Medicines Agency
    (2014) EMA/CHMP/211243/2014: Concept paper on Qualification and Reporting of Physiologically-Based Pharmacokinetic (PBPK) Modelling and Analyses (June 26, 2014), European Medicines Agency, London.
  19. ↵
    1. Evans WE,
    2. Crom WR,
    3. Tsiatis A,
    4. Green AA,
    5. Hayes FA, and
    6. Pratt CB
    (1982) Pharmacokinetic modeling of cisplatin disposition in children and adolescents with cancer. Cancer Chemother Pharmacol 10:22–26.
    OpenUrlPubMed
  20. ↵
    1. Gibbs JP,
    2. Liacouras CA,
    3. Baldassano RN, and
    4. Slattery JT
    (1999) Up-regulation of glutathione S-transferase activity in enterocytes of young children. Drug Metab Dispos 27:1466–1469.
    OpenUrlAbstract/FREE Full Text
  21. ↵
    Hines RN (2009) The ontogeny of drug metabolizing enzymes/pediatric exclusivity, in Handbook of Drug Metabolism (Pearson PG and Wienkers LC eds) 2nd ed, pp 227–239, Informa Healthcare, New York.
  22. ↵
    1. Horton TM and
    2. Berg SL
    (2011) Educational paper. The development of new therapies for pediatric oncology. Eur J Pediatr 170:555–559.
    OpenUrlCrossRefPubMed
  23. ↵
    1. Huang SM,
    2. Abernethy DR,
    3. Wang Y,
    4. Zhao P, and
    5. Zineh I
    (2013) The utility of modeling and simulation in drug development and regulatory review. J Pharm Sci 102:2912–2923.
    OpenUrlCrossRefPubMed
  24. ↵
    1. Jones HM,
    2. Chen Y,
    3. Gibson C,
    4. Heimbach T,
    5. Parrott N,
    6. Peters SA,
    7. Snoeys J,
    8. Upreti VV,
    9. Zheng M, and
    10. Hall SD
    (2015) Physiologically based pharmacokinetic modeling in drug discovery and development: a pharmaceutical industry perspective. Clin Pharmacol Ther 97:247–262.
    OpenUrlCrossRefPubMed
  25. ↵
    1. Jones HM,
    2. Gardner IB, and
    3. Watson KJ
    (2009) Modelling and PBPK simulation in drug discovery. AAPS J 11:155–166.
    OpenUrlCrossRefPubMed
  26. ↵
    1. Kearns GL,
    2. Abdel-Rahman SM,
    3. Alander SW,
    4. Blowey DL,
    5. Leeder JS, and
    6. Kauffman RE
    (2003) Developmental pharmacology--drug disposition, action, and therapy in infants and children. N Engl J Med 349:1157–1167.
    OpenUrlCrossRefPubMed
  27. ↵
    1. Kersting G,
    2. Willmann S,
    3. Würthwein G,
    4. Lippert J,
    5. Boos J, and
    6. Hempel G
    (2012) Physiologically based pharmacokinetic modelling of high- and low-dose etoposide: from adults to children. Cancer Chemother Pharmacol 69:397–405.
    OpenUrlCrossRefPubMed
  28. ↵
    1. Khalil F and
    2. Läer S
    (2011) Physiologically based pharmacokinetic modeling: methodology, applications, and limitations with a focus on its role in pediatric drug development. J Biomed Biotechnol 2011:907461.
    OpenUrlPubMed
  29. ↵
    1. Klaassen CD and
    2. Aleksunes LM
    (2010) Xenobiotic, bile acid, and cholesterol transporters: function and regulation. Pharmacol Rev 62:1–96.
    OpenUrlAbstract/FREE Full Text
  30. ↵
    1. Knutson SK,
    2. Kawano S,
    3. Minoshima Y,
    4. Warholic NM,
    5. Huang KC,
    6. Xiao Y,
    7. Kadowaki T,
    8. Uesugi M,
    9. Kuznetsov G,
    10. Kumar N,
    11. et al.
    (2014) Selective inhibition of EZH2 by EPZ-6438 leads to potent antitumor activity in EZH2-mutant non-Hodgkin lymphoma. Mol Cancer Ther 13:842–854.
    OpenUrlAbstract/FREE Full Text
  31. ↵
    1. Knutson SK,
    2. Warholic NM,
    3. Wigle TJ,
    4. Klaus CR,
    5. Allain CJ,
    6. Raimondi A,
    7. Porter Scott M,
    8. Chesworth R,
    9. Moyer MP,
    10. Copeland RA,
    11. et al.
    (2013) Durable tumor regression in genetically altered malignant rhabdoid tumors by inhibition of methyltransferase EZH2. Proc Natl Acad Sci USA 110:7922–7927.
    OpenUrlAbstract/FREE Full Text
  32. ↵
    1. Kostewicz ES,
    2. Aarons L,
    3. Bergstrand M,
    4. Bolger MB,
    5. Galetin A,
    6. Hatley O,
    7. Jamei M,
    8. Lloyd R,
    9. Pepin X,
    10. Rostami-Hodjegan A,
    11. et al.
    (2014) PBPK models for the prediction of in vivo performance of oral dosage forms. Eur J Pharm Sci 57:300–321.
    OpenUrlCrossRef
  33. ↵
    1. Li J and
    2. Gwilt P
    (2002) The effect of malignant effusions on methotrexate disposition. Cancer Chemother Pharmacol 50:373–382.
    OpenUrlCrossRefPubMed
  34. ↵
    1. Maharaj AR and
    2. Edginton AN
    (2014) Physiologically based pharmacokinetic modeling and simulation in pediatric drug development. CPT Pharmacometrics Syst Pharmacol 3:e150.
    OpenUrlCrossRef
  35. ↵
    1. Mahmood I
    (2006) Prediction of drug clearance in children from adults: a comparison of several allometric methods. Br J Clin Pharmacol 61:545–557.
    OpenUrlCrossRefPubMed
  36. ↵
    1. Mooij MG,
    2. Schwarz UI,
    3. de Koning BA,
    4. Leeder JS,
    5. Gaedigk R,
    6. Samsom JN,
    7. Spaans E,
    8. van Goudoever JB,
    9. Tibboel D,
    10. Kim RB,
    11. et al.
    (2014) Ontogeny of human hepatic and intestinal transporter gene expression during childhood: age matters. Drug Metab Dispos 42:1268–1274.
    OpenUrlAbstract/FREE Full Text
  37. ↵
    1. Ogungbenro K,
    2. Aarons L, and
    3. CRESim & Epi-CRESim Project Groups
    (2014a) Physiologically based pharmacokinetic modelling of methotrexate and 6-mercaptopurine in adults and children. Part 1: methotrexate. J Pharmacokinet Pharmacodyn 41:159–171.
    OpenUrlCrossRef
  38. ↵
    1. Ogungbenro K,
    2. Aarons L, and
    3. CRESim & Epi-CRESim Project Groups
    (2014b) Physiologically based pharmacokinetic modelling of methotrexate and 6-mercaptopurine in adults and children. Part 2: 6-mercaptopurine and its interaction with methotrexate. J Pharmacokinet Pharmacodyn 41:173–185.
    OpenUrlCrossRef
  39. ↵
    1. Poley JR,
    2. Dower JC,
    3. Owen CA Jr., and
    4. Stickler GB
    (1964) Bile acids in infants and children. J Lab Clin Med 63:838–846.
    OpenUrlPubMed
  40. ↵
    1. Pons G,
    2. Carrier O,
    3. Richard MO,
    4. Rey E,
    5. d’Athis P,
    6. Moran C,
    7. Badoual J, and
    8. Olive G
    (1988) Developmental changes of caffeine elimination in infancy. Dev Pharmacol Ther 11:258–264.
    OpenUrlPubMed
  41. ↵
    Ribrag V, Soria JC, Thomson B, Reyderman L, Kraljevic S, Chen R, Salazar P, Kuznetsov G, Keilhack H, Blakemore SJ, et al. (2015) Phase 1 study of EPZ-6438 (E7438), an Enhancer of Zeste Homolog-2 (EZH2) inhibitor: dose determination and preliminary activity in non-Hodgkin lymphoma, in Proceedings of the International Conference on Malignant Lymphoma; 2015 June 17–20; ICML, Bellinzona, Italy.
  42. ↵
    1. Ries LAG,
    2. Smith MA,
    3. Gurney JG,
    4. Linet M,
    5. Tamra T,
    6. Young JL, and
    7. Bunin GR
    , editors (1999) Cancer Incidence and Survival among Children and Adolescents: United States SEER Program 1975–1995, SEER Program Pub No. 99-4649, National Cancer Institute, Bethesda, MD.
  43. ↵
    1. Rioux N,
    2. Keilhack H,
    3. Thomson B,
    4. Chesworth R,
    5. Copeland RA,
    6. Hedrick E,
    7. Ho P, and
    8. Waters NJ
    (2015) Physiologically-based pharmacokinetic modeling to support clinical investigation of the EZH2 inhibitor, tazemetostat (EPZ-6438) in INI1-deficient pediatric tumors (Abstract). Mol Cancer Ther 14:A135.
    OpenUrlCrossRef
  44. ↵
    1. Rossig C,
    2. Juergens H, and
    3. Berdel WE
    (2011) New targets and targeted drugs for the treatment of cancer: an outlook to pediatric oncology. Pediatr Hematol Oncol 28:539–555.
    OpenUrlCrossRefPubMed
  45. ↵
    1. Sager JE,
    2. Yu J,
    3. Ragueneau-Majlessi I, and
    4. Isoherranen N
    (2015) Physiologically based pharmacokinetic (PBPK) modeling and simulation approaches: a systematic review of published models, applications and model verification. Drug Metab Dispos 43:1823–1837.
    OpenUrlAbstract/FREE Full Text
  46. ↵
    1. Salem F,
    2. Johnson TN,
    3. Abduljalil K,
    4. Tucker GT, and
    5. Rostami-Hodjegan A
    (2014) A re-evaluation and validation of ontogeny functions for cytochrome P450 1A2 and 3A4 based on in vivo data. Clin Pharmacokinet 53:625–636.
    OpenUrlCrossRefPubMed
  47. ↵
    1. Salem F,
    2. Johnson TN,
    3. Barter ZE,
    4. Leeder JS, and
    5. Rostami-Hodjegan A
    (2013) Age related changes in fractional elimination pathways for drugs: assessing the impact of variable ontogeny on metabolic drug-drug interactions. J Clin Pharmacol 53:857–865.
    OpenUrlCrossRefPubMed
  48. ↵
    1. Shepard T,
    2. Scott G,
    3. Cole S,
    4. Nordmark A, and
    5. Bouzom F
    (2015) Physiologically based models in regulatory submissions: Output from the ABPI/MHRA forum on physiologically based modeling and simulation. CPT Pharmacometrics Syst Pharmacol 4:221–225.
    OpenUrlCrossRefPubMed
  49. ↵
    Shukla N, O’Brien M, Silverman L, Pauly M, Wetmore C, Loh M, Whitlock J, Brown P, Daigle S, Blakemore SJ, et al. (2015) Preliminary report of the phase 1 study of the DOT1L inhibitor pinometostat (EPZ-5676) in children with MLL-r acute leukemia: safety, exposure and evidence of target inhibition, in Proceedings of the 57th American Society of Hematology Annual Meeting; 2015 December 5–8; Orlando, FL. American Society of Hematology, Washington, DC.
  50. ↵
    1. Skolnik JM,
    2. Barrett JS,
    3. Jayaraman B,
    4. Patel D, and
    5. Adamson PC
    (2008) Shortening the timeline of pediatric phase I trials: the rolling six design. J Clin Oncol 26:190–195.
    OpenUrlAbstract/FREE Full Text
  51. ↵
    1. Stein E,
    2. Garcia-Manero G,
    3. Rizzieri D,
    4. Savona MR,
    5. Tibes R,
    6. Altman J,
    7. Jongen-Lavrencic M,
    8. Döhner H,
    9. Armstrong S,
    10. Pollock RM,
    11. et al.
    (2014) The DOT1L inhibitor EPZ-5676: safety and activity in relapsed/refractory patients with MLL-rearranged leukemia (Abstract). Blood 124:387.
    OpenUrl
  52. ↵
    1. Thai HT,
    2. Mazuir F,
    3. Cartot-Cotton S, and
    4. Veyrat-Follet C
    (2015) Optimizing pharmacokinetic bridging studies in paediatric oncology using physiologically-based pharmacokinetic modelling: application to docetaxel. Br J Clin Pharmacol 80:534–547.
    OpenUrlCrossRef
  53. ↵
    1. Upreti VV and
    2. Wahlstrom JL
    (2016) Meta-analysis of hepatic cytochrome P450 ontogeny to underwrite the prediction of pediatric pharmacokinetics using physiologically based pharmacokinetic modeling. J Clin Pharmacol 56:266–283.
    OpenUrlCrossRef
  54. ↵
    1. Varma MV,
    2. Pang KS,
    3. Isoherranen N, and
    4. Zhao P
    (2015) Dealing with the complex drug-drug interactions: towards mechanistic models. Biopharm Drug Dispos 36:71–92.
    OpenUrlCrossRefPubMed
  55. ↵
    1. Vyhlidal CA,
    2. Gaedigk R, and
    3. Leeder JS
    (2006) Nuclear receptor expression in fetal and pediatric liver: correlation with CYP3A expression. Drug Metab Dispos 34:131–137.
    OpenUrlAbstract/FREE Full Text
  56. ↵
    1. Wagner C,
    2. Zhao P,
    3. Pan Y,
    4. Hsu V,
    5. Grillo J,
    6. Huang SM, and
    7. Sinha V
    (2015) Application of physiologically based pharmacokinetic (PBPK) modeling to support dose selection: report of an FDA public workshop on PBPK. CPT Pharmacometrics Syst Pharmacol 4:226–230.
    OpenUrlCrossRef
  57. ↵
    1. Waters NJ,
    2. Thomson B,
    3. Gardner I,
    4. Johnson TN,
    5. Olhava EJ,
    6. Pollock RM,
    7. Legler M,
    8. Copeland RA, and
    9. Hedrick E
    (2014) Pediatric dose determinations for the phase I study of the DOT1L inhibitor, EPZ-5676, in MLL-r acute leukemia: leveraging early clinical data in adults through physiologically-based pharmacokinetic modeling. Blood 124:3619.
    OpenUrl
  58. ↵
    1. West GB,
    2. Brown JH, and
    3. Enquist BJ
    (1999) The fourth dimension of life: fractal geometry and allometric scaling of organisms. Science 284:1677–1679.
    OpenUrlAbstract/FREE Full Text
  59. ↵
    Zhao P (2014) Application of PBPK modeling and simulations in pediatric drug development, in Proceedings of the American Society for Clinical Pharmacology and Therapeutics Annual Meeting; 2014 March 20; Atlanta, GA. American Society for Clinical Pharmacology and Therapeutics, Alexandria, VA.
  60. ↵
    1. Zhao P,
    2. Rowland M, and
    3. Huang SM
    (2012) Best practice in the use of physiologically based pharmacokinetic modeling and simulation to address clinical pharmacology regulatory questions. Clin Pharmacol Ther 92:17–20.
    OpenUrlCrossRefPubMed
PreviousNext
Back to top

In this issue

Drug Metabolism and Disposition: 44 (7)
Drug Metabolism and Disposition
Vol. 44, Issue 7
1 Jul 2016
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by author
  • Editorial Board (PDF)
  • Front Matter (PDF)
Download PDF
Article Alerts
Sign In to Email Alerts with your Email Address
Email Article

Thank you for sharing this Drug Metabolism & Disposition article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Physiologically Based Pharmacokinetic Modeling in Pediatric Oncology Drug Development
(Your Name) has forwarded a page to you from Drug Metabolism & Disposition
(Your Name) thought you would be interested in this article in Drug Metabolism & Disposition.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Research ArticleSpecial Section on Pediatric Drug Disposition and Pharmacokinetics—Minireview

PBPK Modeling in Pediatric Oncology Drug Development

Nathalie Rioux and Nigel J. Waters
Drug Metabolism and Disposition July 1, 2016, 44 (7) 934-943; DOI: https://doi.org/10.1124/dmd.115.068031

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Share
Research ArticleSpecial Section on Pediatric Drug Disposition and Pharmacokinetics—Minireview

PBPK Modeling in Pediatric Oncology Drug Development

Nathalie Rioux and Nigel J. Waters
Drug Metabolism and Disposition July 1, 2016, 44 (7) 934-943; DOI: https://doi.org/10.1124/dmd.115.068031
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Introduction
    • Case Studies of PBPK Modeling in Pediatric Oncology
    • Pediatric PBPK Modeling In-House at Epizyme
    • Current Challenges and Limitations in Pediatric PBPK Modeling
    • Summary and Conclusions
    • Acknowledgments
    • Authorship Contributions
    • Footnotes
    • Abbreviations
    • References
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF

Related Articles

Cited By...

More in this TOC Section

  • Age-Dependent Hepatic Drug Transport in Pediatrics
Show more Special Section on Pediatric Drug Disposition and Pharmacokinetics—Minireview

Similar Articles

  • Home
  • Alerts
Facebook   Twitter   LinkedIn   RSS

Navigate

  • Current Issue
  • Fast Forward by date
  • Fast Forward by section
  • Latest Articles
  • Archive
  • Search for Articles
  • Feedback
  • ASPET

More Information

  • About DMD
  • Editorial Board
  • Instructions to Authors
  • Submit a Manuscript
  • Customized Alerts
  • RSS Feeds
  • Subscriptions
  • Permissions
  • Terms & Conditions of Use

ASPET's Other Journals

  • Journal of Pharmacology and Experimental Therapeutics
  • Molecular Pharmacology
  • Pharmacological Reviews
  • Pharmacology Research & Perspectives
ISSN 1521-009X (Online)

Copyright © 2021 by the American Society for Pharmacology and Experimental Therapeutics