Exploring the impact of CYP2D6 and UGT2B7 gene-drug interactions, and CYP-mediated DDI on oxycodone and oxymorphone pharmacokinetics using physiologically-based pharmacokinetic modeling and simulation.

Oxycodone is one of the most commonly used opioids to treat moderate to severe pain. It is metabolized mainly by CYP3A4 and CYP2D6, while only a small fraction of the dose is excreted unchanged into the urine. Oxymorphone, the metabolite primarily formed by CYP2D6, has a 40- to 60-fold higher mu-opioid receptor affinity than the parent compound. While CYP2D6-mediated gene-drug-interactions (GDIs) and drug-drug interactions (DDIs) are well-studied, they only account for a portion of the variability in oxycodone and oxymorphone exposure. The combined impact of CYP2D6-mediated GDIs and DDIs, CYP3A4-mediated DDIs, and UGT2B7 GDIs is not fully understood yet and hard to study in head-to-head clinical trials given the relatively large number of scenarios. Instead, we propose the use of a physiologically-based pharmacokinetic model that integrates available information on oxycodone's metabolism to characterize and predict the impact of DDIs and GDIs on the exposure of oxycodone and its major, pharmacologically-active metabolite oxymorphone. To this end, we first developed and verified a PBPK model for oxycodone and its metabolites using published clinical data. The verified model was then applied to determine the dose-exposure relationship of oxycodone and oxymorphone stratified by CYP2D6 and UGT2B7 phenotypes respectively, and administered perpetrators of CYP-based drug interactions. Our simulations demonstrate that the combination of CYP2D6 UM and a UGT2B7Y (268) mutation may lead to a 2.3-fold increase in oxymorphone exposure compared to individuals who are phenotyped as CYP2D6 NM / UGT2B7 NM. The extent of oxymorphone exposure increases up to 3.2-fold in individuals concurrently taking CYP3A4 inhibitors, such as ketoconazole. Inhibition of the CYP3A4 pathway results in a relative increase in the partial metabolic clearance of oxycodone to oxymorphone. Oxymorphone is impacted to a higher extent by GDIs and DDIs than oxycodone. We predict oxymorphone exposure to be highest in CYP2D6 UMs/UGT2B7 PMs in the presence of ketoconazole (strong CYP3A4 index inhibitor) and lowest in CYP2D6 PMs/UGT2B7 NMs in the presence of rifampicin (strong CYP3A4 index inducer) covering a 55-fold exposure range.

Application of a patient-centered reverse translational systems-based approach to understand mechanisms of an adverse drug reaction of immune checkpoint inhibitors.

Immunotherapy became a key pillar of cancer therapeutics with the approvals of ipilimumab, nivolumab, and pembrolizumab, which inhibit either cytotoxic T-lymphocyte antigen-4 (CTLA-4) or programmed death-1 (PD-1) that are negative regulators of T-cell activation. However, boosting T-cell activation is often accompanied by autoimmunity, leading to adverse drug reactions (ADRs), including high grade 3-4 colitis and its severe complications whose prevalence may reach 14% for combination checkpoint inhibitors. In this research, we investigated how mechanistic differences between anti-CTLA-4 (ipilimumab) and anti-PD-1 (nivolumab and pembrolizumab) affect colitis, a general class toxicity. The data analytical platform Molecular Health Effect was utilized to map population ADR data from the US Food and Drug Administration (FDA) Adverse Event Reporting System to chemical and biological databases for hypothesis generation regarding the underlying molecular mechanisms causing colitis. Disproportionality analysis was used to assess the statistical relevance between adverse events of interest and molecular causation. We verified that the anti-CTLA-4 drug is associated with an approximately three-fold higher proportional reporting ratio associated with colitis than those of the anti-PD-1 drugs. The signal of the molecular mechanisms, including signaling pathways of inflammatory cytokines, was statistically insignificant to test the hypothesis that the severer rate of colitis associated with ipilimumab would be due to a greater magnitude of T-cell activation as a result of earlier response of the anti-CTLA-4 drug in the immune response. This patient-centered systems-based approach provides an exploratory process to better understand drug pair adverse events at pathway and target levels through reverse translation from postmarket surveillance safety reports.

Advancing drug safety science by integrating molecular knowledge with post-marketing adverse event reports.

Promising drug development efforts may frequently fail due to unintended adverse reactions. Several methods have been developed to analyze such data, aiming to improve pharmacovigilance and drug safety. In this work, we provide a brief review of key directions to quantitatively analyzing adverse events and explore the potential of augmenting these methods using additional molecular data descriptors. We argue that molecular expansion of adverse event data may provide a path to improving the insights gained through more traditional pharmacovigilance approaches. Examples include the ability to assess statistical relevance with respect to underlying biomolecular mechanisms, the ability to generate plausible causative hypotheses and/or confirmation where possible, the ability to computationally study potential clinical trial designs and/or results, as well as the further provision of advanced features incorporated in innovative methods, such as machine learning. In summary, molecular data expansion provides an elegant way to extend mechanistic modeling, systems pharmacology, and patient-centered approaches for the assessment of drug safety. We anticipate that such advances in real-world data informatics and outcome analytics will help to better inform public health, via the improved ability to prospectively understand and predict various types of drug-induced molecular perturbations and adverse events.

A case study of a patient-centered reverse translational systems-based approach to understand adverse event profiles in drug development.

Adverse drug reactions (ADRs) of targeted therapy drugs (TTDs) are frequently unexpected and long-term toxicities detract from exceptional efficacy of new TTDs. In this proof-of-concept study, we explored how molecular causation involved in trastuzumab-induced cardiotoxicity changes when trastuzumab was given in combination with doxorubicin, tamoxifen, paroxetine, or lapatinib. The data analytical platform Molecular Health Effect was utilized to map population ADR data from the US Food and Drug Administration (FDA) Adverse Event Reporting System to chemical and biological databases (such as UniProt and Reactome), for hypothesis generation regarding the underlying molecular mechanisms causing cardiotoxicity. Disproportionality analysis was used to assess the statistical relevance between adverse events of interest and molecular causation. Literature search was performed to compare the established hypotheses to published experimental findings. We found that the combination therapy of trastuzumab and doxorubicin may affect mitochondrial dysfunction in cardiomyocytes through different molecular pathways such as BCL-X and PGC-1α proteins, leading to a synergistic effect of cardiotoxicity. We found, on the other hand, that trastuzumab-induced cardiotoxicity would be diminished by concomitant use of tamoxifen, paroxetine, and/or lapatinib. Tamoxifen and paroxetine may cause less cardiotoxicity through an increase in antioxidant activities, such as glutathione conjugation. Lapatinib may decrease the apoptotic effects in cardiomyocytes by altering the effects of trastuzumab on BCL-X proteins. This patient-centered systems-based approach provides, based on the trastuzumab-induced ADR cardiotoxicity, an example of how to apply reverse translation to investigate ADRs at the molecular pathway and target level to understand the causality and prevalence during drug development of novel therapeutics.

Dosing Regimen Prediction and Confirmation With Rivaroxaban for Thromboprophylaxis in Children After the Fontan Procedure: Insights From the Phase III UNIVERSE Study.

Thrombosis remains an important complication for children with single-ventricle physiology following the Fontan procedure, and effective thromboprophylaxis is an important unmet medical need. To obviate conventional dose-finding studies and expedite clinical development, a rivaroxaban dose regimen for this indication was determined using a model-informed drug development approach. A physiologically based pharmacokinetic rivaroxaban model was used to predict a pediatric dosing regimen that would produce drug exposures similar to that of 10 mg once daily in adults. This regimen was used in an open-label, multicenter phase III study, which investigated the use of rivaroxaban for thromboprophylaxis in post-Fontan patients 2 to 8 years of age. The pharmacokinetics (PK) of rivaroxaban was assessed in part A (n = 12) and in part B (n = 64) of the UNIVERSE study. The safety and efficacy in the rivaroxaban group were compared to those in the acetylsalicylic acid group for 12 months. Pharmacodynamic end points were assessed in both parts of the study. Rivaroxaban exposures achieved in parts A and B were similar to the adult reference exposures. Prothrombin time also showed similarity to the adult reference. Exposure-response analysis did not identify a quantitative relationship between rivaroxaban exposures and efficacy/safety outcomes within the observed exposure ranges. A body weight-based dose regimen selected by physiologically based pharmacokinetic modeling was shown in the UNIVERSE study to be appropriate for thromboprophylaxis in the post-Fontan pediatric population. Model-based dose selection can support pediatric drug development and bridge adult dose data to pediatrics, thereby obviating the need for dose-finding studies in pediatric programs.

Brexpiprazole Pharmacokinetics in CYP2D6 Poor Metabolizers: Using Physiologically Based Pharmacokinetic Modeling to Optimize Time to Effective Concentrations.

Brexpiprazole is an oral antipsychotic agent indicated for use in patients with schizophrenia, or as adjunctive treatment for major depressive disorder. As cytochrome P450 (CYP) 2D6 contributes significantly to brexpiprazole metabolism, there is a label-recommended 50% reduction in dose among patients with the CYP2D6 poor metabolizer phenotype. This study uses a whole-body physiologically based pharmacokinetic (PBPK) model to compare the pharmacokinetics of brexpiprazole in patients known to be extensive metabolizers (EMs) and poor metabolizers (PMs). A PBPK model was constructed, verified, and validated against brexpiprazole clinical data, and simulations of 500 subjects were performed to establish the median time to effective concentrations in EMs and PMs. The PBPK simulations captured brexpiprazole PK well and demonstrated significant differences in the time to effective concentrations between EMs and PMs according to the label-recommended titration. Additionally, these simulations suggest that CYP2D6 PMs consistently achieve lower minimum concentrations during the dosing interval than CYP2D6 EMs. Simulations using an alternative dosing strategy of twice-daily dosing (as opposed to once daily) in PMs during the first week of brexpiprazole dosing yielded more consistent plasma concentrations between EMs and PMs, without exceeding the area under the plasma concentration-time curve observed in the EMs. Taken together, the results of these PBPK simulations suggest that product labeling for brexpiprazole titration in CYP2D6 PMs likely overcompensates for the decreased clearance seen in this population. We propose an alternative dosing strategy that decreases the time to effective concentrations and recommend a reevaluation of steady-state PK in this population to potentially allow for higher daily doses in CYP2D6 PMs.

Pharmacodynamic Modeling to Evaluate the Impact of Cimetidine, an OCT2 Inhibitor, on the Anticancer Effects of Cisplatin.

Despite potent anticancer activity, the clinical utilization of cisplatin is limited due to nephrotoxicity. As Organic Cation Transporter 2 (OCT2) has been shown to be one of the key transporters involved in the uptake of cisplatin into renal proximal tubules, OCT2 inhibitors such as cimetidine have been explored to suppress cisplatin-induced nephrotoxicity. Nonetheless, the impact of OCT2 inhibition or cimetidine on the anti-cancer effects of cisplatin has not been extensively examined. The main objective of the present study was to quantitatively characterize the anticancer effects of cisplatin and cimetidine and determine their nature of interactions in two cancer cell lines, OCT2-negative hepatocellular carcinoma (HCC) cell line, Huh7, and OCT2-positive breast cancer cell line, MDA-MB-468. First, we determined the static concentration-response curves of cisplatin and cimetidine as single agents. Next, with the help of three-dimensional (3D) response surface analyses and a competitive interaction model, we determined their nature of interactions at static concentrations to be modestly synergistic or additive in Huh7 and antagonistic in MDA-MB-468. These results were consistent with the cell-level pharmacodynamic (PD) modeling analysis which leveraged the time-course effects of drugs as single agents and drug combinations. Our developed PD model can be further used to design future preclinical studies to further investigate the cisplatin and cimetidine combinations in different in vitro and in vivo cancer models.

Perspective on model-informed drug development.

Model-informed drug development (MIDD) is a process intended to expedite drug development, enhance regulatory science, and produce benefits for patients. Quantitative modeling and simulation-principally by population pharmacokinetics (PK), exposure-response, and physiologically based pharmacokinetic (PBPK) analysis-is the technology that provides the capability to deploy MIDD across a range of applications. MIDD was codified in the 2017 Prescription Drug User Fee Act Reauthorization 1 (PDUFA VI, 2018-2022) and a performance goal was a MIDD pilot program to hold 2 to 4 industry-U.S. Food and Drug Administration (FDA) paired meetings quarterly through 2022.

Physiologically-based pharmacokinetics modeling to investigate formulation factors influencing the generic substitution of dabigatran etexilate.

The exposure-response relationship of direct acting oral anti-coagulants (DOACs) for bleeding risk is steep relative to ischemic stroke reduction. As a result, small changes in exposure may lead to bleeding events. The overall goal of this project was to determine the effect of critical formulation parameters on the pharmacokinetics (PKs) and thus safety and efficacy of generic DOACs. In this first installment of our overall finding, we developed and verified a physiologically-based PK (PBPK) model for dabigatran etexilate (DABE) and its metabolites. The model was developed following a middle out approach leveraging available in vitro and in vivo data. External validity of the model was confirmed by overlapping predicted and observed PK profiles for DABE as well as free and total dabigatran for a dataset not used during model development. The verified model was applied to interrogate the impact of modulating the microenvironment pH on DABE systemic exposure. The PBPK exploratory analyses highlighted the high sensitivity of DABE exposure to supersaturation ratio and precipitation kinetics.

Quantitative Systems Pharmacology Model-Based Predictions of Clinical Endpoints to Optimize Warfarin and Rivaroxaban Anti-Thrombosis Therapy.

BACKGROUND: Tight monitoring of efficacy and safety of anticoagulants such as warfarin is imperative to optimize the benefit-risk ratio of anticoagulants in patients. The standard tests used are measurements of prothrombin time (PT), usually expressed as international normalized ratio (INR), and activated partial thromboplastin time (aPTT). OBJECTIVE: To leverage a previously developed quantitative systems pharmacology (QSP) model of the human coagulation network to predict INR and aPTT for warfarin and rivaroxaban, respectively. METHODS: A modeling and simulation approach was used to predict INR and aPTT measurements of patients receiving steady-state anticoagulation therapy. A previously developed QSP model was leveraged for the present analysis. The effect of genetic polymorphisms known to influence dose response of warfarin (CYP2C9, VKORC1) were implemented into the model by modifying warfarin clearance (CYP2C9 *1: 0.2 L/h; *2: 0.14 L/h, *3: 0.04 L/h) and the concentration of available vitamin K (VKORC1 GA: -22% from normal vitamin K concentration; AA: -44% from normal vitamin K concentration). Virtual patient populations were used to assess the ability of the model to accurately predict routine INR and aPTT measurements from patients under long-term anticoagulant therapy. RESULTS: The introduced model accurately described the observed INR of patients receiving long-term warfarin treatment. The model was able to demonstrate the influence of genetic polymorphisms of CYP2C9 and VKORC1 on the INR measurements. Additionally, the model was successfully used to predict aPTT measurements for patients receiving long-term rivaroxaban therapy. CONCLUSION: The QSP model accurately predicted INR and aPTT measurements observed during routine therapeutic drug monitoring. This is an exemplar of how a QSP model can be adapted and used as a model-based precision dosing tool during clinical practice and drug development to predict efficacy and safety of anticoagulants to ultimately help optimize anti-thrombotic therapy.

A clinical population pharmacokinetic/pharmacodynamic model for BIIB059, a monoclonal antibody for the treatment of systemic and cutaneous lupus erythematosus.

A population pharmacokinetic/pharmacodynamic (popPK/PD) model for BIIB059 (anti-blood dendritic cell antigen 2 [anti-BDCA2]), a humanized immunoglobulin G1 monoclonal antibody currently under development for the treatment of SLE and CLE, is presented. BIIB059 binds BDCA2, a plasmacytoid dendritic cell (pDC)-specific receptor that inhibits the production of IFN-I and other inflammatory mediators when ligated. Phase 1 PK and PD data of healthy adult volunteers (HV, n = 87) and SLE subjects (n = 22) were utilized for the development of the popPK/PD model. The data included single and multiple dosing of intravenous and subcutaneous BIIB059. BDCA2 internalization (PD marker) was measured for all subjects by monitoring reduction of BDCA2 on pDC cell surface and used for development of the popPD model. A two-compartment popPK model with linear plus non-linear elimination was found to best describe BIIB059 PK. BDCA2 levels were best captured using an indirect response model with stimulation of the elimination of BDCA2. Clearance in SLE subjects was 25% higher compared to HV (6.87 vs 5.52 mL/h). Bodyweight was identified as only other covariate on clearance and central volume. The estimates of EC50 and Emax were 0.35 µg/mL and 8.92, respectively. No difference in EC50 and Emax was observed between SLE and HV. The popPK/PD model described the data accurately, as evaluated by pcVPCs and bootstrap. The presented popPK/PD model for BIIB059 provides valuable insight into the dynamics and dose-response relationship of BIIB059 for the treatment of SLE and CLE and was used to guide dose selection for the Phase 2 clinical study (NCT02847598).

Semi-mechanistic modelling platform to assess cardiac contractility and haemodynamics in preclinical cardiovascular safety profiling of new molecular entities.

BACKGROUND AND PURPOSE: Cardiovascular safety is one of the most frequent causes of safety-related attrition both preclinically and clinically. Preclinical cardiovascular safety is routinely assessed using dog telemetry monitoring key cardiovascular functions. The present research was to develop a semi-mechanistic modelling platform to simultaneously assess changes in contractility (dPdtmax ), heart rate (HR) and mean arterial pressure (MAP) in preclinical studies. EXPERIMENTAL APPROACH: Data from dPdtmax , HR, preload (left ventricular end-diastolic pressure [LVEDP]) and MAP were available from dog telemetry studies after dosing with atenolol (n = 27), salbutamol (n = 5), L-NG -nitroarginine methyl ester (L-NAME; n = 4), milrinone (n = 4), verapamil (n = 12), dofetilide (n = 8), flecainide (n = 4) and AZ001 (n = 14). Literature model for rat CV function was used for the structural population pharmacodynamic model development. LVEDP was evaluated as covariate to account for the effect of preload on dPdtmax . KEY RESULTS: The model was able to describe drug-induced changes in dPdtmax , HR and MAP for all drugs included in the developed framework adequately, by incorporating appropriate drug effects on dPdtmax , HR and/or total peripheral resistance. Consistent with the Starling's law, incorporation of LVEDP as a covariate on dPdtmax to correct for the preload effect was found to be statistically significant. CONCLUSIONS AND IMPLICATIONS: The contractility and haemodynamics semi-mechanistic modelling platform accounts for diurnal variation, drug-induced changes and inter-animal variation. It can be used to hypothesize and evaluate pharmacological effects and provide a holistic cardiovascular safety profile for new drugs.

Evaluating the Clinical Impact of Formulation Variability: A Metoprolol Extended-Release Case Study.

The objective of this research was to evaluate the impact of changes in the formulation of metoprolol extended-release (ER) tablets on dissolution, pharmacokinetic, and exercise-induced heart rate (EIHR) using a combined physiologically based absorption pharmacokinetic, and population pharmacokinetic/pharmacodynamic modeling and simulation approach. Using a previously developed physiologically based absorption pharmacokinetic model in DDDPlus and GastroPlus, we simulated the changes in drug release and exposure as the result of quantitative changes in the release-controlling excipient, hydroxylpropylmethylcellulose, for 50 and 200 mg. The similarity of dissolution profiles was assessed using the f2 test, and bioequivalence was tested on the simulated pharmacokinetic profiles. We used the simulated concentration-time profiles following formulation changes as pharmacokinetic input into a population pharmacokinetic/pharmacodynamic model newly developed in NONMEM to determine if changes in pharmacokinetics lead to clinically significant changes in pharmacodynamics. Pharmacodynamic assessment was based on the percentage reduction in the EIHR from baseline. Therapeutic effect was considered similar when the model-predicted EIHR was within 50% to 85% of the average maximum EIHR of healthy 30-year-old subjects. A 40% or more increase in the release rate constant resulted in dissimilarity in dissolution profiles and bioINequivalence in pharmacokinetics for both 50 and 200 mg. Formulation-related differences in drug release of metoprolol ER tablets can lead to differences in pharmacokinetics. However, the evaluated pharmacokinetic differences do not lead to clinically meaningful differences in EIHR, suggesting that EIHR may not be sensitive enough to detect changes in pharmacokinetics of metoprolol ER products.

Is Bioequivalence Established Based on the Reference-Scaled Average Bioequivalence Approach Relevant to Chronic Administration of Phenytoin? Perspectives Based on Population Pharmacokinetic Modeling and Simulations.

Phenytoin demonstrates time-dependent and nonlinear pharmacokinetics (PK) within the therapeutic range of 10 to 20 µg/mL. There are discussions on the relevance of bioequivalence (BE) demonstrated in a single-dose BE study in healthy subjects to exposure under chronic use conditions in patients, particularly given that phenytoin has a narrow therapeutic index. The objective of this study was to quantitatively evaluate the appropriateness of single-dose PK BE through simulations for the phenytoin extended-capsule products. A previously published population PK model was updated to account for the interoccasion variability using the dense PK data of the reference listed drug (Dilantin) from 5 single-dose, fully replicated BE studies (n = 124). BE studies with alternative designs were simulated using the developed PK model and subsequently analyzed accordingly: Scenario 1, multiple-dose, 2-period, crossover BE studies with an average BE approach; Scenario 2, single-dose, 4-period, fully replicated BE studies with a reference-scaled average BE approach as recommended in the product-specific guidance. In both scenarios, hypothetical phenytoin capsules with different formulation-related PK parameters, such as relative bioavailability and absorption rates, were included in the simulations. The results showed that the both scenarios provided the same results with respect to BE conclusions.

A pre-clinical quantitative model predicts the pharmacokinetics/pharmacodynamics of an anti-BDCA2 monoclonal antibody in humans.

BIIB059 is a novel humanized monoclonal antibody (mAb) that is currently under development for the treatment of Systemic Lupus Erythematosus and Cutaneous Lupus Erythematosus. BIIB059 is targeted against the blood dendritic cell antigen 2 (BDCA2), a receptor exclusively expressed on the surface of plasmacytoid dendritic cells (pDCs). Herein, we utilized pre-clinical pharmacokinetic (PK) and pharmacodynamic (PD) data to develop a non-human primate (NHP) model and to address whether the NHP model can be successfully scaled to predict the human PK/PD. In particular, PK data from 17 cynomolgus monkeys were utilized for PK model development, wherein BIIB059 was administered intravenously (1 and 10 mg/kg single-dosing and 5 mg/kg multiple-dosing) or subcutaneously (0.2 and 7.5 mg/kg single-dosing). Additionally, PD data (BDCA2 receptor density on pDCs) from 6 cynomolgus monkeys were used for the development of the PD model. The developed NHP two-compartment PK model, linked with an indirect response PD model, was subsequently scaled to humans by combining traditional allometric PK scaling with sensitivity-analysis-driven scaling of the PD. The scaled PK/PD model was then used to simulate the human PK/PD for different dose levels. When clinical data from the BIIB059 Phase I study became available, they were used to evaluate the predictability of the scaled PK/PD model and the model simulations were in agreement with the clinical data. Therefore, the presented approach is suggested to be employed in scaling pre-clinical mAb models to support the selection of safe first-in-human doses and, more broadly, the prediction of PK/PD in the clinic.

Model-Informed Precision Dosing at the Bedside: Scientific Challenges and Opportunities.

The development of model-informed precision dosing (MIPD) tools, especially in the form of native or web-based applications to be used at the bedside, has garnered marked attention in recent years. Their potential clinical benefit can be large, but it should be ensured that such tools make optimal use of available clinical data and have adequate predictive ability. Unique scientific challenges specific to MIPD remain, which may require adaptation of commonly used diagnostics in pharmacometrics.

Establishing a Multidisciplinary Framework to Study Drug-Drug Interactions of Hormonal Contraceptives: An Invitation to Collaborate.

Hormonal contraceptive agents (HCAs) are widely used throughout the world, and women taking HCAs are likely to take other medications. However, little is known about the clinical effect of most drug-drug interactions (DDIs) associated with HCAs. A team of interdisciplinary outcomes and pharmacometric researchers from academia and industry jointly engage in a research project to (i) quantitatively elucidate DDI impacts on unintended pregnancies and breakthrough bleeding, and (ii) establish a DDI-prediction framework to inform optimal use of HCAs.

A Semi-Mechanistic Population Pharmacokinetic Model of Nusinersen: An Antisense Oligonucleotide for the Treatment of Spinal Muscular Atrophy.

A pharmacokinetic (PK) model was developed for nusinersen, an antisense oligonucleotide (ASO) that is the first approved treatment for spinal muscular atrophy (SMA). The model was built with data from 92 nonhuman primates (NHPs) following intrathecal doses (0.3-7 mg) and characterized the PK in cerebrospinal fluid (CSF), plasma, total spinal cord, brain, and pons. The estimated volumes were 13.6, 937, 4.5, 53.8, and 2.11 mL, respectively. Global sensitivity analysis demonstrated that the CSF-to-plasma drug distribution rate (0.09 hour-1 ) is a major determinant of the maximum nusinersen concentration in central nervous system (CNS) tissues. Physiological age-based and body weight-based allometric scaling was implemented with exponent values of -0.08 and 1 for the rate constants and the volume of distribution, respectively. Simulations of the scaled model were in agreement with clinical observations from 52 pediatric phase I PK profiles. The developed model can be used to guide the design of clinical trials with ASOs.